Source author record

Geoffrey A. Blake

Geoffrey A. Blake appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.GA astro-ph.SR astro-ph physics.chem-ph astro-ph.IM physics.optics astro-ph.HE cond-mat.stat-mech

Catalog footprint

What is connected

62works

9topics

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

Protoplanetary disk cavities with JWST-MIRI: a dichotomy in molecular emission

The evolution of planet-forming regions in protoplanetary disks is of fundamental importance to understanding planet formation. Disks with a central deficit in dust emission, a "cavity", have long attracted interest as potential evidence for advanced disk clearing by protoplanets and/or winds. Before JWST, infrared spectra showed that these disks typically lack the strong molecular emission observed in full disks. In this work, we combine a sample of 12 disks with millimeter cavities of a range of sizes ($\sim2$-70 au) and different levels of millimeter and infrared continuum deficits. We analyze their molecular spectra as observed with MIRI on JWST, homogeneously reduced with the new JDISCS pipeline. This analysis demonstrates a stark dichotomy in molecular emission where "molecule-rich" (MR) cavities follow global trends between water, CO, and OH luminosity and accretion luminosity as in full disks, while "molecule-poor" (MP) cavities are significantly sub-luminous in all molecules except sometimes OH. Disk cavities generally show sub-luminous organic emission, higher OH/H$_2$O ratios, and suggest a lower water column density. The sub-thermal excitation of CO and water vibrational lines suggests a decreased gas density in the emitting layer in all cavities, supporting model expectations for C$_2$H$_2$ photodissociation. We discover a bifurcation in infrared index (lower in MR cavities) suggesting that the molecular dichotomy is linked to residual $μ$m-size dust within millimeter disk cavities. Put together, these results suggest a feedback process between dust depletion, gas density decrease, and molecule dissociation. Disk cavities may have a common evolutionary sequence where MR switch into MP over time.

preprint2022arXiv

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz/Raman spectroscopic techniques provide invaluable insight into molecular structure and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

preprint2022arXiv

New Constraints on Protoplanetary Disk Gas Masses in Lupus

Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H$_2$ content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologues of CO. However, without additional constraints, the degeneracy between H$_2$ mass and the elemental composition of the gas leads to large uncertainties in such estimates. Here we explore the gas compositions of seven disks from the Lupus sample representing a range of CO-to-dust ratios. With Band 6 and 7 ALMA observations, we measure line emission for HCO$^+$, HCN, and N$_2$H$^+$. We find a tentative correlation among the line fluxes for these three molecular species across the sample, but no correlation with $^{13}$CO or sub-mm continuum fluxes. For the three disks where N$_2$H$^+$ is detected, we find that a combination of high disk gas masses and sub-interstellar C/H and O/H are needed to reproduce the observed values. We find increases of $\sim$10-100$\times$ previous mass estimates are required to match the observed line fluxes. This study highlights how multi-molecular studies are essential for constraining the physical and chemical properties of the gas in populations of protoplanetary disks and that CO isotopologues alone are not sufficient for determining the mass of many observed disks.

preprint2022arXiv

Retrieving the C and O Abundances of HR 7672~AB: a Solar-Type Primary Star with a Benchmark Brown Dwarf

A benchmark brown dwarf (BD) is a BD whose properties (e.g., mass and chemical composition) are precisely and independently measured. Benchmark BDs are valuable in testing theoretical evolutionary tracks, spectral synthesis, and atmospheric retrievals for sub-stellar objects. Here, we report results of atmospheric retrieval on a synthetic spectrum and a benchmark BD -- HR 7672~B -- with \petit. First, we test the retrieval framework on a synthetic PHOENIX BT-Settl spectrum with a solar composition. We show that the retrieved C and O abundances are consistent with solar values, but the retrieved C/O is overestimated by 0.13-0.18, which is $\sim$4 times higher than the formal error bar. Second, we perform retrieval on HR 7672~B using high spectral resolution data (R=35,000) from the Keck Planet Imager and Characterizer (KPIC) and near infrared photometry. We retrieve [C/H], [O/H], and C/O to be $-0.24\pm0.05$, $-0.19\pm0.04$, and $0.52\pm0.02$. These values are consistent with those of HR 7672~A within 1.5-$σ$. As such, HR 7672~B is among only a few benchmark BDs (along with Gl 570~D and HD 3651~B) that have been demonstrated to have consistent elemental abundances with their primary stars. Our work provides a practical procedure of testing and performing atmospheric retrieval, and sheds light on potential systematics of future retrievals using high- and low-resolution data.

preprint2021arXiv

Mapping physical parameters in Orion KL at high spatial resolution

The Orion Kleinmann-Low nebula (Orion KL) is notoriously complex and exhibits a range of physical and chemical components. We conducted high angular resolution (sub-arcsecond) observations of $^13$CH$_{3}$OH $ν=0$ ($\sim$0.3$^{\prime\prime}$ and $\sim$0.7$^{\prime\prime}$) and CH$_3$CN $ν_8=1$ ($\sim$0.2$^{\prime\prime}$ and $\sim$0.9$^{\prime\prime}$) line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) to investigate Orion KL's structure on small spatial scales (${\le}350$ au). Gas kinematics, excitation temperatures, and column densities were derived from the molecular emission via a pixel-by-pixel spectral line fitting of the image cubes, enabling us to examine the small-scale variation of these parameters. Sub-regions of the Hot Core have a higher excitation temperature in a 0.2$^{\prime\prime}$ beam than a 0.9$^{\prime\prime}$ beam, indicative of possible internal sources of heating. Furthermore, the velocity field includes a bipolar ${\sim}7{-}8$ km s$^{-1}$ feature with a southeast-northwest orientation against the surrounding ${\sim}4{-}5$ km s$^{-1}$ velocity field, which may be due to an outflow. We also find evidence of a possible source of internal heating toward the Northwest Clump, since the excitation temperature there is higher in a smaller beam versus a larger beam. Finally, the region southwest of the Hot Core (Hot Core-SW) presents itself as a particularly heterogeneous region bridging the Hot Core and Compact Ridge. Additional studies to identify the (hidden) sources of luminosity and heating within Orion KL are necessary to better understand the nebula and its chemistry.

preprint2020arXiv

An ALMA Survey of H$_2$CO in Protoplanetary Disks

H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and centrally-depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H$_2$CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H$_2$CO excitation temperatures and column densities for four disks with multiple H$_2$CO line detections. The temperatures are between 20-50K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H$_2$CO emission in disks is generally emerging from the warm molecular layer, with some contributions from the colder midplane. Applying the same H$_2$CO excitation temperatures to all disks in the survey, we find that H$_2$CO column densities span almost three orders of magnitude ($\sim 5 \times 10^{11} - 5 \times 10^{14} \mathrm{cm}^{-2}$). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H$_2$CO compared to T Tauri disks, possibly because of less CO freeze-out. More H$_2$CO observations towards Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H$_2$CO and other oxygen-bearing organics efficiently form during planet formation.

preprint2020arXiv

An evolutionary study of volatile chemistry in protoplanetary disks

The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.

preprint2020arXiv

An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks

The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

preprint2020arXiv

Carbon-grain sublimation: a new top-down component of protostellar chemistry

Earth's carbon deficit has been an outstanding problem in our understanding of the formation of our Solar System. A possible solution would be the sublimation of carbon grains at the so-called soot line (~300 K) early in the planet-formation process. Here, we argue that the most likely signatures of this process are an excess of hydrocarbons and nitriles inside the soot line, and a higher excitation temperature for these molecules compared to oxygen-bearing complex organics that desorb around the water snowline (~100 K). Such characteristics have been reported in the literature, for example, in Orion KL, although not uniformly, potentially due to differences in observational settings and analysis methods of different studies or related to the episodic nature of protostellar accretion. If this process is active, this would mean that there is a heretofore unknown component to the carbon chemistry during the protostellar phase that is acting from the top down - starting from the destruction of larger species - instead of from the bottom up from atoms. In the presence of such a top-down component, the origin of organic molecules needs to be re-explored.

preprint2020arXiv

Contrast and Temperature Dependence of Multi-Epoch High-Resolution Cross-Correlation Exoplanet Spectroscopy

