Search

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., 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

Tartu Observatory telescopes offer unique guaranteed access to objects in the Northern hemisphere. The observational facilities include a 1.5- m and 0.6-m classic Cassegrain reflectors, and 0.31-m remotely controllable telescope . The 1.5 m telescope is currently equipped with a long-slit Cassegrain spectrograph used for stellar characterisation. Historically, the objects of interest have been massive stars but we are now developing new research direction and expanding the list of targets to exoplanet- and disk-hosting stars. We have started evaluating our capabilities to characterise host stars spectroscopically to determine their parameters and composition. In 2020, we carried out a pilot study of a TESS candidate planet host, which we found to have a rare, strong chemical peculiarity [1]. This also allowed us to prepare our tools, workflow, and end-to-end analysis. We are also contributing to the European Space Agency Ariel space mission by offering stellar activity monitoring. The 0.6-m and 0.31-m telescopes are utilised for photometric measurements and the 0.31-m one in particular has been a workhorse for exoplanet transit monitoring. Since 2020, we have made significant preparations to develop and prove our transit observation capabilities: we have observed more than 70 transit light curves. About half of them have been submitted to ExoClock to contribute to Ariel mission planning. Concerning future upgrades, Tartu Observatory will have new instruments by the middle of 2023. The upgrades include procuring a medium resolution echelle spectrograph (projected bandwidth 390 nm to 750 nm, R= 25 000) and new photometer (Johnson-Cousins BVRI and SDSS filters) for the 1.5-m telescope, which will not only enhance our capabilities in both spectroscopic and photometric data retrieval of host stars. In addition, a new remote control system of the telescope will be installed and improvements on instrumentation for the 0.6-m and 0.31-m photometric telescopes will be made. This presentation will give an overview of our facilities, and of current and future spectroscopic and photometric capabilities. References: “A rare phosphorus-rich star in an eclipsing binary from TESS”, Colin P. et al., A&A 658 A105 (2022), DOI: 10.1051/0004-6361/202142124

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. This presentation provides an overall summary of the science and instrument design for Ariel and presents the many activities that the Ariel team have planned to engage the science community at large and the public prior to launch. These include the Ariel Dry-Run program and citizen-science programs such as ExoClock and the Ariel Data Challenges.

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 Ariel mission, and list the desirable theoretical developments as well as the challenges in linking planetary atmospheres, bulk composition and interior structure., Experimental Astronomy, 53, ISSN:0922-6435, ISSN:1572-9508

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made., 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. © The Author(s) 2021., D.T., S.F., S.M., E.S., and A.N. acknowledge the support of the Italian Space Agency (ASI) through the ASI-INAF contract 2018-22-HH.0. D.T., C.C., D.F., and L.P. acknowledge the support of the PRIN-INAF 2016 “The Cradle of Life - GENESIS-SKA (General Conditions in Early Planetary Systems for the rise of life with SKA”. D.T., S.F., S.M. D.F, J.M., F.O., P.W. acknowledge the support of the Italian National Institute of Astrophysics (INAF) through the INAF Main Stream project “Ariel and the astrochemical link between circumstellar discs and planets” (CUP: C54I19000700005). S.M. acknowledges support from the European Research Council via the Horizon 2020 Framework Programme ERC Synergy “ECOGAL” Project GA-855130. M.K. acknowledges funding by the University of Tartu ASTRA project 2014-2020.4.01.16-0029 KOMEET “Benefits for Estonian Society from Space Research and Application”, financed by the EU European Regional Development Fund. J.M.D.K. gratefully acknowledges funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through an Emmy Noether Research Group (grant number KR4801/1-1) and the DFG Sachbeihilfe (grant number KR4801/2-1), as well as from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme via the ERC Starting Grant MUSTANG (grant agreement number 714907). The research of O.P. is funded by the Royal Society, through Royal Society Dorothy Hodgkin Fellowship DH140243. M.P. thanks the support to NuGrid from STFC (through the University of Hull’s Consolidated Grant ST/R000840/1), and access to VIPER, the University of Hull High Performance Computing Facility. M.P. acknowledges the support from the ”Lendulet-2014” Program of the Hungarian Academy of Sciences (Hungary), from the ERC Consolidator Grant (Hungary) funding scheme (Project RADIOSTAR, G.A. n. 724560), by the National Science Foundation (NSF, USA) under grant No. PHY-1430152 (JINA Center for the Evolution of the Elements). M.P. also thanks the UK network BRIDGCE and the ChETEC COST Action (CA16117), supported by COST (European Cooperation in Science and Technology). M.I. thanks the support by JSPS KAKENHI 18H05439. C.D. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709), and the Group project Ref. PID2019-110689RB-I00/AEI/10.13039/501100011033.

