Source author record

E. F. van Dishoeck

E. F. van Dishoeck appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.EP astro-ph astro-ph.IM physics.chem-ph physics.space-ph

Catalog footprint

What is connected

129works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Characterizing the physical and chemical properties of the Class I protostellar system Oph-IRS 44. Binarity, infalling streamers, and accretion shocks

(Abridged) In the low-mass star formation process, theoretical models predict that material from the infalling envelope could be shocked as it encounters the outer regions of the disk. Nevertheless, only a few protostars show evidence of these shocks at the disk-envelope interface, and the main formation path of shocked-related species is still unclear. We present new ALMA observations of IRS 44, a Class I source that has previously been associated with accretion shocks, taken at high angular resolution (0.1"). We target multiple molecular transitions of CO, H2CO, and simple sulfur-bearing species. In continuum emission, the binary nature of IRS 44 is observed for the first time at sub-millimeter wavelengths. Infalling signatures are seen for the CO line and the emission peaks at the edges of the continuum emission around IRS 44 B, the same region where bright SO and SO2 emission is seen. Weak CS and H2CO emission is observed, while OCS, H2S, and H2CS transitions are not detected. IRS 44 B seems to be more embedded than IRS 44 A, indicating a non-coeval formation scenario or the rejuvenation of source B due to late infall. CO emission is tracing the outflow component at large scales, infalling envelope material at intermediate scales, and two infalling streamer candidates are identified at disk scales. Infalling streamers might produce accretion shocks when they encounter the outer regions of the infalling-rotating envelope. These shocks heat the dust and release S-bearing species as well as promoting a lukewarm chemistry in the gas phase. With the majority of carbon locked in CO, there is little free C available to form CS and H2CS in the gas, leaving an oxygen-rich environment. The high column densities of SO and SO2 might be a consequence of two processes: direct thermal desorption from dust grains and gas-phase formation due to the availability of O and S.

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects' orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects' bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects' atmospheres. It also constrains these objects' bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2025arXiv

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) enables the characterisation of young self-luminous gas giants at previously inaccessible wavelengths, revealing physical processes in gas, dust, and clouds. We characterise the young planetary system TWA 27 (2M1207) in the mid-infrared (MIR), studying the atmosphere and disk spectra of the M9 brown dwarf TWA 27A and its L6 planetary-mass companion TWA 27b. We obtained data with the MIRI Medium Resolution Spectrometer (MRS) from 4.9 to 20 um, and MIRI Imaging in the F1000W and F1500W filters. We applied high-contrast imaging methods to extract the companion's spectral energy distribution up to 15 um at 0.78 arcsec separation and a contrast of 60. Combining these with published JWST/NIRSpec spectra, we analysed the 1-20 um range using self-consistent atmospheric grids and 0D slab models for molecular disk emission. The atmosphere of TWA 27A is well fitted by a BT-SETTL model with Teff 2780 K, log g 4.3, plus a 740 K blackbody for the inner disk rim. The disk shows at least 11 organic molecules, with no water or silicate dust emission detected. The atmosphere of TWA 27b is matched by a Teff 1400 K low-gravity model with extinction, best fit by the ExoREM grid. MIRI spectra and photometry for TWA 27b reveal a silicate cloud absorption feature between 8-10 um and significant (>5 sigma) infrared excess at 15 um consistent with circumplanetary disk emission. These MIRI observations provide new insights into TWA 27, revealing diverse features to study the formation and evolution of circumplanetary disks and young dusty atmospheres.

preprint2022arXiv

Gas temperature structure across transition disk cavities

[Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of $^{12}$CO, $^{13}$CO and C$^{18}$O to observationally constrain the gas temperature and surface density. Furthermore, we use the thermochemical code DALI to model the temperature and density structure of a typical transition disk. The $6-5/2-1$ line ratio in LkCa15 constrains the gas temperature in the emitting layers inside the dust cavity to be up to 65 K, warmer than in the outer disk at 20-30 K. For the HD 169142, the peak brightness temperature constrains the gas in the dust cavity of HD 169142 to be 170 K, whereas that in the outer disk is only 100 K. Models also show that a more luminous central star, a lower abundance of PAHs and the absence of a dusty inner disk increase the temperature of the emitting layers and hence the line ratio in the gas cavity. The gas column density in the LkCa15 dust cavity drops by a factor >2 compared to the outer disk, with an additional drop of an order of magnitude inside the gas cavity at 10 AU. In the case of HD 169142, the gas column density drops by a factor of 200$-$500 inside the gas cavity, which could be due to a massive companion of several M$_{\mathrm{J}}$. The broad dust-depleted gas region from 10-68 AU for LkCa15 may imply several lower mass planets. This work demonstrates that knowledge of the gas temperature is important to determine the gas surface density and thus whether planets, and if so what kind of planets, are the most likely carving the dust cavities.

preprint2022arXiv

H$_2$O distribution in the disc of HD 100546 and HD 163296: the role of dust dynamics and planet--disc interaction

[Abridged] Far-infrared observations with Herschel revealed a surprisingly low abundance of cold-water reservoirs in protoplanetary discs. On the other hand, a handful of discs show emission of hot water transitions excited at temperatures above a few hundred Kelvin. In particular, the protoplanetary discs around the Herbig Ae stars HD 100546 and HD 163296 show opposite trends in terms of cold versus hot water emission: in the first case, the ground-state transitions are detected and the high-J lines are undetected, while the trend is opposite in HD 163296. We performed a spectral analysis using the thermo-chemical model DALI. We find that HD 163296 is characterised by a water-rich (abundance $\gtrsim 10^{-5}$) hot inner disc (within the snowline) and a water-poor ($< 10^{-10}$) outer disc: the relative abundance may be due to the thermal desorption of icy grains that have migrated inward. Remarkably, the size of the H$_2$O emitting region corresponds to a narrow dust gap visible in the millimeter continuum at $r=10\,$au with ALMA. The low-J lines detected in HD 100546 instead imply an abundance of a few $10^{-9}$ in the cold outer disc ($> 40$ au). The emitting region of the cold H$_2$O transitions is spatially coincident with that of the H$_2$O ice previously seen in the near-infrared. Notably, millimetre observations with ALMA reveal the presence of a large dust gap between nearly 40 and 150 au, likely opened by a massive embedded protoplanet. In both discs, we find that the warm molecular layer in the outer region (beyond the snow line) is highly depleted of water molecules, implying an oxygen-poor chemical composition of the gas. We speculate that gas-phase oxygen in the outer disc is readily depleted and its distribution in the disc is tightly coupled to the dynamics of the dust grains.

preprint2022arXiv

Importance of source structure on complex organics emission II. Can disks explain lack of methanol emission from some low-mass protostars?

Some protostellar systems show little or no millimetre line emission of complex organics. This can be interpreted as a low abundance of these molecules, alternatively they could be present in the system but are not seen in the gas. The goal is to investigate the latter hypothesis for methanol. We will attempt to answer the question: Does the presence of a disk and optically thick dust reduce methanol emission even if methanol is abundant in the ices and gas? Using the radiative transfer code RADMC-3D, methanol emission lines from an envelope-only model and an envelope-plus-disk model are calculated and compared with each other and the observations. Methanol emission from the envelope-only model is always stronger than from the envelope-plus-disk model by at least a factor 2 as long as the disk radius is larger than 30 au (for L=8 L$_{\odot}$). In most cases, this is due to lower temperatures and, hence, the smaller amount of warm methanol inside the snow surface of the envelope-plus-disk model. The intensities drop by more than an order of magnitude for models including high mm opacity dust grains and disk radii of at least 50 au (for L=8 L$_{\odot}$) due to continuum over-subtraction. The line intensities from the envelope-only models overproduce the observations of protostars with lower methanol emission even with large dust optical depth effects. The envelope-plus-disk models can explain the bulk of the observations. However, they can only reproduce the observations of sources with high luminosities and low methanol emission when dust optical depth effects especially continuum over-subtraction in the disk becomes effective. Therefore, both the effects of disk and dust optical depth should be considered to explain the observations. In conclusion: Absence of methanol emission does not imply absence of methanol molecules in either gas or ice.

preprint2022arXiv

Importance of source structure on complex organics emission. I. Observations of CH$_3$OH from low-mass to high-mass protostars

Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.70\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs.

preprint2022arXiv

Self-consistent modelling of the dust component in protoplanetary and circumplanetary disks: the case of PDS 70

Direct observations of young stellar objects are important to test established theories of planet formation. PDS 70 is one of the few cases where robust evidence favours the presence of two planetary mass companions inside the gap of the transition disk. Those planets are believed to be going through the last stages of accretion from the protoplanetary disk, a process likely mediated by a circumplanetary disk (CPD). We aim to develop a three dimensional radiative transfer model for the dust component of the PDS 70 system which reproduces the system's global features observed at two different wavelengths: 855 $μ\, \mathrm{m}$ with ALMA and 1.25 $μ\, \mathrm{m}$ with VLT/SPHERE. We use this model to investigate the physical properties of the planetary companion PDS 70 c and its potential circumplanetary disk. We select initial values for the physical properties of the planet and CPD through appropriate assumptions about the nature and evolutionary stage of the object. We modify iteratively the properties of the protoplanetary disk until the predictions retrieved from the model are consistent with both data sets. We provide a model that jointly explains the global features of the PDS 70 system seen in submillimeter and polarised-scattered light. Our model suggests that spatial segregation of dust grains is present in the protoplanetary disk. The submillimeter modelling of the PDS 70 c source favours the presence of an optically thick CPD and places an upper limit to its dust mass of 0.7 $M_\oplus$. Furthermore, analysis of the thermal structure of the CPD demonstrates that the planet luminosity is the dominant heating mechanism of dust grains inside 0.6 au from the planet while heating by stellar photons dominates at larger planetocentric distances.

preprint2022arXiv

Survey of Orion Disks with ALMA (SODA) I: Cloud-level demographics of 873 protoplanetary disks

(Abridged) Surveys of protoplanetary disks in nearby star-forming regions (SFRs) have provided important information on their demographics, but due to their sample sizes, they cannot be used to study how disk properties vary with the environment. We conduct a survey of the unresolved millimeter continuum emission of 873 protoplanetary disks identified by Spitzer in the L1641 and L1647 regions of the Orion A cloud. This is the largest such survey yet, allowing us to identify even weak trends in the median disk mass as a function of position in the cloud and cluster membership. The sample detection rates and median masses are also compared to those of nearby (<300 pc) SFRs. The sample was observed with ALMA at 225 GHz, with a median rms of 0.08 mJy/beam (1.5 M$_{\oplus}$). The data were reduced and imaged using an innovative parallel data processing approach. We detect 58% (502/873) of the observed disks. This includes 20 disks with dust masses >100 M$_{\oplus}$, and two objects associated with extended dust emission. We infer a median disk dust mass in the full sample of $2.2^{+0.2}_{-0.2}$ M$_{\oplus}$. In L1641 and L1647 median dust masses are $2.1^{+0.2}_{-0.2}$ M$_{\oplus}$ and $2.6^{+0.4}_{-0.5}$ M$_{\oplus}$, respectively. The disk mass distribution of the full sample is similar to that of nearby low-mass SFRs at similar ages of 1-3 Myr. We find only weak trends in disk (dust) masses with galactic longitude and between the YSO clusters identified in the sample, with median masses varying by $\lesssim$ 50%. Age differences may explain these median disk mass variations. Apart from this, disk masses are essentially constant at scales of ~100 pc. This also suggests that the majority of disks, even in different SFRs, are formed with similar initial masses and evolve at similar rates, assuming no external irradiation, with disk mass loss rates of $\sim 10^{-8}$ M$_{\odot}$/yr.

preprint2022arXiv

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We present the results of the first successful ACA (Atacama Compact Array) [C I] $J$ = 1-0 mini-survey of seven protoplanetary disks. Using tailored azimuthally symmetric DALI (Dust And LInes) thermo-chemical disk models, supported by the [C I] $J$ = 1-0 and resolved CO isotopologue data, we determine the system-averaged elemental volatile carbon abundance in the outer disk of three sources. Six out of seven sources are detected in [C I] $J$ = 1-0 with ACA, four of which show a distinct disk component. Based on the modeling we find severe cold gaseous carbon depletion in the outer disk of DL Tau and moderate depletion in the outer disks of DR Tau and DO Tau. Combining the outer and inner disk carbon abundances, we demonstrate definitive evidence for radial drift in the disk of DL Tau, where the existence of multiple dust rings points to either short lived or leaky dust traps. We find dust locking in the compact and smooth disks of DO Tau and DR Tau, hinting at unresolved dust substructure. Comparing our results with stars of different ages and luminosities, we identify an observational evolutionary trend in gaseous carbon depletion that is consistent with dynamical models of CO depletion processes. Transport efficiency of solids in protoplanetary disks can significantly differ from what we expect based on the current resolved substructure in the continuum observations. This has important implications for our understanding of the impact of radial drift and pebble accretion on planetary compositions.

preprint2021arXiv

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.

preprint2021arXiv

Early planet formation in embedded protostellar disks: Setting the stage for the first generation of planetesimals

(Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally believed to be the product of the streaming instability. Hence the question can be restated: are the physical conditions of embedded disks conducive to the growth of the streaming instability? Here we model the collapse of a `dusty' proto-stellar cloud to show that if there is sufficient drift between the falling gas and dust, regions of the embedded disk can become sufficiently enhanced in dust to drive the streaming instability. We include four models, three with different dust grain sizes and one with a different initial cloud angular momentum to test a variety of collapse trajectories. We find a `sweet spot' for planetesimal formation for grain sizes of a few 10s of micron since they fall sufficiently fast relative to the gas to build a high dust-to-gas ratio along the disk midplane, but have slow enough radial drift speeds in the embedded disk to maintain the high dust-to-gas ratio. Unlike the gas, which is held in hydrostatic equilibrium for a time due to gas pressure, the dust can begin collapsing from all radii at a much earlier time. The streaming instability can produce at least between 7-35 M$_\oplus$ of planetesimals in the Class 0/I phase of our smooth embedded disks, depending on the size of the falling dust grains. This first generation of planetesimals could represent the first step in planet formation, and occurs earlier in the lifetime of the young star than is traditionally thought.

preprint2021arXiv

Probing inner and outer disk misalignments in transition disks

For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and $^{13}$CO molecular line data from ALMA to derive the orientation of the outer disk regions. We fit simple models to the GRAVITY data to derive the inner disks inclination and position angles. The outer disk geometries were derived from Keplerian fits to the ALMA velocity maps and compared to the inner disk constraints. We also predicted the locations of shadows for significantly misaligned systems. Our analysis reveals six disks to exhibit significant misalignments between their inner and outer disks. The predicted shadow positions agree well with the scattered light images of HD100453 and HD142527, and we find supporting evidence for a shadow in the disk around CQ Tau. In the other three targets for which we infer significantly misaligned disks, V1247 Ori, V1366 Ori, and RY Lup, we do not see any evident sign of shadows in the scattered light images. The scattered light shadows observed in DoAr44, HD135344B, and HD139614 are consistent with our observations, yet the underlying morphology is likely too complex to be described by our models and the accuracy achieved by our observations. Whereas we can derive precise constraints on the potential shadow positions for well-resolved inner disks around HAeBe stars, the statistical uncertainties for the marginally resolved inner disks around the TTS of our sample make it difficult to extract conclusive constraints for the presence of shadows in these systems.

preprint2021arXiv

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.

preprint2020arXiv

A cryogenic ice setup to simulate carbon atom reactions in interstellar ices

The design, implementation, and performance of a customized carbon atom beam source for the purpose of investigating solid-state reaction routes in interstellar ices in molecular clouds are discussed. The source is integrated into an existing ultrahigh vacuum setup, SURFace REaction SImulation DEvice (SURFRESIDE$^{2}$), which extends this double atom (H/D, O, and N) beamline apparatus with a third atom (C) beamline to a unique system that is fully suited to explore complex organic molecule solid-state formation under representative interstellar cloud conditions. The parameter space for this system is discussed, which includes the flux of the carbon atoms hitting the ice sample, their temperature, and the potential impact of temperature on ice reactions. Much effort has been put into constraining the beam size to within the limits of the sample size with the aim to reduce carbon pollution inside the setup. How the C-atom beam performs is quantitatively studied through the example experiment, C + $^{18}$O$_2$, and supported by computationally-derived activation barriers. The potential for this source to study the solid-state formation of interstellar complex organic molecules through C-atom reactions is discussed.

preprint2020arXiv

An experimental study of the surface formation of methane in interstellar molecular clouds

Methane is one of the simplest stable molecules that is both abundant and widely distributed across space. It is thought to have partial origin from interstellar molecular clouds, which are near the beginning of the star formation cycle. Observational surveys of CH$_4$ ice towards low- and high-mass young stellar objects showed that much of the CH$_4$ is expected to be formed by the hydrogenation of C on dust grains, and that CH$_4$ ice is strongly correlated with solid H$_2$O. Yet, this has not been investigated under controlled laboratory conditions, as carbon-atom chemistry of interstellar ice analogues has not been experimentally realized. In this study, we successfully demonstrate with a C-atom beam implemented in an ultrahigh vacuum apparatus the formation of CH$_4$ ice in two separate co-deposition experiments: C + H on a 10 K surface to mimic CH$_4$ formation right before H$_2$O ice is formed on the dust grain, and C + H + H$_2$O on a 10 K surface to mimic CH$_4$ formed simultaneously with H$_2$O ice. We confirm that CH$_4$ can be formed by the reaction of atomic C and H, and that the CH$_4$ formation rate is 2 times greater when CH$_4$ is formed within a H$_2$O-rich ice. This is in agreement with the observational finding that interstellar CH$_4$ and H$_2$O form together in the polar ice phase, i.e., when C- and H-atoms simultaneously accrete with O-atoms on dust grains. For the first time, the conditions that lead to interstellar CH$_4$ (and CD$_4$) ice formation are reported, and can be incorporated into astrochemical models to further constrain CH$_4$ chemistry in the interstellar medium and in other regions where CH$_4$ is inherited.

preprint2020arXiv

Chemically tracing the water snowline in protoplanetary disks with HCO$^+$

[Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earth's atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research is to investigate the validity of HCO$^+$ and H$^{13}$CO$^+$, as tracers of the water snowline in protoplanetary disks, as HCO$^+$ is destroyed by gas-phase water. Two small chemical networks are used to predict the HCO$^+$ abundance in a typical Herbig Ae disk. Subsequently, the corresponding emission profiles are modelled for H$^{13}$CO$^+$ and HCO$^+$ $J=2-1$, which provides the best balance between brightness and optical depth effects of the continuum emission. The HCO$^+$ abundance jumps by two orders of magnitude just outside the water snowline at 4.5 AU. We find that the emission of H$^{13}$CO$^+$ and HCO$^+$ is ring-shaped due to three effects: destruction of HCO$^+$ by gas-phase water, continuum optical depth, and molecular excitation effects. The presence of gas-phase water causes an additional drop of only $\sim$13% and 24% in the center of the disk, for H$^{13}$CO$^+$ and HCO$^+$, respectively. For the much more luminous outbursting source V883Ori, our models predict that the effect of dust and excitation are not limiting if the snowline is located outside $\sim$40 AU. Our analysis of ALMA observations of HCO$^+$ $J=3-2$ is consistent with the water snowline located around 100 AU. The HCO$^+$ abundance drops steeply around the water snowline, but dust and excitation can conceal the drop in HCO$^+$ emission due to the water snowline. Therefore, locating the water snowline with HCO$^+$ in Herbig disks is very difficult, but it is possible for outbursting sources like V883Ori.

preprint2020arXiv

Complex organic molecules in low-mass protostars on solar system scales -- I. Oxygen-bearing species

Complex organic molecules (COMs) are thought to form on icy dust grains in the earliest phase of star formation. The evolution of these COMs from the youngest Class 0/I protostellar phases toward the more evolved Class II phase is still not fully understood. Since planet formation seems to start early, and mature disks are too cold for characteristic COM emission lines, studying the inventory of COMs on solar system scales in the Class 0/I stage is relevant. ALMA Band 3 (3 mm) and Band 6 (1 mm) observations are obtained of seven Class 0 protostars in the Perseus and Serpens star-forming regions. By modeling the inner protostellar region using 'LTE' models, the excitation temperature and column densities are determined for several O-bearing COMs. B1-c, B1-bS, and Serpens S68N show COM emission, i.e, three out of the seven sources. No clear correlation seems to exist between the occurrence of COMs and source luminosity. The abundances of several COMs with respect to CH3OH are remarkably similar for the three COM-rich sources, and to IRAS 16293-2422B and HH 212. For other COMs the abundances differ by up to an order of magnitude, indicating that local source conditions are case determining. B1-c hosts a cold ($T_{ex}\approx60$ K), more extended component of COM emission with a column density of typically a few % of the warm/hot ($T_{ex}\sim 200$ K), central component. A D/H ratio of 1-3 % is derived based on the CH2DOH/CH3OH ratio suggesting a temperature of $\sim$15~K during the formation of methanol. This ratio is consistent with other low-mass protostars. Future mid-infrared facilities such as JWST/MIRI will be essential to directly observe COM ices. Combining this with a larger sample of COM-rich sources with ALMA will allow for directly linking ice and gas-phase abundances in order to constrain the routes that produce and maintain chemical complexity during the star formation process.

preprint2020arXiv

Formation of complex molecules in translucent clouds: Acetaldehyde, vinyl alcohol, ketene, and ethanol via nonenergetic processing of C2H2 ice

Complex organic molecules (COMs) have been identified toward high- and low-mass protostars as well as molecular clouds, suggesting that these interstellar species originate from the early stage(s) of starformation. The reaction pathways resulting in COMs described by the formula C$_2$H$_\text{n}$O are still under debate. In this work, we investigate the laboratory possible solid-state reactions that involve simple hydrocarbons and OH-radicals along with H$_2$O ice under translucent cloud conditions (1$\leq$A$_V$$\leq$5 and \textit{n}$_\text{H}$$\sim$10$^3$ cm$^{-3}$). We focus on the interactions of C$_2$H$_2$ with H-atoms and OH-radicals, which are produced along the H$_2$O formation sequence on grain surfaces at 10 K. Ultra-high vacuum (UHV) experiments were performed to study the surface chemistry observed during C$_2$H$_2$ + O$_2$ + H codeposition, where O$_2$ was used for the in-situ generation of OH-radicals. Reflection absorption infrared spectroscopy (RAIRS) was applied to in situ monitor the initial and newly formed species. After that, a temperature-programmed desorption experiment combined with a Quadrupole mass spectrometer (TPD-QMS) was used as a complementary analytical tool. The investigated 10 K surface chemistry of C$_2$H$_2$ with H-atoms and OH-radicals not only results in semi and fully saturated hydrocarbons, such as ethylene (C$_2$H$_4$) and ethane (C$_2$H$_6$), but it also leads to the formation of COMs, such as vinyl alcohol, acetaldehyde, ketene, ethanol, and possibly acetic acid. It is concluded that OH-radical addition reactions to C$_2$H$_2$, acting as a molecular backbone, followed by isomerization (i.e., keto-enol tautomerization) via an intermolecular pathway and successive hydrogenation provides a so far experimentally unreported solid-state route for the formation of these species without the need of energetic input.

preprint2020arXiv

HCN/HNC intensity ratio: a new chemical thermometer for the molecular ISM

Context. The gas kinetic temperature (TK) determines the physical and chemical evolution of the Interestellar Medium (ISM). However, obtaining reliable TK estimates usually requires expensive observations including the combination of multi-line analysis and dedicated radiative transfer calculations. Aims. This work explores the use of HCN and HNC observations, and particularly its I(HCN)/I(HNC) intensity ratio of their J=1-0 lines, as direct probe of the gas kinetic temperature in the molecular ISM. Methods. We obtained a new set of large-scale observations of both HCN and HNC (1-0) lines along the Integral Shape Filament (ISF) in Orion. In combination with ancillary gas and dust temperature measurements, we find a systematic temperature dependence of the observed I(HCN)/I(HNC) intensity ratio across our maps. Additional comparisons with chemical models demonstrate that these observed I(HCN)/I(HNC) variations are driven by the effective destruction and isomerization mechanisms of HNC under low energy barriers. Results. The observed variations of I(HCN)/I(HNC) with TK can be described with a two-part linear function. This empirical calibration is then used to create a temperature map of the entire ISF. Comparisons with similar dust temperature measurements in this cloud, as well as in other regions and galactic surveys, validate this simple technique to obtain direct estimates of the gas kinetic temperature in a wide range of physical conditions and scales with an optimal working range between 15 K and 40 K. Conclusions. Both observations and models demonstrate the strong sensitivity of the I(HCN)/I(HNC) ratio to the gas kinetic temperature. Since these lines are easily obtained in observations of local and extragalactic sources, our results highlight the potential use of this observable as new chemical thermometer for the ISM.

preprint2020arXiv

Molecule formation in dust-poor irradiated jets I. Stationary disk winds

Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations.

preprint2020arXiv

Observed sizes of planet-forming disks trace viscous evolution

The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3).

preprint2020arXiv

Probing UV-Sensitive Pathways for CN and HCN Formation in Protoplanetary Disks with the Hubble Space Telescope

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and sub-mm wavelengths, we present an analysis of new and archival HST, Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure LyA and FUV continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g. HCN, N2). Semi-forbidden C II] 2325 and UV-fluorescent H2 emission were also measured to constrain inner disk populations of C+ and vibrationally excited H2. We find a significant positive correlation between 14 micron HCN emission and fluxes from the FUV continuum and C II] 2325, consistent with model predictions requiring N2 photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between sub-mm CN emission and both C II] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and LyA or UV-H2 emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with JWST will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