While high-resolution cross-correlation spectroscopy (HRCCS) techniques have proven effective at characterizing the atmospheres of transiting and non-transiting hot Jupiters, the limitations of these techniques are not well understood. We present a series of simulations of one HRCCS technique, which combines the cross-correlation functions from multiple epochs, to place temperature and contrast limits on the accessible exoplanet population for the first time. We find that planets approximately Saturn-size and larger within $\sim$0.2 AU of a Sun-like star are likely to be detectable with current instrumentation in the $L$-band, a significant expansion compared with the previously-studied population. Cooler ($ \rm T_{eq} \leq 1000$ K) exoplanets are more detectable than suggested by their photometric contrast alone as a result of chemical changes which increase spectroscopic contrast. The $L$-band CH$_4$ spectrum of cooler exoplanets enables robust constraints on the atmospheric C/O ratio at $\rm T_{eq} \sim 900K$, which have proven difficult to obtain for hot Jupiters. These results suggest that the multi-epoch approach to HRCCS can detect and characterize exoplanet atmospheres throughout the inner regions of Sun-like systems with existing high-resolution spectrographs. We find that many epochs of modest signal-to-noise ($\rm S/N_{epoch} \sim 1500$) yield the clearest detections and constraints on C/O, emphasizing the need for high-precision near-infrared telluric correction with short integration times.

preprint2020arXiv

Simulating the Multi-Epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-Transiting Hot Jupiter HD187123B

We report the 6.5$σ$ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity $K_p$ of 53 $\pm$ 13 km/s. This high confidence detection is made using a multi-epoch, high resolution, cross correlation technique, and corresponds to a planetary mass of 1.4$^{+0.5}_{-0.3}$ $M_J$ and an orbital inclination of 21 $\pm$ 5$^{\circ}$. The technique works by treating the planet/star system as a spectroscopic binary and obtaining high signal-to-noise, high resolution observations at multiple points across the planet's orbit to constrain the system's binary dynamical motion. All together, seven epochs of Keck/NIRSPEC $L$-band observations were obtained, with five before the instrument upgrade and two after. Using high resolution SCARLET planetary and PHOENIX stellar spectral models, along with a line-by-line telluric absorption model, we were able to drastically increase the confidence of the detection by running simulations that could reproduce, and thus remove, the non-random structured noise in the final likelihood space well. The ability to predict multi-epoch results will be extremely useful for furthering the technique. Here, we use these simulations to compare three different approaches to combining the cross correlations of high resolution spectra and find that the Zucker 2003 log(L) approach is least affected by unwanted planet/star correlation for our HD187123 data set. Furthermore, we find that the same total S/N spread across an orbit in many, lower S/N epochs rather than fewer, higher S/N epochs could provide a more efficient detection. This work provides a necessary validation of multi-epoch simulations which can be used to guide future observations and will be key to studying the atmospheres of further separated, non-transiting exoplanets.

preprint2020arXiv

Sum-Frequency Signals in 2D-Terahertz-Terahertz-Raman Spectroscopy

We demonstrate that halogenated methane 2D-Terahertz Terahertz Raman (2D-TTR) spectra are determined by the complicated structure of the instrument response function (IRF) along $f_1$ and by the molecular coherences along $f_2$. Experimental improvements have helped increase the resolution and dynamic range of the measurements, including accurate THz pulse shape characterization. Sum-frequency excitations convolved with the IRF are found to quantitatively reproduce the 2D-TTR signal. A new Reduced Density Matrix model which incorporates sum-frequency pathways, with linear and harmonic operators fully supports this (re)interpretation of the 2D-TTR spectra.

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

CSO and CARMA Observations of L1157. II. Chemical Complexity in the Shocked Outflow

L1157, a molecular dark cloud with an embedded Class 0 protostar possessing a bipolar outflow, is an excellent source for studying shock chemistry, including grain-surface chemistry prior to shocks, and post-shock, gas-phase processing. The L1157-B1 and B2 positions experienced shocks at an estimated ~2000 and 4000 years ago, respectively. Prior to these shock events, temperatures were too low for most complex organic molecules to undergo thermal desorption. Thus, the shocks should have liberated these molecules from the ice grain-surfaces en masse, evidenced by prior observations of SiO and multiple grain mantle species commonly associated with shocks. Grain species, such as OCS, CH3OH, and HNCO, all peak at different positions relative to species that are preferably formed in higher velocity shocks or repeatedly-shocked material, such as SiO and HCN. Here, we present high spatial resolution (~3") maps of CH3OH, HNCO, HCN, and HCO+ in the southern portion of the outflow containing B1 and B2, as observed with CARMA. The HNCO maps are the first interferometric observations of this species in L1157. The maps show distinct differences in the chemistry within the various shocked regions in L1157B. This is further supported through constraints of the molecular abundances using the non-LTE code RADEX (Van der Tak et al. 2007). We find the east/west chemical differentiation in C2 may be explained by the contrast of the shock's interaction with either cold, pristine material or warm, previously-shocked gas, as seen in enhanced HCN abundances. In addition, the enhancement of the HNCO abundance toward the the older shock, B2, suggests the importance of high-temperature O-chemistry in shocked regions.

preprint2016arXiv

Discovery of the Interstellar Chiral Molecule Propylene Oxide (CH$_3$CHCH$_2$O)

Life on Earth relies on chiral molecules, that is, species not superimposable on their mirror images. This manifests itself in the selection of a single molecular handedness, or homochirality, across the biosphere. We present the astronomical detection of a chiral molecule, propylene oxide (CH$_3$CHCH$_2$O), in absorption toward the Galactic Center. Propylene oxide is detected in the gas phase in a cold, extended molecular shell around the embedded, massive protostellar clusters in the Sagittarius B2 star-forming region. This material is representative of the earliest stage of solar system evolution in which a chiral molecule has been found.

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

Hydrocarbon emission rings in protoplanetary disks induced by dust evolution

We report observations of resolved C2H emission rings within the gas-rich protoplanetary disks of TW Hya and DM Tau using the Atacama Large Millimeter Array (ALMA). In each case the emission ring is found to arise at the edge of the observable disk of mm-sized grains (pebbles) traced by (sub)mm-wave continuum emission. In addition, we detect a C3H2 emission ring with an identical spatial distribution to C2H in the TW Hya disk. This suggests that these are hydrocarbon rings (i.e. not limited to C2H). Using a detailed thermo-chemical model we show that reproducing the emission from C2H requires a strong UV field and C/O > 1 in the upper disk atmosphere and outer disk, beyond the edge of the pebble disk. This naturally arises in a disk where the ice-coated dust mass is spatially stratified due to the combined effects of coagulation, gravitational settling and drift. This stratification causes the disk surface and outer disk to have a greater permeability to UV photons. Furthermore the concentration of ices that transport key volatile carriers of oxygen and carbon in the midplane, along with photochemical erosion of CO, leads to an elemental C/O ratio that exceeds unity in the UV-dominated disk. Thus the motions of the grains, and not the gas, lead to a rich hydrocarbon chemistry in disk surface layers and in the outer disk midplane.

preprint2016arXiv

On the Commonality of 10-30AU Sized Axisymmetric Dust Structures in Protoplanetary Disks

An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive mm-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small scale structures. Here we investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau), a sample of disks with the highest spatial resolution ALMA observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ~axi-symmetric but rich in 10-30AU-sized radial structures, possibly due to physical gaps, surface density enhancements or localized dust opacity variations within the disks. These results suggest that small scale axi-symmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snowlines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snowlines.

preprint2016arXiv

The Radial Distribution of H2 and CO in TW Hya as Revealed by Resolved ALMA Observations of CO Isotopologues

CO is widely used as a tracer of molecular gas. However, there is now mounting evidence that gas phase carbon is depleted in the disk around TW Hya. Previous efforts to quantify this depletion have been hampered by uncertainties regarding the radial thermal structure in the disk. Here we present resolved ALMA observations of 13CO 3-2, C18O 3-2, 13CO 6-5, and C18O 6-5 emission in TW Hya, which allow us to derive radial gas temperature and gas surface density profiles, as well as map the CO abundance as a function of radius. These observations provide a measurement of the surface CO snowline at ~30 AU and show evidence for an outer ring of CO emission centered at 53 AU, a feature previously seen only in less abundant species. Further, the derived CO gas temperature profile constrains the freeze-out temperature of CO in the warm molecular layer to < 21 K. Combined with the previous detection of HD 1-0, these data constrain the surface density of the warm H2 gas in the inner ~30 AU. We find that CO is depleted by two orders of magnitude from R=10-60 AU, with the small amount of CO returning to the gas phase inside the surface CO snowline insufficient to explain the overall depletion. Finally, this new data is used in conjunction with previous modeling of the TW Hya disk to constrain the midplane CO snowline to 17-23 AU.

