Source author record

John M. Carpenter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

51 published item(s)

preprint2026arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

New Constraints on Protoplanetary Disk Gas Masses in Lupus

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

preprint2022arXiv

Update on the Systematics in the ALMA Proposal Review Process after Cycle 8

We present an updated analysis of systematics in the Atacama Large Millimeter/submillimeter Array (ALMA) proposal ranks from Carpenter (2020) to include the last two ALMA cycles, when significant changes were introduced in the proposal review process. In Cycle 7, the investigator list on the proposal cover sheet was randomized such that the reviewers were aware of the overall proposal team but did not know the identity of the principal investigator (PI). In Cycle 8, ALMA adopted distributed peer review for most proposals and implemented dual-anonymous review for all proposals, in which the identity of the proposal team was not revealed to the reviewers. The most significant change in the systematics in Cycles 7 and 8 compared to previous cycles is related to the experience of PIs in submitting ALMA proposals. PIs that submit a proposal every cycle tend to have ranks that are consistent with average in Cycles 7 and 8 whereas previously they had the best overall ranks. Also, PIs who submitted a proposal for the second time show improved ranks over previous cycles. These results suggest some biases related to the relative prominence of the PI have been present in the ALMA review process. Systematics related to regional affiliation remain largely unchanged in that PIs from Chile, East Asia, and non-ALMA regions tend to have poorer overall ranks than PIs from Europe and North America. The systematics of how one region ranks proposals from another region are also investigated. No significant differences in the overall ranks based on gender of the PI are observed.

preprint2020arXiv

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

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

preprint2020arXiv

ALMA Band-to-band Phase Referencing: Imaging Capabilities on Long Baselines and High Frequencies

High-frequency long-baseline experiments with the Atacama Large Millimeter/submillimeter Array were organized to test the high angular resolution imaging capabilities in the submillimeter wave regime using baselines up to 16 km. Four experiments were conducted, two Band 7 (289 GHz) and two Band 8 (405 GHz) observations. Phase correction using band-to-band (B2B) phase referencing was used with a phase calibrator only 0.7deg away observed in Band 3 (96 GHz) and Band 4 (135 GHz), respectively. In Band 8, we achieved the highest resolution of 14x11 mas. We compared the synthesis images of the target quasar using 20 and 60 s switching cycle times in the phase referencing. In Band 7, the atmosphere had good stability in phase rms (<0.5 rad over 2 minutes), and there was little difference in image coherence between the 20 and 60 s switching cycle times. One Band 8 experiment was conducted under a worse phase rms condition (>1 rad over 2 minutes), which led to a significantly reduced coherence when using the 60 s switching cycle time. One of our four experiments indicates that the residual phase rms error after phase referencing can be reduced to 0.16 rad at 289 GHz in using the 20 s switching cycle time. Such conditions would meet the phase correction requirement of image coherence of >70% in Band 10, assuming a similar phase calibrator separation angle, emphasizing the need for such B2B phase referencing observing at high frequencies.

preprint2020arXiv

ALMA High-frequency Long Baseline Campaign in 2017: Band-to-band Phase Referencing in Submillimeter Waves

In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequency. For B2B phase referencing, it is expected that a closer phase calibrator to a target can be used, comparing to standard in-band phase referencing. In the first step, it is ensured that an instrumental phase offset difference between low- and high-frequency Bands can be removed using a differential gain calibration in which a phase calibrator is certainly detected while frequency switching. In the next step, comparative experiments are arranged to investigate the image quality between B2B and in-band phase referencing with phase calibrators at various separation angles. In the final step, we conducted long baseline imaging tests for a quasar at 289 GHz in Band 7 and 405 GHz in Band 8 and complex structure sources of HL Tau and VY CMa at ~670 GHz in Band 9. The B2B phase referencing was successfully applied, allowing us to achieve an angular resolution of 14x11 and 10x8 mas for HL Tau and VY CMa, respectively. There is a high probability of finding a low-frequency calibrator within 5.4 deg in B2B phase referencing, bright enough to use an 8 s scan length combined with a 7.5 GHz bandwidth.

preprint2020arXiv

ALMA High-frequency Long-baseline Campaign in 2017: A Comparison of the Band-to-band and In-band Phase Calibration Techniques and Phase-calibrator Separation Angles

The Atacama Large millimeter/submillimeter Array (ALMA) obtains spatial resolutions of 15 to 5 milli-arcsecond (mas) at 275-950GHz (0.87-0.32mm) with 16km baselines. Calibration at higher-frequencies is challenging as ALMA sensitivity and quasar density decrease. The Band-to-Band (B2B) technique observes a detectable quasar at lower frequency that is closer to the target, compared to one at the target high-frequency. Calibration involves a nearly constant instrumental phase offset between the frequencies and the conversion of the temporal phases to the target frequency. The instrumental offsets are solved with a differential-gain-calibration (DGC) sequence, consisting of alternating low and high frequency scans of strong quasar. Here we compare B2B and in-band phase referencing for high-frequencies ($>$289GHz) using 2-15km baselines and calibrator separation angles between $\sim$0.68 and $\sim$11.65$^{\circ}$. The analysis shows that: (1) DGC for B2B produces a coherence loss $<$7% for DGC phase RMS residuals $<$30$^{\circ}$. (2) B2B images using close calibrators ( $<$1.67$^{\circ}$ ) are superior to in-band images using distant ones ( $>$2.42$^{\circ}$ ). (3) For more distant calibrators, B2B is preferred if it provides a calibrator $\sim$2$^{\circ}$ closer than the best in-band calibrator. (4) Decreasing image coherence and poorer image quality occur with increasing phase calibrator separation angle because of uncertainties in the antenna positions and sub-optimal phase referencing. (5) To achieve $>$70% coherence for long-baseline (16 km) band 7 (289GHz) observations, calibrators should be within $\sim$4$^{\circ}$ of the target.

