Search

The origin and evolution of CO gas in debris disks has been debated since its initial detection. The gas could have a primordial origin, as a remnant of the protoplanetary disk or a secondary exocometary origin. This paper investigates the origin of gas in two debris disks, HD110058 and HD131488, using HST observations of CI and CO, which play critical roles in the gas evolution. We fitted several electronic transitions of CI and CO rovibronic bands to derive column densities and temperatures for each system, revealing high CO column densities ($\sim$3-4 orders of magnitude higher than $\beta$ Pictoris), and low CI/CO ratios in both. Using the exogas model, we simulated the radial evolution of the gas in the debris disk assuming a secondary gas origin. We explored a wide range of CO exocometary release rates and $\alpha$ viscosities, which are the key parameters of the model. Additionally, we incorporated photodissociation due to stellar UV to the exogas model and found that it is negligible for typical CO-rich disks and host stars, even at a few au due to the high radial optical depths in the EUV. We find that the current steady-state secondary release model cannot simultaneously reproduce the CO and CI HST-derived column densities, as it predicts larger CI/CO ratios than observed. Our direct UV measurement of low CI/CO ratios agrees with results derived from recent ALMA findings and may point to vertical layering of CI, additional CI removal, CO shielding processes, or different gas origin scenarios., Comment: 21 pages, 16 figures, accepted to MNRAS

As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those gas, dust and stellar populations form what has come to be known as the ecology of our Galaxy. Because we are deeply embedded in the plane of our Galaxy, it takes up a significant fraction of the sky, with complex dust lanes scattered throughout the optically recognisable bands of the Milky Way. These bands become bright at (sub-)millimetre wavelengths, where we can study dust thermal emission and the chemical and kinematic signatures of the gas. To properly study such large-scale environments, requires deep, large area surveys that are not possible with current facilities. Moreover, where stars form, so too do planetary systems, growing from the dust and gas in circumstellar discs, to planets and planetesimal belts. Understanding the evolution of these belts requires deep imaging capable of studying belts around young stellar objects to Kuiper belt analogues around the nearest stars. Here we present a plan for observing the Galactic Plane and circumstellar environments to quantify the physical structure, the magnetic fields, the dynamics, chemistry, star formation, and planetary system evolution of the galaxy in which we live with AtLAST; a concept for a new, 50m single-dish sub-mm telescope with a large field of view which is the only type of facility that will allow us to observe our Galaxy deeply and widely enough to make a leap forward in our understanding of our local ecology., Comment: 27 pages, submitted to Open Research Europe as part of the AtLAST collection: https://open-research-europe.ec.europa.eu/collections/atlast/about

Mounting observational evidence indicates that cold circumstellar gas is present in debris disk systems. This work focuses on various dynamical processes that debris-disk gas may undergo. We review five mechanisms that can transport angular momentum and their applications to debris disks. These include molecular viscosity, hydrodynamic turbulence, magnetohydrodynamic turbulence, magnetized disk winds, and laminar magnetic stress. We find that molecular viscosity can result in $\alpha$ as high as $\lesssim 0.1$ for sufficiently low densities, while the Rossby wave instability is a possible source of hydrodynamic turbulence and structure formation. We argue that the vertical shear instability is unlikely due to the long cooling times. The onset of the magnetorotational instability (MRI) is dichotomous: for low density disks the MRI can be excited at the midplane, while for high mass disks it may only be operating at $z>2-3H$, if at all. The MHD wind and laminar magnetic stress mechanisms rely on the configuration and strength of any background large-scale magnetic field, the existence of which is uncertain and possibly unlikely. We conclude that the dominant mechanism and its efficiency in transporting angular momentum varies from one system to the other, depending especially closely on the gas density. More detailed analyses shall be performed in the future focusing on representative, nearby debris disks., Comment: 16 pages, 9 figures, submitted to MNRAS and revised

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We obtained high-contrast polarimetric images of HD 181327 in the H band with the SPHERE / IRDIS instrument on the Very Large Telescope (ESO). We complemented them with archival data from HST / NICMOS in the F110W filter reprocessed in the context of the Archival Legacy Investigations of Circumstellar Environments (ALICE) project. We developed a combined forward-modelling framework to simultaneously retrieve the scattering phase function in polarisation and intensity. We detected the debris disc around HD 181327 in polarised light and total intensity. We measured the scattering phase function and the degree of linear polarisation of the dust at 1.6 micron in the birth ring. The maximum polarisation is 23.6% +/- 2.6% and occurs between a scattering angle of 70 deg and 82 deg. We show that compact spherical particles made of a highly refractive and relatively absorbing material in a differential power-law size distribution of exponent $-3.5$ can simultaneously reproduce the polarimetric and total intensity scattering properties of the dust. This type of material cannot be obtained with a mixture of silicates, amorphous carbon, water ice, and porosity, and requires a more refracting component such as iron-bearing minerals. We reveal a striking analogy between the near-infrared polarisation of comets and that of HD 181327. The methodology developed here combining VLT/SPHERE and HST/NICMOS may be applicable in the future to combine the polarimetric capabilities of SPHERE with the sensitivity of JWST., Comment: 21 pages, 22 figures, accepted in Astronomy & Astrophysics

