CO Snow Research Articles

Comet 2I/Borisov in Comparison with Comets of the Solar System

The brief review summarizes data on the chemical and mineral composition, as well as on the phys ical properties, of the first extrasolar comet 2I/Borisov, obtained from observations that were carried out from September 2019 to the end of March 2020. It is noted that the qualitative chemical composition of the volatile and mineral components comet 2I/Borisov is similar to the composition of comets in the Solar System, but there are differences that indicate the specific conditions for the formation of its nucleus in a circumstellar gas and dust disk. Different release rates of CO and H2O molecules in the vicinity of perihelion indicate the pos sible heterogeneity of the comet’s nucleus, which was formed from more homogeneous ice blocks, but differ ing in composition. These constituent blocks could have formed over a wide range of radial distances: from the snow line of H2O to the CO snow line. Their accumulation in the comet’s nucleus indicates large-scale mixing of protocometary bodies in the circumstellar disk. No spectra of finely crystalline magnesium silicates were found in cometary coma of 2I/Borisov, which can be interpreted as the absence of a significant amount of gas and dust transfer from the inner hot regions of the disk to the outside, into the zone of formation of protocometary bodies.

A Survey of CO, CO2, and H2O in Comets and Centaurs

CO and CO2 are the two dominant carbon-bearing molecules in comae and have major roles in driving activity. Their relative abundances also provide strong observational constraints to models of solar system formation and evolution but have never been studied together in a large sample of comets. We carefully compiled and analyzed published measurements of simultaneous CO and CO2 production rates for 25 comets. Approximately half of the comae have substantially more CO2 than CO, about one-third are CO-dominated, and about 1/10 produce a comparable amount of both. There may be a heliocentric dependence to this ratio with CO dominating comae beyond 3.5 au. Eight out of nine of the Jupiter Family Comets in our study produce more CO2 than CO. The six dynamically new comets produce more CO2 relative to CO than the eight Oort Cloud comets that have made multiple passes through the inner solar system. This may be explained by long-term cosmic-ray processing of a comet nucleus’s outer layers. We find (Q CO/Q )median = 3 ± 1% and (Q /Q )median = 12 ± 2%. The inorganic volatile carbon budget was estimated to be (Q CO+Q )/Q ∼ 18% for most comets. Between 0.7 and 4.6 au, CO2 outgassing appears to be more intimately tied to the water production in a way that the CO is not. The volatile carbon/oxygen ratio for 18 comets is C/Omedian ∼ 13%, which is consistent with a comet formation environment that is well within the CO snow line.

Broadband spectroscopy of astrophysical ice analogues

Context. Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Such data are still missing from the literature, which can be attributed to the lack of appropriate spectroscopic systems and methods for laboratory studies. Aims. In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are combined to study optical constants of CO and CO2 ices in the broad THz-IR spectral range. Methods. The measured ices were grown at cryogenic temperatures by gas deposition on a cold silicon window. We developed a method to quantify the broadband THz-IR optical constants of ices, based on the direct reconstruction of the complex refractive index of ices in the THz range from the TDS data and the use of the Kramers-Kronig relation in the IR range for the reconstruction from the FTIR data. Uncertainties introduced by the Kramers-Kronig relations were eliminated by merging the THz and IR spectra. Finally, the reconstructed THz-IR response was analyzed using classical models of complex dielectric permittivity. Results. The complex refractive index of CO and CO2 ices deposited at the temperature of 28 K was obtained in the range of 0.312.0 THz and fitted using the analytical Lorentz model. Based on the measured dielectric constants, opacities of the astrophysical dust with CO and CO2 icy mantles were computed. Conclusions. The method developed in this work can be used for a model-independent reconstructions of optical constants of various astrophysical ice analogs in a broad THz-IR range. Such data can provide important benchmarks for interpreting broadband observations from existing and future ground-based facilities and space telescopes. The reported results will be useful in modeling sources that exhibit a drastic molecular freeze-out, such as the central regions of prestellar cores and mid-planes of protoplanetary disks, as well as CO and CO2 snow lines in disks.

Discovery of Line Pressure Broadening and Direct Constraint on Gas Surface Density in a Protoplanetary Disk

