Search

The composition of a planet's atmosphere is intricately linked to the chemical makeup of the protoplanetary disk in which it formed. Determining the elemental abundances from key volatiles within disks is therefore essential for establishing connections between the composition of disks and planets. The disk around the Herbig Ae star HD 169142 is a compelling target for such a study due to its molecule-rich nature and the presence of a newly-forming planet between two prominent dust rings. In this work, we probe the chemistry of the HD 169142 disk at small spatial scales, drawing links between the composition of the disk and the planet-accreted gas. Using thermochemical models and archival data, we constrain the elemental abundances of volatile carbon, oxygen, and sulfur. Carbon and oxygen are only moderately depleted from the gas phase relative to their interstellar abundances, with the inner 60 au appearing enriched in volatile oxygen. The C/O ratio is approximately solar within the inner disk and rises above this in the outer disk, as expected across the H2O snowline. The gas-phase sulfur abundance is depleted by a factor of 1000, consistent with a number of other protoplanetary disks. Interestingly, the observed SiS emission near the HD 169142 b protoplanet vastly exceeds chemical model predictions, supporting previous hypotheses suggesting its origin in shocked gas or a localised outflow. We contextualise our findings in terms of the potential atmospheric composition of the embedded planet, and highlight the utility of sulfur-bearing molecules as probes of protoplanetary disk chemistry., Comment: Accepted for publication in MNRAS

Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three transition disks around Herbig Ae stars all show oxygen-rich gas traced via the unique detections of the molecules SO and CH3OH. We present the first results of an ALMA line survey at 337 to 357 GHz of such disks and focus this paper on the first Herbig Ae disk to exhibit this chemical signature - HD 100546. In these data, we detect 19 different molecules including NO, SO and CH3OCHO (methyl formate). We also make the first tentative detections of H213CO and 34SO in protoplanetary disks. Multiple molecular species are detected in rings, which are, surprisingly, all peaking just beyond the underlying millimeter continuum ring at 200 au. This result demonstrates a clear connection between the large dust distribution and the chemistry in this flat outer disk. We discuss the physical and/or chemical origin of these sub-structures in relation to ongoing planet formation in the HD 100546 disk. We also investigate how similar and/or different the molecular make up of this disk is to other chemically well-characterised Herbig Ae disks. The line-rich data we present motivates the need for more ALMA line surveys to probe the observable chemistry in Herbig Ae systems which offer unique insight into the composition of disk ices, including complex organic molecules., Comment: Accepted to AJ, 25 pages, 11 figures

The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in these disks. The IRS 48 disk shows a particularly rich chemistry that can be directly linked to its asymmetric dust trap. Here, we present ALMA observations of the IRS 48 disk where we detect 16 different molecules and make the first robust detections of H213CO, 34SO, 33SO and c-H2COCH2 (ethylene oxide) in a protoplanetary disk. All of the molecular emissions, aside from CO, are colocated with the dust trap and this includes newly detected simple molecules such as HCO+, HCN and CS. Interestingly, there are spatial offsets between different molecular families, including between the COMs and sulphur-bearing species, with the latter being more azimuthally extended and located radially further from the star. The abundances of the newly detected COMs relative to CH3OH are higher than the expected protostellar ratios, which implies some degree of chemical processing of the inherited ices during the disk lifetime. These data highlight IRS 48 as a unique astrochemical laboratory to unravel the full volatile reservoir at the epoch of planet and comet formation and the role of the disk in (re)setting chemical complexity., Comment: Accepted to AJ, 21 pages, 7 figures

Volatile elements play a crucial role in the formation of planetary systems. Their abundance and distribution in protoplanetary disks provide vital insights into the connection between formation processes and the atmospheric composition of individual planets. Sulfur, being one of the most abundant elements in planet-forming environments, is of great significance, and now observable in exoplanets with JWST. However, planetary formation models currently lack vital knowledge regarding sulfur chemistry in protoplanetary disks. Developing a deeper understanding of the major volatile sulfur carriers in disks is essential to building models that can meaningfully predict planetary atmospheric composition, and reconstruct planetary formation pathways. In this work, we combine archival observations with new data from ALMA and APEX, covering a range of sulfur-bearing species/isotopologs. We interpret this data using the DALI thermo-chemical code, for which our model is highly refined and disk-specific. We find that volatile sulfur is heavily depleted from the cosmic value by a factor of 1000, with a disk-averaged abundance of S/H = 1e-8. We show that the gas-phase sulfur abundance varies radially by 3 orders of magnitude, with the highest abundances inside the inner dust ring and coincident with the outer dust ring at 150 to 230 au. Extracting chemical abundances from our models, we find OCS, H2CS, and CS to be the dominant molecular carriers in the gas phase. We also infer the presence of a substantial OCS ice reservoir. We relate our results to the potential atmospheric composition of planets in HD 100546, and the wider exoplanet population, Comment: Accepted in MNRAS (20/12/2023)

