Source author record

Jonathan P. Williams

Jonathan P. Williams appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

64works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

64 published item(s)

preprint2026arXiv

pintervals: an R package for model-agnostic prediction intervals

The \pkg{pintervals} package aims to provide a unified framework for constructing prediction intervals and calibrating predictions in a model-agnostic setting using set-aside calibration data. It comprises routines to construct conformal as well as parametric and bootstrapped prediction intervals from any model that outputs point predictions. Several R packages and functions already exist for constructing prediction intervals, but they often focus on specific modeling frameworks or types of predictions, or require manual customization for different models or applications. By providing a consistent interface for a variety of prediction interval construction approaches (all model-agnostic), \pkg{pintervals} allows researchers to apply and compare them across different modeling frameworks and applications.

preprint2025arXiv

Bin-Conditional Conformal Prediction of Fatalities from Armed Conflict

Forecasting armed conflicts is a critical area of research with the potential to save lives and mitigate suffering. While existing forecasting models offer valuable point predictions, they often lack individual-level uncertainty estimates, limiting their usefulness for decision-making. Several approaches exist to estimate uncertainty, such as parametric and Bayesian prediction intervals, bootstrapping, quantile regression, but these methods often rely on restrictive assumptions, struggle to provide well-calibrated intervals across the full range of outcomes, or are computationally intensive. Conformal prediction offers a model-agnostic alternative that guarantees a user-specified level of coverage but typically provides only marginal coverage, potentially resulting in non-uniform coverage across different regions of the outcome space. In this paper, we introduce a novel extension called bin-conditional conformal prediction (BCCP), which enhances standard conformal prediction by ensuring consistent coverage rates across user-defined subsets (bins) of the outcome variable. We apply BCCP to simulated data as well as the forecasting of fatalities from armed conflicts, and demonstrate that it provides well-calibrated uncertainty estimates across various ranges of the outcome. Compared to standard conformal prediction, BCCP offers improved local coverage, though this comes at the cost of slightly wider prediction intervals.

preprint2023arXiv

Word Embeddings as Statistical Estimators

Word embeddings are a fundamental tool in natural language processing. Currently, word embedding methods are evaluated on the basis of empirical performance on benchmark data sets, and there is a lack of rigorous understanding of their theoretical properties. This paper studies word embeddings from a statistical theoretical perspective, which is essential for formal inference and uncertainty quantification. We propose a copula-based statistical model for text data and show that under this model, the now-classical Word2Vec method can be interpreted as a statistical estimation method for estimating the theoretical pointwise mutual information (PMI). Next, by building on the work of Levy and Goldberg (2014), we develop a missing value-based estimator as a statistically tractable and interpretable alternative to the Word2Vec approach. The estimation error of this estimator is comparable to Word2Vec and improves upon the truncation-based method proposed by Levy and Goldberg (2014). The proposed estimator also performs comparably to Word2Vec in a benchmark sentiment analysis task on the IMDb Movie Reviews data set.

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

preprint2022arXiv

An APEX search for carbon emission from NGC 1977 proplyds

We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$σ$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments.

preprint2022arXiv

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales.

preprint2022arXiv

High-resolution ALMA observations of transition disk candidates in Lupus

Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged)

preprint2022arXiv

The ODYSSEUS Survey. Motivation and First Results: Accretion, Ejection, and Disk Irradiation of CVSO 109

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

preprint2021arXiv

The ALMA Early Science view of FUor/EXor objects. I. Through the looking-glass of V2775 Ori

As part of an ALMA survey to study the origin of episodic accretion in young eruptive variables, we have observed the circumstellar environment of the star V2775 Ori. This object is a very young, pre-main sequence object which displays a large amplitude outburst characteristic of the FUor class. We present Cycle-2 band 6 observations of V2775 Ori with a continuum and CO (2-1) isotopologue resolution of 0.25\as (103 au). We report the detection of a marginally resolved circumstellar disc in the ALMA continuum with an integrated flux of $106 \pm 2$ mJy, characteristic radius of $\sim$ 30 au, inclination of $14.0^{+7.8}_{-14.5}$ deg, and is oriented nearly face-on with respect to the plane of the sky. The \co~emission is separated into distinct blue and red-shifted regions that appear to be rings or shells of expanding material from quasi-episodic outbursts. The system is oriented in such a way that the disc is seen through the outflow remnant of V2775 Ori, which has an axis along our line-of-sight. The $^{13}$CO emission displays similar structure to that of the \co, while the C$^{18}$O line emission is very weak. We calculated the expansion velocities of the low- and medium-density material with respect to the disc to be of -2.85 km s$^{-1}$ (blue), 4.4 km s$^{-1}$ (red) and -1.35 and 1.15 km s$^{-1}$ (for blue and red) and we derived the mass, momentum and kinetic energy of the expanding gas. The outflow has an hourglass shape where the cavities are not seen. We interpret the shapes that the gas traces as cavities excavated by an ancient outflow. We report a detection of line emission from the circumstellar disc and derive a lower limit of the gas mass of 3 \MJup.

preprint2021arXiv

The Ophiuchus DIsc Survey Employing ALMA (ODISEA). II. The effect of stellar multiplicity on disc properties

We present Adaptive Optics (AO) near infrared (NIR) observations using VLT/NACO and Keck/NIRC2 of ODISEA targets. ODISEA is an ALMA survey of the entire population of circumstellar discs in the Ophiuchus molecular cloud. From the whole sample of ODISEA we select all the discs that are not already observed in the NIR with AO and that are observable with NACO or NIRC2. The NIR-ODISEA survey consists of 147 stars observed in NIR AO imaging for the first time, as well as revisiting almost all the binary systems of Ophiuchus present in the literature (20 out of 21). In total, we detect 20 new binary systems and one triple system. For each of them we calculate the projected separation and position angle of the companion, as well as their NIR and millimeter flux ratios. From the NIR contrast we derived the masses of the secondaries, finding that 9 of them are in the sub-stellar regime (30-50 \MJup). Discs in multiple systems reach a maximum total dust mass of $\sim$ 50 M$_{\oplus}$, while discs in single stars can reach a dust mass of 200 M$_{\oplus}$. Discs with masses above 10 M$_{\oplus}$ are found only around binaries with projected separations larger than $\sim$ 110 au. The maximum disc size is also larger around single star than binaries. However, since most discs in Ophiuchus are very small and low-mass, the effect of visual binaries is relatively weak in the general disc population.

preprint2020arXiv

