BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sean M. Andrews

Sean M. Andrews appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAastro-ph.COastro-ph.IM

Catalog footprint

What is connected

82works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

82 published item(s)

preprint2026arXiv

Tracing Pebble Drift History in Two Protoplanetary Disks with CO Enhancement

Pebble drift is an important mechanism for supplying the materials needed to build planets in the inner region of protoplanetary disks. Thus, constraining pebble drift's timescales and mass flux is essential to understanding planet formation history. Current pebble drift models suggest pebble fluxes can be constrained from the enhancement of gaseous volatile abundances when icy pebbles sublimate after drifting across key snowlines. In this work, we present ALMA observations of spatially resolved $^{13}$C$^{18}$O J=2-1 line emission inside the midplane CO snowline of the HD 163296 and MWC 480 protoplanetary disks. We use radiative transfer and thermochemical models to constrain the spatial distribution of CO gas column density. We find that both disks display centrally peaked CO abundance enhancement of up to ten times of ISM abundance levels. For HD 163296 and MWC 480, the inferred enhancements require 250-350 and 480-660 Earth Masses of pebbles to have drifted across their CO snowlines, respectively. These ranges fall within cumulative pebble mass flux ranges to grow gas giants in the interior to the CO snowline. The centrally peaked CO enhancement is unexpected in current pebble drift models, which predict CO enhancement peaks at the CO snowline or is uniform inside the snowline. We propose two hypotheses to explain the centrally-peaked CO enhancement, including a large CO desorption distance and CO trapped in water ice. By testing both hypotheses with the 1D gas and dust evolution code chemcomp, we find that volatile trapping (about 30\%) best reproduces the centrally peaked CO enhancement observed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A large population study of protoplanetary disks: Explaining the millimeter size-luminosity relation with or without sub-structure

Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $α$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.5 arcsec (200 au) and a median rms noise of 54 mircoJy/beam, which corresponds to a fractional luminosity $L_{dust}/L_{star} \sim 10^{-4}$ for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations.The ALMA non-detections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (HII 1132 and HD 22680) show 24 micron excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gas and Dust Shadows in the TW Hydrae Disk

We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the dust shadow in scattered light from Debes et al. (2017) suggests that an oscillatory motion would provide a better model of the temporal evolution of the dust shadow rather than orbital motion. Both models predict an angular offset between the dust shadow and the gas shadow of up to ~100 deg. We attribute this offset to the finite rate at which dust grains and gas molecules can exchange heat, dominated by the collisional rate between gas molecules and dust grains, $t_{\rm coll}$. The angular offsets derived are equivalent to collisional timescales that range from the near instantaneous up to $t_{\rm coll}$ ~ 10 years, depending on whether a straight or a curved dust shadow, as suggested by HST observations reported by Debes et al. (2017), is adopted. The inferred range of $t_{\rm coll}$ are consistent with those predictions based on representative gas densities, temperatures, gas-to-dust ratios and grain sizes. These results represent the first time empirical constraints can be placed on $t_{\rm coll}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gas Disk Sizes from CO Line Observations: A Test of Angular Momentum Evolution

The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $α$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade

Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λ$ = 450 $μ$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μ$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scattering and sublimation: a multi-scale view of $μ$m-sized dust in the inclined disc of HD 145718

We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim λ/2π$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamical Masses and Stellar Evolutionary Model Predictions of M-Stars

In this era of Gaia and ALMA, dynamical stellar mass measurements provide benchmarks that are independent of observations of stellar characteristics and their uncertainties. These benchmarks can then be used to validate and improve stellar evolutionary models, which can lead to both imprecise and inaccurate mass predictions for pre-main-sequence, low-mass stars. We present the dynamical stellar masses derived from disks around three M-stars (FP Tau, J0432+1827, and J1100-7619) using ALMA observations of $^{12}$CO (J=2--1) and $^{13}$CO (J=2--1) emission. These are the first dynamical stellar mass measurements for J0432+1827 and J1100-7619 and the most precise measurement for FP Tau. Fiducial stellar evolutionary model tracks, which do not include any treatment of magnetic activity, agree with the dynamical measurement of J0432+1827 but underpredict the mass by $\sim$60\% for FP Tau and $\sim$80\% for J1100-7619. Possible explanations for the underpredictions include inaccurate assumptions of stellar effective temperature, undetected binarity for J1100-7619, and that fiducial stellar evolutionary models are not complex enough to represent these stars. In the former case, the stellar effective temperatures would need to be increased by $\sim$40K to $\sim$340K to reconcile the fiducial model predictions with the dynamically-measured masses. In the latter case, we show that the dynamical masses can be reproduced using results from stellar evolutionary models with starspots, which incorporate fractional starspot coverage to represent the manifestation of magnetic activity. Folding in low-mass M-stars from the literature and assuming that the stellar effective temperatures are imprecise but accurate, we find tentative evidence of a relationship between fractional starspot coverage and observed effective temperature for these young, cool stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $α$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ALMA 0.88 mm Survey of Disks around Planetary-mass Companions