preprint2016arXiv

THz Time-Domain Spectroscopy of Mixed CO2-CH3OH Interstellar Ice Analogs

The icy mantles of interstellar dust grains are the birthplaces of the primordial prebiotic molecular inventory that may eventually seed nascent solar systems and the planets and planetesimals that form therein. Here, we present a study of two of the most abundant species in these ices after water: carbon dioxide (CO2) and methanol (CH3OH) using TeraHertz (THz) time-domain spectroscopy and mid-infrared spectroscopy. We study pure and mixed-ices of these species, and demonstrate the power of the THz region of the spectrum to elucidate the long-range structure (i.e. crystalline versus amorphous) of the ice, the degree of segregation of these species within the ice, and the thermal history of the species within the ice. Finally, we comment on the utility of the THz transitions arising from these ices for use in astronomical observations of interstellar ices.

preprint2015arXiv

CSO and CARMA Observations of L1157. I. A Deep Search for Hydroxylamine (NH$_2$OH)

A deep search for the potential glycine precursor hydroxylamine (NH$_2$OH) using the Caltech Submillimeter Observatory (CSO) at $λ= 1.3$ mm and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at $λ= 3$ mm is presented toward the molecular outflow L1157, targeting the B1 and B2 shocked regions. We report non-detections of NH$_2$OH in both sources. We a perform non-LTE analysis of CH$_3$OH observed in our CSO spectra to derive kinetic temperatures and densities in the shocked regions. Using these parameters, we derive upper limit column densities of NH$_2$OH of $\leq1.4 \times 10^{13}$~cm$^{-2}$ and $\leq1.5 \times 10^{13}$~cm$^{-2}$ toward the B1 and B2 shocks, respectively, and upper limit relative abundances of $N_{NH_2OH}/N_{H_2} \leq1.4 \times 10^{-8}$ and $\leq1.5 \times 10^{-8}$, respectively.

preprint2015arXiv

Detection of water vapor in the terrestrial planet forming region of a transition disk

We report a detection of water vapor in the protoplanetary disk around DoAr 44 with the Texas Echelon Cross Echelle Spectrograph --- a visitor instrument on the Gemini north telescope. The DoAr 44 disk consists of an optically thick inner ring and outer disk, separated by a dust-cleared 36 AU gap, and has therefore been termed "pre-transitional". To date, this is the only disk with a large inner gap known to harbor detectable quantities of warm (T=450 K) water vapor. In this work, we detect and spectrally resolve three mid-infrared pure rotational emission lines of water vapor from this source, and use the shapes of the emission lines to constrain the location of the water vapor. We find that the emission originates near 0.3 AU --- the inner disk region. This characteristic region coincides with that inferred for both optically thick and thin thermal infrared dust emission, as well as rovibrational CO emission. The presence of water in the dust-depleted region implies substantial columns of hydrogen (>10^{22} cm-2) as the water vapor would otherwise be destroyed by photodissociation. Combined with the dust modeling, this column implies a gas/small-dust ratio in the optically thin dusty region of >1000. These results demonstrate that DoAr 44 has maintained similar physical and chemical conditions to classical protoplanetary disks in its terrestrial-planet forming regions, in spite of having formed a large gap.

preprint2015arXiv

Dimming and CO absorption toward the AA Tau protoplanetary disk: An infalling flow caused by disk instability?

AA Tau, a classical T Tauri star in the Taurus cloud, has been the subject of intensive photometric monitoring for more than two decades due to its quasi-cyclic variation in optical brightness. Beginning in 2011, AA Tau showed another peculiar variation -- its median optical though near-IR flux dimmed significantly, a drop consistent with a 4-mag increase in visual extinction. It has stayed in the faint state since.Here we present 4.7um CO rovibrational spectra of AA Tau over eight epochs, covering an eleven-year time span, that reveal enhanced 12CO and 13CO absorption features in the $J_{\rm low}\leqslant$13 transitions after the dimming. These newly appeared absorptions require molecular gas along the line of sight with T~500 K and a column density of log (N12CO)~18.5 cm^{-2}, with line centers that show a constant 6 km s$^{-1}$ redshift. The properties of the molecular gas confirm an origin in the circumstellar material. We suggest that the dimming and absorption are caused by gas and dust lifted to large heights by a magnetic buoyancy instability. This material is now propagating inward, and on reaching the star within a few years will be observed as an accretion outburst.

preprint2015arXiv

Evidence of fast pebble growth near condensation fronts in the HL Tau protoplanetary disk

Water and simple organic molecular ices dominate the mass of solid materials available for planetesimal and planet formation beyond the water snow line. Here we analyze ALMA long baseline 2.9, 1.3 and 0.87 mm continuum images of the young star HL Tau, and suggest that the emission dips observed are due to rapid pebble growth around the condensation fronts of abundant volatile species. Specifically, we show that the prominent innermost dip at 13 AU is spatially resolved in the 0.87 mm image, and its center radius is coincident with the expected mid-plane condensation front of water ice. In addition, two other prominent dips, at distances of 32 and 63 AU, cover the mid-plane condensation fronts of pure ammonia or ammonia hydrates and clathrate hydrates (especially with CO and N$_2$) formed from amorphous water ice. The spectral index map of HL Tau between 1.3 and 0.87 mm shows that the flux ratios inside the dips are statistically larger than those of nearby regions in the disk. This variation can be explained by a model with two dust populations, where most of solid mass resides in a component that has grown into decimeter size scales inside the dips. Such growth is in accord with recent numerical simulations of volatile condensation, dust coagulation and settling.

preprint2015arXiv

Herschel observations of EXtra-Ordinary Sources: Analysis of the HIFI 1.2 THz Wide Spectral Survey Toward Orion KL II. Chemical Implications

We present chemical implications arising from spectral models fit to the Herschel/HIFI spectral survey toward the Orion Kleinmann-Low nebula (Orion KL). We focus our discussion on the eight complex organics detected within the HIFI survey utilizing a novel technique to identify those molecules emitting in the hottest gas. In particular, we find the complex nitrogen bearing species CH$_{3}$CN, C$_{2}$H$_{3}$CN, C$_{2}$H$_{5}$CN, and NH$_{2}$CHO systematically trace hotter gas than the oxygen bearing organics CH$_{3}$OH, C$_{2}$H$_{5}$OH, CH$_{3}$OCH$_{3}$, and CH$_{3}$OCHO, which do not contain nitrogen. If these complex species form predominantly on grain surfaces, this may indicate N-bearing organics are more difficult to remove from grain surfaces than O-bearing species. Another possibility is that hot (T$_{\rm kin}$$\sim$300 K) gas phase chemistry naturally produces higher complex cyanide abundances while suppressing the formation of O-bearing complex organics. We compare our derived rotation temperatures and molecular abundances to chemical models, which include gas-phase and grain surface pathways. Abundances for a majority of the detected complex organics can be reproduced over timescales $\gtrsim$ 10$^{5}$ years, with several species being under predicted by less than 3$σ$. Derived rotation temperatures for most organics, furthermore, agree reasonably well with the predicted temperatures at peak abundance. We also find that sulfur bearing molecules which also contain oxygen (i.e. SO, SO$_{2}$, and OCS) tend to probe the hottest gas toward Orion KL indicating the formation pathways for these species are most efficient at high temperatures.

preprint2015arXiv

Measurements of water surface snow lines in classical protoplanetary disks

