Source author record

G. A. Blake

G. A. Blake 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

49works
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

49 published item(s)

preprint2022arXiv

Airfall on Comet 67P/Churyumov-Gerasimenko

We here study the transfer process of material from one hemisphere to the other (deposition of airfall material) on an active comet nucleus, specifically 67P/Churyumov-Gerasimenko. Our goals are to: 1) quantify the thickness of the airfall debris layers and how it depends on the location of the target area, 2) determine the amount of $\mathrm{H_2O}$ and $\mathrm{CO_2}$ ice that are lost from icy dust assemblages of different sizes during transfer through the coma, and 3) estimate the relative amount of vapor loss in airfall material after deposition in order to understand what locations are expected to be more active than others on the following perihelion approach. We use various numerical simulations, that include orbit dynamics, thermophysics of the nucleus and of individual coma aggregates, coma gas kinetics and hydrodynamics, as well as dust dynamics due to gas drag, to address these questions. We find that the thickness of accumulated airfall material varies substantially with location, and typically is of the order $0.1$-$1\,\mathrm{m}$. The airfall material preserves substantial amounts of water ice even in relatively small (cm-sized) coma aggregates after a rather long ($12\,\mathrm{h}$) residence in the coma. However, $\mathrm{CO_2}$ is lost within a couple of hours even in relatively large (dm-sized) aggregates, and is not expected to be an important component in airfall deposits. We introduce reachability and survivability indices to measure the relative capacity of different regions to simultaneously collect airfall and to preserve its water ice until the next perihelion passage, thereby grading their potential of contributing to comet activity during the next perihelion passage.

preprint2021arXiv

Observing Carbon & Oxygen Carriers in Protoplanetary Disks at Mid-infrared Wavelengths

Infrared observations probe the warm gas in the inner regions of planet-forming disks around young sun-like, T Tauri stars. In these systems, H$_2$O, OH, CO, CO$_2$, C$_2$H$_2$, and HCN have been widely observed. However, the potentially abundant carbon carrier CH$_4$ remains largely unconstrained. The James Webb Space Telescope (JWST) will be able to characterize mid-infrared fluxes of CH$_4$ along with several other carriers of carbon and oxygen. In anticipation of the JWST mission, we model the physical and chemical structure of a T Tauri disk to predict the abundances and mid-infrared fluxes of observable molecules. A range of compositional scenarios are explored involving the destruction of refractory carbon materials and alterations to the total elemental (volatile and refractory) C/O ratio. Photon-driven chemistry in the inner disk surface layers largely destroys the initial carbon and oxygen carriers. This causes models with the same physical structure and C/O ratio to have similar steady state surface compositions, regardless of the initial chemical abundances. Initial disk compositions are better preserved in the shielded inner disk midplane. The degree of similarity between the surface and midplane compositions in the inner disk will depend on the characteristics of vertical mixing at these radii. Our modeled fluxes of observable molecules respond sensitively to changes in the disk gas temperature, inner radius, and the total elemental C/O ratio. As a result, mid-infrared observations of disks will be useful probes of these fundamental disk parameters, including the C/O ratio, which can be compared to values determined for planetary atmospheres.

preprint2014arXiv

Dynamics of CO in Amorphous Water Ice Environments

The long-timescale behavior of adsorbed carbon monoxide on the surface of amorphous water ice is studied under dense cloud conditions by means of off-lattice, on-the-fly, kinetic Monte Carlo simula- tions. It is found that the CO mobility is strongly influenced by the morphology of the ice substrate. Nanopores on the surface provide strong binding sites which can effectively immobilize the adsorbates at low coverage. As the coverage increases, these strong binding sites are gradually occupied leav- ing a number of admolecules with the ability to diffuse over the surface. Binding energies, and the energy barrier for diffusion are extracted for various coverages. Additionally, the mobility of CO is determined from isothermal desorption experiments. Reasonable agreement on the diffusivity of CO is found with the simulations. Analysis of the 2152 cm$^{-1}$, polar CO band supports the computational findings that the pores in the water ice provide the strongest binding sites and dominate diffusion at low temperatures.

preprint2014arXiv

