Source author record

D. C. Lis

D. C. Lis 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

76works
5topics
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

76 published item(s)

preprint2026arXiv

A Multi-Wavelength Study of Comet C/2022 E3 (ZTF): Complementary ALMA and JWST Investigations of Water and Methanol in Cometary Comae

Long-period comets, which are often considered to be representative of material in the protoplanetary disk that formed the Solar System, are ideal to investigate the question of chemical inheritance in astronomy. Determining the chemistry of comets, both individually and as a population, has become of great importance in comparative studies against sources representative of evolutionary precursors to planetary systems. Contemporaneous observations of long-period comet C/2022 E3 (ZTF) were obtained with the JWST and the Atacama Large Millimeter/submillimeter Array (ALMA) in early 2023 March. This work focuses on \ce{CH3OH} measurements from both ALMA and JWST as well as \ce{H2O} measurements from JWST. Radiative transfer modeling of \ce{CH3OH} and \ce{H2O} was performed to investigate spatial variations in rotational temperature, column density, and production rates, as well as a comparison of derived values between the two telescopes. Most of the spatial distributions of the modeled values are centrally peaked, and the modeled values from JWST are all within the error bars of the average values from ALMA. C/2022 E3 (ZTF) also displays an enhancement in modeled rotational temperature in the anti-Sunward direction that is shown to be statistically significant. Based on non-LTE radiative transfer modeling, the declining \ce{H2O} rotational temperatures as a function of nucleocentric distance observed by JWST can be explained primarily as a result of rotational line cooling. The values derived in this work are in general agreement with single-dish millimeter-wave observations.

preprint2025arXiv

ALMAGAL V. Relations between the core populations and the parent clump physical properties

Context. The fragmentation of massive molecular clumps into smaller, potentially star-forming cores plays a key role in the processes of high-mass star formation. The ALMAGAL project offers high-resolution data to investigate these processes across various evolutionary stages in the Galactic plane. Aims. This study aims at correlating the fragmentation properties of massive clumps, obtained from ALMA observations, with their global physical parameters (e.g., mass, surface density, and temperature) and evolutionary indicators (such as luminosity-to-mass ratio and bolometric temperature) obtained from Herschel observations. It seeks to assess whether the cores evolve in number and mass in tandem with their host clumps, and to determine the possible factors influencing the formation of massive cores (M > 24M_\odot). Methods. We analyzed the masses of 6348 fragments, estimated from 1.4 mm continuum data for 1007 ALMAGAL clumps. Leveraging this unprecedentedly large data set, we evaluated statistical relationships between clump parameters, estimated over about 0.1 pc scales, and fragment properties, corresponding to scales of a few 1000 au, while accounting for potential biases related to distance and observational resolution. Our results were further compared with predictions from numerical simulations. Results. The fragmentation level correlates preferentially with clump surface density, supporting a scenario of density-driven fragmentation, whereas it does not show any clear dependence on total clump mass. Both the mass of the most massive core and the core formation efficiency show a broad range and increase on average by an order of magnitude in the intervals spanned by evolutionary indicators such as clump dust temperature and the luminosity-to-mass ratio. This suggests that core growth continues throughout the clump evolution, favoring clump-fed over core-fed theoretical scenarios.

preprint2022arXiv

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies.

preprint2022arXiv

Water, hydrogen cyanide, carbon monoxide, and dust production from distant comet 29P/Schwassmann-Wachmann 1

29P/Schwassmann-Wachmann 1 is a distant Centaur/comet, showing persistent CO-driven activity and frequent outbursts. We used the Herschel space observatory in 2010, 2011, and 2013 to observe H$_2$O and NH$_3$ and to image the dust coma. Observations with the IRAM 30 m were undertaken in 2007, 2010, 2011, and 2021 to monitor the CO production rate and to search for HCN. Modeling was performed to constrain the size of the sublimating icy grains and to derive the dust production rate. HCN is detected for the first time in comet 29P (at 5$σ$ in the line area). H$_2$O is detected as well, but not NH$_3$. H$_2$O and HCN line shapes differ strongly from the CO line shape, indicating that these two species are released from icy grains. CO production rates are in the range (2.9-5.6) $\times$ 10$^{28}$ s$^{-1}$ (1400--2600 kg s$^{-1}$). A correlation between the CO production rate and coma brightness is observed, as is a correlation between CO and H$_2$O production. The correlation obtained between the excess of CO production and excess of dust brightness with respect to the quiescent state is similar to that established for the continuous activity of comet Hale-Bopp. The measured $Q$(H$_2$O)/$Q$(CO) and $Q$(HCN)/$Q$(CO) production rate ratios are 10.0 $\pm$ 1.5 % and 0.12 $\pm$ 0.03 %, respectively, averaging the April-May 2010 measurements ($Q$(H$_2$O) = (4.1 $\pm$ 0.6) $\times$ 10$^{27}$ s$^{-1}$, $Q$(HCN) = (4.8 $\pm$ 1.1) $\times$ 10$^{25}$ s$^{-1}$). We derive three independent and similar values of the effective radius of the nucleus, $\sim$ 31 $\pm$ 3 km. The inferred dust mass-loss rates during quiescent phases are in the range 30-120 kg s$^{-1}$, indicating a dust-to-gas mass ratio $<$ 0.1 during quiescent activity. We conclude that strong local heterogeneities exist on the surface of 29P, with quenched dust activity from most of the surface, but not in outbursting regions.

preprint2016arXiv

Collisional excitation of doubly and triply deuterated ammonia ND$_2$H and ND$_3$ by H$_2$

The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios.

preprint2016arXiv

Isotopic ratios of H, C, N, O, and S in comets C/2012 F6 (Lemmon) and C/2014 Q2 (Lovejoy)

