Source author record

F. Lique

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

18works
4topics
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

18 published item(s)

preprint2022arXiv

A SUBLIME 3D Model for Cometary Coma Emission: the Hypervolatile-Rich Comet C/2016 R2 (PanSTARRS)

The coma of comet C/2016 R2 (PanSTARRS) is one of the most chemically peculiar ever observed, in particular due to its extremely high CO/H2O and N2+/H2O ratios}, and unusual trace volatile abundances. However, the complex shape of its CO emission lines, as well as uncertainties in the coma structure and excitation, has lead to ambiguities in the total CO production rate. We performed high resolution, spatially, spectrally and temporally resolved CO observations using the James Clerk Maxwell Telescope (JCMT) and Submillimeter Array (SMA) to elucidate the outgassing behaviour of C/2016 R2. Results are analyzed using a new, time-dependent, three dimensional radiative transfer code (SUBLIME), incorporating for the first time, accurate state-to-state collisional rate coefficients for the CO--CO system. The total CO production rate was found to be in the range $(3.8-7.6)\times10^{28}$ s$^{-1}$ between 2018-01-13 and 2018-02-01, with a mean value of $(5.3\pm0.6)\times10^{28}$ s$^{-1}$ at r_H = 2.8-2.9 au. The emission is concentrated in a near-sunward jet, with an outflow velocity $0.51\pm0.01$ km/s, compared to $0.25\pm0.01$ km/s in the ambient (and night-side) coma. Evidence was also found for an extended source of CO emission, possibly due to icy grain sublimation around $1.2\times10^5$ km from the nucleus. Based on the coma molecular abundances, we propose that the nucleus ices of C/2016 R2 can be divided into a rapidly sublimating apolar phase, rich in CO, CO2, N2 and CH3OH, and a predominantly frozen (or less abundant), polar phase containing more H2O, CH4, H2CO and HCN.

preprint2022arXiv

An improved study of HCO+ and He system: interaction potential, collisional relaxation and pressure broadening

In light of its ubiquitous presence in the interstellar gas, the chemistry and reactivity of the HCO+ ion requires special attention. The availability of up-to-date collisional data between this ion and the most abundant perturbing species in the interstellar medium is a critical resource in order to derive reliable values of its molecular abundance from astronomical observations. This work intends to provide improved scattering parameters for the HCO+ and He collisional system. We have tested the accuracy of explicitly correlated coupled-cluster methods for mapping the short- and long-range multi-dimensional potential energy surface of atom-ion systems. A validation of the methodology employed for the calculation of the potential well has been obtained from the comparison with experimentally derived bound-state spectroscopic parameters. Finally, by solving the close-coupling scattering equations, we have derived the pressure broadening and shift coefficients for the first six rotational transitions of HCO+ as well as inelastic state-to-state transition rates up to j = 5 in the 5-100 K temperature interval.

preprint2022arXiv

Collisional excitation and non-LTE modelling of interstellar chiral propylene oxide