We present deep Herschel-PACS spectroscopy of far-infrared water lines from a sample of four protoplanetary disks around solar-mass stars, selected to have strong water emission at mid-infrared wavelengths. By combining the new Herschel spectra with archival Spitzer-IRS spectroscopy, we retrieve a parameterized radial surface water vapor distribution from 0.1-100 AU using two-dimensional dust and line radiative transfer modeling. The surface water distribution is modeled with a step model comprising of a constant inner and outer relative water abundance and a critical radius at which the surface water abundance is allowed to change. We find that the four disks have critical radii of $\sim 3-11$ AU, at which the surface water abundance decreases by at least 5 orders of magnitude. The measured values for the critical radius are consistently smaller than the location of the surface snow line, as predicted by the observed spectral energy distribution. This suggests that the sharp drop-off of the surface water abundance is not solely due to the local gas-solid balance, but may also be driven by the de-activation of gas-phase chemical pathways to water below 300 K. Assuming a canonical gas-to-dust ratio of 100, as well as coupled gas and dust temperatures $T_{\rm gas}=T_{\rm dust}$, the best-fit inner water abundances become implausibly high (0.01-1.0 ${\rm H_{2}}^{-1}$). Conversely, a model in which the gas and dust temperatures are decoupled leads to canonical inner disk water abundances of $\sim 10^{-4} \rm H_{2}^{-1}$, while retaining gas-to-dust ratios of 100. That is, the evidence for gas-dust decoupling in disk surfaces is stronger than for enhanced gas-to-dust ratios.

preprint2015arXiv

Nonlinear TeraHertz Coherent Excitation of Vibrational Modes of Liquids

We report the first coherent excitation of intramolecular vibrational modes via the nonlinear interaction of a TeraHertz (THz) light field with molecular liquids. A terahertz-terahertz-Raman pulse sequence prepares the coherences with a broadband, high-energy, (sub)picosecond terahertz pulse, that are then measured in a terahertz Kerr effect spectrometer via phase-sensitive, heterodyne detection with an optical pulse. The spectrometer reported here has broader terahertz frequency coverage, and an increased sensitivity relative to previously reported terahertz Kerr effect experiments. Vibrational coherences are observed in liquid diiodomethane at 3.66 THz (122 cm$^{-1}$), and in carbon tetrachloride at 6.50 THz (217 cm$^{-1}$), in exact agreement with literature values of those intramolecular modes. This work opens the door to 2D spectroscopies, nonlinear in terahertz field, that can study the dynamics of condensed-phase molecular systems, as well as coherent control at terahertz frequencies.

preprint2015arXiv

Tracing the Ingredients for a Habitable Earth from Interstellar Space through Planet Formation

We use the C/N ratio as a monitor of the delivery of key ingredients of life to nascent terrestrial worlds. Total elemental C and N contents, and their ratio, are examined for the interstellar medium, comets, chondritic meteorites and terrestrial planets; we include an updated estimate for the Bulk Silicate Earth (C/N = 49.0 +/- 9.3). Using a kinetic model of disk chemistry, and the sublimation/condensation temperatures of primitive molecules, we suggest that organic ices and macro-molecular (refractory or carbonaceous dust) organic material are the likely initial C and N carriers. Chemical reactions in the disk can produce nebular C/N ratios of ~1-12, comparable to those of comets and the low end estimated for planetesimals. An increase of the C/N ratio is traced between volatile-rich pristine bodies and larger volatile-depleted objects subjected to thermal/accretional metamorphism. The C/N ratios of the dominant materials accreted to terrestrial planets should therefore be higher than those seen in carbonaceous chondrites or comets. During planetary formation, we explore scenarios leading to further volatile loss and associated C/N variations owing to core formation and atmospheric escape. Key processes include relative enrichment of nitrogen in the atmosphere and preferential sequestration of carbon by the core. The high C/N BSE ratio therefore is best satisfied by accretion of thermally processed objects followed by large-scale atmospheric loss. These two effects must be more profound if volatile sequestration in the core is effective. The stochastic nature of these processes hints that the surface/atmospheric abundances of biosphere-essential materials will likely be variable.

preprint2014arXiv

A CSO Search for $l$-C$_3$H$^+$: Detection in the Orion Bar PDR

The results of a Caltech Submillimeter Observatory (CSO) search for $l$-C$_3$H$^+$, first detected by Pety et al. (2012) in observations toward the Horsehead photodissociation region (PDR), are presented. A total of 39 sources were observed in the 1 mm window. Evidence of emission from $l$-C$_3$H$^+$ is found in only a single source - the Orion Bar PDR region, which shows a rotational temperature of 178(13) K and a column density of 7(2) x $10^{11}$ cm$^{-2}$. In the remaining sources, upper limits of ~10$^{11} - 10^{13}$ cm$^{-2}$ are found. These results are discussed in the context of guiding future observational searches for this species.

preprint2014arXiv

ALMA Observations of the T Tauri Binary System AS 205: Evidence for Molecular Winds and/or Binary Interactions

In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind (Pontoppidan et al. 2011; Bast et al. 2011), also displays significant deviations from Keplerian disk emission in the observations presented here. Detections near both AS 205 N and S are obtained in 1.3 mm continuum, 12CO 2-1, 13CO 2-1 and C18O 2-1. The 12CO emission is extended up to 2 arcsec from AS 205N, and both 12CO and 13CO display deviations from Keplerian rotation at all angular scales. Two possible explanations for these observations hold up best to close scrutiny - tidal interaction with AS 205 S or disk winds (or a combination of the two), and we discuss these possibilities in some detail.

preprint2014arXiv

Comparison of the dust and gas radial structure in the transition disk [PZ99] J160421.7-213028

We present ALMA observations of the 880um continuum and CO J= 3-2 line emission from the transition disk around [PZ99] J160421.7-213028, a solar mass star in the Upper Scorpius OB association. Analysis of the continuum data indicates that 80% of the dust mass is concentrated in an annulus extending between 79 and 114AU in radius. Dust is robustly detected inside the annulus, at a mass surface density 100 times lower than that at 80 AU. The CO emission in the inner disk also shows a significantly decreased mass surface density, but we infer a cavity radius of only 31AU for the gas. The large separation of the dust and gas cavity edges, as well as the high radial concentration of millimeter-sized dust grains, is qualitatively consistent with the predictions of pressure trap models that include hydrodynamical disk-planet interactions and dust coagulation/fragmentation processes.

preprint2014arXiv

Near-IR Direct Detection of Water Vapor in Tau Boo b

We use high dynamic range, high-resolution L-band spectroscopy to measure the radial velocity variations of the hot Jupiter in the tau Bootis planetary system. The detection of an exoplanet by the shift in the stellar spectrum alone provides a measure of the planet's minimum mass, with the true mass degenerate with the unknown orbital inclination. Treating the tau Boo system as a high flux ratio double-lined spectroscopic binary permits the direct measurement of the planet's true mass as well as its atmospheric properties. After removing telluric absorption and cross-correlating with a model planetary spectrum dominated by water opacity, we measure a 6-sigma detection of the planet at K_p = 111 +- 5 km/s, with a 1-sigma upper limit on the spectroscopic flux ratio of 10^-4. This radial velocity leads to a planetary orbital inclination of i = 45+3-4degrees and a mass of M_P = 5.90+0.35-0.20 M_ Jup. We report the first detection of water vapor in the atmosphere of a non-transiting hot Jupiter, tau Boo b.

preprint2014arXiv

The Search for a Complex Molecule in a Selected Hot Core Region: A Rigorous Attempt to Confirm trans-Ethyl Methyl Ether toward W51 e1/e2

An extensive search has been conducted to confirm transitions of \textit{trans}-ethyl methyl ether (tEME, C$_2$H$_5$OCH$_3$), toward the high mass star forming region W51 e1/e2 using the 12 m Telescope of the Arizona Radio Observatory (ARO) at wavelengths from 2 mm and 3 mm. In short, we cannot confirm the detection of tEME toward W51 e1/e2 and our results call into question the initial identification of this species by \citet{FuchsSpace}. Additionally, reevaluation of the data from the original detection indicates that tEME is not present toward W51 e1/e2 in the abundance reported by Fuchs and colleagues. Typical peak-to-peak noise levels for the present observations of W51 e1/e2 were between 10 - 30 mK, yielding an upper limit of the tEME column density of $\leq$ 1.5 $\times$ 10$^{15}$ cm$^{-2}$. This would make tEME at least a factor 2 times less abundant than dimethyl ether (CH$_3$OCH$_3$) toward W51 e1/e2. We also performed an extensive search for this species toward the high mass star forming region Sgr B2(N-LMH) with the NRAO 100 m Green Bank Telescope (GBT). No transitions of tEME were detected and we were able to set an upper limit to the tEME column density of $\leq$ 4 $\times$ 10$^{14}$ cm$^{-2}$ toward this source. Thus, we are able to show that tEME is not a new molecular component of the interstellar medium and that an exacting assessment must be carried out when assigning transitions of new molecular species to astronomical spectra to support the identification of large organic interstellar molecules.

