Source author record

Evelyne Roueff

Evelyne Roueff 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

30works
6topics
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

30 published item(s)

preprint2025arXiv

PDRs4All XVII: Formation and excitation of HD in photodissociation regions. Application to the Orion Bar

The James Webb Space Telescope enabled the first detection of several rovibrational emission lines of HD in the Orion Bar, a prototypical photodissociation region. This provides an incentive to examine the physics of HD in dense and strong PDRs. Using the latest data available on HD excitation by collisional, radiative and chemical processes, our goal is to unveil HD formation and excitation processes in PDRs by comparing our state-of-the-art PDR model with observations made in the Orion Bar and discuss if and how HD can be used as a complementary tracer of physical parameters in the emitting region. We compute detailed PDR models, using an upgraded version of the Meudon PDR code, which are compared to NIRSpec data using excitation diagrams and synthetic emission spectra. The models predict that HD is mainly produced in the gas phase via the reaction D + H2 = H + HD at the front edge of the PDR and that the D/HD transition is located slightly closer to the edge than the H/H2 transition. Rovibrational levels are excited by UV pumping. In the observations, HD rovibrational emission is detected close to the H/H2 dissociation fronts of the Orion Bar and peaks where vibrationally excited H2 peaks, rather than at the maximum emission of pure rotational H2 levels. We derive an excitation temperature around Tex ~ 480 - 710 K. Due to high continuum in the Orion Bar, fringes lead to high noise levels beyond 15 $μ$m, no pure rotational lines of HD are detected. The comparison to PDR models shows that a range of thermal pressure P = (3-9)x10$^7$ K cm$^{-3}$ with no strong constraints on the intensity of the UV field are compatible with HD observations. This range of pressure is compatible with previous estimates from H2 observations with JWST. This is the first time that observations of HD emission lines in the near-infrared are used to put constraints on the thermal pressure in the PDR.

preprint2024arXiv

Formation of the Methyl Cation by Photochemistry in a Protoplanetary Disk

Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation.

preprint2023arXiv

Analysis of the first infrared spectrum of quasi-bound H2 line emission in Herbig-Haro 7

Context. Highly excited molecular hydrogen H2 has been observed in many regions of shocked molecular gas. A recently published $K$-band spectrum of Herbig-Haro 7 (HH7) contains several vibration-rotation lines of H2 from highly excited energy levels that have not been detected elsewhere, including a line at 2.179 $μ$m identified as arising from the $v$=2 $J$=29 level, which lies above the dissociation limit of H2. One emission line at 2.104 $μ$m in this spectrum was unidentified. Aims. We aim to complete the analysis of the spectrum of HH7 by including previously missing molecular data that have been recently computed. Methods. We re-analysed the $K$-band spectrum, emphasising the physics of quasi-bound upper levels that can produce infrared emission lines in the $K$ band. Results. We confirm the identification of the $2-1$ $S$(27) line at 2.1785 $μ$m and identify the line at 2.1042 $μ$m as due to the 1-0 $S$(29) transition of H2, whose upper level energy is also higher than the dissociation limit. This latter identification, its column density, and the energy of its upper level further substantiate the existence of a hot thermal component at 5000 K in the HH7 environment.} Conclusion. The presence of the newly identified $1-0$ $S$(29) line, whose quasi-bound upper level ($v$=1, $J$=31) has a significant spontaneous dissociation probability, shows that dissociation of H2 is occurring. The mechanism by which virtually all of the H2 levels with energies from 20,000 K to 53,000 K is maintained in local thermodynamic equilibrium at a single temperature of $\sim$5,000 K remains to be understood.

preprint2022arXiv

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.

preprint2021arXiv

Space and laboratory observation of the deuterated cyanomethyl radical HDCCN

Our observations of TMC-1 with the Yebes 40 m radio telescope in the 31.0-50.3 GHz range allowed us to detect a group of unidentified lines, showing a complex line pattern indicative of an open-shell species. {}The observed frequencies of these lines and the similarity of the spectral pattern with that of the 2$_{0,2}$-1$_{0,1}$ rotational transition of H$_2$CCN indicate that the lines arise from the deuterated cyanomethyl radical, HDCCN. Using Fourier transform microwave spectroscopy experiments combined with electric discharges, we succeeded in producing the radical HDCCN in the laboratory and observed its 1$_{0,1}$-0$_{0,0}$ and 2$_{0,2}$-1$_{0,1}$ rotational transitions. From our observations and assuming a rotational temperature of 5 K, we derive an abundance ratio H$_2$CCN/HDCCN=20$\pm$4. The high abundance of the deuterated form of H$_2$CCN is well accounted for by a standard gas-phase model, in which deuteration is driven by deuteron transfer from the H$_2$D$^+$ molecular ion.

