Source author record

Michiel R. Hogerheijde

Michiel R. Hogerheijde appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

The TW Hya Rosetta Stone Project I: Radial and vertical distributions of DCN and DCO+

Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30~au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. There appears to be an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature, or something peculiar to the old TW Hya disk.

preprint2020arXiv

The TW Hya Rosetta Stone Project II: Spatially resolved emission of formaldehyde hints at low-temperature gas-phase formation

Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pathway dominates, gaseous H$_2$CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and $<3$ linked to cold formation in the solid-state. We present spatially resolved ALMA observations of multiple ortho- and para-H$_2$CO transitions in the TW Hya protoplanetary disk to test H$_2$CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of $33\pm2$ K, ($1.1\pm0.1)\times10^{12}$ cm$^{-2}$ and $25\pm2$ K, $(4.4\pm0.3)\times10^{11}$ cm$^{-2}$ for ortho- and para-H$_2$CO, respectively, and an OPR of $2.49\pm0.23$. A radially resolved analysis shows that the observed H$_2$CO emits mostly at rotational temperatures of 30-40 K, corresponding to a layer with $z/R\ge0.25$. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to $2.0\pm0.5$. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, lead us to speculate that gas-phase formation is responsible for the observed H$_2$CO across the TW Hya disk.

preprint2020arXiv

The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future.

preprint2016arXiv

First detection of gas-phase ammonia in a planet-forming disk

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

preprint2016arXiv

Misaligned Disks in the Binary Protostar IRS 43