(Exo-)planets inherit their budget of chemical elements from a protoplanetary disk. The disk temperature determines the phase of each chemical species, which sets the composition of solids and gas available for planet formation. We investigate how gap structures, which are widely seen by recent disk observations, alter the thermal and chemical structure of a disk. Planet-disk interaction is a leading hypothesis of gap formation and so such changes could present a feedback that planets have on planet-forming material. Both the planet gap-opening process and the disk thermal structure are well studied individually, but how the gap-opening process affects disk thermal structure evolution remains an open question. We develop a new modelling method by iterating hydrodynamical and radiative transfer simulations to explore the gap-opening feedback on disk thermal structure. We carry out parameter studies by considering different planet locations rp and planet masses Mp. We find that for the same rp and Mp, our iteration method predicts a wider and deeper gap than the non-iteration method. We also find that the inner disk and gap temperature from the iteration method can vary strongly from the non-iteration or disk without planets, which can further influence dust-trap conditions, iceline locations, and distribution of various ices, such as H2O, CO2, and CO on large dust grains ("pebbles"). Through that, a gap-opening planet can complicate the canonical picture of the non-planet disk C/O ratio and influence the composition of the next generation of planetesimals and planets., Comment: 17 pages, 15 figures, accepted for publication in MNRAS

The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and between disks. This is mostly tied to the icelines of major volatile carriers such as CO and H2O. Using ALMA observations of CS and SO, we have unearthed evidence for an entirely novel type of C/O variation in the protoplanetary disk around HD 100546: an azimuthal variation from a typical, oxygen-dominated ratio (C/O=0.5) to a carbon-dominated ratio (C/O>1.0). We show that the spatial distribution and peculiar line kinematics of both CS and SO molecules can be well-explained by azimuthal variations in the C/O ratio. We propose a shadowing mechanism that could lead to such a chemical dichotomy. Our results imply that tracing the formation history of giant exoplanets using their atmospheric C/O ratios will need to take into account time-dependent azimuthal C/O variations in a planet's accretion zone., Comment: Accepted in Nature Astronomy

KELT-9b is an ultra-hot Jupiter observed to be undergoing extreme mass loss. Its A0-type host star has a radiative envelope, which makes its surface layers prone to retaining recently accreted material. To search for potential signs of planetary material polluting the stellar surface, we carry out the most comprehensive chemical characterisation of KELT-9 to-date. New element detections include Na and Y, which had previously been detected in the ultra-hot Jupiter but not studied in the star; these detections complete the set of nine elements measured in both star and planet. In comparing KELT-9 with similar open cluster stars we find no strong anomalies. This finding is consistent with calculations of photospheric pollution accounting for stellar mixing and using observationally estimated KELT-9b mass loss rates. We also rule out recent, short-lived intensive mass transfer such as the stellar ingestion of an Earth-mass exomoon., Comment: 7 pages, 7 figures, accepted for publication in MNRAS

The nitrogen chemical evolution during star and planet formation is still not fully understood. Ammonia (NH$_3$) is a key specie in the understanding of the molecular evolution in star-forming clouds and nitrogen isotope fractionation. In this paper, we present high spatial resolution observations of multiple emission lines of NH$_3$ toward the protobinary system NGC1333 IRAS4A with Karl G. Jansky Very Large Array (VLA). We spatially resolved the binary (hereafter 4A1 and 4A2) and detected compact emission of NH$_3$ transitions with high excitation energies ($\gtrsim$100 K) from the vicinity of the protostars, indicating the NH$_3$ ice has sublimated at the inner hot region. The NH$_3$ column density is estimated to be $\sim 10^{17}-10^{18}$ cm$^{-2}$. We also detected two NH$_2$D transitions, allowing us to constrain the deuterium fractionation of ammonia. The NH$_2$D/NH$_3$ ratios are as high as $\sim 0.3-1$ in both 4A1 and 4A2. From the comparisons with the astrochemical models in the literature, the high NH$_2$D/NH$_3$ ratios suggest that the formation of NH$_3$ ices mainly started in the prestellar phase after the formation of bulk water ice finished, and that the primary nitrogen reservoir in the star-forming cloud could be atomic nitrogen (or N atoms) rather than nitrogen-bearing species such as N$_2$ and NH$_3$. The implications on the physical properties of IRAS4A cores are discussed as well., Comment: 20 pages, 7 figures, accepted for publication in ApJ

Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming planets. We report high angular resolution ALMA Band 7 observations of sulphur monoxide (SO) in the nearby planet-hosting disk around Herbig star HD 100546. SO is rarely detected in evolved protoplanetary disks but in other environments, it is most often utilised as a tracer of shocks. The SO emission from the HD 100546 disk is primarily originating from gas within the approx. 20 au mm-dust cavity and shows a clear azimuthal brightness asymmetry of a factor of 2. In addition, we see a significant difference in the line profile shape when comparing these new Cycle 7 data to Cycle 0 data of the same SO transitions. We discuss the different physical/chemical mechanisms that could be responsible for this asymmetry and time variability including disk winds, disk warps, and a shock triggered by a (forming) planet. We propose that the SO is enhanced in the cavity due to the presence of a giant planet. The SO asymmetry complements evidence for hot circumplanetary material around the giant planet HD 100546 c traced via CO ro-vibrational emission. This work sets the stage for further observational and modelling efforts to detect and understand the chemical imprint of a forming planet on its parent disk., Comment: Accepted A&A 24th October