preprint2020arXiv

The ALMA-PILS survey: First detection of the unsaturated 3-carbon molecules Propenal (C$_2$H$_3$CHO) and Propylene (C$_3$H$_6$) towards IRAS 16293$-$2422 B

Complex organic molecules with three carbon atoms are found in the earliest stages of star formation. In particular, propenal (C$_2$H$_3$CHO) is a species of interest due to its implication in the formation of more complex species and even biotic molecules. This study aims to search for the presence of C$_2$H$_3$CHO and other three-carbon species such as propylene (C$_3$H$_6$) in the hot corino region of the low-mass protostellar binary IRAS 16293--2422 to understand their formation pathways. We use ALMA observations in Band 6 and 7 from various surveys to search for the presence of C$_3$H$_6$ and C$_2$H$_3$CHO towards the protostar IRAS 16293--2422 B (IRAS 16293B). We report the detection of both C$_3$H$_6$ and C$_2$H$_3$CHO towards IRAS 16293B, however, no unblended lines were found towards the other component of the binary system, IRAS 16293A. We derive column density upper limits for C$_3$H$_8$, HCCCHO, n-C$_3$H$_7$OH, i-C$_3$H$_7$OH, C$_3$O, and cis-HC(O)CHO towards IRAS 16293B. We then use a three-phase chemical model to simulate the formation of these species in a typical prestellar environment followed by its hydrodynamical collapse until the birth of the central protostar. Different formation paths, such as successive hydrogenation and radical-radical additions on grain surfaces, are tested and compared to the observational results. The simulations reproduce the abundances within one order of magnitude from those observed towards IRAS 16293B, with the best agreement found for a rate of $10^{-12}$ cm$^3$ s$^{-1}$ for the gas-phase reaction C$_3$ + O $\rightarrow$ C$_2$ + CO. Successive hydrogenations of C$_3$, HC(O)CHO, and CH$_3$OCHO on grain surfaces are a major and crucial formation route of complex organics molecules, whereas both successive hydrogenation pathways and radical-radical addition reactions contribute to the formation of C$_2$H$_5$CHO.