Class III stars are those in star forming regions without large non-photospheric infrared emission, suggesting recent dispersal of their protoplanetary disks. We observed 30 class III stars in the 1-3 Myr Lupus region with ALMA at ${\sim}856��$m, resulting in 4 detections that we attribute to circumstellar dust. Inferred dust masses are $0.036{-}0.093M_\oplus$, ${\sim}1$ order of magnitude lower than any previous measurements; one disk is resolved with radius ${\sim}80$ au. Two class II sources in the field of view were also detected, and 11 other sources, consistent with sub-mm galaxy number counts. Stacking non-detections yields a marginal detection with mean dust mass ${\sim}0.0048M_\oplus$. We searched for gas emission from the CO J=3-2 line, and present its detection to NO Lup inferring a gas mass ($4.9 {\pm} 1.1$) ${\times}10^{-5} M_\oplus$ and gas-to-dust ratio $1.0{\pm}0.4$. Combining our survey with class II sources shows a gap in the disk mass distribution from $0.09{-}2M_\oplus$ for ${>}0.7M_\odot$ Lupus stars, evidence of rapid dispersal of mm-sized dust from protoplanetary disks. The class III disk mass distribution is consistent with a population model of planetesimal belts that go on to replenish the debris disks seen around main sequence stars. This suggests that planetesimal belt formation does not require long-lived protoplanetary disks, i.e., planetesimals form within ${\sim}$2 Myr. While all 4 class III disks are consistent with collisional replenishment, for two the gas and/or mid-IR emission could indicate primordial circumstellar material in the final stages of protoplanetary disk dispersal. Two class III stars without sub-mm detections exhibit hot emission that could arise from ongoing planet formation processes inside ${\sim}1$ au., 30 pages, 20 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

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. Submitted to MNRAS, preprint for reference at Exoplanets IV

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

This is a correction for 2021 AJ 161 38

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

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/

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., Accepted for publication in A&A. 15 pages (corrected minor typos, two references modified)

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., subm. to Nature 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., 16 pages, MNRAS, in press

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

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., 15 pages, 8 figures, Accepted in ApJ

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

We observed the K7 class III star NO Lup in an ALMA survey of the 1-3 Myr Lupus association and detected circumstellar dust and CO gas. Here we show that the J = 3-2 CO emission is both spectrally and spatially resolved, with a broad velocity width ${\sim}19$kms$^{-1}$ for its resolved size ${\sim}1''$ (${\sim}130$ au). We model the gas emission as a Keplerian disc, finding consistency, but only with a central mass of ${\sim}11M_{\odot}$, which is implausible given its spectral type and X-Shooter spectrum. A good fit to the data can also be found by modelling the CO emission as outflowing gas with a radial velocity ${\sim}22$kms$^{-1}$. We interpret NO Lup's CO emission as the first imaged class III circumstellar disc with outflowing gas. We conclude that the CO is continually replenished, but cannot say if this is from the break-up of icy planetesimals or from the last remnants of the protoplanetary disc. We suggest further work to explore the origin of this CO, and its higher than expected velocity in comparison to photoevaporative models., 6 pages, 4 figures. Accepted for publication in MNRAS following peer review (this is a pre-copyedited, author-produced PDF version)

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., 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