preprint2020arXiv

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

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

preprint2020arXiv

An evolutionary study of volatile chemistry in protoplanetary disks

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

preprint2020arXiv

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

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

preprint2020arXiv

The CARMA-NRO Orion Survey: Filament Formation via Collision-Induced Magnetic Reconnection -- The Stick in Orion A

A unique filament is identified in the {\it Herschel} maps of the Orion A giant molecular cloud. The filament, which, we name the Stick, is ruler-straight and at an early evolutionary stage. Transverse position-velocity diagrams show two velocity components closing in on the Stick. The filament shows consecutive rings/forks in C$^{18}$O(1-0) channel maps, which is reminiscent of structures generated by magnetic reconnection. We propose that the Stick formed via collision-induced magnetic reconnection (CMR). We use the magnetohydrodynamics (MHD) code Athena++ to simulate the collision between two diffuse molecular clumps, each carrying an anti-parallel magnetic field. The clump collision produces a narrow, straight, dense filament with a factor of $>$200 increase in density. The production of the dense gas is seven times faster than free-fall collapse. The dense filament shows ring/fork-like structures in radiative transfer maps. Cores in the filament are confined by surface magnetic pressure. CMR can be an important dense-gas-producing mechanism in the Galaxy and beyond.

preprint2020arXiv

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

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

preprint2020arXiv

The Formation of a Stellar Association in the NGC 7000/IC 5070 Complex: Results from Kinematic Analysis of Stars and Gas

We examine the clustering and kinematics of young stellar objects (YSOs) in the North America/Pelican Nebulae, as revealed by Gaia astrometry, in relation to the structure and motions of the molecular gas, as indicated in molecular line maps. The Gaia parallaxes and proper motions allow us to significantly refine previously published lists of YSOs, demonstrating that many of the objects previously thought to form a distributed population turn out to be non-members. The members are subdivided into at least 6 spatio-kinematic groups, each of which is associated with its own molecular cloud component or components. Three of the groups are expanding, with velocity gradients of 0.3-0.5 km/s/pc, up to maximum velocities of ~8 km/s away from the groups' centers. The two known O-type stars associated with the region, 2MASS J20555125+4352246 and HD 199579, are rapidly escaping one of these groups, following the same position-velocity relation as the low-mass stars. We calculate that a combination of gas expulsion and tidal forces from the clumpy distribution of molecular gas could impart the observed velocity gradients within the groups. However, on a global scale, the relative motions of the groups do not appear either divergent or convergent. The velocity dispersion of the whole system is consistent with the kinetic energy gained due to gravitational collapse of the complex. Most of the stellar population has ages similar to the free-fall timescales for the natal clouds. Thus, we suggest the nearly free-fall collapse of a turbulent molecular cloud as the most likely scenario for star formation in this complex.

preprint2019arXiv

Systematics in the ALMA Proposal Review Rankings

The results from the ALMA proposal peer review process in Cycles 0-6 are analyzed to identify any systematics in the scientific rankings that may signify bias. Proposal rankings are analyzed with respect to the experience level of a Principal Investigator (PI) in submitting ALMA proposals, regional affiliation (Chile, East Asia, Europe, North America, or Other), and gender. The analysis was conducted for both the Stage 1 rankings, which are based on the preliminary scores from the reviewers, and the Stage 2 rankings, which are based on the final scores from the reviewers after participating in a face-to-face panel discussion. Analysis of the Stage 1 results shows that PIs who submit an ALMA proposal in multiple cycles have systematically better proposal ranks than PIs who have submitted proposals for the first time. In terms of regional affiliation, PIs from Europe and North America have better Stage 1 rankings than PIs from Chile and East Asia. Consistent with Lonsdale et al. (2016), proposals led by men have better Stage 1 rankings than women when averaged over all cycles. This trend was most noticeably present in Cycle 3, but no discernible differences in the Stage 1 rankings are present in recent cycles. Nonetheless, in each cycle to date, women have had a lower proposal acceptance rate than men even after differences in demographics are considered. Comparison of the Stage 1 and Stage 2 rankings reveal no significant changes in the distribution of proposal ranks by experience level, regional affiliation, or gender as a result of the panel discussions, although the proposal ranks for East Asian PIs show a marginally significant improvement from Stage 1 to Stage 2 when averaged over all cycles. Thus any systematics in the proposal rankings are introduced primarily in the Stage 1 process and not from the face-to-face discussions.

preprint2016arXiv

A Neptune-sized transiting planet closely orbiting a 5-10-million-year-old star

Theories of the formation and early evolution of planetary systems postulate that planets are born in circumstellar disks, and undergo radial migration during and after dissipation of the dust and gas disk from which they formed. The precise ages of meteorites indicate that planetesimals - the building blocks of planets - are produced within the first million years of a star's life. A prominent question is: how early can one find fully formed planets like those frequently detected on short orbital periods around mature stars? Some theories suggest the in situ formation of planets close to their host stars is unlikely and the existence of such planets is evidence for large scale migration. Other theories posit that planet assembly at small orbital separations may be common. Here we report on a newly-born, transiting planet orbiting its star every 5.4 days. The planet is 50 per cent larger than Neptune, and its mass is less than 3.6 times Jupiter (at 99.7 per cent confidence), with a true mass likely to be within a factor of several of Neptune's. The 5-10 million year old star has a tenuous dust disk extending outwards from about 2 times the Earth-Sun separation, in addition to the large planet located at less than one-twentieth the Earth-Sun separation.

preprint2016arXiv

ALMA Observations of Circumstellar Disks in the Upper Scorpius OB Association

