Source author record

Javier R. Goicoechea

Javier R. Goicoechea 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

29works
9topics
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

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

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.

preprint2022arXiv

The role of highly vibrationally excited H2 initiating the N chemistry: Quantum study and 3-sigma detection of NH emission in the Orion Bar PDR

The formation of hydrides by gas-phase reactions between H2 and a heavy element atom is a very selective process. Reactions with ground-state neutral carbon, oxygen, nitrogen, and sulfur atoms are very endoergic and have high energy barriers because the H2 molecule has to be fragmented before a hydride bond is formed. In cold interstellar clouds, these barriers exclude the formation of CH, OH, NH, and SH radicals through hydrogen abstraction reactions. Here we study a very energetically unfavorable process, the reaction of N(4S) atoms with H2 molecules. We calculated the reaction rate coefficient for H2 in different vibrational levels, using quantum methods for v=0-7 and quasi-classical methods up to v=12. Owing to the high energy barrier, these rate coefficients increase with v and also with the gas temperature. We implemented the new rates in the Meudon PDR code and studied their effect on models with different ultraviolet (UV) illumination conditions. In strongly UV-irradiated dense gas (Orion Bar conditions), the presence of H2 in highly vibrationally excited levels (v>7) enhances the NH abundance by two orders of magnitude (at the PDR surface) compared to models that use the thermal rate coefficient for reaction N(4S) + H2 -> NH + H. The increase in NH column density across the PDR is a factor of ~25. We explore existing Herschel/HIFI observations of the Orion Bar and Horsehead PDRs. We report a 3-sigma emission feature at the ~974 GHz frequency of the NH N_J=1_2-0_1 line toward the Bar. The emission level implies N(NH)~10^13 cm^-2, which is consistent with PDR models using the new rate coefficients for reactions between N and UV-pumped H2. This formation route dominates over hydrogenation reactions involving the less abundant N+ ion. JWST observations will soon quantify the amount and reactivity of UV-pumped H2 in many interstellar and circumstellar environments.

preprint2020arXiv

Bringing high spatial resolution to the Far-infrared -- A giant leap for astrophysics

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

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

CF$^+$ excitation in the interstellar medium

The detection of CF$^+$ in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF$^+$ abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF$^+$ in collisions with para- and ortho-H2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF$^+$ emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF$^+$ as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

preprint2020arXiv

Distribution of Water Vapor in Molecular Clouds. II

The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.

preprint2020arXiv

Hyperfine excitation of SH$^+$ by H

SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of the physical conditions prevailing in these energetic environments relies on knowing the rate coefficients of inelastic collisions between SH$^+$ molecules and hydrogen atoms, hydrogen molecules, and electrons. In this paper, we derive SH$^+$--H fine and hyperfine-resolved rate coefficients from the recent quantum calculations for the SH$^+$--H collisions, including inelastic, exchange and reactive processes. The method used is based on the infinite order sudden approach. State-to-state rate coefficients between the first 31 fine levels and 61 hyperfine levels of SH$^+$ were obtained for temperatures ranging from 10 to 1000 K. Fine structure-resolved rate coefficients present a strong propensity rule in favour of $Δj = ΔN$ transitions. The $Δj = ΔF$ propensity rule is observed for the hyperfine transitions. {The new rate coefficients will help significantly in the interpretation of SH$^+$ spectra from PDRs and UV-irradiated shocks where the abundance of hydrogen atoms with respect to hydrogen molecules can be significant.

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)

preprint2016arXiv

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

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

preprint2016arXiv

Interstellar Hydrides

Interstellar hydrides -- that is, molecules containing a single heavy element atom with one or more hydrogen atoms -- were among the first molecules detected outside the solar system. They lie at the root of interstellar chemistry, being among the first species to form in initially-atomic gas, along with molecular hydrogen and its associated ions. Because the chemical pathways leading to the formation of interstellar hydrides are relatively simple, the analysis of the observed abundances is relatively straightforward and provides key information about the environments where hydrides are found. Recent years have seen rapid progress in our understanding of interstellar hydrides, thanks largely to far-IR and submillimeter observations performed with the Herschel Space Observatory. In this review, we will discuss observations of interstellar hydrides, along with the advanced modeling approaches that have been used to interpret them, and the unique information that has thereby been obtained.

preprint2015arXiv

Herschel observations of interstellar chloronium. II - Detections toward G29.96-0.02, W49N, W51, and W3(OH), and determinations of the ortho-to-para and $^{35}$Cl/$^{37}$Cl isotopic ratios

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

preprint2014arXiv

Herschel Far-Infrared Spectral-mapping of Orion BN/KL Outflows: Spatial distribution of excited CO, H2O, OH, O and C+ in shocked gas