preprint2020arXiv

C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: An analysis of the precision achievable when modeling spectral line within the Local Thermodynamic Equilibrium approximation

CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

preprint2020arXiv

Quantitative inference of the $H_2$ column densities from 3 mm molecular emission: A case study towards Orion B

Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is non-linear and sensitive to excitation conditions, optical depths, abundance variations due to the underlying physico-chemistry. We aim to use multi-molecule line emission to infer NH2 from radio observations. We propose a data-driven approach to determine NH2 from radio molecular line observations. We use supervised machine learning methods (Random Forests) on wide-field hyperspectral IRAM-30m observations of the Orion B molecular cloud to train a predictor of NH2, using a limited set of molecular lines as input, and the Herschel-based dust-derived NH2 as ground truth output. For conditions similar to the Orion B molecular cloud, we obtain predictions of NH2 within a typical factor of 1.2 from the Herschel-based estimates. An analysis of the contributions of the different lines to the predictions show that the most important lines are $^{13}$CO(1-0), $^{12}$CO(1-0), C$^{18}$O(1-0), and HCO$^+$(1-0). A detailed analysis distinguishing between diffuse, translucent, filamentary, and dense core conditions show that the importance of these four lines depends on the regime, and that it is recommended to add the N$_2$H$^+$(1-0) and CH$_3$OH(20-10) lines for the prediction of NH2 in dense core conditions. This article opens a promising avenue to directly infer important physical parameters from the molecular line emission in the millimeter domain. The next step will be to try to infer several parameters simultaneously (e.g., NH2 and far-UV illumination field) to further test the method. [Abridged]

preprint2020arXiv

Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids

The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models. We build grids of models that sample randomly a large space of physical conditions (unobservable quantities such as gas density, temperature, etc.) and compute the corresponding observables (line intensities, column densities) and the ionization fraction. We estimate the predictive power of each potential tracer by training a Random Forest model to predict the ionization fraction from that tracer, based on these model grids. In both translucent medium and cold dense medium conditions, several observable tracers with very good predictive power for the ionization fraction are found. Several tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.). (abridged)

preprint2019arXiv

The role of C/O in nitrile astrochemistry in PDRs and planet-forming disks

Complex nitriles, such as HC3N, and CH3CN, are observed in a wide variety of astrophysical environments, including at relatively high abundances in photon-dominated regions (PDR) and the UV exposed atmospheres of planet-forming disks. The latter have been inferred to be oxygen-poor, suggesting that these observations may be explained by organic chemistry in C-rich environments. In this study we first explore if the PDR complex nitrile observations can be explained by gas-phase PDR chemistry alone if the elemental C/O ratio is elevated. In the case of the Horsehead PDR, we find that gas-phase chemistry with C/O $\gtrsim$ 0.9 can indeed explain the observed nitrile abundances, increasing predicted abundances by several orders of magnitude compared to standard C/O assumptions. We also find that the nitrile abundances are sensitive to the cosmic ray ionization treatment, and provide constraints on the branching ratios between CH3CN and CH3NC productions. In a fiducial disk model, an elevated C/O ratio increases the CH3CN and HC3N productions by more than an order of magnitude, bringing abundance predictions within an order of magnitude to what has been inferred from observations. The C/O ratio appears to be a key variable in predicting and interpreting complex organic molecule abundances in photon-dominated regions across a range of scales.

preprint2016arXiv

Radiative and mechanical feedback into the molecular gas in the Large Magellanic Cloud. I. N159W

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

preprint2015arXiv

HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions

We present new analytic theory and radiative transfer computations for the atomic to molecular (HI-to-H2) transitions, and the build-up of atomic-hydrogen (HI) gas columns, in optically thick interstellar clouds, irradiated by far-ultraviolet photodissociating radiation fields. We derive analytic expressions for the total HI column densities for (1D) planar slabs, for beamed or isotropic radiation fields, from the weak- to strong-field limits, for gradual or sharp atomic to molecular transitions, and for arbitrary metallicity. Our expressions may be used to evaluate the HI column densities as functions of the radiation field intensity and the H2-dust-limited dissociation flux, the hydrogen gas density, and the metallicity-dependent H2 formation rate-coefficient and far-UV dust-grain absorption cross-section. We make the distinction between "HI-dust" and "H2-dust" opacity, and we present computations for the "universal H2-dust-limited effective dissociation bandwidth". We validate our analytic formulae with Meudon PDR code computations for the HI-to-H2 density profiles, and total HI column densities. We show that our general 1D formulae predict HI columns and H2 mass fractions that are essentially identical to those found in more complicated (and approximate) spherical (shell/core) models. We apply our theory to compute H2 mass fractions and star-formation thresholds for individual clouds in self-regulated galaxy disks, for a wide range of metallicities. Our formulae for the HI columns and H2 mass fractions may be incorporated into hydrodynamics simulations for galaxy evolution.

preprint2015arXiv

HI-to-H2 transitions in the Perseus Molecular Cloud

We apply the Sternberg et al. (2014) theoretical model to analyze HI and H2 observations in the Perseus molecular cloud. We constrain the physical properties of the HI shielding envelopes and the nature of the HI-to-H2 transitions. Our analysis (Bialy et al. 2015) implies that in addition to cold neutral gas (CNM), less dense thermally-unstable gas (UNM) significantly contributes to the shielding of the H2 cores in Perseus.

preprint2015arXiv

HI-to-H2 Transitions in the Perseus Molecular Cloud

We use the Sternberg et al. (2014) theory for interstellar atomic to molecular (HI-to-H$_2$) conversion to analyze HI-to-H$_2$ transitions in five (low-mass) star-forming and dark regions in the Perseus molecular cloud, B1, B1E, B5, IC348, and NGC1333. The observed HI mass surface densities of 6.3 to 9.2 M$_{\odot}$ pc$^{-2}$ are consistent with HI-to-H$_2$ transitions dominated by HI-dust shielding in predominantly atomic envelopes. For each source, we constrain the dimensionless parameter $αG$, and the ratio $I_{\rm UV}/n$, of the FUV intensity to hydrogen gas density. We find $αG$ values from 5.0 to 26.1, implying characteristic atomic hydrogen densities 11.8 to 1.8 cm$^{-3}$, for $I_{\rm UV} \approx 1$ appropriate for Perseus. Our analysis implies that the dusty HI shielding layers are probably multiphased, with thermally unstable UNM gas in addition to cold CNM within the 21 cm kinematic radius.

preprint2015arXiv

Physical conditions in the central molecular zone inferred by H3+

The H3+ molecule has been detected in many lines of sight within the central molecular zone (CMZ) with exceptionally large column densities and unusual excitation properties compared to diffuse local clouds. The detection of the (3,3) metastable level has been suggested to be the signature of warm and diffuse gas in the CMZ. We use the Meudon PDR code to re-examine the relationship between the column density of H3+ and the cosmic-ray ionization rate, $ζ$, up to large values of $ζ$. We study the impact of the various mechanisms that can excite H3+ in its metastable state. We produce grids of PDR models exploring different parameters ($ζ$, size of clouds, metallicity) and infer the physical conditions that best match the observations toward ten lines of sight in the CMZ. For one of them, Herschel observations of HF, OH+, H2O+, and H3O+ can be used as additional constraints. We check that the results found for H3+ also account for the observations of these molecules. We find that the linear relationship between N(H3+) and $ζ$ only holds up to a certain value of the cosmic-ray ionization rate, which depends on the proton density. A value $ζ\sim 1 - 11 \times 10^{-14}$ s$^{-1}$ explains both the large observed H3+ column density and its excitation in the metastable level (3,3) in the CMZ. It also reproduces N(OH+), N(H2O+) and N(H3O+) detected toward Sgr B2(N). We confirm that the CMZ probed by H3+ is diffuse, nH $\lesssim$ 100 cm-3 and warm, T $\sim$ 212-505 K. This warm medium is due to cosmic-ray heating. We also find that the diffuse component probed by H3+ must fill a large fraction of the CMZ. Finally, we suggest the warm gas in the CMZ enables efficient H2 formation via chemisorption sites as in PDRs. This contributes to enhance the abundance of H3+ in this high cosmic-ray flux environment.

