Source author record

David J. Wilner

David J. Wilner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

94works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

94 published item(s)

preprint2026arXiv

Isotopic Ratios in the Disk of HD 163296

Isotopic abundance ratios in protoplanetary disks are critical for understanding volatile inheritance and chemical evolution in planet-forming environments. We present Atacama Large Millimeter/submillimeter Array observations of the rare isotopologue 13C18O(2-1) at approximately 0.3 arcsec resolution from the disk around the Herbig Ae star HD 163296, combined with archival observations of C17O(2-1), C18O(1-0), and C17O(1-0), to empirically constrain carbon and oxygen isotopic ratios without detailed disk modeling. Both the C17O/13C18O(2-1) and C18O/C17O(1-0) flux ratios rise sharply across the CO snowline and flatten beyond 1.5 arcsec (r >= 150 au), where the emission becomes optically thin. This transition, reflecting a steep drop in CO column density set by the disk's thermal structure, makes HD 163296 an optimal case for isotopic analysis. Using beam-averaged intensities of the four transitions measured in this optically thin region, we derive isotopic ratios of 12C/13C = 75.3 (+14.7/-11.4) and 18O/17O = 3.28 (+0.31/-0.26), both consistent with local interstellar medium values. The 16O/18O ratio remains weakly constrained due to moderate optical depth in the C18O(1-0) line and degeneracy with CO column density. These results demonstrate that rare CO isotopologues can provide robust, empirical constraints on isotopic ratios in disks when sharp structural transitions allow for the identification of optically thin regions, and establish HD 163296 as a benchmark for extending such studies to other systems with resolved snowline structures.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IX: Gas-driven origin for the continuum arc in the debris disc of HD 121617

Debris discs were long considered to be largely gas-free environments governed by collisional fragmentation, gravitational stirring, and radiative forces. Recent CO detections show that gas is present, but its abundance and origin remain uncertain. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) revealed a narrow gas and dust ring in the disc HD 121617 with an asymmetric arc 40% brighter than the rest of the ring. We aim to constrain the total gas mass in HD 121617 assuming the dust arc is produced by hydrodynamical gas-dust interactions. We used the Dusty FARGO-ADSG code, modelling dust as Lagrangian particles, including radiation pressure and dust feedback, and varying the total gas mass. Simulations were compared to observations using radiative transfer. An unstable gas ring creates a size-dependent radial and azimuthal dust trap whose efficiency depends on gas mass. Two models, with 50 and 5 Earth masses of gas, reproduce both the ALMA band 7 arc and the outward offset of the VLT/SPHERE scattered-light ring via gas drag and radiation pressure. We infer a conservative gas-mass range of 2.5 to 250 Earth masses. If the ALMA asymmetry is caused by gas drag, the required gas mass compared with the observed CO implies substantial H2, consistent with primordial gas. HD 121617 would then be a hybrid disc between protoplanetary and debris stages. Since a planet could also create an arc, future observations are needed to distinguish these scenarios.

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.

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}}$.

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.

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.

preprint2022arXiv

Mapping NGC 7027 in New Light: CO$^+$ and HCO$^+$ Emission Reveal Its Photon- and X-ray-Dominated Regions

The young and well-studied planetary nebula NGC 7027 harbors significant molecular gas that is irradiated by luminous, point-like UV (central star) and diffuse (shocked nebular) X-ray emission. This nebula represents an excellent subject to investigate the molecular chemistry and physical conditions within photon- and X-ray-dominated regions (PDRs and XDRs). As yet, the exact formation routes of CO$^+$ and HCO$^+$ in PN environments remain uncertain. Here, we present $\sim$2$"$ resolution maps of NGC 7027 in the irradiation tracers CO$^+$ and HCO$^+$, obtained with the IRAM NOEMA interferometer, along with SMA CO and HST 2.12~$μ$m H$_2$ data for context. The CO$^+$ map constitutes the first interferometric map of this molecular ion in any PN. Comparison of CO$^+$ and HCO$^+$ maps reveal strikingly different emission morphologies, as well as a systematic spatial displacement between the two molecules; the regions of brightest HCO$^+$, found along the central waist of the nebula, are radially offset by $\sim$1$"$ ($\sim$900 au) outside the corresponding CO$^+$ emission peaks. The CO$^+$ emission furthermore precisely traces the inner boundaries of the nebula's PDR (as delineated by near-IR H$_2$ emission), suggesting that central star UV emission drives CO$^+$ formation. The displacement of HCO$^+$ radially outward with respect to CO$^+$ is indicative that dust-penetrating soft X-rays are responsible for enhancing the HCO$^+$ abundance in the surrounding molecular envelope, forming an XDR. These interferometric CO$^+$ and HCO$^+$ observations of NGC 7027 thus clearly establish the spatial distinction between the PDR and XDR formed (respectively) by intense UV and X-ray irradiation of molecular gas.

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.

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.

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.

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.