preprint2013arXiv

A Search for l-C3H+ and l-C3H in Sgr B2(N), Sgr B2(OH), and the Dark Cloud TMC-1

Pety et al. (2012) recently reported the detection of several transitions of an unknown carrier in the Horsehead PDR and attribute them to l-C3H+. Here, we have tested the predictive power of their fit by searching for, and identifying, the previously unobserved J=1-0 and J=2-1 transitions of the unknown carrier (B11244) towards Sgr B2(N) in data from the publicly available PRIMOS project. Also presented here are observations of the J=6-5 and J=7-6 transitions towards Sgr B2(N) and Sgr B2(OH) using the Barry E. Turner Legacy Survey and results from the Kaifu et al. (2004) survey of TMC-1. We calculate an excitation temperature and column density of B11244 of ~10 K and ~10^13 cm-2 in Sgr B2(N) and ~79 K with an upper limit of < 1.5 x 10^13 cm-2 in Sgr B2(OH) and find trace evidence for the cation's presence in TMC-1. Finally, we present spectra of the neutral species in both Sgr B2(N) and TMC-1, and comment on the robustness of the assignment of the detected signals to l-C3H+.

preprint2013arXiv

A Significantly Low CO Abundance Toward the TW Hya Protoplanetary Disk: A Path to Active Carbon Chemistry?

In this Letter we report the CO abundance relative to H2 derived toward the circumstellar disk of the T-Tauri star TW Hya from the HD (1-0) and C18O (2-1) emission lines. The HD (1-0) line was observed by the Herschel Space Observatory Photodetector Array Camera and Spectrometer whereas C18O (2-1) observations were carried out with the Submillimeter Array at a spatial resolution of 2.8" x 1.9" (corresponding to 142 x 97 AU). In the disk's warm molecular layer (T>20 K) we measure a disk-averaged gas-phase CO abundance relative to H2 of $χ{\rm(CO)}=(0.1-3)x10^{-5}$, substantially lower than the canonical value of $χ{\rm(CO)}=10^{-4}$. We infer that the best explanation of this low $χ$(CO) is the chemical destruction of CO followed by rapid formation of carbon chains, or perhaps CO2, that can subsequently freeze-out, resulting in the bulk mass of carbon locked up in ice grain mantles and oxygen in water. As a consequence of this likely time-dependent carbon sink mechanism, CO may be an unreliable tracer of H2 gas mass.

preprint2013arXiv

An observational investigation of the identity of B11244 ($l$-C$_3$H$^+$/C$_3$H$^-$)

Pety et al. (2012) have reported the detection of eight transitions of a closed-shell, linear molecule (B11244) in observations toward the Horsehead PDR, which they attribute to the $l$-C$_3$H$^+$ cation. Recent high-level \textit{ab initio} calculations have called this assignment into question; the anionic C$_3$H$^-$ molecule has been suggested as a more likely candidate. Here, we examine observations of the Horsehead PDR, Sgr B2(N), TMC-1, and IRC+10216 in the context of both $l$-C$_3$H$^+$ and C$_3$H$^-$. We find no observational evidence of $K_a = 1$ lines which should be present were the carrier indeed C$_3$H$^-$. Additionally, we find a strong anti-correlation between the presence of known molecular anions and B11244 in these regions. Finally, we discuss the formation and destruction chemistry of C$_3$H$^-$ in the context of the physical conditions in the regions. Based on these results, we conclude there is little evidence to support the carrier is C$_3$H$^-$.

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

Measuring Protoplanetary Disk Accretion with H I Pf-beta

In this work, we introduce the use of H I Pfund-beta (4.6538 micron) as a tracer of mass accretion from protoplanetary disks onto young stars. Pf-beta was serendipitously observed in NIRSPEC and CRIRES surveys of CO fundamental emission, amounting to a sample size of 120 young stars with detected Pf-beta emission. Using a subsample of disks with previously measured accretion luminosities, we show that Pf-beta line luminosity is well correlated with accretion luminosity over a range of at least three orders of magnitude. We use this correlation to derive accretion luminosities for all 120 targets, 65 of which are previously unreported in the literature. We also report accretion rates for 67 targets, 16 previously unmeasured. Our large sample size and our ability to probe high extinction values allow for relatively unbiased comparisons between different types of disks. We find that the transitional disks in our sample have lower than average Pf-beta line luminosities, and thus accretion luminosities, at a marginally significant level. We also show that high Pf-beta equivalent width is a signature of transitional disks with high inner disk gas/dust ratios. In contrast, we find that disks with signatures of slow disk winds have Pf-beta luminosities comparable to those of other disks in our sample. Finally, we investigate accretion rates for stage I disks, including significantly embedded targets. We find that stage I and stage II disks have statistically indistinguishable Pf-beta line luminosities, implying similar accretion rates, and that the accretion rates of stage I disks are too low to be consistent with quiescent accretion. Our results are instead consistent with both observational and theoretical evidence that stage I objects experience episodic, rather than quiescent, accretion.

preprint2013arXiv

VLT-CRIRES Survey of Rovibrational CO Emission from Protoplanetary Disks

We present a large, comprehensive survey of rovibrational CO line emission at 4.7 micron from 69 protoplanetary disks, obtained with CRIRES on the ESO Very Large Telescope at the highest available spectral resolving power (R=95,000, v=3.2 km/s). The CO fundamental band (Delta v=1) is a well-known tracer of warm gas in the inner, planet-forming regions of gas-rich disks around young stars, with the lines formed in the super-heated surfaces of the disks at radii of 0.1-10 AU. Our high spectral resolution data provide new insight into the kinematics of the inner disk gas. Pure double-peaked Keplerian profiles are surprisingly uncommon, beyond the frequency expected based on disk inclination. The majority of the profiles are consistent with emission from a disk plus a slow (few km/s) molecular disk wind. This is evidenced by analysis of different categories as well as an overall tendency for line profiles to have excess blue emission. Weak emission lines from isotopologues and vibrationally excited levels are readily detected. In general, 13CO lines trace cooler gas than the bulk 12CO emission and may arise from further out in the disk, as indicated by narrower line profiles. A high fraction of the sources show vibrationally excited emission (~50%) which is correlated with accretion luminosity, consistent with ultra-violet (UV) fluorescent excitation. Disks around early-type Herbig AeBe stars have narrower lines, on average, than their lower-mass late-type counterparts, due to their increased luminosity. Evolutionary changes in CO are also seen. Removal of the protostellar envelope between class I and II results in the disappearance of the strong absorption lines and CO ice feature characteristic of class I spectra. However, CO emission from class I and II objects is similar in detection frequency, excitation and line shape, indicating that inner disk characteristics are established early.

preprint2012arXiv

Interstellar Carbodiimide (HNCNH) - A New Astronomical Detection from the GBT PRIMOS Survey via Maser Emission Features

In this work, we identify carbodiimide (HNCNH), which is an isomer of the well-known interstellar species cyanamide (NH2CN), in weak maser emission, using data from the GBT PRIMOS survey toward Sgr B2(N). All spectral lines observed are in emission and have energy levels in excess of 170 K, indicating that the molecule likely resides in relatively hot gas that characterizes the denser regions of this star forming region. The anticipated abundance of this molecule from ice mantle experiments is ~10% of the abundance of NH2CN, which in Sgr B2(N) corresponds to ~2 x 10^13 cm-2. Such an abundance results in transition intensities well below the detection limit of any current astronomical facility and, as such, HNCNH could only be detected by those transitions which are amplified by masing.

preprint2011arXiv

CO rovibrational emission as a probe of inner disk structure

We present an analysis of CO emission lines from a sample of T Tauri, Herbig Ae/Be, and transitional disks with known inclinations, in order to study the structure of inner disk molecular gas. We calculate CO inner radii by fitting line profiles with a simple parameterized model. We find that, for optically thick disks, CO inner radii are strongly correlated with the total system luminosity (stellar plus accretion), and consistent with the dust sublimation radius. Transitional disk inner radii show the same trend with luminosity, but are systematically larger. Using rotation diagram fits, we derive, for classical T Tauri disks, emitting areas consistent with a ring of width ~0.15 AU located at the CO inner radius; emitting areas for transitional disks are systematically smaller. We also measure lower rotational temperatures for transitional disks, and disks around Herbig Ae/Be stars, than for those around T Tauri stars. Finally, we find that rotational temperatures are similar to, or slightly lower than, the expected temperature of blackbody grains located at the CO inner radius, in contrast to expectations of thermal decoupling between gas and dust.

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