We present ALMA observations of 106 G-, K-, and M-type stars in the Upper Scorpius OB Association hosting circumstellar disks. With these data, we measure the 0.88 mm continuum and $^{12}$CO $J$ = 3$-$2 line fluxes of disks around low-mass ($0.14-1.66$ $M_{\odot}$) stars at an age of 5-11 Myr. Of the 75 primordial disks in the sample, 53 are detected in the dust continuum and 26 in CO. Of the 31 disks classified as debris/evolved transitional disks, 5 are detected in the continuum and none in CO. The lack of CO emission in approximately half of the disks with detected continuum emission can be explained if CO is optically thick but has a compact emitting area ($\lesssim 40$ au), or if the CO is heavily depleted by a factor of at least $\sim1000$ relative to interstellar medium abundances and is optically thin. The continuum measurements are used to estimate the dust mass of the disks. We find a correlation between disk dust mass and stellar host mass consistent with a power-law relation of $M_{\mathrm{dust}}\propto M_*^{1.67\pm0.37}$. Disk dust masses in Upper Sco are compared to those measured in the younger Taurus star-forming region to constrain the evolution of disk dust mass. We find that the difference in the mean of $\log(M_{\mathrm{dust}}/M_*)$ between Taurus and Upper Sco is $0.64\pm0.09$, such that $M_{\mathrm{dust}}/M_*$ is lower in Upper Sco by a factor of $\sim4.5$.

preprint2016arXiv

Debris Disks in the Scorpius-Centaurus OB Association Resolved by ALMA

We present a CO(2-1) and 1240 um continuum survey of 23 debris disks with spectral types B9-G1, observed at an angular resolution of 0.5-1 arcsec with the Atacama Large Millimeter/Submillimeter Array (ALMA). The sample was selected for large infrared excess and age ~10 Myr, to characterize the prevalence of molecular gas emission in young debris disks. We identify three CO-rich debris disks, plus two additional tentative (3-sigma) CO detections. Twenty disks were detected in the continuum at the >3-sigma level. For the 12 disks in the sample that are spatially resolved by our observations, we perform an independent analysis of the interferometric continuum visibilities to constrain the basic dust disk geometry, as well as a simultaneous analysis of the visibilities and broad-band spectral energy distribution to constrain the characteristic grain size and disk mass. The gas-rich debris disks exhibit preferentially larger outer radii in their dust disks, and a higher prevalence of characteristic grain sizes smaller than the blowout size. The gas-rich disks do not exhibit preferentially larger dust masses, contrary to expectations for a scenario in which a higher cometary destruction rate would be expected to result in a larger mass of both CO and dust. The three debris disks in our sample with strong CO detections are all around A stars: the conditions in disks around intermediate-mass stars appear to be the most conducive to the survival or formation of CO.

preprint2016arXiv

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

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

preprint2016arXiv

Spiral Density Waves in a Young Protoplanetary Disk

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

preprint2015arXiv

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

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

preprint2015arXiv

Discovery of Seven Companions to Intermediate Mass Stars with Extreme Mass Ratios in the Scorpius-Centaurus Association

We report the detection of seven low mass companions to intermediate-mass stars (SpT B/A/F; $M$$\approx$1.5-4.5 solar masses) in the Scorpius-Centaurus Association using nonredundant aperture masking interferometry. Our newly detected objects have contrasts $ΔL'$$\approx$4-6, corresponding to masses as low as $\sim$20 Jupiter masses and mass ratios of $q$$\approx$0.01-0.08, depending on the assumed age of the target stars. With projected separations $ρ$$\approx$10-30 AU, our aperture masking detections sample an orbital region previously unprobed by conventional adaptive optics imaging of intermediate mass Scorpius-Centaurus stars covering much larger orbital radii ($\approx$30-3000 AU). At such orbital separations, these objects resemble higher mass versions of the directly imaged planetary mass companions to the 10-30 Myr, intermediate-mass stars HR 8799, $β$ Pictoris, and HD95086. These newly discovered companions span the brown dwarf desert, and their masses and orbital radii provide a new constraint on models of the formation of low-mass stellar and substellar companions to intermediate-mass stars.

preprint2015arXiv

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

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

preprint2015arXiv

K2 Discovery of Young Eclipsing Binaries in Upper Scorpius: Direct Mass and Radius Determinations for the Lowest Mass Stars and Initial Characterization of an Eclipsing Brown Dwarf Binary

We report the discovery of three low-mass double-lined eclipsing binaries in the pre-main sequence Upper Scorpius association, revealed by $K2$ photometric monitoring of the region over $\sim$ 78 days. The orbital periods of all three systems are $<$5 days. We use the $K2$ photometry plus multiple Keck/HIRES radial velocities and spectroscopic flux ratios to determine fundamental stellar parameters for both the primary and secondary components of each system, along with the orbital parameters. We present tentative evidence that EPIC 203868608 is a hierarchical triple system comprised of an eclipsing pair of $\sim$25 $M_\mathrm{Jup}$ brown dwarfs with a wide M-type companion. If confirmed, it would constitute only the second double-lined eclipsing brown dwarf binary system discovered to date. The double-lined system EPIC 203710387 is composed of nearly identical M4.5-M5 stars with fundamentally determined masses and radii measured to better than 3% precision ($M_1=0.1183\pm0.0028 M_\odot$, $M_2=0.1076\pm0.0031 M_\odot$ and $R_1=0.417\pm0.010 R_\odot$, $R_2=0.450\pm0.012 R_\odot$) from combination of the light curve and radial velocity time series. These stars have the lowest masses of any stellar mass double-lined eclipsing binary to date. Finally, EPIC 203476597 is a compact single-lined system with a G8-KO primary and a likely mid-K secondary whose line are revealed in spectral ratios. Continued measurement of radial velocities and spectroscopic flux ratios will better constrain fundamental parameters and should elevate the objects to benchmark status. We also present revised parameters for the double-lined eclipsing binary UScoCTIO 5 ($M_1=0.3336\pm0.0022 M_\odot$, $M_2=0.3200\pm0.0022 M_\odot$ and $R_1=0.862\pm0.012$, $R_2=0.852\pm0.013 R_\odot$). We discuss the implications of our results on these $\sim$0.1-1.5 $M_\odot$ stars for pre-main-sequence evolutionary models.