We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, most of them rotationally excited lines of 12CO (up to J=48-47), H2O, OH, 13CO, and HCN. Peaks 1/2 are characterized by a very high L(CO)/L(FIR)~5x10^{-3} ratio and a plethora of far-IR H2O emission lines. The high-J CO and OH lines are a factor ~2 brighter toward Peak 1 whereas several excited H2O lines are ~50% brighter toward Peak 2. A simplified non-LTE model allowed us to constrain the dominant gas temperature components. Most of the CO column density arises from Tk~200-500 K gas that we associate with low-velocity shocks that fail to sputter grain ice mantles and show a maximum gas-phase H2O/CO~10^{-2} abundance ratio. In addition, the very excited CO (J>35) and H2O lines reveal a hotter gas component (Tk~2500 K) from faster (v_S>25 km/s) shocks that are able to sputter the frozen-out H2O and lead to high H2O/CO>~1 abundance ratios. The H2O and OH luminosities cannot be reproduced by shock models that assume high (undepleted) abundances of atomic oxygen in the preshock gas and/or neglect the presence of UV radiation in the postshock gas. Although massive outflows are a common feature in other massive star-forming cores, Orion BN/KL seems more peculiar because of its higher molecular luminosities and strong outflows caused by a recent explosive event.

preprint2014arXiv

Herschel HIFI observations of O$_2$ toward Orion: special conditions for shock enhanced emission

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

preprint2014arXiv

Revised spectroscopic parameters of SH$^+$ from ALMA and IRAM 30 m observations

Hydrides represent the first steps of interstellar chemistry. Sulfanylium (SH$^+$), in particular, is a key tracer of energetic processes. We used ALMA and the IRAM 30 m telescope to search for the lowest frequency rotational lines of SH$^+$ toward the Orion Bar, the prototypical photo-dissociation region illuminated by a strong UV radiation field. On the basis of previous $Herschel$/HIFI observations of SH$^+$, we expected to detect emission of the two SH$^+$ hyperfine structure (HFS) components of the $N_J = 1_0 - 0_1$ fine structure (FS) component near 346 GHz. While we did not observe any lines at the frequencies predicted from laboratory data, we detected two emission lines, each $\sim$15 MHz above the SH$^+$ predictions and with relative intensities and HFS splitting expected for SH$^+$. The rest frequencies of the two newly detected lines are more compatible with the remainder of the SH$^+$ laboratory data than the single line measured in the laboratory near 346 GHz and previously attributed to SH$^+$. Therefore, we assign these new features to the two SH$^+$ HFS components of the $N_J = 1_0 - 0_1$ FS component and re-determine its spectroscopic parameters, which will be useful for future observations of SH$^+$, in particular if its lowest frequency FS components are studied. Our observations demonstrate the suitability of these lines for SH$^+$ searches at frequencies easily accessible from the ground.

preprint2013arXiv

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

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

preprint2013arXiv

Probing the role of polycyclic aromatic hydrocarbons in the photoelectric heating within photodissociation regions

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

preprint2012arXiv

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

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

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.

preprint2010arXiv

Interstellar OH+, H2O+ and H3O+ along the sight-line to G10.6-0.4

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

preprint2009arXiv

Far-Infrared detection of neutral atomic oxygen toward the Horsehead Nebula

We present the first detection of neutral atomic oxygen (3P_1-3P_2 fine structure line at ~63um) toward the Horsehead photodissociation region (PDR). The cloud has been mapped with the Spitzer Space Telescope at far-IR (FIR) wavelengths using MIPS in the spectral energy distribution (SED) mode. The [OI]63um line peaks at the illuminated edge of the cloud at AV~0.1-0.5 (inward the gas becomes too cold and outward the gas density drops). The luminosity carried by the [OI]63um line represents a significant fraction of the total FIR dust luminosity (I_63/I_FIR~4x10^-3). We analyze the dust continuum emission and the nonlocal OI excitation and radiative transfer in detail. The observations are reproduced with a gas density of n_H~10^4 cm^-3 and gas and dust temperatures of T_k~100 K and T_d~30 K. We conclude that the determination of the OI 3P_J level populations and emergent line intensities at such ``low'' densities is a complex non-LTE problem. FIR radiative pumping, [OI]63um subthermal emission, [OI]145um suprathermal and even maser emission can occur and decrease the resulting [OI]63/145 intensity ratio. The Herschel Space Observatory, observing from ~55 to 672um, will allow us to exploit the diagnostic power of FIR fine structure lines with unprecedented resolution and sensitivity.

preprint2009arXiv

Protoplanetary gas disks in the far infrared