The apparition of bright comets C/2012 F6 (Lemmon) and C/2014 Q2 (Lovejoy) in March-April 2013 and January 2015, combined with the improved observational capabilities of submillimeter facilities, offered an opportunity to carry out sensitive compositional and isotopic studies of the volatiles in their coma. We observed comet Lovejoy with the IRAM 30m telescope between 13 and 26 January 2015, and with the Odin submillimeter space observatory on 29 January - 3 February 2015. We detected 22 molecules and several isotopologues. The H$_2^{16}$O and H$_2^{18}$O production rates measured with Odin follow a periodic pattern with a period of 0.94 days and an amplitude of ~25%. The inferred isotope ratios in comet Lovejoy are $^{16}$O/$^{18}$O = 499 $\pm$ 24 and D/H = 1.4 $\pm$ 0.4 $\times 10^{-4}$ in water, $^{32}$S/$^{34}$S = 24.7 $\pm$ 3.5 in CS, all compatible with terrestrial values. The ratio $^{12}$C/$^{13}$C = 109 $\pm$ 14 in HCN is marginally higher than terrestrial and $^{14}$N/$^{15}$N = 145 $\pm$ 12 in HCN is half the Earth ratio. Several upper limits for D/H or 12C/13C in other molecules are reported. From our observation of HDO in comet C/2014 Q2 (Lovejoy), we report the first D/H ratio in an Oort Cloud comet that is not larger than the terrestrial value. On the other hand, the observation of the same HDO line in the other Oort-cloud comet, C/2012 F6 (Lemmon), suggests a D/H value four times higher. Given the previous measurements of D/H in cometary water, this illustrates that a diversity in the D/H ratio and in the chemical composition, is present even within the same dynamical group of comets, suggesting that current dynamical groups contain comets formed at very different places or times in the early solar system.

preprint2016arXiv

Measuring the Distribution and Excitation of Cometary CH3OH Using ALMA

The Atacama Large Millimeter/submillimeter Array (ALMA) was used to obtain measurements of spatially and spectrally resolved CH$_3$OH emission from comet C/2012 K1 (PanSTARRS) on 28-29 June 2014. Detection of 12-14 emission lines of CH$_3$OH on each day permitted the derivation of spatially-resolved rotational temperature profiles (averaged along the line of sight), for the innermost 5000 km of the coma. On each day, the CH$_3$OH distribution was centrally peaked and approximately consistent with spherically symmetric, uniform outflow. The azimuthally-averaged CH$_3$OH rotational temperature ($T_{rot}$) as a function of sky-projected nucleocentric distance ($ρ$), fell by about 40 K between $ρ=0$ and 2500 km on 28 June, whereas on 29 June, $T_{rot}$ fell by about 50 K between $ρ=$0 km and 1500 km. A remarkable ($\sim50$ K) rise in $T_{rot}$ at $ρ=$1500-2500 km on 29 June was not present on 28 June. The observed variations in CH$_3$OH rotational temperature are interpreted primarily as a result of variations in the coma kinetic temperature due to adiabatic cooling, and heating through Solar irradiation, but collisional and radiative non-LTE excitation processes also play a role.

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.

preprint2016arXiv

Star Formation and Feedback: A Molecular Outflow-Prestellar Core Interaction in L1689N

We present Herschel, ALMA Compact Array (ACA), and Caltech Submillimeter Observatory (CSO) observations of the prestellar core in L1689N, which has been suggested to be interacting with a molecular outflow driven by the nearby solar type protostar IRAS 16293-2422. This source is characterized by some of the highest deuteration levels seen in the interstellar medium. The change in the NH2D line velocity and width across the core provides clear evidence of an interaction with the outflow, traced by the high-velocity water emission. Quiescent, cold gas, characterized by narrow line widths is seen in the NE part of the core, while broader, more disturbed line profiles are seen in the W/SW part. Strong N2D+ and ND3 emission is detected with the ACA, extending S/SW from the peak of the single-dish NH2D emission. The ACA data also reveal the presence a compact dust continuum source, with a mean size of ~1100 au, a central density of (1-2) 10^7 cm-3, and a mass of 0.2-0.4 Msun. The dust emission peak is displaced ~5" to the south with respect to the N2D+ and ND3 emission, as well as the single-dish dust continuum peak, suggesting that the northern, quiescent part of the core is characterized by spatially extended continuum emission, which is resolved out by the interferometer. We see no clear evidence of fragmentation in this quiescent part of the core, which could lead to a second generation of star formation, although a weak dust continuum source is detected in this region in the ACA data.

preprint2016arXiv

The physical and chemical structure of Sagittarius B2, I. Three-dimensional thermal dust and free-free continuum modeling on 100 au to 45 pc scales

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

preprint2016arXiv

WISH VI. Constraints on UV and X-ray irradiation from a survey of hydrides in low- to high-mass YSOs

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

preprint2015arXiv

Herschel HIFI Observations of the Sgr A +50 km/s Cloud. Deep Searches for O2 in Emission and Foreground Absorption

To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).

preprint2015arXiv

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

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

preprint2014arXiv

[CII] absorption and emission in the diffuse interstellar medium across the Galactic Plane

Ionized carbon is the main gas-phase reservoir of carbon in the neutral diffuse interstellar medium and its 158 micron fine structure transition [CII] is the most important cooling line of the diffuse interstellar medium (ISM). We combine [CII] absorption and emission spectroscopy to gain an improved understanding of physical conditions in the different phases of the ISM. We present high resolution [CII] spectra obtained with the Herschel/HIFI instrument towards bright dust continuum sources regions in the Galactic plane, probing simultaneously the diffuse gas along the line of sight and the background high-mass star forming regions. These data are complemented by observations of the 492 and 809 GHz fine structure lines of atomic carbon and by medium spectral resolution spectral maps of the fine structure lines of atomic oxygen at 63 and 145 microns with Herschel/PACS. We show that the presence of foreground absorption may completely cancel the emission from the background source in medium spectral resolution data and that high spectral resolution spectra are needed to interpret the [CII] and [OI] emission and the [CII]/FIR ratio. This phenomenon may explain part of the [CII]/FIR deficit seen in external luminous infrared galaxies. The C+ and C excitation in the diffuse gas is consistent with a median pressure of 5900 Kcm-3 for a mean TK ~100 K. The knowledge of the gas density allows us to determine the filling factor of the absorbing gas along the selected lines of sight: the median value is 2.4 %, in good agreement with the CNM properties. The mean excitation temperature is used to derive the average cooling due to C+ in the Galactic plane : 9.5 x 10^{-26} erg/s/H. Along the observed lines of sight, the gas phase carbon abundance does not exhibit a strong gradient as a function of Galacto-centric radius and has a weighted average of C/H = 1.5 +/- 0.4 x 10^{-4}.