Measuring the amount of gas and dust in protoplanetary disks is a key challenge in planet formation studies. Here we provide a new set of dust depletion factors and relative mass surface densities of gas and dust for the innermost regions of a sample of protoplanetary disks. We do this by combining stellar theory with observed refractory element abundances in both disk hosts and open cluster stars. Our results are independent of, and complementary to, those obtained from spatially resolved disk observations., Comment: 1 table

The composition of forming planets is strongly affected by the protoplanetary disc's thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps - regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly understood until now. In this work, we use the radiative transfer code MCMax to model the 2D dust thermal structure with individual gaps corresponding to planets with masses of 0.1 M$_J$ - 5 M$_J$ and orbital radii of 3, 5, and 10 AU. Low dust opacity in the gap allows radiation to penetrate deeper into the disc and warm the midplane by up to 16 K, but only for gaps located in the region of the disc where stellar irradiation is the dominant source of heat (here, a$\gtrsim$4 AU). In viscously-heated regions (a$\lesssim$4 AU), the midplane of the gap is relatively cooler by up to 100 K. Outside of the gap, broad radial oscillations in heating and cooling are present due to changes in disc flaring. These thermal features affect local segregation of volatile elements (H$_2$O, CH$_4$, CO2, CO) between the dust and gas. We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 AU towards the star and 4.3 AU towards the midplane. While quantitative predictions of iceline deviations will require more sophisticated models which include transport and sublimation/condensation kinetics, our results provide evidence that planet-induced iceline variations represent a potential feedback from the planet onto the composition of material it is accreting., Comment: 11 pages, 5 figures. Accepted to MNRAS, published 20 April 2023

Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds., Comment: 39 pages, 7 figures; under review for publication as a chapter in Protostars and Planets VII (Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura)

Context: Few exoplanets around hot stars with radiative envelopes have been discovered, although new observations from the TESS mission are improving this. Stars with radiative envelopes have little mixing at their surface, and thus their surface abundances provide a sensitive test case for a variety of processes including potentially star-planet interactions. Atomic diffusion is particularly important in these envelopes, producing chemically peculiar objects such as Am and HgMn stars. Aims: An exoplanet candidate around the B6 star HD 235349 was identified by TESS. Here we determine the nature of this transiting object and identify possible chemical peculiarities in the star. Methods: HD 235349 was observed using the long-slit spectrograph at Tartu Observatory, as well as photometrically by the TESS mission. The spectra were modeled to determine stellar parameters and chemical abundances. The photometric light curve was then analyzed in the context of the stellar parameters to determine properties of the transiting object. Results: We find the transiting object is a low-mass stellar companion, not a planet. However, the primary of this eclipsing binary is a rare type of chemically peculiar star. A strong overabundance of P is found with overabundances of Ne and Nd, and mild overabundances of Ti and Mn, while He is mildly underabundant. There is also clear evidence for vertical stratification of P in the atmosphere of the star. The lack of Hg and weak Mn overabundance suggests that this is not a typical HgMn star. It may be in the class of helium-weak phosphorus-gallium (He-weak PGa) stars, or an intermediate between these two classes. Conclusions: We show that HD 235349 is a rare type of chemically peculiar star (He-weak PGa) in an eclipsing binary system with a low-mass stellar companion. This appears to be the first He-weak PGa star discovered in an eclipsing binary., Comment: Accepted for publication in A&A. 15 pages (corrected minor typos, two references modified)

The goal of the Ariel space mission is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. The planetary bulk and atmospheric compositions bear the marks of the way the planets formed: Ariel's observations will therefore provide an unprecedented wealth of data to advance our understanding of planet formation in our Galaxy. A number of environmental and evolutionary factors, however, can affect the final atmospheric composition. Here we provide a concise overview of which factors and effects of the star and planet formation processes can shape the atmospheric compositions that will be observed by Ariel, and highlight how Ariel's characteristics make this mission optimally suited to address this very complex problem., Comment: 67 pages, 11 figures, 1 table. Postprint of the article to appear as open access in the Ariel Special Issue by the journal Experimental Astronomy