Herschel observations of Extra-Ordinary Sources: H2S as a Probe of Dense Gas and Possibly Hidden Luminosity Toward the Orion KL Hot Core

We present Herschel/HIFI observations of the light hydride H$_{2}$S obtained from the full spectral scan of the Orion Kleinmann-Low nebula (Orion KL) taken as part of the HEXOS GT key program. In total, we observe 52, 24, and 8 unblended or slightly blended features from H$_{2}$$^{32}$S, H$_{2}$$^{34}$S, and H$_{2}$$^{33}$S, respectively. We only analyze emission from the so called hot core, but emission from the plateau, extended ridge, and/or compact ridge are also detected. Rotation diagrams for ortho and para H$_{2}$S follow straight lines given the uncertainties and yield T$_{\rm rot}$=141$\pm$12 K. This indicates H$_{2}$S is in LTE and is well characterized by a single kinetic temperature or an intense far-IR radiation field is redistributing the population to produce the observed trend. We argue the latter scenario is more probable and find that the most highly excited states (E$_{\rm up}$>1000 K) are likely populated primarily by radiation pumping. We derive an H$_{2}$$^{32}$S column density, N$_{\rm tot}$(H$_{2}$$^{32}$S)=9.5$\pm$1.9$\times$10$^{17}$ cm$^{-2}$, gas kinetic temperature, T$_{\rm kin}$=120$\pm^{13}_{10}$ K, and constrain the H$_{2}$ volume density, n$_{H2}$>9$\times$10$^{7}$ cm$^{-3}$, for the H$_{2}$S emitting gas. These results point to an H$_{2}$S origin in markedly dense, heavily embedded gas, possibly in close proximity to a hidden self-luminous source (or sources), which are conceivably responsible for Orion KL's high luminosity. We also derive an H$_{2}$S ortho/para ratio of 1.7$\pm$0.8 and set an upper limit for HDS/H2S of <4.9$\times$10$^{-3}$.

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.

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

Evidence for a Snow Line Beyond the Transitional Radius in the TW Hya Protoplanetary Disk

We present an observational reconstruction of the radial water vapor content near the surface of the TW Hya transitional protoplanetary disk, and report the first localization of the snow line during this phase of disk evolution. The observations are comprised of Spitzer-IRS, Herschel-PACS, and Herschel-HIFI archival spectra. The abundance structure is retrieved by fitting a two-dimensional disk model to the available star+disk photometry and all observed H2O lines, using a simple step-function parameterization of the water vapor content near the disk surface. We find that water vapor is abundant (~10^{-4} per H2) in a narrow ring, located at the disk transition radius some 4AU from the central star, but drops rapidly by several orders of magnitude beyond 4.2 AU over a scale length of no more than 0.5AU. The inner disk (0.5-4AU) is also dry, with an upper limit on the vertically averaged water abundance of 10^{-6} per H2. The water vapor peak occurs at a radius significantly more distant than that expected for a passive continuous disk around a 0.6 Msun star, representing a volatile distribution in the TW Hya disk that bears strong similarities to that of the solar system. This is observational evidence for a snow line that moves outward with time in passive disks, with a dry inner disk that results either from gas giant formation or gas dissipation and a significant ice reservoir at large radii. The amount of water present near the snow line is sufficient to potentially catalyze the (further) formation of planetesimals and planets at distances beyond a few AU.

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

First Detection of Near-Infrared Line Emission from Organics in Young Circumstellar Disks

We present an analysis of high-resolution spectroscopy of several bright T Tauri stars using the VLT/CRIRES and Keck/NIRSPEC spectrographs, revealing the first detections of emission from HCN and C2H2 in circumstellar disks at near-infrared wavelengths. Using advanced data reduction techniques we achieve a dynamic range with respect to the disk continuum of ~500 at 3 microns, revealing multiple emission features of H2O, OH, HCN, and C2H2. We also present stringent upper limits for two other molecules thought to be abundant in the inner disk, CH4 and NH3. Line profiles for the different detected molecules are broad but centrally peaked in most cases, even for disks with previously determined inclinations of greater than 20 degrees, suggesting that the emission has both a Keplerian and non-Keplerian component as observed previously for CO emission. We apply two different modeling strategies to constrain the molecular abundances and temperatures: we use a simplified single-temperature LTE slab model with a Gaussian line profile to make line identifications and determine a best-fit temperature and initial abundance ratios, and we compare these values with constraints derived from a detailed disk radiative transfer model assuming LTE excitation but utilizing a realistic temperature and density structure. Abundance ratios from both sets of models are consistent with each other and consistent with expected values from theoretical chemical models, and analysis of the line shapes suggests the molecular emission originates from within a narrow region in the inner disk (R < 1 AU).

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.