The gas surface density profile of protoplanetary disks is one of the most fundamental physical properties to understanding planet formation. However, it is challenging to determine the surface density profile observationally, because the H2 emission cannot be observed in low-temperature regions. We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the 12CO J = 3 − 2 line toward the protoplanetary disk around TW Hya and discovered extremely broad line wings due to the pressure broadening. In conjunction with a previously reported optically thin CO isotopologue line, the pressure broadened line wings enabled us to directly determine the midplane gas density for the first time. The gas surface density at ∼5 au from the central star reaches ∼103 g cm−2, which suggests that the inner region of the disk has enough mass to form a Jupiter-mass planet. Additionally, the gas surface density drops at the inner cavity by ∼2 orders of magnitude compared to outside the cavity. We also found a low CO abundance of ∼10−6 with respect to H2, even inside the CO snow line, which suggests conversion of CO to less volatile species. Combining our results with previous studies, the gas surface density jumps at r ∼ 20 au, suggesting that the inner region (3 < r < 20 au) might be the magnetorotational instability dead zone. This study sheds light on the direct gas surface density constraint without assuming the CO/H2 ratio using ALMA.

The Volatile Carbon-to-oxygen Ratio as a Tracer for the Formation Locations of Interstellar Comets

Based on the occurrence rates implied by the discoveries of 1I/‘Oumuamua and 2I/Borisov, the forthcoming Rubin Observatory Legacy Survey of Space and Time (LSST) should detect ≥one interstellar object every year. We advocate for future measurements of the production rates of H2O, CO2, and CO in these objects to estimate their carbon-to-oxygen ratios, which trace formation locations within their original protoplanetary disks. We review similar measurements for solar system comets, which indicate formation interior to the CO snow line. By quantifying the relative processing in the interstellar medium and solar system, we estimate that production rates will not be representative of primordial compositions for the majority of interstellar comets. Preferential desorption of CO and CO2 relative to H2O in the interstellar medium implies that measured C/O ratios represent lower limits on the primordial ratios. Specifically, production rate ratios of Q(CO)/Q(H2O) < 0.2 and Q(CO)/Q(H2O) > 1 likely indicate formation interior and exterior to the CO snow line, respectively. The high C/O ratio of 2I/Borisov implies that it formed exterior to the CO snow line. We provide an overview of the currently operational facilities capable of obtaining these measurements that will constrain the fraction of ejected comets that formed exterior to the CO snow line. This fraction will provide key insights into the efficiency of and mechanisms for cometary ejection in exoplanetary systems.

The 12CO/13CO isotopologue ratio of a young, isolated brown dwarf

Context. Linking atmospheric characteristics of planets to their formation pathways is a central theme in the study of extrasolar planets. Although the 12C/13C isotope ratio shows little variation in the Solar System, the atmosphere of a super-Jupiter was recently shown to be rich in 13CO, possibly as a result of dominant ice accretion beyond the CO snow line during its formation. Carbon isotope ratios are therefore suggested to be a potential tracer of formation pathways of planets. Aims. In this work, we aim to measure the 12CO/13CO isotopologue ratio of a young, isolated brown dwarf. While the general atmospheric characteristics of young, low-mass brown dwarfs are expected to be very similar to those of super-Jupiters, their formation pathways may be different, leading to distinct isotopologue ratios. In addition, such objects allow high-dispersion spectroscopy at high signal-to-noise ratios. Methods. We analysed archival K-band spectra of the L dwarf 2MASS J03552337+1133437 taken with NIRSPEC at the Keck telescope. A free retrieval analysis was applied to the data using the radiative transfer code petitRADTRANS coupled with the nested sampling tool PyMultiNest to determine the isotopologue ratio 12CO/13CO in its atmosphere. Results. The isotopologue 13CO is detected in the atmosphere through the cross-correlation method at a signal-to-noise of ~8.4. The detection significance is determined to be ~9.5σ using a Bayesian model comparison between two retrieval models (including or excluding 13CO). We retrieve an isotopologue 12CO/13CO ratio of 97−18+25 (90% uncertainty), marginally higher than the local interstellar standard. Its C/O ratio of ~0.56 is consistent with the solar value. Conclusions. Although only one super-Jupiter and one brown dwarf now have a measured 12CO/13CO ratio, it is intriguing that they are different, possibly hinting to distinct formation pathways. Regardless of spectroscopic similarities, isolated brown dwarfs may experience a top-down formation via gravitational collapse, which resembles star formation, while giant exoplanets favourably form through core accretion, which potentially alters isotopologue ratios in their atmospheres depending on the material they accrete from protoplanetary disks. This further emphasises atmospheric carbon isotopologue ratio as a tracer of the formation history of exoplanets. In the future, analyses such as those presented here should be conducted on a wide range of exoplanets using medium-to-high-resolution spectroscopy to further assess planet formation processes.