Quantifying the composition of the material in protoplanetary disks is paramount to determining the potential for exoplanetary systems to produce and support habitable environments. A key complex organic molecule (COM) to detect is methanol (CH3OH). CH3OH primarily forms at low temperatures via the hydrogenation of CO ice on the surface of icy dust grains and is a necessary basis for the formation of more complex species like amino acids and proteins. We report the detection of CH3OH in a disk around a young, luminous A-type star HD100546. This disk is warm and therefore does not host a significant CO ice reservoir. We argue that the CH3OH cannot form in situ, and hence, this disk has likely inherited COMs rich ice from an earlier cold dark cloud phase. This is strong evidence that at least some of the organic material survives the disk formation process and can then be incorporated into forming planets, moons and comets. Therefore, crucial pre-biotic chemical evolution already takes place in dark star-forming clouds., Comment: subm. to Nature Astronomy

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

The recently adopted Ariel ESA mission will measure the atmospheric composition of a large number of exoplanets. This information will then be used to better constrain planetary bulk compositions. While the connection between the composition of a planetary atmosphere and the bulk interior is still being investigated, the combination of the atmospheric composition with the measured mass and radius of exoplanets will push the field of exoplanet characterisation to the next level, and provide new insights of the nature of planets in our galaxy. In this white paper, we outline the ongoing activities of the interior working group of the {\it Ariel} mission, and list the desirable theoretical developments as well as the challenges in linking planetary atmospheres, bulk composition and interior structure., Comment: Accepted for publication in Experimental Astronomy

White dwarfs that have accreted planetary bodies are a powerful probe of the bulk composition of exoplanetary material. In this paper, we present a Bayesian model to explain the abundances observed in the atmospheres of 202 DZ white dwarfs by considering the heating, geochemical differentiation, and collisional processes experienced by the planetary bodies accreted, as well as gravitational sinking. The majority (>60%) of systems are consistent with the accretion of primitive material. We attribute the small spread in refractory abundances observed to a similar spread in the initial planet-forming material, as seen in the compositions of nearby stars. A range in Na abundances in the pollutant material is attributed to a range in formation temperatures from below 1,000K to higher than 1,400K, suggesting that pollutant material arrives in white dwarf atmospheres from a variety of radial locations. We also find that Solar System-like differentiation is common place in exo-planetary systems. Extreme siderophile (Fe, Ni or Cr) abundances in 8 systems require the accretion of a core-rich fragment of a larger differentiated body to at least a 3sigma significance, whilst one system shows evidence that it accreted a crust-rich fragment. In systems where the abundances suggest that accretion has finished (13/202), the total mass accreted can be calculated. The 13 systems are estimated to have accreted masses ranging from the mass of the Moon to half that of Vesta. Our analysis suggests that accretion continues for 11Myrs on average., Comment: 16 pages, MNRAS, in press

Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at $15-20$ au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ comes from a layer near the surface that extends no deeper than $z/r=0.25$. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry., Comment: 15 pages, 8 figures, Accepted in ApJ

The composition of giant planets is imprinted by their migration history and the compositional structure of their hosting disks. Studies in recent literature investigate how the abundances of C and O can constrain the formation pathways of giant planets forming within few tens of au from the star. New ALMA observations, however, suggest planet-forming regions possibly extending to hundreds of au. We explore the implications of these wider formation environments through n-body simulations of growing and migrating giant planets embedded in planetesimal disks, coupled with a compositional model of the protoplanetary disk where volatiles are inherited from the molecular cloud and refractories are calibrated against extrasolar and Solar System data. We find that the C/O ratio provides limited insight on the formation pathways of giant planets that undergo large-scale migration. This limitation can be overcome thanks to nitrogen and sulphur. Jointly using the C/N, N/O and C/O ratios breaks any degeneracy in the formation and migration tracks of giant planets. The use of elemental ratios normalized to the respective stellar ratios supplies additional information on the nature of giant planets, thanks to the relative volatility of O, C and N in disks. When the planetary metallicity is dominated by the accretion of solids C/N* $>$ C/O* $>$ N/O* (* denoting this normalized scale), otherwise N/O* $>$ C/O* $>$ C/N*. The S/N ratio provides an additional independent probe into the metallicity of giant planets and their accretion of solids., Comment: 27 pages, 10 figures and 7 tables. Accepted version of the article published on The Astrophysical Journal (includes 1 extra figure and additional text with respect to previous preprint)