preprint2021arXiv

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: I. ALMA Observations at $λ= $1.4-0.4 mm and Continuum Analysis

We observed with ALMA three deeply buried nuclei in two galaxies, NGC 4418 and Arp 220, at $\sim$0.2$''$ resolution over a total bandwidth of 67 GHz in $f_{\rm rest}$ = 215 - 697 GHz. Here we (1) introduce our program, (2) describe our data reduction method for wide-band, high-resolution imaging spectroscopy, (3) analyze in visibilities the compact nuclei with line forests, (4) develop a continuum-based estimation method of dust opacity and gas column density in heavily obscured nuclei, which uses the BGN (buried galactic nuclei) model and is sensitive to $\log(N_{\rm H_2}/{\rm cm}^{-2}) \sim $ 25 - 26 at $λ\sim 1$ mm, and (5) present the continuum data and diagnosis of our targets. The three continuum nuclei have major-axis FWHM of $\sim$0.1$''$-0.3$''$ (20-140 pc) aligned to their rotating nuclear disks of molecular gas. However, each nucleus is described better with two or three concentric components than with a single Gaussian. The innermost cores have sizes of 0.05$''$-0.10$''$ (8-40 pc), peak brightness temperatures of ~100-500 K at 350 GHz, and more fractional flux at lower frequencies. The intermediate components correspond to the nuclear disks. They have axial ratios of $\approx$0.5 and hence inclinations $\stackrel{>}{\sim} 60$ deg. The outermost elements include the bipolar outflow from Arp 220W. We estimate 1 mm dust opacity of $τ_{\rm d,1mm} \approx 2.2$, $1.2$, and $\stackrel{<}{\sim} 0.4$ respectively for NGC 4418, Arp 220W, and Arp 220E. The first two correspond to $\log(N_{\rm H}/{\rm cm}^{-2}) \sim 26$ for conventional dust-opacity laws, and hence the nuclei are highly Compton thick.

preprint2021arXiv

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: II. Line Forests at $λ= $1.4--0.4 mm and Circumnuclear Gas Observed with ALMA

We present the line observations in our ALMA imaging spectral scan toward three deeply buried nuclei in NGC 4418 and Arp 220. We cover 67 GHz in $f_{\rm rest}$=215-697 GHz at about 0.2$"$ (30, 80 pc) resolution. All the nuclei show dense line forests; we report our initial line identification using 55 species. The line velocities generally indicate gas rotation around each nucleus, tracing nuclear disks of $\sim$100 pc sizes. We confirmed the counter-rotation of the nuclear disks in Arp 220 and that of the nuclear disk and the galactic disk in NGC 4418. While the brightest lines exceed 100 K, most of the major lines and many $^{13}$C isotopologues show absorption against even brighter continuum cores of the nuclei. The lines with higher upper-level energies, including those from vibrationally-excited molecules, tend to arise from smaller areas, indicating radially varying conditions in these nuclei. The outflows from the two Arp 220 nuclei cause blueshifted line absorption below the continuum level. The absorption mostly has small spatial offsets from the continuum peaks to indicate the outflow orientations. The bipolar outflow from the western nucleus is also imaged in multiple emission lines, showing the extent of $\sim$1$"$ (400 pc). Redshifted line absorption against the nucleus of NGC 4418 indicates either an inward gas motion or a small collimated outflow slanted to the nuclear disk. We also resolved some previous confusions due to line blending and misidentification.

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.

preprint2021arXiv

Hot corino chemistry in the Class I binary source Ser-emb 11