preprint2014arXiv

ALMA measurements of the HNC and HC$_3$N distributions in Titan's atmosphere

We present spectrally and spatially-resolved maps of HNC and HC$_3$N emission from Titan's atmosphere, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 2013 November 17. These maps show anisotropic spatial distributions for both molecules, with resolved emission peaks in Titan's northern and southern hemispheres. The HC$_3$N maps indicate enhanced concentrations of this molecule over the poles, consistent with previous studies of Titan's photochemistry and atmospheric circulation. Differences between the spectrally-integrated flux distributions of HNC and HC$_3$N show that these species are not co-spatial. The observed spectral line shapes are consistent with HNC being concentrated predominantly in the mesosphere and above (at altitudes $z\gtrsim 400$ km), whereas HC$_3$N is abundant at a broader range of altitudes ($z\approx70$-600 km). From spatial variations in the HC$_3$N line profile, the locations of the HC$_3$N emission peaks are shown to be variable as a function of altitude. The peaks in the integrated emission from HNC and the line core (upper-atmosphere) component of HC$_3$N (at $z\gtrsim300$ km) are found to be asymmetric with respect to Titan's polar axis, indicating that the mesosphere may be more longitudinally-variable than previously thought. The spatially-integrated HNC and HC$_3$N spectra are modeled using the NEMESIS planetary atmosphere code and the resulting best-fitting disk-averaged vertical mixing ratio (VMR) profiles are found to be in reasonable agreement with previous measurements for these species. Vertical column densities of the best-fitting gradient models for HNC and HC$_3$N are $1.9\times10^{13}$ cm$^{-2}$ and $2.3\times10^{14}$ cm$^{-2}$, respectively.

preprint2014arXiv

Collisional excitation of singly deuterated ammonia NH$_2$D by H$_2$

The availability of collisional rate coefficients with H$_2$ is a pre-requisite for interpretation of observations of molecules whose energy levels are populated under non local thermodynamical equilibrium conditions. In the current study, we present collisional rate coefficients for the NH$_2$D / para--H$_2$($J_2 = 0,2$) collisional system, for energy levels up to $J_τ= 7_7$ ($E_u$$\sim$735 K) and for gas temperatures in the range $T = 5-300$K. The cross sections are obtained using the essentially exact close--coupling (CC) formalism at low energy and at the highest energies, we used the coupled--states (CS) approximation. For the energy levels up to $J_τ= 4_2$ ($E_u$$\sim$215 K), the cross sections obtained through the CS formalism are scaled according to a few CC reference points. These reference points are subsequently used to estimate the accuracy of the rate coefficients for higher levels, which is mainly limited by the use of the CS formalism. Considering the current potential energy surface, the rate coefficients are thus expected to be accurate to within 5\% for the levels below $J_τ= 4_2$, while we estimate an accuracy of 30\% for higher levels.

preprint2014arXiv

Complex organic molecules in comets C/2012 F6 (Lemmon) and C/2013 R1 (Lovejoy): detection of ethylene glycol and formamide