preprint2011arXiv

A Spitzer Survey of Mid-Infrared Molecular Emission from Protoplanetary Disks II: Correlations and LTE Models

We present an analysis of Spitzer-IRS observations of H2O, OH, HCN, C2H2, and CO2 emission, and Keck-NIRSPEC observations of CO emission, from a diverse sample of T Tauri and Herbig Ae/Be circumstellar disks. We find that detections and strengths of most mid-IR molecular emission features are correlated with each other, suggesting a common origin and similar excitation conditions. We note that the line detection efficiency is anti-correlated with the 13/30 um SED spectral slope, which is a measure of the degree of grain settling in the disk atmosphere. We also note a correlation between detection efficiency and H-alpha equivalent width, and tentatively with accretion rate, suggesting that accretional heating contributes to line excitation. If detected, H2O line fluxes are correlated with the mid-IR continuum flux, and other co-varying system parameters, such as L_star. However, significant sample variation, especially in molecular line ratios, remains. LTE models of the H2O emission show that line strength is primarily related to the best-fit emitting area, and this accounts for most source-to-source variation in H2O emitted flux. Best-fit temperatures and column densities cover only a small range of parameter space, near 10^{18} cm-2 and 450 K for all sources, suggesting a high abundance of H2O in many planet-forming regions. Other molecules have a range of excitation temperatures from ~500-1500 K, also consistent with an origin in planet-forming regions. We find molecular ratios relative to water of ~10^{-3} for all molecules, with the exception of CO, for which n(CO)/n(H2O)~1. However, LTE fitting caveats and differences in the way thermo-chemical modeling results are reported make comparisons with such models difficult, and highlight the need for additional observations coupled with the use of line-generating radiative transfer codes.

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

Ice and Dust in the Quiescent Medium of Isolated Dense Cores

The relation between ices in the envelopes and disks surrounding YSOs and those in the quiescent interstellar medium is investigated. For a sample of 31 stars behind isolated dense cores, ground-based and Spitzer spectra and photometry in the 1-25 um wavelength range are combined. The baseline for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the separation of continuum and feature extinction. The extinction between 13-25 um is ~50% relative to that at 2.2 um. The strengths of the 6.0 and 6.85 um absorption bands are in line with those of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH4+, HCOOH, H2CO and NH3, are readily formed in the dense core phase, before stars form. The 3.53 um C-H stretching mode of solid CH3OH was discovered. The CH3OH/H2O abundance ratios of 5-12% are larger than upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, therefore depends on the environment. Signs of thermal and energetic processing that were found toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of the 9.7 um band of silicates relative to the continuum extinction at 2.2 um is significantly shallower than in the diffuse interstellar medium. This extends the results of Chiar et al. (2007) to a larger sample and higher extinctions.

preprint2011arXiv

Temporal and Spatial Aspects of Gas Release During the 2010 Apparition of Comet 103P/Hartley-2

We report measurements of eight primary volatiles (H2O, HCN, CH4, C2H6, CH3OH, C2H2, H2CO, and NH3) and two product species (OH and NH2) in comet 103P/Hartley-2 using high dispersion infrared spectroscopy. We quantified the long- and short-term behavior of volatile release over a three-month interval that encompassed the comet's close approach to Earth, its perihelion passage, and flyby of the comet by the Deep Impact spacecraft during the EPOXI mission. We present production rates for individual species, their mixing ratios relative to water, and their spatial distributions in the coma on multiple dates. The production rates for water, ethane, HCN, and methanol vary in a manner consistent with independent measures of nucleus rotation, but mixing ratios for HCN, C2H6, & CH3OH are independent of rotational phase. Our results demonstrate that the ensemble average composition of gas released from the nucleus is well defined, and relatively constant over the three-month interval (September 18 through December 17). If individual vents vary in composition, enough diverse vents must be active simultaneously to approximate (in sum) the bulk composition of the nucleus. The released primary volatiles exhibit diverse spatial properties which favor the presence of separate polar and apolar ice phases in the nucleus, establish dust and gas release from icy clumps (and also, directly from the nucleus), and provide insights into the driver for the cyanogen (CN) polar jet. The spatial distributions of C2H6 & HCN along the near-polar jet (UT 19.5 October) and nearly orthogonal to it (UT 22.5 October) are discussed relative to the origin of CN. The ortho-para ratio (OPR) of water was 2.85 \pm 0.20; the lower bound (2.65) defines Tspin > 32 K. These values are consistent with results returned from ISO in 1997.