preprint2020arXiv

The ALMA-PILS survey: Inventory of complex organic molecules towards IRAS 16293-2422 A

Complex organic molecules (COM) are detected in many sources in the warm inner regions of envelopes surrounding deeply embedded protostars. Exactly how these COM form remains an open question. This study aims to constrain the formation of complex organic molecules through comparisons of their abundances towards the Class 0 protostellar binary IRAS 16293-2422 (IRAS16293). We utilised observations from the ALMA Protostellar Interferometric Line Survey of IRAS16293. The species identification and the rotational temperature and column density estimation were derived by fitting the extracted spectra towards IRAS16293 A and IRAS16293 B with synthetic spectra. The majority of the work in this paper pertains to the analysis of IRAS16293 A for a comparison with the results from the other binary component, which have already been published. We detect 15 different COM, as well as 16 isotopologues towards the most luminous companion protostar IRAS16293 A. Tentative detections of an additional 11 isotopologues are reported. We also searched for and report on the first detections of CH3OCH2OH and t-C2H5OCH3 towards IRAS16293 B and the follow-up detection of CH2DCHO and CH3CDO. Twenty-four lines of CHD2OH are also identified. The comparison between the two protostars of the binary system shows significant differences in abundance for some of the species, which are partially correlated to their spatial distribution. The spatial distribution is consistent with the sublimation temperature of the species; those with higher expected sublimation temperatures are located in the most compact region of the hot corino towards IRAS16293 A. This spatial differentiation is not resolved in IRAS16293 B and will require observations at a higher angular resolution. In parallel, the list of identified CHD2OH lines shows the need of accurate spectroscopic data including their line strength.

preprint2020arXiv

The GRAVITY young stellar object survey. II. First spatially resolved observations of the CO bandhead emission in a high-mass YSO

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $μ$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

preprint2019arXiv

Mass constraints for 15 protoplanetary disks from HD 1-0

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $Δ_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.

preprint2019arXiv

Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $β$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $β$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $β$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

preprint2016arXiv

Constraining the physical structure of the inner few 100 AU scales of deeply embedded low-mass protostars

(Abridged) The physical structure of deeply-embedded low-mass protostars (Class 0) on scales of less than 300 AU is still poorly constrained. Determining this is crucial for understanding the physical and chemical evolution from cores to disks. In this study two models of the emission, a Gaussian disk intensity distribution and a parametrized power-law disk model, are fitted to sub-arcsecond resolution interferometric continuum observations of five Class 0 sources, including one source with a confirmed Keplerian disk. For reference, a spherically symmetric single power-law envelope is fitted to the larger scale ($\sim$1000 AU) emission and investigated further for one of the sources on smaller scales. A thin disk model can approximate the emission and physical structure in the inner few 100 AU scales of the studied deeply-embedded low-mass protostars and paves the way for analysis of a larger sample with ALMA. While the disk radii agree with previous estimates the masses are different for some of the sources studied. Assuming a typical temperature distribution, the fractional amount of mass in the disk above 100 K varies in between 7% to 30%. Kinematic data are needed to determine the presence of any Keplerian disk. Using previous observations of p-H$_2^{18}$O, we estimate the relative gas phase water abundances roughly an order of magnitude higher than previously inferred when both warm and cold H$_2$ was used as reference. A spherically symmetric single power-law envelope model fails to simultaneously reproduce both the small and large scale emission.

preprint2016arXiv

Evidence for a correlation between mass accretion rates onto young stars and the mass of their protoplanetary disks

A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/X-Shooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium (ISM) gas-to-dust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mm-sized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.

preprint2016arXiv

H-atom addition and abstraction reactions in mixed CO, H2CO and CH3OH ices: an extended view on complex organic molecule formation

Complex organic molecules (COMs) have been observed not only in the hot cores surrounding low- and high- mass protostars, but also in cold dark clouds. Therefore, it is interesting to understand how such species can be formed without the presence of embedded energy sources. We present new laboratory experiments on the low-temperature solid state formation of three complex molecules: methyl formate (HC(O)OCH3), glycolaldehyde (HC(O)CH2OH) and ethylene glycol (H2C(OH)CH2OH), through recombination of free radicals formed via H-atom addition and abstraction reactions at different stages in the CO-H2CO-CH3OH hydrogenation network at 15 K. The experiments extend previous CO hydrogenation studies and aim at resembling the physical&chemical conditions typical of the CO freeze-out stage in dark molecular clouds, when H2CO and CH3OH form by recombination of accreting CO molecules and H-atoms on ice grains. We confirm that H2CO, once formed through CO hydrogenation, not only yields CH3OH through ongoing H-atom addition reactions, but is also subject to H-atom-induced abstraction reactions, yielding CO again. In a similar way, H2CO is also formed in abstraction reactions involving CH3OH. The dominant methanol H-atom abstraction product is expected to be CH2OH, while H-atom additions to H2CO should at least partially proceed through CH3O intermediate radicals. The occurrence of H-atom abstraction reactions in ice mantles leads to more reactive intermediates (HCO, CH3O and CH2OH) than previously thought, when assuming sequential H-atom addition reactions only. This enhances the probability to form COMs through radical-radical recombination without the need of UV photolysis or cosmic rays as external triggers.