Molecules with ALMA at Planet-forming Scales (MAPS). XIII. HCO+ and Disk Ionization Structure

Abstract We observed HCO+ J = 1 − 0 and H13CO+ J = 1 − 0 emission toward the five protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 as part of the MAPS project. HCO+ is detected and mapped at 0.″3 resolution in all five disks, while H13CO+ is detected (S/N &gt; 6σ) toward GM Aur and HD 163296 and tentatively detected (S/N &gt; 3σ) toward the other disks by a matched filter analysis. Inside a radius of R ∼ 100 au, the HCO+ column density is flat or shows a central dip. At outer radii (≳100 au), the HCO+ column density decreases outward, while the column density ratio of HCO+/CO is mostly in the range of ∼10−5–10−4. We derived the HCO+ abundance in the warm CO-rich layer, where HCO+ is expected to be the dominant molecular ion. At R ≳ 100 au, the HCO+ abundance is ∼3 × 10−11 − 3 × 10−10, which is consistent with a template disk model with X-ray ionization. At the smaller radii, the abundance decreases inward, which indicates that the ionization degree is lower in denser gas, especially inside the CO snow line, where the CO-rich layer is in the midplane. Comparison of template disk models with the column densities of HCO+, N2H+, and N2D+ indicates that the midplane ionization rate is ≳10−18 s−1 for the disks around IM Lup, AS 209, and HD 163296. We also find hints of an increased HCO+ abundance around the location of dust continuum gaps in AS 209, HD 163296, and MWC 480. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

If you like C/O variations, you should have put a ring on it

Context. The C/O ratio as traced with C2H emission in protoplanetary disks is fundamental for constraining the formation mechanisms of exoplanets and for our understanding of volatile depletion in disks, but current C2H observations show an apparent bimodal distribution that is not well understood, indicating that the C/O distribution is not described by a simple radial dependence. Aims. The transport of icy pebbles has been suggested to alter the local elemental abundances in protoplanetary disks through settling, drift, and trapping in pressure bumps, resulting in a depletion of volatiles in the surface layer and an increase in the elemental C/O. Methods. We combine all disks with spatially resolved ALMA C2H observations with high-resolution continuum images and constraints on the CO snow line to determine if the C2H emission is indeed related to the location of the icy pebbles. Results. We report a possible correlation between the presence of a significant CO-ice dust reservoir and high C2H emission, which is only found in disks with dust rings outside the CO snow line. In contrast, compact dust disks (without pressure bumps) and warm transition disks (with their dust ring inside the CO snow line) are not detected in C2H, suggesting that such disks may have never contained a significant CO ice reservoir. Conclusions. This correlation provides evidence for the regulation of the C/O profile by the complex interplay of CO snow line and pressure bump locations in the disk. These results demonstrate the importance of including dust transport in chemical disk models for a proper interpretation of exoplanet atmospheric compositions and a better understanding of volatile depletion in disks, in particular the use of CO isotopologs to determine gas surface densities.

Warm millimetre dust in protoplanetary discs near massive stars

ABSTRACT Dust plays a key role in the formation of planets and its emission also provides one of our most accessible views of protoplanetary discs. If set by radiative equilibrium with the central star, the temperature of dust in the disc plateaus at around 10–20 K in the outer regions. However, sufficiently nearby massive stars can heat the outer disc to substantially higher temperatures. In this paper, we study the radiative equilibrium temperature of discs in the presence of massive external sources and gauge the effect that it has on millimetre dust mass estimates. Since millimetre grains are not entrained in any wind, we focus on geometrically simple 2D axisymmetric disc models using radiative transfer calculations with both the host star and an external source. Recent surveys have searched for evidence of massive stars influencing disc evolution using disc properties as a function of projected separation. In assuming a disc temperature of 20 K for a disc a distance D from a strong radiation source, disc masses are overestimated by a factor that scales with D−1/2 interior to the separation that external heating becomes important. This could significantly alter dust mass estimates of discs in close proximity to θ1C in the Orion Nebular Cluster (ONC). We also make an initial assessment of the effect upon snow lines. Within a parsec of an O star like θ1C a CO snow line no longer exists, though the water snow line is virtually unaffected except for very close separations of $\le 0.01\,$pc.

Using HCO+ isotopologues as tracers of gas depletion in protoplanetary disk gaps