Ices in starless and starforming cores

Icy grain mantles are commonly observed through infrared spectroscopy toward dense clouds, cloud cores, protostellar envelopes and protoplanetary disks. Up to 80% of the available oxygen, carbon and nitrogen are found in such ices; the most common ice constituents - H2O, CO2 and CO - are second in abundance only to H2 in many star forming regions. In addition to being a molecular reservoir, ice chemistry is responsible for much of the chemical evolution from H2O to complex, prebiotic molecules. Combining the existing ISO, Spitzer, VLT and Keck ice data results in a large sample of ice sources (\sime80) that span all stages of star formation and a large range of protostellar luminosities (<0.1-105 L\odot). Here we summarize the different techniques that have been applied to mine this ice data set on information on typical ice compositions in different environments and what this implies about how ices form and evolve during star and planet formation. The focus is on how to maximize the use of empirical constraints from ice observations, followed by the application of information from experiments and models. This strategy is used to identify ice bands and to constrain which ices form early during cloud formation, which form later in the prestellar core and which require protostellar heat and/or UV radiation to form. The utility of statistical tests, survival analysis and ice maps is highlighted; the latter directly reveals that the prestellar ice formation takes place in two phases, associated with H2O and CO ice formation, respectively, and that most protostellar ice variation can be explained by differences in the prestellar CO ice formation stage. Finally, special attention is paid to the difficulty of observing complex ices directly and how gas observations, experiments and models help in constraining this ice chemistry stage.

preprint2011arXiv

The Spitzer ice legacy: Ice evolution from cores to protostars

Ices regulate much of the chemistry during star formation and account for up to 80% of the available oxygen and carbon. In this paper, we use the Spitzer c2d ice survey, complimented with data sets on ices in cloud cores and high-mass protostars, to determine standard ice abundances and to present a coherent picture of the evolution of ices during low- and high-mass star formation. The median ice composition H2O:CO:CO2:CH3OH:NH3:CH4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 toward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. In the low-mass sample, the ice abundances with respect to H2O of CH4, NH3, and the component of CO2 mixed with H2O typically vary by <25%, indicative of co-formation with H2O. In contrast, some CO and CO2 ice components, XCN and CH3OH vary by factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, consistent with its OCN- identification. The origin(s) of the different levels of ice abundance variations are constrained by comparing ice inventories toward different types of protostars and background stars, through ice mapping, analysis of cloud-to-cloud variations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase is driven by hydrogenation of atoms, which results in a H2O-dominated ice. At later prestellar times, CO freezes out and variations in CO freeze-out levels and the subsequent CO-based chemistry can explain most of the observed ice abundance variations. The last important ice evolution stage is thermal and UV processing around protostars, resulting in CO desorption, ice segregation and formation of complex organic molecules. The distribution of cometary ice abundances are consistent with with the idea that most cometary ices have a protostellar origin.

preprint2011arXiv

The structure and dynamics of molecular gas in planet-forming zones: A CRIRES spectro-astrometric survey

We present a spectro-astrometric survey of molecular gas in the inner regions of 16 protoplanetary disks using CRIRES, the high resolution infrared imaging spectrometer on the Very Large Telescope. Spectro-astrometry with CRIRES measures the spatial extent of line emission to sub-milliarcsecond precision, or <0.2 AU at the distance of the observed targets. The sample consists of gas-rich disks surrounding stars with spectral types ranging from K to A. The properties of the spectro-astrometric signals divide the sources into two distinct phenomenological classes: one that shows clear Keplerian astrometric spectra, and one in which the astrometric signatures are dominated by gas with strong non-Keplerian (radial) motions. Similarly to the near-infrared continuum emission, as determined by interferometry, we find that the size of the CO line emitting region in the Keplerian sources obeys a size-luminosity relation as $R_CO L_*^0.5. The non-Keplerian spectro-astrometric signatures are likely indicative of the presence of wide-angle disk winds. The central feature of the winds is a strong sub-Keplerian velocity field due to conservation of angular momentum as the wind pressure drives the gas outwards. We construct a parametrized 2-dimensional disk+wind model that reproduces the observed characteristics the observed CO spectra and astrometry. The modeled winds indicate mass-loss rates of >10^-10 to 10^-8 Msol/yr. We suggest a unifying model in which all disks have slow molecular winds, but where the magnitude of the mass-loss rate determines the degree to which the mid-infrared molecular lines are dominated by the wind relative to the Keplerian disk surface.

preprint2010arXiv

A Spitzer c2d Legacy Survey to Identify and Characterize Disks with Inner Dust Holes

Understanding how disks dissipate is essential to studies of planet formation. However, identifying exactly how dust and gas dissipates is complicated due to difficulty in finding objects clearly in the transition of losing their surrounding material. We use Spitzer IRS spectra to examine 35 photometrically-selected candidate cold disks (disks with large inner dust holes). The infrared spectra are supplemented with optical spectra to determine stellar and accretion properties and 1.3mm photometry to measure disk masses. Based on detailed SED modeling, we identify 15 new cold disks. The remaining 20 objects have IRS spectra that are consistent with disks without holes, disks that are observed close to edge-on, or stars with background emission. Based on these results, we determine reliable criteria for identifying disks with inner holes from Spitzer photometry and examine criteria already in the literature. Applying these criteria to the c2d surveyed star-forming regions gives a frequency of such objects of at least 4% and most likely of order 12% of the YSO population identified by Spitzer. We also examine the properties of these new cold disks in combination with cold disks from the literature. Hole sizes in this sample are generally smaller than for previously discovered disks and reflect a distribution in better agreement with exoplanet orbit radii. We find correlations between hole size and both disk and stellar masses. Silicate features, including crystalline features, are present in the overwhelming majority of the sample although 10 micron feature strength above the continuum declines for holes with radii larger than ~7 AU. In contrast, PAHs are only detected in 2 out of 15 sources. Only a quarter of the cold disk sample shows no signs of accretion, making it unlikely that photoevaporation is the dominant hole forming process in most cases.

preprint2010arXiv

A Spitzer survey of mid-infrared molecular emission from protoplanetary disks I: Detection rates

We present a Spitzer InfraRed Spectrometer search for 10-36 micron molecular emission from a large sample of protoplanetary disks, including lines from H2O, OH, C2H2, HCN and CO2. This paper describes the sample and data processing and derives the detection rate of mid-infrared molecular emission as a function of stellar mass. The sample covers a range of spectral type from early M to A, and is supplemented by archival spectra of disks around A and B stars. It is drawn from a variety of nearby star forming regions, including Ophiuchus, Lupus and Chamaeleon. In total, we identify 22 T Tauri stars with strong mid-infrared H2O emission. Integrated water line luminosities, where water vapor is detected, range from 5x10^-4 to 9x10^-3 Lsun, likely making water the dominant line coolant of inner disk surfaces in classical T Tauri stars. None of the 5 transitional disks in the sample show detectable gaseous molecular emission with Spitzer upper limits at the 1% level in terms of line-to-continuum ratios (apart from H2). We find a strong dependence on detection rate with spectral type; no disks around our sample of 25 A and B stars were found to exhibit water emission, down to 1-2% line-to-continuum ratios, in the mid-infrared, while almost 2/3 of the disks around K stars show sufficiently intense water emission to be detected by Spitzer. Some Herbig Ae/Be stars show tentative H2O/OH emission features beyond 20 micron at the 1-2 level, however, and one of them shows CO2 in emission. We argue that the observed differences between T Tauri disks and Herbig Ae/Be disks is due to a difference in excitation and/or chemistry depending on spectral type and suggest that photochemistry may be playing an important role in the observable characteristics of mid-infrared molecular line emission from protoplanetary disks.

preprint2010arXiv

c2d Spitzer IRS spectra of embedded low-mass young stars: gas-phase emission lines