Recent high angular resolution ($\sim$0.2") ALMA observations of the 1.1 mm continuum and of HCO+ J=3-2 and HCN J=3-2 gas towards the binary protostar IRS 43 reveal multiple Keplerian disks which are significantly misaligned ($\gt$ 60$^\circ$), both in inclination and position angle and also with respect to the binary orbital plane. Each stellar component has an associated circumstellar disk while the binary is surrounded by a circumbinary disk. Together with archival VLA measurements of the stellar positions over 25 years, and assuming a circular orbit, we use our continuum measurements to determine the binary separation, a = 74 $\pm$ 4 AU, and its inclination, i $\lt$ 30$^\circ$. The misalignment in this system suggests that turbulence has likely played a major role in the formation of IRS 43.

preprint2015arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: IV. Mechanical and radiative feedback

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

preprint2015arXiv

Steepening of the 820 micron continuum surface-brightness profile signals dust evolution in TW Hya's disk

Grain growth in planet-forming disks is the first step toward the formation of planets. The growth of grains and their inward drift leaves a distinct imprint on the dust surface-density distribution and the resulting surface-brightness profile of the thermal continuum emission. We determine the surface-brightness profile of the continuum emission using resolved observations at millimeter wavelengths of the disk around TW Hya, and infer the signature of dust evolution on the surface density and dust opacity. Archival ALMA observations at 820 micron on baselines up to 410 kilolambda are compared to parametrized disk models to determine the surface-brightness profile. Under the assumption of a constant dust opacity, a broken radial power law best describes the dust surface density, with a slope of -0.53 +/- 0.01 from the 4.1 au radius of the (already known) inner hole to a turn-over radius of 47.1 +/- 0.2 au, steepening to -8.0 +/- 0.1 at larger radii. The emission drops below the detection limit beyond ~60 au. The shape of the dust surface density is consistent with theoretical expectations for grain growth, fragmentation, and drift, but its total dust content and its turn-over radius are too large for TW Hya's age of 8-10 Myr even when taking into account a radially varying dust opacity. Higher resolution imaging with ALMA of TW Hya and other disks is required to establish if unseen gaps associated with, e.g., embedded planets trap grains at large radii or if locally enhanced grain growth associated with the CO snow line explains the extent of the millimeter-continuum surface brightness profile. In the latter case, population studies should reveal a correlation between the location of the CO snow line and the extent of the millimeter continuum. In the former case, and if CO freeze out promotes planet formation, this correlation should extend to the location of gaps as well.

preprint2015arXiv

Volatile depletion in the TW Hydrae disk atmosphere

An abundance decrease in carbon- and oxygen-bearing species relative to dust has been frequently found in planet-forming disks, which can be attributed to an overall reduction of gas mass. However, in the case of TW Hya, the only disk with gas mass measured directly with HD rotational lines, the inferred gas mass ($\lesssim$0.005 solar mass) is significantly below the directly measured value ($\gtrsim$0.05 solar mass). We show that this apparent conflict can be resolved if the elemental abundances of carbon and oxygen are reduced in the upper layers of the outer disk but are normal elsewhere (except for a possible enhancement of their abundances in the inner disk). The implication is that in the outer disk, the main reservoir of the volatiles (CO, water, ...) resides close to the midplane, locked up inside solid bodies that are too heavy to be transported back to the atmosphere by turbulence. An enhancement in the carbon and oxygen abundances in the inner disk can be caused by inward migration of these solid bodies. This is consistent with estimates based on previous models of dust grain dynamics. Indirect measurements of the disk gas mass and disk structure from species such as CO will thus be intertwined with the evolution of dust grains, and possibly also with the formation of planetesimals.

preprint2014arXiv

ALMA hints at the presence of two companions in the disk around HD 100546

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

preprint2014arXiv

Dynamical structure of the inner 100 AU of the deeply embedded protostar IRAS 16293-2422

A fundamental question about the early evolution of low-mass protostars is when circumstellar disks may form. High angular resolution observations of molecular transitions in the (sub)millimeter wavelength windows make it possible to investigate the kinematics of the gas around newly-formed stars, for example to identify the presence of rotation and infall. IRAS 16293-2422 was observed with the extended Submillimeter Array (eSMA) resulting in subarcsecond resolution (0.46" x 0.29", i.e. $\sim$ 55 $\times$ 35~AU) images of compact emission from the C$^{17}$O (3-2) and C$^{34}$S (7-6) transitions at 337~GHz (0.89~mm). To recover the more extended emission we have combined the eSMA data with SMA observations of the same molecules. The emission of C$^{17}$O (3-2) and C$^{34}$S (7-6) both show a velocity gradient oriented along a northeast-southwest direction with respect to the continuum marking the location of one of the components of the binary, IRAS16293A. Our combined eSMA and SMA observations show that the velocity field on the 50--400~AU scales is consistent with a rotating structure. It cannot be explained by simple Keplerian rotation around a single point mass but rather needs to take into account the enclosed envelope mass at the radii where the observed lines are excited. We suggest that IRAS 16293-2422 could be among the best candidates to observe a pseudo-disk with future high angular resolution observations.

preprint2014arXiv

Peculiar Near-Nucleus Outgassing of Comet 17P/Holmes During Its 2007 Outburst

We present high angular resolution Submillimeter Array observations ofthe outbursting Jupiter family comet 17P/Holmes on 2007 October 26-29, achieving a spatial resolution of 2.5", or ~3000 km at the comet distance. The observations resulted in detections of the rotational lines CO 3-2, HCN 4-3, H$^{13}$CN 4-3, CS 7-6, H$_2$CO 3$_{1,2}$-2$_{1,1}$, H$_2$S 2$_{2,0}$-2$_{1,1}$, and multiple CH$_3$OH lines, along with the associated dust continuum at 221 and 349 GHz. The continuum has a spectral index of 2.7$\pm$0.3, slightly steeper than blackbody emission from large dust particles. From the imaging data, we identify two components in the molecular emission. One component is characterized by a relatively broad line width (~1 km s$^{-1}$ FWHM) exhibiting a symmetric outgassing pattern with respect to the nucleus position. The second component has a narrower line width (<0.5 km s$^{-1}$ FWHM) with the line center red-shifted by 0.1-0.2 km s$^{-1}$ (cometocentric frame), and shows a velocity shift across the nucleus position with the position angle gradually changing from 66 to 30 degrees within the four days of observations. We determine distinctly different CO/HCN ratios for each of the components. For the broad-line component we find CO/HCN <7, while in the narrow-line component, CO/HCN = 40$\pm$5. We hypothesize that the narrow-line component originates from the ice grain halo found in near-nucleus photometry, believed to be created by sublimating recently released ice grains around the nucleus during the outburst. In this interpretation, the high CO/HCN ratio of this component reflects the more pristine volatile composition of nucleus material released in the outburst.

preprint2013arXiv

High-J CO survey of low-mass protostars observed with Herschel-HIFI

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

preprint2013arXiv

Imaging of the CO Snow Line in a Solar Nebula Analog

Planets form in the disks around young stars. Their formation efficiency and composition are intimately linked to the protoplanetary disk locations of "snow lines" of abundant volatiles. We present chemical imaging of the CO snow line in the disk around TW Hya, an analog of the solar nebula, using high spatial and spectral resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of N2H+, a reactive ion present in large abundance only where CO is frozen out. The N2H+ emission is distributed in a large ring, with an inner radius that matches CO snow line model predictions. The extracted CO snow line radius of ~ 30 AU helps to assess models of the formation dynamics of the Solar System, when combined with measurements of the bulk composition of planets and comets.

preprint2013arXiv

Rotationally-supported disks around Class I sources in Taurus: disk formation constraints

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

preprint2012arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: III. NGC 1333 IRAS 4A/4B envelope, outflow and UV heating

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

preprint2012arXiv

Evidence for Multiple Pathways to Deuterium Enhancements in Protoplanetary Disks

The distributions of deuterated molecules in protoplanetary disks are expected to depend on the molecular formation pathways. We use observations of spatially resolved DCN emission from the disk around TW Hya, acquired during ALMA Science verification with a ~3" synthesized beam, together with comparable DCO+ observations from the Submillimeter Array, to investigate differences in the radial distributions of these species and hence differences in their formation chemistry. In contrast to DCO+, which shows an increasing column density with radius, DCN is better fit by a model that is centrally peaked. We infer that DCN forms at a smaller radii and thus at higher temperatures than DCO+. This is consistent with chemical network model predictions of DCO+ formation from H2D+ at T<30 K and DCN formation from additional pathways involving CH2D+ at higher temperatures. We estimate a DCN/HCN abundance ratio of ~0.017, similar to the DCO+/HCO+ abundance ratio. Deuterium fractionation appears to be efficient at a range of temperatures in this protoplanetary disk. These results suggest caution in interpreting the range of deuterium fractions observed in Solar System bodies, as multiple formation pathways should be taken into account.

preprint2011arXiv

Detection of the Water Reservoir in a Forming Planetary System

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

preprint2011arXiv

Energetic processes revealed by spectrally resolved high-J CO lines in low-mass star-forming regions with Herschel-HIFI

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars in NGC1333. The observations show several energetic components including shocked and quiescent gas. Radiative transfer models are used to quantify the C18O envelope abundance which require a jump in the abundance at an evaporation temperature, T_ev ~25 K, providing new direct evidence of a CO ice evaporation zone around protostars. The abundance in the outermost part of the envelope, X_0, is within the canonical value of 2x10^-4; however the inner abundance, X_in, is found around a factor of 3-5 lower than X_0.

preprint2011arXiv

The Young Binary DQ Tau: A Hunt For X-ray Emission From Colliding Magnetospheres

The young high-eccentricity binary DQ Tau exhibits powerful recurring millimeter-band (mm) flaring attributed to collisions between the two stellar magnetospheres near periastron, when the stars are separated by only ~8Rstar. These magnetospheric interactions are expected to have scales and magnetic field strengths comparable to those of large X-ray flares from single pre-main-sequence (PMS) stars observed in the Chandra Orion Ultradeep Project (COUP). To search for X-rays arising from processes associated with colliding magnetospheres, we performed simultaneous X-ray and mm observations of DQ Tau near periastron phase. We report here several results. 1) As anticipated, DQ Tau was caught in a flare state in both mm and X-rays. A single long X-ray flare spanned the entire 16.5 hour Chandra exposure. 2) The inferred morphology, duration, and plasma temperature of the X-ray flare are typical of those of large flares from COUP stars. 3) However, our study provides three lines of evidence that this X-ray flare likely arises from colliding magnetospheres: the chance of capturing a large COUP-like flare within the span of our observation is small; the relative timing of the X-ray and mm flares indicates the Neupert effect and is consistent with a common coronal structure; the size of the emitting coronal structure (4-5Rstar) inferred from our analysis (which is admittedly model-dependent and should be considered with caution) is comparable to half the binary separation. 4) The peak flare X-ray luminosity is in agreement with an estimate of the power dissipated by magnetic reconnection within the framework of a simple model of interacting magnetospheres.