preprint2015arXiv

The interstellar chemistry of H2C3O isomers

We present the detection of two H2C3O isomers, propynal and cyclopropenone, toward various starless cores and molecular clouds, together with upper limits for the third isomer propadienone. We review the processes controlling the abundances of H2C3O isomers in interstellar media showing that the reactions involved are gas-phase ones. We show that the abundances of these species are controlled by kinetic rather than thermodynamic effects.

preprint2014arXiv

Ion-molecule reactions involving HCO$^+$ and N$_2$H$^+$: Isotopologue equilibria from new theoretical calculations and consequences for interstellar isotope fractionation

$Aims$: We revisit with new augmented accuracy the theoretical dynamics of basic isotope exchange reactions involved in the $^{12}$C/$^{13}$C, $^{16}$O/$^{18}$O, and $^{14}$N/$^{15}$N balance because these reactions have already been studied experimentally in great detail. $Methods$: Electronic structure methods were employed to explore potential energy surfaces, full-dimensional rovibrational calculations to compute rovibrational energy levels that are numerically exact, and chemical network models to estimate the abundance ratios under interstellar conditions. $Results$: New exothermicities, derived for HCO$^+$ reacting with CO, provide rate coefficients markedly different from previous theoretical values in particular at low temperatures, resulting in new abundance ratios relevant for carbon chemistry networks. In concrete terms, we obtain a reduction in the abundance of H$^{12}$C$^{18}$O$^+$ and an increase in the abundance of H$^{13}$C$^{16}$O$^+$ and D$^{13}$C$^{16}$O$^+$. In all studied cases, the reaction of the ion with a neutral polarizable molecule proceeds through the intermediate proton-bound complex found to be very stable. For the complexes OCH$^+$...CO, OCH$^+$...OC, COHOC$^+$, N$_2$...HCO$^+$, N$_2$H$^+$...OC, and N$_2$HN$_2^+$, we also calculated vibrational frequencies and dissociation energies. $Conclusions$: The linear proton-bound complexes possess sizeable dipole moments, which may facilitate their detection.

preprint2013arXiv

An observational investigation of the identity of B11244 ($l$-C$_3$H$^+$/C$_3$H$^-$)

Pety et al. (2012) have reported the detection of eight transitions of a closed-shell, linear molecule (B11244) in observations toward the Horsehead PDR, which they attribute to the $l$-C$_3$H$^+$ cation. Recent high-level \textit{ab initio} calculations have called this assignment into question; the anionic C$_3$H$^-$ molecule has been suggested as a more likely candidate. Here, we examine observations of the Horsehead PDR, Sgr B2(N), TMC-1, and IRC+10216 in the context of both $l$-C$_3$H$^+$ and C$_3$H$^-$. We find no observational evidence of $K_a = 1$ lines which should be present were the carrier indeed C$_3$H$^-$. Additionally, we find a strong anti-correlation between the presence of known molecular anions and B11244 in these regions. Finally, we discuss the formation and destruction chemistry of C$_3$H$^-$ in the context of the physical conditions in the regions. Based on these results, we conclude there is little evidence to support the carrier is C$_3$H$^-$.

preprint2012arXiv

Herschel observations of interstellar chloronium

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

preprint2012arXiv

Surface chemistry in the interstellar medium - I - H2 formation by Langmuir-Hinshelwood and Eley-Rideal mechanisms

H2 formation remains a major issue for the understanding of interstellar physics. We investigate H2 formation in the interstellar medium at the light of the most recent experimental and theoretical data. We implemented detailed H2 formation mechanisms on grains surface in the Meudon PDR code : i) Langmuir-Hinshelwood mechanism taking into account the contribution of the different sizes of dust grains in the diffusion processes and ii) the Eley-Rideal mechanism. We show that, thanks to these processes, H2 can be formed even in regions where dust temperature is larger than 25 K. We also show that formation by Eley-Rideal mechanism can be a significant source of heating of the gas. We derive line intensities for various astrophysical conditions. Such an approach results in an enhanced H2 formation rate compared to the standard formation determined by observations in absorption in the UV. We derive some H2 line intensities for isobaric and isochoric models.

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

preprint2011arXiv

NUV Excess in Slowly Accreting T Tauri Stars: Limits Imposed by Chromospheric Emission