We report the detection of more than 120 emission lines corresponding to 8 complex organic molecules (CH3OH, CH3CH2OH, CH3OCH3, CH3OCHO, CH3COCH3, NH2CHO, CH2DCN, and CH3CH2CN) and 3 isotopologues (CH2DOH, 13CH3CN, and CH3C15N) toward the western component of the Ser-emb 11 binary young stellar object (YSO) using observations with the Atacama Large Millimeter/submillimeter Array at ~1 mm. The complex organic emission was unresolved with a ~0.5" beam (~220 au) in a compact region around the central protostar, and a population diagram analysis revealed excitation temperatures above 100 K for all COMs, indicating the presence of a hot corino. The estimated column densities were in the range of 10^17 - 10^18 cm^-2 for the O-bearing COMs, and three orders of magnitude lower for the N-bearing species. We also report the detection of H2CO and CH3OH emission in a nearby millimeter source that had not been previously catalogued. Ser-emb 11 is classified in the literature as a Class I source near the Class 0/I cutoff. The estimated COM relative abundances in Ser-emb 11 W and the other three Class I hot corino sources reported in the literature are consistent with those of Class 0 hot corinos, suggesting a continuity in the chemical composition of hot corinos during protostellar evolution.

preprint2020arXiv

A Dust Trap in the Young Multiple System HD 34700

Millimeter observations of disks around young stars reveal substructures indicative of gas pressure traps that may aid grain growth and planet formation. We present Submillimeter Array observations of HD 34700- two Herbig Ae stars in a close binary system (Aa/Ab, $\sim$0.25 AU), surrounded by a disk presenting a large cavity and spiral arms seen in scattered light, and two distant, lower mass companions. These observations include 1.3 mm continuum emission and the $^{12}$CO 2-1 line at $\sim0.5$" (178 AU) resolution. They resolve a prominent azimuthal asymmetry in the continuum, and Keplerian rotation of a circumbinary disk in the $^{12}$CO line. The asymmetry is located at a radius of $155^{+11}_{-7}$ AU, consistent with the edge of the scattered light cavity, being resolved in both radius ($72 ^{+14}_{-15}$ AU) and azimuth (FWHM = $64 ^{\circ +8}_{-7}$). The strong asymmetry in millimeter continuum emission could be evidence for a dust trap, together with the more symmetric morphology of $^{12}$CO emission and small grains. We hypothesize an unseen circumbinary companion, responsible for the cavity in scattered light and creating a vortex at the cavity edge that manifests in dust trapping. The disk mass has limitations imposed by the detection of $^{12}$CO and non-detection of $^{13}$CO. We discuss its consequences for the potential past gravitational instability of this system, likely accounting for the rapid formation of a circumbinary companion. We also report the discovery of resolved continuum emission associated with HD 34700B (projected separation $\sim1850$AU), which we explain through a circumstellar disk.

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."

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.

preprint2020arXiv

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

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

preprint2020arXiv

An evolutionary study of volatile chemistry in protoplanetary disks

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

preprint2020arXiv

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

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

preprint2020arXiv

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.

preprint2020arXiv

Dust Populations in the Iconic Vega Planetary System Resolved by ALMA

The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a steep inner edge (at 60-80 au). The belt is very broad, extending out to at least 150-200 au. We strongly detect the star and set a stringent upper limit to warm dust emission previously detected in the infrared. We discuss three scenarios that could explain the architecture of Vega's planetary system, including the new {ALMA} constraints: no outer planets, a chain of low-mass planets, and a single giant planet. The planet-less scenario is only feasible if the outer belt was born with the observed sharp inner edge. If instead the inner edge is currently being truncated by a planet, then the planet must be $\gtrsim$6 M$_{\oplus}$ and at $\lesssim71$ au to have cleared its chaotic zone within the system age. In the planet chain scenario, outward planet migration and inward scattering of planetesimals could produce the hot and warm dust observed in the inner regions of the system. In the single giant planet scenario, an asteroid belt could be responsible for the warm dust, and mean motion resonances with the planet could put asteroids on star-grazing orbits, producing the hot dust.

preprint2020arXiv

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

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

preprint2020arXiv

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.

preprint2020arXiv

Multiband GPI Imaging of the HR 4796A Debris Disk

We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2016arXiv

ALMA Observations of the Debris Disk of Solar Analogue Tau Ceti