The physical and chemical conditions in young protoplanetary disks set the boundary conditions for planet formation. Although the dust in disks is relatively easily detected as a far-IR photometric ``excess'' over the expected photospheric emission, much less is known about the gas phase. It seems clear that an abrupt transition from massive optically thick disks (gas-rich structures where only ~1% of the total mass is in the form of dust) to tenuous debris disks almost devoid of gas occurs at ~10^7 years, by which time the majority of at least the giant planets must have formed. Indeed, these planets are largely gaseous and thus they must assemble before the gas disk dissipates. Spectroscopic studies of the disk gas content at different evolutive stages are thus critical. Far-IR water vapor lines and atomic fine structure lines from abundant gas reservoirs (e.g., [OI]63um, [SI]56um, [SiII]34um) are robust tracers of the gas in disks. Spectrometers on board Herschel will detect some of these lines toward the closest, youngest and more massive protoplanetary disks. However, according to models, Herschel will not reach the required sensitivity to (1) detect the gas residual in more evolved and tenuous transational disks that are potentially forming planets and (2) detect the gas emission from less massive protoplanetary disks around the most numerous stars in the Galaxy (M-type and cooler dwarfs). Both are unique goals for SPICA/SAFARI. Besides, SAFARI will be able to detect the far-IR modes of water ice at ~44 and ~62um, and thus allow water ice to be observed in many protoplanetary systems and fully explore its impact on planetary formation and evolution.

preprint2009arXiv

Studies of Expolanets and Solar Systems with SPICA

The SPace Infrared telescope for Cosmology and Astrophysics (SPICA) is a proposed mid-to-far infrared (4-200 um) astronomy mission, scheduled for launch in 2017. A single, 3.5m aperture telescope would provide superior image quality at 5-200 um, and its very cold (~5 K) instrumentation would provide superior sensitivity in the 25-200 um wavelength regimes. This would provide a breakthrough opportunity for studies of exoplanets, protoplanetary and debris disk, and small solar system bodies. This paper summarizes the potential scientific impacts for the proposed instrumentation.

preprint2009arXiv

The interstellar gas seen in the mid- and far-infrared: The promise of SPICA Space Telescope

The mid- and far-IR spectral ranges are critical windows to characterize the physical and chemical processes that transform the interstellar gas and dust into stars and planets. Sources in the earliest phases of star formation and in the latest stages of stellar evolution release most of their energy at these wavelengths. Besides, the mid- and far-IR ranges provide key spectral diagnostics of the gas chemistry (water, light hydrides, organic species ...), of the prevailing physical conditions (H2, atomic fine structure lines...), and of the dust mineral and ice composition that can not be observed from ground-based telescopes. With the launch of JAXA's SPICA telescope, uninterrupted studies in the mid- and far-IR will be possible since ESA's Infrared Space Observatory (1995). In particular, SAFARI will provide full access to the 34-210um waveband through several detector arrays and flexible observing modes (from broadband photometry to medium resolution spectroscopy with R~3,000 at 63um), and reaching very high line sensitivities (~10^-19 Wm^-2, 5sigma-1hr) within a large FOV (~2'x2'). Compared to previous far-IR instruments (ISO/LWS, Akari/FIS, Spitzer/MIPS and Herschel/PACS), SAFARI will provide a superior way to obtain fully-sampled spectro-images and continuous SEDs of very faint and extended ISM sources in a wavelength domain not accessible to JWST or ALMA. The much increased sensitivity of SPICA will allow us to step forward and reveal not only the chemical complexity in the local ISM, but also in the extragalactic ISM routinely.

preprint2008arXiv

The far-IR view of Sgr B2 and Orion KL

We summarize the main highlights from ISO observations towards Sgr B2 and Orion KL in the far-IR domain (~43 to 197 um). Both Star-Forming Regions are among the best sources to construct a template for more distant and unresolved regions (e.g., extragalactic). We stress some peculiarities in the interpretation (excitation and radiative transfer) of far-IR spectral lines and dust continuum emission.

preprint2007arXiv

Deuterium fractionation in the Horsehead edge

Deuterium fractionation is known to enhance the [DCO+]/[HCO+] abundance ratio over the D/H elemental ratio of about 1e-5 in the cold and dense gas typically found in pre-stellar cores. We report the first detection and mapping of very bright DCO+ J=3-2 and J=2-1 lines (3 and 4 K respectively) towards the Horsehead photodissociation region (PDR) observed with the IRAM-30m telescope. The DCO+ emission peaks close to the illuminated warm edge of the nebula (< 50" or about 0.1 pc away). Detailed nonlocal, non-LTE excitation and radiative transfer analyses have been used to determine the prevailing physical conditions and to estimate the DCO+ and H13CO+ abundances from their line intensities. A large [DCO+]/[HCO+] abundance ratio (>= 0.02) is inferred at the DCO+ emission peak, a condensation shielded from the illuminating far-UV radiation field where the gas must be cold (10-20 K) and dense (>= 2x10^5 cm-3). DCO+ is not detected in the warmer photodissociation front, implying a lower [DCO+]/[HCO+] ratio (< 1e-3). According to our gas phase chemical predictions, such a high deuterium fractionation of HCO+ can only be explained if the gas temperature is below 20 K, in good agreement with DCO+ excitation calculations.

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.