Debris discs are our best means to probe the outer regions of planetary systems. Many studies assume that planets lie at the inner edges of debris discs, akin to Neptune and the Kuiper Belt, and use the disc morphologies to constrain those otherwise-undetectable planets. However, this produces a degeneracy in planet mass and semimajor axis. We investigate the effect of a sculpting planet on the radial surface-density profile at the disc inner edge, and show that this degeneracy can be broken by considering the steepness of the edge profile. Like previous studies, we show that a planet on a circular orbit ejects unstable debris and excites surviving material through mean-motion resonances. For a non-migrating, circular-orbit planet, in the case where collisions are negligible, the steepness of the disc inner edge depends on the planet-to-star mass ratio and the initial-disc excitation level. We provide a simple analytic model to infer planet properties from the steepness of ALMA-resolved disc edges. We also perform a collisional analysis, showing that a purely planet-sculpted disc would be distinguishable from a purely collisional disc and that, whilst collisions flatten planet-sculpted edges, they are unlikely to fully erase a planet's signature. Finally, we apply our results to ALMA-resolved debris discs and show that, whilst many inner edges are too steep to be explained by collisions alone, they are too flat to arise through completed sculpting by non-migrating, circular-orbit planets. We discuss implications of this for the architectures, histories and dynamics in the outer regions of planetary systems., Comment: 24 pages, 18 figures, accepted for publication in MNRAS

We show the first SPHERE/IRDIS and IFS data of the CO-rich debris disc around HD 131488. We use N-body simulations to model both the scattered light images and the SED of the disc in a self-consistent way. We apply the Henyey-Greenstein approximation, Mie theory, and the Discrete Dipole Approximation to model the emission of individual dust grains. Our study shows that only when gas drag is taken into account can we find a model that is consistent with scattered light as well as thermal emission data of the disc. The models suggest a gas surface density of $2\times10^{-5}\,M_\oplus/$au$^2$ which is in agreement with estimates from ALMA observations. Thus, our modelling procedure allows us to roughly constrain the expected amount of gas in a debris disc without actual gas measurements. We also show that the shallow size distribution of the dust leads to a significant contribution of large particles to the overall amount of scattered light. The scattering phase function indicates a dust porosity of $\sim0.2\ldots 0.6$ which is in agreement with a pebble pile scenario for planetesimal growth., Comment: 27 pages, 25 figures, 2 tables, accepted for publication in MNRAS

While most debris disks consist of dust with little or no gas, a fraction has significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semi-analytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplanetary disks that have become depleted of small grains compared to the interstellar medium. Photoelectric heating by the A-stars' FUV radiation is then inefficient, while the stars' EUV and X-ray emissions are weak owing to a lack of surface convective zones capable of driving magnetic activity. In this picture, stars outside the range of spectral types from A through early B are relatively hard to have such long-lived gas disks. Less-massive stars have stronger magnetic activity in the chromosphere, transition region, and corona with resulting EUV and X-ray emission, while more-massive stars have photospheres hot enough to produce strong EUV radiation. In both cases, primordial disk gas is likely to photoevaporate well before 50 Myr. These results come from 0D disk evolution models where we incorporate internal accretion stresses, MHD winds, and photoevaporation by EUV and X-ray photons with luminosities that are functions of the stellar mass and age. A key issue this work leaves open is how some disks become depleted in small dust so that FUV photoevaporation slows. Candidates include grains' growth, settling, radial drift, radiation force, and incorporation into planetary systems., Comment: Submitted to ApJL

Previous surveys in the far-infrared have found very few, if any, M-dwarf debris discs among their samples. It has been questioned whether M-dwarf discs are simply less common than earlier types, or whether the low detection rate derives from the wavelengths and sensitivities available to those studies. The highly sensitive, long wavelength Atacama Large Millimetre/submillimetre Array can shed light on the problem. This paper presents a survey of M-dwarf stars in the young and nearby Beta Pictoris Moving Group with ALMA at Band 7 (880\,$\mu$m). From the observational sample we detect two new sub-mm excesses that likely constitute unresolved debris discs around GJ\,2006\,A and AT\,Mic\,A and model distributions of the disc fractional luminosities and temperatures. From the science sample of 36 M-dwarfs including AU\,Mic we find a disc detection rate of 4/36 or 11.1$^{+7.4}_{-3.3}$\% that rises to 23.1$^{+8.3}_{-5.5}$\% when adjusted for completeness. We conclude that this detection rate is consistent with the detection rate of discs around G and K type stars and that the disc properties are also likely consistent with earlier type stars. We additionally conclude that M-dwarf stars are not less likely to host debris discs, but instead their detection requires longer wavelength and higher sensitivity observations than have previously been employed., Comment: Accepted to MNRAS

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the JWST Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of $0.5\lambda/D$ for an interferometer), which are inaccessible with the classical inner working angles of the JWST coronagraphs. When combined with JWST's unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a $5\sigma$ contrast of $\Delta m{\sim}7.62{\pm}0.13$ mag relative to the host star at separations ${\gtrsim}0.07{"}$, and the contrast deteriorates steeply at separations ${\lesssim}0.07{"}$. However, we detect no additional companions interior to the known companion HIP 65426 b (at separation ${\sim}0.82{"}$ or, $87^{+108}_{-31}\,\rm{au}$). Our observations thus rule out companions more massive than $10{-}12\,\rm{M_{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP 65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on JWST is sensitive to planetary mass companions at close-in separations (${\gtrsim}0.07{"}$), even for thousands of more distant stars at $\sim$100 pc, in addition to the stars in the nearby young moving groups as stated in previous works. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening an essentially unexplored parameter space., Comment: 20 pages, 12 figures, submitted to ApJL

We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim \lambda/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy., Comment: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

Radial velocity surveys suggest that the Solar System may be unusual and that Jupiter-like planets have a frequency <20% around solar-type stars. However, they may be much more common in one of the closest associations in the solar neighbourhood. Young moving stellar groups are the best targets for direct imaging of exoplanets and four massive Jupiter-like planets have been already discovered in the nearby young beta Pic Moving Group (BPMG) via high-contrast imaging, and four others were suggested via high precision astrometry by the European Space Agency's Gaia satellite. Here we analyze 30 stars in BPMG and show that 20 of them might potentially host a Jupiter-like planet as their orbits would be stable. Considering incompleteness in observations, our results suggest that Jupiter-like planets may be more common than previously found. The next Gaia data release will likely confirm our prediction., Comment: 21 pages. Final edited version published on Nature Communications available at URL: https://www.nature.com/articles/s41467-023-41665-0