Context. The widespread rings and gaps seen in the dust continuum in protoplanetary disks are sometimes accompanied by similar substructures seen in molecular line emission. One example is the outer gap at ~100 au in AS 209, which shows that the H13CO+ and C18O emission intensities decrease along with the continuum in the gap, while the DCO+ emission increases inside the gap. Aims. We aim to study the behavior of DCO+/H13CO+ and DCO+/HCO+ ratios in protoplanetary disk gaps assuming the two scenarios: (A) the gas depletion follows the dust depletion and (B) only the dust is depleted. Methods. We first modeled the physical disk structure using the thermo-chemical model ANDES. This 1+1D steady-state disk model calculates the thermal balance of gas and dust and includes the far ultraviolet, X-rays, cosmic rays, and other ionization sources together with the reduced chemical network for molecular coolants. Afterward, this physical structure was adopted for calculations of molecular abundances with the extended gas-grain chemical network with deuterium fractionation. Ideal synthetic spectra and 0th-moment maps were produced with the LIne Modeling Engine. Results. We are able to qualitatively reproduce the increase in the DCO+ intensity and the decrease in the H13CO+ and C18O intensities inside the disk gap, which is qualitatively similar to what is observed in the outer AS 209 gap. The corresponding disk model (A) assumes that both the gas and dust are depleted in the gap. The model (B) with the gas-rich gap, where only the dust is depleted, produces emission that is too bright in all HCO+ isotopologues and C18O. Conclusions. The DCO+/H13CO+ line ratio can be used to probe gas depletion in dust continuum gaps outside of the CO snow line. The DCO+/C18O line ratio shows a similar, albeit weaker, effect; however, these species can be observed simultaneously with a single (sub)mm interferometer setup.

CO destruction in protoplanetary disk midplanes: Inside versus outside the CO snow surface

Context. The total gas mass is one of the most fundamental properties of disks around young stars, because it controls their evolution and their potential to form planets. To measure disk gas masses, CO has long been thought to be the best tracer as it is readily detected at (sub)mm wavelengths in many disks. However, inferred gas masses from CO in recent ALMA observations of large samples of disks in the 1–5 Myr age range seem inconsistent with their inferred dust masses. The derived gas-to-dust mass ratios from CO are between one and two orders of magnitude lower than the ISM value of ~100 even if photodissociation and freeze-out are included. In contrast, Herschel measurements of hydrogen deuteride line emission of a few disks imply gas masses in line with gas-to-dust mass ratios of 100. This suggests that at least one additional mechanism is removing CO from the gas phase. Aims. Here we test the suggestion that the bulk of the CO is chemically processed and that the carbon is sequestered into less volatile species such as CO2, CH3OH, and CH4 in the dense, shielded midplane regions of the disk. This study therefore also addresses the carbon reservoir of the material which ultimately becomes incorporated into planetesimals. Methods. Using our gas-grain chemical code, we performed a parameter exploration and follow the CO abundance evolution over a range of conditions representative of shielded disk midplanes. Results. Consistent with previous studies, we find that no chemical processing of CO takes place on 1–3 Myr timescales for low cosmic-ray ionisation rates, &lt;5 × 10−18 s−1. Assuming an ionisation rate of 10−17 s−1, more than 90% of the CO is converted into other species, but only in the cold parts of the disk below 30 K. This order of magnitude destruction of CO is robust against the choice of grain-surface reaction rate parameters, such as the tunnelling efficiency and diffusion barrier height, for temperatures between 20 and 30 K. Below 20 K there is a strong dependence on the assumed efficiency of H tunnelling. Conclusions. The low temperatures needed for CO chemical processing indicate that the exact disk temperature structure is important, with warm disks around luminous Herbig stars expected to have little to no CO conversion. In contrast, for cold disks around sun-like T Tauri stars, a large fraction of the emitting CO layer is affected unless the disks are young (&lt;1 Myr). This can lead to inferred gas masses that are up to two orders of magnitude lower. Moreover, unless CO is locked up early in large grains, the volatile carbon composition of the icy pebbles and planetesimals forming in the midplane and drifting to the inner disk will be dominated by CH3OH, CO2 and/or hydrocarbons.

An Empirical Planetesimal Belt Radius–Stellar Luminosity Relation