ALMA Observations of Young Eruptive Stars: continuum disk sizes and molecular outflows

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of four young, eruptive star-disk systems at 0.4" resolution: two FUors (V582 Aur and V900 Mon), one EXor (UZ Tau E) and one source with an ambiguous FU/EXor classification (GM Cha). The disks around GM Cha, V900 Mon and UZ Tau E are resolved. These observations increase the sample of FU/EXors observed at sub-arcsecond resolution by 15%. The disk sizes and masses of FU/EXors objects observed by ALMA so far suggest that FUor disks are more massive than Class 0/I disks in Orion and Class II disks in Lupus of similar size. EXor disks in contrast do not seem to be distinguishable from these two populations. We reach similar conclusions when comparing the FU/EXor sample to the Class I and Class II disks in Ophiuchus. FUor disks around binaries are host to more compact disks than those in single-star systems, similar to non-eruptive young disks. We detect a wide-angle outflow around GM Cha in $^{12}$CO emission, wider than typical Class I objects and more similar to those found around some FUor objects. We use radiative transfer models to fit the continuum and line data of the well-studied disk around UZ Tau E. The line data is well described by a keplerian disk, with no evidence of outflow activity (similar to other EXors). The detection of wide-angle outflows in FUors and not in EXors support to the current picture in which FUors are more likely to represent an accretion burst in the protostellar phase (Class I), while EXors are smaller accretion events in the protoplanetary (Class II) phase.

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

Keck/NIRC2 $L$'-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

preprint2020arXiv

Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup

While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3''$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like'' structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup's peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup's extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems.

preprint2020arXiv

Protoplanetary disk masses in NGC 2024: Evidence for two populations

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.

preprint2019arXiv

Protostellar Evolution in Serpens Main: Possible Origin of Disk-Size Diversity

We have observed the submillimeter continuum condensations SMM2, SMM4, SMM9, and SMM11 in the star forming cluster Serpens Main using the Atacama Large Millimeter/submillimeter Array during Cycle 3 in the 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at an angular resolution of ~0.55 (240 au). Sixteen sources have been detected in the 1.3 mm continuum, which can be classified into three groups. Group 1 consists of six sources showing extended continuum emission and bipolar/monopolar 12CO outflows. Although all the Group 1 members are classified as Class 0 protostars, our observations suggest evolutionary trends among them in terms of 12CO outflow dynamical time, SO emission distribution, C18O fractional abundance, and continuum morphology. Group 2 consists of four sources associated with a continuum filamentary structure and no 12CO outflows. Central densities estimated from the 1.3 mm continuum intensity suggest that they are prestellar sources in a marginally Jeans unstable state. Group 3 consists of six Spitzer sources showing point-like 1.3 mm continuum emission and clumpy 12CO outflows. These features of Group 3 suggest envelope dissipation, preventing disk growth from the present size, r <~ 60 au. The Group 3 members are protostars that may be precursors to the T Tauri stars associated with small disks at tens-au radii identified in recent surveys.

preprint2016arXiv

A VLA Survey For Faint Compact Radio Sources in the Orion Nebula Cluster

We present Karl G. Janksy Very Large Array (VLA) 1.3 cm, 3.6 cm, and 6 cm continuum maps of compact radio sources in the Orion Nebular Cluster. We mosaicked 34 square arcminutes at 1.3 cm, 70 square arcminutes at 3.6 cm and 109 square arcminutes at 6 cm, containing 778 near-infrared detected YSOs and 190 HST-identified proplyds (with significant overlap between those characterizations). We detected radio emission from 175 compact radio sources in the ONC, including 26 sources that were detected for the first time at these wavelengths. For each detected source we fit a simple free-free and dust emission model to characterize the radio emission. We extrapolate the free-free emission spectrum model for each source to ALMA bands to illustrate how these measurements could be used to correctly measure protoplanetary disk dust masses from sub-millimeter flux measurements. Finally, we compare the fluxes measured in this survey with previously measured fluxes for our targets, as well as four separate epochs of 1.3 cm data, to search for and quantify variability of our sources.

preprint2016arXiv

Imaging the water-snow line during a protostellar outburst

A snow-line is the region of a protoplanetary disk at which a major volatile, such as water or carbon monoxide, reaches its condensation temperature. Snow-lines play a crucial role in disk evolution by promoting the rapid growth of ice-covered grains. Signatures of the carbon monoxide snow-line (at temperatures of around 20 kelvin) have recently been imaged in the disks surrounding the pre-main-sequence stars TW Hydra and HD163296, at distances of about 30 astronomical units (au) from the star. But the water snow-line of a protoplanetary disk (at temperatures of more than 100 kelvin) has not hitherto been seen, as it generally lies very close to the star (less than 5 au away for solar-type stars). Water-ice is important because it regulates the efficiency of dust and planetesimal coagulation, and the formation of comets, ice giants and the cores of gas giants. Here we report ALMA images at 0.03-arcsec resolution (12 au) of the protoplanetary disk around V883 Ori, a protostar of 1.3 solar masses that is undergoing an outburst in luminosity arising from a temporary increase in the accretion rate. We find an intensity break corresponding to an abrupt change in the optical depth at about 42 au, where the elevated disk temperature approaches the condensation point of water, from which we conclude that the outburst has moved the water snow-line. The spectral behaviour across the snow-line confirms recent model predictions: dust fragmentation and the inhibition of grain growth at higher temperatures results in soaring grain number densities and optical depths. As most planetary systems are expected to experience outbursts caused by accretion during their formation our results imply that highly dynamical water snow-lines must be considered when developing models of disk evolution and planet formation.

preprint2016arXiv

Measuring protoplanetary disk gas surface density profiles with ALMA

The gas and dust are spatially segregated in protoplanetary disks due to the vertical settling and radial drift of large grains. A fuller accounting of the mass content and distribution in disks therefore requires spectral line observations. We extend the modeling approach presented in Williams & Best (2014) to show that gas surface density profiles can be measured from high fidelity 13CO integrated intensity images. We demonstrate the methodology by fitting ALMA observations of the HD 163296 disk to determine a gas mass, Mgas = 0.048 solar masse, and accretion disk characteristic size Rc = 213au and gradient gamma = 0.39. The same parameters match the C18O 2--1 image and indicates an abundance ratio [13CO]/[C18O] of 700 independent of radius. To test how well this methodology can be applied to future line surveys of smaller, lower mass T Tauri disks, we create a large 13CO 2--1 image library and fit simulated data. For disks with gas masses 3-10 Jupiter masses at 150pc, ALMA observations with a resolution of 0.2-0.3 arcseconds and integration times of about 20 minutes allow reliable estimates of Rc to within about 10au and gamma to within about 0.2. Economic gas imaging surveys are therefore feasible and offer the opportunity to open up a new dimension for studying disk structure and its evolution toward planet formation.

preprint2015arXiv

A SCUBA-2 850-$μ$m survey of circumstellar disks in the $λ$ Orionis cluster

We present results from an 850-$μ$m survey of the $\sim$ 5 Myr old $λ$ Orionis star-forming region. We used the SCUBA-2 camera on the James Clerk Maxwell Telescope to survey a $\sim$0.5-diameter circular region containing 36 (out of 59) cluster members with infrared excesses indicative of circumstellar disks. We detected only one object at $>3σ$ significance, the Herbig Ae star HD 245185, with a flux density of $\sim74$ mJy beam$^{-1}$ corresponding to a dust mass of $\sim150$ M$_{\oplus}$. Stacking the individually undetected sources did not produce a significant mean signal but gives an upper limit on the average dust mass for $λ$ Orionis disks of $\sim3$ M$_{\oplus}$. Our follow-up observations of HD 245185 with the Submillimeter Array found weak CO 2--1 line emission with an integrated flux of $\sim170$ mJy km s$^{-1}$ but no $^{13}$CO or C$^{18}$O isotopologue emission at 30 mJy km s$^{-1}$ sensitivity, suggesting a gas mass of $\lesssim1$ M$_{\rm Jup}$. The implied gas-to-dust ratio is thus $\gtrsim50$ times lower than the canonical interstellar medium value, setting HD 245185 apart from other Herbig Ae disks of similar age, which have been found to be gas rich; as HD 245185 also shows signs of accretion, we may be catching it in the final phases of disk clearing. Our study of the $λ$ Orionis cluster places quantitative constraints on planet formation timescales, indicating that at $\sim5$ Myr the average disk no longer has sufficient dust and gas to form giant planets and perhaps even super Earths; the bulk material has been mostly dispersed or is locked in pebbles/planetesimals larger than a few mm in size.