A spectral survey in the 1 mm wavelength range was undertaken in the long-period comets C/2012 F6 (Lemmon) and C/2013 R1 (Lovejoy) using the 30 m telescope of the Institut de radioastronomie millimétrique (IRAM) in April and November-December 2013. We report the detection of ethylene glycol (CH$_2$OH)$_2$ (aGg' conformer) and formamide (NH$_2$CHO) in the two comets. The abundances relative to water of ethylene glycol and formamide are 0.2-0.3% and 0.02% in the two comets, similar to the values measured in comet C/1995 O1 (Hale-Bopp). We also report the detection of HCOOH and CH$_3$CHO in comet C/2013 R1 (Lovejoy), and a search for other complex species (methyl formate, glycolaldehyde).

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 observations of gas and dust in comet C/2006 W3 (Christensen) at 5 AU from the Sun

We aimed to measure the H2O and dust production rates in C/2006 W3 (Christensen) with the Herschel Space Observatory at a heliocentric distance of ~ 5 AU. We have searched for emission in the H2O and NH3 ground-state rotational transitions at 557 GHz and 572 GHz, simultaneously, with HIFI onboard Herschel on UT 1.5 September 2010. Photometric observations of the dust coma in the 70 and 160 μm channels were acquired with the PACS instrument on UT 26.5 August 2010. A tentative 4-σ H2O line emission feature was found in the spectra obtained with the HIFI wide-band and high-resolution spectrometers, from which we derive a water production rate of $2.0(5) \times 10^{27}$ molec. s$^{-1}$. A 3-σ upper limit for the ammonia production rate of <$1.5 \times 10^{27}$ molec. s$^{-1}$ is obtained taking into account the contribution from all hyperfine components. The blueshift of the water line detected by HIFI suggests preferential emission from the subsolar point. However, it is also possible that water sublimation occurs in small ice-bearing grains that are emitted from an active region on the nucleus surface at a speed of ~ 0.2 km s$^{-1}$. The dust thermal emission was detected in the 70 and 160 μm filters, with a more extended emission in the blue channel. The dust production rates, obtained for a dust size distribution index that explains the fluxes at the photocenters of the PACS images, lie in the range 70-110 kg s$^{-1}$. Scaling the CO production rate measured post-perihelion at 3.20 and 3.32 AU, these values correspond to a dust-to-gas production rate ratio in the range 0.3-0.4. The dust production rates derived in August 2010 are roughly one order of magnitude lower than in September 2009, suggesting that the dust-to-gas production rate ratio remained approximately constant during the period when the activity became increasingly dominated by CO outgassing.

preprint2014arXiv

Herschel Survey of Galactic OH+, H2O+, and H3O+: Probing the Molecular Hydrogen Fraction and Cosmic-Ray Ionization Rate

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

preprint2014arXiv

Hydrogen Fluoride toward Luminous Nearby Galaxies: NGC 253 and NGC 4945

We present the detection of hydrogen fluoride, HF, in two luminous nearby galaxies NGC 253 and NGC 4945 using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory. The HF line toward NGC 253 has a P-Cygni profile, while an asymmetric absorption profile is seen toward NGC 4945. The P-Cygni profile in NGC 253 suggests an outflow of molecular gas with a mass of M(H$_2$)$_{out}$ $\sim$ 1 $\times$ 10$^7$ M$_\odot$ and an outflow rate as large as M $\sim$ 6.4 M$_\odot$ yr$^{-1}$. In the case of NGC 4945, the axisymmetric velocity components in the HF line profile is compatible with the interpretation of a fast-rotating nuclear ring surrounding the nucleus and the presence of inflowing gas. The gas falls into the nucleus with an inflow rate of $\le$ 1.2 M$_\odot$ yr$^{-1}$, inside a inner radius of $\le$ 200 pc. The gas accretion rate to the central AGN is much smaller, suggesting that the inflow can be triggering a nuclear starburst. From these results, the HF $J = 1-0$ line is seen to provide an important probe of the kinematics of absorbing material along the sight-line to nearby galaxies with bright dust continuum and a promising new tracer of molecular gas in high-redshift galaxies.

preprint2014arXiv

Ionization toward the high-mass star-forming region NGC 6334 I

Context. Ionization plays a central role in the gas-phase chemistry of molecular clouds. Since ions are coupled with magnetic fields, which can in turn counteract gravitational collapse, it is of paramount importance to measure their abundance in star-forming regions. Aims. We use spectral line observations of the high-mass star-forming region NGC 6334 I to derive the abundance of two of the most abundant molecular ions, HCO+ and N2H+, and consequently, the cosmic ray ionization rate. In addition, the line profiles provide information about the kinematics of this region. Methods. We present high-resolution spectral line observations conducted with the HIFI instrument on board the Herschel Space Observatory of the rotational transitions with Jup > 5 of the molecular species C17O, C18O, HCO+, H13CO+, and N2H+. Results. The HCO+ and N2H+ line profiles display a redshifted asymmetry consistent with a region of expanding gas. We identify two emission components in the spectra, each with a different excitation, associated with the envelope of NGC 6334 I. The physical parameters obtained for the envelope are in agreement with previous models of the radial structure of NGC 6334 I based on submillimeter continuum observations. Based on our new Herschel/HIFI observations, combined with the predictions from a chemical model, we derive a cosmic ray ionization rate that is an order of magnitude higher than the canonical value of 10^(-17) s-1. Conclusions. We find evidence of an expansion of the envelope surrounding the hot core of NGC 6334 I, which is mainly driven by thermal pressure from the hot ionized gas in the region. The ionization rate seems to be dominated by cosmic rays originating from outside the source, although X-ray emission from the NGC 6334 I core could contribute to the ionization in the inner part of the envelope.

preprint2014arXiv

Mapping the release of volatiles in the inner comae of comets C/2012 F6 (Lemmon) and C/2012 S1 (ISON) using the Atacama Large Millimeter/Submillimeter Array

Results are presented from the first cometary observations using the Atacama Large Millimeter/Submillimeter Array (ALMA), including measurements of the spatially-resolved distributions of HCN, HNC, H$_2$CO and dust within the comae of two comets: C/2012 F6 (Lemmon) and C/2012 S1 (ISON), observed at heliocentric distances of 1.5 AU and 0.54 AU, respectively. These observations (with angular resolution $\approx0.5''$), reveal an unprecedented level of detail in the distributions of these fundamental cometary molecules, and demonstrate the power of ALMA for quantitative measurements of the distributions of molecules and dust in the inner comae of typical bright comets. In both comets, HCN is found to originate from (or within a few hundred km of) the nucleus, with a spatial distribution largely consistent with spherically-symmetric, uniform outflow. By contrast, the HNC distributions are clumpy and asymmetrical, with peaks at cometocentric radii $\sim$500-1000~km, consistent with release of HNC in collimated outflow(s). Compared to HCN, the H$_2$CO distribution in comet Lemmon is very extended. The interferometric visibility amplitudes are consistent with coma production of H$_2$CO and HNC from unidentified precursor material(s) in both comets. Adopting a Haser model, the H$_2$CO parent scale-length is found to be a few thousand km in Lemmon and only a few hundred km in ISON, consistent with destruction of the precursor by photolysis or thermal degradation at a rate which scales in proportion to the Solar radiation flux.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Widespread Rotationally-Hot Hydronium Ion in the Galactic Interstellar Medium

We present new observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ~ 500 K, rotational temperatures characterizing the population of the highly-excited metastable H3O+ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone, but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ~ 380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic ray fluxes, shocks and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly-excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation to the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process, not limited to the active environments associated with galactic nuclei.

preprint2013arXiv

A Herschel Study of D/H in Water in the Jupiter-Family Comet 45P/Honda-Mrkos-Pajdusakova and Prospects for D/H Measurements with CCAT

We present Herschel observations of water isotopologues in the atmosphere of the Jupiter-family comet 45P/Honda-Mrkos-Pajdusakova. No HDO emission is detected, with a 3 sigma upper limit of 2.0 10-4 for the D/H ratio. This value is consistent with the earlier Herschel measurement in the Jupiter-family comet 103P/Hartley 2. The canonical value of 3 10-4 measured pre-Herschel in a sample of Oort-cloud comets can be excluded at a 4.5 sigma level. The observations presented here further confirm that a diversity of D/H ratios exists in the comet population and emphasize the need for additional measurements with future ground-based facilities, such as CCAT, in the post-Herschel era.

preprint2013arXiv

ALMA Observations of the Galactic Center: SiO Outflows and High Mass Star Formation near Sgr A*

ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FWZI) $\sim110-147$ \kms. The remaining clumps, distributed mainly to the NE of the ionized mini-spiral, have narrow FWZI ($\sim18-56$ \kms). Using CARMA SiO (2-1) data, LVG modeling of the the SiO line ratios for the broad velocity clumps, constrains the column density N(SiO) $\sim10^{14}$ cm$^{-2}$, and the H$_2$ gas density n$_{\rm H_2}=(3-9)\times10^5$ cm$^{-3}$ for an assumed kinetic temperature 100-200K. The SiO clumps are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last $10^4-10^5$ years. Support for this interpretation is provided by the SiO (5-4) line luminosities and velocity widths which lie in the range measured for protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources shows two YSO candidates near SiO clumps, supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhances the gas density, before the gas cloud become gravitationally unstable near Sgr A*. Alternatively, collisions between clumps in the ring may trigger gravitational collapse.