A survey of mid-IR gas-phase emission lines of H2, H2O and various atoms toward a sample of 43 embedded low-mass young stars in nearby star-forming regions is presented. The sources are selected from the Spitzer "Cores to Disks" (c2d) legacy program. The environment of embedded protostars is complex both in its physical structure (envelopes, outflows, jets, protostellar disks) and the physical processes (accretion, irradiation by UV and/or X-rays, excitation through slow and fast shocks) which take place. A key point is to spatially resolve the emission in the Spitzer-IRS spectra. An optimal extraction method is used to separate both spatially unresolved (compact, up to a few 100 AU) and spatially resolved (extended, 1000 AU or more) emission from the IRS spectra. The results are compared with the c2d disk sample and literature PDR and shock models to address the physical nature of the sources. Both compact and extended emission features are observed. Warm (Tex few 100 K) H2, observed through the pure rotational H2 S(0), S(1) and S(2) lines, and [S I] 25 mu emission is observed primarily in the extended component. [S I] is observed uniquely toward truly embedded sources and not toward disks. On the other hand hot (Tex>700 K) H2, observed primarily through the S(4) line, and [Ne II] emission is seen mostly in the compact component. [Fe II] and [Si II] lines are observed in both spatial components. Hot H2O emission is found in the compact component of some sources. The observed emission on >=1000 AU scales is characteristic of PDR emission and likely originates in the outflow cavities in the remnant envelope created by the stellar wind and jets from the embedded young stars. Weak shocks along the outflow wall can also contribute. The compact emission is likely of mixed origin, comprised of optically thick circumstellar disk and/or jet/outflow emission from the protostellar object.

preprint2010arXiv

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

preprint2010arXiv

Investigating the Cosmic-Ray Ionization Rate Near the Supernova Remnant IC 443 Through H3+ Observations

Observational and theoretical evidence suggests that high-energy Galactic cosmic rays are primarily accelerated by supernova remnants. If also true for low-energy cosmic rays, the ionization rate near a supernova remnant should be higher than in the general Galactic interstellar medium (ISM). We have searched for H3+ absorption features in 6 sight lines which pass through molecular material near IC 443---a well-studied case of a supernova remnant interacting with its surrounding molecular material---for the purpose of inferring the cosmic-ray ionization rate in the region. In 2 of the sight lines (toward ALS 8828 and HD 254577) we find large H3+ column densities, N(H3+)~3*10^14 cm^-2, and deduce ionization rates of zeta_2~2*10^-15 s^-1, about 5 times larger than inferred toward average diffuse molecular cloud sight lines. However, the 3 sigma upper limits found for the other 4 sight lines are consistent with typical Galactic values. This wide range of ionization rates is likely the result of particle acceleration and propagation effects, which predict that the cosmic-ray spectrum and thus ionization rate should vary in and around the remnant. While we cannot determine if the H3+ absorption arises in post-shock (interior) or pre-shock (exterior) gas, the large inferred ionization rates suggest that IC 443 is in fact accelerating a large population of low-energy cosmic rays. Still, it is unclear whether this population can propagate far enough into the ISM to account for the ionization rate inferred in diffuse Galactic sight lines.

preprint2010arXiv

Spectrally resolved pure rotational lines of water in protoplanetary disks

We present ground-based high resolution N-band spectra (Δv = 15 km/s) of pure rotational lines of water vapor in two protoplanetary disks surrounding the pre-main sequence stars AS 205N and RNO 90, selected based on detections of rotational water lines by the Spitzer IRS. Using VISIR on the Very Large Telescope, we spectrally resolve individual lines and show that they have widths of 30-60 km/s, consistent with an origin in Keplerian disks at radii of ~1 AU. The water lines have similar widths to those of the CO at 4.67 micron, indicating that the mid-infrared water lines trace similar radii. The rotational temperatures of the water are 540 and 600K in the two disks, respectively. However, the lines ratios show evidence of non-LTE excitation, with low-excitation line fluxes being over-predicted by 2-dimensional disk LTE models. Due to the limited number of observed lines and the non-LTE line ratios, an accurate measure of the water ortho/para ratio is not available, but a best estimate for AS 205N is ortho/para = 4.5 +/- 1.0, apparently ruling out a low-temperature origin of the water. The spectra demonstrate that high resolution spectroscopy of rotational water lines is feasible from the ground, and further that ground-based high resolution spectroscopy is likely to significantly improve our understanding of the inner disk chemistry recently revealed by recent Spitzer observations.

preprint2010arXiv

The c2d Spitzer Spectroscopic Survey of Ices Around Low-Mass Young Stellar Objects. IV. NH3 and CH3OH

NH3 and CH3OH are key molecules in astrochemical networks leading to the formation of more complex N- and O-bearing molecules, such as CH3CN and HCOOCH3. Despite a number of recent studies, little is known about their abundances in the solid state. (...) In this work, we investigate the ~ 8-10 micron region in the Spitzer IRS (InfraRed Spectrograph) spectra of 41 low-mass young stellar objects (YSOs). These data are part of a survey of interstellar ices in a sample of low-mass YSOs studied in earlier papers in this series. We used both an empirical and a local continuum method to correct for the contribution from the 10 micron silicate absorption in the recorded spectra. In addition, we conducted a systematic laboratory study of NH3- and CH3OH-containing ices to help interpret the astronomical spectra. We clearly detect a feature at ~9 micron in 24 low-mass YSOs. Within the uncertainty in continuum determination, we identify this feature with the NH3 nu_2 umbrella mode, and derive abundances with respect to water between ~2 and 15%. Simultaneously, we also revisited the case of CH3OH ice by studying the nu_4 C-O stretch mode of this molecule at ~9.7 micron in 16 objects, yielding abundances consistent with those derived by Boogert et al. 2008 (hereafter paper I) based on a simultaneous 9.75 and 3.53 micron data analysis. Our study indicates that NH3 is present primarily in H2O-rich ices, but that in some cases, such ices are insufficient to explain the observed narrow FWHM. The laboratory data point to CH3OH being in an almost pure methanol ice, or mixed mainly with CO or CO2, consistent with its formation through hydrogenation on grains. Finally, we use our derived NH3 abundances in combination with previously published abundances of other solid N-bearing species to find that up to 10-20 % of nitrogen is locked up in known ices.

preprint2009arXiv

A new raytracer for modeling AU-scale imaging of lines from protoplanetary disks

The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.

preprint2009arXiv

High Resolution Spectroscopy of [NeII] Emission from AA Tau and GM Aur

We present high resolution (R=80,000) spectroscopy of [NeII] emission from two young stars, GM Aur and AA Tau, which have moderate to high inclinations. The emission from both sources appears centered near the stellar velocity and is broader than the [NeII] emission measured previously for the face-on disk system TW Hya. These properties are consistent with a disk origin for the [NeII] emission we detect, with disk rotation (rather than photoevaporation or turbulence in a hot disk atmosphere) playing the dominant role in the origin of the line width. In the non-face-on systems, the [NeII] emission is narrower than the CO fundamental emission from the same sources. If the widths of both diagnostics are dominated by Keplerian rotation, this suggests that the [NeII] emission arises from larger disk radii on average than does the CO emission. The equivalent width of the [NeII] emission we detect is less than that of the spectrally unresolved [NeII] feature in the Spitzer spectra of the same sources. Variability in the [NeII] emission or the mid-infrared continuum, a spatially extended [NeII] component, or a very (spectrally) broad [NeII] component might account for the difference in the equivalent widths.

preprint2008arXiv

The c2d Spitzer spectroscopy survey of ices around low-mass young stellar objects, III: CH4

CH4 is proposed to be the starting point of a rich organic chemistry. Solid CH4 abundances have previously been determined mostly toward high mass star forming regions. Spitzer/IRS now provides a unique opportunity to probe solid CH4 toward low mass star forming regions as well. Infrared spectra from the Spitzer Space Telescope are presented to determine the solid CH4 abundance toward a large sample of low mass young stellar objects. 25 out of 52 ice sources in the $c2d$ (cores to disks) legacy have an absorption feature at 7.7 um, attributed to the bending mode of solid CH4. The solid CH4 / H2O abundances are 2-8%, except for three sources with abundances as high as 11-13%. These latter sources have relatively large uncertainties due to small total ice column densities. Toward sources with H2O column densities above 2E18 cm-2, the CH4 abundances (20 out of 25) are nearly constant at 4.7+/-1.6%. Correlation plots with solid H2O, CH3OH, CO2 and CO column densities and abundances relative to H2O reveal a closer relationship of solid CH4 with CO2 and H2O than with solid CO and CH3OH. The inferred solid CH4 abundances are consistent with models where CH4 is formed through sequential hydrogenation of C on grain surfaces. Finally the equal or higher abundances toward low mass young stellar objects compared with high mass objects and the correlation studies support this formation pathway as well, but not the two competing theories: formation from CH3OH and formation in gas phase with subsequent freeze-out.