preprint2015arXiv

ALMA Observations of the Largest Proto-Planetary Disk in the Orion Nebula, 114-426: A CO Silhouette

We present ALMA observations of the largest protoplanetary disk in the Orion Nebula, 114-426. Detectable 345 GHz (856 micron) dust continuum is produced only in the 350 AU central region of the ~1000 AU diameter silhouette seen against the bright H-alpha background in HST images. Assuming optically thin dust emission at 345 GHz, a gas-to-dust ratio of 100, and a grain temperature of 20 K, the disk gas-mass is estimated to be 3.1 +/- 0.6 Jupiter masses. If most solids and ices have have been incorporated into large grains, however, this value is a lower limit. The disk is not detected in dense-gas tracers such as HCO+ J=4-3, HCN J=4-3, or CS =7-6. These results may indicate that the 114-426 disk is evolved and depleted in some light organic compounds found in molecular clouds. The CO J=3-2 line is seen in absorption against the bright 50 to 80 K background of the Orion A molecular cloud over the full spatial extent and a little beyond the dust continuum emission. The CO absorption reaches a depth of 27 K below the background CO emission at VLSR ~6.7 km/s about 0.52 arcseconds (210 AU) northeast and 12 K below the background CO emission at VLSR ~ 9.7 km/s about 0.34 arcseconds (140 AU) southwest of the suspected location of the central star, implying that the embedded star has a mass less than 1 Solar mass .

preprint2015arXiv

Molecule sublimation as a tracer of protostellar accretion: Evidence for accretion bursts from high angular resolution C18O images

The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such as intermittent bursts? We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current luminosities. Using observations from the Submillimeter Array we measure the extents of the emission from the C18O isotopologue toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current luminosities of the sources through dust and line radiative transfer calculations. Eight out of sixteen sources show more extended C18O emission than predicted by the models. The modeling shows that the likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. For four of those sources the increase would have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C18O only sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H2O in the grain ice-mantles. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the chemical structures of embedded protostars.

preprint2015arXiv

On the nature of the tertiary companion to FW Tau: ALMA CO observations and SED modeling

It is thought that planetary mass companions may form through gravitational disk instabilities or core accretion. Identifying such objects in the process of formation would provide the most direct test for the competing formation theories. One of the most promising candidates for a planetary mass object still in formation is the third object in the FW Tau system. We here present ALMA cycle 1 observations confirming the recently published 1.3 mm detection of a dust disk around this third object and present for the first time a clear detection of a single peak $^{12}$CO (2-1) line, providing direct evidence for the simultaneous existence of a gas disk. We perform radiative transfer modeling of the third object in FW Tau and find that current observations are consistent with either a brown dwarf embedded in an edge-on disk or a planet embedded in a low inclination disk, which is externally irradiated by the binary companion. Further observations with ALMA, aiming for high SNR detections of non contaminated gas lines, are required to conclusively unveil the nature of the third object in FW Tau.

preprint2015arXiv

Protoplanetary Disk Masses in the Young NGC 2024 Cluster

We present the results from a Submillimeter Array survey of the 887 micron continuum emission from the protoplanetary disks around 95 young stars in the young cluster NGC 2024. Emission was detected from 22 infrared sources, with flux densities from ~5 to 330 mJy; upper limits (at 3sigma) for the other 73 sources range from 3 to 24 mJy. For standard assumptions, the corresponding disk masses range from ~0.003 to 0.2Msolar, with upper limits at 0.002--0.01Msolar. The NGC 2024 sample has a slightly more populated tail at the high end of its disk mass distribution compared to other clusters, but without more information on the nature of the sample hosts it remains unclear if this difference is statistically significant or a superficial selection effect. Unlike in the Orion Trapezium, there is no evidence for a disk mass dependence on the (projected) separation from the massive star IRS2b in the NGC 2024 cluster. We suggest that this is due to either the cluster youth or a comparatively weaker photoionizing radiation field.

preprint2014arXiv

A Parametric Modeling Approach to Measuring the Gas Masses of Circumstellar Disks

The disks that surround young stars are mostly composed of molecular gas, which is harder to detect and interpret than the accompanying dust. Disk mass measurements have therefore relied on large and uncertain extrapolations from the dust to the gas. We have developed a grid of models to study the dependencies of isotopologue CO line strengths on disk structure and temperature parameters and find that a combination of 13CO and C18O observations provides a robust measure of the gas mass. We apply this technique to Submillimeter Array observations of nine circumstellar disks and published measurements of six well studied disks. We find evidence for selective photodissociation of C18O and determine masses to within a factor of about three. The inferred masses for the nine disks in our survey range from 0.7-6 M_Jup, and all are well below the extrapolation from the interstellar medium gas-to-dust ratio of 100. This is consistent with the low masses of planets found around such stars, and may be due to accretion or photoevaporation of a dust-poor upper atmosphere. However, the masses may be underestimated if there are more efficient CO depletion pathways than those known in molecular clouds and cold cores.

preprint2014arXiv

ALMA observations of a misaligned binary protoplanetary disk system in Orion

We present ALMA observations of a wide binary system in Orion, with projected separation 440 AU, in which we detect submillimeter emission from the protoplanetary disks around each star. Both disks appear moderately massive and have strong line emission in CO 3-2, HCO+ 4-3, and HCN 3-2. In addition, CS 7-6 is detected in one disk. The line-to-continuum ratios are similar for the two disks in each of the lines. From the resolved velocity gradients across each disk, we constrain the masses of the central stars, and show consistency with optical-infrared spectroscopy, both indicative of a high mass ratio ~9. The small difference between the systemic velocities indicates that the binary orbital plane is close to face-on. The angle between the projected disk rotation axes is very high, ~72 degrees, showing that the system did not form from a single massive disk or a rigidly rotating cloud core. This finding, which adds to related evidence from disk geometries in other systems, protostellar outflows, stellar rotation, and similar recent ALMA results, demonstrates that turbulence or dynamical interactions act on small scales well below that of molecular cores during the early stages of star formation.

preprint2014arXiv

ALMA Observations of the Orion Proplyds