Recent observations have revealed that most proto-planetary discs show a pattern of bright rings and dark gaps. However, most of the high-resolution observations have focused only on the continuum emission. In this Paper we present high-resolution ALMA band 7 (0.89mm) observations of the disc around the star CI Tau in the $^{12}$CO & $^{13}$CO $J=3$-2 and CS $J=7$-6 emission lines. Our recent work demonstrated that the disc around CI Tau contains three gaps and rings in continuum emission, and we look for their counterparts in the gas emission. While we find no counterpart of the third gap and ring in $^{13}$CO, the disc has a gap in emission at the location of the second continuum ring (rather than gap). We demonstrate that this is mostly an artefact of the continuum subtraction, although a residual gap still remains after accounting for this effect. Through radiative transfer modelling we propose this is due to the inner disc shadowing the outer parts of the disc and making them colder. This raises a note of caution in mapping high-resolution gas emission lines observations to the gas surface density - while possible, this needs to be done carefully. In contrast to $^{13}$CO, CS emission shows instead a ring morphology, most likely due to chemical effects. Finally, we note that $^{12}$CO is heavily absorbed by the foreground preventing any morphological study using this line., Comment: 17 pages, 17 figures; accepted for publication in MNRAS

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., Comment: 28 pages, 12 main figures, 14 appendix figures. Accepted to ApJ

HD\,139614 is known to be a $\sim$14-Myr-old, possibly pre-main-sequence star in the Sco-Cen OB association in the Upper Centaurus-Lupus subgroup, with a slightly warped circumstellar disc containing ring structures hinting at one or more planets. The star's chemical abundance pattern is metal-deficient except for volatile elements, which places it in the $\lambda$ Boo class and suggests it has recently accreted gas-rich but dust-poor material. We identify seven dipole and four radial pulsation modes among its $\delta$ Sct pulsations using the TESS light curve and an \'echelle diagram. Precision modelling with the MESA stellar evolution and GYRE stellar oscillation programs confirms it is on the pre-main sequence. Asteroseismic, grid-based modelling suggests an age of $10.75\pm0.77$ Myr, a mass of $1.52\pm0.02$ M$_{\odot}$, and a global metal abundance of $Z=0.0100\pm0.0010$. This represents the first asteroseismic determination of the bulk metallicity of a $\lambda$ Boo star. The precise age and metallicity offer a benchmark for age estimates in Upper Centaurus--Lupus, and for understanding disc retention and planet formation around intermediate-mass stars., Comment: Figures make use of colour. Accepted in MNRAS

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's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" 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.

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".2-0".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., Comment: Accepted for publication in The Astronomical Journal

We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent radius of ~100 au. The VLA observations detect both the disk in dust continuum and free-free emission extended along the jet direction. The ALMA observations of molecular lines (H$_2$CO, SO, CH$_3$OH, $^{13}$CO, C$^{18}$O, NS, and H$^{13}$CN) reveal rotation of the apparent disk surrounding HOPS-370 orthogonal to the jet/outflow direction. We fit radiative transfer models to both the dust continuum structure of the disk and molecular line kinematics of the inner envelope and disk for the H$_2$CO, CH$_3$OH, NS, and SO lines. The central protostar mass is determined to be $\sim$2.5 M_sun with a disk radius of $\sim$94~au, when fit using combinations of the H$_2$CO, CH$_3$OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution (SED) yields a disk mass of 0.035~M$_{\odot}$ (inferred dust+gas) and a dust disk radius of 62~au, thus the dust disk may have a smaller radius than the gas disk, similar to Class II disks. In order to explain the observed luminosity with the measured protostar mass, HOPS-370 must be accreting at a rate between 1.7 and 3.2$\times$10$^{-5}$~M$_{\odot}$~yr$^{-1}$., Comment: Accepted to ApJ; 51 pages, 12 Figures, 7 Tables

Exocomets are small bodies releasing gas and dust which orbit stars other than the Sun. Their existence was first inferred from the detection of variable absorption features in stellar spectra in the late 1980s using spectroscopy. More recently, they have been detected through photometric transits from space, and through far-IR/mm gas emission within debris disks. As (exo)comets are considered to contain the most pristine material accessible in stellar systems, they hold the potential to give us information about early stage formation and evolution conditions of extra Solar Systems. In the Solar System, comets carry the physical and chemical memory of the protoplanetary disk environment where they formed, providing relevant information on processes in the primordial solar nebula. The aim of this paper is to compare essential compositional properties between Solar System comets and exocomets. The paper aims to highlight commonalities and to discuss differences which may aid the communication between the involved research communities and perhaps also avoid misconceptions. Exocomets likely vary in their composition depending on their formation environment like Solar System comets do, and since exocomets are not resolved spatially, they pose a challenge when comparing them to high fidelity observations of Solar System comets. Observations of gas around main sequence stars, spectroscopic observations of "polluted" white dwarf atmospheres and spectroscopic observations of transiting exocomets suggest that exocomets may show compositional similarities with Solar System comets. The recent interstellar visitor 2I/Borisov showed gas, dust and nuclear properties similar to that of Solar System comets. This raises the tantalising prospect that observations of interstellar comets may help bridge the fields of exocomet and Solar System comets., Comment: 25 pages, 3 figures. To be published in PASP. This paper is the product of a workshop at the Lorentz Centre in Leiden, the Netherlands

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase., Comment: 107 pages, 21 Figures, 11 Tables, accepted to ApJ. Version with all source figures: https://www.cv.nrao.edu/~jtobin/Orion-disks-fullfigs.pdf Reduced data available from https://dataverse.harvard.edu/dataverse/VANDAMOrion

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them., Comment: 32 Pages, 10 Figures, 5 Tables, accepted to ApJ