Formaldehyde (H2CO) is an important precursor to organics like methanol (CH3OH). It is important to understand the conditions that produce H2CO and prebiotic molecules during star and planet formation. H2CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold (≲20 K) dust grains. To understand which pathway dominates, gaseous H2CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and 12 cm−2 and 25 ± 2 K, (4.4 ± 0.3) × 1011 cm−2 for ortho- and para-H2CO, respectively, and an OPR of 2.49 ± 0.23. A radially resolved analysis shows that the observed H2CO emits mostly at rotational temperatures of 30–40 K, corresponding to a layer with z/R ≥ 0.25. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to 2.0 ± 0.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, led us to speculate that gas-phase formation is responsible for the observed H2CO across the TW Hya disk.

We present ALMA (0.87 and 1.3 mm) and VLA (9 mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L bol ∼ 314 L o˙) at ∼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 and 1.3 mm resolves 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 (H2CO, SO, CH3OH, 13CO, C18O, NS, and H13CN) 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 H2CO, CH3OH, NS, and SO lines. The central protostar mass is determined to be ∼2.5 M o˙ with a disk radius of ∼94 au, when fit using combinations of the H2CO, CH3OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution yields a disk mass of 0.035 M o˙ (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 of (1.7-3.2) × 10-5 M o˙ yr-1. © 2020. The American Astronomical Society. All rights reserved., The authors thank the anonymous referee for a constructive report that improved the quality of the manuscript. The authors wish to acknowledge useful discussions with L. Hartmann and J. Bae about the results. J.J.T. acknowledges support from grant AST-1814762 from the National Science Foundation and past support from the Homer L. Dodge Endowed Chair at the University of Oklahoma. P.D.S is supported by a National Science Foundation Astronomy & Astrophysics Postdoctoral Fellowship under Award No. 2001830. G.A., M.O., and A.K.D.-R. acknowledge financial support from the State Agency for Research of the Spanish MCIU through the AYA2017-84390-C2-1-R grant (co-funded by FEDER) and through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). Z.Y.L. is supported in part by NASA 80NSSC18K1095 and NSF AST-1716259. M.K. gratefully acknowledges funding by the University of Tartu ASTRA project 2014-2020.4.01.16-0029 KOMEET "Benefits for Estonian Society from Space Research and Application," financed by the EU European Regional Development Fund. M.L.R.H. acknowledges support from the Michigan Society of Fellows. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00041.S and ADS/JAO.ALMA#2017.1.00419.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. These results are based on observations obtained with the Apache Point Observatory 3.5m telescope, which is owned and operated by the Astrophysical Research Consortium (ARC). We wish to thank the APO 3.5m telescope operators in their assistance in obtaining these data. Access to the APO 3.5m telescope is funded by the University of Oklahoma and the Homer L. Dodge Endowed Chair. This research made use of APLpy, an open-source plotting package for Python hosted at.http://aplpy.github.com.This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration, 2013) http://www.astropy.org.

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)., 17 pages, 9 figures. Published in MNRAS. Updated with discussion of WD pollution

Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the field centred on J1407 in the 340 GHz (Band 7) continuum in order to determine the flux and astrometric location of the ring system relative to the star. We used the VLT/NACO camera to observe the J1407 system in March 2019 and to search for the central planetary mass object at thermal infrared wavelengths. We detect no point source at the expected location of J1407, and derive an upper limit $3\sigma$ level of $57.6~\mu\rm{Jy}$. There is a point source detected at an angular separation consistent with the expected location for a free-floating ring system that occulted J1407 in May 2007, with a flux of $89~\mu\rm{Jy}$ consistent with optically thin dust surrounding a massive substellar companion. At 3.8 microns with the NACO camera, we detect the star J1407 but no other additional point sources within 1.3 arcseconds of the star, with a lower bound on the sensitivity of $6M_{Jup}$ at the location of the ALMA source, and down to $4M_{Jup}$ in the sky background limit. The ALMA upper limit at the location of J1407 implies that a hypothesised bound ring system is composed of dust smaller than $1\rm{~mm}$ in size, implying a young ring structure. The detected ALMA source has multiple interpretations, including: (i) it is an unbound substellar object surrounded by warm dust in Sco-Cen with an upper mass limit of $6M_{Jup}$, or (ii) it is a background galaxy., Comment: Accepted for publication in A&A (6 pages, 5 figures). Reduced data and reduction scripts on GitHub at https://github.com/mkenworthy/j1407_alma_detection