Warps have often been used to explain disc properties, but well characterised examples are important due to their role in disc evolution. Scattered light images of discs with central gaps have revealed sharp warps, such that the outer rings are shadowed by tilted inner discs. The near-IR intensity drops along the ring around TTauri star DoAr44 have been interpreted in terms of a central warp. We report new ALMA observations of DoAr44 in the continuum at 230 GHz and 350 GHz (at ~10 au), along with a new epoch of SPHERE/IRDIS differential polarised imaging taken during excellent weather conditions. The ALMA observations resolve the ring and confirm the decrements proposed from deconvolution of coarse 336 GHz data. The scattered light image constrains the dips, which correspond to a misaligned inner disc with a relative inclination $\xi$ = 21.4 $^{+6.7}_{-8.3}$ deg. The SPHERE intensity profile shows a morphological change compared to a previous epoch that may be interpreted as a variable orientation of the inner disc, from $\xi$ ~30 deg to $\xi$ ~20 deg. The intensity dips probably correspond to temperature decrements, as their mm-spectral index, $\alpha^{230 GHz}_{350 GHz}$ ~2.0 $\pm$ 0.1, is indicative of optically thick emission. The azimuth of the two temperature decrements are leading clockwise relative to the IR-dips, by $\eta$ = 14.95 deg and $\eta$ = 7.92 deg. For a retrograde disc, such shifts are expected from a thermal lag and imply gas surface densities of $\Sigma_g$ = 117 $\pm$ 10 g/cm$^2$ and $\Sigma_g$ = 48 $\pm$ 10 g/cm$^2$. A lopsided disc, with contrast ratio $f_r$=2.4 $\pm$ 0.5, is also consistent with the large continuum crescent., Comment: accepted in MNRAS

Kuiper-like belts of planetesimals orbiting stars other than the Sun are most commonly detected from the thermal emission of small dust produced in collisions. Emission from gas, most notably CO, highlights the cometary nature of these planetesimals. Here we present models for the release of gas from comet-like bodies in these belts, both due to their thermophysical evolution, most notably the decay of long-lived radioactive nuclides and collisional evolution, including catastrophic and gentler resurfacing collisions. We show that the rate of gas release is not proportional to the rate of dust release, if non-catastrophic collisions or thermal evolution dominate the release of CO gas. In this case, care must be taken when inferring the composition of comets. Non-catastrophic collisions dominate the gas production at earlier times than catastrophic collisions, depending on the properties of the planetesimal belt. We highlight the importance of the thermal evolution of comets, including crucially the decay of long-lived radioactive nuclides, as a source of CO gas around young (<50Myr) planetary systems, if large (10-100s kms) planetesimals are present., Comment: Accepted by MNRAS, 16 pages. Updated following review comments

Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brightness profile of circumstellar discs, can be used to deproject their emission at any inclination as long as they are optically thin and axisymmetric. Furthermore, we extend Frankenstein to account for the vertical thickness of an optically thin disc ($H(r)$) and show how it can be constrained by sampling its posterior probability distribution and assuming a functional form (e.g. constant $h=H/r$), while fitting the radial profile non-parametrically. We use this new method to determine the radial and vertical structure of 16 highly inclined debris discs observed by ALMA. We find a wide range of vertical aspect ratios, $h$, ranging from $0.020\pm0.002$ (AU Mic) to $0.20\pm0.03$ (HD 110058), which are consistent with parametric models. We find a tentative correlation between $h$ and the disc fractional width, as expected if wide discs were more stirred. Assuming discs are self-stirred, the thinnest discs would require the presence of at least 500 km-sized planetesimals. The thickest discs would likely require the presence of planets. We also recover previously inferred and new radial structures, including a potential gap in the radial distribution of HD 61005. Finally, our new extension of Frankenstein also allows constraining how $h$ varies as a function of radius, which we test on 49 Ceti, finding that $h$ is consistent with being constant., Comment: Accepted for publication in MNRAS. 17 pages. 16 figures

The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA observations of CO and CI emission, complemented by CII data from Herschel, for a sample of 14 debris disks. This expands the number of disks with ALMA measurements of both CO and CI by ten disks. We present new detections of CI emission towards three disks: HD 21997, HD 121191 and HD 121617. We use a simple disk model to derive gas masses and column densities. We find that current state-of-the-art models of secondary gas production overpredict the neutral carbon content of debris disk gas. This does not rule out a secondary origin, but might indicate that the models require an additional C removal process. Alternatively, the gas might be produced in transient events rather than a steady-state collisional cascade. We also test a primordial gas origin by comparing our results to a simplified thermo-chemical model. This yields promising results, but more detailed work is required before a conclusion can be reached. Our work demonstrates that the combination of C and CO data is a powerful tool to advance our understanding of debris disk gas., Comment: 90 pages, 60 figures. Accepted for publication in ApJ. version 2: additional acknowledgement. versions 3, 4: minor edits