preprint2015arXiv

Protoplanetary Disk Masses in the Young NGC 2024 Cluster

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

preprint2014arXiv

An ALMA Continuum Survey of Circumstellar Disks in the Upper Scorpius OB Association

We present ALMA 880 micron continuum observations of 20 K and M-type stars in the Upper Scorpius OB association that are surrounded by protoplanetary disks. These data are used to measure the dust content in disks around low mass stars (0.1-1.6 Msun) at a stellar age of 5-11 Myr. Thirteen sources were detected in the 880 micron dust continuum at >3 sigma with inferred dust masses between 0.3 and 52 Mearth. The dust masses tend to be higher around the more massive stars, but the significance is marginal in that the probability of no correlation is p ~ 0.03. The evolution in the dust content in disks was assessed by comparing the Upper Sco observations with published continuum measurements of disks around ~ 1-2 Myr stars in the Class II stage in the Taurus molecular cloud. While the dust masses in the Upper Sco disks are on average lower than in Taurus, any difference in the dust mass distributions is significant at less than 3sigma. For stellar masses between 0.49 Msun and 1.6 Msun, the mean dust mass in disks is lower in Upper Sco relative to Taurus by Delta log Mdust = 0.44 +/-0.26.

preprint2014arXiv

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

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

preprint2014arXiv

Large-Scale Asymmetries in the Transitional Disks of SAO 206462 and SR 21

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations in the dust continuum (690 GHz, 0.45 mm) and 12CO J=6-5 spectral line emission, of the transitional disks surrounding the stars SAO 206462 and SR 21. These ALMA observations resolve the dust-depleted disk cavities and extended gaseous disks, revealing large-scale asymmetries in the dust emission of both disks. We modeled these disks structures with a ring and an azimuthal gaussian, where the azimuthal gaussian is motivated by the steady-state vortex solution from Lyra & Lin (2013). Compared to recent observations of HD 142527, Oph IRS 48, and LkHa 330, these are low-contrast (< 2) asymmetries. Nevertheless, a ring alone is not a good fit, and the addition of a vortex prescription describes these data much better. The asymmetric component encompasses 15% and 28% of the total disk emission in SAO 206462 and SR 21 respectively, which corresponds to a lower limit of 2 MJup of material within the asymmetry for both disks. Although the contrast in the dust asymmetry is low, we find that the turbulent velocity inside it must be large (~20% of the sound speed) in order to drive these azimuthally wide and radially narrow vortex-like structures. We obtain residuals from the ring and vortex fitting that are still significant, tracing non-axisymmetric emission in both disks. We compared these submillimeter observations with recently published H-band scattered light observations. For SR 21 the scattered light emission is distributed quite differently from the submillimeter continuum emission, while for SAO 206462 the submillimeter residuals are suggestive of spiral-like structure similar to the near-IR emission.

preprint2014arXiv

Resolved Multifrequency Radio Observations of GG Tau

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

preprint2014arXiv

Searching for circumplanetary disks around LkCa 15

We present Karl G. Jansky Very Large Array (VLA) observations of the 7 mm continuum emission from the disk surrounding the young star LkCa 15. The observations achieve an angular resolution of 70 mas and spatially resolve the circumstellar emission on a spatial scale of 9 AU. The continuum emission traces a dusty annulus of 45 AU in radius that is consistent with the dust morphology observed at shorter wavelengths. The VLA observations also reveal a compact source at the center of the disk, possibly due to thermal emission from hot dust or ionized gas located within a few AU from the central star. No emission is observed between the star and the dusty ring, and, in particular, at the position of the candidate protoplanet LkCa 15 b. By comparing the observations with theoretical models for circumplanetary disk emission, we find that if LkCa~15~b is a massive planet (>5 M_J) accreting at a rate greater than 1.e-6 M_J yr^{-1}, then its circumplanetary disk is less massive than 0.1 M_J, or smaller than 0.4 Hill radii. Similar constraints are derived for any possible circumplanetary disk orbiting within 45 AU from the central star. The mass estimate are uncertain by at least one order of magnitude due to the uncertainties on the mass opacity. Future ALMA observations of this system might be able to detect circumplanetary disks down to a mass of 5.e-4 M_J and as small as 0.2 AU, providing crucial constraints on the presence of giant planets in the act of forming around this young star.

preprint2014arXiv

TADPOL: A 1.3 mm Survey of Dust Polarization in Star-forming Cores and Regions

We present λ1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

preprint2014arXiv

The CARMA Paired Antenna Calibration System: Atmospheric Phase Correction for Millimeter Wave Interferometry and its Application to Mapping the Ultraluminous Galaxy Arp 193