We present 1.3 mm observations of the Sun-like star $τ$ Ceti with the Atacama Large Millimeter/submillimeter Array (ALMA) that probe angular scales of $\sim1$'' (4 AU). This first interferometric image of the $τ$ Ceti system, which hosts both a debris disk and possible multiplanet system, shows emission from a nearly face-on belt of cold dust with a position angle of $90^\circ$ surrounding an unresolved central source at the stellar position. To characterize this emission structure, we fit parametric models to the millimeter visibilities. The resulting best-fit model yields an inner belt edge of $6.2^{+9.8}_{-4.6}$ AU, consistent with inferences from lower resolution, far-infrared Herschel observations. While the limited data at sufficiently short baselines preclude us from placing stronger constraints on the belt properties and its relation to the proposed five planet system, the observations do provide a strong lower limit on the fractional width of the belt, $ΔR/R > 0.75$ with $99\%$ confidence. This fractional width is more similar to broad disks such as HD 107146 than narrow belts such as the Kuiper Belt and Fomalhaut. The unresolved central source has a higher flux density than the predicted flux of the stellar photosphere at 1.3 mm. Given previous measurements of an excess by a factor of $\sim2$ at 8.7 mm, this emission is likely due to a hot stellar chromosphere.

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.

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.

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.

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.

preprint2015arXiv

A Sub-arcsecond Survey Toward Class 0 Protostars in Perseus: Searching for Signatures of Protostellar Disks

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

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.

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.

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.

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.

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.

preprint2015arXiv

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

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

preprint2015arXiv

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.

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.

preprint2015arXiv

The Inner Debris Structure in the Fomalhaut Planetary System

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 μm (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 μJy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3σ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3σ upper limit for such a belt is <1.3 mJy at 870 μm. We also interpret the 2 and 8-13 μm interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

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.

preprint2014arXiv

Cold Molecular Gas in Merger Remnants: I. Formation of Molecular Gas Disks

We present < 1 kpc resolution CO imaging study of 37 optically-selected local merger remnants using new and archival interferometric maps obtained with ALMA, CARMA, SMA and PdBI. We supplement a sub-sample with single-dish measurements obtained at the NRO 45 m telescope for estimating the molecular gas mass (10^7 - 10^11 M_sun), and evaluating the missing flux of the interferometric measurements. Among the sources with robust CO detections, we find that 80 % (24/30) of the sample show kinematical signatures of rotating molecular gas disks (including nuclear rings) in their velocity fields, and the sizes of these disks vary significantly from 1.1 kpc to 9.3 kpc. The size of the molecular gas disks in 54 % of the sources is more compact than the K-band effective radius. These small gas disks may have formed from a past gas inflow that was triggered by a dynamical instability during a potential merging event. On the other hand, the rest (46 %) of the sources have gas disks which are extended relative to the stellar component, possibly forming a late-type galaxy with a central stellar bulge. Our new compilation of observational data suggests that nuclear and extended molecular gas disks are common in the final stages of mergers. This finding is consistent with recent major-merger simulations of gas rich progenitor disks. Finally, we suggest that some of the rotation-supported turbulent disks observed at high redshifts may result from galaxies that have experienced a recent major merger.

preprint2014arXiv

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

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

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.

preprint2014arXiv

The State of the Warm and Cold Gas in the Extreme Starburst at the Core of the Phoenix Galaxy Cluster (SPT-CLJ2344-4243)

[Abridged] We present new optical integral field spectroscopy (Gemini South) and submillimeter spectroscopy (Submillimeter Array) of the central galaxy in the Phoenix cluster (SPT-CLJ2344-4243). This cluster was previously reported to have a massive starburst (~800 Msun/yr) in the central, brightest cluster galaxy, most likely fueled by the rapidly-cooling intracluster medium. These new data reveal a complex emission-line nebula, extending for >30 kpc from the central galaxy. The total Halpha luminosity, assuming Halpha/Hbeta = 2.85, is L_Ha = 7.6 +/- 0.4 x10^43 erg/s, making this the most luminous emission line nebula detected in the center of a cool core cluster. Overall, the relative fluxes of the low-ionization lines (e.g., [O II], Hbeta) to the UV continuum are consistent with photoionization by young stars. In both the center of the galaxy and in a newly-discovered highly-ionized plume to the north of the galaxy, the ionization ratios are consistent with both shocks and AGN photoionization. We speculate that this extended plume may be a galactic wind, driven and partially photoionized by both the starburst and central AGN. We find evidence for shocks throughout the ISM of the central galaxy, most likely driven by a combination of stellar winds from massive young stars, core-collapse supernovae, and the central AGN. In addition to the warm, ionized gas, we detect a substantial amount of cold, molecular gas via the CO(3-2) transition, coincident in position with the galaxy center. We infer a molecular gas mass of M_H2 = 2.2 +/- 0.6 x10^10 Msun, which implies that the starburst will consume its fuel in ~30 Myr if it is not replenished. The combination of the high level of turbulence in the warm phase and the high L_IR/M_H2 ratio suggests that this violent starburst may be in the process of quenching itself.

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%).