Low mass stars might offer today the best opportunities to detect and characterise planetary systems, especially those harbouring close-in low mass temperate planets. Among those stars, TRAPPIST-1 is exceptional since it has seven Earth-sized planets, of which three could sustain liquid water on their surfaces. Here we present new and deep ALMA observations of TRAPPIST-1 to look for an exo-Kuiper belt which can provide clues about the formation and architecture of this system. Our observations at 0.88 mm did not detect dust emission, but can place an upper limit of 23$��$Jy if the belt is smaller than 4 au, and 0.15 mJy if resolved and 100 au in radius. These limits correspond to low dust masses of $\sim10^{-5}-10^{-2}$ $M_\oplus$, which are expected after 8 Gyr of collisional evolution unless the system was born with a $>20$ $M_\oplus$ belt of 100 km-sized planetesimals beyond 40 au or suffered a dynamical instability. This $20$ $M_\oplus$ mass upper limit is comparable to the combined mass in TRAPPIST-1 planets, thus it is possible that most of the available solid mass in this system was used to form the known planets. A similar analysis of the ALMA data on Proxima Cen leads us to conclude that a belt born with a mass $\gtrsim1$ $M_\oplus$ in 100 km-sized planetesimals could explain its putative outer belt at 30 au. We recommend that future characterisations of debris discs around low mass stars should focus on nearby and young systems if possible., accepted for publication in MNRAS, 7 pages, 5 figures

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time., Comment: Accepted for publication in A&A

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 report new interferometric images of cyclopropenylidene, c-C$_3$H$_2$, towards the young protocluster OMC-2 FIR\,4. The observations were performed at 82 and 85 GHz with the NOrthern Extended Millimeter Array (NOEMA) as part of the project Seeds Of Life In Space (SOLIS). In addition, IRAM-30m data observations were used to investigate the physical structure of OMC-2 FIR\,4. We find that the c-C$_3$H$_2$ gas emits from the same region where previous SOLIS observations showed bright HC$_5$N emission. From a non-LTE analysis of the IRAM-30m data, the c-C$_3$H$_2$ gas has an average temperature of $\sim$40K, a H$_2$ density of $\sim$3$\times$10$^{5}$~cm$^{-3}$, and a c-C$_3$H$_2$ abundance relative to H$_2$ of ($7\pm1$)$\times$10$^{-12}$. In addition, the NOEMA observations provide no sign of significant c-C$_3$H$_2$ excitation temperature gradients across the region (about 3-4 beams), with T$_{ex}$ in the range 8$\pm$3 up to 16$\pm$7K. We thus infer that our observations are inconsistent with a physical interaction of the OMC-2 FIR\,4 envelope with the outflow arising from OMC-2 FIR\,3, as claimed by previous studies. The comparison of the measured c-C$_3$H$_2$ abundance with the predictions from an astrochemical PDR model indicates that OMC-2 FIR\,4 is irradiated by a FUV field $\sim$1000 times larger than the interstellar one, and by a flux of ionising particles $\sim$4000 times larger than the canonical value of $1\times10^{-17}$~s$^{-1}$ from the Galaxy cosmic rays, which is consistent with our previous HC$_5$N observations. This provides an important and independent confirmation of other studies that one or more sources inside the OMC-2 FIR\,4 region emit energetic ($\geq10$~MeV) particles., Accepted for publication in ApJ - 14 pages, 10 figures, 4 tables