Phase fluctuations introduced by the atmosphere are the main limiting factor in attaining diffraction limited performance in extended interferometric arrays at millimeter and submillimeter wavelengths. We report the results of C-PACS, the Combined Array for Research in Millimeter-Wave Astronomy Paired Antenna Calibration System. We present a systematic study of several hundred test observations taken during the 2009-2010 winter observing season where we utilize CARMA's eight 3.5-m antennas to monitor an atmospheric calibrator while simultaneously acquiring science observations with 6.1-m and 10.4-m antennas on baselines ranging from a few hundred meters to ~2 km. We find that C-PACS is systematically successful at improving coherence on long baselines under a variety of atmospheric conditions. We find that the angular separation between the atmospheric calibrator and target source is the most important consideration, with consistently successful phase correction at CARMA requiring a suitable calibrator located $\lesssim$6$^\circ$ away from the science target. We show that cloud cover does not affect the success of C-PACS. We demonstrate C-PACS in typical use by applying it to the observations of the nearby very luminous infrared galaxy Arp 193 in $^{12}$CO(2-1) at a linear resolution of ~70 pc (0.12" x 0.18"), 3 times better than previously published molecular maps of this galaxy. We resolve the molecular disk rotation kinematics and the molecular gas distribution and measure the gas surface densities and masses on 90 pc scales. We find that molecular gas constitutes $\sim30\%$ of the dynamical mass in the inner 700 pc of this object with a surface density $\sim10^4 M_\odot$ pc$^{-2}$; we compare these properties to those of the starburst region of NGC 253.

preprint2013arXiv

A joint analysis of AMI and CARMA observations of the recently discovered SZ galaxy cluster system AMI-CL J0300+2613

We present CARMA observations of a massive galaxy cluster discovered in the AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian analysis of the AMI, CARMA and joint AMI & CARMA uv-data is used to quantify the detection significance and parameterise both the physical and observational properties of the cluster whilst accounting for the statistics of primary CMB anisotropies, receiver noise and radio sources. The joint analysis of the AMI & CARMA uv-data was performed with two parametric physical cluster models: the β-model; and the model described in Olamaie et al. 2012 with the pressure profile fixed according to Arnaud et al. 2010. The cluster mass derived from these different models is comparable but our Bayesian evidences indicate a preference for the β-profile which we, therefore, use throughout our analysis. From the CARMA data alone we obtain a Bayesian probability of detection ratio of 12.8:1 when assuming that a cluster exists within our search area; alternatively assuming that Jenkins et al. 2001 accurately predicts the number of clusters as a function of mass and redshift, the Bayesian probability of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis of the AMI & CARMA data with a physical cluster model we derive the total mass internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological β-model to quantify the temperature decrement as a function of angular distance we find a central SZ temperature decrement of 170μK in the AMI & CARMA data. The SZ decrement in the CARMA data is weaker than expected and we speculate that this is a consequence of the cluster morphology. In a forthcoming study we will assess the impact of cluster morphology on the SZ decrements that are observed with interferometers such as AMI and CARMA.

preprint2013arXiv

An azimuthal asymmetry in the LkHa 330 disk

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.

preprint2013arXiv

Misalignment of Magnetic Fields and Outflows in Protostellar Cores

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

preprint2012arXiv

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

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

preprint2012arXiv

Highly Variable Extinction and Accretion in the Jet-driving Class I Type Young Star PTF 10nvg (V2492 Cyg, IRAS 20496+4354)

We report extensive new photometry and spectroscopy of the highly variable young stellar object PTF 10nvg including optical and near-infrared time series data as well as mid-infrared and millimeter data. Following the previously reported 2010 rise, during 2011 and 2012 the source underwent additional episodes of brightening and dimming events including prolonged faint states. The observed high-amplitude variations are largely consistent with extinction changes having a 220 day quasi-periodic signal. Spectral evolution includes not only changes in the spectral slope but correlated variation in the prominence of TiO/VO/CO bands and atomic line emission, as well as anticorrelated variation in forbidden line emission which, along with H_2, dominates optical and infrared spectra at faint epochs. Neutral and singly-ionized atomic species are likely formed in an accretion flow and/or impact while the origin of zero-velocity atomic LiI 6707 in emission is unknown. Forbidden lines, including several rare species, exhibit blueshifted emission profiles and likely arise from an outflow/jet. Several of these lines are also seen spatially offset from the continuum source position, presumably in a shocked region of an extended jet. CARMA maps resolve on larger scales a spatially extended outflow in mm-wavelength CO. We attribute the observed photometric and spectroscopic behavior in terms of occultation of the central star as well as the bright inner disk and the accretion/outflow zones that renders shocked gas in the inner part of the jet amenable to observation at the faint epochs. We discuss PTF 10nvg as a source exhibiting both accretion-driven (perhaps analogous to V1647 Ori) and extinction-driven (perhaps analogous to UX Ori or GM Cep) high-amplitude variability phenomena.

preprint2012arXiv

On the nature of the transition disk around LkCa 15

We present CARMA 1.3 mm continuum observations of the T Tauri star LkCa 15,which resolve the circumstellar dust continuum emission on angular scales between 0.2-3 arcsec, corresponding to 28-420 AU at the distance of the star. The observations resolve the inner gap in the dust emission and reveal an asymmetric dust distribution in the outer disk. (Abridge) We calculate that 90% of the dust emission arises from an azimuthally symmetric ring that contains about 5x10^{-4} M_sun of dust. A low surface-brightness tail that extends to the northwest out to a radius of about 300 AU contains the remaining 10% of the observed continuum emission. The ring is modeled with a rather flat surface density profile between 40 and 120 AU, while the inner cavity is consistent with either a sharp drop of the 1.3 mm dust optical depth at about 42 AU or a smooth inward decrease between 3 and 85 AU. (Abridge). Within 40 AU, the observations constrain the amount of dust between 10^{-6} and 7 Earth masses, where the minimum and maximum limits are set by the near-IR SED modeling and by the mm-wave observations of the dust emission respectively. In addition, we confirm the discrepancy in the outer disk radius inferred from the dust and gas, which corresponds to 150 AU and 900 AU respectively. We cannot reconcile this difference by adopting an exponentially tapered surface density profile as suggested for other systems, but we instead suggest that the gas surface density in the outer disk decreases less steeply than that predicted by model fits to the dust continuum emission. The lack of continuum emission at radii lager than 120 AU suggests a drop of at least a factor of 5 in the dust-to-gas ratio, or in the dust opacity. We show that a sharp dust opacity drop of this magnitude is consistent with a radial variation of the grain size distribution as predicted by existing grain growth models.