The radial structure of debris discs can encode important information about their dynamical and collisional history. In this paper we present a 3-phase analytical model to analyse the collisional evolution of solids in debris discs, focusing on their joint radial and temporal dependence. Consistent with previous models, we find that as the largest planetesimals reach collisional equilibrium in the inner regions, the surface density of dust and solids becomes proportional to $\sim r^{2}$ within a certain critical radius. We present simple equations to estimate the critical radius and surface density of dust as a function of the maximum planetesimal size and initial surface density in solids (and vice versa). We apply this model to ALMA observations of 7 wide debris discs. We use both parametric and non-parametric modelling to test if their inner edges are shallow and consistent with collisional evolution. We find that 4 out of 7 have inner edges consistent with collisional evolution. Three of these would require small maximum planetesimal sizes below 10 km, with HR 8799's disc potentially lacking solids larger than a few centimeters. The remaining systems have inner edges that are much sharper, which requires maximum planetesimal sizes $\gtrsim10$ km. Their sharp inner edges suggest they could have been truncated by planets, which JWST could detect. In the context of our model, we find that the 7 discs require surface densities below a Minimum Mass Solar Nebula, avoiding the so-called disc mass problem. Finally, during the modelling of HD 107146 we discover that its wide gap is split into two narrower ones, which could be due to two low-mass planets formed within the disc., Comment: Accepted for publication in MNRAS, 21 pages, 11 figures

We present Atacama Large Millimeter Array (ALMA) observations of the young, gas-rich debris disk around HD110058 at 0.3-0.6\arcsec resolution. The disk is detected in the 0.85 and 1.3~mm continuum, as well as in the J=2-1 and J=3-2 transitions of $^{12}$CO and $^{13}$CO. The observations resolve the dust and gas distributions and reveal that this is the smallest debris disk around stars of similar luminosity observed by ALMA. The new ALMA data confirm the disk is very close to edge-on, as shown previously in scattered light images. We use radiative transfer modeling to constrain the physical properties of dust and gas disks. The dust density peaks at around 31~au and has a smooth outer edge that extends out to $\sim70$~au. Interestingly, the dust emission is marginally resolved along the minor axis, which indicates that it is vertically thick if truly close to edge-on with an aspect ratio between 0.13 and 0.28. We also find that the CO gas distribution is more compact than the dust \ah{(similarly to the disk around 49 Ceti)}, which could be due to a low viscosity and a higher gas release rate at small radii. Using simulations of the gas evolution taking into account the CO photodissociation, shielding, and viscous evolution, we find that HD~110058's CO gas mass and distribution are consistent with a secondary origin scenario. Finally, we find that the gas densities may be high enough to cause the outward drift of small dust grains in the disk.

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $\mu$m to 20 $\mu$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion., Comment: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $\mu$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $\mu$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $\mu$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$\sigma$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $\mu$m and MIRI at 11.3 $\mu$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $\mu$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail., Comment: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$\mu$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$\mu$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles., Comment: 21 pages, 7 figures. Accepted for Publication in PASP

We present new ALMA Band 7 observations of the edge-on debris disc around the M1V star GSC 07396-00759. At ~20 Myr old and in the beta Pictoris Moving Group along with AU Mic, GSC 07396-00759 joins it in the handful of low mass M-dwarf discs to be resolved in the sub-mm. With previous VLT/SPHERE scattered light observations we present a multi-wavelength view of the dust distribution within the system under the effects of stellar wind forces. We find the mm dust grains to be well described by a Gaussian torus at 70 au with a FWHM of 48 au and we do not detect the presence of CO in the system. Our ALMA model radius is significantly smaller than the radius derived from polarimetric scattered light observations, implying complex behaviour in the scattering phase function. The brightness asymmetry in the disc observed in scattered light is not recovered in the ALMA observations, implying that the physical mechanism only affects smaller grain sizes. High resolution follow-up observations of the system would allow investigation into its unique dust features as well as provide a true coeval comparison for its smaller sibling AU Mic, singularly well observed amongst M-dwarfs systems., Comment: 14 pages, 11 figures, 1 table, Accepted for publication in MNRAS

This review chapter for young researchers presents our current understanding of debris discs. It introduces some of their basic properties and observables, and describes how we think they form and collisionally evolve. Special emphasis is dedicated to ALMA observations of the dust and gas, which constrains the distribution of planetesimals, their volatile composition, and potential volatile delivery to planetary atmospheres., Comment: Review chapter of the book "Planetary Systems Now", to be published by World Scientific

The nearby V4046 Sgr spectroscopic binary hosts a gas-rich disc known for its wide cavity and dusty ring. We present high resolution ($\sim$20 mas or 1.4 au) ALMA observations of the 1.3mm continuum of V4046 Sgr which, combined with SPHERE--IRDIS polarised images and a well-sampled spectral energy distribution (SED), allow us to propose a physical model using radiative transfer (RT) predictions. The ALMA data reveal a thin ring at a radius of 13.15$\pm$0.42 au (Ring13), with a radial width of 2.46$\pm$0.56 au. Ring13 is surrounded by a $\sim$10 au-wide gap, and it is flanked by a mm-bright outer ring (Ring24) with a sharp inner edge at 24 au. Between 25 and $\sim$35 au the brightness of Ring24 is relatively flat and then breaks into a steep tail that reaches out to $\sim$60 au. In addition, central emission is detected close to the star which we interpret as a tight circumbinary ring made of dust grains with a lower size limit of 0.8 mm at 1.1 au. In order to reproduce the SED, the model also requires an inner ring at $\sim$5 au (Ring5) composed mainly of small dust grains, hiding under the IRDIS coronagraph, and surrounding the inner circumbinary disc. The surprisingly thin Ring13 is nonetheless roughly 10 times wider than its expected vertical extent. The strong near-far disc asymmetry at 1.65 $\mu$m points at a very forward-scattering phase function and requires grain radii of no less than 0.4 $\mu$m., Comment: 11 pages, 7 figures