Cometary ices are formed during star and planet formation, and their molecular and elemental makeup can be related to the early solar system via the study of inter- and protostellar material. The first cometary abundance of the halogen element bromine (Br) was recently made available by the Rosetta mission. Its abundance in protostellar gas is thus far unconstrained, however. We set out to place the first observational constraints on the interstellar gas-phase abundance of bromine (Br). We further aim to compare the protostellar Br abundance with that measured by Rosetta in the ices of comet 67P/Churyumov-Gerasimenko. Archival Herschel data of Orion KL, Sgr B2(N), and NGC 6334I are examined for the presence of HBr and HBr$^{+}$ emission or absorption lines. A chemical network for modelling HBr in protostellar molecular gas is compiled to aid in the interpretation. HBr and HBr$^{+}$ were not detected towards any of our targets. However, in the Orion KL Hot Core, our upper limit on HBr/H$_{2}$O is a factor of ten below the ratio measured in comet 67P. This result is consistent with the chemical network prediction that HBr is not a dominant gas-phase Br carrier. Cometary HBr is likely predominantly formed in icy grain mantles which lock up nearly all elemental Br., Comment: Accepted for publication in A&A. 9 pages, 6 figures

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

We report on the discovery of correlations between dust and CO gas tracers of the 0.1-10 au region in planet-forming disks around young intermediate-mass stars. The abundance of refractory elements on stellar photospheres decreases as the location of hot CO gas emission recedes to larger disk radii, and as the near-infrared excess emission from hot dust in the inner disk decreases. The linked behavior between these observables demonstrates that the recession of infrared CO emission to larger disk radii traces an inner disk region where dust is being depleted. We also find that Herbig disk cavities have either low (~ 5-10 %) or high (~ 20-35 %) near-infrared excess, a dichotomy that has not been captured by the classic definition of "pre-transitional" disks., Accepted for publication in A&A Letters; v2 includes language edits as in the A&A publication

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$��$ 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., 18 pages, 13 figures, accepted in A&A

Of the roughly 3000 neutron stars known, only a handful have sub-stellar companions. The most famous of these are the low-mass planets around the millisecond pulsar B1257+12. New evidence indicates that observational biases could still hide a wide variety of planetary systems around most neutron stars. We consider the environment and physical processes relevant to neutron star planets, in particular the effect of X-ray irradiation and the relativistic pulsar wind on the planetary atmosphere. We discuss the survival time of planet atmospheres and the planetary surface conditions around different classes of neutron stars, and define a neutron star habitable zone based on the presence of liquid water and retention of an atmosphere. Depending on as-yet poorly constrained aspects of the pulsar wind, both Super-Earths around B1257+12 could lie within its habitable zone.

The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $\zeta$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process., Comment: Accepted for publication in A&A 18 pages, 12 figures, 15 tables

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 $, 11 pages, 3 figures, accepted for publication in ApJ, correspondence welcome

The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing techniques and their user communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together in a collective effort to understand how disks evolve and disperse at the time of planet formation. This paper was born at the "Protoplanetary Discussions" meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations of disks converge in unveiling disk structure and evolution, and where they diverge and challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions based on limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized by multi-wavelength observations. Specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how the links between them, as well as the apparent contradictions, can help us to disentangle the disk physics and to learn about disk evolution., Comment: Accepted for publication in PASA. 37 pages, 9 figures, 1 table. Revised version: corrected problem in Fig 2

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