preprint2010arXiv

Empirical Constraints on Turbulence in Protoplanetary Accretion Disks

We present arcsecond-scale Submillimeter Array observations of the CO(3-2) line emission from the disks around the young stars HD 163296 and TW Hya at a spectral resolution of 44 m/s. These observations probe below the ~100 m/s turbulent linewidth inferred from lower-resolution observations, and allow us to place constraints on the turbulent linewidth in the disk atmospheres. We reproduce the observed CO(3-2) emission using two physical models of disk structure: (1) a power-law temperature distribution with a tapered density distribution following a simple functional form for an evolving accretion disk, and (2) the radiative transfer models developed by D'Alessio et al. that can reproduce the dust emission probed by the spectral energy distribution. Both types of models yield a low upper limit on the turbulent linewidth (Doppler b-parameter) in the TW Hya system (<40 m/s), and a tentative (3-sigma) detection of a ~300 m/s turbulent linewidth in the upper layers of the HD 163296 disk. These correspond to roughly <10% and 40% of the sound speed at size scales commensurate with the resolution of the data. The derived linewidths imply a turbulent viscosity coefficient, alpha, of order 0.01 and provide observational support for theoretical predictions of subsonic turbulence in protoplanetary accretion disks.

preprint2009arXiv

A break in the gas and dust surface density of the disc around the T Tauri star IM Lup