preprint2011arXiv

Warm dust resolved in the cold disk around T Cha with VLTI/AMBER

The transition between massive Class II circumstellar disks and Class III debris disks, with dust residuals, has not yet been clearly understood. Disks are expected to dissipate with time, and dust clearing in the inner regions can be the consequence of several mechanisms. Planetary formation is one of them that will possibly open a gap inside the disk. According to recent models based on photometric observations, T Cha is expected to present a large gap within its disk, meaning that an inner dusty disk is supposed to have survived close to the star. We investigate this scenario with new near-infrared interferometric observations. We observed T Cha in the H and K bands using the AMBER instrument at VLTI and used the MCFOST radiative transfer code to model the SED of T Cha and the interferometric observations simultaneously and to test the scenario of an inner dusty structure. We also used a toy model of a binary to check that a companion close to the star can reproduce our observations. The scenario of a close (few mas) companion cannot satisfactorily reproduce the visibilities and SED, while a disk model with a large gap and an inner ring producing the bulk of the emission (in H and K-bands) close to 0.1 AU is able to account for all the observations. With this study, the presence of an optically thick inner dusty disk close to the star and dominating the H and K- bands emission is confirmed. According to our model, the large gap extends up to ~ 7.5 AU. This points toward a companion (located at several AU) gap-opening scenario to explain the morphology of T Cha.

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

Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

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 of the Herschel Key Program 'Dust, Ice and Gas in Time': Dust and Gas Spectroscopy of HD 100546

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 μ\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 μ\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 μ\m and [OI] at 63 and 145 μ\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

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

Fundamental Vibrational Transition of CO during the Outburst of EX Lupi in 2008

We report monitoring observations of the T Tauri star EX Lupi during its outburst in 2008 in the CO fundamental band at 4.6-5.0 um. The observations were carried out at the VLT and the Subaru Telescope at six epochs from April to August 2008, covering the plateau of the outburst and the fading phase to a quiescent state. The line flux of CO emission declines with the visual brightness of the star and the continuum flux at 5 um, but composed of two subcomponents that decay with different rates. The narrow line emission (50 km s-1 in FWHM) is near the systemic velocity of EX Lupi. These emission lines appear exclusively in v=1-0. The line widths translate to a characteristic orbiting radius of 0.4 AU. The broad line component (FWZI ~ 150 km s-1) is highly excited upto v<=6. The line flux of the component decreases faster than the narrow line emission. Simple modeling of the line profiles implies that the broad-line emitting gas is orbiting around the star at 0.04-0.4 AU. The excitation state, the decay speed of the line flux, and the line profile, indicate that the broad-line emission component is physically distinct from the narrow-line emission component, and more tightly related to the outburst event.

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 observations of high-J CO lines in the NGC 1333 low-mass star-forming region

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

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

Hunting for millimeter flares from magnetic reconnection in pre-main sequence spectroscopic binaries