Young stars surrounded by disks with very low mass accretion rates are likely in the final stages of inner disk evolution and therefore particularly interesting to study. We present ultraviolet (UV) observations of the ~5-9 Myr old stars RECX-1 and RECX-11, obtained with the Cosmic Origins Spectrograph (COS) and Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST), as well as optical and near infrared spectroscopic observations. The two stars have similar levels of near UV emission, although spectroscopic evidence indicates that RECX-11 is accreting and RECX-1 is not. The line profiles of Halpha and He I 10830 in RECX-11 show both broad and narrow redshifted absorption components that vary with time, revealing the complexity of the accretion flows. We show that accretion indicators commonly used to measure mass accretion rates, e.g. U band excess luminosity or the Ca II triplet line luminosity are unreliable for low accretors, at least in the middle K spectral range. Using RECX-1 as a template for the intrinsic level of photospheric and chromospheric emission, we determine an upper limit of 3 x 10^-10 solar masses/ year for RECX-11. At this low accretion rate, recent photoevaporation models predict that an inner hole should have developed in the disk. However, the spectral energy distribution of RECX-11 shows fluxes comparable to the median of Taurus in the near infrared, indicating that substantial dust remains. Fluorescent H_2 emission lines formed in the innermost disk are observed in RECX-11, showing that gas is present in the inner disk, along with the dust.

preprint2011arXiv

Ortho-H2 and the Age of Interstellar Dark Clouds

Interstellar dark clouds are the sites of star formation. Their main component, dihydrogen, exists under two states, ortho and para. H2 is supposed to form in the ortho:para ratio (OPR) of 3:1 and to subsequently decay to almost pure para-H2 (OPR <= 0.001). Only if the H2 OPR is low enough, will deuteration enrichment, as observed in the cores of these clouds, be efficient. The second condition for strong deuteration enrichment is the local disappearance of CO, which freezes out onto grains in the core formation process. We show that this latter condition does not apply to DCO+, which, therefore, should be present all over the cloud. We find that an OPR >= 0.1 is necessary to prevent DCO+ large-scale apparition. We conclude that the inevitable decay of ortho-H2 sets an upper limit of ~6 million years to the age of starless molecular clouds under usual conditions.

preprint2010arXiv

Collisional excitation of doubly deuterated ammonia NHD2 by para-H2

Collisional de-excitation rates of partially deuterated molecules are different from the fully hydrogenated species because of lowering of symmetry. We compute the collisional (de)excitation rates of ND2H by ground state para-H2, extending the previous results for He- lium. We describe the changes in the potential energy surface of NH3- H2 involved by the pres- ence of two deuterium nuclei. Cross sections are calculated within the full close-coupling ap- proach and augmented with coupled-state calculations. Collisional rate coefficients are given between 5 and 35 K, a range of temperatures which is relevant to cold interstellar conditions. We find that the collisional rates of ND2H by H2 are about one order of magnitude higher than those obtained with Helium as perturber. These results are essential to radiative transfer modelling and will allow to interpret the millimeter and submillimeter detections of ND2H with better constraints than previously.

preprint2009arXiv

Detection of 15NH2D in dense cores: A new tool for measuring the 14N/15N ratio in the cold ISM

Ammonia is one of the best tracers of cold dense cores. It is also a minor constituent of interstellar ices and, as such, one of the important nitrogen reservoirs in the protosolar nebula, together with the gas phase nitrogen, in the form of N2 and N. An important diagnostic of the various nitrogen sources and reservoirs of nitrogen in the Solar System is the 14N/15N isotopic ratio. While good data exist for the Solar System, corresponding measurements in the interstellar medium are scarce and of low quality. Following the successful detection of the singly, doubly, and triply deuterated isotopologues of ammonia, we have searched for 15NH2D in dense cores, as a new tool for investigating the 14N/15N ratio in dense molecular gas. With the IRAM-30m telescope, we have obtained deep integrations of the ortho 15NH2D (1(1,1)-1(0,1)) line at 86.4 GHz, simultaneously with the corresponding ortho NH2D line at 85.9 GHz. o-15NH2D is detected in Barnard-1b, NGC1333-DCO+, and L1689N, while we obtained upper limits towards LDN1544 and NGC1333-IRAS4A, and a tentative detection towards L134N(S). The 14N/15N abundance ratio in NH2D ranges between 350 and 850, similar to the protosolar value of ~ 424, and likely higher than the terrestrial ratio of 270.