preprint2013arXiv

Constraining a Model of Turbulent Coronal Heating for AU Microscopii with X-Ray, Radio, and Millimeter Observations

Many low-mass pre-main-sequence stars exhibit strong magnetic activity and coronal X-ray emission. Even after the primordial accretion disk has been cleared out, the star's high-energy radiation continues to affect the formation and evolution of dust, planetesimals, and large planets. Young stars with debris disks are thus ideal environments for studying the earliest stages of non-accretion-driven coronae. In this paper we simulate the corona of AU Mic, a nearby active M dwarf with an edge-on debris disk. We apply a self-consistent model of coronal loop heating that was derived from numerical simulations of solar field-line tangling and magnetohydrodynamic turbulence. We also synthesize the modeled star's X-ray luminosity and thermal radio/millimeter continuum emission. A realistic set of parameter choices for AU Mic produces simulated observations that agree with all existing measurements and upper limits. This coronal model thus represents an alternative explanation for a recently discovered ALMA central emission peak that was suggested to be the result of an inner "asteroid belt" within 3 AU of the star. However, it is also possible that the central 1.3 mm peak is caused by a combination of active coronal emission and a bright inner source of dusty debris. Additional observations of this source's spatial extent and spectral energy distribution at millimeter and radio wavelengths will better constrain the relative contributions of the proposed mechanisms.

preprint2013arXiv

Imaging of the CO Snow Line in a Solar Nebula Analog

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

preprint2013arXiv

Modeling the Resolved Disk Around the Class 0 Protostar L1527

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at 870 micron and 3.4 mm from the Submillimeter Array (SMA) and Combined Array for Research in Millimeter Astronomy (CARMA). We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M_sun, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed; alpha ~ 2, implying a dust opacity spectral index beta ~ 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L\arcmin\ image, spectral energy distribution, and mid-infrared spectrum. The best fitting model has a disk radius of R = 125 AU, is highly flared (H ~ R^1.3), has a radial density profile rho ~ R^-2.5, and has a mass of 0.0075 M_sun. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 is comparable to disks around pre-main sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of ~1 micron-sized dust grains.

preprint2013arXiv

Resolving the gap and AU-scale asymmetries in the pre-transitional disk of V1247 Orionis

Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating the presence of young planets in these systems. In order to explore the structure of these objects and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This allows us spatially resolve the AU-scale disk structure from near- to mid-infrared wavelengths (1.5 to 13 μm), tracing material at different temperatures and over a wide range of stellocentric radii. Our observations reveal a narrow, optically-thick inner-disk component (located at 0.18 AU from the star) that is separated from the optically thick outer disk (radii >46 AU), providing unambiguous evidence for the existence of a gap in this pre-transitional disk. Surprisingly, we find that the gap region is filled with significant amounts of optically thin material with a carbon-dominated dust mineralogy. The presence of this optically thin gap material cannot be deduced solely from the spectral energy distribution, yet it is the dominant contributor at mid-infrared wavelengths. Furthermore, using Keck/NIRC2 aperture masking observations in the H, K', and L' band, we detect asymmetries in the brightness distribution on scales of about 15-40 AU, i.e. within the gap region. The detected asymmetries are highly significant, yet their amplitude and direction changes with wavelength, which is not consistent with a companion interpretation but indicates an inhomogeneous distribution of the gap material. We interpret this as strong evidence for the presence of complex density structures, possibly reflecting the dynamical interaction of the disk material with sub-stellar mass bodies that are responsible for the gap clearing.

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.

preprint2013arXiv

Submillimeter Interferometry of the Luminous Infrared Galaxy NGC 4418: A Hidden Hot Nucleus with an Inflow and an Outflow

We have observed the nucleus of the nearby luminous infrared galaxy NGC 4418 with subarcsec resolution at 860 and 450 micron for the first time to characterize its hidden power source. A ~20 pc (0.1") hot dusty core was found inside a 100 pc scale concentration of molecular gas at the galactic center. The 860 micron continuum core has a deconvolved (peak) brightness temperature of 120-210 K. The CO(3-2) peak brightness temperature there is as high as 90 K at 50 pc resolution. The core has a bolometric luminosity of about 10^{11} Lsun, which accounts for most of the galaxy luminosity. It is Compton thick (N_H >~ 10^{25} cm^{-2}) and has a high luminosity-to-mass ratio ~500 Lsun/Msun as well as a high luminosity surface density 10^{8.5+-0.5} Lsun pc^{-2}. These parameters are consistent with an AGN to be the main luminosity source (with an Eddington ratio about 0.3) while they can be also due to a young starburst near its maximum L/M. We also found an optical color (reddening) feature that we attribute to an outflow cone emanating from the nucleus. The hidden hot nucleus thus shows evidence of both an inflow, previously seen with absorption lines, and the new outflow reported here in a different direction. The nucleus must be rapidly evolving with these gas flows.

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).

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.