Recent observations of the low-mass pre-main sequence, eccentric spectroscopic binaries DQ Tau and V773 Tau A reveal that their millimeter spectrum is occasionally dominated by flares from non-thermal emission processes. The transient activity is believed to be synchrotron in nature, resulting from powerful magnetic reconnection events when the separate magnetic structures of the binary components are capable of interacting and forced to reorganize, typically near periastron. We conducted the first systematic study of the millimeter variability toward a sample of 12 PMS spectroscopic binaries with the aim to characterize the proliferation of flares amongst sources likely to experience similar interbinary reconnection events. The source sample consists of short-period, close-separation binaries that possess either a high orbital eccentricity or a circular orbit. Using the MAMBO2 array on the IRAM 30m telescope, we carried out continuous monitoring at 1.25 mm over a 4-night period during which all of the high-eccentricity binaries approached periastron. We also obtained simultaneous optical VRI measurements, since a strong link is often observed between stellar reconnection events and optical brightenings. UZ Tau E is the only source to be detected at millimeter wavelengths: it exhibited significant variation; it is also the only source to undergo strong simultaneous optical variability. The binary possesses the largest orbital eccentricity in the current sample, a predicted factor in star-star magnetic interaction events. With orbital parameters and variable accretion activity similar to DQ Tau, the millimeter behavior of UZ Tau E draws many parallels to the DQ Tau model for colliding magnetospheres. However, on the basis of our observations alone, we cannot determine whether the variability is repetitive, or if it could also be due to variable free-free emission in an ionized wind.

preprint2010arXiv

Origin of the hot gas in low-mass protostars: Herschel-PACS spectroscopy of HH 46

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

preprint2010arXiv

Recurring millimeter flares as evidence for star-star magnetic reconnection events in the DQ Tauri PMS binary system

Observations of the T Tauri spectroscopic binary DQ Tau in April 2008 captured an unusual flare at 3 mm, which peaked at an observed max flux of 0.5 Jy (about 27x the quiescent value). Here we present follow-up mm observations that demonstrate a periodicity to the phenomenon. While monitoring 3 new periastron encounters, we detect flares within 17.5 hrs (or 4.6%) of the orbital phase of the first reported flare, and we constrain the main emitting region to a stellar height of 3.7-6.8 Rstar. The recorded activity is consistent with the proposed picture for synchrotron emission initiated by a magnetic reconnection event when the two stellar magnetospheres of the highly eccentric (e=0.556) binary are believed to collide near periastron as the stars approach a minimum separation of 8 Rstar (~13 Rsolar). The similar light curve decay profiles allow us to estimate an average flare duration of 30 hrs. Assuming one mm flare per orbit, DQ Tau could spend approximately 8% of its 15.8-d orbital period in an elevated flux state. Our analysis of the mm emission provides an upper limit of 5% on the linear polarization. We discuss the extent to which a severely entangled magnetic field structure and Faraday rotation effects are likely to reduce the observed polarization fraction. We also predict that, for the current picture, the stellar magnetospheres must be misaligned at a significant angle or, alternatively, that the topologies of the outer magnetospheres are poorly described by a well-ordered dipole inside a radius of 7 Rstar. Finally, to investigate whether reorganization of the magnetic field during the interaction affects mass accretion, we also present simultaneous optical (VRI) monitoring, as an established tracer of accretion activity in this system. We find that an accretion event can occur coincident in both time and duration with the synchrotron fallout of a magnetic reconnection event.

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

The methanol lines and hot core of OMC2-FIR4, an intermediate-mass protostar, with Herschel-HIFI

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

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

C2D Spitzer-IRS spectra of disks around T Tauri stars. IV. Crystalline silicates

Dust grains in the planet forming regions around young stars are expected to be heavily processed due to coagulation, fragmentation and crystallization. This paper focuses on the crystalline silicate dust grains in protoplanetary disks. As part of the Cores to Disks Legacy Program, we obtained more than a hundred Spitzer/IRS spectra of TTauri stars. More than 3/4 of our objects show at least one crystalline silicate emission feature that can be essentially attributed to Mg-rich silicates. Observational properties of the crystalline features seen at lambda > 20 mu correlate with each other, while they are largely uncorrelated with the properties of the amorphous silicate 10 mu feature. This supports the idea that the IRS spectra essentially probe two independent disk regions: a warm zone (< 1 AU) emitting at lambda ~ 10 mu and a much colder region emitting at lambda > 20 mu (< 10 AU). We identify a crystallinity paradox, as the long-wavelength crystalline silicate features are 3.5 times more frequently detected (~55 % vs. ~15%) than the crystalline features arising from much warmer disk regions. This suggests that the disk has an inhomogeneous dust composition within ~10 AU. The abundant crystalline silicates found far from their presumed formation regions suggests efficient outward radial transport mechanisms in the disks. The analysis of the shape and strength of both the amorphous 10 mu feature and the crystalline feature around 23 mu provides evidence for the prevalence of micron-sized grains in upper layers of disks. Their presence in disk atmospheres suggests efficient vertical diffusion, likely accompanied by grain-grain fragmentation to balance the efficient growth expected. Finally, the depletion of submicron-sized grains points toward removal mechanisms such as stellar winds or radiation pressure.