preprint2009arXiv

Far-Ultraviolet H2 Emission from Circumstellar Disks

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 Myr old Taurus-Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion OB1a complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly to electron impact excitation of \h2 to study the evolution of the gas in the inner disk. We find that the \h2 feature is absent in non-accreting sources, but is detected in the spectra of CTTS and correlates with accretion luminosity. Since all young stars have active chromospheres which produce strong X-ray and UV emission capable of exciting \h2 in the disk, the fact that the non-accreting sources show no \h2 emission implies that the \h2 gas in the inner disk has dissipated in the non-accreting sources, although dust (and possibly gas) remains at larger radii. Using the flux at 1600 Å, we estimate that the column density of \h2 left in the inner regions of the debris disks in our sample is less than ~ 3x10^-6 g cm^-2, nine orders of magnitude below the surface density of the minimum mass solar nebula at 1 AU.

preprint2006arXiv

A model for atomic and molecular interstellar gas: The Meudon PDR code

We present the revised ``Meudon'' model of Photon Dominated Region (PDR code), presently available on the web under the Gnu Public Licence at: http://aristote.obspm.fr/MIS. General organisation of the code is described down to a level that should allow most observers to use it as an interpretation tool with minimal help from our part. Two grids of models, one for low excitation diffuse clouds and one for dense highly illuminated clouds, are discussed, and some new results on PDR modelisation highlighted.

preprint2006arXiv

Benchmarking PDR models against the Horsehead edge

To prepare for the unprecedented spatial and spectral resolution provided by ALMA and Herschel/HIFI, chemical models are being benchmarked against each other. It is obvious that chemical models also need well-constrained observations that can serve as references. Photo-dissociation regions (PDRs) are particularly well suited to serve as references because they make the link between diffuse and molecular clouds, thus enabling astronomers to probe a large variety of physical and chemical processes. At a distance of 400 pc (1" corresponding to 0.002 pc), the Horsehead PDR is very close to the prototypical kind of source (i.e. 1D, edge-on) needed to serve as a reference to models.

preprint2006arXiv

The distribution of ND2H in LDN1689N

Finding tracers of the innermost regions of prestellar cores is important for understanding their chemical and dynamical evolution before the onset of gravitational collapse. While classical molecular tracers, such as CO and CS, have been shown to be strongly depleted in cold, dense gas by condensation on grain mantles, it has been a subject of discussion to what extent nitrogen-bearing species, such as ammonia, are affected by this process. As deuterium fractionation is efficient in cold, dense gas, deuterated species are excellent tracers of prestellar cores. A comparison of the spatial distribution of neutral and ionized deuterated species with the dust continuum emission can thus provide important insights into the physical and chemical structure of such regions. We study the spatial distribution of the ground-state 335.5 GHz line of ND2H in LDN1689N, using APEX, and compare it with the distribution of the DCO+(3--2) line, as well as the 350 micron dust continuum emission observed with the SHARC~II bolometer camera at CSO. While the distribution of the ND2H emission in LDN1689N is generally similar to that of the 350 microns dust continuum emission, the peak of the ND2H emission is offset by ~10'' to the East from the dust continuum and DCO+ emission peak. ND2H and ND3 share the same spatial distribution. The observed offset between the ND2H and DCO+ emission is consistent with the hypothesis that the deuterium peak in LDN1689N is an interaction region between the outflow shock from IRAS16293--2422 and the dense ambient gas. We detect the J = 4 - 3 line of H13CO+ at 346.998 GHz in the image side band serendipitously. This line shows the same spatial distribution as DCO+(3--2), and peaks close to the 350 mic emission maximum which provides further support for the shock interaction scenario.

preprint2001arXiv

Small scale structure in diffuse molecular gas from repeated FUSE and visible spectra of HD 34078

We present preliminary results from an ongoing program devoted to a study of small scale structure in the spatial distribution of molecular gas. Our work is based on multi-epoch FUSE and visible observations of HD34078. A detailed comparison of H2, CH and CH+ absorption lines is performed. No short term variations are seen (except for highly excited H2) but long-term changes in N(CH) are clearly detected when comparing our data to spectra taken about 10 years ago.