Sulfur is one of the most abundant elements in the Universe, with important roles in astro-, geo-, and biochemistry. Its main reservoirs in planet-forming disks have previously eluded detection: gaseous molecules only account for $<1\,$\% of total elemental sulfur, with the rest likely in either ices or refractory minerals. Mechanisms such as giant planets can filter out dust from gas accreting onto disk-hosting stars. For stars above 1.4 solar masses, this leaves a chemical signature on the stellar photosphere that can be used to determine the fraction of each element that is locked in dust. Here, we present an application of this method to sulfur, zinc, and sodium. We analyse the accretion-contaminated photospheres of a sample of young stars and find $(89\pm8)\,$\% of elemental sulfur is in refractory form in their disks. The main carrier is much more refractory than water ice, consistent with sulfide minerals such as FeS., Comment: 11 pages, 3 figures, accepted for publication in ApJ, correspondence welcome

Gas giants are the most readily detectable exoplanets but fundamental questions about their system architectures, formation, migration, and atmospheres have been unanswerable with the current generation of ground- and space-based facilities. The dominant techniques to detect and characterize giant planets $-$ radial velocities, transits, direct imaging, microlensing, and astrometry $-$ are each isolated to a limited range of planet masses, separations, ages, and temperatures. These windows into the arrangement and physical properties of giant planets have spawned new questions about the timescale and location of their assembly; the distributions of planet mass and orbital separation at young and old ages; the composition and structure of their atmospheres; and their orbital and rotational angular momentum architectures. The ELTs will address these questions by building bridges between these islands of mass, orbital distance, and age. The angular resolution, collecting area, all-sky coverage, and novel instrumentation suite of these facilities are needed to provide a complete map of the orbits and atmospheric evolution of gas giant planets (0.3$-$10 $M_\mathrm{Jup}$) across space (0.1$-$100 AU) and time (1 Myr to 10 Gyr). This white paper highlights the scientific potential of the GMT and TMT to address these outstanding questions, with a particular focus on the role of direct imaging and spectroscopy of large samples of giant planets that will soon be made available with $Gaia$., Comment: White paper for the Astro2020 decadal survey

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade., Comment: White paper submitted to the Astro2020 decadal survey

Despite being $>10$Myr, there are $\sim$10 debris discs with as much CO gas as in protoplanetary discs. Such discs have been assumed to be "hybrid", i.e., with secondary dust but primordial gas. Here we show that both the dust and gas in such systems could instead be secondary, with the high CO content caused by accumulation of neutral carbon (C$^0$) that shields CO from photodissociating; i.e., these could be "shielded secondary discs". New ALMA observations are presented of HD131835 that detect $\sim 3 \times 10^{-3}$ M$_\oplus$ of C$^0$, the majority 40-200au from the star, in sufficient quantity to shield the previously detected CO. A simple semi-analytic model for the evolution of CO, C and O originating in a volatile-rich planetesimal belt shows how CO shielding becomes important when the viscous evolution is slow (low $\alpha$ parameter) and/or the CO production rate is high. Shielding by C$^0$ may also cause the CO content to reach levels at which CO self-shields, and the gas disc may become massive enough to affect the dust evolution. Application to the HD 131835 observations shows these can be explained if $\alpha \sim 10^{-3}$; an inner cavity in C$^0$ and CO may also mean the system has yet to reach steady state. Application to other debris discs with high CO content finds general agreement for $\alpha=10^{-3}$ to $0.1$. The shielded secondary nature of these gas discs can be tested by searching for C$^0$, as well as CN, N$_2$ and CH$^{+}$, which are also expected to be shielded by C$^0$., Comment: accepted to MNRAS

We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disc; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exo-planets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius.Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage., Comment: Accepted on Astrophysical Journal Letters; 7 pages, 2 figures