preprint2013arXiv

VLA and CARMA Observations of Protostars in the Cepheus Clouds: Sub-arcsecond Proto-Binaries Formed via Disk Fragmentation

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $λ= 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $λ= 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $λ= 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

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.

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.

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.

preprint2012arXiv

ALMA and VLA observations of the outflows in IRAS 16293-2422

We present ALMA and VLA observations of the molecular and ionized gas at 0.1-0.3 arcsec resolution in the Class 0 protostellar system IRAS 16293-2422. These data clarify the origins of the protostellar outflows from the deeply embedded sources in this complex region. Source A2 is confirmed to be at the origin of the well known large scale north-east--south-west flow. The most recent VLA observations reveal a new ejection from that protostar, demonstrating that it drives an episodic jet. The central compact part of the other known large scale flow in the system, oriented roughly east-west, is well delineated by the CO(6-5) emission imaged with ALMA and is confirmed to be driven from within component A. Finally, a one-sided blueshifted bubble-like outflow structure is detected here for the first time from source B to the north-west of the system. Its very short dynamical timescale (~ 200 yr), low velocity, and moderate collimation support the idea that source B is the youngest object in the system, and possibly one of the youngest protostars known.

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.

preprint2012arXiv

Evidence for Multiple Pathways to Deuterium Enhancements in Protoplanetary Disks

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

preprint2012arXiv

Grain growth signatures in the protoplanetary discs of Chamaeleon and Lupus

We present ATCA results of a 3 and 7 mm continuum survey of 20 T Tauri stars in the Chamaeleon and Lupus star forming regions. This survey aims to identify protoplanetary discs with signs of grain growth. We detected 90% of the sources at 3 and 7 mm, and determined the spectral slopes, dust opacity indices and dust disc masses. We also present temporal monitoring results of a small sub-set of sources at 7, 15 mm and 3+6 cm to investigate grain growth to cm sizes and constrain emission mechanisms in these sources. Additionally, we investigated the potential correlation between grain growth signatures in the infrared (10 μm silicate feature) and millimetre (1-3 mm spectral slope, α). Eleven sources at 3 and 7 mm have dominant thermal dust emission up to 7 mm, with 7 of these having a 1-3 mm dust opacity index less than unity, suggesting grain growth up to at least mm sizes. The Chamaeleon sources observed at 15 mm and beyond show the presence of excess emission from an ionised wind and/or chromo- spheric emission. Long-timescale monitoring at 7 mm indicated that cm-sized pebbles are present in at least four sources. Short-timescale monitoring at 15 mm suggests the excess emission is from thermal free-free emission. Finally, a weak correlation was found between the strength of the 10 \mum feature and α, suggesting simultaneous dust evolution of the inner and outer parts of the disc. This survey shows that grain growth up to cm-sized pebbles and the presence of excess emission at 15 mm and beyond are common in these systems, and that temporal monitoring is required to disentangle these emission mechanisms.

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.

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.

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.

preprint2012arXiv

The Magnetized Environment of the W3(H2O) Protostars

We present the first interferometric polarization map of the W3(OH) massive star-forming region observed with the Submillimeter Array (SMA) at 878 mum with an angular resolution of 1.5 (about 3 \times 10 AU). Polarization is detected in the W3(H2O) hot core, an extended emission structure in the north-west of W3(H2O), and part of the W3(OH) ultracompact HII region. The W3(H2O) hot core is known to be associated with a synchrotron jet along the east-west direction. In this core, the inferred magnetic field orientation is well aligned with the synchrotron jet and close to the plane of sky. Using the Chandrasekhar-Fermi method with the observed dispersion in polarization angle, we estimate a plane-of-sky magnetic field strength of 17.0 mG. Combined with water maser Zeeman measurements, the total magnetic field strength is estimated to be 17.1 mG, comparable to the field strength estimated from the synchrotron model. The magnetic field energy dominates over turbulence in this core. In addition, the depolarization effect is discerned in both SMA and JCMT measurements. Despite the great difference in angular resolutions and map extents, the polarization percentage shows a similar power-law dependence with the beam averaged column density. We suggest that the column density may be an important factor to consider when interpreting the depolarization effect.

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.

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.