We present ALMA observations of protoplanetary disks ("proplyds") in the Orion Nebula Cluster. We imaged 5 individual fields at 856um containing 22 HST-identified proplyds and detected 21 of them. Eight of those disks were detected for the first time at submillimeter wavelengths, including the most prominent, well-known proplyd in the entire Orion Nebula, 114-426. Thermal dust emission in excess of any free-free component was measured in all but one of the detected disks, and ranged between 1-163 mJy, with resulting disk masses of 0.3-79 Mjup. An additional 26 stars with no prior evidence of associated disks in HST observations were also imaged within the 5 fields, but only 2 were detected. The disk mass upper limits for the undetected targets, which include OB stars, theta1Ori C and theta1Ori F, range from 0.1-0.6 Mjup. Combining these ALMA data with previous SMA observations, we find a lack of massive (>3 Mjup) disks in the extreme-UV dominated region of Orion, within 0.03 pc of O-star theta1Ori C. At larger separations from theta1Ori C, in the far-UV dominated region, there is a wide range of disk masses, similar to what is found in low-mass star forming regions. Taken together, these results suggest that a rapid dissipation of disk masses likely inhibits potential planet formation in the extreme-UV dominated regions of OB associations, but leaves disks in the far-UV dominated regions relatively unaffected.

preprint2014arXiv

Molecular Outflows Driven by Low-Mass Protostars. I. Correcting for Underestimates When Measuring Outflow Masses and Dynamical Properties

We present a survey of 28 molecular outflows driven by low-mass protostars, all of which are sufficiently isolated spatially and/or kinematically to fully separate into individual outflows. Using a combination of new and archival data from several single-dish telescopes, 17 outflows are mapped in CO (2-1) and 17 are mapped in CO (3-2), with 6 mapped in both transitions. For each outflow, we calculate and tabulate the mass, momentum, kinetic energy, mechanical luminosity, and force assuming optically thin emission in LTE at an excitation temperature of 50 K. We show that all of the calculated properties are underestimated when calculated under these assumptions. Taken together, the effects of opacity, outflow emission at low velocities confused with ambient cloud emission, and emission below the sensitivities of the observations increase outflow masses and dynamical properties by an order of magnitude, on average, and factors of 50-90 in the most extreme cases. Different (and non-uniform) excitation temperatures, inclination effects, and dissociation of molecular gas will all work to further increase outflow properties. Molecular outflows are thus almost certainly more massive and energetic than commonly reported. Additionally, outflow properties are lower, on average, by almost an order of magnitude when calculated from the CO (3-2) maps compared to the CO (2-1) maps, even after accounting for different opacities, map sensitivities, and possible excitation temperature variations. It has recently been argued in the literature that the CO (3-2) line is subthermally excited in outflows, and our results support this finding.

preprint2014arXiv

SMA observations on faint Submillimeter Galaxies with S$_{850}$ < 2 mJy: Ultra Dusty Low-Luminosity Galaxies at High Redshift

We obtained SMA observations of eight faint (intrinsic 850 $μ$m fluxes $<$ 2 mJy) submillimeter galaxies (SMGs) discovered in SCUBA images of the massive lensing cluster fields A370, A2390, and A1689 and detected five. In total, we obtain 5 SMA detections, all of which have de-lensed fluxes $<$1 mJy with estimated total infrared luminosities 10$^{10}-10^{12}$ $L_\odot$, comparable to luminous infrared galaxies (LIRGs) and normal star-forming galaxies. Based on the latest number counts, these galaxies contribute $\sim$70% of the 850 $μ$m extragalactic background light and represent the dominant star-forming galaxy population in the dusty universe. However, only 40$^{+30}_{-16}$% of our faint SMGs would be detected in deep optical or near-infrared surveys, which suggests many of these sources are at high redshifts ($z \gtrsim 3$) or extremely dusty, and they are not included in current star formation history estimates.

preprint2013arXiv

A Herschel PACS survey of the dust and gas in Upper Scorpius disks

We present results of far-infrared photometric observations with Herschel PACS of a sample of Upper Scorpius stars, with a detection rate of previously known disk-bearing K and M stars at 70, 100, and 160 micron of 71%, 56%, and 50%, respectively. We fit power-law disk models to the spectral energy distributions of K & M stars with infrared excesses, and have found that while many disks extend in to the sublimation radius, the dust has settled to lower scale heights than in disks of the less evolved Taurus-Auriga population, and have much reduced dust masses. We also conducted Herschel PACS observations for far-infrared line emission and JCMT observations for millimeter CO lines. Among B and A stars, 0 of 5 debris disk hosts exhibit gas line emission, and among K and M stars, only 2 of 14 dusty disk hosts are detected. The OI 63 micron and CII 157 micron lines are detected toward [PZ99] J160421.7-213028 and [PBB2002] J161420.3-190648, which were found in millimeter photometry to host two of the most massive dust disks remaining in the region. Comparison of the OI line emission and 63 micron continuum to that of Taurus sources suggests the emission in the former source is dominated by the disk, while in the other there is a significant contribution from a jet. The low dust masses found by disk modeling and low number of gas line detections suggest that few stars in Upper Scorpius retain sufficient quantities of material for giant planet formation. By the age of Upper Scorpius, giant planet formation is essentially complete.

preprint2013arXiv

A SCUBA-2 850 micron survey of protoplanetary discs in the sigma Orionis cluster

We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution.

preprint2013arXiv

Characterisation of SCUBA-2 450um and 850um-selected Galaxies in the COSMOS Field

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

preprint2013arXiv

Herschel/PACS Survey of protoplanetary disks in Taurus/Auriga -- Observations of [OI] and [CII], and far infrared continuum

The Herschel Space Observatory was used to observe ~ 120 pre-main-sequence stars in Taurus as part of the GASPS Open Time Key project. PACS was used to measure the continuum as well as several gas tracers such as [OI] 63 μm, [OI] 145 μm, [CII] 158 μm, OH, H2O and CO. The strongest line seen is [OI] at 63 μm. We find a clear correlation between the strength of the [OI] 63 μm line and the 63 μm continuum for disk sources. In outflow sources, the line emission can be up to 20 times stronger than in disk sources, suggesting that the line emission is dominated by the outflow. The tight correlation seen for disk sources suggests that the emission arises from the inner disk ($<$ 50 AU) and lower surface layers of the disk where the gas and dust are coupled. The [OI] 63 μm is fainter in transitional stars than in normal Class II disks. Simple SED models indicate that the dust responsible for the continuum emission is colder in these disks, leading to weaker line emission. [CII] 158 μm emission is only detected in strong outflow sources. The observed line ratios of [OI] 63 μm to [OI] 145 μm are in the regime where we are insensitive to the gas-to-dust ratio, neither can we discriminate between shock or PDR emission. We detect no Class III object in [OI] 63 μm and only three in continuum, at least one of which is a candidate debris disk.

preprint2013arXiv

Outflows and Massive Stars in the protocluster IRAS 05358+3543

We present new near-IR H2, CO J=2-1, and CO J = 3-2 observations to study outflows in the massive star forming region IRAS 05358+3543. The Canada-France-Hawaii Telescope H2 images and James Clerk Maxwell Telescope CO data cubes of the IRAS 05358 region reveal several new outflows, most of which emerge from the dense cluster of sub-mm cores associated with the Sh 2-233IR NE cluster to the northeast of IRAS 05358. We used Apache Point Observatory (APO) JHK spectra to determine line of sight velocities of the outflowing material. Analysis of archival VLA cm continuum data and previously published VLBI observations reveal a massive star binary as a probable source of one or two of the outflows. We have identified probable sources for 6 outflows and candidate counterflows for 7 out of a total of 11 seen to be originating from the IRAS 05358 clusters. We classify the clumps within Sh 2-233IR NE as an early protocluster and Sh 2-233IR SW as a young cluster, and conclude that the outflow energy injection rate approximately matches the turbulent decay rate in Sh 2-233IR NE.