preprint2016arXiv

Observations and modelling of CO and [CI] in disks. First detections of [CI] and constraints on the carbon abundance

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

preprint2016arXiv

Probing the 2D temperature structure of protoplanetary disks with Herschel observations of high-J CO lines

The gas temperature structure of protoplanetary disks is a key ingredient for interpreting various disk observations and for quantifying the subsequent evolution of these systems. The comparison of low- and mid-$J$ CO rotational lines is a powerful tool to assess the temperature gradient in the warm molecular layer of disks. Spectrally resolved high-$J$ ($J_{\rm u} > 14$) CO lines probe intermediate distances and heights from the star that are not sampled by (sub-)millimeter CO spectroscopy. This paper presents new {\it Herschel}/HIFI and archival PACS observations of $^{12}$CO, $^{13}$CO and \cii \ emission in 4 Herbig AeBe (HD 100546, HD 97048, IRS 48, HD 163296) and 3 T Tauri (AS 205, S CrA, TW Hya) disks. In the case of the T Tauri systems AS 205 and S CrA, the CO emission has a single-peaked profile, likely due to a slow wind. For all other systems, the {\it Herschel} CO spectra are consistent with pure disk emission and the spectrally-resolved lines (HIFI) and the CO rotational ladder (PACS) are analyzed simultaneously assuming power-law temperature and column density profiles, using the velocity profile to locate the emission in the disk. The temperature profile varies substantially from disk to disk. In particular, $T_{\rm gas}$ in the disk surface layers can differ by up to an order of magnitude among the 4 Herbig AeBe systems with HD 100546 being the hottest and HD 163296 the coldest disk of the sample. Clear evidence of a warm disk layer where $T_{\rm gas} > T_{\rm dust}$ is found in all the Herbig Ae disks. The observed CO fluxes and line profiles are compared to predictions of physical-chemical models. The primary parameters affecting the disk temperature structure are the flaring angle, the gas-to-dust mass ratio the scale height and the dust settling.

preprint2016arXiv

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

preprint2016arXiv

The ALMA Protostellar Interferometric Line Survey (PILS): First results from an unbiased submillimeter wavelength line survey of the Class 0 protostellar binary IRAS 16293-2422 with ALMA

The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array (ALMA). The survey covers the full frequency range from 329 to 363 GHz (0.8 mm) with additional targeted observations at 3.0 and 1.3 mm. More than 10,000 features are detected toward one component in the protostellar binary. Glycolaldehyde, its isomers, methyl formate and acetic acid, and its reduced alcohol, ethylene glycol, are clearly detected. For ethylene glycol both lowest state conformers, aGg' and gGg', are detected, the latter for the first time in the ISM. The abundance of glycolaldehyde is comparable to or slightly larger than that of ethylene glycol. In comparison to the Galactic Center, these two species are over-abundant relative to methanol, possibly an indication of formation at low temperatures in CO-rich ices. Both 13C and deuterated isotopologues of glycolaldehyde are detected, also for the first time ever in the ISM. For the deuterated species, a D/H ratio of approximately 5% is found with no differences between the deuteration in the different functional groups of glycolaldehyde. Measurements of the 13C-species lead to a 12C:13C ratio of approximately 30, lower than the typical ISM value. This low ratio may reflect an enhancement of 13CO in the ice due to either ion-molecule reactions in the gas before freeze-out or differences in the temperatures where 12CO and 13CO ices sublimate. The results reinforce the importance of low-temperature grain surface chemistry for the formation of prebiotic molecules seen here in the gas after sublimation of the entire ice mantle.

preprint2016arXiv

The Herschel-HIFI view of mid-IR quiet massive protostellar objects

We present Herschel/HIFI observations (WISH KP) of 14 water lines in a small sample of galactic massive protostellar objects: NGC6334I(N), DR21(OH), IRAS16272-4837, and IRAS05358+3543. We analyze the gas dynamics from the line profiles. Through modeling of the observations using RATRAN, we estimate outflow, infall, turbulent velocities, molecular abundances, and investigate any correlation with the evolutionary status of each source. The molecular line profiles exhibit a broad component coming from the shocks along the cavity walls associated with the protostars, and an infalling (or expansion for IRAS05358+3543) and passively heated envelope component, with highly supersonic turbulence likely increasing with the distance from the center. Accretion rates between 6.3 10^{-5} and 5.6 10^{-4} \msun yr^{-1} are derived from the infall observed in three of our sources. The outer water abundance is estimated to be at the typical value of a few 10^{-8} while the inner abundance varies from 1.7 10^{-6} to 1.4 10^{-4} with respect to H2 depending on the source. We confirm that regions of massive star formation are highly turbulent and that the turbulence likely increases in the envelope with the distance to the star. The inner abundances are lower than the expected 10^{-4} perhaps because our observed lines do not probe deep enough into the inner envelope, or because photodissociation through protostellar UV photons is more efficient than expected. We show that the higher the infall/expansion velocity in the protostellar envelope, the higher is the inner abundance, maybe indicating that larger infall/expansion velocities generate shocks that will sputter water from the ice mantles of dust grains in the inner region. High-velocity water must be formed in the gas-phase from shocked material.

preprint2016arXiv

Volatile carbon locking and release in protoplanetary disks. A study of TW Hya and HD 100546

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $λ\,$Boötis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

preprint2016arXiv

Water around IRAS15398-3359 observed with ALMA

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

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

[OI]63micron jets in class 0 sources detected by Herschel

We present Herschel PACS mapping observations of the [OI]63 micron line towards protostellar outflows in the L1448, NGC1333-IRAS4, HH46, BHR71 and VLA1623 star forming regions. We detect emission spatially resolved along the outflow direction, which can be associated with a low excitation atomic jet. In the L1448-C, HH46 IRS and BHR71 IRS1 outflows this emission is kinematically resolved into blue- and red-shifted jet lobes, having radial velocities up to 200 km/s. In the L1448-C atomic jet the velocity increases with the distance from the protostar, similarly to what observed in the SiO jet associated with this source. This suggests that [OI] and molecular gas are kinematically connected and that this latter could represent the colder cocoon of a jet at higher excitation. Mass flux rates (M$_{jet}$(OI)) have been measured from the [OI]63micron luminosity adopting two independent methods. We find values in the range 1-4 10$^{-7}$ Mo/yr for all sources but HH46, for which an order of magnitude higher value is estimated. M$_{jet}$(OI) are compared with mass accretion rates (M$_{acc}$) onto the protostar and with M$_{jet}$ derived from ground-based CO observations. M$_{jet}$(OI)/M$_{acc}$ ratios are in the range 0.05-0.5, similar to the values for more evolved sources. M$_{jet}$(OI) in HH46 IRS and IRAS4A are comparable to M$_{jet}$(CO), while those of the remaining sources are significantly lower than the corresponding M$_{jet}$(CO). We speculate that for these three sources most of the mass flux is carried out by a molecular jet, while the warm atomic gas does not significantly contribute to the dynamics of the system.

preprint2015arXiv

Classifying the embedded young stellar population in Perseus and Taurus & the LOMASS database