preprint2011arXiv

Initial conditions for star formation in clusters: physical and kinematical structure of the starless core Oph A-N6

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

preprint2011arXiv

Ionized Gas Kinematics and Morphology in Sgr B2 Main on 1000 AU Scales

We have imaged the Sgr B2 Main region with the Very Large Array in the BnA configuration ($θ_{beam}$ = 0\farcs13) in both the H52$α$ (45.453 GHz) radio recombination line (RRL) and 7 mm continuum emission. At a distance of 8500 pc, this spatial resolution corresponds to a physical scale of 0.005 pc ($\sim$1100 AU). The current observations detect H52$α$ emission in 12 individual ultracompact (UC) and hypercompact (HC) HII regions. Two of the sources with detected H52 $α$ emission have broad ($Δ$V$_{FWHM}\sim$50 \kms) recombination lines, and two of the sources show lines with peaks at more than one velocity. We use line parameters from the H52$α$ lines and our previous H66$α$ line observations to determine the relative contribution of thermal, pressure and kinematic broadening, and electron density. These new observations suggest that pressure broadening can account for the broad lines in some of the sources, but that gas motions (e.g. turbulence, accretion or outflow) contribute significantly to the broad lines in at least one of the sources (Sgr B2 F3).

preprint2011arXiv

Millimeter Observations of the Transition Disk around HD135344B (SAO 206462)

We present \sim1" resolution 1.3 millimeter dust continuum and spectral line (12CO and 13CO J=2-1) observations of the transitional disk system HD 135344B obtained with the Submillimeter Array. The disk shows a Keplerian rotation pattern with an inclination of \sim11°, based on the spatially and spectrally resolved 12CO and 13CO emission. The data show clear evidence for both dust and gas surface density reductions in the inner region of the disk (radius \lesssim 50 AU) from the continuum and 13CO J=2-1 data, respectively. The presence of this inner cavity in both the dust and gas is more consistent with clearing by giant planet formation than by photoevaporation or by grain growth. There is also an indication of global CO gas depletion in the disk, as the mass estimated from 13CO emission (\sim3.8 \times 10^(-4) M\odot) is about two orders of magnitude lower than that derived from the 1.3 mm continuum (\sim2.8 \times 10^(-2) M\odot).

preprint2011arXiv

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

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

preprint2011arXiv

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.

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.

preprint2011arXiv

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

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

preprint2010arXiv

[CII] line emission in BRI1335-0417 at z=4.4

Using the 12m APEX telescope, we have detected redshifted emission from the 157.74micron [CII] line in the z=4.4074 quasar BRI1335-0417. The linewidth and redshift are in good agreement with previous observations of high-J CO line emission. We measure a [CII] line luminosity, L_[CII] = (16.4 +/- 2.6)x10^9 Lsun, making BRI~1335-0417 the most luminous, unlensed [CII] line emitter known at high-redshift. The [CII]-to-FIR luminosity ratio of (5.3+/-0.8)x10^-4 is ~3x higher than expected for an average object with a FIR luminosity L_FIR = 3.1x10^13 Lsun, if this ratio were to follow the trend observed in other FIR-bright galaxies that have been detected in [CII] line emission. These new data suggest that the scatter in the [CII]-to-FIR luminosity ratio could be larger than previously expected for high luminosity objects. BR1335-0417 has a similar FIR luminosity and [CII]/CO luminosity compared to local ULIRGS and appears to be a gas-rich merger forming stars at a rate of a few thousand solar masses per year.

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.

preprint2010arXiv

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

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

preprint2010arXiv

Imaging Carbon Monoxide Emission in the Starburst Galaxy NGC 6000