preprint2013arXiv

Resolving the Cosmic Far-Infrared Background at 450 and 850 Microns with SCUBA-2

We use the SCUBA-2 submillimeter camera mounted on the JCMT to obtain extremely deep number counts at 450 and 850um. We combine data on two cluster lensing fields, A1689 and A370, and three blank fields, CDF-N, CDF-S, and COSMOS, to measure the counts over a wide flux range at each wavelength. We use statistical fits to broken power law representations to determine the number counts. This allows us to probe to the deepest possible level in the data. At both wavelengths our results agree well with the literature in the flux range over which they have been measured, with the exception of the 850um counts in CDF-S, where we do not observe the counts deficit found by previous single-dish observations. At 450um, we detect significant counts down to ~1mJy, an unprecedented depth at this wavelength. By integrating the number counts above this flux limit, we measure 113.9^{+49.7}_{-28.4} Jydeg^{-2} of the 450um extragalactic background light (EBL). The majority of this contribution is from sources with S_450um between 1-10mJy, and these sources are likely to be the ones that are analogous to the local luminous infrared galaxies (LIRGs). At 850um, we measure 37.3^{+21.1}_{-12.9} Jydeg^{-2} of the EBL. Because of the large systematic uncertainties on the COBE measurements, the percentage of the EBL we resolve could range from 48%-153% (44%-178%) at 450 (850)um. Based on high-resolution SMA observations of around half of the 4sigma 850um sample in CDF-N, we find that 12.5^{+12.1}_{-6.8}% of the sources are blends of multiple fainter sources. This is a low multiple fraction, and we find no significant difference between our original SCUBA-2 850um counts and the multiplicity corrected counts.

preprint2013arXiv

Resolving The Moth at Millimeter Wavelengths

HD 61005, also known as "The Moth," is one of only a handful of debris disks that exhibit swept-back "wings" thought to be caused by interaction with the ambient interstellar medium (ISM). We present 1.3 mm Submillimeter Array (SMA) observations of the debris disk around HD 61005 at a spatial resolution of 1.9 arcsec that resolve the emission from large grains for the first time. The disk exhibits a double-peaked morphology at millimeter wavelengths, consistent with an optically thin ring viewed close to edge-on. To investigate the disk structure and the properties of the dust grains we simultaneously model the spatially resolved 1.3 mm visibilities and the unresolved spectral energy distribution. The temperatures indicated by the SED are consistent with expected temperatures for grains close to the blowout size located at radii commensurate with the millimeter and scattered light data. We also perform a visibility-domain analysis of the spatial distribution of millimeter-wavelength flux, incorporating constraints on the disk geometry from scattered light imaging, and find suggestive evidence of wavelength-dependent structure. The millimeter-wavelength emission apparently originates predominantly from the thin ring component rather than tracing the "wings" observed in scattered light. The implied segregation of large dust grains in the ring is consistent with an ISM-driven origin for the scattered light wings.

preprint2012arXiv

870 micron Imaging of a Transitional Disk in Upper Scorpius: Holdover from the Era of Giant Planet Formation?

We present 880 micron images of the transition disk around the star [PZ99] J160421.7-213028, a solar-mass star in the nearby Upper Scorpius association. With a resolution down to 0.34 arcsec, we resolve the inner hole in this disk, and via model fitting to the visibilities and spectral energy distribution we determine both the structure of the outer region and the presence of sparse dust within the cavity. The disk contains about 0.1 Jupiter masses of mm-emitting grains, with an inner disk edge of about 70 AU. The inner cavity contains a small amount of dust with a depleted surface density in a region extending from about 20-70 AU. Taking into account prior observations indicating little to no stellar accretion, the lack of a binary companion, and the presence of dust near 0.1 AU, we determine that the most likely mechanism for the formation of this inner hole is the presence of one or more giant planets.

preprint2012arXiv

Astronomical Evidence for the Rapid Growth of Millimeter Sized Particles in Protoplanetary Disks

I summarize recent surveys of protoplanetary disks at millimeter wavelengths and show that the distribution of luminosity, equivalent to the mass in small dust grains, declines rapidly. This contrasts with statistics on the lifetime of disks from infrared observations and the high occurrence of planets from radial velocity and transit surveys. I suggest that these disparate results can be reconciled if most of the dust in a disk is locked up in millimeter and larger sized particles within about 2 Myr. This statistical result on disk evolution agrees with detailed modeling of a small number of individual disks and with cosmochemical measurements of rapid planetesimal formation.

preprint2012arXiv

Faint Submillimter Galaxy Counts at 450 micron

We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge]

preprint2012arXiv

HII Regions, Embedded Protostars, and Starless Cores in Sharpless 2-157

We present arcsecond resolution 1.4mm observations of the high mass star forming region, Sharpless 2-157, that reveal the cool dust associated with the first stages of star formation. These data are compared with archival images at optical, infrared, and radio wavelengths, and complemented with new arcsecond resolution mid-infrared data. We identify a dusty young HII region, numerous infrared sources within the cluster envelope, and four starless condensations. Three of the cores lie in a line to the south of the cluster peak, but the most massive one is right at the center and associated with a jumble of bright radio and infrared sources. This presents an interesting juxtaposition of high and low mass star formation within the same cluster which we compare with similar observations of other high mass star forming regions and discuss in the context of cluster formation theory.

preprint2012arXiv

Submillimeter Array Observations of the RX J1633.9-2442 Transition Disk: Evidence for Multiple Planets in the Making

We present continuum high resolution Submillimeter Array (SMA) observations of the transition disk object RX J1633.9-2442, which is located in the Ophiuchus molecular cloud and has recently been identified as a likely site of ongoing giant planet formation. The observations were taken at 340 GHz (880 micron) with the SMA in its most extended configuration, resulting in an angular resolution of 0.3" (35 AU at the distance of the target). We find that the disk is highly inclined (i ~50 deg) and has an inner cavity ~25 AU in radius, which is clearly resolved by our observations. We simultaneously model the entire optical to millimeter wavelength spectral energy distribution (SED) and SMA visibilities of RX J1633.9-2442 in order to constrain the structure of its disk. We find that an empty cavity ~25 AU in radius is inconsistent with the excess emission observed at 12, 22, and 24 micron. Instead, the mid-IR excess can be modeled by either a narrow, optically thick ring at ~10 AU or an optically thin region extending from ~7 AU to ~25 AU. The inner disk (r < 5 AU) is mostly depleted of small dust grains as attested by the lack of detectable near-IR excess. We also present deep Keck aperture masking observations in the near-IR, which rule out the presence of a companion up to 500 times fainter than the primary star (in K-band) for projected separations in the 5-20 AU range. We argue that the complex structure of the RX J1633.9-2442 disk is best explained by multiple planets embedded within the disk. We also suggest that the properties and incidence of objects such as RX J1633.9-2442, T Cha, and LkCa 15 (and those of the companions recently identified to these two latter objects) are most consistent with the runaway gas accretion phase of the core accretion model, when giant planets gain their envelopes and suddenly become massive enough to open wide gaps in the disk.