preprint2009arXiv

Evidence for dust clearing through resolved submillimeter imaging

Mid-infrared spectrophotometric observations have revealed a small sub-class of circumstellar disks with spectral energy distributions (SEDs) suggestive of large inner gaps with low dust content. However, such data provide only an indirect and model-dependent method of finding central holes. Imaging of protoplanetry disks provides an independent check of SED modeling. We present here the direct characterization of three 33-47 AU radii inner gaps, in the disks around LkHa 330, SR 21N and HD 135344B, via 340 GHz (880 micron) dust continuum aperture synthesis observations obtained with the Submillimeter Array (SMA). The large gaps are fully resolved at ~0\farcs3 by the SMA observations and mostly empty of dust, with less than 1 - 7.5 x 10^-6 Msolar of fine grained solids inside the holes. Gas (as traced by atomic accretion markers and CO 4.7 micron rovibrational emission) is still present in the inner regions of all three disks. For each, the inner hole exhibits a relatively steep rise in dust emission to the outer disk, a feature more likely to originate from the gravitational influence of a companion body than from a process expected to show a more shallow gradient like grain growth. Importantly, the good agreement of the spatially resolved data and spectrophotometry-based models lends confidence to current interpretations of SEDs, wherein the significant dust emission deficits arise from disks with inner gaps or holes. Further SED-based searches can therefore be expected to yield numerous additional candidates that can be examined at high spatial resolution.

preprint2009arXiv

Radiative transfer models of mid-infrared H2O lines in the Planet-forming Region of Circumstellar Disks

The study of warm molecular gas in the inner regions of protoplanetary disks is of key importance for the study of planet formation and especially for the transport of H2O and organic molecules to the surfaces of rocky planets/satellites. Recent Spitzer observations have shown that the mid-infrared spectra of protoplanetary disks are covered in emission lines due to water and other molecules. Here, we present a non-LTE 2D radiative transfer model of water lines in the 10-36 mum range that can be used to constrain the abundance structure of water vapor, given an observed spectrum, and show that an assumption of local thermodynamic equilibrium (LTE) does not accurately estimate the physical conditions of the water vapor emission zones. By applying the model to published Spitzer spectra we find that: 1) most water lines are subthermally excited, 2) the gas-to-dust ratio must be one to two orders of magnitude higher than the canonical interstellar medium ratio of 100-200, and 3) the gas temperature must be higher than the dust temperature, and 4) the water vapor abundance in the disk surface must be truncated beyond ~ 1 AU. A low efficiency of water formation below ~ 300 K may naturally result in a lower water abundance beyond a certain radius. However, we find that chemistry, may not be sufficient to produce an abundance drop of many orders of magnitude and speculate that the depletion may also be caused by vertical turbulent diffusion of water vapor from the superheated surface to regions below the snow line, where the water can freeze out and be transported to the midplane as part of the general dust settling. Such a vertical cold finger effect is likely to be efficient due to the lack of a replenishment mechanism of large, water-ice coated dust grains to the disk surface.

preprint2008arXiv

Probing dust grain evolution in IM Lupi's circumstellar disc. Multi-wavelength observations and modelling of the dust disc