We present measurements of carbon monoxide emission in the central region of the nearby starburst NGC 6000 taken with the Submillimeter Array. The J=2-1 transition of 12CO, 13CO, and C18O were imaged at a resolution of ~3''x2'' (450x300 pc). We accurately determine the dynamical center of NGC 6000 at R.A(J2000.0)=15h49m49.5s and dec(J2000.0)=-29d23'13'' which agrees with the peak of molecular emission position. The observed CO dynamics could be explained in the context of the presence of a bar potential affecting the molecular material, likely responsible for the strong nuclear concentration where more than 85% of the gas is located. We detect a kinematically detached component of dense molecular gas at relatively high velocity which might be fueling the star formation. A total nuclear dynamical mass of 7x10^9 Msun is derived and a total mass of gas of 4.6x10^8 Msun, yielding a Mgas/Mdyn~6%, similar to other previously studied barred galaxies with central starbursts. We determined the mass of molecular gas with the optically thin isotopologue C18O and we estimate a CO-to-H2 conversion factor X(CO)=0.4x10^20 cm-2/(K km s-1) in agreement with that determined in other starburst galaxies.

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.

preprint2010arXiv

Observational Constraints on Submillimeter Dust Opacity

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (\kaprat) and the submillimeter opacity power-law index ($κ\propto λ^{-β}$). Using the average value of theoretical dust opacity models at 2.2 \micron, we constrain the dust opacity at 850 and 450 \micron . Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are $\frac{κ_{850}}{κ_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{κ_{450}}{κ_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $β_{smm} = (2.18 - 2.58)^{+0.30}_{-0.30}$. The range of quoted values are determined from the uncertainty in the physical model for B335. For an average 2.2 \micron\ opacity of $3800 \pm 700$ cm$^2$g$^{-1}$, we find a dust opacity at 850 and 450 \micron\ of $κ_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $κ_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) κ^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron.

preprint2010arXiv

Vibrationally Excited HCN in the Luminous Infrared Galaxy NGC 4418

Infrared pumping and its effect on the excitation of HCN molecules can be important when using rotational lines of HCN to probe dense molecular gas in galaxy nuclei. We report the first extragalactic detection of (sub)millimeter rotational lines of vibrationally excited HCN, in the dust-enshrouded nucleus of the luminous infrared galaxy NGC 4418. We estimate the excitation temperature of T_vib ~ 230 K between the vibrational ground and excited (v_2=1) states. This excitation is most likely due to infrared radiation. At this high vibrational temperature the path through the v_2=1 state must have a strong impact on the rotational excitation in the vibrational ground level, although it may not be dominant for all rotational levels. Our observations also revealed nearly confusion limited lines of CO, HCN, HCO+, H13CN, HC15N, CS, N2H+, and HC3N at lambda ~ 1 mm. Their relative intensities may also be affected by the infrared pumping.

preprint2009arXiv

A Spatially Resolved Inner Hole in the Disk around GM Aurigae

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

preprint2009arXiv

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.

preprint2009arXiv

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

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

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.

preprint2008arXiv

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

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

preprint2008arXiv

SMA CO(J=6-5) and 435 micron interferometric imaging of the nuclear region of Arp 220

We have used the Submillimeter Array (SMA) to make the first interferometric observations (beam size ~1") of the 12CO J=6-5 line and 435 micron (690 GHz) continuum emission toward the central region of the nearby ULIRG Arp 220. These observations resolve the eastern and western nuclei from each other, in both the molecular line and dust continuum emission. At 435 micron, the peak intensity of the western nucleus is stronger than the eastern nucleus, and the difference in peak intensities is less than at longer wavelengths. Fitting a simple model to the dust emission observed between 1.3 mm and 435 micron suggests that dust emissivity power law index in the western nucleus is near unity and steeper in the eastern nucleus, about 2, and that the dust emission is optically thick at the shorter wavelength. Comparison with single dish measurements indicate that the interferometer observations are missing ~60% of the dust emission, most likely from a spatially extended component to which these observations are not sensitive. The 12CO J=6-5 line observations clearly resolve kinematically the two nuclei. The distribution and kinematics of the 12CO J=6-5 line appear to be very similar to lower J CO lies observed at similar resolution. Analysis of multiple 12CO line intensities indicates that the molecular gas in both nuclei have similar excitation conditions, although the western nucleus is warmer and denser. The excitation conditions are similar to those found in other extreme environments, including M82, Mrk 231, and BR 1202-0725. Simultaneous lower resolution observations of the 12CO, 13CO, and C18O J=2-1 lines show that the 13CO and C18O lines have similar intensities, which suggests that both of these lines are optically thick, or possibly that extreme high mass star formation has produced in an overabundance of C18O.

preprint2007arXiv

PROSAC: A Submillimeter Array Survey of Low-Mass Protostars. I. Overview of Program: Envelopes, Disks, Outflows and Hot Cores

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.