sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Andrews, Sean M.

Normalized to: Andrews, S.

198 article(s) in total. 684 co-authors, from 1 to 116 common article(s). Median position in authors list is 4,0.

[1]  oai:arXiv.org:2007.02974  [pdf] - 2129498
Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup
Comments: 32 pages, 13 figures, accepted by ApJ. Self-calibrated visibilities, images, and reduction scripts are available at https://zenodo.org/record/3928862
Submitted: 2020-07-06
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.
[2]  oai:arXiv.org:2006.16187  [pdf] - 2123771
An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks
Comments: 20 pages, 11 figures
Submitted: 2020-06-29
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.
[3]  oai:arXiv.org:2006.12584  [pdf] - 2120122
An evolutionary study of volatile chemistry in protoplanetary disks
Comments: Accepted to ApJ
Submitted: 2020-06-22
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.
[4]  oai:arXiv.org:2006.03120  [pdf] - 2109248
Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks
Comments: accepted for publication in ApJ
Submitted: 2020-06-04
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 ($\alpha_{\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 ($\tau\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.
[5]  oai:arXiv.org:2004.13551  [pdf] - 2085331
Protoplanetary disk masses in NGC 2024: Evidence for two populations
Comments: Accepted for publication in A&A; 19 pages, 13 figures. Key results shown in Fig. 8 & 10
Submitted: 2020-04-28
Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.
[6]  oai:arXiv.org:2004.12176  [pdf] - 2110030
Measuring turbulent motion in planet-forming disks with ALMA: A detection around DM Tau and non-detections around MWC 480 and V4046 Sgr
Comments: 21 pages, 10 figures, accepted to ApJ
Submitted: 2020-04-25, last modified: 2020-04-28
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 ($\alpha\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.
[7]  oai:arXiv.org:2004.01204  [pdf] - 2075055
Disk tearing in a young triple star system with misaligned disk/orbit planes
Comments: 74 pages, 3+19 Figures, 6 Tables, paper under review at 'Science' (version prior to peer review)
Submitted: 2020-04-02
In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.
[8]  oai:arXiv.org:2003.08658  [pdf] - 2089156
ALMA 0.88 mm Survey of Disks around Planetary-mass Companions
Comments: Accepted for publication in AJ
Submitted: 2020-03-19, last modified: 2020-03-20
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$\sigma$ limits of 120-210 $\mu$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 $\mu$Jy.
[9]  oai:arXiv.org:2002.12525  [pdf] - 2056378
An ALMA Survey of H$_2$CO in Protoplanetary Disks
Comments: 20 pages, 14 figures, 3 figure sets, 9 tables. Published in ApJ (February 2020)
Submitted: 2020-02-27
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.
[10]  oai:arXiv.org:2001.11040  [pdf] - 2061705
A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk
Comments: 31 pages, 15 figures, accepted by ApJ. Self-calibration and imaging scripts, images, and self-calibrated visibilities available at https://zenodo.org/record/3628656#.XjEXyRNKiYV
Submitted: 2020-01-29
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 $\alpha$ 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."
[11]  oai:arXiv.org:2001.07720  [pdf] - 2050390
Nine localised deviations from Keplerian rotation in the DSHARP circumstellar disks: Kinematic evidence for protoplanets carving the gaps
Comments: Accepted for publication in ApJ letters, 10 pages, 3 figures
Submitted: 2020-01-21
We present evidence for localised deviations from Keplerian rotation, i.e., velocity "kinks", in 8 of 18 circumstellar disks observed by the DSHARP program: DoAr 25, Elias 2-27, GW Lup, HD 143006, HD 163296, IM Lup, Sz 129 and WaOph 6. Most of the kinks are detected over a small range in both radial extent and velocity, suggesting a planetary origin, but for some of them foreground contamination prevents us from measuring their spatial and velocity extent. Because of the DSHARP limited spectral resolution and signal-to-noise in the 12CO J=2-1 line, as well as cloud contamination, the kinks are usually detected in only one spectral channel, and will require confirmation. The strongest circumstantial evidence for protoplanets in the absence of higher spectral resolution data and additional tracers is that, upon deprojection, we find that all of the candidate planets lie within a gap and/or at the end of a spiral detected in dust continuum emission. This suggests that a significant fraction of the dust gaps and spirals observed by ALMA in disks are caused by embedded protoplanets.
[12]  oai:arXiv.org:2001.05007  [pdf] - 2032031
Observations of Protoplanetary Disk Structures
Comments: invited review to be published in Annual Reviews of Astronomy and Astrophysics; 43 pages, 13 figures
Submitted: 2020-01-14
The disks that orbit young stars are the essential conduits and reservoirs of material for star and planet formation. Their structures, meaning the spatial variations of the disk physical conditions, reflect the underlying mechanisms that drive those formation processes. Observations of the solids and gas in these disks, particularly at high resolution, provide fundamental insights on their mass distributions, dynamical states, and evolutionary behaviors. Over the past decade, rapid developments in these areas have largely been driven by observations with the Atacama Large Millimeter/submillimeter Array (ALMA). This review highlights the state of observational research on disk structures, emphasizing three key conclusions that reflect the main branches of the field: (1) Relationships among disk structure properties are also linked to the masses, environments, and evolutionary states of their stellar hosts; (2) There is clear, qualitative evidence for the growth and migration of disk solids, although the implied evolutionary timescales suggest the classical assumption of a smooth gas disk is inappropriate; and (3) Small-scale substructures with a variety of morphologies, locations, scales, and amplitudes -- presumably tracing local gas pressure maxima -- broadly influence the physical and observational properties of disks. The last point especially is reshaping the field, with the recognition that these disk substructures likely trace active sites of planetesimal growth or are the hallmarks of planetary systems at their formation epoch.
[13]  oai:arXiv.org:1912.10542  [pdf] - 2050320
The Nucleosynthetic Yields of Core-Collapse Supernovae, prospects for the Next Generation of Gamma-Ray Astronomy
Comments: 17 pages, 19 figures, submitted to ApJ
Submitted: 2019-12-22
Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants.
[14]  oai:arXiv.org:1911.12434  [pdf] - 2042280
Midplane temperature and outer edge of the protoplanetary disk around HD 163296
Comments: Accepted for publication in Astronomy and Astrophysics
Submitted: 2019-11-27
Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, also in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2-1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be an indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.
[15]  oai:arXiv.org:1911.04214  [pdf] - 2026299
Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?
Comments:
Submitted: 2019-11-11
We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.
[16]  oai:arXiv.org:1909.04674  [pdf] - 1975444
The DSHARP Rings: Evidence of Ongoing Planetesimal Formation?
Comments: Accepted for publication in ApJL
Submitted: 2019-09-10
Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process. We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations. We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.
[17]  oai:arXiv.org:1907.10647  [pdf] - 1966895
Probing CO and N$_2$ Snow Surfaces in Protoplanetary Disks with N$_2$H$^+$ Emission
Comments: Accepted for publication in ApJ
Submitted: 2019-07-24
Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$''$-0.4$''$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks.
[18]  oai:arXiv.org:1906.03269  [pdf] - 1966797
The Degree of Alignment Between Circumbinary Disks and Their Binary Hosts
Comments: ApJ accepted. 30 pages, 16 figures. Transcription error in Table 4 fixed
Submitted: 2019-06-07, last modified: 2019-06-17
All four circumbinary (CB) protoplanetary disks orbiting short-period ($P < 20$ day) double-lined spectroscopic binaries (SB2s)---a group that includes UZ Tau E, for which we present new ALMA data---exhibit sky-plane inclinations $i_{\rm disk}$ which match, to within a few degrees, the sky-plane inclinations $i_\star$ of their stellar hosts. Although for these systems the true mutual inclinations $\theta$ between disk and binary cannot be directly measured because relative nodal angles are unknown, the near-coincidence of $i_{\rm disk}$ and $i_\star$ suggests that $\theta$ is small for these most compact of systems. We confirm this hypothesis using a hierarchical Bayesian analysis, showing that 68% of CB disks around short-period SB2s have $\theta < 3.0^\circ$. Near co-planarity of CB disks implies near co-planarity of CB planets discovered by Kepler, which in turn implies that the occurrence rate of close-in CB planets is similar to that around single stars. By contrast, at longer periods ranging from $30-10^5$ days (where the nodal degeneracy can be broken via, e.g., binary astrometry), CB disks exhibit a wide range of mutual inclinations, from co-planar to polar. Many of these long-period binaries are eccentric, as their component stars are too far separated to be tidally circularized. We discuss how theories of binary formation and disk-binary gravitational interactions can accommodate all these observations.
[19]  oai:arXiv.org:1904.02127  [pdf] - 1890422
One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering
Comments: 15 pages, 11 figures, Accepted for publication in ApJL
Submitted: 2019-04-03, last modified: 2019-05-07
ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.
[20]  oai:arXiv.org:1905.00909  [pdf] - 1877423
A dust and gas cavity in the disc around CQ Tau revealed by ALMA
Comments:
Submitted: 2019-05-02
The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.
[21]  oai:arXiv.org:1904.09866  [pdf] - 1890482
Constraints on a Putative Planet Sculpting the V4046 Sagittarii Circumbinary Disk
Comments: 13 pages, 5 figures, Accepted
Submitted: 2019-04-18
We analyze the highest-resolution millimeter continuum and near-infrared (NIR) scattered-light images presented to date of the circumbinary disk orbiting V4046 Sgr, a ~20 Myr old actively accreting, close binary T Tauri star system located a mere 72.4 pc from Earth. We observed the disk with the Atacama Large Millimeter/submillimeter Array (ALMA) at 870-micron during Cycle 4, and we analyze these data in conjunction with archival NIR (H band) polarimetric images obtained with SPHERE/IRDIS on the ESO Very Large Telescope. At 0.3'' (20 au) resolution, the 870-micron image reveals a marginally resolved ring that peaks at ~32 au and has an extension of ~ 90 au. We infer a lower limit on dust mass of ~ 60.0 M_earth within the 870-micron ring, and confirm that the ring is well aligned with the larger-scale gaseous disk. A second, inner dust ring is also tentatively detected in the ALMA observations; its position appears coincident with the inner (~14 au radius) ring detected in scattered light. Using synthetic 870 micron and H-band images obtained from disk-planet interaction simulations, we attempt to constrain the mass of the putative planet orbiting at 20 au. Our trials suggest that a circumbinary Jovian-mass planet may be responsible for generating the dust ring and gap structures detected within the disk. We discuss the longevity of the gas-rich disk orbiting V4046 Sgr in the context of the binary nature of the system.
[22]  oai:arXiv.org:1904.00027  [pdf] - 1871561
The Mass of Stirring Bodies in the AU Mic Debris Disk Inferred from Resolved Vertical Structure
Comments: 21 pages, 8 figures, 5 tables. Accepted for publication in ApJ
Submitted: 2019-03-29, last modified: 2019-04-08
The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.
[23]  oai:arXiv.org:1903.11692  [pdf] - 1857802
Realizing the Unique Potential of ALMA to Probe the Gas Reservoir of Planet Formation
Comments: Astro2020 White Paper
Submitted: 2019-03-27
Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.
[24]  oai:arXiv.org:1903.06527  [pdf] - 1850019
On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits
Comments: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584
Submitted: 2019-03-13
Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.
[25]  oai:arXiv.org:1903.05077  [pdf] - 1847123
Protoplanetary Disk Science Enabled by Extremely Large Telescopes
Comments: Science white paper submitted to the Astro2020 Decadal Survey
Submitted: 2019-03-12
The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.
[26]  oai:arXiv.org:1812.04041  [pdf] - 1814209
The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures
Comments: 33 pages, 12 figures, minor typos corrected
Submitted: 2018-12-10, last modified: 2019-01-14
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.
[27]  oai:arXiv.org:1812.04193  [pdf] - 1815119
The Disk Substructures at High Angular Resolution Project (DSHARP): III. Spiral Structures in the Millimeter Continuum of the Elias 27, IM Lup, and WaOph 6 Disks
Comments: 19 pages, 9 figures, minor typos corrected
Submitted: 2018-12-10, last modified: 2019-01-14
We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks.
[28]  oai:arXiv.org:1901.02467  [pdf] - 1838320
Multiple spiral arms in the disk around intermediate-mass binary HD 34700A
Comments: Accepted by the Astrophysical Journal on Dec 15, 2018
Submitted: 2019-01-08
We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.
[29]  oai:arXiv.org:1812.04536  [pdf] - 1811254
The Disk Substructures at High Angular Resolution Project (DSHARP): IV. Characterizing substructures and interactions in disks around multiple star systems
Comments: 15 pages, 10 figures, accepted to ApJ Letters
Submitted: 2018-12-11
To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symmetric spiral arms in the disk around AS 205 N, with pitch angle of $14^\circ$, while the southern component AS 205 S, itself a spectroscopic binary, is surrounded by a compact inner disk and a bright ring at a radius of $34\,$au. The $^{12}$CO line exhibits clear signatures of tidal interactions, with spiral arms, extended arc-like emission, and high velocity gas, possible evidence of a recent close encounter between the disks in the AS 205 system, as these features are predicted by hydrodynamic simulations of fly-by encounters. In the HT Lup system, we detect continuum emission from all three components. The primary disk, HT Lup A, also shows two-armed symmetric spiral structure with a pitch angle of $4^\circ$, while HT Lup B and C, located at $25$ and $434\,$au in projected separation from HT Lup A, are barely resolved with $\sim5$ and $\sim10\,$au in diameter, respectively. The gas kinematics for the closest pair indicates a different sense of rotation for each disk, which could be explained by either a counter rotation of the two disks in different, close to parallel, planes, or by a projection effect of these disks with a close to $90^\circ$ misalignment between them.
[30]  oai:arXiv.org:1812.04043  [pdf] - 1811243
The Disk Substructures at High Angular Resolution Project (DSHARP): V. Interpreting ALMA maps of protoplanetary disks in terms of a dust model
Comments: ApJ Letters, in press
Submitted: 2018-12-10
The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_\nu$ to dust surface density $\Sigma_d$.
[31]  oai:arXiv.org:1812.04040  [pdf] - 1811241
The Disk Substructures at High Angular Resolution Project (DSHARP): I. Motivation, Sample, Calibration, and Overview
Comments: ApJ Letters, in press; some figure resolutions reduced, full-res version at https://www.cfa.harvard.edu/~sandrews/DSHARP_I.pdf ; Data Release at https://almascience.org/alma-data/lp/DSHARP
Submitted: 2018-12-10
We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($<$10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data -- including visibilities, images, calibration scripts, and more -- are released for community use at https://almascience.org/alma-data/lp/DSHARP.
[32]  oai:arXiv.org:1812.04047  [pdf] - 1811247
The Disk Substructures at High Angular Resolution Project (DSHARP) - IX. A high definition study of the HD 163296 planet forming disk
Comments: ApJ Letters, in press. 19 pages, 11 figures
Submitted: 2018-12-10
ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets.
[33]  oai:arXiv.org:1812.04049  [pdf] - 1811248
The Disk Substructures at High Angular Resolution Project (DSHARP): X. Multiple rings, a misaligned inner disk, and a bright arc in the disk around the T Tauri star HD 143006
Comments: 15 pages, 8 figures, accepted for publication at ApJ Letters
Submitted: 2018-12-10
We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006.
[34]  oai:arXiv.org:1812.04044  [pdf] - 1811244
The Disk Substructures at High Angular Resolution Project (DSHARP) VI: Dust trapping in thin-ringed protoplanetary disks
Comments: Accepted for publication in ApJ Letters
Submitted: 2018-12-10
A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper.
[35]  oai:arXiv.org:1812.04045  [pdf] - 1811245
The Disk Substructures at High Angular Resolution Project (DSHARP): VII. The Planet-Disk Interactions Interpretation
Comments: 36 pages, 22 figures, ApJ Letters, in press
Submitted: 2018-12-10
The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.
[36]  oai:arXiv.org:1812.04046  [pdf] - 1811246
Disk Substructures at High Angular Resolution Program (DSHARP): VIII. The Rich Ringed Substructures in the AS 209 Disk
Comments: ApJ Letters, in press
Submitted: 2018-12-10
We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content.
[37]  oai:arXiv.org:1811.12036  [pdf] - 1818788
The Flared Gas Structure of the Transitional Disk around Sz 91
Comments: 18 pages, 13 figures; Accepted for publication in The Astrophysical Journal
Submitted: 2018-11-29
We report $0.14"$ resolution observations of the dust continuum at band 7, and the CO(3--2) and HCO$^{+}$(4--3) line emissions toward the transitional disk around Sz 91 with Atacama Large Millimeter/submillimeter Array (ALMA). The dust disk appears to be an axisymmetric ring, peaking a radius of $\sim$95~au from a Gaussian fit. The Gaussian fit widths of the dust ring are 24.6 and 23.7~au for the major and the minor axes, respectively, indicating that the dust ring is not geometrically thin. The gas disk extends out to $\sim$320~au and is also detected in the inner hole of the dust ring. A twin-line pattern is found in the channel maps of CO, which can be interpreted as the emission from the front and rear of the flared gas disk. We perform the radiative transfer calculations using RADMC-3D, to check whether the twin-line pattern can be reproduced under the assumption that the flared gas disk has a power-law form for the column density and $T_\mathrm{gas}=T_\mathrm{dust}$. The thermal Monte Carlo calculation in RADMC-3D shows that the disk temperature has a gradient along the vertical direction beyond the dust ring, as it blocks the stellar radiation, and thus the twin-line pattern can be naturally explained by the flared gas disk in combination with the dust ring. In addition, no significant depletion of the CO molecules in the cold midplane achieves a reasonable agreement with the observed twin-line pattern. This result indicates that the CO emission from the rear surface must be heavily absorbed in the cold midplane.
[38]  oai:arXiv.org:1811.06880  [pdf] - 1799985
Imaging the disc rim and a moving close-in companion candidate in the pre-transitional disc of V1247 Orionis
Comments: 17 pages, 14 figures
Submitted: 2018-11-16
V1247 Orionis harbours a pre-transitional disc with a partially cleared gap. Earlier interferometric and polarimetric observations revealed strong asymmetries both in the gap region and in the outer disc. The presence of a companion was inferred to explain these asymmetric structures and the ongoing disc clearing. Using an extensive set of multi-wavelength and multi-epoch observations we aimed to identify the origin of the previously detected asymmetries. We have observed V1247 Ori at three epochs spanning $\sim678$ days using sparse aperture masking interferometry with Keck/NIRC2 and VLT/NACO. In addition, we search for signs of accretion through VLT/SPHERE-ZIMPOL spectral differential imaging in H$\alpha$ and R-band continuum. Our SMA sub-millimetre interferometry in 880 $\mu$m continuum and in the CO(3-2) line allows us to constrain the orientation and direction of rotation of the outer disc. We find the L'-band emission to be dominated by static features which trace forward-scattered dust emission from the inner edge of the outer disc located to the north-east. In H- and K-band, we see evidence for a companion candidate that moved systematically by 45$^{\circ}$ within the first $\sim$345 days. The separation of the companion candidate is not well constrained, but the observed position angle change is consistent with Keplerian motion of a body located on a 6 au orbit. From the SMA CO moment map, the location of the disc rim, and the detected orbital motion, we deduced the three-dimensional orientation of the disc. We see no indication of accretion in H$\alpha$ and set upper limits for an accreting companion. The measured contrast of the companion candidate in H and K is consistent with an actively accreting protoplanet. Hence, we identify V1247 Ori as a unique laboratory for studying companion-disc interactions and disc clearing.
[39]  oai:arXiv.org:1810.11150  [pdf] - 1774136
Disk Winds and the Evolution of Planet-Forming Disks
Comments: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-25
Disk winds are thought to play an important role in the evolution and dispersal of planet-forming disks. While high-resolution optical and infrared spectroscopy has identified several disk wind diagnostics, wind mass loss rates remain largely unconstrained mostly due to the lack of spatial resolution to measure the extent of the wind emitting region. Here, we show that the ngVLA will have the sensitivity and resolution to detect and spatially resolve the free-free emission from the fully or partially ionized component of disk winds. Hence, ngVLA observations will be critical to estimate mass loss rates and clarify the role of disk winds in the evolution and dispersal of disk mass.
[40]  oai:arXiv.org:1810.06604  [pdf] - 1767802
Science with an ngVLA: Tracing the Water Snowline in Protoplanetary disks with the ngVLA
Comments: 5 pages, 2 figures. To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-15
The water snowline in protoplanetary disks is one of the most pivotal locations for planet formation: it sets a critical boundary of chemical composition in the disk and likely serves as a favorable site of planetesimal growth and planet formation. However, the water snowline at the mid-plane cannot be directly traced by water line observations as the dust emission around the water snowline is highly optically thick at these line frequencies. Alternatively, the water snowline maybe traced through the sharp changes in physical properties (surface density and dust size distribution) at this transition or through the snowlines of other molecules that sublimate co-spatially with water. In this chapter, we discuss the ngVLA's capability of tracing the mid-plane water snowline through mapping sharp transitions in the continuum spectral index caused by changes of dust size distributions and by using the NH3 snowline as a proxy. We found that with the current design of the ngVLA ability, a sharp transition in spectral index across the water snowline can be readily traced in relatively bright disks. A detection of the NH3 snowline is challenging under the current design, but may still be possible if the NH3-to-water abundance in the disk is close to that of the ISM ices. In summary, the proposed ngVLA is the only facility can provide sufficient sensitivity and spatial resolution to trace the evolution of water snowline at the disk mid-plane and to reveal its role in accelerating the planetesimal formation and in building the architecture of planetary systems.
[41]  oai:arXiv.org:1810.06598  [pdf] - 1767801
Science with an ngVLA: Resolved Substructures in Protoplanetary Disks
Comments: To be published in the ASP Monograph Series, Science with a Next-Generation VLA, ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-15
Terrestrial planets and the cores of giant planets are thought to be built by the collisional agglomeration of solids spanning over 20 orders of magnitude in size within a few million years. However, there is tension between this basic picture of planet formation and standard theoretical assumptions associated with the migration of "pebbles" ($\sim$mm/cm-sized particles) in gas-rich disks and the presumably much longer timescales necessary to assemble ($\sim$km-scale) "planetesimals". To confront these potential theoretical discrepancies with observational constraints, the ideal tracer of the solids concentrated in protoplanetary disk substructures is the 30-100 GHz continuum, which strikes the best balance in sensitivity (emission still bright), optical depth (low enough to reliably estimate densities), and angular resolution (high enough to resolve fine-scale features at disk radii as small as 1 au). With its combination of sensitivity, frequency coverage, and angular resolution, the next-generation VLA will be the only facility that has the capabilities to open up this new window into the physics of planetesimal formation.
[42]  oai:arXiv.org:1809.01082  [pdf] - 1789496
Shadows and asymmetries in the T Tauri disk HD 143006: Evidence for a misaligned inner disk
Comments: Accepted for publication in A&A
Submitted: 2018-09-04
While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, reaching an angular resolution of ~0.037" (~6 au). We obtained two datasets, one with a 145 mas diameter coronagraph, and the other without, enabling us to probe the disk structure down to an angular separation of ~0.06" (~10 au). In our observations, the disk of HD 143006 is clearly resolved up to ~0.5" and shows a clear large-scale asymmetry with the eastern side brighter than the western side. We detect a number of additional features, including two gaps and a ring. The ring shows an overbrightness at a position angle (PA) of ~140 deg, extending over a range in position angle of ~60 deg, and two narrow dark regions. The two narrow dark lanes and the overall large-scale asymmetry are indicative of shadowing effects, likely due to a misaligned inner disk. We demonstrate the remarkable resemblance between the scattered light image of HD 143006 and a model prediction of a warped disk due to an inclined binary companion. The warped disk model, based on the hydrodynamic simulations combined with 3D radiative transfer calculations, reproduces all major morphological features. However, it does not account for the observed overbrightness at PA~140 deg. Shadows have been detected in several protoplanetary disks, suggesting that misalignment in disks is not uncommon. However, the origin of the misalignment is not clear. As-yet-undetected stellar or massive planetary companions could be responsible for them, and naturally account for the presence of depleted inner cavities.
[43]  oai:arXiv.org:1808.10682  [pdf] - 1767625
Constraining Gas-Phase Carbon, Oxygen, and Nitrogen in the IM Lup Protoplanetary Disk
Comments: 16 pages, 8 figures, accepted for publication in ApJ
Submitted: 2018-08-31
We present new constraints on gas-phase C, N, and O abundances in the molecular layer of the IM Lup protoplanetary disk. Building on previous physical and chemical modeling of this disk, we use new ALMA observations of C$_2$H to constrain the C/O ratio in the molecular layer to be $\sim0.8$, i.e., higher than the solar value of $\sim0.54$. We use archival ALMA observations of HCN and H$^{13}$CN to show that no depletion of N is required (assuming an interstellar abundance of $7.5\times10^{-5}$ per H). These results suggest that an appreciable fraction of O is sequestered in water ice in large grains settled to the disk mid-plane. Similarly, a fraction of the available C is locked up in less volatile molecules. By contrast, N remains largely unprocessed, likely as N$_2$. This pattern of depletion suggests the presence of true abundance variations in this disk, and not a simple overall depletion of gas mass. If these results hold more generally, then combined CO, C$_2$H, and HCN observations of disks may provide a promising path for constraining gas-phase C/O and N/O during planet-formation. Together, these tracers offer the opportunity to link the volatile compositions of disks to the atmospheres of planets formed from them.
[44]  oai:arXiv.org:1808.10510  [pdf] - 1767623
Scaling Relations Associated with Millimeter Continuum Sizes in Protoplanetary Disks
Comments: Astrophysical Journal, in press
Submitted: 2018-08-30
We present a combined, homogenized analysis of archival Submillimeter Array (SMA) and Atacama Large Millimeter/submillimeter Array (ALMA) observations of the spatially resolved 340 GHz (870 $\mu$m) continuum emission from 105 nearby protoplanetary disks. Building on the previous SMA survey, we infer surface brightness profiles using a simple model of the observed visibilities to derive the luminosities ($L_{\rm mm}$) and effective sizes ($R_{\rm eff}$) of the continuum emission. With this sample, we confirm the shapes, normalizations, and dispersions for the strong correlations between $L_{\rm mm}$, $M_\ast$ (or $L_\ast$), and $\dot{M}_\ast$ found in previous studies. We also verify the continuum size--luminosity relation determined from the SMA survey alone (extending to an order of magnitude lower $L_{\rm mm}$), demonstrating that the amount of emission scales linearly with the emitting surface area. Moreover, we identify new, although weaker, relationships between $R_{\rm eff}$ and the host and accretion properties, such that disks are larger around more massive hosts with higher accretion rates. We explore these inter-related demographic properties with some highly simplified approximations. These multi-dimensional relationships can be explained if the emission is optically thick with a filling factor of $\sim$0.3, or if the emission is optically thin and disks have roughly the same optical depth profile shapes and normalizations independent of host properties. In both scenarios, we require the dust disk sizes to have a slightly sub-linear relationship with the host mass and a non-negligible dispersion ($\sim$0.2 dex at a given $M_\ast$).
[45]  oai:arXiv.org:1807.01325  [pdf] - 1717237
The structure of post-starburst galaxies at $0.5 < z < 2$: evidence for two distinct quenching routes at different epochs
Comments: 21 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2018-07-03
We present an analysis of the structure of post-starburst (PSB) galaxies in the redshift range $0.5 < z < 2$, using a photometrically-selected sample identified in the Ultra Deep Survey (UDS) field. We examine the structure of $\sim80$ of these transient galaxies using radial light $\mu(r)$ profiles obtained from CANDELS $\textit{Hubble Space Telescope}$ near-infrared/optical imaging, and compare to a large sample of $\sim2000$ passive and star-forming galaxies. For each population, we determine their typical structural properties (effective radius $r_{\rm e}$, S\'ersic index $n$) and find significant differences in PSB structure at different epochs. At high redshift ($z > 1$), PSBs are typically massive ($M_* > 10^{10}\rm\,M_{\odot}$), very compact and exhibit high S\'ersic indices, with structures that differ significantly from their star-forming progenitors but are similar to massive passive galaxies. In contrast, at lower redshift ($0.5 < z < 1$), PSBs are generally of low mass ($M_* < 10^{10}\rm\,M_{\odot}$) and exhibit compact but less concentrated profiles (i.e. lower S\'ersic indices), with structures similar to low-mass passive discs. Furthermore, for both epochs we find remarkably consistent PSB structure across the optical/near-infrared wavebands (which largely trace different stellar populations), suggesting that any preceding starburst and/or quenching in PSBs was not strongly centralized. Taken together, these results imply that PSBs at $z > 1$ have been recently quenched during a major disruptive event (e.g. merger or protogalactic collapse) which formed a compact remnant, while at $z < 1$ an alternative less disruptive process is primarily responsible. Our results suggest that high-$z$ PSBs are an intrinsically different population to those at lower redshifts, and indicate different quenching routes are active at different epochs.
[46]  oai:arXiv.org:1806.10553  [pdf] - 1736426
A Subarcsecond ALMA Molecular Line Imaging Survey of the Circumbinary, Protoplanetary Disk Orbiting V4046 Sgr
Comments: 32 pages, 10 figures; accepted 26 June 2018 by The Astrophysical Journal
Submitted: 2018-06-27
We present a suite of ALMA interferometric molecular line and continuum images that elucidate, on linear size scales of $\sim$30--40 AU, the chemical structure of the nearby, evolved, protoplanetary disk orbiting the close binary system V4046 Sgr. The observations were undertaken in the 1.1--1.4 mm wavelength range (ALMA Bands 6 and 7) with antenna configurations involving maximum baselines of several hundred meters, yielding subarcsecond-resolution images in more than a dozen molecular species and isotopologues. Isotopologues of CO and HCN display centrally peaked morphologies of integrated emission line intensity, whereas the line emission from complex nitrile group molecules (HC$_3$N, CH$_3$CN), deuterated molecules (DCN, DCO$^+$), hydrocarbons (as traced by C$_2$H), and potential CO ice line tracers (N$_2$H$^+$, and H$_2$CO) appears as a sequence of sharp and diffuse rings. The dimensions and morphologies of HC$_3$N and CH$_3$CN emission are suggestive of photodesorption of organic ices from the surfaces of dust grains, while the sequence of increasing radius of peak intensity represented by DCN (smallest), DCO$^+$, N$_2$H$^+$, and H$_2$CO (largest) is qualitatively consistent with the expected decline of midplane gas temperature with increasing disk radius. Empirical modeling indicates that the sharp-edged C$_2$H emission ring lies at relatively deep disk layers, leaving open the question of the origin of C$_2$H abundance enhancements in evolved disks. This study of the "molecular anatomy" of V4046 Sgr should serve as motivation for additional subarcsecond ALMA molecular line imaging surveys of nearby, evolved protoplanetary disks aimed at addressing major uncertainties in protoplanetary disk physical and chemical structure and molecular production pathways.
[47]  oai:arXiv.org:1805.12141  [pdf] - 1709508
The Eccentric Cavity, Triple Rings, Two-Armed Spirals, and Double Clumps of the MWC 758 Disk
Comments: accepted by ApJ; the fits file of the ALMA image will be made available by ApJ
Submitted: 2018-05-30
Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.
[48]  oai:arXiv.org:1805.07590  [pdf] - 1694183
Alma Survey Of Circumstellar Disks In The Young Stellar Cluster IC 348
Comments: 19 pages, 12 figures
Submitted: 2018-05-19
We present a 1.3 mm continuum survey of the young (2-3 Myr) stellar cluster IC 348, which lies at a distance of 310 pc, and is dominated by low-mass stars (M$_{\star}$ $\sim$ 0.1-0.6 M$_{\odot}$). We observed 136 Class II sources (disks that are optically thick in the infrared) at 0.8$''$ (200 au) resolution with a 3$\sigma$ sensitivity of $\sim$ 0.45 mJy (M$_{\rm dust}$ $\sim$ 1.3 M$_{\oplus}$). We detect 40 of the targets and construct a mm-continuum luminosity function. We compare the disk mass distribution in IC 348 to those of younger and older regions, taking into account the dependence on stellar mass. We find a clear evolution in disk masses from 1 to 5-10 Myr. The disk masses in IC 348 are significantly lower than those in Taurus (1-3 Myr) and Lupus (1-3 Myr), similar to those of Chamaleon~I, (2-3 Myr) and $\sigma$ Ori (3-5 Myr) and significantly higher than in Upper Scorpius (5$-$10 Myr). About 20 disks in our sample ($\sim$5$\%$ of the cluster members) have estimated masses (dust $+$ gas) $>$1 M$_{\rm Jup}$ and hence might be the precursors of giant planets in the cluster. Some of the most massive disks include transition objects with inner opacity holes based on their infrared SEDs. From a stacking analysis of the 96 non-detections, we find that these disks have a typical dust mass of just $\lesssim$ 0.4 M$_{\oplus}$, even though the vast majority of their infrared SEDs remain optically thick and show little signs of evolution. Such low-mass disks may be the precursors of the small rocky planets found by \emph{Kepler} around M-type stars.
[49]  oai:arXiv.org:1805.06457  [pdf] - 1717138
The Millimeter Continuum Size-Frequency Relationship in the UZ Tau E Disk
Comments: Accepted for publication in The Astrophysical Journal (13 pages, 8 figures)
Submitted: 2018-05-16
We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.
[50]  oai:arXiv.org:1805.03669  [pdf] - 1705313
Protoplanetary Disk Properties in the Orion Nebula Cluster: Initial Results from Deep, High-Resolution ALMA Observations
Comments: Accepted for publication in ApJ
Submitted: 2018-05-09
We present ALMA 850 $\mu$m continuum observations of the Orion Nebula Cluster that provide the highest angular resolution ($\sim 0\rlap{.}''1 \approx 40$ AU) and deepest sensitivity ($\sim 0.1$ mJy) of the region to date. We mosaicked a field containing $\sim 225$ optical or near-IR-identified young stars, $\sim 60$ of which are also optically-identified "proplyds". We detect continuum emission at 850 $\mu$m towards $\sim 80$% of the proplyd sample, and $\sim 50$% of the larger sample of previously-identified cluster members. Detected objects have fluxes of $\sim 0.5$-80 mJy. We remove sub-mm flux due to free-free emission in some objects, leaving a sample of sources detected in dust emission. Under standard assumptions of isothermal, optically thin disks, sub-mm fluxes correspond to dust masses of $\sim 0.5$ to 80 Earth masses. We measure the distribution of disk sizes, and find that disks in this region are particularly compact. Such compact disks are likely to be significantly optically thick. The distributions of sub-mm flux and inferred disk size indicate smaller, lower-flux disks than in lower-density star-forming regions of similar age. Measured disk flux is correlated weakly with stellar mass, contrary to studies in other star forming regions that found steeper correlations. We find a correlation between disk flux and distance from the massive star $\theta^1$ Ori C, suggesting that disk properties in this region are influenced strongly by the rich cluster environment.
[51]  oai:arXiv.org:1804.06269  [pdf] - 1705261
ALMA observations of polarization from dust scattering in the IM Lup protoplanetary disk
Comments: 11 pages, 4 figures, 1 table. Accepted for publication in the Astrophysical Journal. Materials accessible in the online version of the (open-access) ApJ article include the FITS files used to make the ALMA image in Figure 1, and a full, machine-readable version of Table 1
Submitted: 2018-04-17
We present 870 $\mu$m ALMA observations of polarized dust emission toward the Class II protoplanetary disk IM Lup. We find that the orientation of the polarized emission is along the minor axis of the disk, and that the value of the polarization fraction increases steadily toward the center of the disk, reaching a peak value of ~1.1%. All of these characteristics are consistent with models of self-scattering of submillimeter-wave emission from an optically thin inclined disk. The distribution of the polarization position angles across the disk reveals that while the average orientation is along the minor axis, the polarization orientations show a significant spread in angles; this can also be explained by models of pure scattering. We compare the polarization with that of the Class I/II source HL Tau. A comparison of cuts of the polarization fraction across the major and minor axes of both sources reveals that IM Lup has a substantially higher polarization fraction than HL Tau toward the center of the disk. This enhanced polarization fraction could be due a number of factors, including higher optical depth in HL Tau, or scattering by larger dust grains in the more evolved IM Lup disk. However, models yield similar maximum grain sizes for both HL Tau (72 $\mu$m) and IM Lup (61 $\mu$m, this work). This reveals continued tension between grain-size estimates from scattering models and from models of the dust emission spectrum, which find that the bulk of the (unpolarized) emission in disks is most likely due to millimeter (or even centimeter) sized grains.
[52]  oai:arXiv.org:1803.04987  [pdf] - 1670661
Detecting Weak Spectral Lines in Interferometric Data through Matched Filtering
Comments: 20 pages, 8 figures, accepted for publication in The Astronomical Journal
Submitted: 2018-03-13
Modern radio interferometers enable observations of spectral lines with unprecedented spatial resolution and sensitivity. In spite of these technical advances, many lines of interest are still at best weakly detected and therefore necessitate detection and analysis techniques specialized for the low signal-to-noise ratio (SNR) regime. Matched filters can leverage knowledge of the source structure and kinematics to increase sensitivity of spectral line observations. Application of the filter in the native Fourier domain improves SNR while simultaneously avoiding the computational cost and ambiguities associated with imaging, making matched filtering a fast and robust method for weak spectral line detection. We demonstrate how an approximate matched filter can be constructed from a previously observed line or from a model of the source, and we show how this filter can be used to robustly infer a detection significance for weak spectral lines. When applied to ALMA Cycle 2 observations of CH3OH in the protoplanetary disk around TW Hya, the technique yields a ~53% SNR boost over aperture-based spectral extraction methods, and we show that an even higher boost will be achieved for observations at higher spatial resolution. A Python-based open-source implementation of this technique is available under the MIT license at https://github.com/AstroChem/VISIBLE.
[53]  oai:arXiv.org:1803.04467  [pdf] - 1648981
Witnessing Planetary Systems in the Making with the Next Generation Very Large Array
Comments: 6 pages, 3 figures, White Paper submitted to the National Academies of Science, Engineering, and Medicine's Exoplanet Science Strategy Call for Papers
Submitted: 2018-03-12
The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.
[54]  oai:arXiv.org:1803.03842  [pdf] - 1663380
Turbulence in the TW Hya Disk
Comments: 14 pages, 4 figures, accepted to ApJ
Submitted: 2018-03-10
Turbulence is a fundamental parameter in models of grain growth during the early stages of planet formation. As such, observational constraints on its magnitude are crucial. Here we self-consistently analyze ALMA CO(2-1), SMA CO(3-2), and SMA CO(6-5) observations of the disk around TW Hya and find an upper limit on the turbulent broadening of $<$0.08c$_s$ ($\alpha<$0.007 for $\alpha$ defined only within 2-3 pressure scale heights above the midplane), lower than the tentative detection previously found from an analysis of the CO(2-1) data. We examine in detail the challenges of image plane fitting vs directly fitting the visibilities, while also considering the role of the vertical temperature gradient, systematic uncertainty in the amplitude calibration, and assumptions about the CO abundance, as potential sources of the discrepancy in the turbulence measurements. These tests result in variations of the turbulence limit between $<$0.04c$_s$ and $<$0.13c$_s$, consistently lower than the 0.2-0.4c$_s$ found previously. Having ruled out numerous factors, we restrict the source of the discrepancy to our assumed coupling between temperature and density through hydrostatic equilibrium in the presence of a vertical temperature gradient and/or the confinement of CO to a thin molecular layer above the midplane, although further work is needed to quantify the influence of these prescriptions. Assumptions about hydrostatic equilibrium and the CO distribution are physically motivated, and may have a small influence on measuring the kinematics of the gas, but they become important when constraining small effects such as the strength of the turbulence within a protoplanetary disk.
[55]  oai:arXiv.org:1803.00054  [pdf] - 1652524
Resolved Millimeter Observations of the HR 8799 Debris Disk
Comments: 18 pages, 7 figures, accepted by ApJ
Submitted: 2018-02-20
We present 1.3 millimeter observations of the debris disk surrounding the HR 8799 multi-planet system from the Submillimeter Array to complement archival ALMA observations that spatially filtered away the bulk of the emission. The image morphology at $3.8$ arcsecond (150 AU) resolution indicates an optically thin circumstellar belt, which we associate with a population of dust-producing planetesimals within the debris disk. The interferometric visibilities are fit well by an axisymmetric radial power-law model characterized by a broad width, $\Delta R/R\gtrsim 1$. The belt inclination and orientation parameters are consistent with the planet orbital parameters within the mutual uncertainties. The models constrain the radial location of the inner edge of the belt to $R_\text{in}= 104_{-12}^{+8}$ AU. In a simple scenario where the chaotic zone of the outermost planet b truncates the planetesimal distribution, this inner edge location translates into a constraint on the planet~b mass of $M_\text{pl} = 5.8_{-3.1}^{+7.9}$ M$_{\rm Jup}$. This mass estimate is consistent with infrared observations of the planet luminosity and standard hot-start evolutionary models, with the uncertainties allowing for a range of initial conditions. We also present new 9 millimeter observations of the debris disk from the Very Large Array and determine a millimeter spectral index of $2.41\pm0.17$. This value is typical of debris disks and indicates a power-law index of the grain size distribution $q=3.27\pm0.10$, close to predictions for a classical collisional cascade.
[56]  oai:arXiv.org:1801.05812  [pdf] - 1630364
Rings and gaps in the disc around Elias 24 revealed by ALMA
Comments: 17 pages, 11 figures. Accepted for publication in MNRAS
Submitted: 2018-01-17
We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.
[57]  oai:arXiv.org:1801.03948  [pdf] - 1626504
CO and dust properties in the TW Hya disk from high-resolution ALMA observations
Comments: 27 pages, 11 figures, accepted by ApJ; FITS image files available at https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/PXDKBC
Submitted: 2018-01-11
We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the $^{12}$CO $J = 3-2$ line at a spatial resolution of 8 AU (0\farcs14). The CO emission exhibits a bright inner core, a shoulder at $r\approx70$ AU, and a prominent break in slope at $r\approx90$ AU. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a $^{12}$CO column density profile featuring a steep decrease at $r\approx15$ AU and a secondary bump peaking at $r\approx70$ AU. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disk's vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 $\mu$m continuum observations to produce a spectral index map at a spatial resolution of 2 AU. The spectral index rises sharply at the continuum emission gaps at radii of 25, 41, and 47 AU. This behavior suggests that the grains within the gaps are no larger than a few millimeters. Outside the continuum gaps, the low spectral index values of $\alpha\approx 2$ indicate either that grains up to centimeter size are present, or that the bright continuum rings are marginally optically thick at millimeter wavelengths.
[58]  oai:arXiv.org:1801.03945  [pdf] - 1634316
Molecular reconnaissance of the $\beta$ Pictoris gas disk with the SMA: a low HCN/(CO+CO$_2$) outgassing ratio and predictions for future surveys
Comments: Accepted for publication in ApJ. 19 pages, 7 figures
Submitted: 2018-01-11
The exocometary origin of CO gas has been confirmed in several extrasolar Kuiper belts, with CO ice abundances consistent with Solar System comets. We here present a molecular survey of the $\beta$ Pictoris belt with the Submillimeter Array (SMA), reporting upper limits for CN, HCN, HCO$^+$, N$_2$H$^+$ and H$_2$CO, as well as for H$_2$S, CH$_3$OH, SiO and DCN from archival ALMA data. Non-detections can be attributed to rapid molecular photodissociation due to the A-star's strong UV flux. CN is the longest-lasting and most easily detectable molecule after CO in this environment. We update our NLTE excitation model to include UV fluorescence, finding it plays a key role in CO and CN excitation, and use it to turn the SMA CN/CO flux ratio constraint into an upper limit of $<2.5$% on the HCN/(CO+CO$_2$) ratio of outgassing rates. This value is consistent with, but at the low end of, the broad range observed in Solar System comets. If sublimation dominates outgassing, then this low value may be caused by decreased outgassing for the less volatile molecule HCN compared to CO. If instead UV photodesorption or collisional vaporization of unbound grains dominates outgassing, then this low ratio of rates would imply a low ice abundance ratio, which would in turn indicate a variation in cometary cyanide abundances across planetary systems. To conclude, we make predictions for future molecular surveys and show that CN and HCN should be readily detectable with ALMA around $\beta$ Pictoris for Solar-System-like exocometary compositions.
[59]  oai:arXiv.org:1801.03584  [pdf] - 1641414
On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits
Comments: 53 pages, 5 figures, Accepted for publication in ApJS
Submitted: 2018-01-10
We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $\mu$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade.
[60]  oai:arXiv.org:1712.03415  [pdf] - 1663248
Parameterizing the Supernova Engine and its Effects on Remnants and Basic Yields
Comments: 14 pages, 10 figures, submitted to ApJ
Submitted: 2017-12-09
Core-collapse supernova science is now entering an era where engine models are beginning to make both qualitative and quantitative predictions. Although the evidence in support of the convective engine for core-collapse supernova continues to grow, it is difficult to place quantitative constraints on this engine. Some studies have made specific predictions for the remnant distribution from the convective engine, but the results differ between different groups. Here we use a broad parameterization for the supernova engine to understand the differences between distinct studies. With this broader set of models, we place error bars on the remnant mass and basic yields from the uncertainties in the explosive engine. We find that, even with only 3 progenitors and a narrow range of explosion energies, we can produce a wide range of remnant masses and nucleosynthetic yields.
[61]  oai:arXiv.org:1710.03153  [pdf] - 1608494
The Architecture of the GW Ori Young Triple Star System and Its Disk: Dynamical Masses, Mutual Inclinations, and Recurrent Eclipses
Comments: 26 pages, 15 figures, accepted to ApJ
Submitted: 2017-10-09, last modified: 2017-11-20
We present spatially and spectrally resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of gas and dust orbiting the pre-main sequence hierarchical triple star system GW Ori. A forward-modeling of the ${}^{13}$CO and C${}^{18}$O $J$=2-1 transitions permits a measurement of the total stellar mass in this system, $5.29 \pm 0.09\,M_\odot$, and the circum-triple disk inclination, $137.6 \pm 2.0^\circ$. Optical spectra spanning a 35 year period were used to derive new radial velocities and, coupled with a spectroscopic disentangling technique, revealed that the A and B components of GW Ori form a double-lined spectroscopic binary with a $241.50\pm0.05$ day period; a tertiary companion orbits that inner pair with a $4218\pm50$ day period. Combining the results from the ALMA data and the optical spectra with three epochs of astrometry in the literature, we constrain the individual stellar masses in the system ($M_\mathrm{A} \approx 2.7\,M_\odot$, $M_\mathrm{B} \approx 1.7\,M_\odot$, $M_\mathrm{C} \approx 0.9\,M_\odot$) and find strong evidence that at least one (and likely both) stellar orbital planes are misaligned with the disk plane by as much as $45^\circ$. A $V$-band light curve spanning 30 years reveals several new $\sim$30 day eclipse events 0.1-0.7~mag in depth and a 0.2 mag sinusoidal oscillation that is clearly phased with the AB-C orbital period. Taken together, these features suggest that the A-B pair may be partially obscured by material in the inner disk as the pair approaches apoastron in the hierarchical orbit. Lastly, we conclude that stellar evolutionary models are consistent with our measurements of the masses and basic photospheric properties if the GW Ori system is $\sim$1 Myr old.
[62]  oai:arXiv.org:1711.01266  [pdf] - 1602761
An SMA Continuum Survey of Circumstellar Disks in the Serpens Star-Forming Region
Comments: 20 pages, 12 figures, accept for publication in AJ
Submitted: 2017-11-03
We present observations with the Submillimeter Array of the continuum emission at $\lambda = 1.3$ mm from 62 young stars surrounded by a protoplanetary disk in the Serpens star-forming region. The typical angular resolution for the survey in terms of beam size is $3.5^{\prime \prime}\times2.5^{\prime \prime}$ with a median rms noise level of 1.6 mJy beam$^{-1}$. These data are used to infer the dust content in disks around low-mass stars $(0.1$-$2.5\,M_{\odot})$ at a median stellar age of $1$-$3$ Myr. Thirteen sources were detected in the 1.3 mm dust continuum with inferred dust masses of ${\approx} 10$-$260\,M_{\oplus}$ and an upper limit to the median dust mass of $5.1_{-4.3}^{+6.1}\,M_{\oplus}$, derived using survival analysis. Comparing the protoplanetary disk population in Serpens to those of other nearby star-forming regions, we find that the populations of dust disks in Serpens and Taurus, which have a similar age, are statistically indistinguishable. This is potentially surprising since Serpens has a stellar surface density two orders of magnitude in excess of Taurus. Hence, we find no evidence that dust disks in Serpens have been dispersed as a result of more frequent and/or stronger tidal interactions due its elevated stellar density. We also report that the fraction of Serpens disks with $M_{\rm{dust}} \geq 10\,M_{\oplus}$ is less than 20%, which supports the notion that the formation of giant planets is likely inherently rare or has substantially progressed by a few Myrs.
[63]  oai:arXiv.org:1710.11329  [pdf] - 1605051
Modelling the cosmic spectral energy distribution and extragalactic background light over all time
Comments: Accepted by MNRAS
Submitted: 2017-10-31
We present a phenomological model of the Cosmic Spectral Energy Distribution (CSED) and the integrated galactic light (IGL) over all cosmic time. This model, based on an earlier model by Driver et al. (2013), attributes the cosmic star formation history to two processes -- firstly, chaotic clump accretion and major mergers, resulting in the early-time formation of bulges and secondly, cold gas accretion, resulting in late-time disc formation. Under the assumption of a Universal Chabrier initial mass function, we combine the Bruzual & Charlot (2003) stellar libraries, the Charlot & Fall (2000) dust attenuation prescription and template spectra for emission by dust and active galactic nuclei to predict the CSED -- pre- and post-dust attenuation -- and the IGL throughout cosmic time. The phenomological model, as constructed, adopts a number of basic axioms and empirical results and has minimal free parameters. We compare the model output, as well as predictions from the semi-analytic model GALFORM to recent estimates of the CSED out to $z=1$. By construction, our empirical model reproduces the full energy output of the Universe from the ultraviolet to the far-infrared extremely well. We use the model to derive predictions of the stellar and dust mass densities, again finding good agreement. We find that GALFORM predicts the CSED for $z < 0.3$ in good agreement with the observations. This agreement becomes increasingly poor towards $z = 1$, when the model CSED is $\sim$50 per cent fainter. The latter is consistent with the model underpredicting the cosmic star formation history. As a consequence, GALFORM predicts a $\sim$30 per cent fainter IGL.
[64]  oai:arXiv.org:1710.08686  [pdf] - 1670498
A 1.3 mm SMA Survey of 29 Variable Young Stellar Objects
Comments: 23 pages, 17 figures, accepted to A&A
Submitted: 2017-10-24
Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.
[65]  oai:arXiv.org:1710.08426  [pdf] - 1590159
Far-infrared to millimeter data of protoplanetary disks: dust growth in the Taurus, Ophiuchus, and Chamaeleon I star-forming regions
Comments: ApJ, accepted. 38 pages, 11 figures, 16 tables, 3 appendices. The associated data and tables can be found in the Zenodo repository https://doi.org/10.5281/zenodo.889053
Submitted: 2017-10-23
Far-infrared and (sub)millimeter fluxes can be used to study dust in protoplanetary disks, the building blocks of planets. Here, we combine observations from the Herschel Space Observatory with ancillary data of 284 protoplanetary disks in the Taurus, Chamaeleon I, and Ophiuchus star-forming regions, covering from the optical to mm/cm wavelengths. We analyze their spectral indices as a function of wavelength and determine their (sub)millimeter slopes when possible. Most disks display observational evidence of grain growth, in agreement with previous studies. No correlation is found between other tracers of disk evolution and the millimeter spectral indices. A simple disk model is used to fit these sources, and we derive posterior distributions for the optical depth at 1.3 mm and 10 au, the disk temperature at this same radius, and the dust opacity spectral index. We find the fluxes at 70 microns to correlate strongly with disk temperatures at 10 au, as derived from these simple models. We find tentative evidence for spectral indices in Chamaeleon I being steeper than those of disks in Taurus/Ophiuchus, although more millimeter observations are needed to confirm this trend and identify its possible origin. Additionally, we determine the median spectral energy distribution of each region and find them to be similar across the entire wavelength range studied, possibly due to the large scatter in disk properties and morphologies.
[66]  oai:arXiv.org:1710.06628  [pdf] - 1637549
GAMA/G10-COSMOS/3D-HST: The 0<z<5 cosmic star-formation history, stellar- and dust-mass densities
Comments: Accepted for publication in MNRAS
Submitted: 2017-10-18, last modified: 2017-10-19
We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.
[67]  oai:arXiv.org:1708.07287  [pdf] - 1705134
On the Radio Detectability of Circumplanetary Discs
Comments: 15 pages, 10 figures, 2 tables, submitted to MNRAS
Submitted: 2017-08-24
Discs around young planets, so-called circumplanetary discs (CPDs), are essential for planet growth, satellite formation, and planet detection. We study the millimetre and centimetre emission from accreting CPDs by using the simple $\alpha$ disc model. We find that it is easier to detect CPDs at shorter radio wavelengths (e.g. $\lambda\lesssim$ 1 mm). For example, if the system is 140 pc away from us, deep observations (e.g. 5 hours) at ALMA Band 7 (0.87 mm) are sensitive to as small as 0.03 lunar mass of dust in CPDs. If the CPD is around a Jupiter mass planet 20 AU away from the host star and has $\alpha\lesssim 0.001$, ALMA can detect this disc when it accretes faster than $10^{-10} M_{\odot}/yr$. ALMA can also detect the "minimum mass sub-nebulae" disc if such a disc exists around a young planet in YSOs. However, to distinguish the embedded compact CPD from the circumstellar disc material, we should observe circumstellar discs with large gaps/cavities using the highest resolution possible. We also calculate the CPD fluxes at VLA bands, and discuss the possibility of detecting radio emission from jets/winds launched in CPDs. Finally we argue that, if the radial drift of dust particles is considered, the drifting timescale for millimetre dust in CPDs can be extremely short. It only takes 10$^2$-10$^{3}$ years for CPDs to lose millimetre dust. Thus, for CPDs to be detectable at radio wavelengths, mm-sized dust in CPDs needs to be replenished continuously, or the disc has a significant fraction of micron-sized dust or a high gas surface density so that the particle drifting timescale is long, or the radial drift is prevented by other means (e.g. pressure traps).
[68]  oai:arXiv.org:1706.08977  [pdf] - 1585220
A Millimeter Continuum Size-Luminosity Relationship for Protoplanetary Disks
Comments: Accepted for publication in ApJ
Submitted: 2017-06-27
We present a sub-arcsecond resolution survey of the 340 GHz dust continuum emission from 50 nearby protoplanetary disks, based on new and archival observations with the Submillimeter Array. The observed visibility data were modeled with a simple prescription for the radial surface brightness profile. The results were used to extract intuitive, empirical estimates of the emission "size" for each disk, $R_{\rm eff}$, defined as the radius that encircles a fixed fraction of the total continuum luminosity, $L_{\rm mm}$. We find a significant correlation between the sizes and luminosities, such that $R_{\rm eff} \propto L_{\rm mm}^{0.5}$, providing a confirmation and quantitative characterization of a putative trend that was noted previously. This correlation suggests that these disks have roughly the same average surface brightness interior to their given effective radius, ~0.2 Jy arcsec$^{-2}$ (or 8 K in brightness temperature). The same trend remains, but the 0.2dex of dispersion perpendicular to this relation essentially disappears, when we account for the irradiation environment of each disk with a crude approximation of the dust temperatures based on the stellar host luminosities. We consider two (not mutually exclusive) explanations for the origin of this size-luminosity relationship. Simple models of the growth and migration of disk solids can account for the observed trend for a reasonable range of initial conditions, but only on timescales that are much shorter than the nominal ages present in the sample. An alternative scenario invokes optically thick emission concentrated on unresolved scales, with filling factors of a few tens of percent, that are perhaps manifestations of localized particle traps.
[69]  oai:arXiv.org:1706.04504  [pdf] - 1584678
A Three-Dimensional View of Turbulence: Constraints on Turbulent Motions in the HD 163296 Protoplanetary Disk using DCO$^+$
Comments: Accepted to ApJ
Submitted: 2017-06-14
Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $\alpha\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment.
[70]  oai:arXiv.org:1706.03708  [pdf] - 1584555
ALMA Observations of the Young Substellar Binary System 2M1207
Comments: 10 pages, 6 figures, accepted for publication in AJ
Submitted: 2017-06-12
We present ALMA observations of the 2M1207 system, a young binary made of a brown dwarf with a planetary-mass companion at a projected separation of about 40 au. We detect emission from dust continuum at 0.89 mm and from the $J = 3 - 2$ rotational transition of CO from a very compact disk around the young brown dwarf. The small radius found for this brown dwarf disk may be due to truncation from the tidal interaction with the planetary-mass companion. Under the assumption of optically thin dust emission, we estimated a dust mass of 0.1 $M_{\oplus}$ for the 2M1207A disk, and a 3$\sigma$ upper limit of $\sim 1~M_{\rm{Moon}}$ for dust surrounding 2M1207b, which is the tightest upper limit obtained so far for the mass of dust particles surrounding a young planetary-mass companion. We discuss the impact of this and other non-detections of young planetary-mass companions for models of planet formation, which predict the presence of circum-planetary material surrounding these objects.
[71]  oai:arXiv.org:1706.00833  [pdf] - 1584202
Variable H$^{13}$CO$^+$ Emission in the IM Lup Disk: X-ray Driven Time-Dependent Chemistry?
Comments: 7 pages, 4 figures, accepted for publication in ApJL
Submitted: 2017-06-02
We report the first detection of a substantial brightening event in an isotopologue of a key molecular ion, HCO$^+$, within a protoplanetary disk of a T Tauri star. The H$^{13}$CO$^+$ $J=3-2$ rotational transition was observed three times toward IM Lup between July 2014 and May 2015 with the Atacama Large Millimeter Array. The first two observations show similar spectrally integrated line and continuum fluxes, while the third observation shows a doubling in the disk integrated $J=3-2$ line flux compared to the continuum, which does not change between the three epochs. We explore models of an X-ray active star irradiating the disk via stellar flares, and find that the optically thin H$^{13}$CO$^+$ emission variation can potentially be explained via X-ray driven chemistry temporarily enhancing the HCO$^+$ abundance in the upper layers of the disk atmosphere during large or prolonged flaring events. If the HCO$^+$ enhancement is indeed caused by a X-ray flare, future observations should be able to spatially resolve these events and potentially enable us to watch the chemical aftermath of the high-energy stellar radiation propagating across the face of protoplanetary disks, providing a new pathway to explore ionization physics and chemistry, including electron density, in disks.
[72]  oai:arXiv.org:1705.07596  [pdf] - 1583685
Galaxy And Mass Assembly: the evolution of the cosmic spectral energy distribution from z = 1 to z = 0
Comments: 19 pages, 14 figures, accepted by MNRAS
Submitted: 2017-05-22
We present the evolution of the Cosmic Spectral Energy Distribution (CSED) from $z = 1 - 0$. Our CSEDs originate from stacking individual spectral energy distribution fits based on panchromatic photometry from the Galaxy and Mass Assembly (GAMA) and COSMOS datasets in ten redshift intervals with completeness corrections applied. Below $z = 0.45$, we have credible SED fits from 100 nm to 1 mm. Due to the relatively low sensitivity of the far-infrared data, our far-infrared CSEDs contain a mix of predicted and measured fluxes above $z = 0.45$. Our results include appropriate errors to highlight the impact of these corrections. We show that the bolometric energy output of the Universe has declined by a factor of roughly four -- from $5.1 \pm 1.0$ at $z \sim 1$ to $1.3 \pm 0.3 \times 10^{35}~h_{70}$~W~Mpc$^{-3}$ at the current epoch. We show that this decrease is robust to cosmic variance, SED modelling and other various types of error. Our CSEDs are also consistent with an increase in the mean age of stellar populations. We also show that dust attenuation has decreased over the same period, with the photon escape fraction at 150~nm increasing from $16 \pm 3$ at $z \sim 1$ to $24 \pm 5$ per cent at the current epoch, equivalent to a decrease in $A_\mathrm{FUV}$ of 0.4~mag. Our CSEDs account for $68 \pm 12$ and $61 \pm 13$ per cent of the cosmic optical and infrared backgrounds respectively as defined from integrated galaxy counts and are consistent with previous estimates of the cosmic infrared background with redshift.
[73]  oai:arXiv.org:1705.04074  [pdf] - 1583328
Galaxy And Mass Assembly (GAMA): The galaxy stellar mass function to $z=0.1$ from the r-band selected equatorial regions
Comments: 22 pages; 15 figures; Accepted for publication in MNRAS
Submitted: 2017-05-11
We derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties.
[74]  oai:arXiv.org:1704.05133  [pdf] - 1567776
H2CO distribution and formation in the TW Hya disk
Comments: Accepted for publication in ApJ
Submitted: 2017-04-17
H2CO is one of the most readily detected organic molecules in protoplanetary disks. Yet its distribution and dominant formation pathway(s) remain largely unconstrained. To address these issues, we present ALMA observations of two H2CO lines (3_{12}-2_{11} and 5_{15}-4_{14}) at 0".5 (~30 au) spatial resolution toward the disk around the nearby T Tauri star TW Hya. Emission from both lines is spatially resolved, showing a central depression, a peak at 0".4 radius, and a radial decline at larger radii with a bump at ~1", near the millimeter continuum edge. We adopt a physical model for the disk and use toy models to explore the radial and vertical H2CO abundance structure. We find that the observed emission implies the presence of at least two distinct H2CO gas reservoirs: (1) a warm and unresolved inner component (<10 au), and (2) an outer component that extends from ~15 au to beyond the millimeter continuum edge. The outer component is further constrained by the line ratio to arise in a more elevated disk layer at larger radii. The inferred H2CO abundance structure agrees well with disk chemistry models, which predict efficient H2CO gas-phase formation close to the star, and cold H2CO grain surface formation, through H additions to condensed CO, followed by non-thermal desorption in the outer disk. The implied presence of active grain surface chemistry in the TW Hya disk is consistent with the recent detection of CH3OH emission, and suggests that more complex organic molecules are formed in disks, as well.
[75]  oai:arXiv.org:1704.01972  [pdf] - 1567764
Radial Surface Density Profiles of Gas and Dust in the Debris Disk around 49 Ceti
Comments: 20 pages, 8 figures, accepted for publication in ApJ March 31, 2017 (submitted Nov 2016)
Submitted: 2017-04-06
We present ~0.4 resolution images of CO(3-2) and associated continuum emission from the gas-bearing debris disk around the nearby A star 49 Ceti, observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We analyze the ALMA visibilities in tandem with the broad-band spectral energy distribution to measure the radial surface density profiles of dust and gas emission from the system. The dust surface density decreases with radius between ~100 and 310 au, with a marginally significant enhancement of surface density at a radius of ~110 au. The SED requires an inner disk of small grains in addition to the outer disk of larger grains resolved by ALMA. The gas disk exhibits a surface density profile that increases with radius, contrary to most previous spatially resolved observations of circumstellar gas disks. While ~80% of the CO flux is well described by an axisymmetric power-law disk in Keplerian rotation about the central star, residuals at ~20% of the peak flux exhibit a departure from axisymmetry suggestive of spiral arms or a warp in the gas disk. The radial extent of the gas disk (~220 au) is smaller than that of the dust disk (~300 au), consistent with recent observations of other gas-bearing debris disks. While there are so far only three broad debris disks with well characterized radial dust profiles at millimeter wavelengths, 49 Ceti's disk shows a markedly different structure from two radially resolved gas-poor debris disks, implying that the physical processes generating and sculpting the gas and dust are fundamentally different.
[76]  oai:arXiv.org:1704.02006  [pdf] - 1574722
A Multi-Ringed, Modestly-Inclined Protoplanetary Disk around AA Tau
Comments: Accepted for publication in ApJ, 9 pages, 7 figures
Submitted: 2017-04-06
AA Tau is the archetype for a class of stars with a peculiar periodic photometric variability thought to be related to a warped inner disk structure with a nearly edge-on viewing geometry. We present high resolution ($\sim$0.2") ALMA observations of the 0.87 and 1.3~mm dust continuum emission from the disk around AA Tau. These data reveal an evenly spaced three-ringed emission structure, with distinct peaks at 0.34", 0.66", and 0.99", all viewed at a modest inclination of 59.1$^{\circ}\pm$0.3$^{\circ}$ (decidedly not edge-on). In addition to this ringed substructure, we find non-axisymmetric features including a `bridge' of emission that connects opposite sides of the innermost ring. We speculate on the nature of this `bridge' in light of accompanying observations of HCO$^+$ and $^{13}$CO (J=3--2) line emission. The HCO$^+$ emission is bright interior to the innermost dust ring, with a projected velocity field that appears rotated with respect to the resolved disk geometry, indicating the presence of a warp or inward radial flow. We suggest that the continuum bridge and HCO$^+$ line kinematics could originate from gap-crossing accretion streams, which may be responsible for the long-duration dimming of optical light from AA Tau.
[77]  oai:arXiv.org:1703.09227  [pdf] - 1567743
A Multi-Wavelength Analysis of Dust and Gas in the SR 24S Transition Disk
Comments: Accepted for publication in ApJ
Submitted: 2017-03-27
We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the $^{13}$CO and C$^{18}$O (J=2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is $\lesssim 0.12\,M_{\bigoplus}$, which gives a disk mass ratio in dust between the two components of $M_{\mathrm{dust, SR\,24S}}/M_{\mathrm{dust, SR\,24N}}\gtrsim840$. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to mm-sizes is inhibited there and that only warm gas, as seen by ro-vibrational CO emission inside the truncation radii of the binary, is present.
[78]  oai:arXiv.org:1702.05652  [pdf] - 1574641
Disentangling Time Series Spectra with Gaussian Processes: Applications to Radial Velocity Analysis
Comments: submitted to ApJ
Submitted: 2017-02-18
Measurements of radial velocity variations from the spectroscopic monitoring of stars and their companions are essential for a broad swath of astrophysics, providing access to the fundamental physical properties that dictate all phases of stellar evolution and facilitating the quantitative study of planetary systems. The conversion of those measurements into both constraints on the orbital architecture and individual component spectra can be a serious challenge, however, especially for extreme flux ratio systems and observations with relatively low sensitivity. Gaussian processes define sampling distributions of flexible, continuous functions that are well-motivated for modeling stellar spectra, enabling proficient search for companion lines in time-series spectra. We introduce a new technique for spectral disentangling, where the posterior distributions of the orbital parameters and intrinsic, rest-frame stellar spectra are explored simultaneously without needing to invoke cross-correlation templates. To demonstrate its potential, this technique is deployed on red-optical time-series spectra of the mid-M dwarf binary LP661-13. We report orbital parameters with improved precision compared to traditional radial velocity analysis and successfully reconstruct the primary and secondary spectra. We discuss potential applications for other stellar and exoplanet radial velocity techniques and extensions to time-variable spectra. The code used in this analysis is freely available as an open source Python package.
[79]  oai:arXiv.org:1702.04780  [pdf] - 1554031
Polarized disk emission from Herbig Ae/Be stars observed using Gemini Planet Imager: HD 144432, HD 150193, HD 163296, and HD 169142
Comments: Accepted Astrophysical Journal
Submitted: 2017-02-15
In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.
[80]  oai:arXiv.org:1701.01735  [pdf] - 1534242
An ALMA survey of DCN/H$^{13}$CN and DCO$^+$/H$^{13}$CO$^+$ in protoplanetary disks
Comments: 36 pages, 14 figures, updated to match figure order of published version in ApJ
Submitted: 2017-01-06, last modified: 2017-02-10
The deuterium enrichment of molecules is sensitive to their formation environment. Constraining patterns of deuterium chemistry in protoplanetary disks is therefore useful for probing how material is inherited or reprocessed throughout the stages of star and planet formation. We present ALMA observations at $\sim0.6"$ resolution of DCO$^+$, H$^{13}$CO$^+$, DCN, and H$^{13}$CN in the full disks around T Tauri stars AS 209 and IM Lup, the transition disks around T Tauri stars V4046 Sgr and LkCa 15, and the full disks around Herbig Ae stars MWC 480 and HD 163296. We also present ALMA observations of HCN in the IM Lup disk. DCN, DCO$^+$, and H$^{13}$CO$^+$ are detected in all disks, and H$^{13}$CN in all but the IM Lup disk. We find efficient deuterium fractionation for the sample, with estimates of disk-averaged DCO$^+$/HCO$^+$ and DCN/HCN abundance ratios ranging from $\sim0.02-0.06$ and $\sim0.005-0.08$, respectively, which is comparable to values reported for other ISM environments. The relative distributions of DCN and DCO$^+$ vary between disks, suggesting that multiple formation pathways may be needed to explain the diverse emission morphologies. In addition, gaps and rings observed in both H$^{13}$CO$^+$ and DCO$^+$ emission provide new evidence that DCO$^+$ bears a complex relationship with the location of the midplane CO snowline.
[81]  oai:arXiv.org:1611.06229  [pdf] - 1532715
ALMA Measurements of Circumstellar Material in the GQ Lup System
Comments: 11 pages, 4 figures
Submitted: 2016-11-18
We present ALMA observations of the GQ Lup system, a young Sun-like star with a substellar mass companion in a wide-separation orbit. These observations of 870 $\mu$m continuum and CO J=3-2 line emission with beam size $\sim0.3''$ ($\sim45$ AU) resolve the disk of dust and gas surrounding the primary star, GQ Lup A, and provide deep limits on any circumplanetary disk surrounding the companion, GQ Lup b. The circumprimary dust disk is compact with a FWHM of $59\pm12$ AU, while the gas has a larger extent with a characteristic radius of $46.5\pm1.8$ AU. By forward-modeling the velocity field of the circumprimary disk based on the CO emission, we constrain the mass of GQ Lup A to be $M_* = (1.03\pm0.05)*(d/156\text{ pc})$ $M_\odot$, where $d$ is a known distance, and determine that we view the disk at an inclination angle of $60.5^\circ\pm0.5^\circ$ and a position angle of $346^\circ \pm1^\circ$. The $3\sigma$ upper limit on the 870 $\mu$m flux density of any circumplanetary disk associated with GQ Lup b of $<0.15$ mJy implies an upper limit on the dust disk mass of $<0.04$ $M_\oplus$ for standard assumptions about optically thin emission. We discuss proposed mechanisms for the formation of wide-separation substellar companions given the non-detection of circumplanetary disks around GQ Lup b and other similar systems.
[82]  oai:arXiv.org:1610.05139  [pdf] - 1497661
Spiral Density Waves in a Young Protoplanetary Disk
Comments: This is our own version of the manuscript, the definitive version was published in Science (DOI: 10.1126/science.aaf8296) on September 30, 2016. Posted to the arxiv for non-commercial use
Submitted: 2016-10-17
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.
[83]  oai:arXiv.org:1610.00715  [pdf] - 1528262
The Coupled Physical Structure of Gas and Dust in the IM Lup Protoplanetary Disk
Comments: 19 pages, 16 figures, accepted to ApJ
Submitted: 2016-10-03
The spatial distribution of gas and solids in protoplanetary disks determines the composition and formation efficiency of planetary systems. A number of disks show starkly different distributions for the gas and small grains compared to millimeter-centimeter sized dust. We present new Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the dust continuum, CO, $^{13}$CO, and C$^{18}$O in the IM Lup protoplanetary disk, one of the first systems where this dust-gas dichotomy was clearly seen. The $^{12}$CO is detected out to a radius of 970 AU, while the millimeter continuum emission is truncated at just 313 AU. Based upon this data, we have built a comprehensive physical and chemical model for the disk structure, which takes into account the complex, coupled nature of the gas and dust and the interplay between the local and external environment. We constrain the distributions of gas and dust, the gas temperatures, the CO abundances, the CO optical depths, and the incident external radiation field. We find that the reduction/removal of dust from the outer disk exposes this region to higher stellar and external radiation and decreases the rate of freeze-out, allowing CO to remain in the gas out to large radial distances. We estimate a gas-phase CO abundance of $5\%$ of the ISM value and a low external radiation field ($G_0\lesssim4$). The latter is consistent with that expected from the local stellar population. We additionally find tentative evidence for ring-like continuum substructure, suggestions of isotope-selective photodissociation, and a diffuse gas halo.
[84]  oai:arXiv.org:1609.09579  [pdf] - 1513726
G10/COSMOS: 38 band (far-UV to far-IR) panchromatic photometry using LAMBDAR
Comments: Accepted by MNRAS
Submitted: 2016-09-29
We present a consistent total flux catalogue for a $\sim$1 deg$^2$ subset of the COSMOS region (R.A. $\in [149.55\degr, 150.65\degr]$, DEC $\in [1.80\degr, 2.73\degr]$) with near-complete coverage in 38 bands from the far-ultraviolet to the far-infrared. We produce aperture matched photometry for 128,304 objects with i < 24.5 in a manner that is equivalent to the Wright et al. (2016) catalogue from the low-redshift (z < 0.4) Galaxy and Mass Assembly (GAMA) survey. This catalogue is based on publicly available imaging from GALEX, CFHT, Subaru, VISTA, Spitzer and Herschel, contains a robust total flux measurement or upper limit for every object in every waveband and complements our re-reduction of publicly available spectra in the same region. We perform a number of consistency checks, demonstrating that our catalogue is comparable to existing data sets, including the recent COSMOS2015 catalogue (Laigle et al. 2016). We also release an updated Davies et al. (2015) spectroscopic catalogue that folds in new spectroscopic and photometric redshift data sets. The catalogues are available for download at http://cutout.icrar.org/G10/dataRelease.php. Our analysis is optimised for both panchromatic analysis over the full wavelength range and for direct comparison to GAMA, thus permitting measurements of galaxy evolution for 0 < z < 1 while minimising the systematic error resulting from disparate data reduction methods.
[85]  oai:arXiv.org:1609.01800  [pdf] - 1483537
Galaxy And Mass Assembly (GAMA): The absence of stellar mass segregation in galaxy groups and consistent predictions from GALFORM and EAGLE simulations
Comments: 17 pages, 10 figures, Accepted for publication in MNRAS
Submitted: 2016-09-06
We investigate the contentious issue of the presence, or lack thereof, of satellites mass segregation in galaxy groups using the Galaxy And Mass Assembly (GAMA) survey, the GALFORM semi-analytic and the EAGLE cosmological hydrodynamical simulation catalogues of galaxy groups. We select groups with halo mass $12 \leqslant \log(M_{\text{halo}}/h^{-1}M_\odot) <14.5$ and redshift $z \leqslant 0.32$ and probe the radial distribution of stellar mass out to twice the group virial radius. All the samples are carefully constructed to be complete in stellar mass at each redshift range and efforts are made to regularise the analysis for all the data. Our study shows negligible mass segregation in galaxy group environments with absolute gradients of $\lesssim0.08$ dex and also shows a lack of any redshift evolution. Moreover, we find that our results at least for the GAMA data are robust to different halo mass and group centre estimates. Furthermore, the EAGLE data allows us to probe much fainter luminosities ($r$-band magnitude of 22) as well as investigate the three-dimensional spatial distribution with intrinsic halo properties, beyond what the current observational data can offer. In both cases we find that the fainter EAGLE data show a very mild spatial mass segregation at $z \leqslant 0.22$, which is again not apparent at higher redshift. Interestingly, our results are in contrast to some earlier findings using the Sloan Digital Sky Survey. We investigate the source of the disagreement and suggest that subtle differences between the group finding algorithms could be the root cause.
[86]  oai:arXiv.org:1608.05526  [pdf] - 1498409
Galaxy And Mass Assembly (GAMA): the Stellar Mass Budget of Galaxy Spheroids and Disks
Comments: 14 pages; 10 figures; MNRAS accepted
Submitted: 2016-08-19
We build on a recent photometric decomposition analysis of 7506 Galaxy and Mass Assembly (GAMA) survey galaxies to derive stellar mass function fits to individual spheroid and disk component populations down to a lower mass limit of log(M_*/M_sun)= 8. We find that the spheroid/disk mass distributions for individual galaxy morphological types are well described by single Schechter function forms. We derive estimates of the total stellar mass densities in spheroids (rho_spheroid = 1.24+/-0.49 * 10^8 M_sun Mpc^-3 h_0.7) and disks (rho_disk = 1.20+/-0.45 * 10^8 M_sun Mpc^-3 h_0.7), which translates to approximately 50% of the local stellar mass density in spheroids and 48% in disks. The remaining stellar mass is found in the dwarf "little blue spheroid" class, which is not obviously similar in structure to either classical spheroid or disk populations. We also examine the variation of component mass ratios across galaxy mass and group halo mass regimes, finding the transition from spheroid to disk mass dominance occurs near galaxy stellar mass ~10^11 M_sun and group halo mass ~10^12.5 M_sun/h. We further quantify the variation in spheroid-to-total mass ratio with group halo mass for central and satellite populations as well as the radial variation of this ratio within groups.
[87]  oai:arXiv.org:1606.06299  [pdf] - 1426467
GAMA/H-ATLAS: A meta-analysis of SFR indicators - comprehensive measures of the SFR-M* relation and Cosmic Star Formation History at z < 0.4
Comments: 31 pages, 15 Figures
Submitted: 2016-06-20
We present a meta-analysis of star-formation rate (SFR) indicators in the GAMA survey, producing 12 different SFR metrics and determining the SFR-M* relation for each. We compare and contrast published methods to extract the SFR from each indicator, using a well-defined local sample of morphologically-selected spiral galaxies, which excludes sources which potentially have large recent changes to their SFR. The different methods are found to yield SFR-M* relations with inconsistent slopes and normalisations, suggesting differences between calibration methods. The recovered SFR-M* relations also have a large range in scatter which, as SFRs of the targets may be considered constant over the different timescales, suggests differences in the accuracy by which methods correct for attenuation in individual targets. We then recalibrate all SFR indicators to provide new, robust and consistent luminosity-to-SFR calibrations, finding that the most consistent slopes and normalisations of the SFR-M* relations are obtained when recalibrated using the radiation transfer method of Popescu et al. These new calibrations can be used to directly compare SFRs across different observations, epochs and galaxy populations. We then apply our calibrations to the GAMA II equatorial dataset and explore the evolution of star-formation in the local Universe. We determine the evolution of the normalisation to the SFR-M* relation from 0 < z < 0.35 - finding consistent trends with previous estimates at 0.3 < z < 1.2. We then provide the definitive z < 0.35 Cosmic Star Formation History, SFR-M* relation and its evolution over the last 3 billion years.
[88]  oai:arXiv.org:1605.06516  [pdf] - 1415239
Evidence for a CO desorption front in the outer AS 209 disk
Comments: 8 pages, 4 figures, published in ApJL
Submitted: 2016-05-20
Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=2-1$ emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the C$^{18}$O emission consists of a central peak and a ring at a radius of $\sim1''$ (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.
[89]  oai:arXiv.org:1605.01523  [pdf] - 1460876
Extra-galactic background light measurements from the far-UV to the far-IR from deep ground and space-based galaxy counts
Comments: Accepted for publication in The Astrophysical Journal
Submitted: 2016-05-05
We combine wide and deep galaxy number-count data from GAMA, COSMOS/G10, HST ERS, HST UVUDF and various near-, mid- and far- IR datasets from ESO, Spitzer and Herschel. The combined data range from the far-UV (0.15microns) to far-IR (500microns), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated-IGL in all bands. We argue undetected low surface brightness galaxies and intra-cluster/group light is modest, and that our extrapolated-IGL measurements are an accurate representation of the extra-galactic background light. Our data agree with most earlier IGL estimates and with direct measurements in the far-IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC), suggest that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted Zodiacal light model. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the Cosmic Optical Background (COB) and Cosmic Infrared Background (CIB) of: $24^{+4}_{-4}$nWm$^{-2}$sr$^{-1}$ and $26^{+5}_{-5}$nWm$^{-2}$sr$^{-1}$ respectively (48:52%). Over the next decade, upcoming space missions such as Euclid and WFIRST, have the capacity to reduce the COB error to $<1\%$, at which point comparisons to the very high energy data could have the potential to provide a direct detection and measurement of the reionisation field.
[90]  oai:arXiv.org:1604.01923  [pdf] - 1396995
Galaxy And Mass Assembly (GAMA): Accurate Panchromatic Photometry from Optical Priors using LAMBDAR
Comments: 38 pages, 27 figures, 3 tables, 3 appendices, accepted for publication in MNRAS
Submitted: 2016-04-07, last modified: 2016-04-08
We present the Lambda Adaptive Multi-Band Deblending Algorithm in R (LAMBDAR), a novel code for calculating matched aperture photometry across images that are neither pixel- nor PSF-matched, using prior aperture definitions derived from high resolution optical imaging. The development of this program is motivated by the desire for consistent photometry and uncertainties across large ranges of photometric imaging, for use in calculating spectral energy distributions. We describe the program, specifically key features required for robust determination of panchromatic photometry: propagation of apertures to images with arbitrary resolution, local background estimation, aperture normalisation, uncertainty determination and propagation, and object deblending. Using simulated images, we demonstrate that the program is able to recover accurate photometric measurements in both high-resolution, low-confusion, and low-resolution, high-confusion, regimes. We apply the program to the 21-band photometric dataset from the Galaxy And Mass Assembly (GAMA) Panchromatic Data Release (PDR; Driver et al. 2016), which contains imaging spanning the far-UV to the far-IR. We compare photometry derived from LAMBDAR with that presented in Driver et al. (2016), finding broad agreement between the datasets. Nonetheless, we demonstrate that the photometry from LAMBDAR is superior to that from the GAMA PDR, as determined by a reduction in the outlier rate and intrinsic scatter of colours in the LAMBDAR dataset. We similarly find a decrease in the outlier rate of stellar masses and star formation rates using LAMBDAR photometry. Finally, we note an exceptional increase in the number of UV and mid-IR sources able to be constrained, which is accompanied by a significant increase in the mid-IR colour-colour parameter-space able to be explored.
[91]  oai:arXiv.org:1603.09352  [pdf] - 1389005
Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk
Comments: published in Astrophysical Journal Letters; image files are available at http://www.cfa.harvard.edu/~sandrews/data/twhya/
Submitted: 2016-03-30
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.
[92]  oai:arXiv.org:1603.05644  [pdf] - 1415131
Constraints on Planetesimal Collision Models in Debris Disks
Comments: 34 pages, 5 figures
Submitted: 2016-03-17
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$ $\mu$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.
[93]  oai:arXiv.org:1601.07542  [pdf] - 1382118
Dust properties across the CO snowline in the HD 163296 disk from ALMA and VLA observations
Comments: 12 pages, 17 figures. Accepted for publication in A\&A
Submitted: 2016-01-27
To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.
[94]  oai:arXiv.org:1601.03806  [pdf] - 1370345
A Disk-based Dynamical Constraint on the Mass of the Young Binary DQ Tau
Comments: Accepted to The Astrophysical Journal. 9 pages, 5 figures
Submitted: 2016-01-14
We present new ALMA observations of CO $J$=2$-$1 line emission from the DQ Tau circumbinary disk. These data are used to tomographically reconstruct the Keplerian disk velocity field in a forward-modeling inference framework, and thereby provide a dynamical constraint on the mass of the DQ Tau binary of $M_\ast = 1.27_{-0.27}^{+0.46} \,M_\odot$. Those results are compared with an updated and improved orbital solution for this double-lined system based on long-term monitoring of its stellar radial velocities. Both of these independent dynamical constraints on the binary mass are in excellent agreement: taken together, they demonstrate that the DQ Tau system mass is $1.21\pm0.26\,M_\odot$ and that the disk and binary orbital planes are aligned within $3^\circ$ (at 3$\sigma$ confidence). The predictions of various theoretical models for pre-main sequence stellar evolution are also consistent with these masses, although more detailed comparisons are difficult due to lingering uncertainties in the photospheric properties of the individual components. DQ Tau is the third nearly equal-mass double-lined spectroscopic binary with a circumbinary disk that has been dynamically "weighed" with these two independent techniques: all show consistent results, validating the overall accuracy of the disk-based approach and demonstrating that it can be robustly applied to large samples of young, single stars as ALMA ramps up to operations at full capacity.
[95]  oai:arXiv.org:1601.02642  [pdf] - 1359200
Resolved Gas Interior to the Dust Rings of the HD 141569 Disk
Comments: Accepted to ApJ
Submitted: 2016-01-11
The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks.
[96]  oai:arXiv.org:1512.05679  [pdf] - 1378751
A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure
Comments: 19 pages, 15 figures, accepted for publication in A&A
Submitted: 2015-12-17
Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.
[97]  oai:arXiv.org:1512.05553  [pdf] - 1351651
Steepening of the 820 micron continuum surface-brightness profile signals dust evolution in TW Hya's disk
Comments: (8 pages, 4 figures, accepted by Astronomy & Astrophysics)
Submitted: 2015-12-17
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.
[98]  oai:arXiv.org:1511.07149  [pdf] - 1327590
Resolved gas cavities in transitional disks inferred from CO isotopologues with ALMA
Comments: Accepted by A&A; version after language-editing
Submitted: 2015-11-23, last modified: 2015-12-01
Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.
[99]  oai:arXiv.org:1508.02076  [pdf] - 1323446
Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget
Comments: 31 pages and 30 figures. Accepted in MNRAS. High-resolution copy available from our release site: http://gama-psi.icrar.org/ or directly via http://www.simondriver.org/mwavev05.pdf
Submitted: 2015-08-09, last modified: 2015-11-20
We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/
[100]  oai:arXiv.org:1510.08890  [pdf] - 1339006
Resolved Millimeter-Wavelength Observations of Debris Disks around Solar-Type Stars
Comments: 17 pages, 5 figures, 7 tables
Submitted: 2015-10-29
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.
[101]  oai:arXiv.org:1510.05660  [pdf] - 1304218
Dust Evolution Can Produce Scattered Light Gaps in Protoplanetary Disks
Comments: 6 pages, 4 figures, accepted for publication in ApJL, fixed typos
Submitted: 2015-10-19, last modified: 2015-10-21
Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps." The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths.
[102]  oai:arXiv.org:1509.03040  [pdf] - 1312007
Testing particle trapping in transition disks with ALMA
Comments: Minor changes after language edition. Accepted for publication in A&A (abstract slightly shortened for arXiv)
Submitted: 2015-09-10, last modified: 2015-10-13
We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336GHz of two transition disks, SR21 and HD135344B. In combination with previous ALMA observations from Cycle 0 at 689GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($\alpha_{\rm{mm}}$). The observations of SR21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR21 with values of $\alpha_{\rm{mm}}{\sim}3.8-4.2$ in the inner region ($r<35$ AU) and $\alpha_{\rm{mm}}{\sim}2.6-3.0$ outside. An axisymmetric ring (which we call the ring model) or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex model), is assumed for fitting the disk morphology. For SR21, the ring model better fits the emission at 336GHz, conversely the vortex model better fits the 689GHz emission. For HD135344B, neither a significant shift in the null of the visibilities nor radial variations of $\alpha_{\rm{mm}}$ are detected. Furthermore, for HD135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $\alpha_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD135344B, although it cannot be excluded either.
[103]  oai:arXiv.org:1510.01375  [pdf] - 1308349
Weak Turbulence in the HD 163296 Protoplanetary Disk Revealed by ALMA CO Observations
Comments: Accepted to ApJ
Submitted: 2015-10-05
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.
[104]  oai:arXiv.org:1510.00968  [pdf] - 1308341
Chemical Imaging of the CO Snow Line in the HD 163296 Disk
Comments: Accepted for publication in ApJ
Submitted: 2015-10-04
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.
[105]  oai:arXiv.org:1510.00739  [pdf] - 1351509
Rings of C2H in the Molecular Disks Orbiting TW Hya and V4046 Sgr
Comments: 2 pages, 2 figures, to appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314 (Cambridge University Press), J.H. Kastner, B. Stelzer, S.A. Metchev, eds
Submitted: 2015-10-02
We have used the Submillimeter Array to image, at ~1" resolution, C2H(3-2) emission from the molecule-rich circumstellar disks orbiting the nearby, classical T Tauri star systems TW Hya and V4046 Sgr. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology within each disk, the inner hole radius of the C2H ring within the V4046 Sgr disk (~70 AU) is somewhat larger than than of its counterpart within the TW Hya disk (~45 AU). We suggest that, in each case, the C2H emission likely traces irradiation of the tenuous surface layers of the outer disks by high-energy photons from the central stars.
[106]  oai:arXiv.org:1509.07520  [pdf] - 1304167
Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25
Comments: 18 pages, 15 figures, accepted for publication at ApJ
Submitted: 2015-09-24
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 {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 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, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} 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.
[107]  oai:arXiv.org:1412.5177  [pdf] - 1296100
Constructing A Flexible Likelihood Function For Spectroscopic Inference
Comments: Accepted to ApJ. Incorporated referees' comments. New figures 1, 8, 10, 12, and 14. Supplemental website: http://iancze.github.io/Starfish/
Submitted: 2014-12-16, last modified: 2015-09-15
We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically address the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line "outliers." By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high resolution $V$-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate resolution $K$-band spectrum of Gliese 51, an M5 field dwarf.
[108]  oai:arXiv.org:1508.01003  [pdf] - 1275506
A concentration of centimeter-sized grains in the Oph IRS 48 dust trap
Comments: Language editing and addition reference ALMA data
Submitted: 2015-08-05, last modified: 2015-09-11
Azimuthally asymmetric dust distributions observed with ALMA in transition disks have been interpreted as dust traps. We present VLA Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 GHz or 0.45 mm) observations at 0.2" resolution of the Oph IRS 48 disk, which suggest that larger particles could be more azimuthally concentrated than smaller dust grains, assuming an axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 $\pm0.9$ times as wide as the centimeter emission, marginally consistent with the particle trapping mechanism under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized gas emission from the star. Both the morphology and the spectral index variations are consistent with an increase of large particles in the center of the trap, but uncertainties remain due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet formation models to allow dust particles to grow beyond millimeter sizes in the outer regions of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the dust trapping mechanism for centimeter-sized grains, although additional data is required for definitive confirmation.
[109]  oai:arXiv.org:1509.02543  [pdf] - 1300297
The Early ALMA View of the FU Ori Outburst System
Comments: ApJ Accepted
Submitted: 2015-09-08
We have obtained ALMA Band 7 observations of the FU Ori outburst system at 0.6"x0.5" resolution to measure the link between the inner disk instability and the outer disk through sub-mm continuum and molecular line observations. Our observations detect continuum emission which can be well modeled by two unresolved sources located at the position of each binary component. The interferometric observations recover the entire flux reported in previous single-dish studies, ruling out the presence of a large envelope. Assuming that the dust is optically thin, we derive disk dust masses of $2\times 10^{-4}$M$_{\odot}$ and $8\times 10^{-5}$M$_{\odot}$, for the north and south components respectively. We place limits on the disks' radii of $r<$45 AU. We report the detection of molecular emission from $^{12}$CO(3-2), HCO$^{+}$(4-3) and from HCN(4-3). The $^{12}$CO appears widespread across the two binary components, and is slightly more extended than the continuum emission. The denser gas tracer HCO$^{+}$ peaks close to the position of the southern binary component, while HCN appears peaked at the position of the northern component. This suggests that the southern binary component is embedded in denser molecular material, consistent with previous studies that indicate a heavily reddened object. At this angular resolution any interaction between the two unresolved disk components cannot be disentangled. Higher resolution images are vital to understanding the process of star formation via rapid accretion FU Ori-type episodes.
[110]  oai:arXiv.org:1508.07296  [pdf] - 1281028
Double DCO+ rings reveal CO ice desorption in the outer disk around IM Lup
Comments: Accepted for publication in ApJ
Submitted: 2015-08-28
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.
[111]  oai:arXiv.org:1508.07004  [pdf] - 1267976
The Distribution and Chemistry of H$_2$CO in the DM Tau Protoplanetary Disk
Comments: 6 pages, 4 figures, 1 table
Submitted: 2015-08-27
H$_2$CO ice on dust grains is an important precursor of complex organic molecules (COMs). H$_2$CO gas can be readily observed in protoplanetary disks and may be used to trace COM chemistry. However, its utility as a COM probe is currently limited by a lack of constraints on the relative contributions of two different formation pathways: on icy grain-surfaces and in the gas-phase. We use archival ALMA observations of the resolved distribution of H$_2$CO emission in the disk around the young low-mass star DM Tau to assess the relative importance of these formation routes. The observed H$_2$CO emission has a centrally peaked and radially broad brightness profile (extending out to 500 AU). We compare these observations with disk chemistry models with and without grain-surface formation reactions, and find that both gas and grain-surface chemistry are necessary to explain the spatial distribution of the emission. Gas-phase H$_2$CO production is responsible for the observed central peak, while grain-surface chemistry is required to reproduce the emission exterior to the CO snowline (where H$_2$CO mainly forms through the hydrogenation of CO ice before being non-thermally desorbed). These observations demonstrate that both gas and grain-surface pathways contribute to the observed H$_2$CO in disks, and that their relative contributions depend strongly on distance from the host star.
[112]  oai:arXiv.org:1507.05574  [pdf] - 1277111
A Combined Spitzer and Herschel Infrared Study of Gas and Dust in the Circumbinary Disk Orbiting V4046 Sgr
Comments: 33 pages, 9 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2015-07-20
We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.
[113]  oai:arXiv.org:1507.04758  [pdf] - 1370188
Observations of Solids in Protoplanetary Disks
Comments: invited pedagogical review for PASP (appearing in October issue); 29 pages, 20 figures
Submitted: 2015-07-16
This review addresses the state of research that employs astronomical (remote sensing) observations of solids ("dust") in young circumstellar disks to learn about planet formation. The intention is for it to serve as an accessible, introductory, pedagogical resource for junior scientists interested in the subject. After some historical background and a basic observational primer, the focus is shifted to the three fundamental topics that broadly define the field: (1) demographics -- the relationships between disk properties and the characteristics of their environments and hosts; (2) structure -- the spatial distribution of disk material and its associated physical conditions and composition; and (3) evolution -- the signposts of key changes in disk properties, including the growth and migration of solids and the impact of dynamical interactions with young planetary systems. Based on the state of the art results in these areas, suggestions are made for potentially fruitful lines of work in the near future.
[114]  oai:arXiv.org:1507.01642  [pdf] - 1264106
The Epsilon Eridani System Resolved by Millimeter Interferometry
Comments: 27 pages, 6 figures
Submitted: 2015-07-06
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 $\Delta 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.
[115]  oai:arXiv.org:1506.08222  [pdf] - 1259091
Galaxy And Mass Assembly (GAMA): end of survey report and data release 2
Comments: Accepted for publication in MMRAS, 40 pages, 33 figures
Submitted: 2015-06-26
The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component S\'ersic fits, stellar masses, H$\alpha$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.
[116]  oai:arXiv.org:1506.03391  [pdf] - 1258994
ALMA Observations of the Largest Proto-Planetary Disk in the Orion Nebula, 114-426: A CO Silhouette
Comments: 20 pages, 4 figures
Submitted: 2015-06-10
We present ALMA observations of the largest protoplanetary disk in the Orion Nebula, 114-426. Detectable 345 GHz (856 micron) dust continuum is produced only in the 350 AU central region of the ~1000 AU diameter silhouette seen against the bright H-alpha background in HST images. Assuming optically thin dust emission at 345 GHz, a gas-to-dust ratio of 100, and a grain temperature of 20 K, the disk gas-mass is estimated to be 3.1 +/- 0.6 Jupiter masses. If most solids and ices have have been incorporated into large grains, however, this value is a lower limit. The disk is not detected in dense-gas tracers such as HCO+ J=4-3, HCN J=4-3, or CS =7-6. These results may indicate that the 114-426 disk is evolved and depleted in some light organic compounds found in molecular clouds. The CO J=3-2 line is seen in absorption against the bright 50 to 80 K background of the Orion A molecular cloud over the full spatial extent and a little beyond the dust continuum emission. The CO absorption reaches a depth of 27 K below the background CO emission at VLSR ~6.7 km/s about 0.52 arcseconds (210 AU) northeast and 12 K below the background CO emission at VLSR ~ 9.7 km/s about 0.34 arcseconds (140 AU) southwest of the suspected location of the central star, implying that the embedded star has a mass less than 1 Solar mass .
[117]  oai:arXiv.org:1506.01007  [pdf] - 1245833
The Transitional Disk around IRAS 04125+2902
Comments: accepted to ApJ
Submitted: 2015-06-02
Resolved submillimeter imaging of transitional disks is increasingly revealing the complexity of disk structure. Here we present the first high-resolution submillimeter image of a recently identified transitional disk around IRAS 04125+2902 in the Taurus star-forming region. We measure an inner disk hole of ~20 AU around IRAS 04125+2902 by simultaneously modeling new 880 micron Submillimeter Array (SMA) data along with an existing spectral energy distribution supplemented by new Discovery Channel Telescope (DCT) photometry. We also constrain the outer radius of the dust disk in IRAS~04125+2902 to ~50-60 AU. Such a small dust disk could be attributed to initial formation conditions, outward truncation by an unseen companion, or dust evolution in the disk. Notably, the dust distribution of IRAS 04125+2902 resembles a narrow ring (delta R ~ 35 AU) composed of large dust grains at the location of the disk wall. Such narrow dust rings are also seen in other transitional disks and may be evidence of dust trapping in pressure bumps, possibly produced by planetary companions. More sensitive submillimeter observations of the gas are necessary to further probe the physical mechanisms at work in shaping the spatial distribution of large dust in this disk. Interestingly, the IRAS 04125+2902 disk is significantly fainter than other transitional disks that have been resolved at submillimeter wavelengths, hinting that more objects with large disk holes may exist at the faint end of the submillimeter luminosity distribution that await detection with more sensitive imaging telescopes.
[118]  oai:arXiv.org:1505.06347  [pdf] - 1095730
The cometary composition of a protoplanetary disk as revealed by complex cyanides
Comments: Definitive version of the manuscript is published in Nature, 520, 7546, 198, 2015. This is the author's version
Submitted: 2015-05-23
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.
[119]  oai:arXiv.org:1505.05518  [pdf] - 1045068
Galaxy And Mass Assembly (GAMA): Trends in galaxy colours, morphology, and stellar populations with large scale structure, group, and pair environments
Comments: 22 pages, 18 figures, accepted for publication in MNRAS on May 20, 2015
Submitted: 2015-05-20
We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of S\'{e}rsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.
[120]  oai:arXiv.org:1505.01850  [pdf] - 991725
A Disk-based Dynamical Mass Estimate for the Young Binary AK Sco
Comments: Accepted to The Astrophysical Journal
Submitted: 2015-05-07
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{\deg}, 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\,\sigma$), 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.
[121]  oai:arXiv.org:1505.01483  [pdf] - 1095717
An ALMA Constraint on the GSC 6214-210 B Circum-Substellar Accretion Disk Mass
Comments: Accepted to ApJL
Submitted: 2015-05-06
We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of GSC 6214-210 A and B, a solar-mass member of the 5-10 Myr Upper Scorpius association with a 15 $\pm$ 2 Mjup companion orbiting at $\approx$330 AU (2.2"). Previous photometry and spectroscopy spanning 0.3-5 $\mu$m revealed optical and thermal excess as well as strong H$\alpha$ and Pa~$\beta$ emission originating from a circum-substellar accretion disk around GSC 6214-210 B, making it the lowest mass companion with unambiguous evidence of a subdisk. Despite ALMA's unprecedented sensitivity and angular resolution, neither component was detected in our 880 $\mu$m (341 GHz) continuum observations down to a 3-$\sigma$ limit of 0.22 mJy/beam. The corresponding constraints on the dust mass and total mass are <0.15 Mearth and <0.05 Mjup, respectively, or <0.003% and <0.3% of the mass of GSC 6214-210 B itself assuming a 100:1 gas-to-dust ratio and characteristic dust temperature of 10-20 K. If the host star possesses a putative circum-stellar disk then at most it is a meager 0.0015% of the primary mass, implying that giant planet formation has certainly ceased in this system. Considering these limits and its current accretion rate, GSC 6214-210 B appears to be at the end stages of assembly and is not expected to gain any appreciable mass over the next few Myr.
[122]  oai:arXiv.org:1504.06040  [pdf] - 1292175
A SCUBA-2 850 micron Survey of Protoplanetary Discs in the IC 348 Cluster
Comments: Submitted to MNRAS
Submitted: 2015-04-23, last modified: 2015-04-25
We present 850 micron observations of the 2-3 Myr cluster IC 348 in the Perseus molecular cloud using the SCUBA-2 camera on the James Clerk Maxwell Telescope. Our SCUBA-2 map has a diameter of 30 arcmin and contains ~370 cluster members, including ~200 objects with IR excesses. We detect a total of 13 discs. Assuming standard dust properties and a gas to dust mass ratio of 100, we derive disc masses ranging from 1.5 to 16 M_JUP . We also detect 8 Class 0/I protostars. We find that the most massive discs (M_Disc > 3 M_JUP ; 850 micron flux > 10 mJy) in IC 348 tend to be transition objects according to the characteristic "dip" in their infrared Spectral Energy Distributions (SEDs). This trend is also seen in other regions. We speculate that this could be an initial conditions effect (e.g., more massive discs tend to form giant planets that result in transition disc SEDs) and/or a disc evolution effect (the formation of one or more massive planets results in both a transition disc SED and a reduction of the accretion rate, increasing the lifetime of the outer disc). A stacking analysis of the discs that remain undetected in our SCUBA-2 observations suggests that their median 850 micron flux should be ~1 mJy, corresponding to a disc mass ~0.3 M_JUP (gas plus dust) or ~1 M_Earth of dust. While the available data are not deep enough to allow a meaningful comparison of the disc luminosity functions between IC 348 and other young stellar clusters, our results imply that disc masses exceeding the Minimum Mass Solar Nebula are very rare (~1%) at the age of IC 348, especially around very low-mass stars.
[123]  oai:arXiv.org:1504.05980  [pdf] - 1182983
A Ring of C2H in the Molecular Disk Orbiting TW Hya
Comments: 17 pages, 4 figures; to appear in The Astrophysical Journal
Submitted: 2015-04-22
We have used the Submillimeter Array to image, at ~1.5" resolution, C2H (3-2) emission from the circumstellar disk orbiting the nearby (D = 54 pc), ~8 Myr-old, ~0.8 Msun classical T Tauri star TW Hya. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology. Based on a model in which the C2H column density follows a truncated radial power-law distribution, we find that the inner edge of the ring lies at ~45 AU, and that the ring extends to at least ~120 AU. Comparison with previous (single-dish) observations of C2H (4-3) emission indicates that the C2H molecules are subthermally excited and, hence, that the emission arises from the relatively warm, tenuous upper atmosphere of the disk. We propose that the C2H emission most likely traces particularly efficient photo-destruction of small grains and/or photodesorption and photodissociation of hydrocarbons derived from grain ice mantles in the surface layers of the outer disk. The presence of a C2H ring in the TW Hya disk hence likely serves as a marker of dust grain processing and radial and vertical grain size segregation within the disk.
[124]  oai:arXiv.org:1504.05198  [pdf] - 982181
Demographics of Transition Discs in Ophiuchus and Taurus
Comments: 9 pages, 2 figures, accepted for publication in MNRAS
Submitted: 2015-04-20
Transition disc systems are young stars that appear to be on the verge of dispersing their protoplanetary discs. We explore the nature of these systems by comparing the stellar accretion rates and disc masses of transition discs and normal T Tauri stars in Taurus and Ophiuchus. After controlling for the known dependencies of stellar accretion rate and disc mass and on age, stellar accretion rate on stellar mass, and disc mass on the presence of stellar or sub-stellar companions, we find that the normal T Tauri stars show a trend of stellar accretion rate increasing with disc mass. The transition discs tend to have higher average disc masses than normal T Tauri stars as well as lower accretion rates than normal T Tauri stars of the same disc mass. These results are most consistent with the interpretation that the transition discs have formed objects massive enough to alter the accretion flow, i.e., single or multiple giant planets. Several Ophiuchus T Tauri stars that are not known transition disc systems also have very low accretion rates for their disc masses. We speculate on the possible nature of these sources.
[125]  oai:arXiv.org:1503.06192  [pdf] - 1232084
Scattered Light from Dust in the Cavity of the V4046 Sgr Transition Disk
Comments: 13 pages, 4 figures, Accepted to ApJ Letters
Submitted: 2015-03-20
We report the presence of scattered light from dust grains located in the giant planet formation region of the circumbinary disk orbiting the ~20-Myr-old close (~0.045 AU separation) binary system V4046 Sgr AB based on observations with the new Gemini Planet Imager (GPI) instrument. These GPI images probe to within ~7 AU of the central binary with linear spatial resolution of ~3 AU, and are thereby capable of revealing dust disk structure within a region corresponding to the giant planets in our solar system. The GPI imaging reveals a relatively narrow (FWHM ~10 AU) ring of polarized near-infrared flux whose brightness peaks at ~14 AU. This ~14 AU radius ring is surrounded by a fainter outer halo of scattered light extending to ~45 AU, which coincides with previously detected mm-wave thermal dust emission. The presence of small grains that efficiently scatter starlight well inside the mm-wavelength disk cavity supports current models of planet formation that suggest planet-disk interactions can generate pressure traps that impose strong radial variations in the particle size distribution throughout the disk.
[126]  oai:arXiv.org:1501.06512  [pdf] - 1224183
Protoplanetary Disk Masses in the Young NGC 2024 Cluster
Comments: ApJ, in press
Submitted: 2015-01-26
We present the results from a Submillimeter Array survey of the 887 micron continuum emission from the protoplanetary disks around 95 young stars in the young cluster NGC 2024. Emission was detected from 22 infrared sources, with flux densities from ~5 to 330 mJy; upper limits (at 3sigma) for the other 73 sources range from 3 to 24 mJy. For standard assumptions, the corresponding disk masses range from ~0.003 to 0.2Msolar, with upper limits at 0.002--0.01Msolar. The NGC 2024 sample has a slightly more populated tail at the high end of its disk mass distribution compared to other clusters, but without more information on the nature of the sample hosts it remains unclear if this difference is statistically significant or a superficial selection effect. Unlike in the Orion Trapezium, there is no evidence for a disk mass dependence on the (projected) separation from the massive star IRS2b in the NGC 2024 cluster. We suggest that this is due to either the cluster youth or a comparatively weaker photoionizing radiation field.
[127]  oai:arXiv.org:1501.05962  [pdf] - 1224158
Resolved Millimeter Emission from the HD 15115 Debris Disk
Comments: 19 pages, 3 figures. Accepted for publication in ApJ
Submitted: 2015-01-23
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.
[128]  oai:arXiv.org:1412.2175  [pdf] - 1223519
An ALMA Disk Mass for the Candidate Protoplanetary Companion to FW Tau
Comments: 5 pages, 2 figures; accepted to ApJL
Submitted: 2014-12-05
We present ALMA observations of the FW Tau system, a close binary pair of M5 stars with a wide-orbit (300 AU projected separation) substellar companion. The companion is extremely faint and red in the optical and near-infrared, but boasts a weak far-infrared excess and optical/near-infrared emission lines indicative of a primordial accretion disk of gas and dust. The component-resolved 1.3 mm continuum emission is found to be associated only with the companion, with a flux (1.78 +/- 0.03 mJy) that indicates a dust mass of 1-2 M_Earth. While this mass reservoir is insufficient to form a giant planet, it is more than sufficient to produce an analog of the Kepler-42 exoplanetary system or the Galilean satellites. The mass and geometry of the disk-bearing FW Tau companion remains unclear. Near-infrared spectroscopy shows deep water bands that indicate a spectral type later than M5, but substantial veiling prevents a more accurate determination of the effective temperature (and hence mass). Both a disk-bearing "planetary-mass" companion seen in direct light or a brown dwarf tertiary viewed in light scattered by an edge-on disk or envelope remain possibilities.
[129]  oai:arXiv.org:1411.6355  [pdf] - 1451820
Galaxy And Mass Assembly (GAMA): mass-size relations of z$<$0.1 galaxies subdivided by S\'ersic index, colour and morphology
Comments:
Submitted: 2014-11-24
We use data from the Galaxy And Mass Assembly (GAMA) survey in the redshift range 0.01$<$z$<$0.1 (8399 galaxies in $g$ to $K_s$ bands) to derive the stellar mass $-$ half-light radius relations for various divisions of 'early' and 'late'-type samples. We find the choice of division between early and late (i.e., colour, shape, morphology) is not particularly critical, however, the adopted mass limits and sample selections (i.e., the careful rejection of outliers and use of robust fitting methods) are important. In particular we note that for samples extending to low stellar mass limits ($<10^{10}\mathcal{M_{\odot}}$) the S\'ersic index bimodality, evident for high mass systems, becomes less distinct and no-longer acts as a reliable separator of early- and late-type systems. The final set of stellar mass $-$ half-light radius relations are reported for a variety of galaxy population subsets in 10 bands ($ugrizZYJHK_s$) and are intended to provide a comprehensive low-z benchmark for the many ongoing high-z studies. Exploring the variation of the stellar mass $-$ half-light radius relations with wavelength we confirm earlier findings that galaxies appear more compact at longer wavelengths albeit at a smaller level than previously noted: at $10^{10}\mathcal{M_{\odot}}$ both spiral systems and ellipticals show a decrease in size of 13% from $g$ to $K_s$ (which is near linear in log wavelength). Finally we note that the sizes used in this work are derived from 2D S\'ersic light profile fitting (using GALFIT3), i.e., elliptical semi-major half light radii, improving on earlier low-z benchmarks based on circular apertures.
[130]  oai:arXiv.org:1410.3570  [pdf] - 1222559
ALMA observations of a misaligned binary protoplanetary disk system in Orion
Comments: 6 pages, 5 figures, accepted for publication in ApJ
Submitted: 2014-10-14
We present ALMA observations of a wide binary system in Orion, with projected separation 440 AU, in which we detect submillimeter emission from the protoplanetary disks around each star. Both disks appear moderately massive and have strong line emission in CO 3-2, HCO+ 4-3, and HCN 3-2. In addition, CS 7-6 is detected in one disk. The line-to-continuum ratios are similar for the two disks in each of the lines. From the resolved velocity gradients across each disk, we constrain the masses of the central stars, and show consistency with optical-infrared spectroscopy, both indicative of a high mass ratio ~9. The small difference between the systemic velocities indicates that the binary orbital plane is close to face-on. The angle between the projected disk rotation axes is very high, ~72 degrees, showing that the system did not form from a single massive disk or a rigidly rotating cloud core. This finding, which adds to related evidence from disk geometries in other systems, protostellar outflows, stellar rotation, and similar recent ALMA results, demonstrates that turbulence or dynamical interactions act on small scales well below that of molecular cores during the early stages of star formation.
[131]  oai:arXiv.org:1404.5652  [pdf] - 1209105
Resolved Multifrequency Radio Observations of GG Tau
Comments: ApJ, in press
Submitted: 2014-04-22
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.
[132]  oai:arXiv.org:1404.0709  [pdf] - 1208776
Herschel evidence for disk flattening or gas depletion in transitional disks
Comments: Accepted for publication in ApJ; 52 pages, 16 figures, 8 tables
Submitted: 2014-04-02
Transitional disks are protoplanetary disks characterized by reduced near- and mid-infrared emission with respect to full disks. This characteristic spectral energy distribution indicates the presence of an optically thin inner cavity within the dust disk believed to mark the disappearance of the primordial massive disk. We present new Herschel Space Observatory PACS spectra of [OI] 63 micron for 21 transitional disks. Our survey complements the larger Herschel GASPS program "Gas in Protoplanetary Systems" (Dent et al. 2013) by quadrupling the number of transitional disks observed with PACS at this wavelength. [OI] 63 micron traces material in the outer regions of the disk, beyond the inner cavity of most transitional disks. We find that transitional disks have [OI] 63 micron line luminosities two times fainter than their full disk counterparts. We self consistently determine various stellar properties (e.g. bolometric luminosity, FUV excess, etc.) and disk properties (e.g. disk dust mass, etc.) that could influence the [OI] 63 micron line luminosity and we find no correlations that can explain the lower [OI] 63 micron line luminosities in transitional disks. Using a grid of thermo-chemical protoplanetary disk models, we conclude that either transitional disks are less flared than full disks or they possess lower gas-to-dust ratios due to a depletion of gas mass. This result suggests that transitional disks are more evolved than their full disk counterparts, possibly even at large radii.
[133]  oai:arXiv.org:1403.2026  [pdf] - 1208269
ALMA Observations of the Orion Proplyds
Comments: ApJ, in press
Submitted: 2014-03-08
We present ALMA observations of protoplanetary disks ("proplyds") in the Orion Nebula Cluster. We imaged 5 individual fields at 856um containing 22 HST-identified proplyds and detected 21 of them. Eight of those disks were detected for the first time at submillimeter wavelengths, including the most prominent, well-known proplyd in the entire Orion Nebula, 114-426. Thermal dust emission in excess of any free-free component was measured in all but one of the detected disks, and ranged between 1-163 mJy, with resulting disk masses of 0.3-79 Mjup. An additional 26 stars with no prior evidence of associated disks in HST observations were also imaged within the 5 fields, but only 2 were detected. The disk mass upper limits for the undetected targets, which include OB stars, theta1Ori C and theta1Ori F, range from 0.1-0.6 Mjup. Combining these ALMA data with previous SMA observations, we find a lack of massive (>3 Mjup) disks in the extreme-UV dominated region of Orion, within 0.03 pc of O-star theta1Ori C. At larger separations from theta1Ori C, in the far-UV dominated region, there is a wide range of disk masses, similar to what is found in low-mass star forming regions. Taken together, these results suggest that a rapid dissipation of disk masses likely inhibits potential planet formation in the extreme-UV dominated regions of OB associations, but leaves disks in the far-UV dominated regions relatively unaffected.
[134]  oai:arXiv.org:1402.7103  [pdf] - 1203630
An Observational Perspective of Transitional Disks
Comments: 24 pages, 13 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning
Submitted: 2014-02-27
Transitional disks are objects whose inner disk regions have undergone substantial clearing. The Spitzer Space Telescope produced detailed spectral energy distributions (SEDs) of transitional disks that allowed us to infer their radial dust disk structure in some detail, revealing the diversity of this class of disks. The growing sample of transitional disks also opened up the possibility of demographic studies, which provided unique insights. There now exist (sub)millimeter and infrared images that confirm the presence of large clearings of dust in transitional disks. In addition, protoplanet candidates have been detected within some of these clearings. Transitional disks are thought to be a strong link to planet formation around young stars and are a key area to study if further progress is to be made on understanding the initial stages of planet formation. Here we provide a review and synthesis of transitional disk observations to date with the aim of providing timely direction to the field, which is about to undergo its next burst of growth as ALMA reaches its full potential. We discuss what we have learned about transitional disks from SEDs, color-color diagrams, and imaging in the (sub)mm and infrared. We then distill the observations into constraints for the main disk clearing mechanisms proposed to date (i.e., photoevaporation, grain growth, and companions) and explore how the expected observational signatures from these mechanisms, particularly planet-induced disk clearing, compare to actual observations. Lastly, we discuss future avenues of inquiry to be pursued with ALMA, JWST, and next generation of ground-based telescopes.
[135]  oai:arXiv.org:1402.6597  [pdf] - 820798
On the structure of the transition disk around TW Hya
Comments: 22 pages, 12 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2014-02-26
For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.
[136]  oai:arXiv.org:1402.1354  [pdf] - 1203212
Dust Evolution in Protoplanetary Disks
Comments: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning
Submitted: 2014-02-06
(abridged) In the core accretion scenario for the formation of planetary rocky cores, the first step toward planet formation is the growth of dust grains into larger and larger aggregates and eventually planetesimals. Although dust grains are thought to grow from the submicron sizes typical of interstellar dust to micron size particles in the dense regions of molecular clouds and cores, the growth from micron size particles to pebbles and kilometre size bodies must occur in protoplanetary disks. This step in the formation of planetary systems is the last stage of solids evolution that can be observed directly in young extrasolar systems. In this chapter we review the constraints on the physics of grain-grain collisions as they have emerged from laboratory experiments and numerical computations. We then review the current theoretical understanding of the global processes governing the evolution of solids in protoplanetary disks, including dust settling, growth, and radial transport. The predicted observational signatures are summarized. We discuss recent developments in the study of grain growth in molecular cloud cores and in collapsing envelopes of protostars as these provide the initial conditions for the dust in disks. We discuss the observational evidence for the growth of grains in young disks from mm surveys, as well as the recent evidence of radial variations of the dust properties in disks. We include a brief discussion of the constraints on the small end of the grain size distribution and on dust settling as derived from optical and IR observations. The observations are discussed in the context of global dust evolution models, in particular we focus on the emerging evidence for a very efficient early growth of grains and the radial distribution of grain sizes in disks. We also highlight the limits of current models, including the need to slow the radial drift of grains.
[137]  oai:arXiv.org:1402.1538  [pdf] - 1203230
High-Resolution Submillimeter and Near-Infrared Studies of the Transition Disk around Sz 91
Comments:
Submitted: 2014-02-06
To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.
[138]  oai:arXiv.org:1312.3817  [pdf] - 1202234
Fast Radial Flows in Transition Disk Holes
Comments: 21 pages, 8 figures, Accepted to ApJ
Submitted: 2013-12-13
Protoplanetary "transition" disks have large, mass-depleted central cavities, yet also deliver gas onto their host stars at rates comparable to disks without holes. The paradox of simultaneous transparency and accretion can be explained if gas flows inward at much higher radial speeds inside the cavity than outside the cavity, since surface density (and by extension optical depth) varies inversely with inflow velocity at fixed accretion rate. Radial speeds within the cavity might even have to approach free-fall values to explain the huge surface density contrasts inferred for transition disks. We identify observational diagnostics of fast radial inflow in channel maps made in optically thick spectral lines. Signatures include (1) twisted isophotes in maps made at low systemic velocities and (2) rotation of structures observed between maps made in high-velocity line wings. As a test case, we apply our new diagnostic tools to archival ALMA data on the transition disk HD 142527, and uncover evidence for free-fall radial velocities inside its cavity. Although the observed kinematics are also consistent with a disk warp, the radial inflow scenario is preferred because it predicts low surface densities that appear consistent with recent observations of optically thin CO isotopologues in this disk. How material in the disk cavity sheds its angular momentum wholesale to fall freely onto the star is an unsolved problem; gravitational torques exerted by giant planets or brown dwarfs are briefly discussed as a candidate mechanism.
[139]  oai:arXiv.org:1311.5222  [pdf] - 1201767
On the Outer Edges of Protoplanetary Dust Disks
Comments: 12 pages, 5 figures, accepted to ApJ
Submitted: 2013-11-20
The expectation that aerodynamic drag will force the solids in a gas-rich protoplanetary disk to spiral in toward the host star on short timescales is one of the fundamental problems in planet formation theory. The nominal efficiency of this radial drift process is in conflict with observations, suggesting that an empirical calibration of solid transport mechanisms in a disk is highly desirable. However, the fact that both radial drift and grain growth produce a similar particle size segregation in a disk (such that larger particles are preferentially concentrated closer to the star) makes it difficult to disentangle a clear signature of drift alone. We highlight a new approach, by showing that radial drift leaves a distinctive "fingerprint" in the dust surface density profile that is directly accessible to current observational facilities. Using an analytical framework for dust evolution, we demonstrate that the combined effects of drift and (viscous) gas drag naturally produce a sharp outer edge in the dust distribution (or, equivalently, a sharp decrease in the dust-to-gas mass ratio). This edge feature forms during the earliest phase in the evolution of disk solids, before grain growth in the outer disk has made much progress, and is preserved over longer timescales when both growth and transport effects are more substantial. The key features of these analytical models are reproduced in detailed numerical simulations, and are qualitatively consistent with recent millimeter-wave observations that find gas/dust size discrepancies and steep declines in dust continuum emission in the outer regions of protoplanetary disks.
[140]  oai:arXiv.org:1311.1819  [pdf] - 1180576
The Dispersal of Protoplanetary Disks
Comments: 23 pages, 6 figures. Refereed review chapter, accepted for publication in Protostars & Planets VI, University of Arizona Press (2014), eds. H.Beuther, C.Dullemond, Th.Henning, R.Klessen
Submitted: 2013-11-07
Protoplanetary disks are the sites of planet formation, and the evolution and eventual dispersal of these disks strongly influences the formation of planetary systems. Disk evolution during the planet-forming epoch is driven by accretion and mass-loss due to winds, and in typical environments photoevaporation by high-energy radiation from the central star is likely to dominate final gas disk dispersal. We present a critical review of current theoretical models, and discuss the observations that are used to test these models and inform our understanding of the underlying physics. We also discuss the role disk dispersal plays in shaping planetary systems, considering its influence on both the process(es) of planet formation and the architectures of planetary systems. We conclude by presenting a schematic picture of protoplanetary disk evolution and dispersal, and discussing prospects for future work.
[141]  oai:arXiv.org:1311.1569  [pdf] - 1180555
Detection, size, measurement and structural analysis limits for the 2MASS, UKIDSS-LAS & VISTA VIKING surveys
Comments: 16 pages, 15 figures, 7 tables. Accepted for publication in Publications of the Astronomical Society of Australia
Submitted: 2013-11-06
The 2MASS, UKIDSS-LAS and VISTA VIKING surveys have all now observed the GAMA 9hr region in the $K_s$ band. Here we compare the detection rates, photometry, basic size measurements, and single-component GALFIT structural measurements for a sample of 37,591 galaxies. We explore the sensitivity limits where the data agree for a variety of issues including: detection, star-galaxy separation, photometric measurements, size and ellipticity measurements, and S\'{e}rsic measurements. We find that 2MASS fails to detect at least 20% of the galaxy population within all magnitude bins, however for those that are detected we find photometry is robust ($\pm$ 0.2 mag) to 14.7 AB mag and star-galaxy separation to 14.8 AB mag. For UKIDSS-LAS we find incompleteness starts to enter at a flux limit of 18.9 AB mag, star-galaxy separation is robust to 16.3 AB mag and structural measurements are robust to 17.7 AB mag. VISTA VIKING data is complete to approximately 20.0 AB mag and structural measurements appear robust to 18.8 AB mag.
[142]  oai:arXiv.org:1308.6575  [pdf] - 1178882
Dust Properties of Local Dust-Obscured Galaxies with the Submillimeter Array
Comments: 8 pages, 4 figures. To appear in ApJ
Submitted: 2013-08-29
We report Submillimeter Array (SMA) observations of the 880 \mu m dust continuum emission for four dust-obscured galaxies (DOGs) in the local universe. Two DOGs are clearly detected with S_\nu(880 \mu m)=10-13 mJy and S/N>5, but the other two are not detected with 3{\sigma} upper limits of S_\nu(880 \mu m)=5-9 mJy. Including an additional two local DOGs with submillimeter data from the literature, we determine the dust masses and temperatures for six local DOGs. The infrared luminosities and dust masses for these DOGs are in the range 1.2-4.9x10^{11} (L_sun) and 4-14x10^{7} (M_sun), respectively. The dust temperatures derived from a two-component modified blackbody function are 23-26 K and 60-124 K for the cold and warm dust components, respectively. Comparison of local DOGs with other infrared luminous galaxies with submillimeter detections shows that the dust temperatures and masses do not differ significantly among these objects. Thus, as argued previously, local DOGs are not a distinctive population among dusty galaxies, but simply represent the high-end tail of the dust obscuration distribution.
[143]  oai:arXiv.org:1308.4358  [pdf] - 1173510
The Structure of the Evolved Circumbinary Disk around V4046 Sgr
Comments: 31 pages, 8 figures, accepted for publication in ApJ
Submitted: 2013-08-20
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.
[144]  oai:arXiv.org:1307.7439  [pdf] - 699304
Imaging of the CO Snow Line in a Solar Nebula Analog
Comments: 21 pages, 7 figures. This is the author's pre-editorial version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science Express 18 July 2013. Full version can be found at http://dx.doi.org/10.1126/science.1239560
Submitted: 2013-07-28
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.
[145]  oai:arXiv.org:1307.7174  [pdf] - 1173015
A SCUBA-2 850 micron survey of protoplanetary discs in the sigma Orionis cluster
Comments: 10 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2013-07-26
We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution.
[146]  oai:arXiv.org:1307.5848  [pdf] - 1172890
An azimuthal asymmetry in the LkHa 330 disk
Comments: ApJ, in press, 12 pages, 6 figures
Submitted: 2013-07-22
Theory predicts that giant planets and low mass stellar companions shape circumstellar disks by opening annular gaps in the gas and dust spatial distribution. For more than a decade it has been debated whether this is the dominant process that leads to the formation of transitional disks. In this paper, we present millimeter-wave interferometric observations of the transitional disk around the young intermediate mass star LkHa330. These observations reveal a lopsided ring in the 1.3 mm dust thermal emission characterized by a radius of about 100 AU and an azimuthal intensity variation of a factor of 2. By comparing the observations with a Gaussian parametric model, we find that the observed asymmetry is consistent with a circular arc, that extends azimuthally by about 90 deg and emits about 1/3 of the total continuum flux at 1.3 mm. Hydrodynamic simulations show that this structure is similar to the azimuthal asymmetries in the disk surface density that might be produced by the dynamical interaction with unseen low mass companions orbiting within 70 AU from the central star. We argue that such asymmetries might lead to azimuthal variations in the millimeter-wave dust opacity and in the dust temperature, which will also affect the millimeter-wave continuum emission. Alternative explanations for the observed asymmetry that do not require the presence of companions cannot be ruled out with the existing data. Further observations of both the dust and molecular gas emission are required to derive firm conclusions on the origin of the asymmetry observed in the LkHa330 disk.
[147]  oai:arXiv.org:1307.3560  [pdf] - 1172705
Resolving The Moth at Millimeter Wavelengths
Comments: 10 pages, 6 figures
Submitted: 2013-07-12
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.
[148]  oai:arXiv.org:1306.6475  [pdf] - 1172332
A Spatially Resolved Vertical Temperature Gradient in the HD 163296 Disk
Comments: 31 pages, 13 figures, accepted for publication in ApJ
Submitted: 2013-06-27
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%).
[149]  oai:arXiv.org:1305.5262  [pdf] - 1171534
The Mass Dependence Between Protoplanetary Disks and their Stellar Hosts
Comments: ApJ, in press: electronic tables are available before publication, by request (typo in Table 2 now corrected)
Submitted: 2013-05-22, last modified: 2013-06-13
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).
[150]  oai:arXiv.org:1306.0275  [pdf] - 1171736
GASPS - a Herschel survey of gas and dust in Protoplanetary Disks: Summary and Initial Statistics
Comments: Accepted for PASP
Submitted: 2013-06-02
GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.
[151]  oai:arXiv.org:1303.1107  [pdf] - 1944130
An Old Disk That Can Still Form a Planetary System
Comments: 16 pages, 3 figures. Pre-editorial version of paper published in Nature. Full version can be found at http://www.nature.com.proxy.lib.umich.edu/nature/journal/v493/n7434/full/nature11805.html
Submitted: 2013-03-05
From the masses of planets orbiting our Sun, and relative elemental abundances, it is estimated that at birth our Solar System required a minimum disk mass of ~0.01 solar masses within ~100 AU of the star. The main constituent, gaseous molecular hydrogen, does not emit from the disk mass reservoir, so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally probes the disk surface, while the conversion from dust emission to gas mass requires knowledge of the grain properties and gas-to-dust mass ratio, which likely differ from their interstellar values. Thus, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3--10 Myr) star TW Hya, with estimates ranging from 0.0005 to 0.06 solar masses. Here we report the detection the fundamental rotational transition of hydrogen deuteride, HD, toward TW Hya. HD is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The HD detection, combined with existing observations and detailed models, implies a disk mass >0.05 solar masses, enough to form a planetary system like our own.
[152]  oai:arXiv.org:1302.5705  [pdf] - 709517
Mapping H-band Scattered Light Emission in the Mysterious SR21 Transitional Disk
Comments:
Submitted: 2013-02-22
We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.
[153]  oai:arXiv.org:1212.1466  [pdf] - 1158341
Spiral Arms in the Asymmetrically Illuminated Disk of MWC 758 and Constraints on Giant Planets
Comments: 45 pages, 9 figures; ApJ in press
Submitted: 2012-12-06
We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.
[154]  oai:arXiv.org:1211.5148  [pdf] - 1157970
Millimeter Emission Structure in the first ALMA Image of the AU Mic Debris Disk
Comments: 13 pages, 2 figures, 2 tables, accepted by ApJ Letters
Submitted: 2012-11-21
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.
[155]  oai:arXiv.org:1210.5252  [pdf] - 1152270
Constraints on the Radial Variation of Grain Growth in the AS 209 Circumstellar Disk
Comments: 7 pages, 5 figures, Accepted to ApJL
Submitted: 2012-10-18
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 \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} 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.
[156]  oai:arXiv.org:1209.4407  [pdf] - 1151516
A Disk-based Dynamical Mass Estimate for the Young Binary V4046 Sgr
Comments: 24 pages, 7 figures, accepted for publication in ApJ
Submitted: 2012-09-19
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.
[157]  oai:arXiv.org:1209.1641  [pdf] - 1151240
Matryoshka Holes: Nested Emission Rings in the Transitional Disk Oph IRS 48
Comments: 5 pages, 2 figures, accepted to ApJ Letters
Submitted: 2012-09-07
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.
[158]  oai:arXiv.org:1208.1285  [pdf] - 1150553
Kinematics of the CO Gas in the Inner Regions of the TW Hya Disk
Comments: 30 pages, 7 figures, accepted for publication in ApJ
Submitted: 2012-08-06
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.
[159]  oai:arXiv.org:1206.4688  [pdf] - 1124261
Keck/NIRC2 Imaging of the Warped, Asymmetric Debris Disk around HD 32297
Comments: 19 pages, 12 figures; ApJ in press (no text changes from previous version)
Submitted: 2012-06-20, last modified: 2012-07-17
We present Keck/NIRC2 $K_{s}$ band high-contrast coronagraphic imaging of the luminous debris disk around the nearby, young A star HD 32297 resolved at a projected separation of $r$ = 0.3-2.5\arcsec{} ($\approx$ 35-280 AU). The disk is highly warped to the north and exhibits a complex, "wavy" surface brightness profile interior to $r$ $\approx$ 110 AU, where the peaks/plateaus in the profiles are shifted between the NE and SW disk lobes. The SW side of the disk is 50--100% brighter at $r$ = 35-80 AU, and the location of its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution modeling suggests that HD 32297 has at least two dust populations that may originate from two separate belts likely at different locations, possibly at distances coinciding with the surface brightness peaks. A disk model for a single dust belt including a phase function with two components and a 5-10 AU pericenter offset explains the disk's warped structure and reproduces some of the surface brightness profile's shape (e.g. the overall "wavy" profile, the SB peak/plateau shifts) but more poorly reproduces the disk's brightness asymmetry. Although there may be alternate explanations, agreement between the SW disk brightness peak and disk's peak mm emission is consistent with an overdensity of very small, sub-blowout-sized dust and large, 0.1-1 mm-sized grains at $\approx$ 45 AU tracing the same parent population of planetesimals. New near-IR and submm observations may be able to clarify whether even more complex grain scattering properties or dynamical sculpting by an unseen planet are required to explain HD 32297's disk structure.
[160]  oai:arXiv.org:1206.5802  [pdf] - 530623
Can grain growth explain transition disks?
Comments: 11 pages, 5 figures, accepted to A&A
Submitted: 2012-06-25
Aims: Grain growth has been suggested as one possible explanation for the diminished dust optical depths in the inner regions of protoplanetary "transition" disks. In this work, we directly test this hypothesis in the context of current models of grain growth and transport. Methods: A set of dust evolution models with different disk shapes, masses, turbulence parameters, and drift efficiencies is combined with radiative transfer calculations in order to derive theoretical spectral energy distributions (SEDs) and images. Results: We find that grain growth and transport effects can indeed produce dips in the infrared SED, as typically found in observations of transition disks. Our models achieve the necessary reduction of mass in small dust by producing larger grains, yet not large enough to be fragmenting efficiently. However, this population of large grains is still detectable at millimeter wavelengths. Even if perfect sticking is assumed and radial drift is neglected, a large population of dust grains is left behind because the time scales on which they are swept up by the larger grains are too long. This mechanism thus fails to reproduce the large emission cavities observed in recent millimeter-wave interferometric images of accreting transition disks.
[161]  oai:arXiv.org:1203.6353  [pdf] - 1117620
A Resolved Census of Millimeter Emission from Taurus Multiple Star Systems
Comments: Accepted for publication in ApJ. 40 pages, 11 figures, 9 tables
Submitted: 2012-03-28
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.
[162]  oai:arXiv.org:1203.1896  [pdf] - 1117171
A Resolved Millimeter Emission Belt in the AU Mic Debris Disk
Comments: 10 pages, 2 figures. Accepted for publication in ApJ Letters
Submitted: 2012-03-08
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.
[163]  oai:arXiv.org:1203.0318  [pdf] - 1117026
Confirming the Primarily Smooth Structure of the Vega Debris Disk at Millimeter Wavelengths
Comments: 11 pages, 7 figures, accepted for publication in ApJ
Submitted: 2012-03-01
Clumpy structure in the debris disk around Vega has been previously reported at millimeter wavelengths and attributed to concentrations of dust grains trapped in resonances with an unseen planet. However, recent imaging at similar wavelengths with higher sensitivity has disputed the observed structure. We present three new millimeter-wavelength observations that help to resolve the puzzling and contradictory observations. We have observed the Vega system with the Submillimeter Array (SMA) at a wavelength of 880 um and angular resolution of 5"; with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) at a wavelength of 1.3 mm and angular resolution of 5"; and with the Green Bank Telescope (GBT) at a wavelength of 3.3 mm and angular resolution of 10". Despite high sensitivity and short baselines, we do not detect the Vega debris disk in either of the interferometric data sets (SMA and CARMA), which should be sensitive at high significance to clumpy structure based on previously reported observations. We obtain a marginal (3-sigma) detection of disk emission in the GBT data; the spatial distribution of the emission is not well constrained. We analyze the observations in the context of several different models, demonstrating that the observations are consistent with a smooth, broad, axisymmetric disk with inner radius 20-100 AU and width >50 AU. The interferometric data require that at least half of the 860 um emission detected by previous single-dish observations with the James Clerk Maxwell Telescope be distributed axisymmetrically, ruling out strong contributions from flux concentrations on spatial scales of <100 AU. These observations support recent results from the Plateau de Bure Interferometer indicating that previous detections of clumpy structure in the Vega debris disk were spurious.
[164]  oai:arXiv.org:1201.1518  [pdf] - 1092782
On the Transitional Disk Class: Linking Observations of T Tauri Stars & Physical Disk Models
Comments: accepted to ApJ
Submitted: 2012-01-06
Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the disk and we call these objects "transitional" and pre-transitional" disks, respectively. For the latter group, we can fit these SEDs with full disk models and find that millimeter data are necessary to break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only be applied to objects that display evidence for a radical change in the disk's radial structure. Using this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. These reduced accretion rates onto the star could be linked to forming planets. Future observations of transitional and pre-transitional disks will allow us to better quantify the signatures of planet formation in young disks.
[165]  oai:arXiv.org:1111.5037  [pdf] - 1091822
The TW Hya Disk at 870 microns: Comparison of CO and Dust Radial Structures
Comments: ApJ, in press (fixed typo in Equation 4)
Submitted: 2011-11-21, last modified: 2011-12-16
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.
[166]  oai:arXiv.org:1110.3865  [pdf] - 1084925
A Closer Look at the LkCa 15 Protoplanetary Disk
Comments: in press at ApJ Letters
Submitted: 2011-10-17
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.
[167]  oai:arXiv.org:1109.2578  [pdf] - 1084008
The Ionization Fraction in the DM Tau Protoplanetary Disk
Comments: Accepted for publication in ApJ. 20 pages, including 5 figures
Submitted: 2011-09-12
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.
[168]  oai:arXiv.org:1107.5061  [pdf] - 1078207
Resolving the CO Snow Line in the Disk around HD 163296
Comments: 40 pages, 13 figures, accepted for publication in ApJ
Submitted: 2011-07-25
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.
[169]  oai:arXiv.org:1107.3153  [pdf] - 1078012
Resolved Submillimeter Observations of the HR 8799 and HD 107146 Debris Disks
Comments: 9 pages, 3 figures, accepted for publication in ApJ
Submitted: 2011-07-15
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.
[170]  oai:arXiv.org:1104.1236  [pdf] - 1053224
Disk Imaging Survey of Chemistry with SMA: II. Southern Sky Protoplanetary Disk Data and Full Sample Statistics
Comments: accepted for publication in ApJ, 41 pages including 7 figures
Submitted: 2011-04-06
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.
[171]  oai:arXiv.org:1103.0284  [pdf] - 1052438
Resolved Images of Large Cavities in Protoplanetary Transition Disks
Comments: ApJ, in press
Submitted: 2011-03-01
Circumstellar disks are thought to experience a rapid "transition" phase in their evolution that can have a considerable impact on the formation and early development of planetary systems. We present new and archival high angular resolution (0.3" = 40-75 AU) Submillimeter Array (SMA) observations of the 880 micron dust continuum emission from 12 such transition disks in nearby star-forming regions. In each case, we directly resolve a dust-depleted disk cavity around the central star. Using radiative transfer calculations, we interpret these dust disk structures in a homogeneous, parametric model framework by reproducing their SMA visibilities and SEDs. The cavities in these disks are large (R_cav = 15-73 AU) and substantially depleted of small (~um-sized) dust grains, although their mass contents are still uncertain. The structures of the remnant material at larger radii are comparable to normal disks. We demonstrate that these large cavities are common among the millimeter-bright disk population, comprising at least 20% of the disks in the bright half of the millimeter luminosity (disk mass) distribution. Utilizing these results, we assess some of the physical mechanisms proposed to account for transition disk structures. As has been shown before, photoevaporation models do not produce the large cavity sizes, accretion rates, and disk masses representative of this sample. It would be difficult to achieve a sufficient decrease of the dust optical depths in these cavities by particle growth alone: substantial growth (to meter sizes or beyond) must occur in large (tens of AU) regions of low turbulence without also producing an abundance of small particles. Given those challenges, we suggest instead that the observations are most commensurate with dynamical clearing due to tidal interactions with low-mass companions --young brown dwarfs or giant planets on long-period orbits.
[172]  oai:arXiv.org:1012.2313  [pdf] - 1042493
Millimeter Imaging of the beta Pictoris Debris Disk: Evidence for a Planetesimal Belt
Comments: 9 pages, 2 figures, accepted by ApJL
Submitted: 2010-12-10
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.
[173]  oai:arXiv.org:1011.3826  [pdf] - 1042010
Empirical Constraints on Turbulence in Protoplanetary Accretion Disks
Comments: 18 pages, 9 figures, accepted for publication in ApJ
Submitted: 2010-11-16
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.
[174]  oai:arXiv.org:1007.5070  [pdf] - 1033983
Protoplanetary Disk Structures in Ophiuchus II: Extension to Fainter Sources
Comments: accepted in ApJ, 39 pages, 10 figures
Submitted: 2010-07-28
We present new results from a significant extension of our previous high angular resolution (0.3" = 40 AU) Submillimeter Array survey of the 880 um continuum emission from dusty circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. An expanded sample is constructed to probe disk structures that emit significantly lower millimeter luminosities (hence dust masses), down to the median value for T Tauri stars. Using a Monte Carlo radiative transfer code, the millimeter visibilities and broadband spectral energy distribution for each disk are simultaneously reproduced with a two-dimensional parametric model for a viscous accretion disk. We find wide ranges of characteristic radii (14-198 AU) and disk masses (0.004-0.143 M_sun), but a narrow distribution of surface density gradients (0.4-1.1) that is consistent with a uniform value $\gamma$ = 0.9 +/- 0.2 and independent of mass (or millimeter luminosity). In this sample, we find a correlation between the disk luminosity/mass and characteristic radius, such that fainter disks are both smaller and less massive. We suggest that this relationship is an imprint of the initial conditions inherited by the disks at their formation epoch, compare their angular momenta with those of molecular cloud cores, and speculate on how future observations can help constrain the distribution of viscous evolution timescales. No other correlations between disk and star properties are found. The inferred disk structures are briefly compared with theoretical models for giant planet formation, although resolution limitations do not permit us to directly comment on material inside R = 20 AU. However, there is some compelling evidence for dust evolution in the planet formation region: 4/17 disks in the sample show resolved regions of significantly reduced optical depths within ~20-40 AU of their central stars.
[175]  oai:arXiv.org:1007.3267  [pdf] - 1033773
Structure and Composition of Two Transitional Circumstellar Disks in Corona Australis
Comments: 11 pages, 5 figures, accepted for publication in AJ
Submitted: 2010-07-19
The late stages of evolution of the primordial circumstellar disks surrounding young stars are poorly understood, yet vital to constrain theories of planet formation. We consider basic structural models for the disks around two ~10 Myr-old members of the nearby RCrA association, RX J1842.9-3532 and RX J1852.3-3700. We present new arcsecond-resolution maps of their 230 GHz continuum emission from the Submillimeter Array and unresolved CO(3-2) spectra from the Atacama Submillimeter Telescope Experiment. By combining these data with broadband fluxes from the literature and infrared fluxes and spectra from the catalog of the Formation and Evolution of Planetary Systems (FEPS) Legacy program on the Spitzer Space Telescope, we assemble a multiwavelength data set probing the gas and dust disks. Using the Monte Carlo radiative transfer code RADMC to model simultaneously the SED and millimeter continuum visibilities, we derive basic dust disk properties and identify an inner cavity of radius 16 AU in the disk around RX J1852.3-3700. We also identify an optically thin 5 AU cavity in the disk around RX J1842.9-3532, with a small amount of optically thick material close to the star. The molecular line observations suggest an intermediate disk inclination in RX J1842.9-3532, consistent with the continuum emission. In combination with the dust models, the molecular data allow us to derive a lower CO content than expected, suggesting that the process of gas clearing is likely underway in both systems, perhaps simultaneously with planet formation.
[176]  oai:arXiv.org:1007.1476  [pdf] - 1033588
Disk Imaging Survey of Chemistry with SMA (DISCS): I. Taurus Protoplanetary Disk Data
Comments: 29 pages, 4 figures, accepted for publication in ApJ
Submitted: 2010-07-08
Chemistry plays an important role in the structure and evolution of protoplanetary disks, with implications for the composition of comets and planets. This is the first of a series of papers based on data from DISCS, a Submillimeter Array survey of the chemical composition of protoplanetary disks. The six Taurus sources in the program (DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) range in stellar spectral type from M1 to A4 and offer an opportunity to test the effects of stellar luminosity on the disk chemistry. The disks were observed in 10 different lines at ~3" resolution and an rms of ~100 mJy beam-1 at ~0.5 km s-1. The four brightest lines are CO 2-1, HCO+ 3-2, CN 2_3-1_2 and HCN 3-2 and these are detected toward all sources (except for HCN toward CQ Tau). The weaker lines of CN 2_2-1_1, DCO+ 3-2, N2H+ 3-2, H2CO 3_03-2_02 and 4_14-3_13 are detected toward two to three disks each, and DCN 3-2 only toward LkCa 15. CH3OH 4_21-3_12 and c-C3H2 are not detected. There is no obvious difference between the T Tauri and Herbig Ae sources with regard to CN and HCN intensities. In contrast, DCO+, DCN, N2H+ and H2CO are detected only toward the T Tauri stars, suggesting that the disks around Herbig Ae stars lack cold regions for long enough timescales to allow for efficient deuterium chemistry, CO freeze-out, and grain chemistry.
[177]  oai:arXiv.org:1005.3828  [pdf] - 1032654
Gas in the protoplanetary disc of HD 169142: Herschel's view
Comments: 5 pages, 3 Postscript figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.
[178]  oai:arXiv.org:1005.3864  [pdf] - 1032661
Gas in Protoplanetary Discs (GASPS) 1. First results
Comments: Comments: 5 pages, 2 Postscript figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.
[179]  oai:arXiv.org:1005.3826  [pdf] - 1032652
The Herschel view of Gas in Protoplanetary Systems (GASPS). First comparisons with a large grid of models
Comments: 5 pages, 4 figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.
[180]  oai:arXiv.org:1005.3830  [pdf] - 1032656
Herschel-PACS observation of the 10 Myr old T Tauri disk TW Hya: Constraining the disk gas mass
Comments: 6 pages, 6 figures, accepted for publication in the A&A Herschel special issue
Submitted: 2010-05-20
Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.
[181]  oai:arXiv.org:0912.3537  [pdf] - 902791
Truncated Disks in TW Hya Association Multiple Star Systems
Comments: ApJ, in press
Submitted: 2009-12-17
We present high angular resolution (down to 0.3" = 13 AU in diameter) Submillimeter Array observations of the 880 micron (340 GHz) thermal continuum emission from circumstellar dust disks in the nearby HD 98800 and Hen 3-600 multiple star systems. In both cases, the dust emission is resolved and localized around one stellar component - the HD 98800 B and Hen 3-600 A spectroscopic binaries - with no evidence for circum-system material. Using two-dimensional Monte Carlo radiative transfer calculations, we compare the SMA visibilities and broadband spectral energy distributions with truncated disk models to empirically locate the inner and outer edges of both disks. The HD 98800 B disk appears to be aligned with the spectroscopic binary orbit, is internally truncated at a radius of 3.5 AU, and extends to only 10-15 AU from the central stars. The Hen 3-600 A disk is slightly larger, with an inner edge at 1 AU and an outer radius of 15-25 AU. These inferred disk structures compare favorably with theoretical predictions of their truncation due to tidal interactions with the stellar companions.
[182]  oai:arXiv.org:0912.2433  [pdf] - 1018745
Millimeter Dust Emission in the GQ Lup System
Comments: 16 pages, 2 figures, accepted to AJ
Submitted: 2009-12-12
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.
[183]  oai:arXiv.org:0909.1345  [pdf] - 901646
Stringent Limits on the Polarized Submillimeter Emission from Protoplanetary Disks
Comments: 15 pages, 6 figures, accepted for publication in ApJ
Submitted: 2009-09-07
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.
[184]  oai:arXiv.org:0906.0730  [pdf] - 24906
Protoplanetary Disk Structures in Ophiuchus
Comments: ApJ in press: 51 pages, 13 figures
Submitted: 2009-06-03
We present a high angular resolution (0.3" = 40 AU) SMA survey of the 870 micron thermal continuum emission from 9 of the brightest, and therefore most massive, circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. Using 2-D radiative transfer calculations, we simultaneously fit the observed continuum visibilities and broadband spectral energy distribution for each disk with a parametric structure model. Compared to previous millimeter studies, this survey includes significant upgrades in modeling, data quality, and angular resolution that provide improved constraints on key structure parameters, particularly those that characterize the spatial distribution of mass in the disks. In the context of a surface density profile motivated by similarity solutions for viscous accretion disks, the best-fit models for the sample disks have characteristic radii R_c = 20-200 AU, high disk masses M_d = 0.005-0.14 M_sun, and a narrow range of radial surface density gradients around a median $\gamma$ = 0.9. These density structures are used in conjunction with accretion rate estimates from the literature to help characterize the viscous evolution of the disk material. Using the standard prescription for disk viscosities, those combined constraints indicate that $\alpha$ = 0.0005-0.08. Three of the sample disks show large (R = 20-40 AU) central cavities in their continuum emission morphologies, marking extensive zones where dust has been physically removed and/or has significantly diminished opacities. Based on the current requirements of planet formation models, these emission cavities and the structure constraints for the sample as a whole suggest that these young disks may eventually produce planetary systems, and have perhaps already started. (abridged)
[185]  oai:arXiv.org:0903.4455  [pdf] - 315756
A Spatially Resolved Inner Hole in the Disk around GM Aurigae
Comments: 12 pages, 6 figures, accepted for publication in ApJ
Submitted: 2009-03-25
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.
[186]  oai:arXiv.org:0809.2276  [pdf] - 16272
The LkHa 101 Cluster
Comments: to appear in the Handbook of Star Forming Regions Vol. I (ASP Monograph Series, ed. B. Reipurth)
Submitted: 2008-09-12
In the infrared, the heavily reddened LkH$\alpha$ 101 is one of the brightest young stars in the sky. Situated just north of the Taurus-Auriga complex in the L1482 dark cloud, it appears to be an early B-type star that has been serendipitously exposed during a rarely observed stage of early evolution, revealing a remarkable spectrum and a directly-imaged circumstellar disk. While detailed studies of this star and its circumstellar environment have become increasingly sophisticated in the 50 years since Herbig (1956) first pointed it out, the true nature of the object still remains a mystery. Recent work has renewed focus on the young cluster of stars surrounding LkH$\alpha$ 101, and what it can tell us about the enigmatic source at its center (e.g., massive star formation timescales, clustered formation mechanisms). This latter effort certainly deserves more intensive study. We describe the current knowledge of this region and point out interesting work that could be done in the future.
[187]  oai:arXiv.org:0806.4163  [pdf] - 13897
Submillimeter Observations of the Young Low-Mass Object IRAS 04158+2805
Comments: accepted for publication in ApJ; 18 pages, 5 figures
Submitted: 2008-06-25
We present high spatial resolution Submillimeter Array observations and supplementary single-dish photometry of the molecular gas and dust around IRAS 04158+2805, a young source with spectral type M5-M6 in the Taurus star-forming region. A bright, highly elongated dust structure that extends 8" (~1120 AU) in diameter is revealed in a 883 micron thermal continuum image. This emission geometry is in good agreement with optical observations that show a similar structure in absorption, aligned perpendicular to bipolar scattered light nebulae. However, the interferometric data also clearly demonstrate that the submillimeter continuum emission is not centrally concentrated, but rather appears to have a toroidal geometry with substantially lower intensities inside a radius of ~250-300 AU. Spatially resolved emission from the CO J=3-2 transition exhibits a velocity gradient along the major axis of the dust structure. If this kinematic pattern is interpreted as the signature of rotation around a central object, a relatively low mass is inferred (M_star = 0.3 M_sun, with a ~50% uncertainty). We discuss several possible explanations for the observed gas and dust environment around IRAS 04158+2805, including a flattened envelope with an outflow cavity and a large circumbinary ring. This source offers unique views of the gas and dust environment surrounding a young low-mass stellar system. Its properties are generally not commensurate with formation scenarios for such low-mass objects that rely on dynamical ejection, but rather confirms that a single mechanism - molecular cloud core collapse and fragmentation - can produce stars over a wide range of stellar masses (at least an order of magnitude).
[188]  oai:arXiv.org:0804.0437  [pdf] - 11451
The Structure of the DoAr 25 Circumstellar Disk
Comments: ApJL in press
Submitted: 2008-04-02
We present high spatial resolution (< 0.3" = 40$ AU) Submillimeter Array observations of the 865 micron continuum emission from the circumstellar disk around the young star DoAr 25. Despite its bright millimeter emission, this source exhibits only a comparatively small infrared excess and low accretion rate, suggesting that the material and structural properties of the inner disk may be in an advanced state of evolution. A simple model of the physical conditions in the disk is derived from the submillimeter visibilities and the complete spectral energy distribution using a Monte Carlo radiative transfer code. For the standard assumption of a homogeneous grain size distribution at all disk radii, the results indicate a shallow surface density profile, $\Sigma \propto r^{-p}$ with p = 0.34, significantly less steep than a steady-state accretion disk (p = 1) or the often adopted minimum mass solar nebula (p = 1.5). Even though the total mass of material is large (M_d = 0.10 M_sun), the densities inferred in the inner disk for such a model may be too low to facilitate any mode of planet formation. However, alternative models with steeper density gradients (p = 1) can explain the observations equally well if substantial grain growth in the planet formation region (r < 40 AU) has occurred. We discuss these data in the context of such models with dust properties that vary with radius and highlight their implications for understanding disk evolution and the early stages of planet formation.
[189]  oai:arXiv.org:0708.4185  [pdf] - 4404
A Submillimeter View of Circumstellar Dust Disks in $\rho$ Ophiuchus
Comments: accepted in ApJ; 38 pages, 11 figures
Submitted: 2007-08-30
We present new multiwavelength submillimeter continuum measurements of the circumstellar dust around 48 young stars in the $\rho$ Ophiuchus dark clouds. Supplemented with previous 1.3 mm observations of an additional 99 objects from the literature, the statistical distributions of disk masses and submillimeter colors are calculated and compared to those in the Taurus-Auriga region. These basic submillimeter properties of young stellar objects in both environments are shown to be essentially identical. As with their Taurus counterparts, the $\rho$ Oph circumstellar dust properties are shown to evolve along an empirical evolution sequence based on the infrared spectral energy distribution. The combined $\rho$ Oph and Taurus Class II samples (173 sources) are used to set benchmark values for basic outer disk characteristics: M_disk ~ 0.005 solar masses, M_disk/M_star ~ 1%, and $\alpha$ ~ 2 (where $F_{\nu} \propto \nu^{\alpha}$ between 350 microns and 1.3 mm). The precision of these numbers are addressed in the context of substantial solid particle growth in the earliest stages of the planet formation process. There is some circumstantial evidence that disk masses inferred from submillimeter emission may be under-estimated by up to an order of magnitude.
[190]  oai:arXiv.org:astro-ph/0610813  [pdf] - 316584
High Resolution Submillimeter Constraints on Circumstellar Disk Structure
Comments: submitted to ApJ (8/22/06)
Submitted: 2006-10-26
We present a high spatial resolution submillimeter continuum survey of 24 circumstellar disks in the Tau-Aur and Oph-Sco star formation regions using the SMA. In the context of a simple model, we use broadband SEDs and submillimeter visibilities to derive constraints on some basic parameters that describe the structure of these disks. For the typical disk in the sample we infer a radial surface density distribution \Sigma ~ r^-p with a median p ~ 0.5, although consideration of the systematic effects of some of our assumptions suggest that steeper distributions with p ~ 0.7-1.0 are more reasonable. The distribution of the outer radii of these disks shows a distinct peak at R_d = 200 AU, with only a few cases where the disk emission is completely unresolved. Based on these disk structure measurements, the mass accretion rates, and the typical spectral and spatial distributions of submillimeter emission, we show that the observations are in good agreement with similarity solutions for steady accretion disks that have a viscosity parameter alpha ~ 0.01. We provide new estimates of the spectral dependence of the disk opacity with a median spectral index of ~0.7, corrected for optically thick emission. This typical value is consistent with model predictions for the collisional growth of solids to millimeter size scales in the outer disk. Although direct constraints on planet formation in these disks are not currently available, the extrapolated density distributions inferred here are substantially shallower than those calculated based on the solar system or extrasolar planets and typically used in planet formation models. It is possible that we are substantially underestimating disk densities due to an incomplete submillimeter opacity prescription.
[191]  oai:arXiv.org:astro-ph/0609122  [pdf] - 84714
The Dust Properties of Eight Debris Disk Candidates as Determined by Submillimeter Photometry
Comments: accepted by ApJ
Submitted: 2006-09-05
The nature of far-infrared dust emission toward main sequence stars, whether interstellar or circumstellar, can be deduced from submillimeter photometry. We present JCMT/SCUBA flux measurements at 850 microns toward 8 stars with large photospheric excesses at 60-100 microns. 5 sources were detected at 3-sigma or greater significance and one was marginally detected at 2.5-sigma. The inferred dust masses and temperatures range from 0.033 to 0.24 Earth masses and 43-65 K respectively. The frequency behavior of the opacity, tau_nu ~ nu^beta, is relatively shallow, beta < 1. These dust properties are characteristic of circumstellar material, most likely the debris from planetesimal collisions. The 2 non-detections have lower temperatures, 35-38 K and steeper opacity indices, beta > 1.5, that are more typical of interstellar cirrus. The confirmed disks all have inferred diameters > 2'', most lie near the upper envelope of the debris disk mass distribution, and 4 are bright enough to be feasible for high resolution imaging.
[192]  oai:arXiv.org:astro-ph/0511589  [pdf] - 77947
Submillimeter studies of circumstellar disks in Taurus and Orion
Comments: Conference proceedings for "Revealing the Molecular Universe: One Telescope is Never Enough"; informal writeup of recently published work for a meeting celebrating the career of Jack Welch
Submitted: 2005-11-18
We highlight two recent studies of circumstellar disks in the Taurus and Orion star forming regions. Using the JCMT and CSO, we measure disk fluxes in Taurus over a wide range of submillimeter wavelengths and determine the frequency dependence of the dust opacity. We find clear evidence for a systematic change in its behavior with time, most readily explained by grain growth. Using the SMA, we observed the protoplanetary disks (proplyds) in the Orion Trapezium cluster. The combination of high resolution, high frequency, and high sensitivity that this instrument provides made it possible to resolve disks from one another, distinguish their emission from background cloud material and surrounding ionized gas, and to detect thermal dust emission. This allowed us to make the first mass measurements of the proplyds and to assess their viability for planet formation.
[193]  oai:arXiv.org:astro-ph/0506225  [pdf] - 73625
The Masses of the Orion Proplyds from Submillimeter Dust Emission
Comments: 15 pages, 3 figures, accepted for publication in the Astrophysical Journal
Submitted: 2005-06-09, last modified: 2005-07-01
We have imaged the 880 micron continuum emission from the "proplyds" in the center of the Trapezium Cluster in Orion using the Submillimeter Array with a beam size 1.5'' FWHM and an rms of 2.7 mJy. Five sources are detected with fluxes in the range 18 to 38 mJy, which includes dust emission from four proplyds and ionized gas from theta1 Ori G. The total masses of the detected proplyds derived from their dust emission range from 0.013 to 0.024 Msun assuming a dust temperature of 20 K and mass opacity of 0.03 cm^2/g. The eighteen other proplyds within the field-of-view were not formally detected individually, but by combining the fluxes measured at their locations, an average flux of 1.1mJy was determined for them on a statistical basis, corresponding to a mass of 8e-4 Msun. The four detected proplyds have sufficient disk mass bound to their central stars to form planetary systems on the scale of our Solar System.
[194]  oai:arXiv.org:astro-ph/0506187  [pdf] - 73587
Circumstellar Dust Disks in Taurus-Auriga: The Submillimeter Perspective
Comments: accepted by ApJ
Submitted: 2005-06-08
We present a sensitive, multiwavelength submillimeter continuum survey of 153 young stellar objects in the Taurus-Auriga star formation region. The submillimeter detection rate is 61% to a completeness limit of ~10 mJy (3-sigma) at 850 microns. The inferred circumstellar disk masses are log-normally distributed with a mean mass of ~0.005 solar masses and a large dispersion (0.5 dex). Roughly one third of the submillimeter sources have disk masses larger than the minimal nebula from which the solar system formed. The median disk to star mass ratio is 0.5%. The empirical behavior of the submillimeter continuum is best described as F_nu ~ nu^(2.0 +/- 0.5) between 350 microns and 1.3 mm, which we argue is due to the combined effects of the fraction of optically thick emission and a flatter frequency behavior of the opacity compared to the ISM. This latter effect could be due to a substantial population of large dust grains, which presumably would have grown through collisional agglomeration. In this sample, the only stellar property that is correlated with the outer disk is the presence of a companion. We find evidence for significant decreases in submillimeter flux densities, disk masses, and submillimeter continuum slopes along the canonical infrared spectral energy distribution evolution sequence for young stellar objects. The fraction of objects detected in the submillimeter is essentially identical to the fraction with excess near-infrared emission, suggesting that dust in the inner and outer disk are removed nearly simultaneously.
[195]  oai:arXiv.org:astro-ph/0411131  [pdf] - 68730
Submillimeter Array Observations of Disks in the SR 24 Multiple Star System
Comments: 10 pages, 2 figures, to appear in ApJ Letters
Submitted: 2004-11-05
We present high-resolution aperture synthesis images from the Submillimeter Array of the 225 GHz (1.3 mm) continuum and CO J=2-1 line emission from the disks around components of the hierarchical triple system SR 24, located in the Ophiuchus star-forming region. The most widely separated component, SR 24 S (with a projected semimajor axis of 832 AU), has a circumstellar disk with properties typical of those around single T Tauri stars. The binary SR 24 N (with a projected semimajor axis of 32 AU) is undetected in the continuum, but has strong, resolved CO emission which likely originates in a circumbinary disk with a central gap. The data constrain the total disk mass in the SR 24 N system to be < 10^-3 M_sun and indicate that the depletion of CO onto dust grains is not more than 100 times larger than the mean value in the interstellar medium. The SR 24 N disk is unusual in that it is only detected in line emission. It is possible that other low mass disks around binaries and single stars may have been missed in single-dish continuum surveys.
[196]  oai:arXiv.org:astro-ph/0406089  [pdf] - 65246
Mid-Infrared and Submillimeter Observations of the Illuminating Source of McNeil's Variable Nebula
Comments: submitted to ApJL
Submitted: 2004-06-03
We present post-outburst observations of the mid-infrared spectrum and submillimeter continuum of the illuminating source of the newly-discovered McNeil's Nebula in the L1630 region of Orion. The 12 micron flux of this source has increased by a factor of 25 after the outburst, whereas the submillimeter continuum remains at its pre-outburst level. The bolometric luminosity has increased by at least an order of magnitude, to ~34 L_sun, and is likely less than 90 L_sun. The mid-infrared spectrum exhibits a strong and red continuum with no emission or absorption features. The infrared slope of the spectral energy distribution characterizes the illuminating source as a flat-spectrum protostar, in both its active and quiescent states. New CO spectral line observations show no evidence for a molecular outflow.
[197]  oai:arXiv.org:astro-ph/0404260  [pdf] - 64168
X-Ray and Infrared Observations of Embedded Young Stars in L1630
Comments: to appear in ApJ. 32 pages, 8 figures. High-resolution figures available at http://www.ifa.hawaii.edu/~andrews/L1630.ps
Submitted: 2004-04-13
The HH 24-26 star-forming region within the L1630 dark cloud in Orion contains a remarkable collection of rare Class 0 and Class I protostars, collimated molecular and ionized jets, and a luminous but spatially unresolved ASCA X-ray source. To study the X-ray properties of the embedded protostar population of that region, we have obtained a deep X-ray image with the ACIS-S camera on board the Chandra X-Ray Observatory. A number of H-alpha emission-line objects were detected in the areas surrounding HH 24-26, of which the weak-line T Tauri star SSV 61 was the brightest source, at a steady luminosity of 10^31.9 erg/s (0.3-10 keV). Two Class I protostars aligned with optical jets in HH 24, SSV 63E and SSV 63W, were also detected, as was the continuum radio source SSV 63NE, which is very likely an extreme Class I or Class 0 object. We observed no X rays from the Class 0 protostars HH 24-MMS and HH 25-MMS, nor any from regions of the cloud bounded by HH 25 and HH 26, at a 2-sigma upper limit of 10^30.0 erg/s. HH 26-IR, the Class I object thought to be the origin of the HH 26 flow, was not detected. Near-infrared spectroscopy obtained at the NASA IRTF reveals 3 micron ice bands in the spectra of SSV 59, 63E, 63W, and HH 26-IR, and 2.3 micron CO overtone absorption bands for SSV 61. SSV 60, which lies astride one end of the great arc of nebulosity forming HH 25, exhibits a deep infrared ice band and CO absorption, but is not an X-ray source, and is most likely a distant background giant of late spectral type.
[198]  oai:arXiv.org:astro-ph/0104028  [pdf] - 41767
Small-Scale Interstellar Na I Structure Toward M92
Comments: 9 pages, 2 figures, accepted for publication in ApJ Letters
Submitted: 2001-04-02
We have used integral field echelle spectroscopy with the DensePak fiber-optic array on the KPNO WIYN telescope to observe the central 27" x 43" of the globular cluster M92 in the Na I D wavelength region at a spatial resolution of 4". Two interstellar Na I absorption components are evident in the spectra at LSR velocities of 0 km/s (Cloud 1) and -19 km/s (Cloud 2). Substantial strength variations in both components are apparent down to scales limited by the fiber-to-fiber separations. The derived Na I column densities differ by a factor of 4 across the Cloud 1 absorption map and by a factor of 7 across the Cloud 2 map. Using distance upper limits of 400 and 800 pc for Cloud 1 and Cloud 2, respectively, the absorption maps indicate structure in the ISM down to scales of 1600 and 3200 AU. The fiber-to-fiber Na I column density differences toward M92 are comparable to those found in a similar study of the ISM toward the globular cluster M15. Overall, the structures in the interstellar components toward M92 have significantly lower column densities than those toward M15. We interpret these low column density structures as small-scale turbulent variations in the gas and compare them to the larger-scale, higher column density variations toward M15, which may be the hallmarks of actual H I structures.