preprint2012arXiv

The Nature of Transition Circumstellar Disks III. Perseus, Taurus, and Auriga

As part of an ongoing program aiming to characterize a large number of Spitzer-selected transition disks (disks with reduced levels of near-IR and/or mid- IR excess emission), we have obtained (sub)millimeter wavelength photometry, high-resolution optical spectroscopy, and adaptive optics near-infrared imaging for a sample of 31 transition objects located in the Perseus, Taurus, and Auriga molecular clouds. We use these ground-based data to estimate disk masses, multiplicity, and accretion rates in order to investigate the mechanisms potentially responsible for their inner holes. Following our previous studies in other regions, we combine disk masses, accretion rates and multiplicity data with other information, such as SED morphology and fractional disk luminosity to classify the disks as strong candidates for the following categories: grain-growth dominated disks (7 objects), giant planet-forming disks (6 objects), photoevaporating disks (7 objects), debris disks (11 objects), and cicumbinary disks (1 object, which was also classified as a photoeavaporating disk). Combining our sample of 31 transition disks with those from our previous studies results in a sample of 74 transition objects that have been selected, characterized, and classified in an homogenous way. We discuss this combined high-quality sample in the context of the current paradigm of the evolution and dissipation of protoplanetary disks and use its properties to constrain different aspects of the key processes driving their evolution. We find that the age distribution of disks that are likely to harbor recently formed giant planets favors core accretion as the main planet formation mechanism and a ~2-3 Myr formation timescale

preprint2011arXiv

Classifying structures in the ISM with Support Vector Machines: the G16.05-0.57 supernova remnant

We apply Support Vector Machines -- a machine learning algorithm -- to the task of classifying structures in the Interstellar Medium. As a case study, we present a position-position velocity data cube of 12 CO J=3--2 emission towards G16.05-0.57, a supernova remnant that lies behind the M17 molecular cloud. Despite the fact that these two objects partially overlap in position-position-velocity space, the two structures can easily be distinguished by eye based on their distinct morphologies. The Support Vector Machine algorithm is able to infer these morphological distinctions, and associate individual pixels with each object at >90% accuracy. This case study suggests that similar techniques may be applicable to classifying other structures in the ISM -- a task that has thus far proven difficult to automate.

preprint2011arXiv

Protoplanetary Disk Masses in IC348: A Rapid Decline in the Population of Small Dust Grains After 1 Myr

We present a 1.3 mm continuum survey of protoplanetary disks in the 2-3 Myr old cluster, IC348, with the Submillimeter Array. We observed 85 young stellar objects and detected 10 with 1.3 mm fluxes greater than 2 mJy. The brightest source is a young embedded protostar driving a molecular outflow. The other 9 detections are dusty disks around optically visible stars. Our millimeter flux measurements translate into total disk masses ranging from 2 to 6 Jupiter masses. Each detected disk has strong mid-infrared emission in excess of the stellar photosphere and has Halpha equivalent widths larger than the average in the cluster and indicative of ongoing gas accretion. The disk mass distribution, however, is shifted by about a factor of 20 to lower masses, compared to that in the ~1 Myr old Taurus and Ophiuchus regions. These observations reveal the rapid decline in the number of small dust grains in disks with time, and probably their concomitant growth beyond millimeter sizes. Moreover, if IC348 is to form planets in the same proportion as detected in the field, these faint millimeter detections may represent the best candidates in the cluster to study the progression from planetesimals to planets.

preprint2011arXiv

Protoplanetary Disks and Their Evolution

Flattened, rotating disks of cool dust and gas extending for tens to hundreds of AU are found around almost all low mass stars shortly after their birth. These disks generally persist for several Myr, during which time some material accretes onto the star, some is lost through outflows and photoevaporation, and some condenses into centimeter- and larger-sized bodies or planetesimals. Through observations mainly at infrared through millimeter wavelengths, we can determine how common disks are at different ages, measure basic properties including mass, size, structure, and composition, and follow their varied evolutionary pathways. In this way, we see the first steps toward exoplanet formation and learn about the origins of the Solar System. This review addresses observations of the outer parts, beyond 1 AU, of protoplanetary disks with a focus on recent infrared and (sub-)millimeter results and an eye to the promise of new facilities in the immediate future.

preprint2011arXiv

Resolved Submillimeter Observations of the HR 8799 and HD 107146 Debris Disks

We present 880 um Submillimeter Array observations of the debris disks around the young solar analogue HD 107146 and the multiple-planet host star HR 8799, at an angular resolution of 3" and 6", respectively. We spatially resolve the inner edge of the disk around HR 8799 for the first time. While the data are not sensitive enough (with rms noise of 1 mJy) to constrain the system geometry, we demonstrate that a model by Su et al. (2009) based on the spectral energy distribution (SED) with an inner radius of 150 AU predicts well the spatially resolved data. Furthermore, by modeling simultaneously the SED and visibilities, we demonstrate that the dust is distributed in a broad (of order 100 AU) annulus rather than a narrow ring. We also model the observed SED and visibilities for the HD 107146 debris disk and generate a model of the dust emission that extends in a broad band between 50 and 170 AU from the star. We perform an a posteriori comparison with existing 1.3 mm CARMA observations and demonstrate that a smooth, axisymmetric model reproduces well all of the available millimeter-wavelength data.

preprint2011arXiv

Submillimeter sources behind the massive lensing clusters A370 and A2390

We report 850 micron Submillimeter Array (SMA) observations of four gravitationally lensed submillimeter galaxies (SMGs), A370-2, A2390-1, A2390-3 and A2390-4, which were originally discovered with the Submillimeter Common-User Bolometer Array (SCUBA). Our SMA detection of A370-2 with a submillimeter flux of 7.95 +/- 0.60 mJy unambiguously identifies the counterparts to this source at optical and radio wavelengths. A2390-1 is an ultraluminous infrared galaxies with a submillimeter flux of 5.55 +/- 0.92 mJy and a redshift of 1.8 +/- 0.2 computed from submillimeter/radio flux ratio analysis. We resolve A2390-3 into two components, A2390-3a and A2390-3b, with fluxes of 3.15 +/- 0.63 mJy and 1.92 +/- 0.60 mJy, respectively. The structure of the system could be consistent with morphological distortion by gravitational lensing. The lack of counterparts in the optical and infrared indicates a heavily dust-enshrouded nature, and a non-detection in the radio implies that these two sources probably lie at z > 4.7, which would make them among the most distant SMGs known to date. Our non-detection of A2390-4 suggests either that there are multiple fainter submillimeter sources within the SCUBA beam or that the SCUBA detection may have been false. Our precise positions allow us to determine accurate amplifications and fluxes for all our detected sources. Our new results give a shallower power-law fit (-1.10) to the faint-end 850 micron cumulative number counts than previous work. We emphasize the need for high-resolution observations of single-dish detected SMGs in order to measure accurately the faint end of the 850 micron counts.