Context. The classification of young stellar objects (YSOs) is typically done using the infrared spectral slope or bolometric temperature, but either can result in contamination of samples. More accurate methods to determine the evolutionary stage of YSOs will improve the reliability of statistics for the embedded YSO population and provide more robust stage lifetimes. Aims. We aim to separate the truly embedded YSOs from more evolved sources. Methods. Maps of HCO+ J=4-3 and C18O J=3-2 were observed with HARP on the James Clerk Maxwell Telescope (JCMT) for a sample of 56 candidate YSOs in Perseus and Taurus in order to characterize emission from high (column) density gas. These are supplemented with archival dust continuum maps observed with SCUBA on the JCMT and Herschel PACS to compare the morphology of the gas and dust in the protostellar envelopes. The spatial concentration of HCO+ J=4-3 and 850 micron dust emission are used to classify the embedded nature of YSOs. Results. Approximately 30% of Class 0+I sources in Perseus and Taurus are not Stage I, but are likely to be more evolved Stage II pre-main sequence (PMS) stars with disks. An additional 16% are confused sources with an uncertain evolutionary stage. Conclusions. Separating classifications by cloud reveals that a high percentage of the Class 0+I sources in the Perseus star forming region are truly embedded Stage I sources (71%), while the Taurus cloud hosts a majority of evolved PMS stars with disks (68%). The concentration factor method is useful to correct misidentified embedded YSOs, yielding higher accuracy for YSO population statistics and Stage timescales. Current estimates (0.54 Myr) may overpredict the Stage I lifetime on the order of 30%, resulting in timescales of 0.38 Myr for the embedded phase.

preprint2015arXiv

Disentangling the jet emission from protostellar systems. The ALMA view of VLA1623

Context: High-resolution studies of class 0 protostars represent the key to constraining protostar formation models. VLA16234-2417 represents the prototype of class 0 protostars, and it has been recently identified as a triple non-coeval system. Aim: We aim at deriving the physical properties of the jets in VLA16234-2417 using tracers of shocked gas. Methods: ALMA Cycle 0 Early Science observations of CO(2-1) in the extended configuration are presented in comparison with previous SMA CO(3-2) and Herschel-PACS [OI}] 63 micron observations. Gas morphology and kinematics were analysed to constrain the physical structure and origin of the protostellar outflows. Results: We reveal a collimated jet component associated with the [OI] 63 micron emission at about 8'' (about 960 AU) from source B. This newly detected jet component is inversely oriented with respect to the large-scale outflow driven by source A, and it is aligned with compact and fast jet emission very close to source B (about 0.3'') rather than with the direction perpendicular to the A disk. We also detect a cavity-like structure at low projected velocities, which surrounds the [OI] 63 micron emission and is possibly associated with the outflow driven by source A. Finally, no compact outflow emission is associated with source W. Conclusions: Our high-resolution ALMA observations seem to suggest there is a fast and collimated jet component associated with source B. This scenario would confirm that source B is younger than A, that it is in a very early stage of evolution, and that it drives a faster, more collimated, and more compact jet with respect to the large-scale slower outflow driven by A. However, a different scenario of a precessing jet driven by A cannot be firmly excluded from the present observations.

preprint2015arXiv

Outflow forces in intermediate mass star formation

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

preprint2015arXiv

Quantifying the gas inside dust cavities in transitional disks: implications for young planets

ALMA observations of a small sample of transitional disks with large dust cavities observed in Cycle 0 and 1 are summarized. The gas and dust surface density structures are inferred from the continuum and 12CO, 13CO and C18O line data using the DALI physical-chemical code. Thanks to its ability to self-shield, CO can survive inside dust cavities in spite of being exposed to intense UV radiation and can thus be used as a probe of the gas structure. Modeling of the existing data shows that gas is present inside the dust cavities in all cases, but at a reduced level compared with the gas surface density profile of the outer disk. The gas density decrease inside the dust cavity radius by factors of up to 10^4 suggests clearing by one or more planetary-mass companions. The accompanying pressure bumps naturally lead to trapping of the mm-sized dust grains observed in the ALMA images.

preprint2015arXiv

Reconstructing the history of water ice formation from HDO/H2O and D2O/HDO ratios in protostellar cores

Recent interferometer observations have found that the D2O/HDO abundance ratio is higher than that of HDO/H2O by about one order of magnitude in the vicinity of low-mass protostar NGC 1333-IRAS 2A, where water ice has sublimated. Previous laboratory and theoretical studies show that the D2O/HDO ice ratio should be lower than the HDO/H2O ice ratio, if HDO and D2O ices are formed simultaneously with H2O ice. In this work, we propose that the observed feature, D2O/HDO > HDO/H2O, is a natural consequence of chemical evolution in the early cold stages of low-mass star formation: 1) majority of oxygen is locked up in water ice and other molecules in molecular clouds, where water deuteration is not efficient, and 2) water ice formation continues with much reduced efficiency in cold prestellar/protostellar cores, where deuteration processes are highly enhanced due to the drop of the ortho-para ratio of H2, the weaker UV radiation field, etc. Using a simple analytical model and gas-ice astrochemical simulations tracing the evolution from the formation of molecular clouds to protostellar cores, we show that the proposed scenario can quantitatively explain the observed HDO/H2O and D2O/HDO ratios. We also find that the majority of HDO and D2O ices are likely formed in cold prestellar/protostellar cores rather than in molecular clouds, where the majority of H2O ice is formed. This work demonstrates the power of the combination of the HDO/H2O and D2O/HDO ratios as a tool to reveal the past history of water ice formation in the early cold stages of star formation and when the enrichment of deuterium in the bulk of water occurred. Further observations are needed to explore if the relation, D2O/HDO > HDO/H2O, is common in low-mass protostellar sources.

preprint2015arXiv

Search for methylamine in high mass hot cores

We aim to detect methylamine, CH$_{3}$NH$_{2}$, in a variety of hot cores and use it as a test for the importance of photon-induced chemistry in ice mantles and mobility of radicals. Specifically, CH$_3$NH$_2$ cannot be formed from atom addition to CO whereas other NH$_2$-containing molecules such as formamide, NH$_2$CHO, can. Submillimeter spectra of several massive hot core regions were taken with the James Clerk Maxwell Telescope. Abundances are determined with the rotational diagram method where possible. Methylamine is not detected, giving upper limit column densities between 1.9 $-$ 6.4 $\times$ 10$^{16}$ cm$^{-2}$ for source sizes corresponding to the 100 K envelope radius. Combined with previously obtained JCMT data analyzed in the same way, abundance ratios of CH$_{3}$NH$_{2}$, NH$_{2}$CHO and CH$_{3}$CN with respect to each other and to CH$_{3}$OH are determined. These ratios are compared with Sagittarius B2 observations, where all species are detected, and to hot core models. The observed ratios suggest that both methylamine and formamide are overproduced by up to an order of magnitude in hot core models. Acetonitrile is however underproduced. The proposed chemical schemes leading to these molecules are discussed and reactions that need further laboratory studies are identified. The upper limits obtained in this paper can be used to guide future observations, especially with ALMA.

preprint2015arXiv

Testing particle trapping in transition disks with ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336GHz of two transition disks, SR21 and HD135344B. In combination with previous ALMA observations from Cycle 0 at 689GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($α_{\rm{mm}}$). The observations of SR21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR21 with values of $α_{\rm{mm}}{\sim}3.8-4.2$ in the inner region ($r<35$ AU) and $α_{\rm{mm}}{\sim}2.6-3.0$ outside. An axisymmetric ring (which we call the ring model) or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex model), is assumed for fitting the disk morphology. For SR21, the ring model better fits the emission at 336GHz, conversely the vortex model better fits the 689GHz emission. For HD135344B, neither a significant shift in the null of the visibilities nor radial variations of $α_{\rm{mm}}$ are detected. Furthermore, for HD135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $α_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD135344B, although it cannot be excluded either.

preprint2015arXiv

The ALMA Band 9 receiver - Design, construction, characterization, and first light

We describe the design, construction, and characterization of the Band 9 heterodyne receivers (600-720 GHz) for the Atacama Large Millimeter / submillimeter Array (ALMA). The ALMA Band 9 receiver units ("cartridges"), which are installed in the telescope's front end, have been designed to detect and down-convert two orthogonal linear polarization components of the light collected by the ALMA antennas. The light entering the front end is refocused with a compact arrangement of mirrors, which is fully contained within the cartridge. The arrangement contains a grid to separate the polarizations and two beam splitters to combine each resulting beam with a local oscillator signal. The combined beams are fed into independent double-sideband mixers, each with a corrugated feedhorn coupling the radiation by way of a waveguide with backshort cavity into an impedance-tuned SIS junction that performs the heterodyne down-conversion. Finally, the generated intermediate frequency signals are amplified by cryogenic and room-temperature HEMT amplifiers and exported to the telescope's back end for further processing and, finally, correlation. The receivers have been constructed and tested in the laboratory and they show excellent performance, complying with ALMA requirements. Performance statistics on all 73 Band 9 receivers are reported. On-sky characterization and tests of the performance of the Band 9 cartridges are presented using commissioning data.

