Researcher profile

Richard Teague

Richard Teague contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 published item(s)

preprint2026arXiv

Tracing Pebble Drift History in Two Protoplanetary Disks with CO Enhancement

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

preprint2024arXiv

Mapping the Vertical Gas Structure of the Planet-hosting PDS 70 Disk

PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (${\approx}$0."1;10 au) observations in a set of CO isotopologue lines and HCO$^+$ J=4-3 to map the full 2D $(r,z)$ disk structure from the disk atmosphere, as traced by $^{12}$CO, to closer to the midplane, as probed by less abundant isotopologues and HCO$^+$. In the PDS 70 disk, $^{12}$CO traces a height of $z/r\approx0.3$, $^{13}$CO is found at $z/r\approx0.1$, and C$^{18}$O originates at, or near, the midplane. The HCO$^+$ surface arises from $z/r\approx0.2$ and is one of the few non-CO emission surfaces constrained with high fidelity in disks to date. In the $^{12}$CO J=3-2 line, we resolve a vertical dip and steep rise in height at the cavity wall, making PDS 70 the first transition disk where this effect is directly seen in line emitting heights. In the outer disk, the CO emission heights of PDS 70 appear typical for its stellar mass and disk size and are not substantially altered by the two inner embedded planets. By combining CO isotopologue and HCO$^+$ lines, we derive the 2D gas temperature structure and estimate a midplane CO snowline of ${\approx}$56-85 au. This implies that both PDS 70b and 70c are located interior to the CO snowline and are likely accreting gas with a high C/O ratio of ${\approx}$1.0, which provides context for future planetary atmospheric measurements from, e.g., JWST, and for properly modeling their formation histories.

preprint2022arXiv

A localized kinematic structure detected in atomic carbon emission spatially coincident with a proposed protoplanet in the HD 163296 disk

Over the last five years, studies of the kinematics in protoplanetary disks have led to the discovery of new protoplanet candidates and several structures linked to possible planet-disk interactions. We detect a localized kinematic bipolar structure in the HD 163296 disk present inside the deepest dust gap at 48 au from atomic carbon line emission. HD 163296's stellar jet and molecular winds have been described in detail in the literature; however, the kinematic anomaly in C I emission is not associated with either of them. Further, the velocity of the kinematic structure points indicates a component fast enough to differentiate it from the Keplerian profile of the disk; and its atomic nature hints at a localized UV source strong enough to dissociate CO and launch a C I outflow, or a strong polar flow from the upper layers of the disk. By discarding the stellar jet and previously observed molecular winds, we explore different sources for this kinematic feature in C I emission that could be associated with a protoplanet inflow/outflow, or disk winds.

preprint2022arXiv

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

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

preprint2022arXiv

Gas and Dust Shadows in the TW Hydrae Disk

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $χ$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $β$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies.

preprint2022arXiv

Unveiling the outer dust disc of TW Hya with deep ALMA observations

The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs.

preprint2021arXiv

Dynamical Masses and Stellar Evolutionary Model Predictions of M-Stars

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

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

A Three Dimensional View of Gomez's Hamburger

Unraveling the 3D physical structure, the temperature and density distribution, of protoplanetary discs is an essential step if we are to confront simulations of embedded planets or dynamical instabilities. In this paper we focus on Submillimeter Array observations of the edge-on source, Gomez&#39;s Hamburger, believed to host an over-density hypothesised to be a product of gravitational instability in the disc, GoHam b. We demonstrate that, by leveraging the well characterised rotation of a Keplerian disc to deproject observations of molecular lines in position-position-velocity space into disc-centric coordinates, we are able to map out the emission distribution in the (r, z) plane and (x, y, z) space. We show that 12CO traces an elevated layer of $z\,/\,r \sim 0.3$, while 13CO traces deeper in the disc at $z\,/\,r \lesssim 0.2$. We localize emission associated with GoHam b, finding it at deprojected radius of approximately 500 au and at polar angle of +\- 30 degrees from the disc major axis. At the spatial resolution of $\sim 1.5^{\prime\prime}$, GoHam b is spatially unresolved, with an upper limit to its radius of $<190$~au.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) V: Sample, overview, and demography of disk molecular emission

We present an overview of the ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT), a campaign devoted to the characterization of the molecular emission from partly embedded, young stars. The project aims at better understanding the gaseous products that are delivered to planets by means of high-resolution maps of assorted lines probing disks at the time of the planet formation (less than 1 Myr). Nine different molecules are surveyed by our observations of six Class I/flat-spectrum sources. A series of accompanying articles analyze specific targets and molecules. Here we describe the sample and provide a general overview of the results, focusing on the spatial distribution, column densities, and abundance ratios of H$_2$CO, CS, and CN. The results of this work are a first step toward the characterization of the disk chemical evolution that need to be complemented by further observations of less exceptional disks and customized thermo-chemical modeling.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT): I. CO, CS, CN, and H2CO around DG Tau B