preprint2012arXiv

The Dynamical Mass and Three-Dimensional Orbit of HR7672B: A Benchmark Brown Dwarf with High Eccentricity

The companion to the G0V star HR7672 directly imaged by Liu et al. (2002) has moved measurably along its orbit since the discovery epoch, making it possible to determine its dynamical properties. Originally targeted with adaptive optics because it showed a long-term radial velocity acceleration (trend), we have monitored this star with precise Doppler measurements and have now established a 24 year time baseline. The radial velocity variations show significant curvature (change in the acceleration) including an inflection point. We have also obtained a recent image of HR7672B with NIRC2 at Keck. The astrometry also shows curvature. In this paper, we use jointly-fitted Doppler and astrometric models to calculate the three-dimensional orbit and dynamical mass of the companion. The mass of the host star is determined using a direct radius measurement from CHARA interferometry in combination with high resolution spectroscopic modeling. We find that HR7672B has a highly eccentric, $e=0.50^{+0.01}_{-0.01}$, near edge-on, $i=97.3^{+0.4}_{-0.5}$ deg, orbit with semimajor axis, $a=18.3^{+0.4}_{-0.5}$ AU. The mass of the companion is $m=68.7^{+2.4}_{-3.1}M_J$ at the 68.2% confidence level. HR7672B thus resides near the substellar boundary, just below the hydrogen-fusing limit. These measurements of the companion mass are independent of its brightness and spectrum and establish HR7672B as a rare and precious "benchmark" brown dwarf with a well-determined mass, age, and metallicity essential for testing theoretical evolutionary models and synthetic spectral models. It is presently the only directly imaged L,T,Y-dwarf known to produce an RV trend around a solar-type star.

preprint2011arXiv

Mapping the Shores of the Brown Dwarf Desert III: Young Moving Groups

We present the results of an aperture masking interferometry survey for substellar companions around 67 members of the young (~8-200Myr) nearby (~5-86pc) AB Doradus, Beta Pictoris, Hercules-Lyra, TW Hya, and Tucana-Horologium stellar associations. Observations were made at near infrared wavelengths between 1.2-3.8 microns using the adaptive optics facilities of the Keck II, VLT UT4, and Palomar Hale Telescopes. Typical contrast ratios of ~100-200 were achieved at angular separations between ~40-320mas, with our survey being 100% complete for companions with masses below 0.25\msolar across this range. We report the discovery of a $0.52 \pm 0.09$\msolar companion to HIP14807, as well as the detections and orbits of previously known stellar companions to HD16760, HD113449, and HD160934. We show that the companion to HD16760 is in a face-on orbit, resulting in an upward revision of its mass from $M_2 \sin i \sim 14$\mjupiter to $M_2 = 0.28 \pm 0.04$\msolar. No substellar companions were detected around any of our sample members, despite our ability to detect companions with masses below 80\mjupiter for 50 of our targets: of these, our sensitivity extended down to 40\mjupiter around 30 targets, with a subset of 22 subject to the still more stringent limit of 20\mjupiter. A statistical analysis of our non-detection of substellar companions allows us to place constraints on their frequency around ~0.2-1.5\msolar stars. In particular, considering companion mass distributions that have been proposed in the literature, we obtain an upper limit estimate of ~9-11% for the frequency of 20-80\mjupiter companions between 3-30AU at 95% confidence, assuming that their semimajor axes are distributed according to $d\mathcal{N}/da \propto a^{-1}$ in this range.

preprint2011arXiv

Multi-wavelength modeling of the spatially resolved debris disk of HD 107146

(abridged) We aim to constrain the location, composition, and dynamical state of planetesimal populations and dust around the young, sun-like (G2V) star HD 107146}. We consider coronagraphic observations obtained with the Advanced Camera for Surveys (HST/ACS) onboard the HST in broad V and broad I filters, a resolved 1.3mm map obtained with the Combined Array for Research in Millimeter-Wave Astronomy (CARMA), Spitzer/IRS low resolution spectra, and the spectral energy distribution (SED) of the object at wavelengths ranging from 3.5micron to 3.1mm. We complement these data with new coronagraphic high resolution observations of the debris disk using the Near Infrared Camera and Multi-Object Spectrometer (HST/NICMOS) aboard the HST in the F110W filter. The SED and images of the disk in scattered light as well as in thermal reemission are combined in our modeling using a parameterized model for the disk density distribution and optical properties of the dust. A detailed analytical model of the debris disk around HD 107146 is presented that allows us to reproduce the almost entire set of spatially resolved and unresolved multi-wavelength observations. Considering the variety of complementary observational data, we are able to break the degeneracies produced by modeling SED data alone. We find the disk to be an extended ring with a peak surface density at 131AU. Furthermore, we find evidence for an additional, inner disk probably composed of small grains released at the inner edge of the outer disk and moving inwards due to Poynting-Robertson drag. A birth ring scenario (i.e., a more or less broad ring of planetesimals creating the dust disk trough collisions) is found to be the most likely explanation of the ringlike shape of the disk.

preprint2011arXiv

Observational Constraints on Companions inside of 10 AU in the HR 8799 Planetary System