Protoplanetary disks, debris disks, and disrupted or evaporating planets can all feed accretion onto stars. The photospheric abundances of such stars may then reveal the composition of the accreted material. This is especially likely in B to mid-F type stars, which have radiative envelopes and hence less bulk--photosphere mixing. We present a theoretical framework (\texttt{CAM}) considering diffusion, rotation, and other stellar mixing mechanisms, to describe how the accreted material interacts with the bulk of the star. This allows the abundance pattern of the circumstellar material to be calculated from measured stellar abundances and parameters ($v_{\rm rot}$, $T_{\rm eff}$). We discuss the \lboo\ phenomenon and the application of \texttt{CAM} on stars hosting protoplanetary disks (HD~100546, HD~163296), debris disks (HD~141569, HD~21997), and evaporating planets (HD~195689/KELT-9)., Comment: 17 pages, 9 figures. Published in MNRAS. Updated with discussion of WD pollution

Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5$\sigma$ detection in the stacked line profile. The emission is compact but is spectrally resolved and the line profile has two components. One component peaks at the source velocity and the other is blue-shifted by 5 km s$^{-1}$. The kinematics and spatial distribution of the SO emission are not consistent with that expected from a purely Keplerian disk. We detect additional blue-shifted emission that we attribute to a disk wind. The disk component was simulated using LIME and a physical disk structure. The disk emission is asymmetric and best fit by a wedge of emission in the north east region of the disk coincident with a `hot-spot' observed in the CO $J=3-2$ line. The favoured hypothesis is that a possible inner disk warp (seen in CO emission) directly exposes the north-east side of the disk to heating by the central star, creating locally the conditions to launch a disk wind. Chemical models of a disk wind will help to elucidate why the wind is particularly highlighted in SO emission and whether a refractory source of sulphur is needed. An alternative explanation is that the SO is tracing an accretion shock from a circumplanetary disk associated with the proposed protoplanet embedded in the disk at 50 au. We also report a non-detection of SO in the protoplanetary disk around HD 97048., Comment: Accepted for publication in A&A on 12/12/17

We present a complete census of the 14N/15N isotopic ratio in the most abundant N-bearing molecules towards the cold envelope of the protocluster OMC-2 FIR4, the best known Sun progenitor. To this scope, we analysed the unbiased spectral survey obtained with the IRAM-30m telescope at 3mm, 2mm and 1mm. We detected several lines of CN, HCN, HNC, HC3N, N2H+, and their respective 13C and 15N isotopologues. The lines relative fluxes are compatible with LTE conditions and moderate line opacities have been corrected via a Population Diagram method or theoretical relative intensity ratios of the hyperfine structures. The five species lead to very similar 14N/15N isotopic ratios, without any systematic difference between amine and nitrile bearing species as previously found in other protostellar sources. The weighted average of the 14N/15N isotopic ratio is 270 +/- 30. This 14N/15N value is remarkably consistent with the [250-350] range measured for the local galactic ratio but significantly differs from the ratio measured in comets (around 140). High-angular resolution observations are needed to examine whether this discrepancy is maintained at smaller scales. In addition, using the CN, HCN and HC3N lines, we derived a 12C/13C isotopic ratio of 50 +/- 5., Comment: Accepted for publication in ApJ ; 19 pages, 5 tables, 12 figures

Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tracer of the disk gas mass, specifically the effect of gas mass on the HD FIR emission and its sensitivity to the vertical structure. Deuterium chemistry reactions relevant for HD were implemented in the thermochemical code DALI and models were run for a range of disk masses and vertical structures. The HD J=1-0 line intensity depends directly on the gas mass through a sublinear power law relation with a slope of ~0.8. Assuming no prior knowledge about the vertical structure of a disk and using only the HD 1-0 flux, gas masses can be estimated to within a factor of 2 for low mass disks (M$_{\rm disk} < 10^{-3}$ M$_\odot$). For more massive disks, this uncertainty increases to more than an order of magnitude. Adding the HD 2-1 line or independent information about the vertical structure can reduce this uncertainty to a factor of ~3 for all disk masses. For TW Hya, using the radial and vertical structure from Kama et al. 2016b the observations constrain the gas mass to $6\cdot10^{-3}$ M$_\odot$ < M$_{\rm disk} < 9\cdot10^{-3}$ M$_\odot$. Future observations require a 5$\sigma$ sensitivity of $1.8\cdot10^{-20}$ W m$^{-2}$ ($2.5\cdot10^{-20}$ W m$^{-2}$) and a spectral resolving power R > 300 (1000) to detect HD 1-0 (HD 2-1) for all disk masses above $10^{-5}$ M$_\odot$ with a line-to-continuum ratio > 0.01. These results show that HD can be used as an independent gas mass tracer with a relatively low uncertainty and should be considered as an important science goal for future FIR missions., Comment: 18 pages, 13 figures, accepted in A&A

High-resolution imaging reveals a large morphological variety of protoplanetary disks. To date, no constraints on their global evolution have been found from this census. An evolutionary classification of disks was proposed based on their IR spectral energy distribution, with the Group I sources showing a prominent cold component ascribed to an earlier stage of evolution than Group II. Disk evolution can be constrained from the comparison of disks with different properties. A first attempt of disk taxonomy is now possible thanks to the increasing number of high-resolution images of Herbig Ae/Be stars becoming available. Near-IR images of six Group II disks in scattered light were obtained with VLT/NACO in Polarimetric Differential Imaging, which is the most efficient technique to image the light scattered by the disk material close to the stars. We compare the stellar/disk properties of this sample with those of well-studied Group I sources available from the literature. Three Group II disks are detected. The brightness distribution in the disk of HD163296 indicates the presence of a persistent ring-like structure with a possible connection with the CO snowline. A rather compact (less than 100 AU) disk is detected around HD142666 and AK Sco. A taxonomic analysis of 17 Herbig Ae/Be sources reveals that the difference between Group I and Group II is due to the presence or absence of a large disk cavity (larger than 5 AU). There is no evidence supporting the evolution from Group I to Group II. Group II are not evolved version of the Group I. Within the Group II disks, very different geometries (both self-shadowed and compact) exist. HD163296 could be the primordial version of a typical Group I. Other Group II, like AK Sco and HD142666, could be smaller counterpart of Group I unable to open cavities as large as those of Group I., Comment: 16 pages, 7 figures, published by A&A

[Abridged] Snowlines in protoplanetary disks play an important role in planet formation and composition. Since the CO snowline is difficult to observe directly with CO emission, its location has been inferred in several disks from spatially resolved ALMA observations of DCO+ and N2H+. N2H+ is considered to be a good tracer of the CO snowline based on astrochemical considerations predicting an anti-correlation between N2H+ and gas-phase CO. In this work, the robustness of N2H+ as a tracer of the CO snowline is investigated. A simple chemical network is used in combination with the radiative transfer code LIME to model the N2H+ distribution and corresponding emission in the disk around TW Hya. The assumed CO and N2 abundances, corresponding binding energies, cosmic ray ionization rate, and degree of large-grain settling are varied to determine the effects on the N2H+ emission and its relation to the CO snowline. For the adopted physical structure of the TW Hya disk and molecular binding energies for pure ices, the balance between freeze-out and thermal desorption predicts a CO snowline at 19 AU, corresponding to a CO midplane freeze-out temperature of 20 K. A model with a total, i.e. gas plus ice, CO abundance of 3e-6 with respect to H2 fits the position of the emission peak observed by Qi et al. 2013 for the TW Hya disk. However, the relationship between N2H+ and the CO snowline is more complicated than generally assumed: for the investigated parameters, the N2H+ column density peaks at least 5 AU outside the CO snowline. Moreover, the N2H+ emission can peak much further out, as far as ~50 AU beyond the snowline. Hence, chemical modeling, as done here, is necessary to derive a CO snowline location from N2H+ observations., Comment: 16 pages, 3 tables, 15 figures. Accepted for publication in A&A

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 \sigma$), reaching a peak signal-to-noise of $5.5\sigma$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation., Comment: 14 pages, 4 figures, 1 table, published online in ApJL on 13th May 2016

Despite intensive studies of protoplanetary disks, there is still no reliable way to determine their total mass and their surface density distribution, quantities that are crucial for describing both the structure and the evolution of disks up to the formation of planets. The goal of this work is to use less abundant CO isotopologues, whose detection is routine for ALMA, to infer the gas mass of disks. Isotope-selective effects need to be taken into account in the analysis, because they can significantly modify CO isotopologues line intensities. CO isotope-selective photodissociation has been implemented in the physical-chemical code DALI and 800 disk models have been run for a range of disk and stellar parameters. Dust and gas temperature structures have been computed self-consistently, together with a chemical calculation of the main species. Both disk structure and stellar parameters have been investigated. Total fluxes have been ray-traced for different CO isotopologues and for various transitions for different inclinations. A combination of 13CO and C18O total intensities allows inference of the total disk mass, although with non-negligible uncertainties. These can be overcome by employing spatially resolved observations, i.e. the disk's radial extent and inclination. Comparison with parametric models shows differences at the factor of a few level, especially for extremely low and high disk masses. Finally, total line intensities for different CO isotopologue and for various low-J transitions are provided and are fitted to simple formulae. The effects of a lower gas-phase carbon abundance and different gas/dust ratios are investigated as well, and comparison with other tracers is made. Disk masses can be determined within a factor of a few by comparing CO isotopologue lines observations with the simulated line fluxes, modulo the uncertainties in the volatile elemental abundances., Comment: 20 pages, 15 figures, A&A accepted

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., Comment: (8 pages, 4 figures, accepted by Astronomy & Astrophysics)

Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps." The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths., Comment: 6 pages, 4 figures, accepted for publication in ApJL, fixed typos