preprint2015arXiv

The Mid-Infrared Instrument for JWST, II: Design and Build

The Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) provides measurements over the wavelength range 5 to 28.5 microns. MIRI has, within a single 'package', four key scientific functions: photometric imaging, coronagraphy, single-source low-spectral resolving power (R ~ 100) spectroscopy, and medium-resolving power (R ~ 1500 to 3500) integral field spectroscopy. An associated cooler system maintains MIRI at its operating temperature of < 6.7 K. This paper describes the driving principles behind the design of MIRI, the primary design parameters, and their realization in terms of the 'as-built' instrument. It also describes the test program that led to delivery of the tested and calibrated Flight Model to NASA in 2012, and the confirmation after delivery of the key interface requirements.

preprint2015arXiv

The Mid-Infrared Instrument for the James Webb Space Telescope, VI: The Medium Resolution Spectrometer

We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral field unit (IFU): an 'image slicer' that reformats the input field for presentation to a grating spectrometer. Two 1024 X 1024 focal plane arrays record the output spectral images with an instantaneous spectral coverage of approximately one third of the full wavelength range of each channel. The full 5 to 28.5 micron spectrum is then obtained by making three exposures using gratings and pass-band-determining filters that are selected using just two three-position mechanisms. The expected on-orbit optical performance is presented, based on testing of the MIRI Flight Model and including spectral and spatial coverage and resolution. The point spread function of the reconstructed images is shown to be diffraction limited and the optical transmission is shown to be consistent with the design expectations.

preprint2015arXiv

Water in star-forming regions with Herschel (WISH) V. The physical conditions in low-mass protostellar outflows revealed by multi-transition water observations

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I.

preprint2014arXiv

First spectrally-resolved H$_2$ observations towards HH 54 / Low H$_2$O abundance in shocks

Context: Herschel observations suggest that the H$_2$O distribution in outflows from low-mass stars resembles the H$_2$ emission. It is still unclear which of the different excitation components that characterise the mid- and near-IR H$_2$ distribution is associated with H$_2$O. Aim: The aim is to spectrally resolve the different excitation components observed in the H$_2$ emission. This will allow us to identify the H$_2$ counterpart associated with H$_2$O and finally derive directly an H$_2$O abundance estimate with respect to H$_2$. Methods: We present new high spectral resolution observations of H$_2$ 0-0 S(4), 0-0 S(9), and 1-0 S(1) towards HH 54, a bright nearby shock region in the southern sky. In addition, new Herschel-HIFI H$_2$O (2$_{12}$$-$1$_{01}$) observations at 1670~GHz are presented. Results: Our observations show for the first time a clear separation in velocity of the different H$_2$ lines: the 0-0 S(4) line at the lowest excitation peaks at $-$7~km~s$^{-1}$, while the more excited 0-0 S(9) and 1-0 S(1) lines peak at $-$15~km~s$^{-1}$. H$_2$O and high-$J$ CO appear to be associated with the H$_2$ 0-0 S(4) emission, which traces a gas component with a temperature of 700$-$1000 K. The H$_2$O abundance with respect to H$_2$ 0-0 S(4) is estimated to be $X$(H$_2$O)$<$1.4$\times$10$^{-5}$ in the shocked gas over an area of 13$^{\prime\prime}$. Conclusions: We resolve two distinct gas components associated with the HH 54 shock region at different velocities and excitations. This allows us to constrain the temperature of the H$_2$O emitting gas ($\leq$1000 K) and to derive correct estimates of H$_2$O abundance in the shocked gas, which is lower than what is expected from shock model predictions.

preprint2014arXiv

Gas and dust structures in protoplanetary disks hosting multiple planets

Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. We study the gas and dust distributions of a disk hosting two massive planets as function of different disk and dust parameters. Observational signatures, such as the spectral energy distribution, sub-millimetre, and polarised images are simulated for the various parameters. We confirm that planets can lead to particle trapping, although for a disk with high viscosity ($α_{\rm{turb}}=10^{-2}$), the planet should be more massive than $5 M_{\rm{Jup}}$ and dust fragmentation should occur with low efficiency ($v_{f}\sim30\rm{m s}^{-1}$). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity ($α_{\rm{turb}}=10^{-2}$) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity ($α_{\rm{turb}}=10^{-3}$), the planet mass is constrained to be less than $5 M_{\rm{Jup}}$. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.

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

Observations of water with Herschel/HIFI toward the high-mass protostar AFGL 2591

Water is a sensitive tracer of physical conditions in star-forming regions because of its large abundance variations between hot and cold regions. We use spectrally resolved observations of rotational lines of H$_2$O and its isotopologs with Herschel/HIFI to constrain the physical conditions of the water emitting region toward the high-mass protostar AFGL2591. We use analytical estimates and rotation diagrams to estimate T$_{ex}$ and column densities of H$_2$O of the envelope, the outflow, and a foreground cloud. Furthermore, we use the non-LTE radiative transfer code to estimate the temperature and volume density of the H$_2$O emitting gas. Assuming LTE, we estimate an T$_{ex}$ of 42 K and a column density of 2$\times$10$^{14}$ cm$^{-2}$ for the envelope and 45 K and 4$\times$10$^{13}$ cm$^{-2}$ for the outflow, in beams of 4" and 30", respectively. Non-LTE models indicate a kinetic temperature of 60-230 K and a volume density of 7$\times$10$^6$-10$^8$ cm$^{-3}$ for the envelope, and a kinetic temperature of 70-90 K and a gas density of 10$^7$-10$^8$ cm$^{-3}$ for the outflow. The o/p ratio of the foreground absorption is 1.9$\pm$0.4, suggesting a low temperature. In contrast, the o/p ratio seen in absorption by the outflow is 3.5$\pm$1.0, as expected for warm gas. The water abundance in the envelope is 10$^{-9}$, similar to the low values found for other high- and low-mass protostars, suggesting that this abundance is constant during the embedded phase of high-mass star formation. The water abundance in the outflow is 10$^{-10}$, which is 10$\times$ lower than in the envelope and in the outflows of high- and low-mass protostars. Since beam size effects can only increase this estimate by a factor of 2, we suggest that the water in the outflow is affected by dissociating UV radiation due to the low extinction in the outflow lobe.

preprint2014arXiv

Shockingly low water abundances in Herschel / PACS observations of low-mass protostars in Perseus

Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.

preprint2014arXiv

The deuterium fractionation of water on solar-system scales in deeply-embedded low-mass protostars

(Abridged) The water deuterium fractionation (HDO/H$_2$O abundance ratio) has traditionally been used to infer the amount of water brought to Earth by comets. Measuring this ratio in deeply-embedded low-mass protostars makes it possible to probe the critical stage when water is transported from clouds to disks in which icy bodies are formed. We present sub-arcsecond resolution observations of HDO in combination with H$_2^{18}$O from the PdBI toward the three low-mass protostars NGC 1333-IRAS 2A, IRAS 4A-NW, and IRAS 4B. The resulting HDO/H$_2$O ratio is $7.4\pm2.1\times10^{-4}$ for IRAS 2A, $19.1\pm5.4\times10^{-4}$ for IRAS 4A-NW, and $5.9\pm1.7\times10^{-4}$ for IRAS 4B. Derived ratios agree with radiative transfer models within a factor of 2-4 depending on the source. Our HDO/H$_2$O ratios for the inner regions (where $T>100$ K) of four young protostars are only a factor of 2 higher than those found for pristine, solar system comets. These small differences suggest that little processing of water occurs between the deeply embedded stage and the formation of planetesimals and comets.

preprint2014arXiv

Water distribution in shocked regions of the NGC1333-IRAS4A protostellar outflow