We present a panchromatic study, involving a multiple technique approach, of the circumstellar disc surrounding the T Tauri star IM Lupi (Sz 82). We have undertaken a comprehensive observational study of IM Lupi using photometry, spectroscopy, millimetre interferometry and multi-wavelength imaging. For the first time, the disc is resolved from optical and near-infrared wavelengths in scattered light, to the millimetre regime in thermal emission. Our data-set, in conjunction with existing photometric data, provides an extensive coverage of the spectral energy distribution, including a detailed spectrum of the silicate emission bands. We have performed a simultaneous modelling of the various observations, using the radiative transfer code MCFOST, and analysed a grid of models over a large fraction of the parameter space via Bayesian inference. We have constructed a model that can reproduce all of the observations of the disc. Our analysis illustrates the importance of combining a wide range of observations in order to fully constrain the disc model, with each observation providing a strong constraint only on some aspects of the disc structure and dust content. Quantitative evidence of dust evolution in the disc is obtained: grain growth up to millimetre-sized particles, vertical stratification of dust grains with micrometric grains close to the disc surface and larger grains which have settled towards the disc midplane, and possibly the formation of fluffy aggregates and/or ice mantles around grains.

preprint2007arXiv

Abundant crystalline silicates in the disk of a very low mass star

We announce the discovery of SST-Lup3-1, a very low mass star close to the brown dwarf boundary in Lupus III with a circum(sub)stellar disk, discovered by the `Cores to Disks' Spitzer Legacy Program from mid-, near-infrared and optical data, with very conspicuous crystalline silicate features in its spectrum. It is the first of such objects with a full 5 to 35 micron spectrum taken with the IRS and it shows strong 10 and 20 micron silicate features with high feature to continuum ratios and clear crystalline features out to 33 micron. The dust in the disk upper layer has a crystalline silicate grain fraction between 15% and 33%, depending on the assumed dust continuum. The availability of the full Spitzer infrared spectrum allows an analysis of the dust composition as a function of temperature and position in the disk. The hot (~ 300 K) dust responsible for the 10 micron feature consists of a roughly equal mix of small (~ 0.1 micron) and large (~ 1.5 micron) grains, whereas the cold (~ 70 K) dust responsible for the longer wavelength silicate features contains primarily large grains (> 1 micron). Since the cold dust emission arises from deeper layers in the inner (< 3 AU) disk as well as from the surface layers of the outer (3-5 AU) disk, this provides direct evidence for combined grain growth and settling in the disk. The inferred crystalline mass fractions in the two components are comparable. Since only the inner 0.02 AU of the disk is warm enough to anneal the amorphous silicate grains, even the lowest fraction of 15% of crystalline material requires either very efficient mixing or other formation mechanisms.

preprint2007arXiv

Cold Disks: Spitzer Spectroscopy of Disks around Young Stars with Large Gaps

We have identified four circumstellar disks with a deficit of dust emission from their inner 15-50 AU. All four stars have F-G spectral type, and were uncovered as part of the Spitzer Space Telescope ``Cores to Disks'' Legacy Program Infrared Spectrograph (IRS) first look survey of ~100 pre-main sequence stars. Modeling of the spectral energy distributions indicates a reduction in dust density by factors of 100-1000 from disk radii between ~0.4 and 15-50 AU, but with massive gas-rich disks at larger radii. This large contrast between the inner and outer disk has led us to use the term `cold disks' to distinguish these unusual systems. However, hot dust [0.02-0.2 Mmoon] is still present close to the central star (R ~0.8 AU). We introduce the 30/13 micron, flux density ratio as a new diagnostic for identifying cold disks. The mechanisms for dust clearing over such large gaps are discussed. Though rare, cold disks are likely in transition from an optically thick to an optically thin state, and so offer excellent laboratories for the study of planet formation.

preprint2007arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds: VI. Perseus Observed with MIPS

We present observations of 10.6 square degrees of the Perseus molecular cloud at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). The image mosaics show prominent, complex extended emission dominated by illuminating B stars on the East side of the cloud, and by cold filaments of 160 micron emission on the West side. Of 3950 point sources identified at 24 microns, 1141 have 2MASS counterparts. A quarter of these populate regions of the Ks vs. Ks-[24] diagram that are distinct from stellar photospheres and background galaxies, and thus are likely to be cloud members with infrared excess. Nearly half (46%) of these 24 micron excess sources are distributed outside the IC 348 and NGC 1333 clusters. NGC 1333 shows the highest fraction of stars with flat or rising spectral energy distributions (28%), while Class II SEDs are most common in IC 348. These results are consistent with previous relative age determinations for the two clusters. The intercluster region contains several tightly clumped (r~0.1 pc) young stellar aggregates whose members exhibit a wide variety of infrared spectral energy distributions characteristic of different circumstellar environments. One possible explanation is a significant age spread among the aggregate members, such that some have had time to evolve more than others. Alternatively, if the aggregate members all formed at roughly the same time, then remarkably rapid circumstellar evolution would be required to account for the association of Class I and Class III sources at ages <~1 Myr. We highlight important results for several other objects as well (full abstract in the paper).