preprint2008arXiv

The Spitzer c2d Legacy Results: Star Formation Rates and Efficiencies; Evolution and Lifetimes

(Abridged) The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instruments for five large, nearby molecular clouds. This paper combines information drawn from studies of individual clouds into a combined and updated statistical analysis of star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. Current star formation efficiencies range from 3% to 6%. Taken together, the five clouds are producing about 260 solar masses of stars per Myr. The star formation surface density is more than an order of magnitude larger than would be predicted from the Kennicutt relation used in extragalactic studies. Measured against the dense gas probed by the maps of dust continuum emission, the efficiencies are much higher, and the current stock of dense cores would be exhausted in 1.8 Myr on average. The derived lifetime for the Class I phase is 0.44 to 0.54 Myr, considerably longer than some estimates. Similarly, the lifetime for the Class 0 SED class, 0.10 to 0.16 Myr, is longer than early estimates. The great majority (90%) of young stars lie within loose clusters with at least 35 members and a stellar density of 1 solar mass per cubic pc. Accretion at the sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star of about 0.25 solar masses, given an efficiency of 0.3. Our data confirm and aggravate the "luminosity problem" for protostars. Our results strongly suggest that accretion is time variable, with prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight most luminous objects.

preprint2007arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds VIII. Serpens Observed with MIPS

We present maps of 1.5 square degrees of the Serpens dark cloud at 24, 70, and 160\micron observed with the Spitzer Space Telescope MIPS Camera. More than 2400 compact sources have been extracted at 24um, nearly 100 at 70um, and 4 at 160um. We estimate completeness limits for our 24um survey from Monte Carlo tests with artificial sources inserted into the Spitzer maps. We compare source counts, colors, and magnitudes in the Serpens cloud to two reference data sets, a 0.50 deg^2 set on a low-extinction region near the dark cloud, and a 5.3 deg^2 subset of the SWIRE ELAIS N1 data that was processed through our pipeline. These results show that there is an easily identifiable population of young stellar object candidates in the Serpens Cloud that is not present in either of the reference data sets. We also show a comparison of visual extinction and cool dust emission illustrating a close correlation between the two, and find that the most embedded YSO candidates are located in the areas of highest visual extinction.

preprint2006arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. III. Perseus Observed with IRAC

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

preprint2006arXiv

The Spitzer c2d survey of Large, Nearby, Interstellar Clouds. V. Chamaeleon II Observed with IRAC

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

preprint2003arXiv

Interferometric Observations of Formaldehyde in the Protoplanetary Disk around LkCa15

Emission from the $2_{12}-1_{11}$ line of H$_2$CO has been detected and marginally resolved toward LkCa15 by the Nobeyama Millimeter Array. The column density of H$_2$CO is higher than that observed in DM Tau and than predicted by theoretical models of disk chemistry; also the line-intensity profile is less centrally peaked than that for CO. A similar behavior is observed in other organic gaseous molecules in the LkCa 15 disk.

Geoffrey A. Blake

What is connected

Connect this record

See the researcher in context

Building this map preview

62 published item(s)

Protoplanetary disk cavities with JWST-MIRI: a dichotomy in molecular emission

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

New Constraints on Protoplanetary Disk Gas Masses in Lupus

Retrieving the C and O Abundances of HR 7672~AB: a Solar-Type Primary Star with a Benchmark Brown Dwarf

Mapping physical parameters in Orion KL at high spatial resolution

An ALMA Survey of H$_2$CO in Protoplanetary Disks

An evolutionary study of volatile chemistry in protoplanetary disks

An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks

Carbon-grain sublimation: a new top-down component of protostellar chemistry

Contrast and Temperature Dependence of Multi-Epoch High-Resolution Cross-Correlation Exoplanet Spectroscopy

Simulating the Multi-Epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-Transiting Hot Jupiter HD187123B

Sum-Frequency Signals in 2D-Terahertz-Terahertz-Raman Spectroscopy

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

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

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

CSO and CARMA Observations of L1157. II. Chemical Complexity in the Shocked Outflow

Discovery of the Interstellar Chiral Molecule Propylene Oxide (CH$_3$CHCH$_2$O)

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

Hydrocarbon emission rings in protoplanetary disks induced by dust evolution

On the Commonality of 10-30AU Sized Axisymmetric Dust Structures in Protoplanetary Disks

The Radial Distribution of H2 and CO in TW Hya as Revealed by Resolved ALMA Observations of CO Isotopologues

THz Time-Domain Spectroscopy of Mixed CO2-CH3OH Interstellar Ice Analogs

CSO and CARMA Observations of L1157. I. A Deep Search for Hydroxylamine (NH$_2$OH)

Detection of water vapor in the terrestrial planet forming region of a transition disk

Dimming and CO absorption toward the AA Tau protoplanetary disk: An infalling flow caused by disk instability?

Evidence of fast pebble growth near condensation fronts in the HL Tau protoplanetary disk

Herschel observations of EXtra-Ordinary Sources: Analysis of the HIFI 1.2 THz Wide Spectral Survey Toward Orion KL II. Chemical Implications

Measurements of water surface snow lines in classical protoplanetary disks

Nonlinear TeraHertz Coherent Excitation of Vibrational Modes of Liquids

Tracing the Ingredients for a Habitable Earth from Interstellar Space through Planet Formation

A CSO Search for $l$-C$_3$H$^+$: Detection in the Orion Bar PDR

ALMA Observations of the T Tauri Binary System AS 205: Evidence for Molecular Winds and/or Binary Interactions

Comparison of the dust and gas radial structure in the transition disk [PZ99] J160421.7-213028

Near-IR Direct Detection of Water Vapor in Tau Boo b

The Search for a Complex Molecule in a Selected Hot Core Region: A Rigorous Attempt to Confirm trans-Ethyl Methyl Ether toward W51 e1/e2

A Search for l-C3H+ and l-C3H in Sgr B2(N), Sgr B2(OH), and the Dark Cloud TMC-1

A Significantly Low CO Abundance Toward the TW Hya Protoplanetary Disk: A Path to Active Carbon Chemistry?

An observational investigation of the identity of B11244 ($l$-C$_3$H$^+$/C$_3$H$^-$)

Imaging of the CO Snow Line in a Solar Nebula Analog

Measuring Protoplanetary Disk Accretion with H I Pf-beta

VLT-CRIRES Survey of Rovibrational CO Emission from Protoplanetary Disks

Interstellar Carbodiimide (HNCNH) - A New Astronomical Detection from the GBT PRIMOS Survey via Maser Emission Features

CO rovibrational emission as a probe of inner disk structure

Detection of the Water Reservoir in a Forming Planetary System

Ices in starless and starforming cores

The Spitzer ice legacy: Ice evolution from cores to protostars

The structure and dynamics of molecular gas in planet-forming zones: A CRIRES spectro-astrometric survey

A Spitzer c2d Legacy Survey to Identify and Characterize Disks with Inner Dust Holes

A Spitzer survey of mid-infrared molecular emission from protoplanetary disks I: Detection rates

c2d Spitzer IRS spectra of embedded low-mass young stars: gas-phase emission lines

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

Investigating the Cosmic-Ray Ionization Rate Near the Supernova Remnant IC 443 Through H3+ Observations

Spectrally resolved pure rotational lines of water in protoplanetary disks

The c2d Spitzer Spectroscopic Survey of Ices Around Low-Mass Young Stellar Objects. IV. NH3 and CH3OH

A new raytracer for modeling AU-scale imaging of lines from protoplanetary disks

High Resolution Spectroscopy of [NeII] Emission from AA Tau and GM Aur

The c2d Spitzer spectroscopy survey of ices around low-mass young stellar objects, III: CH4

The Spitzer c2d Legacy Results: Star Formation Rates and Efficiencies; Evolution and Lifetimes

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds VIII. Serpens Observed with MIPS

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. III. Perseus Observed with IRAC

The Spitzer c2d survey of Large, Nearby, Interstellar Clouds. V. Chamaeleon II Observed with IRAC

Interferometric Observations of Formaldehyde in the Protoplanetary Disk around LkCa15