We report the results of Keck L'-band non-redundant aperture masking of HR 8799, a system with four confirmed planetary mass companions at projected orbital separations of 14 to 68 AU. We use these observations to place constraints on the presence of planets and brown dwarfs at projected orbital separations inside of 10 AU---separations out of reach to more conventional direct imaging methods. No companions were detected at better than 99% confidence between orbital separations of 0.8 to 10 AU. Assuming an age of 30 Myr and adopting the Baraffe models, we place upper limits to planetary mass companions of 80, 60, and 11 Jupiter Masses at projected orbital separations of 0.8, 1, and 3-10 AU respectively. Our constraints on massive companions to HR 8799 will help clarify ongoing studies of the orbital stability of this multi-planet system, and may illuminate future work dedicated to understanding the dust-free hole interior to ~6 AU.

preprint2010arXiv

A Spitzer IRS Study of Debris Disks Around Planet-Host Stars

Since giant planets scatter planetesimals within a few tidal radii of their orbits, the locations of existing planetesimal belts indicate regions where giant planet formation failed in bygone protostellar disks. Infrared observations of circumstellar dust produced by colliding planetesimals are therefore powerful probes of the formation histories of known planets. Here we present new Spitzer IRS spectrophotometry of 111 Solar-type stars, including 105 planet hosts. Our observations reveal 11 debris disks, including two previously undetected debris disks orbiting HD 108874 and HD 130322. Combining our 32 micron spectrophotometry with previously published MIPS photometry, we find that the majority of debris disks around planet hosts have temperatures in the range 60 < T < 100 K. Assuming a dust temperature T = 70 K, which is representative of the nine debris disks detected by both IRS and MIPS, we find that debris rings surrounding Sunlike stars orbit between 15 and 240 AU, depending on the mean particle size. Our observations imply that the planets detected by radial-velocity searches formed within 240 AU of their parent stars. If any of the debris disks studied here have mostly large, blackbody emitting grains, their companion giant planets must have formed in a narrow region between the ice line and 15 AU.

preprint2010arXiv

Atmospheric Phase Correction using CARMA-PACS: High Angular Resolution Observations of the FU-Orionis star PP 13S*

We present 0.15" resolution observations of the 227 GHz continuum emission from the circumstellar disk around the FU-Orionis star PP 13S*. The data were obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) Paired Antenna Calibration System (C-PACS), which measures and corrects the atmospheric delay fluctuations on the longest baselines of the array in order to improve the sensitivity and angular resolution of the observations. A description of the C-PACS technique and the data reduction procedures are presented. C-PACS was applied to CARMA observations of PP 13S*, which led to a factor of 1.6 increase in the observed peak flux of the source, a 36% reduction in the noise of the image, and a 52% decrease in the measured size of the source major axis. The calibrated complex visibilities were fitted with a theoretical disk model to constrain the disk surface density. The total disk mass from the best fit model corresponds to 0.06 \msun, which is larger than the median mass of a disk around a classical T Tauri star. The disk is optically thick at a wavelength of 1.3 mm for orbital radii less than 48 AU. At larger radii, the inferred surface density of the PP 13S* disk is an order of magnitude lower than that needed to develop a gravitational instability.

preprint2010arXiv

Investigating planet formation in circumstellar disks: CARMA observations of RY Tau and DG Tau

(Abridged) We present CARMA observations of the thermal dust emission from the circumstellar disks around the young stars RYTau and DGTau at wavelengths of 1.3mm and 2.8mm. The angular resolution of the maps is as high as 0.15arcsec, or 20AU at the distance of the Taurus cloud, which is a factor of 2 higher than has been achieved to date at these wavelengths. The unprecedented detail of the resulting disk images enables us to address three important questions related to the formation of planets. (1) What is the radial distribution of the circumstellar dust? (2) Does the dust emission show any indication of gaps that might signify the presence of (proto-)planets? (3) Do the dust properties depend on the orbital radius? We find that modeling the disk surface density in terms of either a classical power law or the similarity solution for viscous disk evolution, reproduces the observations well. The 1.3mm image from RYTau shows two peaks separated by 0.2arcsec with a decline in the dust emission toward the stellar position, which is significant at about 2-4sigma. For both RYTau and DGTau, the dust emission at radii larger than 15 AU displays no significant deviation from an unperturbed viscous disk model. In particular, no radial gaps in the dust distribution are detected. Under reasonable assumptions, we exclude the presence of planets more massive than 5 Jupiter masses orbiting either star at distances between about 10 and 60 AU. The radial variation of the dust opacity slope, beta, was investigated by comparing the 1.3mm and 2.8mm observations. We find mean values of beta of 0.5 and 0.7 for DGTau and RYTau respectively. Variations in beta are smaller than 0.7 between 20 and 70 AU. These results confirm that the circumstellar dust throughout these disks differs significantly from dust in the interstellar medium.

preprint2010arXiv

Millimeter imaging of MWC 758: probing the disk structure and kinematics

We investigate the structure and kinematics of the circumstellar disk around the Herbig Ae star MWC 758 using high resolution observations of the 12CO (3-2) and dust continuum emission at the wavelengths of 0.87 and 3.3 mm. We find that the dust emission peaks at an orbital radius of about 100 AU, while the CO intensity has a central peak coincident with the position of the star. The CO emission is in agreement with a disk in keplerian rotation around a 2.0 Msun star, confirming that MWC758 is indeed an intermediate mass star. By comparing the observation with theoretical disk models, we derive that the disk surface density Sigma(r) steeply increases from 40 to 100 AU, and decreases exponentially outward. Within 40 AU, the disk has to be optically thin in the continuum emission at millimeter wavelengths to explain the observed dust morphology, though our observations lack the angular resolution and sensitivity required to constrain the surface density on these spatial scales. The surface density distribution in MWC 758 disk is similar to that of ``transition'' disks, though no disk clearing has been previously inferred from the analysis of the spectral energy distribution (SED). Moreover, the asymmetries observed in the dust and CO emission suggest that the disk may be gravitationally perturbed by a low mass companion orbiting within a radius of 30 AU. Our results emphasize that SEDs alone do not provide a complete picture of disk structure and that high resolution millimeter-wave images are essential to reveal the structure of the cool disk mid plane.