preprint2013arXiv

CH2D+, the Search for the Holy Grail

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

preprint2013arXiv

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

Gas and dust productions of comet 103P/Hartley 2 from millimetre observations: interpreting rotation-induced time variations

Comet 103P/Hartley 2 made a close approach to the Earth in October 2010. It was the target of an extensive observing campaign and was visited by the Deep Impact spacecraft (mission EPOXI). We present observations of HCN and CH3OH emission lines conducted with the IRAM Plateau de Bure interferometer on 22-23, 28 October and 4, 5 November 2010 at 1.1, 1.9 and 3.4 mm wavelengths. The thermal emission from the dust coma and nucleus is detected simultaneously. Interferometric images with unprecedented spatial resolution are obtained. A sine-wave variation of the thermal continuum is observed in the 23 October data, that we associate with the nucleus thermal light curve. The nucleus contributes up to 30-55 % of the observed continuum. The large dust-to-gas ratio (in the range 2-6) can be explained by the unusual activity of the comet for its size, which allows decimeter size particles and large boulders to be entrained by the gas. The rotational temperature of CH3OH is measured. We attribute the increase from 35 to 46 K with increasing beam size (from 150 to 1500 km) to radiative processes. The HCN production rate displays strong rotation-induced variations. The HCN production curve, as well as those of CO2 and H2O measured by EPOXI, are interpreted with a geometric model which takes into account the rotation and the shape of the comet. The HCN and H2O production curves are in phase, showing common sources. The 1.7h delay, in average, of HCN and H2O with respect to the CO2 production curve suggests that HCN and H2O are mainly produced by subliming icy grains. The scale length of production of HCN is determined to be on the order of 500-1000 km, implying a mean velocity of 100-200 m/s for the icy grains producing HCN. The modulation of the CO2 prouction and of the velocity offset of the HCN lines are interpreted in terms of localized sources of gas on the nucleus surface.

preprint2013arXiv

Hydrogen chloride in diffuse interstellar clouds along the line of sight to W31C (G10.6-0.4)

We report the detection of hydrogen chloride, HCl, in diffuse molecular clouds on the line of sight towards the star-forming region W31C (G10.6-0.4). The J = 1-0 lines of the two stable HCl isotopologues, H35Cl and H37Cl, are observed using the 1b receiver of the Heterodyne Instrument for the Far-Infrared (HIFI) aboard the Herschel Space Observatory. The HCl line is detected in absorption, over a wide range of velocities associated with diffuse clouds along the line of sight to W31C. The analysis of the absorption strength yields a total HCl column density of few 10^13 cm^-2, implying that HCl accounts for ~0.6 % of the total gasphase chlorine, which exceeds by a factor of ~6 the theoretical model predictions. This result is comparable to those obtained from the chemically-related species H2Cl+ and HCl+, for which large column densities have also been reported on the same line of sight. The source of discrepancy between models and observations is still unknown; however, the detection of these Cl-bearing molecules, provides key constraints for the chlorine chemistry in the diffuse 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.

preprint2012arXiv

A multi-wavelength view of the Galactic center dust ridge reveals little star formation

The Galactic center dust ridge consists of a narrow string of massive condensations identified in submillimeter dust continuum emission. To determine whether new high-mass stars are forming in this region, we performed new observations at 870 $μ$m with the Atacama Pathfinder Experiment telescope and at 8.4 GHz with the Very Large Array. We complement our data with recent maser and mid-infrared results. The ridge's clouds are dark at mid-infrared wavelengths, indicating the presence of cold, high column density material. In combination with existing temperature measurements in the dust ridge, we determined the masses of the largest clouds. The results show that the dust ridge contains a very massive reservoir of molecular material. We find five radio sources at 8.4 GHz in the general dust ridge vicinity but outside of the dust ridge clouds, which are probably all excited by massive young stars, whose properties we constrain. Our observations exclude the existence of zero age main sequence stars with spectral types earlier than B0.5 within the dust ridge clouds. The only indication of ongoing high-mass star formation inside the clouds are class II methanol masers that are found in two of the clouds. Except for a weak water maser, found in previous observations, no signs of star formation are detected in the massive cloud M0.25+0.012.

preprint2012arXiv

Ammonia and other parent molecules in comet 10P/Tempel 2 from Herschel/HIFI and ground-based radio observations

The Jupiter-family comet 10P/Tempel 2 was observed during its 2010 return with the Herschel Space Observatory. We present here the observation of the (J, K) = (1, 0)-(0, 0) transition of ammonia at 572 GHz in this comet with the Heterodyne Instrument for the Far Infrared (HIFI) of Herschel. We also report on radio observations of other molecules (HCN, CH3OH, H2S and CS) obtained during the 1999 return of the comet with the CSO telescope and the JCMT, and during its 2010 return with the IRAM 30-m telescope. Molecular abundances relative to water are 0.09%, 1.8%, 0.4%, and 0.08% for HCN, CH3OH, H2S, and CS, respectively. An abundance of 0.5% for NH3 is obtained, which is similar to the values measured in other comets. The hyperfine structure of the ammonia line is resolved for the first time in an astronomical source. Strong anisotropy in the outgassing is present in all observations from 1999 to 2010 and is modelled to derive the production rates.

preprint2012arXiv

An upper limit for the water outgassing rate of the main-belt comet 176P/LINEAR observed with Herschel/HIFI

176P/LINEAR is a member of the new cometary class known as main-belt comets (MBCs). It displayed cometary activity shortly during its 2005 perihelion passage that may be driven by the sublimation of sub-surface ices. We have therefore searched for emission of the H2O 110-101 ground state rotational line at 557 GHz toward 176P/LINEAR with the Heterodyne Instrument for the Far Infrared (HIFI) on board the Herschel Space Observatory on UT 8.78 August 2011, about 40 days after its most recent perihelion passage, when the object was at a heliocentric distance of 2.58 AU. No H2O line emission was detected in our observations, from which we derive sensitive 3-sigma upper limits for the water production rate and column density of < 4e25 molec/s and of < 3e10 cm^{-2}, respectively. From the peak brightness measured during the object's active period in 2005, this upper limit is lower than predicted by the relation between production rates and visual magnitudes observed for a sample of comets by Jorda et al. (2008) at this heliocentric distance. Thus, 176P/LINEAR was likely less active at the time of our observation than during its previous perihelion passage. The retrieved upper limit is lower than most values derived for the H2O production rate from the spectroscopic search for CN emission in MBCs.