Characterizing the physical properties and compositions of circumplanetary disks can provide important insights into the formation of giant planets and satellites. We report ALMA 0.88 mm (Band 7) continuum observations of six planetary-mass (10-20 $M_{\rm Jup}$) companions: CT Cha b, 1RXS 1609 b, ROXs 12 b, ROXs 42B b, DH Tau b, and FU Tau b. No continuum sources are detected at the locations of the companions down to 3$σ$ limits of 120-210 $μ$Jy. Given these non-detections, it is not clear whether disks around planetary-mass companions indeed follow the disk flux-host mass trend in the stellar regime. The faint radio brightness of these companion disks may result from a combination of fast radial drift and a lack of dust traps. Alternatively, as disks in binary systems are known to have significantly lower millimeter fluxes due to tidal interactions, these companion disks may instead follow the relationship of moderate-separation binary stars. This scenario can be tested with sensitive continuum imaging at rms levels of $\lesssim$10 $μ$Jy.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($α_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($τ\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring turbulent motion in planet-forming disks with ALMA: A detection around DM Tau and non-detections around MWC 480 and V4046 Sgr

Turbulence is a crucial factor in many models of planet formation, but it has only been directly constrained among a small number of planet forming disks. Building on the upper limits on turbulence placed in disks around HD 163296 and TW Hya, we present ALMA CO J=2-1 line observations at $\sim0.3"$ (20-50 au) resolution and 80 ms$^{-1}$ channel spacing of the disks around DM Tau, MWC 480, and V4046 Sgr. Using parametric models of disk structure, we robustly detect non-thermal gas motions around DM Tau of between 0.25 c$_s$ and 0.33 c$_s$, with the range dominated by systematic effects, making this one of the only systems with directly measured non-zero turbulence. Using the same methodology, we place stringent upper limits on the non-thermal gas motion around MWC 480 ($<$0.08c$_s$) and V4046 Sgr ($<$0.12 c$_s$). The preponderance of upper limits in this small sample, and the modest turbulence levels consistent with dust studies, suggest that weak turbulence ($α\lesssim10^{-3}$) may be a common, albeit not universal, feature of planet-forming disks. We explore the particular physical conditions around DM Tau that could lead this system to be more turbulent than the others.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on Planetesimal Collision Models in Debris Disks

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and $\sim5$ $μ$Jy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, $n(a)\propto a^{-q}$, of $\langle q \rangle = 3.36\pm0.02$, with a possible trend of decreasing $q$ for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution ($q=3.5$) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence for a CO desorption front in the outer AS 209 disk

Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=2-1$ emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the C$^{18}$O emission consists of a central peak and a ring at a radius of $\sim1''$ (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resolved Gas Interior to the Dust Rings of the HD 141569 Disk

The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk

We present long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870 micron continuum emission from the nearest gas-rich protoplanetary disk, around TW Hya, that trace millimeter-sized particles down to spatial scales as small as 1 AU (20 mas). These data reveal a series of concentric ring-shaped substructures in the form of bright zones and narrow dark annuli (1-6 AU) with modest contrasts (5-30%). We associate these features with concentrations of solids that have had their inward radial drift slowed or stopped, presumably at local gas pressure maxima. No significant non-axisymmetric structures are detected. Some of the observed features occur near temperatures that may be associated with the condensation fronts of major volatile species, but the relatively small brightness contrasts may also be a consequence of magnetized disk evolution (the so-called zonal flows). Other features, particularly a narrow dark annulus located only 1 AU from the star, could indicate interactions between the disk and young planets. These data signal that ordered substructures on ~AU scales can be common, fundamental factors in disk evolution, and that high resolution microwave imaging can help characterize them during the epoch of planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spiral Density Waves in a Young Protoplanetary Disk

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Disk-based Dynamical Mass Estimate for the Young Binary AK Sco

We present spatially and spectrally resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of gas and dust in the disk orbiting the pre-main sequence binary AK Sco. By forward-modeling the disk velocity field traced by CO J=2-1 line emission, we infer the mass of the central binary, $M_\ast = 2.49 \pm 0.10~M_\odot$, a new dynamical measurement that is independent of stellar evolutionary models. Assuming the disk and binary are co-planar within $\sim$2°, this disk-based binary mass measurement is in excellent agreement with constraints from radial velocity monitoring of the combined stellar spectra. These ALMA results are also compared with the standard approach of estimating masses from the location of the binary in the Hertzsprung-Russell diagram, using several common pre-main sequence model grids. These models predict stellar masses that are marginally consistent with our dynamical measurement (at $\sim 2\,σ$), but are systematically high (by $\sim$10%). These same models consistently predict an age of $18\pm1$ Myr for AK Sco, in line with its membership in the Upper Centaurus-Lupus association but surprisingly old for it to still host a gas-rich disk. As ALMA accumulates comparable data for large samples of pre-main sequence stars, the methodology employed here to extract a dynamical mass from the disk rotation curve should prove extraordinarily useful for efforts to characterize the fundamental parameters of early stellar evolution.

0 citations0 reviewsNo code linkedOpen access
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 .

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An ALMA Constraint on the GSC 6214-210 B Circum-Substellar Accretion Disk Mass

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of GSC 6214-210 A and B, a solar-mass member of the 5-10 Myr Upper Scorpius association with a 15 $\pm$ 2 Mjup companion orbiting at $\approx$330 AU (2.2"). Previous photometry and spectroscopy spanning 0.3-5 $μ$m revealed optical and thermal excess as well as strong H$α$ and Pa~$β$ emission originating from a circum-substellar accretion disk around GSC 6214-210 B, making it the lowest mass companion with unambiguous evidence of a subdisk. Despite ALMA's unprecedented sensitivity and angular resolution, neither component was detected in our 880 $μ$m (341 GHz) continuum observations down to a 3-$σ$ limit of 0.22 mJy/beam. The corresponding constraints on the dust mass and total mass are <0.15 Mearth and <0.05 Mjup, respectively, or <0.003% and <0.3% of the mass of GSC 6214-210 B itself assuming a 100:1 gas-to-dust ratio and characteristic dust temperature of 10-20 K. If the host star possesses a putative circum-stellar disk then at most it is a meager 0.0015% of the primary mass, implying that giant planet formation has certainly ceased in this system. Considering these limits and its current accretion rate, GSC 6214-210 B appears to be at the end stages of assembly and is not expected to gain any appreciable mass over the next few Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chemical Imaging of the CO Snow Line in the HD 163296 Disk

The condensation fronts (snow lines) of H2O, CO and other abundant volatiles in the midplane of a protoplanetary disk affect several aspects of planet formation. Locating the CO snow line, where the CO gas column density is expected to drop substantially, based solely on CO emission profiles is challenging. This has prompted an exploration of chemical signatures of CO freeze-out. We present ALMA Cycle 1 observations of the N2H+ J=3-2 and DCO+ J=4-3 emission lines toward the disk around the Herbig Ae star HD~163296 at ~0.5" (60 AU) resolution, and evaluate their utility as tracers of the CO snow line location. The N2H+ emission is distributed in a ring with an inner radius at 90 AU, corresponding to a midplane temperature of 25 K. This result is consistent with a new analysis of optically thin C18O data, which implies a sharp drop in CO abundance at 90 AU. Thus N2H+ appears to be a robust tracer of the midplane CO snow line. The DCO+ emission also has a ring morphology, but neither the inner nor the outer radius coincides with the CO snow line location of 90 AU, indicative of a complex relationship between DCO+ emission and CO freeze-out in the disk midplane. Compared to TW Hya, CO freezes out at a higher temperature in the disk around HD 163296 (25 vs. 17 K in the TW Hya disk), perhaps due to different ice compositions. This highlights the importance of actually measuring the CO snow line location, rather than assuming a constant CO freeze-out temperature for all disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constructing A Flexible Likelihood Function For Spectroscopic Inference

We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically address the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line "outliers." By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high resolution $V$-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate resolution $K$-band spectrum of Gliese 51, an M5 field dwarf.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Demographics of Transition Discs in Ophiuchus and Taurus

Transition disc systems are young stars that appear to be on the verge of dispersing their protoplanetary discs. We explore the nature of these systems by comparing the stellar accretion rates and disc masses of transition discs and normal T Tauri stars in Taurus and Ophiuchus. After controlling for the known dependencies of stellar accretion rate and disc mass and on age, stellar accretion rate on stellar mass, and disc mass on the presence of stellar or sub-stellar companions, we find that the normal T Tauri stars show a trend of stellar accretion rate increasing with disc mass. The transition discs tend to have higher average disc masses than normal T Tauri stars as well as lower accretion rates than normal T Tauri stars of the same disc mass. These results are most consistent with the interpretation that the transition discs have formed objects massive enough to alter the accretion flow, i.e., single or multiple giant planets. Several Ophiuchus T Tauri stars that are not known transition disc systems also have very low accretion rates for their disc masses. We speculate on the possible nature of these sources.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Double DCO+ rings reveal CO ice desorption in the outer disk around IM Lup

In a protoplanetary disk, a combination of thermal and non-thermal desorption processes regulate where volatiles are liberated from icy grain mantles into the gas phase. Non-thermal desorption should result in volatile-enriched gas in disk-regions where complete freeze-out is otherwise expected. We present ALMA observations of the disk around the young star IM Lup in 1.4 mm continuum, C18O 2-1, H13CO+ 3-2 and DCO+ 3-2 emission at ~0".5 resolution. The images of these dust and gas tracers are clearly resolved. The DCO+ line exhibits a striking pair of concentric rings of emission that peak at radii of ~0".6 and 2" (~90 and 300 AU, respectively). Based on disk chemistry model comparison, the inner DCO+ ring is associated with the balance of CO freeze-out and thermal desorption due to a radial decrease in disk temperature. The outer DCO+ ring is explained by non-thermal desorption of CO ice in the low-column-density outer disk, repopulating the disk midplane with cold CO gas. The CO gas then reacts with abundant H2D+ to form the observed DCO+ outer ring. These observations demonstrate that spatially resolved DCO+ emission can be used to trace otherwise hidden cold gas reservoirs in the outmost disk regions, opening a new window onto their chemistry and kinematics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dust Evolution Can Produce Scattered Light Gaps in Protoplanetary Disks

Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps." The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at λ = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at λ = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, β, is consistent with the observed dust emission in both disks, with low-β in the inner disk and high-β in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observations of Solids in Protoplanetary Disks

This review addresses the state of research that employs astronomical (remote sensing) observations of solids ("dust") in young circumstellar disks to learn about planet formation. The intention is for it to serve as an accessible, introductory, pedagogical resource for junior scientists interested in the subject. After some historical background and a basic observational primer, the focus is shifted to the three fundamental topics that broadly define the field: (1) demographics -- the relationships between disk properties and the characteristics of their environments and hosts; (2) structure -- the spatial distribution of disk material and its associated physical conditions and composition; and (3) evolution -- the signposts of key changes in disk properties, including the growth and migration of solids and the impact of dynamical interactions with young planetary systems. Based on the state of the art results in these areas, suggestions are made for potentially fruitful lines of work in the near future.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved gas cavities in transitional disks inferred from CO isotopologues with ALMA

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved Millimeter Emission from the HD 15115 Debris Disk

We have used the Submillimeter Array (SMA) to make 1.3 millimeter observations of the debris disk surrounding HD 15115, an F-type star with a putative membership in the beta Pictoris moving group. This nearly edge-on debris disk shows an extreme asymmetry in optical scattered light, with an extent almost two times larger to the west of the star than to the east (originally dubbed the "Blue Needle"). The SMA observations reveal resolved emission that we model as a circumstellar belt of thermal dust emission. This belt extends to a radius of ~110 AU, coincident with the break in the scattered light profile convincingly seen on the western side of the disk. This outer edge location is consistent with the presence of an underlying population of dust-producing planetesimals undergoing a collisional cascade, as hypothesized in "birth ring" theory. In addition, the millimeter emission shows a ~3 sigma feature aligned with the asymmetric western extension of the scattered light disk. If this millimeter extension is real, then mechanisms for asymmetry that affect only small grains, such as interactions with interstellar gas, are disfavored. This tentative feature might be explained by secular perturbations to grain orbits introduced by neutral gas drag, as previously invoked to explain asymmetric morphologies of other, similar debris disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved Millimeter-Wavelength Observations of Debris Disks around Solar-Type Stars

The presence of debris disks around young main sequence stars hints at the existence and structure of planetary systems. Millimeter-wavelength observations probe large grains that trace the location of planetesimal belts. The FEPS (Formation and Evolution of Planetary Systems) $Spitzer$ Legacy survey of nearby young solar analogues yielded a sample of five debris disk-hosting stars with millimeter flux suitable for interferometric follow-up. We present observations with the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at ~2" resolution that spatially resolve the debris disks around these nearby ($d\sim$50 pc) stars. Two of the five disks (HD 377, HD 8907) are spatially resolved for the first time and one (HD 104860) is resolved at millimeter wavelengths for the first time. We combine our new observations with archival SMA and Atacama Large Millimeter/Submillimeter Array (ALMA) data to enable a uniform analysis of the full five-object sample. We simultaneously model the broad-band photometric data and resolved millimeter visibilities to constrain the dust temperatures and disk morphologies, and perform an MCMC analysis to fit for basic structural parameters. We find that the radii and widths of the cold outer belts exhibit properties consistent with scaled-up versions of the Solar System's Kuiper Belt. All the disks exhibit characteristic grain sizes comparable to the blowout size, and all the resolved observations of emission from large dust grains are consistent with an axisymmetric dust distribution to within the uncertainties. These results are consistent with comparable studies carried out at infrared wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Scattered Light from Dust in the Cavity of the V4046 Sgr Transition Disk

We report the presence of scattered light from dust grains located in the giant planet formation region of the circumbinary disk orbiting the ~20-Myr-old close (~0.045 AU separation) binary system V4046 Sgr AB based on observations with the new Gemini Planet Imager (GPI) instrument. These GPI images probe to within ~7 AU of the central binary with linear spatial resolution of ~3 AU, and are thereby capable of revealing dust disk structure within a region corresponding to the giant planets in our solar system. The GPI imaging reveals a relatively narrow (FWHM ~10 AU) ring of polarized near-infrared flux whose brightness peaks at ~14 AU. This ~14 AU radius ring is surrounded by a fainter outer halo of scattered light extending to ~45 AU, which coincides with previously detected mm-wave thermal dust emission. The presence of small grains that efficiently scatter starlight well inside the mm-wavelength disk cavity supports current models of planet formation that suggest planet-disk interactions can generate pressure traps that impose strong radial variations in the particle size distribution throughout the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The cometary composition of a protoplanetary disk as revealed by complex cyanides

Observations of comets and asteroids show that the Solar Nebula that spawned our planetary system was rich in water and organic molecules. Bombardment brought these organics to the young Earth's surface, seeding its early chemistry. Unlike asteroids, comets preserve a nearly pristine record of the Solar Nebula composition. The presence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect to water, is of special interest because of the importance of C-N bonds for abiotic amino acid synthesis. Comet-like compositions of simple and complex volatiles are found in protostars, and can be readily explained by a combination of gas-phase chemistry to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances complexity[3]. Simple volatiles, including water and HCN, have been detected previously in Solar Nebula analogues - protoplanetary disks around young stars - indicating that they survive disk formation or are reformed in situ. It has been hitherto unclear whether the same holds for more complex organic molecules outside of the Solar Nebula, since recent observations show a dramatic change in the chemistry at the boundary between nascent envelopes and young disks due to accretion shocks[8]. Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary disk around the young star MWC 480. We find abundance ratios of these N-bearing organics in the gas-phase similar to comets, which suggests an even higher relative abundance of complex cyanides in the disk ice. This implies that complex organics accompany simpler volatiles in protoplanetary disks, and that the rich organic chemistry of the Solar Nebula was not unique.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Distribution and Chemistry of H$_2$CO in the DM Tau Protoplanetary Disk

H$_2$CO ice on dust grains is an important precursor of complex organic molecules (COMs). H$_2$CO gas can be readily observed in protoplanetary disks and may be used to trace COM chemistry. However, its utility as a COM probe is currently limited by a lack of constraints on the relative contributions of two different formation pathways: on icy grain-surfaces and in the gas-phase. We use archival ALMA observations of the resolved distribution of H$_2$CO emission in the disk around the young low-mass star DM Tau to assess the relative importance of these formation routes. The observed H$_2$CO emission has a centrally peaked and radially broad brightness profile (extending out to 500 AU). We compare these observations with disk chemistry models with and without grain-surface formation reactions, and find that both gas and grain-surface chemistry are necessary to explain the spatial distribution of the emission. Gas-phase H$_2$CO production is responsible for the observed central peak, while grain-surface chemistry is required to reproduce the emission exterior to the CO snowline (where H$_2$CO mainly forms through the hydrogenation of CO ice before being non-thermally desorbed). These observations demonstrate that both gas and grain-surface pathways contribute to the observed H$_2$CO in disks, and that their relative contributions depend strongly on distance from the host star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Epsilon Eridani System Resolved by Millimeter Interferometry

We present observations of Epsilon Eridani from the Submillimeter Array (SMA) at 1.3 millimeters and from the Australia Telescope Compact Array (ATCA) at 7 millimeters that reach an angular resolution of ~4" (13 AU). These first millimeter interferometer observations of Epsilon Eridani, which hosts the closest debris disk to the Sun, reveal two distinct emission components: (1) the well-known outer dust belt, which, although patchy, is clearly resolved in the radial direction, and (2) an unresolved source coincident with the position of the star. We use direct model-fitting of the millimeter visibilities to constrain the basic properties of these two components. A simple Gaussian shape for the outer belt fit to the SMA data results in a radial location of $64.4^{+2.4}_{-3.0}$ AU and FWHM of $20.2^{+6.0}_{-8.2}$ AU (fractional width $ΔR/R = 0.3$. Similar results are obtained taking a power law radial emission profile for the belt, though the power law index cannot be usefully constrained. Within the noise obtained (0.2 mJy/beam), these data are consistent with an axisymmetric belt model and show no significant azimuthal structure that might be introduced by unseen planets in the system. These data also limit any stellocentric offset of the belt to $<9$ AU, which disfavors the presence of giant planets on highly eccentric ($>0.1$) and wide (10's of AU) orbits. The flux density of the unresolved central component exceeds predictions for the stellar photosphere at these long wavelengths, by a marginally significant amount at 1.3 millimeters but by a factor of a few at 7 millimeters (with brightness temperature $13000 \pm 1600$ K for a source size of the optical stellar radius). We attribute this excess emission to ionized plasma from a stellar corona or chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weak Turbulence in the HD 163296 Protoplanetary Disk Revealed by ALMA CO Observations

Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An ALMA Disk Mass for the Candidate Protoplanetary Companion to FW Tau

We present ALMA observations of the FW Tau system, a close binary pair of M5 stars with a wide-orbit (300 AU projected separation) substellar companion. The companion is extremely faint and red in the optical and near-infrared, but boasts a weak far-infrared excess and optical/near-infrared emission lines indicative of a primordial accretion disk of gas and dust. The component-resolved 1.3 mm continuum emission is found to be associated only with the companion, with a flux (1.78 +/- 0.03 mJy) that indicates a dust mass of 1-2 M_Earth. While this mass reservoir is insufficient to form a giant planet, it is more than sufficient to produce an analog of the Kepler-42 exoplanetary system or the Galilean satellites. The mass and geometry of the disk-bearing FW Tau companion remains unclear. Near-infrared spectroscopy shows deep water bands that indicate a spectral type later than M5, but substantial veiling prevents a more accurate determination of the effective temperature (and hence mass). Both a disk-bearing "planetary-mass" companion seen in direct light or a brown dwarf tertiary viewed in light scattered by an edge-on disk or envelope remain possibilities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Resolved Multifrequency Radio Observations of GG Tau

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Spatially Resolved Vertical Temperature Gradient in the HD 163296 Disk

We analyze sensitive, sub-arcsecond resolution ALMA Science Verification observations of CO emission lines in the protoplanetary disk hosted by the young, isolated Ae star HD 163296. The observed spatial morphology of the 12CO J=3-2 emission line is asymmetric across the major axis of the disk; the 12CO J=2-1 line features a much less pronounced, but similar, asymmetry. The J=2-1 emission from 12CO and its main isotopologues have no resolved spatial asymmetry. We associate this behavior as the direct signature of a vertical temperature gradient and layered molecular structure in the disk. This is demonstrated using both toy models and more sophisticated calculations assuming non-local thermodynamic equilibrium (non LTE) conditions. A model disk structure is developed to reproduce both the distinctive spatial morphology of the 12CO J=3-2 line as well as the J=2-1 emission from the CO isotopologues assuming relative abundances consistent with the interstellar medium. This model disk structure has tau=1 emitting surfaces for the 12CO emission lines that make an angle of about 15 degrees with respect to the disk midplane. Furthermore, we show that the spatial and spectral sensitivity of these data can distinguish between models that have sub-Keplerian gas velocities due to the vertical extent of the disk and its associated radial pressure gradient (a fractional difference in the bulk gas velocity field of approximately greater than 5%).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dust Properties of Local Dust-Obscured Galaxies with the Submillimeter Array

We report Submillimeter Array (SMA) observations of the 880 μm dust continuum emission for four dust-obscured galaxies (DOGs) in the local universe. Two DOGs are clearly detected with S_ν(880 μm)=10-13 mJy and S/N>5, but the other two are not detected with 3σ upper limits of S_ν(880 μm)=5-9 mJy. Including an additional two local DOGs with submillimeter data from the literature, we determine the dust masses and temperatures for six local DOGs. The infrared luminosities and dust masses for these DOGs are in the range 1.2-4.9x10^{11} (L_sun) and 4-14x10^{7} (M_sun), respectively. The dust temperatures derived from a two-component modified blackbody function are 23-26 K and 60-124 K for the cold and warm dust components, respectively. Comparison of local DOGs with other infrared luminous galaxies with submillimeter detections shows that the dust temperatures and masses do not differ significantly among these objects. Thus, as argued previously, local DOGs are not a distinctive population among dusty galaxies, but simply represent the high-end tail of the dust obscuration distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast Radial Flows in Transition Disk Holes

Protoplanetary "transition" disks have large, mass-depleted central cavities, yet also deliver gas onto their host stars at rates comparable to disks without holes. The paradox of simultaneous transparency and accretion can be explained if gas flows inward at much higher radial speeds inside the cavity than outside the cavity, since surface density (and by extension optical depth) varies inversely with inflow velocity at fixed accretion rate. Radial speeds within the cavity might even have to approach free-fall values to explain the huge surface density contrasts inferred for transition disks. We identify observational diagnostics of fast radial inflow in channel maps made in optically thick spectral lines. Signatures include (1) twisted isophotes in maps made at low systemic velocities and (2) rotation of structures observed between maps made in high-velocity line wings. As a test case, we apply our new diagnostic tools to archival ALMA data on the transition disk HD 142527, and uncover evidence for free-fall radial velocities inside its cavity. Although the observed kinematics are also consistent with a disk warp, the radial inflow scenario is preferred because it predicts low surface densities that appear consistent with recent observations of optically thin CO isotopologues in this disk. How material in the disk cavity sheds its angular momentum wholesale to fall freely onto the star is an unsolved problem; gravitational torques exerted by giant planets or brown dwarfs are briefly discussed as a candidate mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imaging of the CO Snow Line in a Solar Nebula Analog

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mapping H-band Scattered Light Emission in the Mysterious SR21 Transitional Disk

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Outer Edges of Protoplanetary Dust Disks

The expectation that aerodynamic drag will force the solids in a gas-rich protoplanetary disk to spiral in toward the host star on short timescales is one of the fundamental problems in planet formation theory. The nominal efficiency of this radial drift process is in conflict with observations, suggesting that an empirical calibration of solid transport mechanisms in a disk is highly desirable. However, the fact that both radial drift and grain growth produce a similar particle size segregation in a disk (such that larger particles are preferentially concentrated closer to the star) makes it difficult to disentangle a clear signature of drift alone. We highlight a new approach, by showing that radial drift leaves a distinctive "fingerprint" in the dust surface density profile that is directly accessible to current observational facilities. Using an analytical framework for dust evolution, we demonstrate that the combined effects of drift and (viscous) gas drag naturally produce a sharp outer edge in the dust distribution (or, equivalently, a sharp decrease in the dust-to-gas mass ratio). This edge feature forms during the earliest phase in the evolution of disk solids, before grain growth in the outer disk has made much progress, and is preserved over longer timescales when both growth and transport effects are more substantial. The key features of these analytical models are reproduced in detailed numerical simulations, and are qualitatively consistent with recent millimeter-wave observations that find gas/dust size discrepancies and steep declines in dust continuum emission in the outer regions of protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Mass Dependence Between Protoplanetary Disks and their Stellar Hosts

We present a substantial extension of the mm-wave continuum photometry catalog for Taurus circumstellar dust disks. Combining new Submillimeter Array data with measurements in the literature, we construct a mm-wave luminosity distribution for Class II disks that is statistically complete for stellar hosts with spectral types earlier than M8.5 and has a (3-sigma) depth of ~3 mJy. The resulting census eliminates a longstanding bias against disks with late-type hosts, and thereby reveals a strong correlation between L_mm and the host spectral type. We confirm that this corresponds to a statistically robust relationship between the masses of dust disks and the stars that host them. A Bayesian regression technique is used to characterize these relationships: the results indicate a typical 1.3 mm flux density of 25 mJy for solar mass hosts and a power-law scaling L_mm \propto M_star^1.5-2.0. We suggest that a reasonable treatment of dust temperature in the conversion from L_mm to M_disk favors an inherently linear M_disk \propto M_star scaling, with a typical disk-to-star mass ratio of $\sim$0.2--0.6%. The RMS dispersion around this regression is 0.7 dex, suggesting that the combined effects of diverse evolutionary states, dust opacities, and temperatures in these disks imprint a FWHM range of a factor of 40 on the inferred M_disk (or L_mm) at any given host mass. We argue that this relationship between M_disk and M_star likely represents the origin of the inferred correlation between giant planet frequency and host star mass in the exoplanet population, and provides some basic support for the core accretion model for planet formation. Moreover, we caution that selection bias must be considered in comparative studies of disk evolution, and illustrate that fact with statistical comparisons of L_mm between Taurus and other clusters (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Structure of the Evolved Circumbinary Disk around V4046 Sgr

We present sensitive, sub-arcsecond resolution Submillimeter Array observations of the protoplanetary disk around the nearby, pre-main sequence spectroscopic binary V4046 Sgr. We report for the first time a large inner hole (r=29 AU) spatially resolved in the 1.3 mm continuum emission and study the structure of this disk using radiative transfer calculations to model the spectral energy distribution (SED), continuum visibilities, and spectral line emission of CO and its main isotopologues. Our modeling scheme demonstrates that the majority of the dust mass is distributed in a narrow ring (centered at 37 AU with a FWHM of 16 AU) that is ~5 times more compact than the gas disk. This structure implies that the dust-to-gas mass ratio has a strong spatial variation, ranging from a value much larger than typical of the interstellar medium (ISM) at the ring to much smaller than that of the ISM at larger disk radii. We suggest that these basic structural features are potentially observational signatures of the accumulation of solids at a local gas pressure maximum. These models also require a substantial population of ~micron-sized grains inside the central disk cavity. We suggest that this structure is likely the result of dynamical interactions with a low-mass companion, although photoevaporation may also play a secondary role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Disk-based Dynamical Mass Estimate for the Young Binary V4046 Sgr

We present sensitive, arcsecond-resolution Submillimeter Array observations of the 12CO J=2-1 line emission from the circumstellar disk orbiting the double-lined spectroscopic binary star V4046 Sgr. Based on a simple model of the disk structure, we use a novel Monte Carlo Markov Chain technique to extract the Keplerian velocity field of the disk from these data and estimate the total mass of the central binary. Assuming the distance inferred from kinematic parallax measurements in the literature (d is approximately 73 pc), we determine a total stellar mass M_star = 1.75^{+0.09}_{-0.06} solar masses and a disk inclination i_d = 33.5^{+0.7}_{-1.4} degrees from face-on. These measurements are in excellent agreement with independent dynamical constraints made from multi-epoch monitoring of the stellar radial velocities, confirming the absolute accuracy of this precise (~ few percent uncertainties) disk-based method for estimating stellar masses and reaffirming previous assertions that the disk and binary orbital planes are well aligned (with |i_d - i_star| \approx 0.1\pm1 degree). Using these results as a reference, we demonstrate that various pre-main sequence evolution models make consistent and accurate predictions for the masses of the individual components of the binary, and uniformly imply an advanced age of ~5-30 Myr. Taken together, these results verify that V4046 Sgr is one of the precious few nearby and relatively evolved pre-main sequence systems that still hosts a gas-rich accretion disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Resolved Census of Millimeter Emission from Taurus Multiple Star Systems

We present a high angular resolution millimeter-wave dust continuum imaging survey of circumstellar material associated with individual components of 23 multiple star systems in the Taurus-Auriga young cluster. Combined with previous measurements, these new data permit a comprehensive look at how millimeter luminosity (a tracer of disk mass) relates to the separation and mass of a stellar companion. Approximately one third (28-37%) of individual stars in multiples have detectable millimeter emission, a rate half that for single stars (~62%). There is a strong correlation between the luminosity and projected separation (a_p) of a stellar pair. Wide pairs (a_p > 300 AU) have a similar luminosity distribution as single stars, medium pairs (a_p ~ 30-300 AU) are a factor of 5 fainter, and close pairs (a_p < 30 AU) are ~ 5 times fainter yet (aside from a small population of bright circumbinary disks). In most cases, the emission is dominated by a disk around the primary (or a wide tertiary in triples), but there is no clear relationship between luminosity and stellar mass ratio. A direct comparison of resolved disk sizes with predictions from tidal truncation models yields mixed results; some disks are larger than expected given their companion separations. We suggest that the presence of a stellar companion impacts disk properties at a level comparable to the internal evolution mechanisms operating in isolated systems, with both the multiple star formation process itself and star-disk tidal interactions likely playing important roles in the evolution of disk material. From the perspective of the mass content of the disk, we expect that (giant) planet formation is inhibited around the components of close pairs or secondaries, but should be as likely as for single stars around the primaries (or wide tertiaries in hierarchical triples) in more widely-separated multiple star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Resolved Millimeter Emission Belt in the AU Mic Debris Disk

We present imaging observations at 1.3 millimeters of the debris disk surrounding the nearby M-type flare star AU Mic with beam size 3 arcsec (30 AU) from the Submillimeter Array. These data reveal a belt of thermal dust emission surrounding the star with the same edge-on geometry as the more extended scattered light disk detected at optical wavelengths. Simple modeling indicates a central radius of ~35 AU for the emission belt. This location is consistent with the reservoir of planetesimals previously invoked to explain the shape of the scattered light surface brightness profile through size-dependent dust dynamics. The identification of this belt further strengthens the kinship between the debris disks around AU Mic and its more massive sister star beta Pic, members of the same ~10 Myr-old moving group.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on the Radial Variation of Grain Growth in the AS 209 Circumstellar Disk

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity κ_ν. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (β), we find that β(R) increases from β<0.5 at \sim20 AU to β>1.5 for R>80 AU, inconsistent with a constant value of β across the disk (at the 10σ level). Furthermore, if radial variations of κ_ν are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinematics of the CO Gas in the Inner Regions of the TW Hya Disk

We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Matryoshka Holes: Nested Emission Rings in the Transitional Disk Oph IRS 48

The processes that form transition disks - disks with depleted inner regions - are not well understood; possible scenarios include planet formation, grain growth and photoevaporation. Disks with spatially resolved dust holes are rare, but, in general, even less is known about the gas structure. The disk surrounding A0 star Oph IRS 48 in the nearby Rho Ophiuchus region has a 30 AU radius hole previously detected in the 18.7 micron dust continuum and in warm CO in the 5 micron fundamental ro-vibrational band. We present here Submillimeter Array 880 micron continuum imaging resolving an inner hole. However, the radius of the hole in the millimeter dust is only 13 AU, significantly smaller than measured at other wavelengths. The nesting structure of the disk is counter-intuitive, with increasingly large radii rings of emission seen in the millimeter dust (12.9 +1.7/-3.4 AU), 5 micron CO (30 AU) and 18.7 micron dust (peaking at 55 AU). We discuss possible explanations for this structure, including self-shadowing that cools the disk surface layers, photodissociation of CO, and photoevaporation. However, understanding this unusual disk within the stringent multi-wavelength spatial constraints will require further observations to search for cold atomic and molecular gas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Millimeter Emission Structure in the first ALMA Image of the AU Mic Debris Disk

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Closer Look at the LkCa 15 Protoplanetary Disk

We present 870 micron observations of dust continuum emission from the LkCa 15 protoplanetary disk at high angular resolution (with a characteristic scale of 0.25" = 35 AU), obtained with the IRAM Plateau de Bure interferometer and supplemented by slightly lower resolution observations from the Submillimeter Array. We fit these data with simple morphological models to characterize the spectacular ring-like emission structure of this disk. Our analysis indicates that a small amount of 870 micron dust emission (~5 mJy) originates inside a large (40-50 AU radius) low optical depth cavity. This result can be interpreted either in the context of an abrupt decrease by a factor of ~5 in the radial distribution of millimeter-sized dust grains or as indirect evidence for a gap in the disk, in agreement with previous inferences from the unresolved infrared spectrum and scattered light images. A preliminary model focused on the latter possibility suggests the presence of a low-mass (planetary) companion, having properties commensurate with those inferred from the recent discovery of LkCa 15b.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Disk Imaging Survey of Chemistry with SMA: II. Southern Sky Protoplanetary Disk Data and Full Sample Statistics

This is the second in a series of papers based on data from DISCS, a Submillimeter Array observing program aimed at spatially and spectrally resolving the chemical composition of 12 protoplanetary disks. We present data on six Southern sky sources - IM Lup, SAO 206462 (HD 135344b), HD 142527, AS 209, AS 205 and V4046 Sgr - which complement the six sources in the Taurus star forming region reported previously. CO 2-1 and HCO+ 3-2 emission are detected and resolved in all disks and show velocity patterns consistent with Keplerian rotation. Where detected, the emission from DCO+ 3-2, N2H+ 3-2, H2CO 3-2 and 4-3,HCN 3-2 and CN 2-1 are also generally spatially resolved. The detection rates are highest toward the M and K stars, while the F star SAO 206462 has only weak CN and HCN emission, and H2CO alone is detected toward HD 142527. These findings together with the statistics from the previous Taurus disks, support the hypothesis that high detection rates of many small molecules depend on the presence of a cold and protected disk midplane, which is less common around F and A stars compared to M and K stars. Disk-averaged variations in the proposed radiation tracer CN/HCN are found to be small, despite two orders of magnitude range of spectral types and accretion rates. In contrast, the resolved images suggest that the CN/HCN emission ratio varies with disk radius in at least two of the systems. There are no clear observational differences in the disk chemistry between the classical/full T Tauri disks and transitional disks. Furthermore, the observed line emission does not depend on measured accretion luminosities or the number of infrared lines detected, which suggests that the chemistry outside of 100 AU is not coupled to the physical processes that drive the chemistry in the innermost few AU.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolving the CO Snow Line in the Disk around HD 163296

We report Submillimeter Array (SMA) observations of CO (J=2--1, 3--2 and 6--5) and its isotopologues (13CO J=2--1, C18O J=2--1 and C17O J=3--2) in the disk around the Herbig Ae star HD 163296 at ~2" (250 AU) resolution, and interpret these data in the framework of a model that constrains the radial and vertical location of the line emission regions. First, we develop a physically self-consistent accretion disk model with an exponentially tapered edge that matches the spectral energy distribution and spatially resolved millimeter dust continuum emission. Then, we refine the vertical structure of the model using wide range of excitation conditions sampled by the CO lines, in particular the rarely observed J=6--5 transition. By fitting 13CO data in this structure, we further constrain the vertical distribution of CO to lie between a lower boundary below which CO freezes out onto dust grains (T ~ 19 K) and an upper boundary above which CO can be photodissociated (the hydrogen column density from the disk surface is ~ 10^{21} cm-2). The freeze-out at 19 K leads to a significant drop in the gas-phase CO column density beyond a radius of ~155 AU, a "CO snow line" that we directly resolve. By fitting the abundances of all CO isotopologues, we derive isotopic ratios of 12C/13C, 16O/18O and 18O/17O that are consistent with quiescent interstellar gas-phase values. This detailed model of the HD 163296 disk demonstrates the potential of a staged, parametric technique for constructing unified gas and dust structure models and constraining the distribution of molecular abundances using resolved multi-transition, multi-isotope observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Ionization Fraction in the DM Tau Protoplanetary Disk

We present millimeter-wave observations of several molecular ions in the disk around the pre-main-sequence star DM Tau and use these to investigate the ionization fraction in different regions of the disk. New Submillimeter Array (SMA) observations of H2D+ J=1_10 - 1_11, N2H+ J=4-3 and CO J=3-2 are presented. H2D+ and N2H+ are not detected and using the CO 3-2 disk size the observations result in an upper limit of <0.47 K km s-1 for both lines, a factor of 2.5 below previous single-dish H2D+ observations. Assuming LTE, a disk midplane temperature of 10-20 K and estimates of the H2D+ o/p ratio, the observed limit corresponds to NH2D+ < 4 - 21 \times 1012 cm-2. We adopt a parametric model for the disk structure from the literature and use new IRAM 30 meter telescope observations of the H13CO+ J=3-2 line and previously published SMA observations of the N2H+ J=3-2, HCO+ J=3-2 and DCO+ J=3-2 lines to constrain the ionization fraction, xi, in three temperature regions in the disk where theoretical considerations suggest different ions should dominate: (1) a warm, upper layer with T>20 K where CO is in the gas-phase and HCO+ is most abundant, where we estimate xi \simeq 4 \times 10-10, (2) a cooler molecular layer with T = 16-20 K where N2H+ and DCO+ abundances are predicted to peak, with xi \simeq 3\times10-11, and (3) the cold, dense midplane with T<16 K where H3+ and its deuterated isotopologues are the main carriers of positive charge, with xi < 3\times10-10. While there are considerable uncertainties, these estimates are consistent with a decreasing ionization fraction into the deeper, colder, and denser disk layers. Stronger constraints on the ionization fraction in the disk midplane will require not only substantially more sensitive observations of the H2D+ 1_10 - 1_11 line, but also robust determinations of the o/p ratio, observations of D2H+ and stronger constraints on where N2 is present in the gas phase.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Empirical Constraints on Turbulence in Protoplanetary Accretion Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Millimeter Imaging of the beta Pictoris Debris Disk: Evidence for a Planetesimal Belt

We present observations at 1.3 millimeters wavelength of the beta Pictoris debris disk with beam size 4.3 x 2.6 arcsec (83 x 50 AU) from the Submillimeter Array. The emission shows two peaks separated by ~7 arsec along the disk plane, which we interpret as a highly inclined dust ring or belt. A simple model constrains the belt center to 94+/-8 AU, close to the prominent break in slope of the optical scattered light. We identify this region as the location as the main reservoir of dust producing planetesimals in the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protoplanetary Disk Structures in Ophiuchus II: Extension to Fainter Sources

We present new results from a significant extension of our previous high angular resolution (0.3" = 40 AU) Submillimeter Array survey of the 880 um continuum emission from dusty circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. An expanded sample is constructed to probe disk structures that emit significantly lower millimeter luminosities (hence dust masses), down to the median value for T Tauri stars. Using a Monte Carlo radiative transfer code, the millimeter visibilities and broadband spectral energy distribution for each disk are simultaneously reproduced with a two-dimensional parametric model for a viscous accretion disk. We find wide ranges of characteristic radii (14-198 AU) and disk masses (0.004-0.143 M_sun), but a narrow distribution of surface density gradients (0.4-1.1) that is consistent with a uniform value $γ$ = 0.9 +/- 0.2 and independent of mass (or millimeter luminosity). In this sample, we find a correlation between the disk luminosity/mass and characteristic radius, such that fainter disks are both smaller and less massive. We suggest that this relationship is an imprint of the initial conditions inherited by the disks at their formation epoch, compare their angular momenta with those of molecular cloud cores, and speculate on how future observations can help constrain the distribution of viscous evolution timescales. No other correlations between disk and star properties are found. The inferred disk structures are briefly compared with theoretical models for giant planet formation, although resolution limitations do not permit us to directly comment on material inside R = 20 AU. However, there is some compelling evidence for dust evolution in the planet formation region: 4/17 disks in the sample show resolved regions of significantly reduced optical depths within ~20-40 AU of their central stars.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Millimeter Dust Emission in the GQ Lup System

We present Submillimeter Array observations of the GQ Lup system at 1.3 millimeters wavelength with $0\farcs4$ ($\sim$60 AU) resolution. Emission is detected from the position of the primary star, GQ Lup A, and is marginally resolved. No emission is detected from the substellar companion, GQ Lup B, $0\farcs7$ away. These data, together with models of the spectral energy distribution, suggest a compact disk around GQ Lup A with mass $\sim 3$ M$_{Jup}$, perhaps truncated by tidal forces. There is no evidence for a gap or hole in the disk that might be the signature of an additional inner companion body capable of scattering GQ Lup B out to $\sim100$ AU separation from GQ Lup A. For GQ Lup B to have formed {\it in situ}, the disk would have to have been much more massive and extended.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stringent Limits on the Polarized Submillimeter Emission from Protoplanetary Disks

We present arcsecond-resolution Submillimeter Array (SMA) polarimetric observations of the 880 um continuum emission from the protoplanetary disks around two nearby stars, HD 163296 and TW Hydrae. Although previous observations and theoretical work have suggested that a 2-3% polarization fraction should be common for the millimeter continuum emission from such disks, we detect no polarized continuum emission above a 3-sigma upper limit of 7 mJy in each arcsecond-scale beam, or <1% in integrated continuum emission. We compare the SMA upper limits with the predictions from the exploratory Cho & Lazarian (2007) model of polarized emission from T Tauri disks threaded by toroidal magnetic fields, and rule out their fiducial model at the ~10-sigma level. We explore some potential causes for this discrepancy, focusing on model parameters that describe the shape, magnetic field alignment, and size distribution of grains in the disk. We also investigate related effects like the magnetic field strength and geometry, scattering off of large grains, and the efficiency of grain alignment, including recent advances in grain alignment theory, which are not considered in the fiducial model. We discuss the impact each parameter would have on the data and determine that the suppression of polarized emission plausibly arises from rounding of large grains, reduced efficiency of grain alignment with the magnetic field, and/or some degree of magnetic field tangling (perhaps due to turbulence). A poloidal magnetic field geometry could also reduce the polarization signal, particularly for a face-on viewing geometry like the TW Hya disk. The data provided here offer the most stringent limits to date on the polarized millimeter-wavelength emission from disks around young stars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Truncated Disks in TW Hya Association Multiple Star Systems

We present high angular resolution (down to 0.3" = 13 AU in diameter) Submillimeter Array observations of the 880 micron (340 GHz) thermal continuum emission from circumstellar dust disks in the nearby HD 98800 and Hen 3-600 multiple star systems. In both cases, the dust emission is resolved and localized around one stellar component - the HD 98800 B and Hen 3-600 A spectroscopic binaries - with no evidence for circum-system material. Using two-dimensional Monte Carlo radiative transfer calculations, we compare the SMA visibilities and broadband spectral energy distributions with truncated disk models to empirically locate the inner and outer edges of both disks. The HD 98800 B disk appears to be aligned with the spectroscopic binary orbit, is internally truncated at a radius of 3.5 AU, and extends to only 10-15 AU from the central stars. The Hen 3-600 A disk is slightly larger, with an inner edge at 1 AU and an outer radius of 15-25 AU. These inferred disk structures compare favorably with theoretical predictions of their truncation due to tidal interactions with the stellar companions.

0 citations0 reviewsNo code linkedOpen access