Models of terrestrial planet formation predict that the final stages of planetary assembly, lasting tens of millions of years beyond the dispersal of young protoplanetary disks, are dominated by planetary collisions. It is through these giant impacts that planets like the young Earth grow to their final mass and achieve long-term stable orbital configurations. A key prediction is that these impacts produce debris. To date, the most compelling observational evidence for post-impact debris comes from the planetary system around the nearby 23 Myr-old A star HD 172555. This system shows large amounts of fine dust with an unusually steep size distribution and atypical dust composition, previously attributed to either a hypervelocity impact or a massive asteroid belt. Here, we report the spectrally resolved detection of a CO gas ring co-orbiting with dusty debris between ~6-9 au - a region analogous to the outer terrestrial planet region of our Solar System. Taken together, the dust and CO detections favor a giant impact between large, volatile-rich bodies. This suggests that planetary-scale collisions, analogous to the Moon-forming impact, can release large amounts of gas as well as debris, and that this gas is observable, providing a window into the composition of young planets., Comment: 21 pages (double spaced, including references and figures), 3 figures, 1 table, published in Nature 20/10/2021, published version available online through SharedIt initiative at https://rdcu.be/cAlEy

Context. Dynamical and turbulent motions of gas in a protoplanetary disk are crucial for their evolution and affect planet formation. Recent observations suggest weak turbulence in the disk's outer regions. However, the physical mechanism of turbulence in these outer regions remains uncertain. The vertical shear instability (VSI) is a promising mechanism to produce turbulence in disks. Aims. Our aim is to study the observability of the gas velocity structure produced by the VSI via CO kinematics with ALMA. Methods. We perform global 3D hydrodynamical simulations of a VSI-unstable disk. We post-process the simulation results with radiative transfer calculations, and produce synthetic predictions of CO rotational emission lines. Following, we compute the line of sight velocity map, and its deviations from a sub-Keplerian equilibrium solution. We explore the detectability of the VSI by identifying kinematic signatures using realistic simulated observations. Results. Our 3D simulations of the VSI show the steady state dynamics of the gas in great detail. From the velocity structure we infer a turbulent stress value of $\alpha_{r\phi}=1.4 \times 10^{-4}$. On large scales, we observe velocity deviations of 50 m s$^{-1}$ as axisymmetric rings. We find optimal conditions at $i \lesssim 20^{\circ}$ to trace for the kinematic structures of the VSI. We found that current diagnostics to constrain gas turbulence from non-thermal broadening of the line emission are not applicable to anisotropic VSI turbulence. Conclusions. The detection of kinematic signatures produced by the VSI is possible with ALMA. Observations including an extended antenna configuration combined with the highest spectral resolution available are needed for robust detection. The characterization of the large-scale velocity perturbations is required to constrain the turbulence level produced by the VSI from gas observations., Comment: Accepted for publication in A&A. 17 pages with 10 figures and 1 table

Self-assembly is one of the crucial mechanisms allowing to design multifunctional materials. Soft hybrid materials contain components of different nature and exhibit competitive interactions which drive self-organisation into structures of a particular function. Here we demonstrate a novel type of a magnetic hybrid material where the molecular tilt can be manipulated through a delicate balance between the topologically-assisted colloidal self-assembly of \rev{magnetic nanoparticles} and the anisotropic molecular interactions in a liquid crystal matrix.

The dust production in debris discs by grinding collisions of planetesimals requires their orbits to be stirred. However, stirring levels remain largely unconstrained, and consequently the stirring mechanisms as well. This work shows how the sharpness of the outer edge of discs can be used to constrain the stirring levels. Namely, the sharper the edge is the lower the eccentricity dispersion must be. For a Rayleigh distribution of eccentricities ($e$), I find that the disc surface density near the outer edge can be parametrised as $\tanh[(r_{\max}-r)/l_{\rm out}]$, where $r_{\max}$ approximates the maximum semi-major axis and $l_{\rm out}$ defines the edge smoothness. If the semi-major axis distribution has sharp edges $e_\mathrm{rms}$ is roughly $1.2 l_{\rm out}/r_{\max}$, or $e_\mathrm{rms}=0.77 l_{\rm out}/r_{\max}$ if semi-major axes have diffused due to self-stirring. This model is fitted to ALMA data of five wide discs: HD 107146, HD 92945, HD 206893, AU Mic and HR 8799. The results show that HD 107146, HD 92945 and AU Mic have the sharpest outer edges, corresponding to $e_\mathrm{rms}$ values of $0.121\pm0.05$, $0.15^{+0.07}_{-0.05}$ and $0.10\pm0.02$ if their discs are self-stirred, suggesting the presence of Pluto-sized objects embedded in the disc. Although these stirring values are larger than typically assumed, the radial stirring of HD 92945 is in good agreement with its vertical stirring constrained by the disc height. HD 206893 and HR~8799, on the other hand, have smooth outer edges that are indicative of scattered discs since both systems have massive inner companions., Comment: Accepted for publication in MNRAS. 15 ppages, 11 figures

The exoplanetary system of HR 8799 is one of the rare systems in which multiple planets have been directly imaged. Its architecture is strikingly similar to that of the Solar System, with the four imaged giant planets surrounding a warm dust belt analogous to the Asteroid Belt, and themselves being surrounded by a cold dust belt analogue to the Kuiper Belt. Previous observations of this cold belt with ALMA in Band 6 (1.3 mm) revealed its inner edge, but analyses of the data differ on its precise location. It was therefore unclear whether the outermost planet HR 8799 b was dynamically sculpting it or not. We present here new ALMA observations of this debris disk in Band 7 (340 GHz, 880 micron). These are the most detailed observations of this disk obtained so far, with a resolution of 1" (40 au) and sensitivity of $9.8\,\mu\mathrm{Jy\,beam^{-1}}$, which allowed us to recover the disk structure with high confidence. In order to constrain the disk morphology, we fit its emission using radiative transfer models combined with a MCMC procedure. We find that this disk cannot be adequately represented by a single power law with sharp edges. It exhibits a smoothly rising inner edge and smoothly falling outer edge, with a peak in between, as expected from a disk that contains a high eccentricity component, hence confirming previous findings. Whether this excited population and inner edge shape stem from the presence of an additional planet remains, however, an open question., Comment: 30 pages, 12 figures, 7 tables. Accepted in AJ

Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and wavelength of 1.3 mm with the Atacama Large Millimeter/submillimeter Array (ALMA). We observe a broad disk extending from a radius of <51 au to 194^{+13}_{-2} au. We model the disk with a continuous, gapped, and double power-law model of the surface density profile, and find strong evidence for a local minimum in the surface density distribution near a radius of 70 au, consistent with a gap in the disk with an inner radius of 63^{+8}_{-16} au and width 31^{+11}_{-7} au. Gapped structure has been observed in four other debris disks -- essentially every other radially resolved debris disk observed with sufficient angular resolution and sensitivity with ALMA -- and could be suggestive of the presence of an additional planetary-mass companion., Comment: 24 pages, 12 figures, Accepted for publication in ApJ