preprint2012arXiv

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

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

preprint2012arXiv

Discovery of Hydrogen Fluoride in the Cloverleaf Quasar at z = 2.56

We report the first detection of hydrogen fluoride (HF) toward a high redshift quasar. Using the Caltech Submillimeter Observatory (CSO) we detect the HF J = 1 - 0 transition in absorption toward the Cloverleaf, a broad absorption line (BAL) quasi-stellar object (QSO) at z=2.56. The detection is statistically significant at the ~ 6 sigma level. We estimate a lower limit of 4 \times 1014 cm-2 for the HF column density and using a previous estimate of the hydrogen column density, we obtain a lower limit of 1.7 \times 10-9 for the HF abundance. This value suggests that, assuming a Galactic N(HF)/NH ratio, HF accounts for at least ~10% of the fluorine in the gas phase along the line of sight to the Cloverleaf quasar. This observation corroborates the prediction that HF should be a good probe of the molecular gas at high redshift. Measurements of the HF abundance as a function of redshift are urgently needed to better constrain the fluorine nucleosynthesis mechanism(s).

preprint2012arXiv

Herschel measurements of the D/H and 16O/18O ratios in water in the Oort-cloud comet C/2009 P1 (Garradd)

The D/H ratio in cometary water is believed to be an important indicator of the conditions under which icy planetesimals formed and can provide clues to the contribution of comets to the delivery of water and other volatiles to Earth. Available measurements suggest that there is isotopic diversity in the comet population. The Herschel Space Observatory revealed an ocean-like ratio in the Jupiter-family comet 103P/Hartley 2, whereas most values measured in Oort-cloud comets are twice as high as the ocean D/H ratio. We present here a new measurement of the D/H ratio in the water of an Oort-cloud comet. HDO, H_2O, and H_2^18O lines were observed with high signal-to-noise ratio in comet C/2009 P1 (Garradd) using the Herschel HIFI instrument. Spectral maps of two water lines were obtained to constrain the water excitation. The D/H ratio derived from the measured H_2^16O and HDO production rates is 2.06+/-0.22 X 10**-4. This result shows that the D/H in the water of Oort-cloud comets is not as high as previously thought, at least for a fraction of the population, hence the paradigm of a single, archetypal D/H ratio for all Oort-cloud comets is no longer tenable. Nevertheless, the value measured in C/2009 P1 (Garradd) is significantly higher than the Earth's ocean value of 1.558 X 10**-4. The measured H_2^16O/H_2^18O ratio of 523+/-32 is, however, consistent with the terrestrial value.

preprint2012arXiv

Herschel/HIFI observations of [CII] and [13CII] in PDRs

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

preprint2012arXiv

Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the SMA

We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions. In the framework of the Herschel guaranteed time key program CHESS, NGC 6334I is investigated by using HIFI aboard the Herschel Space Observatory. A spectral line survey is carried out in the frequency range 480-1907 GHz, and auxiliary interferometric data from the SMA in the 230 GHz band provide spatial information for disentangling the different physical components contributing to the HIFI spectrum. The spectral lines are identified with the aid of former surveys and spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum with XCLASS, assuming local thermal equilibrium, and best fit parameters are derived using the model optimization package MAGIX. A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High- energy levels of the dominant emitter methanol and vibrationally excited HCN are detected. The number of unidentified lines remains low with 75, or less than 2 percent of the lines detected. The modeling suggests that several spectral features need two or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible in the SiO emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules identified in SMA1 and O-bearing molecules in SMA2. Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances provide further constraints for astrochemical models.

preprint2012arXiv

Multi-line detection of O2 toward rho Oph A

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

preprint2012arXiv

The Abundance, Ortho/Para Ratio, and Deuteration of Water in the High-Mass Star Forming Region NGC 6334 I

We present Herschel/HIFI observations of 30 transitions of water isotopologues toward the high-mass star forming region NGC 6334 I. The line profiles of H_2^{16}O, H_2^{17}O, H_2^{18}O, and HDO show a complex pattern of emission and absorption components associated with the embedded hot cores, a lower-density envelope, two outflow components, and several foreground clouds, some associated with the NGC 6334 complex, others seen in projection against the strong continuum background of the source. Our analysis reveals an H2O ortho/para ratio of 3 +/- 0.5 in the foreground clouds, as well as the outflow. The water abundance varies from ~10^{-8} in the foreground clouds and the outer envelope to ~10^{-6} in the hot core. The hot core abundance is two orders of magnitude below the chemical model predictions for dense, warm gas, but within the range of values found in other Herschel/HIFI studies of hot cores and hot corinos. This may be related to the relatively low gas and dust temperature (~100 K), or time dependent effects, resulting in a significant fraction of water molecules still locked up in dust grain mantles. The HDO/H_2O ratio in NGC 6334 I, ~2 10^{-4}, is also relatively low, but within the range found in other high-mass star forming regions.

preprint2012arXiv

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

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

preprint2012arXiv

Upper limit for the D2H+ ortho-to-para ratio in the prestellar core 16293E (CHESS)

The H3+ ion plays a key role in the chemistry of dense interstellar gas clouds where stars and planets are forming. The low temperatures and high extinctions of such clouds make direct observations of H3+ impossible, but lead to large abundances of H2D+ and D2H+, which are very useful probes of the early stages of star and planet formation. The ground-state rotational ortho-D2H+ 111-000 transition at 1476.6 GHz in the prestellar core 16293E has been searched for with the Herschel/HIFI instrument, within the CHESS (Chemical HErschel Surveys of Star forming regions) Key Program. The line has not been detected at the 21 mK km/s level (3 sigma integrated line intensity). We used the ortho-H2D+ 110-111 transition and para-D2H+ 110-101 transition detected in this source to determine an upper limit on the ortho-to-para D2H+ ratio as well as the para-D2H+/ortho-H2D+ ratio from a non-LTE analysis. The comparison between our chemical modeling and the observations suggests that the CO depletion must be high (larger than 100), with a density between 5e5 and 1e6 cm-3. Also the upper limit on the ortho-D2H+ line is consistent with a low gas temperature (~ 11 K) with a ortho-to-para ratio of 6 to 9, i.e. 2 to 3 times higher than the value estimated from the chemical modeling, making it impossible to detect this high frequency transition with the present state of the art receivers.