We study the distribution and physical properties of molecular gas in the disc around the T Tauri star IM Lup on scales close to 200 AU. We investigate how well the gas and dust distributions compare and work towards a unified disc model that can explain both gas and dust emission. 12CO, 13CO, and C18O J=2-1 line emission, as well as the dust continuum at 1.3 mm, is observed at 1.8" resolution towards IM Lup using the Submillimeter Array. A detailed disc model based on the dust emission is tested against these observations with the aid of a molecular excitation and radiative transfer code. Apparent discrepancies between the gas and dust distribution are investigated by adopting simple modifications to the existing model. The disc is seen at an inclination of 54+/-3 degrees and is in Keplerian rotation around a 0.8-1.6 Msun star. The outer disc radius traced by molecular gas emission is 900 AU, while the dust continuum emission and scattered light images limit the amount of dust present beyond 400 AU and are consistent with the existing model that assumes a 400 AU radius. Our observations require a drastic density decrease close to 400 AU with the vertical gas column density at 900 AU in the range of 5.d20 - 1.d22 cm-2. We derive a gas-to-dust mass ratio of 100 or higher in disc regions beyond 400 AU. Within 400 AU from the star our observations are consistent with a gas-to-dust ratio of 100 but other values are not ruled out.

preprint2009arXiv

Hot and cool water in Herbig Ae protoplanetary disks. A challenge for Herschel

The spatial origin and detectability of rotational H2O emission lines from Herbig Ae type protoplanetary disks beyond 70 micron is discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: (1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond the snow-line r>20AU, and (3) a layer of irradiated hot water at high altitudes z/r=0.1...0.3, extending from about 1AU to 30AU, where the kinetic gas temperature ranges from 200K to 1500K. Although region 3 contains only little amounts of water vapour (~3x10^-5 M_Earth), we find this warm layer to be almost entirely responsible for the rotational water emission lines, execpt for the 3 lowest excitation lines. Thus, Herschel will probe first and foremost the conditions and radial extension of region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust T_gas>T_dust. The observations do not allow for a determination of the snow-line, because the snow-line truncates the radial extension of region 1, whereas the lines originate from region 3. Different line transfer approximations (LTE, escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the escape probability method are found to underestimate the Monte Carlo results by 2%...45%.

preprint2009arXiv

PROSAC: A Submillimeter Array survey of low-mass protostars. II. The mass evolution of envelopes, disks, and stars from the Class 0 through I stages

The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through the circumstellar disk onto the central star. A sample of 20 Class 0 and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular resolution with the Submillimeter Array. Using detailed dust radiative transfer models, we have developed a framework for disentangling the continuum emission from the envelopes and disks, and from that estimated their masses. For the Class I sources in the sample, HCO+ 3-2 line emission has furthermore been observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, constraining the masses of the central stars. Both Class 0 and I protostars are surrounded by disks with typical masses of about 0.05 M_sun. No evidence is found for a correlation between the disk mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from 1 M_sun in the Class 0 stage to <0.1 M_sun in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses in the Class 0 stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the embedded protostellar stages. Theoretical models tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution and rapidly process material accreted from the larger scale envelope onto the central star.

preprint2008arXiv

Detection of CI in absorption towards PKS 1830-211 with the eSMA

We report the first science observations and results obtained with the "extended" SMA (eSMA), which is composed of the SMA (Submillimeter Array), JCMT (James Clerk Maxwell Telescope) and CSO (Caltech Submillimeter Observatory). Redshifted absorptions at z=0.886 of CI (^3P_1 - ^3P_0) were observed with the eSMA with an angular resolution of 0.55"x0.22" at 1.1 mm toward the southwestern image of the remarkable lensed quasar PKS 1830-211, but not toward the northeastern component at a separation of ~1". Additionally, SMA observations of CO, 13CO and C18O (all J=4-3) were obtained toward this object: CO was also detected toward the SW component, but none of the isotopologues were. This is the first time [CI] is detected in this object, allowing the first direct determination of relative abundances of neutral atomic carbon to CO in the molecular clouds of a spiral galaxy at z>0.1. The [CI] and CO profiles can be decomposed into two and three velocity components respectively. We derive C/CO column density ratios ranging from <0.5 (representative of dense cores) to ~2.5 (close to translucent clouds values). This could indicate that we are seeing environments with different physical conditions or that we are witnessing chemical evolution of regions where C has not completely been converted into CO.