We present the study of the H2O spatial distribution at two bright shocked regions along IRAS4A, one of the strongest H2O emitters among the Class 0 outflows. We obtained Herschel-PACS maps of the IRAS4A outflow and HIFI observations of two shocked positions. The largest HIFI beam of 38 arcsec at 557 GHz was mapped in several key water lines with different upper energy levels, to reveal possible spatial variations of the line profiles. We detect four H2O lines and CO (16-15) at the two selected positions. In addition, transitions from related outflow and envelope tracers are detected. Different gas components associated with the shock are identified in the H2O emission. In particular, at the head of the red lobe of the outflow, two distinct gas components with different excitation conditions are distinguished in the HIFI emission maps: a compact component, detected in the ground-state water lines, and a more extended one. Assuming that these two components correspond to two different temperature components observed in previous H2O and CO studies, the excitation analysis of the H2O emission suggests that the compact (about 3 arcsec) component is associated with a hot (T~1000 K) gas with densities ~(1-4)x10^5 cm^{-3}, whereas the extended one (10-17 arcsec) traces a warm (T~300-500 K) and dense gas (~(3-5)x10^7 cm^{-3}). Finally, using the CO (16-15) emission observed at R2, we estimate the H2O/H2 abundance of the warm and hot components to be (7-10)x10^{-7} and (3-7)x10^{-5}. Our data allowed us, for the first time, to resolve spatially the two temperature components previously observed with HIFI and PACS. We propose that the compact hot component may be associated with the jet that impacts the surrounding material, whereas the warm, dense, and extended component originates from the compression of the ambient gas by the propagating flow.

preprint2014arXiv

Water in Low-Mass Star-Forming Regions with Herschel: The Link Between Water Gas and Ice in Protostellar Envelopes

Aims: Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods: A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., FUV field, initial abundances, timescales, and kinematics). Results: The observed water ice abundances with respect to hydrogen nuclei in our sample are 30-80ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T(dust)>15 K. This requires short prestellar core lifetimes of less than about 0.1Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions: The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short pre-collapse time of less than 0.1 Myr at densities of n(H)~1e4 cm-3, or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.

E. F. van Dishoeck

What is connected

Connect this record

See the researcher in context

Building this map preview

129 published item(s)

Characterizing the physical and chemical properties of the Class I protostellar system Oph-IRS 44. Binarity, infalling streamers, and accretion shocks

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

Gas temperature structure across transition disk cavities

H$_2$O distribution in the disc of HD 100546 and HD 163296: the role of dust dynamics and planet--disc interaction

Importance of source structure on complex organics emission II. Can disks explain lack of methanol emission from some low-mass protostars?

Importance of source structure on complex organics emission. I. Observations of CH$_3$OH from low-mass to high-mass protostars

Self-consistent modelling of the dust component in protoplanetary and circumplanetary disks: the case of PDS 70

Survey of Orion Disks with ALMA (SODA) I: Cloud-level demographics of 873 protoplanetary disks

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

Early planet formation in embedded protostellar disks: Setting the stage for the first generation of planetesimals

Probing inner and outer disk misalignments in transition disks

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

A cryogenic ice setup to simulate carbon atom reactions in interstellar ices

An experimental study of the surface formation of methane in interstellar molecular clouds

Chemically tracing the water snowline in protoplanetary disks with HCO$^+$

Complex organic molecules in low-mass protostars on solar system scales -- I. Oxygen-bearing species

Formation of complex molecules in translucent clouds: Acetaldehyde, vinyl alcohol, ketene, and ethanol via nonenergetic processing of C2H2 ice

HCN/HNC intensity ratio: a new chemical thermometer for the molecular ISM

Molecule formation in dust-poor irradiated jets I. Stationary disk winds

Observed sizes of planet-forming disks trace viscous evolution

Probing UV-Sensitive Pathways for CN and HCN Formation in Protoplanetary Disks with the Hubble Space Telescope

The ALMA-PILS survey: First detection of the unsaturated 3-carbon molecules Propenal (C$_2$H$_3$CHO) and Propylene (C$_3$H$_6$) towards IRAS 16293$-$2422 B

The ALMA-PILS survey: Inventory of complex organic molecules towards IRAS 16293-2422 A

The GRAVITY young stellar object survey. II. First spatially resolved observations of the CO bandhead emission in a high-mass YSO

Mass constraints for 15 protoplanetary disks from HD 1-0

Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry

Constraining the physical structure of the inner few 100 AU scales of deeply embedded low-mass protostars

Evidence for a correlation between mass accretion rates onto young stars and the mass of their protoplanetary disks

H-atom addition and abstraction reactions in mixed CO, H2CO and CH3OH ices: an extended view on complex organic molecule formation

Observations and modelling of CO and [CI] in disks. First detections of [CI] and constraints on the carbon abundance

Probing the 2D temperature structure of protoplanetary disks with Herschel observations of high-J CO lines

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

The ALMA Protostellar Interferometric Line Survey (PILS): First results from an unbiased submillimeter wavelength line survey of the Class 0 protostellar binary IRAS 16293-2422 with ALMA

The Herschel-HIFI view of mid-IR quiet massive protostellar objects

Volatile carbon locking and release in protoplanetary disks. A study of TW Hya and HD 100546

Water around IRAS15398-3359 observed with ALMA

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

[OI]63micron jets in class 0 sources detected by Herschel

Classifying the embedded young stellar population in Perseus and Taurus & the LOMASS database

Disentangling the jet emission from protostellar systems. The ALMA view of VLA1623

Outflow forces in intermediate mass star formation

Quantifying the gas inside dust cavities in transitional disks: implications for young planets

Reconstructing the history of water ice formation from HDO/H2O and D2O/HDO ratios in protostellar cores

Search for methylamine in high mass hot cores

Testing particle trapping in transition disks with ALMA

The ALMA Band 9 receiver - Design, construction, characterization, and first light

The Mid-Infrared Instrument for JWST, II: Design and Build

The Mid-Infrared Instrument for the James Webb Space Telescope, VI: The Medium Resolution Spectrometer

Water in star-forming regions with Herschel (WISH) V. The physical conditions in low-mass protostellar outflows revealed by multi-transition water observations

First spectrally-resolved H$_2$ observations towards HH 54 / Low H$_2$O abundance in shocks

Gas and dust structures in protoplanetary disks hosting multiple planets

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

Observations of water with Herschel/HIFI toward the high-mass protostar AFGL 2591

Shockingly low water abundances in Herschel / PACS observations of low-mass protostars in Perseus

The deuterium fractionation of water on solar-system scales in deeply-embedded low-mass protostars

Water distribution in shocked regions of the NGC1333-IRAS4A protostellar outflow

Water in Low-Mass Star-Forming Regions with Herschel: The Link Between Water Gas and Ice in Protostellar Envelopes

ALMA detection of the rotating molecular disk wind from the young star HD 163296

ALMA imaging of the CO snowline of the HD 163296 disk with DCO+

Deuterated water in the solar-type protostars NGC 1333 IRAS 4A and IRAS 4B

DIGIT survey of far-infrared lines from protoplanetary disks I

DIGIT: Herschel and Spitzer spectro-imaging of SMM3 and SMM4 in Serpens

Far-infrared molecular lines from Low- to High-Mass Star Forming Regions observed with Herschel

Herschel PACS observations of shocked gas associated with the jets of L1448 and L1157

Herschel-HIFI observations of high-J CO and isotopologues in star-forming regions: from low- to high-mass

High-pressure, low-abundance water in bipolar outflows. Results from a Herschel-WISH survey

Infrared Spectroscopic Survey of the Quiescent Medium of Nearby Clouds: I. Ice Formation and Grain Growth in Lupus

Observational evidence for dissociative shocks in the inner 100 AU of low-mass protostars using Herschel-HIFI

Probing the radial temperature structure of protoplanetary disks with Herschel/HIFI

The 69 micron forsterite band in spectra of protoplanetary disks - Results from the Herschel DIGIT programme

Water in star forming regions with Herschel (WISH) III. Far-infrared cooling lines in low-mass young stellar objects

Waterfalls around protostars: Infall motions towards Class 0/I envelopes as probed by water