preprint2012arXiv

Water absorption in Galactic translucent clouds: conditions and history of the gas derived from Herschel/HIFI PRISMAS observations

We present Herschel/HIFI observations of nine transitions of \hho and \hheo towards six high-mass star-forming regions, obtained as part of the PRISMAS Key Program. Water vapor in translucent clouds is detected in absorption along every sightline. We derive the column density of \hho or \hheo for the lower energy level of each transition observed. The total water column density is about a few $10^{13} \rm{cm^{-2}}$. We find that the abundance of water relative to hydrogen nuclei is $1\times10^{-8}$ in agreement with models for oxygen chemistry with high cosmic ray ionization rates. Relative to \hh, the abundance of water is remarkably constant at $5\times10^{-8}$. The abundance of water in excited levels is at most 15%, implying that the excitation temperature $T_{ex}$ in the ground state transitions is below 10 K. The column densities derived from the two ortho ground state transitions indicates that $T_{ex}\simeq5$ K and that the density $n($\hh$)$ in the clouds is $\le10^4 \rm{cm^{-3}}$. For most clouds we derive a water ortho-to-para ratio consistent with the value of 3 expected in thermodynamic equilibrium in the high temperature limit. Two clouds with large column densities exhibit a ratio significantly below 3. This may argue that the history of water molecules includes a cold phase, either when the molecules were formed on cold grains, or when they later become at least partially thermalized with the cold gas ($\sim25$ K) in the shielded, low temperature regions of the clouds; evidently, they have not yet fully thermalized with the warmer ($\sim50$ K) translucent portions of the clouds.

preprint2011arXiv

Discovery of Water Vapor in the High-Redshift Quasar APM 08279+5255 at Z=3.91

We report a detection of the excited 220-211 rotational transition of para-H2O in APM 08279+5255 using the IRAM Plateau de Bure interferometer. At z = 3.91, this is the highest-redshift detection of interstellar water to date. From LVG modeling, we conclude that this transition is predominantly radiatively pumped and on its own does not provide a good estimate of the water abundance. However, additional water transitions are predicted to be detectable in this source, which would lead to an improved excitation model. We also present a sensitive upper limit for the HF J = 1 - 0 absorption toward APM 08279+5255. While the face-on geometry of this source is not favorable for absorption studies, the lack of HF absorption is still puzzling and may be indicative of a lower fluorine abundance at z = 3.91 compared with the Galactic ISM.

preprint2011arXiv

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

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

preprint2011arXiv

Herschel/HIFI Observations of Hydrogen Fluoride Toward Sagittarius B2(M)

Herschel/HIFI observations have revealed the presence of widespread absorption by hydrogen fluoride (HF) J = 1-0 rotational transition, toward a number of Galactic sources. We present observations of HF J = 1-0 toward the high-mass star-forming region Sagittarius B2(M). The spectrum obtained shows a complex pattern of absorption, with numerous features covering a wide range of local standard of rest velocities (-130 to 100 km -1). An analysis of this absorption yields HF abundances relative to H2 of ~1.3 {\times}10-8, in most velocity intervals. This result is in good agreement with estimates from chemical models, which predict that HF should be the main reservoir of gas-phase fluorine under a wide variety of interstellar conditions. Interestingly, we also find velocity intervals in which the HF spectrum shows strong absorption features that are not present, or are very weak, in spectra of other molecules, such as 13CO (1-0) and CS (2-1). HF absorption reveals components of diffuse clouds with small extinction that can be studied for the first time. Another interesting observation is that water is significantly more abundant than hydrogen fluoride over a wide range of velocities toward Sagittarius B2(M), in contrast to the remarkably constant H2O/HF abundance ratio with average value close to unity measured toward other Galactic sources.

preprint2011arXiv

Submillimeter continuum observations of Sagittarius B2 at subarcsecond spatial resolution

We report the first high spatial resolution submillimeter continuum observations of the Sagittarius B2 cloud complex using the Submillimeter Array (SMA). With the subarcsecond resolution provided by the SMA, the two massive star-forming clumps Sgr B2(N) and Sgr B2(M) are resolved into multiple compact sources. In total, twelve submillimeter cores are identified in the Sgr B2(M) region, while only two components are observed in the Sgr B2(N) clump. The gas mass and column density are estimated from the dust continuum emission. We find that most of the cores have gas masses in excess of 100 M$_{\odot}$ and column densities above 10$^{25}$ cm$^{-2}$. The very fragmented appearance of Sgr B2(M), in contrast to the monolithic structure of Sgr B2 (N), suggests that the former is more evolved. The density profile of the Sgr B2(N)-SMA1 core is well fitted by a Plummer density distribution. This would lead one to believe that in the evolutionary sequence of the Sgr B2 cloud complex, a massive star forms first in an homogeneous core, and the rest of the cluster forms subsequently in the then fragmenting structure.

preprint2011arXiv

The Influence of Deuteration and Turbulent Diffusion on the Observed D/H Ratio

The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mixing acts to reduce the efficiency of deuteration in these clouds by increasing the ionization fraction and reducing freeze-out of heavy molecules. This leads to lower abundances for many deuterated species. We also examine the influence of turbulent mixing on the transition from atomic hydrogen to H2 and from atomic deuterium to HD near the cloud edge. We find that including turbulent diffusion in our models serves to push these transitions deeper into the cloud and helps maintain a higher atomic fraction throughout the cloud envelope. Based on these findings, we propose a new process to account for the significant scatter in the observed atomic D/H ratio for galactic sightlines extending beyond the Local Bubble. Although several mechanisms have been put forward to explain this scatter, they are unable to fully account for the range in D/H values. We suggest a scenario in which turbulent mixing of atomic and molecular gas at the edges of molecular clouds causes the observed atomic D/H ratio to vary by a factor of ~2.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

First results of Herschel/PACS observations of Neptune

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

preprint2010arXiv

First results on Martian carbon monoxide from Herschel/HIFI observations

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

