Search

New observations with powerful infrared facilities, most notably JWST, will reveal the spectral signatures of ices in protoplanetary disks in unprecedented detail. The diversity, abundance, and distribution of disk ices play vital roles in planet formation and the delivery of biocritical volatiles to planets. Disks are physically and chemically heterogeneous, so inferring the ice properties from observations is non-trivial. I will present results from an ongoing program to forward-model disk ice spectral features, providing a powerful tool to fully understand the ice reservoirs probed by observations. The models use a kinetics-based gas-grain chemical evolution code to simulate the distribution of ices in a disk, followed by radiative transfer using a subset of key ice species to simulate the observations. I will discuss which ice species should be detectable and how the observable features vary with system properties, including stellar type, disk inclination, inheritance vs. reset initial chemical conditions, and subsequent chemical evolution. I will also highlight the value of obtaining spatially resolved spectra of edge-on disks (possible with, e.g., JWST's integral field units) to constrain the vertical distribution of ices.

We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets., Comment: This paper has been accepted for publication in the Astrophysical Journal (ApJ)

The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon-and nitrogen-rich chemistry at cold temperatures ( 1) to form carbon-bearing radicals and ions. This further results in gas phase formation of organic molecules, which then would be accreted by any actively forming planets present in the evolved disk., Accepted to Nature Astronomy, Published Dec 8th 2022

[Abridged] Planets and their atmospheres are built from gas and solid material in protoplanetary disks. This solid material grows from smaller, micron-sized grains to larger sizes in the disks, during the process of planet formation. Our goal is to model the compositional evolution of volatile ices on grains of different sizes, assuming both time-dependent grain growth and several constant grain sizes. The state-of-the-art Walsh chemical kinetics code is utilised for modeling chemical evolution. This code has been upgraded to account for the time-evolving sizes of solids. Chemical evolution is modelled locally at four different radii in a protoplanetary disk midplane for up to 10Myr. The evolution is modelled for five different constant grain sizes, and one model where the grain size changes with time according to a grain growth model appropriate for the disk midplane. Local grain growth, with conservation of total grain mass and the assumption of spherical grains, acts to reduced the total grain-surface area that is available for ice-phase reactions. This reduces these reactions efficiency compared to a chemical scenario with a conventional grain-size choice of 0.1$\mu$m. The modelled chemical evolution with grain growth leads to increased abundances of H$_{2}$O ice. For carbon in the inner disk, grain growth leads CO gas to overtake CO$_{2}$ ice as dominant carrier, and in the outer disk, CH$_{4}$ ice to become the dominant carrier. Overall, a constant grain size adopted from a grain evolution model leads to almost identical chemical evolution, when compared with chemical evolution with evolving grain sizes. A constant grain size choice, albeit larger than 0.1$\mu$m, may therefore be an appropriate simplification when approximating the impact of grain growth on chemical evolution., Comment: Accepted by Astronomy & Astrophysics. 21 pages

M-stars are the most common hosts of planetary systems in the Galaxy. Protoplanetary disks around M-stars thus offer a prime opportunity to study the chemistry of planet-forming environments. We present an ALMA survey of molecular line emission toward a sample of five protoplanetary disks around M4-M5 stars (FP Tau, J0432+1827, J1100-7619, J1545-3417, and Sz 69). These observations can resolve chemical structures down to tens of AU. Molecular lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, C$_2$H, and HCN are detected toward all five disks. Lines of H$_2$CO and DCN are detected toward 2/5 and 1/5 disks, respectively. For disks with resolved C$^{18}$O, C$_2$H, HCN, and H$_2$CO emission, we observe substructures similar to those previously found in disks around solar-type stars (e.g., rings, holes, and plateaus). C$_2$H and HCN excitation conditions estimated interior to the pebble disk edge for the bright disk J1100-7619 are consistent with previous measurements around solar-type stars. The correlation previously found between C$_2$H and HCN fluxes for solar-type disks extends to our M4-M5 disk sample, but the typical C$_2$H/HCN ratio is higher for the M4-M5 disk sample. This latter finding is reminiscent of the hydrocarbon enhancements found by previous observational infrared surveys in the innermost ($, Comment: 46 pages (20 pages in the main document, 26 pages in the appendix), 48 figures, 7 tables. Accepted in ApJ (February 2021)

H₂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 Atacama Large Millimeter/submillimeter Array survey of H₂CO toward 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H₂CO is detected toward 13 disks and tentatively detected toward a fourteenth. 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₂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₂CO excitation temperatures and column densities for four disks with multiple H₂CO line detections. The temperatures are between 20–50 K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H₂CO emission in disks generally emerges from the warm molecular layer, with some contributions from the colder midplane. Applying the same H₂CO excitation temperatures to all disks in the survey, we find that H₂CO column densities span almost three orders of magnitude (~5 × 10¹¹–5 × 10¹⁴ cm⁻²). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H₂CO compared to T Tauri disks, possibly because of less CO freeze-out. More H₂CO observations toward Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H₂CO and other oxygen-bearing organics efficiently form during planet formation.

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era., Comment: Astro2020 White Paper

One of the primary mechanisms for inferring the dynamical history of planets in our Solar System and in exoplanetary systems is through observation of elemental ratios (i.e. C/O). The ability to effectively use these observations relies critically on a robust understanding of the chemistry and evolutionary history of the observed abundances. Significant efforts have been devoted to this area from within astrochemistry circles, and these efforts should be supported going forward by the larger exoplanetary science community. In addition, the construction of a next-generation radio interferometer will be required to test many of these predictive models in situ, while simultaneously providing the resolution necessary to pinpoint the location of planets in formation., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

In this chapter we outline some of the basic understanding of the chemistry that accompanies planet formation. We discuss the basic physical environment which dictates the dominant chemical kinetic pathways for molecule formation. We focus on three zones from both observational and theoretical perspectives: (1) the planet forming midplane and ice/vapor transition zones (snow-lines), (2) the warm disk surface that is shielded from radiation, which can be readily accessed by todays observational facilities, and (3) the surface photodissociation layers where stellar radiation dominates. We end with a discussion of how chemistry influences planet formation along with how to probe the link between formation and ultimate atmospheric composition for gas giants and terrestrial worlds., Comment: Invited review, accepted for publication in the 'Handbook of Exoplanets', eds. H.J. Deeg and J.A. Belmonte, Springer (2018); 29 pages, 12 figures

Planets assemble in the midplanes of protoplanetary disks. The compositions of dust and gas in the disk midplane region determine the compositions of nascent planets, including their chemical hospitality to life. In this context, the distributions of volatile organic material across the planet and comet forming zones is of special interest. These are difficult to access in the disk midplane at IR and even millimeter wavelengths due to dust opacity, which can veil the midplane, low intrinsic molecular abundances due to efficient freeze-out, and, in the case of mid-sized organics, a mismatch between expected excitation temperatures and accessible line upper energy levels. At ngVLA wavelengths, the dust is optically thin, enabling observations into the planet forming disk midplane. ngVLA also has the requisite sensitivity. Using TW Hya as a case study, we show that ngVLA will be able to map out the distributions of diagnostic organics, such as CH3CN, in nearby protoplanetary disks., Comment: To be published in the ASP Monograph Series, 'Science with a Next-Generation VLA', ed. E. J. Murphy (ASP, San Francisco, CA)

This data set consists ofimage cubes, images, radial profiles and spectra of two H2CO emission lines towards the protoplanetary disk TW Hya as presented in Öberg et al. (2017,ApJ 839:43 'H2CO Distribution and Formation in the TW HYA Disk'). The data was observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in 2014as a part of observing program 2013.1.00114.S.In particular the data set corresponds to the channel maps presented in Fig. 1, and the Moment 0 maps, radial profiles and spectra presented in Fig. 2. The image cubes from which all other data products were derived were produced using CASA following standard calibration and imaging procedures as outlined in the publication. The raw data is archived by the ALMA observatory and can be found here: http://almascience.nrao.edu/aq/ All publications making use of this data must include the following acknowledgement: This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00114.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $

From the masses of planets orbiting our Sun, and relative elemental abundances, it is estimated that at birth our Solar System required a minimum disk mass of ~0.01 solar masses within ~100 AU of the star. The main constituent, gaseous molecular hydrogen, does not emit from the disk mass reservoir, so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally probes the disk surface, while the conversion from dust emission to gas mass requires knowledge of the grain properties and gas-to-dust mass ratio, which likely differ from their interstellar values. Thus, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3--10 Myr) star TW Hya, with estimates ranging from 0.0005 to 0.06 solar masses. Here we report the detection the fundamental rotational transition of hydrogen deuteride, HD, toward TW Hya. HD is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The HD detection, combined with existing observations and detailed models, implies a disk mass >0.05 solar masses, enough to form a planetary system like our own., 16 pages, 3 figures. Pre-editorial version of paper published in Nature. Full version can be found at http://www.nature.com.proxy.lib.umich.edu/nature/journal/v493/n7434/full/nature11805.html

We model the radiatively driven flow from IM Lup - a large protoplanetary disc expected to be irradiated by only a weak external radiation field (at least 104 times lower than the ultraviolet field irradiating the Orion Nebula Cluster proplyds). We find that material at large radii (>400 au) in this disc is sufficiently weakly gravitationally bound that significant mass-loss can be induced. Given the estimated values of the disc mass and accretion rate, the viscous time-scale is long (˜10 Myr) so the main evolutionary behaviour for the first Myr of the disc's lifetime is truncation of the disc by photoevaporation, with only modest changes effected by viscosity. We also produce approximate synthetic observations of our models, finding substantial emission from the flow that can explain the CO halo observed about IM Lup out to ≥1000 au. Solutions that are consistent with the extent of the observed CO emission generally imply that IM Lup is still in the process of having its disc outer radius truncated. We conclude that IM Lup is subject to substantial external photoevaporation, which raises the more general possibility that external irradiation of the largest discs can be of significant importance even in low mass star forming regions.