preprint2011arXiv

The Late Stages of Protoplanetary Disk Evolution: A Millimeter Survey of Upper Scorpius

We present deep 1.2 millimeter photometry of 37 stars in the young (5 Myr) Upper Scorpius OB association, sensitive to ~4 x 10^-3 Mjup of cool millimeter dust. Disks around four low- and solar-mass stars are detected, as well as one debris disk around an intermediate mass star, with dust masses ranging from 3.6 x 10^-3 -- 1.0 x 10^-1 Mjup. The source with the most massive disk exhibits a transition-disk spectral energy distribution. Combining our results with previous studies, we find the millimeter-detection fraction of Class II sources has significantly decreased from younger ages, and comparison with near-infrared and Halpha measurements indicates the present disks have undergone significant evolution in composition or structure at all radii. The disks of Upper Scorpius represent the tail-end of the depletion of primordial disks; while a few near-solar mass stars may still sustain giant planet formation, this process has finished around higher mass stars

preprint2011arXiv

The TW Hya Disk at 870 microns: Comparison of CO and Dust Radial Structures

We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

preprint2010arXiv

A Submillimeter Array Survey of Protoplanetary Disks in the Orion Nebula Cluster

We present the full results of our 3-year long Submillimeter Array survey of protoplanetary disks in the Orion Nebula Cluster. We imaged 23 fields at 880 microns and 2 fields at 1330 microns, covering an area of ~6.5 arcmin^2 and containing 67 disks. We detected 42 disks with fluxes between 6-135 mJy and at rms noise levels between 0.6 to 5.3 mJy/beam. Thermal dust emission above any free-free component was measured in 40 of the 42 detections, and the inferred disk masses range from 0.003-0.07 Msolar. We find that disks located within 0.3 pc of theta^1 Ori C have a truncated mass distribution, while disks located beyond 0.3 pc have masses more comparable to those found in low-mass star forming regions. The disk mass distribution in Orion has a distance dependence, with a derived relationship max(M_(disk)) = 0.046Msolar(d/0.3pc)^0.33 for the maximum disk masses. We found evidence of grain growth in disk 197-427, the only disk detected at both 880 microns and 1330 microns with the SMA. Despite the rapid erosion of the outer parts of the Orion disks by photoevaporation, the potential for planet formation remains high in this massive star forming region, with approximately 18% of the surveyed disks having masses greater than or equal to 0.01 Msolar within 60 AU.

preprint2010arXiv

A submillimetre survey of the kinematics of the Perseus molecular cloud - II. Molecular outflows

We present a census of molecular outflows across four active regions of star formation in the Perseus molecular cloud (NGC 1333, IC348/HH211, L1448 and L1455), totalling an area of over 1000 sq arcmin. This is one of the largest surveys of outflow evolution in a single molecular cloud published to date. We analyse large-scale, sensitive CO J=3-2 datasets from the James Clerk Maxwell Telescope, including new data towards NGC 1333. Where possible we make use of our complementary 13CO and C18O data to correct for the 12CO optical depth and measure ambient cloud properties. Of the 65 submillimetre cores in our fields, we detect outflows towards 45. We compare various parameters between the outflows from Class 0 and I protostars, including their mass, momentum, energy and momentum flux. Class 0 outflows are longer, faster, more massive and have more energy than Class I outflows. The dynamical time-scales we derive from these outflows are uncorrelated to the age of the outflow driving source, computed from the protostar's bolometric temperature. We confirm the results of Bontemps et al., that outflows decrease in force as they age, suggesting a decline in the mass accretion rate. If F_rad=L_bol/c is the flux expected in radiation from the central source, then F_CO(Class I)~100F_rad and F_CO(Class 0)~1000F_rad. Finally, we note that the total energy contained in outflows in NGC 1333, L1448 and L1455 is greater than the estimated turbulent energy in the respective regions, which may have implications for the regions' evolution.

preprint2010arXiv

The Astrophysical Environment of the Solar Birthplace

Our Sun, like all stars, formed within a cold molecular cloud. Astronomical observations and theory provide considerable detail into this process. Yet cosmochemical observations of short lived radionuclides in primitive meteorites, in particular 60Fe, provide unequivocal evidence that the early solar system inherited fresh nucleosynthetic material from the core of a hot, massive star, almost certainly ejected in a supernova explosion. I give a short introduction to the fields of star formation and meteoritics and discuss how the reconciliation of their disparate clues to our origin places strong constraints on the environment of the Solar birthplace. Direct injection of supernova ejecta into a protoplanetary disk or a dense molecular core is unlikely since their small sizes require placement unusually close to the massive star. Lower density molecular cloud clumps can capture more ejecta but the radionuclides decay during the slow gravitational collapse. The most likely scenario is on the largest scales via the formation of enriched molecular clouds at the intersection of colliding supernova bubbles in spiral arms.

preprint2010arXiv

The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

preprint2010arXiv

The Bolocam Galactic Plane Survey: Survey Description and Data Reduction

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

preprint2010arXiv

The Nature of Transition Circumstellar Disks I. The Ophiuchus Molecular Cloud

We have obtained millimeter wavelength photometry, high-resolution optical spectroscopy and adaptive optics near-infrared imaging for a sample of 26 Spitzer-selected transition circumstellar disks. All of our targets are located in the Ophiuchus molecular cloud (d ~125 pc) and have Spectral Energy Distributions (SEDs) suggesting the presence of inner opacity holes. We use these ground-based data to estimate the disk mass, multiplicity, and accretion rate for each object in our sample in order to investigate the mechanisms potentially responsible for their inner holes. We find that transition disks are a heterogeneous group of objects, with disk masses ranging from < 0.6 to 40 M_JUP and accretion rates ranging from <1E-11 1E-7 M_solar/yr, but most tend to have much lower masses and accretion rates than "full disks" (i.e., disks without opacity holes). Eight of our targets have stellar companions: 6 of them are binaries and the other 2 are triple systems. In four cases, the stellar companions are close enough to suspect they are responsible for the inferred inner holes. We find that 9 of our 26 targets have low disk mass (< 2.5 M_JUP) and negligible accretion (< 1E-11) M_solar/yr, and are thus consistent with photoevaporating (or photoevaporated) disks. Four of these 9 non-accreting objects have fractional disk luminosities < 1E-3 and could already be in a debris disk stage. Seventeen of our transition disks are accreting. Thirteen of these accreting objects are consistent with grain growth. The remaining 4 accreting objects have SEDs suggesting the presence of sharp inner holes, and thus are excellent candidates for harboring giant planets.

preprint2009arXiv

A Spatially Resolved Inner Hole in the Disk around GM Aurigae