preprint2009arXiv

Anachronistic Grain Growth and Global Structure of the Protoplanetary Disk Associated with the Mature Classical T Tauri Star, PDS 66

We present ATCA interferometric observations of the old (13 Myr), nearby (86pc) classical T Tauri star, PDS 66. Unresolved 3 and 12 mm continuum emission is detected towards PDS 66, and upper limits are derived for the 3 and 6 cm flux densities. The mm-wave data show a spectral slope flatter than that expected for ISM-sized dust particles, which is evidence of grain growth. We also present HST/NICMOS 1.1 micron PSF-subtracted coronagraphic imaging of PDS 66. The HST observations reveal a bilaterally symmetric circumstellar region of dust scattering about 0.32% of the central starlight, declining radially in surface brightness. The light-scattering disk of material is inclined 32 degrees from face-on, and extends to a radius of 170 AU. These data are combined with published optical and longer wavelength observations to make qualitative comparisons between the median Taurus and PDS 66 spectral energy distributions (SEDs). By comparing the near-infrared emission to a simple model, we determine that the location of the inner disk radius is consistent with the dust sublimation radius (1400 K at 0.1 AU). We place constraints on the total disk mass using a flat-disk model and find that it is probably too low to form gas giant planets according to current models. Despite the fact that PDS 66 is much older than a typical classical T Tauri star (< 5 Myr), its physical properties are not much different.

preprint2009arXiv

Debris Disks in the Upper Scorpius OB Association

We present MIPS 24um and 70um photometry for 205 members of the Upper Scorpius OB Association. These data are combined with published MIPS photometry for 15 additional association members to assess the frequency of circumstellar disks around 5 Myr old stars with spectral types between B0 and M5. Twelve stars have a detectable 70um excess, each of which also has a detectable 24um excess. A total of 54 stars are identified with a 24um excess more than 32% above the stellar photosphere. The MIPS observations reveal 19 excess sources -- 8 A/F/G stars and 11 K/M stars -- that were not previously identified with an 8um or 16um excess. The lack of short-wavelength emission and the weak 24um excess suggests that these sources are debris systems or the remnants of optically thick primordial disks with inner holes. Despite the wide range of luminosities of the stars hosting apparent debris systems, the excess characteristics are consistent with all stars having dust at similar orbital radii after factoring in variations in the radiation blowout particle size with spectral type. The results for Upper Sco are compared to similar photometric surveys from the literature to re-evaluate the evolution of debris emission. After considering the completeness limits of published surveys and the effects of stellar evolution on the debris luminosity, we find that the magnitude of the 24um excess around F-type stars increases between ages of 5 and 17 Myr as found by previous studies, but at < 2.6 sigma confidence. For B7-A9 and G0-K5 stars, any variations in the observed 24um excess emission over this age range are significant at less than 2 sigma confidence.

preprint2006arXiv

Are Debris Disks and Massive Planets Correlated?

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

preprint2006arXiv

Evidence for Mass-dependent Circumstellar Disk Evolution in the 5 Myr-old Upper Scorpius OB Association

We present 4.5, 8, and 16um photometry from the Spitzer Space Telescope for 204 stars in the Upper Scorpius OB Association. The data are used to investigate the frequency and properties of circumstellar disks around stars with masses between ~ 0.1 and 20 Msun at an age of ~ 5 Myr. We identify 35 stars that have emission at 8um or 16um in excess of the stellar photosphere. The lower mass stars (~ 0.1-1.2 Msun) appear surrounded by primordial optically thick disks based on the excess emission characteristics. Stars more massive than ~ 1.8 Msun have lower fractional excess luminosities suggesting that the inner ~ 10 AU of the disk has been largely depleted of primordial material. None of the G and F stars (~ 1.2-1.8 Msun) in our sample have an infrared excess at wavelengths <= 16um. These results indicate that the mechanisms for dispersing primordial optically thick disks operate less efficiently on average for low mass stars, and that longer time scales are available for the buildup of planetary systems in the terrestrial zone for stars with masses < 1 Msun.

preprint2000arXiv

2MASS Observations of the Perseus, Orion A, Orion B, and Monoceros R2 Molecular Clouds

We use the 2MASS Second Incremental Release Point Source Catalog to investigate the spatial distribution of young stars in the Perseus, Orion A, Orion B, and MonR2 molecular clouds. After subtracting a semi-empirical model of the field star contamination from the observed star counts, stellar surface density maps are used to identify compact clusters and any stellar population found more uniformly distributed over the molecular cloud. Each cloud contains between 2 to 7 clusters, with at least half of the cluster population found in a single, rich cluster. In addition, a distributed stellar population is inferred in the Orion A and MonR2 molecular clouds within the uncertainties of the field star subtraction with a surface density between 0.013 - 0.083 arcmin**-2. The fraction of the total stellar population contained in clusters for the nominal extinction model ranges from ~50-100% if the distributed population is relatively young (< 10 Myr), to ~25%-70% if it is relatively old (~100 Myr). The relatively high fraction of stars contained in clusters regardless of the age of the distributed population, in conjunction with the young ages generally inferred for embedded clusters in nearby molecular clouds, indicates that a substantial fraction of the total stellar population in these regions has formed within the past few million years in dense clusters. This suggests that either the star formation rate in each these clouds has recently peaked if one assumes clouds have ages > 10 Myr, or molecular clouds are younger than typically thought if one assumes that the star formation rate has been approximately constant in time.