The first set of theoretical cross sections for propylene oxide (CH3CHCH2O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H2 are deduced in the temperature range 5-30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made towards the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (~ 10 K) and translucent (nH ~ 2000 cm-3) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is Ntot ~ 3e12 cm-2, corresponding to a fractional abundance relative to total hydrogen of ~ 2.5e-11. The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.

preprint2021arXiv

Anomalous HCN emission from warm giant molecular clouds

HCN is considered a good tracer of the dense molecular gas that serves as fuel for star formation. However, recent large-scale surveys of giant molecular clouds (GMCs) have detected extended HCN line emission. Such observations often resolve the HCN J=1-0 hyperfine structure (HFS). A precise determination of the physical conditions of the gas requires treating the HFS line overlap effects. Here, we study the HCN HFS excitation and line emission using nonlocal radiative transfer models that include line overlaps and new HFS-resolved collisional rate coefficients for inelastic collisions of HCN with both para-H2 and ortho-H2 (computed via the scaled-IOS approximation up to Tk=500 K). In addition, we account for the role of electron collisions in the HFS level excitation. We find that line overlap and opacity effects frequently produce anomalous HCN J=1-0 HFS line intensity ratios (inconsistent with the common assumption of the same Tex for all HFS lines) as well as anomalous HFS line width ratios. Line overlap and electron collisions also enhance the excitation of the higher J rotational lines. Electron excitation becomes important for molecular gas with H2 densities below a few 10^5 cm-3 and electron abundances above ~10^-5. In particular, electron excitation can produce low-surface-brightness HCN emission from very extended but low-density gas in GMCs. The existence of such a widespread HCN emission component may affect the interpretation of the extragalactic relationship HCN luminosity versus star-formation rate. Alternatively, extended HCN emission may arise from dense star-forming cores and become resonantly scattered by large envelopes of lower density gas. There are two scenarios - namely, electron-assisted (weakly) collisionally excited versus scattering - that lead to different HCN J=1-0 HFS intensity ratios, which can be tested on the basis of observations.

preprint2021arXiv

Sulfur gas-phase abundance in dense cores

The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. However, atomic sulfur, which could be the main gas-phase carrier, cannot be observed directly in cold cores. We present measurements of the nitrogen sulfide (NS) radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that $n(\rm S)/n_{\rm H}$ is close, or even equal, to the cosmic abundance of sulfur, 14$\times 10^{-6}$, demonstrating that sulfur is not depleted and is atomic, which is in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime also need to be explored from their chemical and dynamical perspectives.

preprint2020arXiv

The effect of CO-H2O collisions in the rotational excitation of cometary CO

We present the first accurate rate coefficients for the rotational excitation of CO by H2O in the kinetic temperature range 5-100 K. The statistical adiabatic channel method (SACM) is combined with a high-level rigid-rotor CO-H2O intermolecular potential surface. Transitions among the first 11 rotational levels of CO and the first 8 rotational levels of both para-H2O and ortho-H2O are considered. Our rate coefficients are compared to previous data from the literature and they are also incorporated in a simple non-LTE model of cometary coma including collision-induced transitions, solar radiative pumping and radiative decay. We find that the uncertainties in the collision data have significant influence on the CO population distribution for H2O densities in the range 10^3-10^8 cm^-3. We also show that the rotational distribution of H2O plays an important role in CO excitation (owing to correlated energy transfer in both CO and H2O), while the impact of the ortho-to-para ratio of H2O is found to be negligible.

preprint2019arXiv

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

preprint2016arXiv

Collisional excitation of HC3N by para- and ortho-H2

New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks.

preprint2016arXiv

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

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

preprint2016arXiv

N$_2$H$^+$ and N$^{15}$NH$^+$ towards the prestellar core 16293E in L1689N

Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statistically significant studies. In particular, the N$_2$H$^+$ molecule currently shows the largest spread of $^{14}$N/$^{15}$N ratios in high-mass star forming regions. However, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ was obtained in only two low-mass star forming regions (L1544 and B1b). The current work extends this sample to a third dark cloud. We targeted the 16293E prestellar core, where the N$^{15}$NH$^+$ $J$=1-0 line was detected. Using a model previously developed for the physical structure of the source, we solved the molecular excitation with a non-local radiative transfer code. For that purpose, we computed specific collisional rate coefficients for the N$^{15}$NH$^+$-H$_2$ collisional system. As a first step of the analysis, the N$_2$H$^+$ abundance profile was constrained by reproducing the N$_2$H$^+$ $J$=1-0 and 3-2 maps. A scaling factor was then applied to this profile to match the N$^{15}$NH$^+$ $J$=1-0 spectrum. We derive a column density ratio N$_2$H$^+$ / N$^{15}$NH$^+$ = $330^{+170}_{-100}$. The current estimate $\sim$330 agrees with the value typical of the elemental isotopic ratio in the local ISM. It is however lower than in some other cores, where values as high as 1300 have been reported.

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

Stratified NH and ND emission in the prestellar core 16293E in L1689N

High degrees of deuterium fractionation are commonly found in cold prestellar cores and in the envelopes around young protostars. As it brings strong constraints to chemical models, deuterium chemistry is often used to infer core history or molecule formation pathways. Whereas a large number of observations is available regarding interstellar deuterated stable molecules, relatively little is known about the deuteration of hydride radicals, as their fundamental rotational transitions are at high frequencies where the atmosphere is mostly opaque. Nitrogen hydride radicals are important species in nitrogen chemistry, as they are thought to be related to ammonia formation. Observations have shown that ammonia is strongly deuterated, with [NH2D]/[NH3] ~ 10%. Models predict similarly high [ND]/[NH] ratios, but so far only one observational determination of this ratio is available, towards the envelope of the protostar IRAS16293-2422. In order to test model predictions, we aim here at determining [ND]/[NH] in a dense, starless core. We observed NH and ND in 16293E with the HIFI spectrometer on board the Herschel Space Observatory as part of the CHESS guaranteed time key programme, and derived the abundances of these two species using a non-LTE non-local radiative transfer model. Both NH and ND are detected in the source, with ND in emission and NH in absorption against the continuum arising from the cold dust emission. Our model shows however that the ND emission and the NH absorption originate from different layers in the cloud, as further evidenced by their different velocities. In the central region of the core, we can set a lower limit to the [ND]/[NH] ratio of ~2%. This estimate is consistent with recent pure gas-phase models of nitrogen chemistry

preprint2014arXiv

Collisional excitation of water by hydrogen atoms

We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of both ortho- and para-H$_2$O and for temperatures in the range T = 5-1500 K. These rate coefficients are subsequently compared to the values previously published for the H$_2$O / He and H$_2$O / H$_2$ collisional systems. It is shown that no simple relation exists between the three systems and that specific calculations are thus mandatory.

preprint2014arXiv

Depletion of chlorine into HCl ice in a protostellar core

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

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

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.

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.

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.