We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in more evolved sources with Class II SEDs. Our sample includes 5 objects with previously known large dust cavities (r $>$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques., Comment: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

We present ALMA observations of the $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1 transitions and the 1.3\,mm continuum emission for the circumbinary disc around HD 142527, at an angular resolution of $\approx$\,0\farcs3. We observe multiple spiral structures in intensity, velocity and velocity dispersion for the $^{12}$CO and $^{13}$CO gas tracers. A newly detected $^{12}$CO spiral originates from the dust horseshoe, and is rotating at super-Keplerian velocity or vertically ascending, whilst the inter-spiral gas is rotating at sub-Keplerian velocities. This new spiral possibly connects to a previously identified spiral, thus spanning > 360$^\circ$. A spatial offset of ~30 au is observed between the $^{12}$CO and $^{13}$CO spirals, to which we hypothesize that the gas layers are propagating at different speeds (``surfing'') due to a non-zero vertical temperature gradient. Leveraging the varying optical depths between the CO isotopologues, we reconstruct temperature and column density maps of the outer disc. Gas surface density peaks at r\,$\approx$\,180\,au, coincident with the peak of continuum emission. Here the dust grains have a Stokes number of $\approx$\,1, confirming radial and azimuthal trapping in the horseshoe. We measure a cavity radius at half-maximum surface density of $\approx$\,100\,au, and a cavity eccentricity between 0.3 and 0.45.

HD 38206 is an A0V star in the Columba association, hosting a debris disc first discovered by IRAS. Further observations by Spitzer and Herschel showed that the disc has two components, likely analogous to the asteroid and Kuiper belts of the Solar System. The young age of this star makes it a prime target for direct imaging planet searches. Possible planets in the system can be constrained using the debris disc. Here we present the first ALMA observations of the system's Kuiper belt and fit them using a forward modelling MCMC approach. We detect an extended disc of dust peaking at around 180 au with a width of 140 au. The disc is close to edge on and shows tentative signs of an asymmetry best fit by an eccentricity of $0.25^{+0.10}_{-0.09}$. We use the fitted parameters to determine limits on the masses of planets interior to the cold belt. We determine that a minimum of four planets are required, each with a minimum mass of 0.64 M$_J$, in order to clear the gap between the asteroid and Kuiper belts of the system. If we make the assumption that the outermost planet is responsible for the stirring of the disc, the location of its inner edge and the eccentricity of the disc, then we can more tightly predict its eccentricity, mass and semimajor axis to be $e_{\rm{p}}=0.34^{+0.20}_{-0.13}$, $m_{\rm{p}}=0.7^{+0.5}_{-0.3}\,\rm{M}_{\rm{J}}$ and $a_{\rm{p}}=76^{+12}_{-13}\,\rm{au}$., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a steep inner edge (at 60-80 au). The belt is very broad, extending out to at least 150-200 au. We strongly detect the star and set a stringent upper limit to warm dust emission previously detected in the infrared. We discuss three scenarios that could explain the architecture of Vega's planetary system, including the new {ALMA} constraints: no outer planets, a chain of low-mass planets, and a single giant planet. The planet-less scenario is only feasible if the outer belt was born with the observed sharp inner edge. If instead the inner edge is currently being truncated by a planet, then the planet must be $\gtrsim$6 M$_{\oplus}$ and at $\lesssim71$ au to have cleared its chaotic zone within the system age. In the planet chain scenario, outward planet migration and inward scattering of planetesimals could produce the hot and warm dust observed in the inner regions of the system. In the single giant planet scenario, an asteroid belt could be responsible for the warm dust, and mean motion resonances with the planet could put asteroids on star-grazing orbits, producing the hot dust., Comment: 18 pages, 3 figures, Accepted for publication in ApJ

Gas detection around main sequence stars is becoming more common with around 20 systems showing the presence of CO. However, more detections are needed, especially around later spectral type stars to better understand the origin of this gas and refine our models. To do so, we carried out a survey of 10 stars with predicted high likelihoods of secondary CO detection using ALMA in band 6. We looked for continuum emission of mm-dust as well as gas emission (CO and CN transitions). The continuum emission was detected in 9/10 systems for which we derived the discs' dust masses and geometrical properties, providing the first mm-wave detection of the disc around HD 106906, the first mm-wave radius for HD 114082, 117214, HD 15745, HD 191089 and the first radius at all for HD 121191. A crucial finding of our paper is that we detect CO for the first time around the young 10-16 Myr old G1V star HD 129590, similar to our early Sun. The gas seems colocated with its planetesimal belt and its total mass is likely between $2-10 \times 10^{-5}$ M$_\oplus$. This first gas detection around a G-type main-sequence star raises questions as to whether gas may have been released in the Solar System as well in its youth, which could potentially have affected planet formation. We also detected CO gas around HD 121191 at a higher S/N than previously and find that the CO lies much closer-in than the planetesimals in the system, which could be evidence for the previously suspected CO viscous spreading owing to shielding preventing its photodissociation. Finally, we make estimates for the CO content in planetesimals and the HCN/CO outgassing rate (from CN upper limits), which we find are below the level seen in Solar System comets in some systems., Comment: Accepted for publication in MNRAS. 22 pages, 13 figures

Demographic surveys of protoplanetary disks, carried out mainly with ALMA, have provided access to a large range of disk dust masses ($M_{\rm{dust}}$) around stars with different stellar types and in different star-forming regions. These surveys found a power-law relation between $M_{\rm{dust}}$ and $M_{\star}$ that steepens in time, but which is also flatter for transition disks (TDs). We performed dust evolution models, which included perturbations to the gas surface density with different amplitudes to investigate the effect of particle trapping on the $M_{\rm{dust}}-M_{\star}$ relation. These perturbations were aimed at mimicking pressure bumps that originated from planets. We focused on the effect caused by different stellar and disk masses based on exoplanet statistics that demonstrate a dependence of planet mass on stellar mass and metallicity. Models of dust evolution can reproduce the observed $M_{\rm{dust}}-M_{\star}$ relation in different star-forming regions when strong pressure bumps are included and when the disk mass scales with stellar mass (case of $M_{\rm{disk}}=0.05\,M_\star$ in our models). This result arises from dust trapping and dust growth beyond centimeter-sized grains inside pressure bumps. However, the flatter relation of $M_{\rm{dust}}-M_{\star}$ for TDs and disks with substructures cannot be reproduced by the models unless the formation of boulders is inhibited inside pressure bumps. In the context of pressure bumps originating from planets, our results agree with current exoplanet statistics on giant planet occurrence increasing with stellar mass, but we cannot draw a conclusion about the type of planets needed in the case of low-mass stars. This is attributed to the fact that for $M_\star<1\,M_\odot$, the observed $M_{\rm{dust}}$ obtained from models is very low due to the efficient growth of dust particles beyond centimeter-sizes inside pressure bumps., Comment: Accepted for publication in A&A, minor changes after language edition

The Solar System gas giants are each surrounded by many moons, with at least 50 prograde satellites thought to have formed from circumplanetary material. Just like the Sun is not the only star surrounded by planets, extrasolar gas giants are likely surrounded by satellite systems. Here, we report on ALMA observations of four <40 Myr old stars with directly-imaged companions: PZ Tel, AB Pic, 51 Eri, and $\kappa$ And. Continuum emission at 1.3 mm is undetected for any of the systems. Since these are directly-imaged companions, there is knowledge of their temperatures, masses and locations. These allow for upper limits on the amount of circumplanetary dust to be derived from detailed radiative transfer models. These protolunar disc models consider two disc sizes: 0.4 and 0.04 times the exoplanet's Hill radius. The former is representative of hydrodynamic simulations of circumplanetary discs while the latter a case with significant radial drift of solids. The more compact case is also motivated by the semi-major axis of Callisto, enclosing Jupiter's Galilean satellites. All upper limits fall below the expected amount of dust required to explain regular satellite systems ($10^{-4}$ times the mass of their central planet). Upper limits are compared with viscous evolution and debris disc models. Our analysis suggests that the non detections can be interpreted as evidence of dust growth beyond metre sizes to form moonetesimals in timescales <10 Myr. This sample increases by 50% the number of ALMA non-detections of young companions available in the literature., Comment: 6 pages, 3 figures, accepted for publication in MNRAS

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts., Comment: 14 pages, 6 figures, accepted to ApJ

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening., Comment: 10 pages, 6 figures. Accepted for publication in ApJ Letters

The mass of the gaseous reservoir in young circumstellar disks is a crucial initial condition for the formation of planetary systems, but estimates vary by orders of magnitude. In some disks with resolvable cavities, sharp inner disk warps cast two-sided shadows on the outer rings; can the cooling of the gas as it crosses the shadows bring constraints on its mass? The finite cooling timescale should result in dust temperature decrements shifted ahead of the optical/IR shadows in the direction of rotation. However, some systems show temperature drops, while others do not. The depth of the drops and the amplitude of the shift depend on the outer disk surface density Sigma through the extent of cooling during the shadow crossing time, and also on the efficiency of radiative diffusion. These phenomena may bear observational counterparts, which we describe with a simple one-dimensional model. An application to the HD142527 disk suggests an asymmetry in its shadows, and predicts a >~10deg shift for a massive gaseous disk, with peak Sigma > 8.3 g/cm2. Another application to the DoAr44 disk limits the peak surface density to Sigma < 13g/cm2, Comment: accepted to MNRAS Letters

Evidence for exocomets, icy bodies in extrasolar planetary systems, has rapidly increased over the past decade. Volatiles are detected through the gas that exocomets release as they collide and grind down within their natal belts, or as they sublimate once scattered inwards to the regions closest to their host star. Most detections are in young, 10 to a few 100 Myr-old systems that are undergoing the final stages of terrestrial planet formation. This opens the exciting possibility to study exocomets at the epoch of volatile delivery to the inner regions of planetary systems. Detection of molecular and atomic gas in exocometary belts allows us to estimate molecular ice abundances and overall elemental abundances, enabling comparison with the Solar Nebula and Solar System comets. At the same time, observing star-grazing exocomets transiting in front of their star (for planetary systems viewed edge-on) and exozodiacal dust in the systems' innermost regions gives unique dynamical insights into the inward scattering process producing delivery to inner rocky planets. The rapid advances of this budding subfield of exoplanetary science will continue in the short term with the upcoming JWST, WFIRST and PLATO missions. In the longer term, the priority should be to explore the full composition of exocomets, including species crucial for delivery and later prebiotic synthesis. Doing so around an increasingly large population of exoplanetary systems is equally important, to enable comparative studies of young exocomets at the epoch of volatile delivery. We identify the proposed LUVOIR and Origins flagship missions as the most promising for a large-scale exploration of exocometary gas, a crucial component of the chemical heritage of young exo-Earths., Comment: White Paper submitted to the US National Academy of Sciences Astro2020 Decadal Survey (8 pages, 3 figures)

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario., Comment: 18 pages, 8 figures, accepted for publication in MNRAS

The inclination distribution of the Kuiper belt provides unique constraints on its origin and dynamical evolution, motivating vertically resolved observations of extrasolar planetesimal belts. We present ALMA observations of millimeter emission in the near edge-on planetesimal belt around $\beta$ Pictoris, finding that the vertical distribution is significantly better described by the sum of two Gaussians compared to a single Gaussian. This indicates that, as for the Kuiper belt, the inclination distribution of $\beta$ Pic's belt is better described by the sum of dynamically hot and cold populations rather than a single component. The hot and cold populations have RMS inclinations of 8.9$^{+0.7}_{-0.5}$ and 1.1$^{+0.5}_{-0.5}$ degrees. We also report that an axisymmetric belt model provides a good fit to new and archival ALMA visibilities, and confirm that the midplane is misaligned with respect to $\beta$ Pic b's orbital plane. However, we find no significant evidence for either the inner disk tilt observed in scattered light and CO emission or the South-West/North-East (SW/NE) asymmetry previously reported for millimeter emission. Finally, we consider the origin of the belt's inclination distribution. Secular perturbations from $\beta$ Pic b are unlikely to provide sufficient dynamical heating to explain the hot population throughout the belt's radial extent, and viscous stirring from large bodies within the belt alone cannot reproduce the two populations observed. This argues for an alternative or additional scenario, such as planetesimals being born with high inclinations, or the presence of a `$\beta$ Pic c' planet, potentially migrating outwards near the belt's inner edge., Comment: 19 pages, 12 figures; Accepted for publication in AJ

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

Exocomets scattered by planets have been invoked to explain observations in multiple contexts, including the frequently found near- and mid-infrared excess around nearby stars arising from exozodiacal dust. Here we investigate how the process of inward scattering of comets originating in an outer belt, is affected by the architecture of a planetary system, to determine whether this could lead to observable exozodi levels or deliver volatiles to inner planets. Using N-body simulations, we model systems with different planet mass and orbital spacing distributions in the 1-50 AU region. We find that tightly packed ($\Delta a_\mathrm{p}<20 R_\mathrm{Hill,m}$) low mass planets are the most efficient at delivering material to exozodi regions (5-7% of scattered exocomets end up within 0.5 AU at some point), although the exozodi levels do not vary by more than a factor of ~7 for the architectures studied here. We suggest that emission from scattered dusty material in between the planets could provide a potential test for this delivery mechanism. We show that the surface density of scattered material can vary by two orders of magnitude (being highest for systems of low mass planets with medium spacing), whilst the exozodi delivery rate stays roughly constant, and that future instruments such as JWST could detect it. In fact for $\eta$ Corvi, the current Herschel upper limit rules out the scattering scenario by a chain of $\lesssim$30 M$_\oplus$ planets. Finally, we show that exocomets could be efficient at delivering cometary material to inner planets (0.1-1% of scattered comets are accreted per inner planet). Overall, the best systems at delivering comets to inner planets are the ones that have low mass outer planets and medium spacing ($\sim20 R_\mathrm{Hill,m}$)., Comment: Accepted for publication in MNRAS. 23 pages, 11 figures

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains., Comment: accepted to MNRAS

Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring in the submm radio continuum, DoAr44 is likely also a warped system. We constrain the warp geometry by comparing radiative transfer predictions with the DPI data in H band (Q_\phi(H)) and with a re-processing of archival 336GHz ALMA observations. The observed DPI shadows have coincident radio counterparts, but the intensity drops are much deeper in Q_\phi(H) (~88%), compared to the shallow drops at 336GHz (~24%). Radiative transfer predictions with an inner disk tilt of ~30+-5deg approximately account for the observations. ALMA long-baseline observations should allow the observation of the warped gas kinematics inside the cavity of DoAr44., Comment: 11 pages, 10 figures, 3 tables, accepted for publication in MNRAS

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions., Comment: 16 pages, 12 figures, accepted to MNRAS. Movies at http://users.monash.edu.au/~dprice/pubs/HD142527/

The TRAPPIST-1 system is unique in that it has a chain of seven terrestrial Earth-like planets located close to or in its habitable zone. In this paper, we study the effect of potential cometary impacts on the TRAPPIST-1 planets and how they would affect the primordial atmospheres of these planets. We consider both atmospheric mass loss and volatile delivery with a view to assessing whether any sort of life has a chance to develop. We ran N-body simulations to investigate the orbital evolution of potential impacting comets, to determine which planets are more likely to be impacted and the distributions of impact velocities. We consider three scenarios that could potentially throw comets into the inner region (i.e within 0.1au where the seven planets are located) from an (as yet undetected) outer belt similar to the Kuiper belt or an Oort cloud: Planet scattering, the Kozai-Lidov mechanism and Galactic tides. For the different scenarios, we quantify, for each planet, how much atmospheric mass is lost and what mass of volatiles can be delivered over the age of the system depending on the mass scattered out of the outer belt. We find that the resulting high velocity impacts can easily destroy the primordial atmospheres of all seven planets, even if the mass scattered from the outer belt is as low as that of the Kuiper belt. However, we find that the atmospheres of the outermost planets f, g and h can also easily be replenished with cometary volatiles (e.g. $\sim$ an Earth ocean mass of water could be delivered). These scenarios would thus imply that the atmospheres of these outermost planets could be more massive than those of the innermost planets, and have volatiles-enriched composition., Comment: accepted to MNRAS. We now quantify the amount of water lost owing to impacts