We present 0.3 arcsec resolution observations of the disk around GM Aurigae with the Submillimeter Array (SMA) at a wavelength of 860 um and with the Plateau de Bure Interferometer at a wavelength of 1.3 mm. These observations probe the distribution of disk material on spatial scales commensurate with the size of the inner hole predicted by models of the spectral energy distribution. The data clearly indicate a sharp decrease in millimeter optical depth at the disk center, consistent with a deficit of material at distances less than ~20 AU from the star. We refine the accretion disk model of Calvet et al. (2005) based on the unresolved spectral energy distribution (SED) and demonstrate that it reproduces well the spatially resolved millimeter continuum data at both available wavelengths. We also present complementary SMA observations of CO J=3-2 and J=2-1 emission from the disk at 2" resolution. The observed CO morphology is consistent with the continuum model prediction, with two significant deviations: (1) the emission displays a larger CO J=3-2/J=2-1 line ratio than predicted, which may indicate additional heating of gas in the upper disk layers; and (2) the position angle of the kinematic rotation pattern differs by 11 +/- 2 degrees from that measured at smaller scales from the dust continuum, which may indicate the presence of a warp. We note that photoevaporation, grain growth, and binarity are unlikely mechanisms for inducing the observed sharp decrease in opacity or surface density at the disk center. The inner hole plausibly results from the dynamical influence of a planet on the disk material. Warping induced by a planet could also potentially explain the difference in position angle between the continuum and CO data sets.

preprint2009arXiv

Al-26 and the formation of the Solar System from a molecular cloud contaminated by Wolf-Rayet winds

In agreement with previous work, we show that the presence of the short-lived radionuclide Al-26 in the early Solar System was unlikely (<2% a priori probability) to be the result of direct introduction of supernova ejecta into the gaseous disk during the Class II stage of protosolar evolution. We also show that any Bondi-Hoyle accretion of contaminated residual gas from the natal star cluster made a negligible contribution to the primordial Al-26 inventory of the Solar System. These results are consistent with the absence of the oxygen isotopic signature expected with any late introduction of supernova ejecta into the protoplanetary disk. Instead, the presence of Al-26 in the oldest Solar System solids (calcium-aluminum-rich inclusions or CAIs) and its apparent uniform distribution with the inferred canonical Al-26/Al-27 ratio of (4.5-5)E-5 support the inheritance of Al-26 from the parent giant molecular cloud. We propose that this radionuclide originated in a prior generation of massive stars that formed in the same host molecular cloud as the Sun and contaminated that cloud by Wolf-Rayet winds. We calculated the Galactic distribution of Al-26/Al-27 ratios that arise from such contamination using the established embedded cluster mass and stellar initial mass functions, published nucleosynthetic yields from the winds of massive stars, and by assuming rapid and uniform mixing into the cloud. Although our model predicts that the majority of stellar systems contain no Al-26 from massive stars, and that the a priori probability that the Al-26/Al-27 ratio will reach or exceed the canonical Solar System value is only ~6%, the maximum in the distribution of non-zero values is close to the canonical ratio.

preprint2009arXiv

Massive Protoplanetary Disks in Orion Beyond the Trapezium Cluster

We present Submillimeter Array observations of the 880 micron continuum emission from three circumstellar disks around young stars in Orion that lie several arcminutes (> 1-pc) north of the Trapezium cluster. Two of the three disks are in the binary system 253-1536. Silhouette disks 216-0939 and 253-1536a are found to be more massive than any previously observed Orion disks, with dust masses derived from their submillimeter emission of 0.045 Msun and 0.066 Msun, respectively. The existence of these massive disks reveals the disk mass distribution in Orion does extend to high masses, and that the truncation observed in the central Trapezium cluster is a result of photoevaporation due to the proximity of O-stars. 253-1536b has a disk mass of 0.018 Msun, making the 253-1536 system the first optical binary in which each protoplanetary disk is massive enough to potentially form Solar systems.

preprint2009arXiv

Molecular Rings and the Thickness of Star-Forming Clouds

The winds and radiation from massive stars clear out large cavities in the interstellar medium. These bubbles, as they have been called, impact their surrounding molecular clouds and may influence the formation of stars therein. Here we present JCMT observations of the J=3-2 line of CO in 43 bubbles identified with Spitzer Space Telescope observations. These spectroscopic data reveal the three-dimensional structure of the bubbles. In particular, we show that the cold gas lies in a ring, not a sphere, around the bubbles indicating that the parent molecular clouds are flattened with a typical thickness of a few parsecs. We also mapped 7 bubbles in the J=4-3 line of HCO+ and find that the column densities inferred from the CO and HCO+ line intensities are below that necessary for "collect and collapse" models of induced star formation. We hypothesize that the flattened molecular clouds are not greatly compressed by expanding shock fronts, which may hinder the formation of new stars.

preprint2009arXiv

The Bolocam Galactic Plane Survey -- II. Catalog of the Image Data

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

preprint2006arXiv

Hot corinos in NGC1333-IRAS4B and IRAS2A

(Abridged) We aim to enlarge the number of known hot corinos and carry out a first comparative study with hot cores. The ultimate goal is to understand whether complex organic molecules form in the gas phase or on grain surfaces, and what the possible key parameters are. We observed millimeter rotational transitions of HCOOH, HCOOCH3, CH3OCH3, CH3CN, and C2H5CN in a sample of low-mass protostars with the IRAM-30m. Using the rotational diagram method coupled with the information about the sources' structure, we calculate the abundances of the observed molecules. To interpret these abundances, we review the proposed formation processes of the above molecules. We report the detection of HCOOCH3 and/or CH3CN towards NGC1333-IRAS4B and NGC1333-IRAS2A. We find that abundance ratios of O-bearing molecules to methanol or formaldehyde in hot corinos are comparable and about unity, and are relatively (depending on how the ratios are determined) higher than those in hot cores and in Galactic center clouds. So far, complex organic molecules were detected in all the hot corinos where they were searched for, suggesting that it is a common phase for low-mass protostars. While some evidence points to grain-surface synthesis (either in the cold or warm-up phase) of these molecules (in particular for HCOOH and HCOOCH3), the present data do not allow us to disregard gas-phase formation. More observational, laboratory, and theoretical studies are required to improve our understanding of hot corinos.

preprint2003arXiv

High resolution imaging of CO outflows in OMC-2 and OMC-3

A large scale, high resolution map of CO(1-0) emission toward the OMC-2 and OMC-3 star forming regions is presented. The map is a mosaic of 46 fields using the Berkeley-Illinois-Maryland Array (BIMA) and covers $\sim 10'\times 15'$ at $\sim 10" $ resolution. These data are combined with singledish FCRAO observations and analyzed to identify and determine the properties of nine protostellar outflows. The BIMA data alone almost completely resolve out the cloud emission at central velocities and only recover 1/20 of the flux in the high velocity gas showing that outflows are generally broadly dispersed over $\sim 1'$ angular scales. All nine identified outflows emanate from known Class 0 or borderline Class 0/I sources, are associated with knots of shocked \h2\ emission, and have short dynamical times. It is suggested that only the youngest, most spatially compact, and energetic outflows have been found and that more distributed high velocity gas undetected by BIMA is due to older outflows continuing through the Class I phase of protostellar evolution. The mechanical energy injection rate into the cloud is estimated to be $\sim 1.5 L_\odot$ which is comparable to the turbulent energy dissipation rate. Outf`lows appear capable, therefore, of sustaining cloud turbulence but a high starformation rate is required implying a short cloud lifetime $\leq 5$ `Myr.