preprint2006arXiv

Probing Protoplanetary Disks with Silicate Emission: Where is the silicate emission zone?

Recent results indicate that the grain size and crystallinity inferred from observations of silicate features may be correlated with spectral type of the central star and/or disk geometry. In this paper, we show that grain size, as probed by the 10 um silicate feature peak-to-continuum and 11.3-to-9.8 um flux ratios, is inversely proportional to log L_star. These trends can be understood using a simple two-layer disk model for passive irradiated flaring disks, CGPLUS. We find that the radius, R_10, of the 10 um silicate emission zone in the disk goes as (L_star/L_sun)^0.56, with slight variations depending on disk geometry (flaring angle, inner disk radius). The observed correlations, combined with simulated emission spectra of olivine and pyroxene mixtures, imply a grain size dependence on luminosity. Combined with the fact that R_10 is smaller for less luminous stars, this implies that the apparent grain size of the emitting dust is larger for low-luminositysources. In contrast, our models suggest that the crystallinity is only marginally affected, because for increasing luminosity, the zone for thermal annealing (assumed to be at T>800 K) is enlarged by roughly the same factor as the silicate emission zone. The observed crystallinity is affected by disk geometry, however, with increased crystallinity in flat disks. The apparent crystallinity may also increase with grain growth due to a corresponding increase in contrast between crystalline and amorphous silicate emission bands.

preprint2005arXiv

Protostellar holes: Spitzer Space Telescope observations of the protostellar binary IRAS16293-2422

Mid-infrared (23-35 micron) emission from the deeply embedded "Class 0" protostar IRAS16293-2422 is detected with the Spitzer Space Telescope infrared spectrograph. A detailed radiative transfer model reproducing the full spectral energy distribution (SED) from 23 micron to 1.3 mm requires a large inner cavity of radius 600 AU in the envelope to avoid quenching the emission from the central sources. This is consistent with a previous suggestion based on high angular resolution millimeter interferometric data. An alternative interpretation using a 2D model of the envelope with an outflow cavity can reproduce the SED but not the interferometer visibilities. The cavity size is comparable to the centrifugal radius of the envelope and therefore appears to be a natural consequence of the rotation of the protostellar core, which has also caused the fragmentation leading to the central protostellar binary. With a large cavity such as required by the data, the average temperature at a given radius does not increase above 60-80 K and although hot spots with higher temperatures may be present close to each protostar, these constitute a small fraction of the material in the inner envelope. The proposed cavity will also have consequences for the interpretation of molecular line data, especially of complex species probing high temperatures in the inner regions of the envelope.

preprint2004arXiv

Methane Abundance Variations toward the Massive Protostar NGC 7538 : IRS9

Absorption and emission lines originating from the nu3 C-H stretching manifold of gas phase CH4 were discovered in the high resolution (R=25,000) infrared L band spectrum along the line of sight toward NGC 7538 : IRS9. These observations provide a diagnostic of the complex dynamics and chemistry in a massive star forming region. The line shapes resemble P Cygni profiles with the absorption and emission components shifted by ~7 km/s with respect to the systemic velocity. Similar velocity components were observed in CO at 4.7 um, but in contrast to CH4, the CO shows deep absorption due to a high velocity outflow as well as absorption at the systemic velocity due to the cold outer envelope. It is concluded that the gas phase CH4 abundance varies by an order of magnitude in this line of sight: it is low in the envelope and the outflow (X[CH4]<0.4e-6), and at least a factor of 10 larger in the central core. The discovery of solid CH4 in independent ground and space based data sets shows that methane is nearly entirely frozen onto grains in the envelope. It thus appears that CH4 is formed by grain surface reactions, evaporates into the gas phase in the warm inner regions of protostellar cores and is efficiently destroyed in shocks related to outflows.