Evaporating rocky exoplanets, such as KIC 12557548b, eject large amounts of dust grains, which can trail the planet in a comet-like tail. When such objects occult their host star, the resulting transit signal contains information about the dust in the tail. We aim to use the detailed shape of the Kepler light curve of KIC 12557548b to constrain the size and composition of the dust grains that make up the tail, as well as the mass loss rate of the planet. Using a self-consistent numerical model of the dust dynamics and sublimation, we calculate the shape of the tail by following dust grains from their ejection from the planet to their destruction due to sublimation. From this dust cloud shape, we generate synthetic light curves (incorporating the effects of extinction and angle-dependent scattering), which are then compared with the phase-folded Kepler light curve. We explore the free-parameter space thoroughly using a Markov chain Monte Carlo method. Our physics-based model is capable of reproducing the observed light curve in detail. Good fits are found for initial grain sizes between 0.2 and 5.6 micron and dust mass loss rates of 0.6 to 15.6 M_earth/Gyr (2-sigma ranges). We find that only certain combinations of material parameters yield the correct tail length. These constraints are consistent with dust made of corundum (Al2O3), but do not agree with a range of carbonaceous, silicate, or iron compositions. Using a detailed, physically motivated model, it is possible to constrain the composition of the dust in the tails of evaporating rocky exoplanets. This provides a unique opportunity to probe to interior composition of the smallest known exoplanets., Comment: 18 pages, 11 figures, A&A accepted

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

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $\lambda\,$Bo\"{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems., Comment: Accepted for publication in A&A

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)

Spiral structures are a common feature in scattered-light images of protoplanetary disks, and of great interest as possible tracers of the presence of planets. However, other mechanisms have been put foward to explain them, including self-gravity, disk-envelope interactions, and dead zone boundaries. These mechanisms explain many spirals very well, but are unable to easily account for very loosely wound spirals and single spiral arms. We study the effect of light travel time on the shape of a shadow cast by a clump orbiting close (within ${\sim}1\,$au) of the central star, where there can be significant orbital motion during the light travel time from the clump to the outer disk and then to the sky plane. This delay in light rays reaching the sky plane gives rise to a variety of spiral- and arc-shaped shadows, which we describe with a general fitting formula for a flared, inclined disk., Accepted for publication in A&A Letters. Videos available at dl.dropboxusercontent.com/u/3526708/spiralmovies.zip

The gas temperature structure of protoplanetary disks is a key ingredient for interpreting various disk observations and for quantifying the subsequent evolution of these systems. The comparison of low- and mid-$J$ CO rotational lines is a powerful tool to assess the temperature gradient in the warm molecular layer of disks. Spectrally resolved high-$J$ ($J_{\rm u} > 14$) CO lines probe intermediate distances and heights from the star that are not sampled by (sub-)millimeter CO spectroscopy. This paper presents new {\it Herschel}/HIFI and archival PACS observations of $^{12}$CO, $^{13}$CO and \cii \ emission in 4 Herbig AeBe (HD 100546, HD 97048, IRS 48, HD 163296) and 3 T Tauri (AS 205, S CrA, TW Hya) disks. In the case of the T Tauri systems AS 205 and S CrA, the CO emission has a single-peaked profile, likely due to a slow wind. For all other systems, the {\it Herschel} CO spectra are consistent with pure disk emission and the spectrally-resolved lines (HIFI) and the CO rotational ladder (PACS) are analyzed simultaneously assuming power-law temperature and column density profiles, using the velocity profile to locate the emission in the disk. The temperature profile varies substantially from disk to disk. In particular, $T_{\rm gas}$ in the disk surface layers can differ by up to an order of magnitude among the 4 Herbig AeBe systems with HD 100546 being the hottest and HD 163296 the coldest disk of the sample. Clear evidence of a warm disk layer where $T_{\rm gas} > T_{\rm dust}$ is found in all the Herbig Ae disks. The observed CO fluxes and line profiles are compared to predictions of physical-chemical models. The primary parameters affecting the disk temperature structure are the flaring angle, the gas-to-dust mass ratio the scale height and the dust settling., Accepted for publication in A&A