Abstract Resolved observations of millimeter-sized dust, tracing larger planetesimals, have pinpointed the location of 26 Edgeworth–Kuiper Belt analogs. We report that a belt’s distance R to its host star correlates with the star’s luminosity L ⋆, following with a low intrinsic scatter of ∼17%. Remarkably, our Edgeworth–Kuiper Belt in the solar system and the two CO snow lines imaged in protoplanetary disks lie close to this R–L ⋆ relation, suggestive of an intrinsic relationship between protoplanetary disk structures and belt locations. To test the effect of bias on the relation, we use a Monte Carlo approach and simulate uncorrelated model populations of belts. We find that observational bias could produce the slope and intercept of the R–L ⋆ relation but is unable to reproduce its low scatter. We then repeat the simulation taking into account the collisional evolution of belts, following the steady-state model that fits the belt population as observed through infrared excesses. This significantly improves the fit by lowering the scatter of the simulated R–L ⋆ relation; however, this scatter remains only marginally consistent with the one observed. The inability of observational bias and collisional evolution alone to reproduce the tight relationship between belt radius and stellar luminosity could indicate that planetesimal belts form at preferential locations within protoplanetary disks. The similar trend for CO snow line locations would then indicate that the formation of planetesimals or planets in the outer regions of planetary systems is linked to the volatility of their building blocks, as postulated by planet formation models.

Probing Episodic Accretion in Very Low Luminosity Objects

Abstract Episodic accretion has been proposed as a solution to the long-standing luminosity problem in star formation; however, the process remains poorly understood. We present observations of line emission from N2H+ and CO isotopologues using the Atacama Large Millimeter/submillimeter Array (ALMA) in the envelopes of eight very low luminosity objects (VeLLOs). In five of the sources the spatial distribution of emission from N2H+ and CO isotopologues shows a clear anticorrelation. It is proposed that this is tracing the CO snow line in the envelopes: N2H+ emission is depleted toward the center of these sources, in contrast to the CO isotopologue emission, which exhibits a peak. The positions of the CO snow lines traced by the N2H+ emission are located at much larger radii than those calculated using the current luminosities of the central sources. This implies that these five sources have experienced a recent accretion burst because the CO snow line would have been pushed outward during the burst because of the increased luminosity of the central star. The N2H+ and CO isotopologue emission from DCE161, one of the other three sources, is most likely tracing a transition disk at a later evolutionary stage. Excluding DCE161, five out of seven sources (i.e., ∼70%) show signatures of a recent accretion burst. This fraction is larger than that of the Class 0/I sources studied by Jørgensen et al. and Frimann et al., suggesting that the interval between accretion episodes in VeLLOs is shorter than that in Class 0/I sources.

Direct mapping of the temperature and velocity gradients in discs

Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of 21 ± 2 K. In the outer disc (&gt;300 au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to ≈15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.

Disk Evolution and the Fate of Water

We review the general theoretical concepts and observational constraints on the distribution and evolution of water vapor and ice in protoplanetary disks, with a focus on the Solar System. Water is expected to freeze out at distances greater than 1-3 AU from solar-type central stars; more precise estimates are difficult to obtain due to uncertainties in the complex processes involved in disk evolution, including dust growth, settling, and radial drift, and the level of turbulence and viscous dissipation within disks. Interferometric observations are now providing constraints on the positions of CO snow lines, but extrapolation to the unresolved regions where water ice sublimates will require much better theoretical understanding of mass and angular momentum transport in disks as well as more refined comparison of observations with sophisticated disk models.

Polarized Disk Emission from Herbig Ae/Be Stars Observed Using Gemini Planet Imager: HD 144432, HD 150193, HD 163296, and HD 169142

Abstract In order to look for signs of ongoing planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager, along with new H band observations of HD 144432. We confirm the complex “double ring” structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 au above the midplane at a radial distance of 77 au, co-spatial with a ring seen at 1.3 mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

Probing the CO and methanol snow lines in young protostars

Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)

ALMA OBSERVATIONS OF A GAP AND A RING IN THE PROTOPLANETARY DISK AROUND TW HYA

ABSTRACT We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ∼30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ∼15 au. In addition, the 13CO and C18O J = 3 &minus; 2 ?> lines show a decrement in CO line emission throughout the disk, down to ∼10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2 M Neptune ?> ), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

ON THE COMMONALITY OF 10–30 AU SIZED AXISYMMETRIC DUST STRUCTURES IN PROTOPLANETARY DISKS

ABSTRACT An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive millimeter-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small-scale structures. Here, we investigate dust structures in the (sub)millimeter continuum emission of four disks (TW Hya, HL Tau, HD 163296, and DM Tau), a sample of disks with the highest spatial resolution Atacama Large Millimeter/submillimeter Array observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ∼axisymmetric but rich in 10–30 AU-sized radial structures, possibly due to physical gaps, surface density enhancements, or localized dust opacity variations within the disks. These results suggest that small-scale axisymmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snow lines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snow lines.

Steepening of the 820μm continuum surface brightness profile signals dust evolution in TW Hydrae’s disk

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