preprint2010arXiv

Herschel observations of deuterated water towards Sgr B2(M)

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

Herschel/HIFI observations of interstellar OH+ and H2O+ towards W49N: a probe of diffuse clouds with a small molecular fraction

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

preprint2010arXiv

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

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

preprint2010arXiv

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

preprint2010arXiv

HIFI Observations of Water in the Atmosphere of Comet C/2008 Q3 (Garradd)

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

preprint2010arXiv

Nitrogen Isotopic Fractionation in Interstellar Ammonia

Using the Green Bank Telescope (GBT), we have obtained accurate measurements of the $^{14}$N/$^{15}$N isotopic ratio in ammonia in two nearby cold, dense molecular clouds, Barnard~1 and NGC 1333. The $^{14}$N/$^{15}$N ratio in Barnard~1, $334 \pm 50$ (3$σ$), is particularly well constrained and falls in between the local interstellar medium/proto-solar value of $\sim 450$ and the terrestrial atmospheric value of 272. The NGC 1333 measurement is consistent with the Barnard~1 result, but has a larger uncertainty. We do not see evidence for the very high $^{15}$N enhancements seen in cometary CN. Sensitive observations of a larger, carefully selected sample of prestellar cores with varying temperatures and gas densities can significantly improve our understanding of the nitrogen fractionation in the local interstellar medium and its relation to the isotopic ratios measured in various solar system reservoirs.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

preprint2010arXiv

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

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

preprint2009arXiv

Deuterium chemistry in the Orion Bar PDR - "warm" chemistry starring CH2D+

High levels of deuterium fractionation in gas-phase molecules are usually associated with cold regions, such as prestellar cores. Significant fractionation ratios are also observed in hot environments such as hot cores or hot corinos, where they are believed to be produced by the evaporation of the icy mantles surrounding dust grains, and thus are remnants of a previous cold (either gas-phase or grain surface) chemistry. The recent detection of DCN towards the Orion Bar, in a clump at a characteristic temperature of 70K, has shown that high deuterium fractionation can also be detected in PDRs. The Orion Bar clumps thus appear as a good environment for the observational study of deuterium fractionation in luke-warm gas, allowing to validate chemistry models in a different temperature range, where dominating fractionation processes are predicted to be different than in cold gas (< 20K). We aimed at studying observationally in detail the chemistry at work in the Orion Bar PDR, to understand if DCN is produced by ice mantle evaporation, or is the result of warm gas-phase chemistry, involving the CH2D+ precursor ion (which survives higher temperatures than the usual H2D+ precursor). Using the APEX and the IRAM 30m telescopes, we targetted selected deuterated species towards two clumps in the Orion Bar. We confirmed the detection of DCN and detected two new deuterated molecules (DCO+ and HDCO) towards one clump in the Orion Bar PDR. Significant deuterium fractionations are found for HCN and H2CO, but a low fractionation in HCO+. We also give upper limits for other molecules relevant for the deuterium chemistry. (...) We show evidence that warm deuterium chemistry driven by CH2D+ is at work in the clumps.

preprint2009arXiv

The chemical diversity of comets

A fundamental question in cometary science is whether the different dynamical classes of comets have different chemical compositions, which would reflect different initial conditions. From the ground or Earth orbit, radio and infrared spectroscopic observations of a now significant sample of comets indeed reveal deep differences in the relative abundances of cometary ices. However, no obvious correlation with dynamical classes is found. Further results come, or are expected, from space exploration. Such investigations, by nature limited to a small number of objects, are unfortunately focussed on short-period comets (mainly Jupiter-family). But these in situ studies provide "ground truth" for remote sensing. We discuss the chemical differences in comets from our database of spectroscopic radio observations, which has been recently enriched by several Jupiter-family and Halley-type comets.

preprint2004arXiv

Radial Distribution of Dust Grains Around HR 4796A

We present high-dynamic-range images of circumstellar dust around HR 4796A that were obtained with MIRLIN at the Keck II telescope at lambda = 7.9, 10.3, 12.5 and 24.5 um. We also present a new continuum measurement at 350 um obtained at the Caltech Submillimeter Observatory. Emission is resolved in Keck images at 12.5 and 24.5 um with PSF FWHM's of 0.37" and 0.55", respectively, and confirms the presence of an outer ring centered at 70 AU. Unresolved excess infrared emission is also detected at the stellar position and must originate well within 13 AU of the star. A model of dust emission fit to flux densities at 12.5, 20.8, and 24.5 um indicates dust grains are located 4(+3/-2) AU from the star with effective size, 28+/-6 um, and an associated temperature of 260+/-40 K. We simulate all extant data with a simple model of exozodiacal dust and an outer exo-Kuiper ring. A two-component outer ring is necessary to fit both Keck thermal infrared and HST scattered-light images. Bayesian parameter estimates yield a total cross-sectional area of 0.055 AU^2 for grains roughly 4 AU from the star and an outer-dust disk composed of a narrow large-grain ring embedded within a wider ring of smaller grains. The narrow ring is 14+/-1 AU wide with inner radius 66+/-1 AU and total cross-sectional area 245 AU^2. The outer ring is 80+/-15 AU wide with inner radius 45+/-5 AU and total cross-sectional area 90 AU^2. Dust grains in the narrow ring are about 10 times larger and have lower albedos than those in the wider ring. These properties are consistent with a picture in which radiation pressure dominates the dispersal of an exo-Kuiper belt.

preprint2003arXiv

The line-of-sight distribution of water in the SgrB2 complex

We report the detection, with the Caltech Submillimeter Observatory, of the 894-GHz HDO(111-000) transition, observed in absorption against the background continuum emission of the SgrB2 cores M and N. Radiative transfer modeling of this feature, together with the published data set of mm and submm HDO and H2O-18 transitions, suggests that ground-state absorption features from deuterated and non-deuterated water trace different gas components along the line of sight. In particular, while the HDO line seems to be produced by the large column densities of gas located in the SgrB2 warm envelope, the H2O-18 ground-state transition detected by SWAS and KAO at 548 GHz (Neufeld et al. 2000; Zmuidzinas et al. 1995a) is instead a product of the hot, diffuse, thin gas layer lying in the foreground of the SgrB2 complex.