The chemical composition of planets is inherited by the distribution of the various molecular species in the protoplanetary disk at thetime of their formation. As of today, only a handful of disks has been imaged in multiple spectral lines with high spatial resolution. As part of a small campaign devoted to the chemical characterization of disk-outflow sources in Taurus, we report on new ALMA Band 6 observations with 20 au resolution toward the embedded young star DG Tau B. Images of the continuum emission reveals a dust disk with rings and, putatively, a leading spiral arm. The disk, as well as the prominent outflow cavities, are detected in CO, H2CO, CS, and CN while they remain undetected in SO2, HDO, and CH3OH. From the absorption of the back-side outflow, we inferred that the disk emission is optically thick in the inner 50 au. This morphology explains why no line emission is detected from this inner region and poses some limitations toward the calculation of the dust mass and the characterization of the inner gaseous disk. The H2CO and CS emission from the inner 200 au is mostly from the disk and their morphology is very similar. The CN emission significantly differs from the other two molecules as it is observed only beyond 150 au. This ring-like morphology is consistent with previous observations and the predictions of thermochemical disk models. Finally, we constrained the disk-integrated column density of all molecules. In particular, we found that the CH3OH/H2CO ratio must be smaller than 2, making the methanol non-detection still consistent with the only such a ratio available from the literature (1.27 in TW Hya).

preprint2020arXiv

Hints of a Population of Solar System Analog Planets from ALMA

The recent ALMA DSHARP survey provided illuminating results on the diversity of substructures in planet forming disks. These substructures trace pebble-sized grains accumulated at local pressure maxima, possibly due to planet-disk interactions or other planet formation processes. DSHARP sources are heavily biased to large and massive disks that only represent the high (dust flux) tail end of the disk population. Thus it is unclear whether similar substructures and corresponding physical processes also occur in the majority of disks which are fainter and more compact. Here we explore the presence and characteristics of features in a compact disk around GQ Lup A, the effective radius of which is 1.5 to 10 times smaller than those of DSHARP disks. We present our analysis of ALMA 1.3mm continuum observations of the GQ Lup system. By fitting visibility profiles of the continuum emission, we find substructures including a gap at ~ 10 au. The compact disk around GQ Lup exhibits similar substructures to those in the DSHARP sample, suggesting that mechanisms of trapping pebble-sized grains are at work in small disks as well. Characteristics of the feature at ~ 10 au, if due to a hidden planet, are evidence of planet formation at Saturnian distances. Our results hint at a rich world of substructures to be identified within the common population of compact disks, and subsequently a population of solar system analogs within these disks. Such study is critical to understanding the formation mechanisms and planet populations in the majority of protoplanetary disks.

preprint2020arXiv

The efficiency of dust trapping in ringed proto-planetary discs

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $α/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($α/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $α$ turbulent parameter ($α\sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

preprint2020arXiv

The Excitation Conditions of CN in TW Hya

We report observations of the cyanide anion, CN, in the disk around TW~Hya covering the $N=1-0$, $N=2-1$ and $N=3-2$ transitions. Using line stacking techniques, 24 hyperfine transitions are detected out of the 30 within the observed frequency ranges. Exploiting the super-spectral resolution from the line stacking method reveals the splitting of hyperfine components previously unresolved by laboratory spectroscopy. All transitions display a similar emission morphology, characterized by an azimuthally symmetric ring, peaking at $\approx 45$~au (0.75&#34;), and a diffuse outer tail extending out to the disk edge at $\approx 200$~au. Excitation analyses assuming local thermodynamic equilibrium (LTE) yield excitation temperatures in excess of the derived kinetic temperatures based on the local line widths for all fine structure groups, suggesting assumptions of LTE are invalid. Using the 0D radiative transfer code RADEX, we demonstrate that such non-LTE effects may be present when the local H$_2$ density drops to $10^{7}~{\rm cm^{-3}}$ and below. Comparison with models of TW~Hya find similar densities at elevated regions in the disk, typically $z \, / \, r \gtrsim 0.2$, consistent with model predictions where CN is formed via vibrationally excited H$_2$ in the disk atmospheric layers where UV irradiation is less attenuated.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future.

preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the &#39;Visualizing the Kinematics of Planet Formation&#39; workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.