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.

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.

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.

The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long timescales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded onto an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity ($\alpha\sim0.1$), indicating that gas is lost on timescales $\sim1-10$ Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesise that the observed carbon cavities could be due to radiation pressure or accreting planets.

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$\mu$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.

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.

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.

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.

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

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.

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.

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.

In the last few years, multiwavelength observations have revealed the ubiquity of gaps/rings in circumstellar discs. Here we report the first ALMA observations of HD 92945 at 0.86 mm, that reveal a gap at about 73$\pm$3 au within a broad disc of planetesimals that extends from 50 to 140 au. We find that the gap is $20^{+10}_{-8}$ au wide. If cleared by a planet in situ, this planet must be less massive than 0.6 $M_\mathrm{Jup}$, or even lower if the gap was cleared by a planet that formed early in the protoplanetary disc and prevented planetesimal formation at that radius. By comparing opposite sides of the disc we also find that the disc could be asymmetric. Motivated by the asymmetry and the fact that planets might be more frequent closer to the star in exoplanetary systems, we show that the gap and asymmetry could be produced by two planets interior to the disc through secular resonances. These planets excite the eccentricity of bodies at specific disc locations, opening radial gaps in the planetesimal distribution. New observations are necessary to confirm if the disc is truly asymmetric, thus favouring the secular resonance model, or if the apparent asymmetry is due to a background galaxy, favouring the in-situ planet scenario. Finally, we also report the non-detection of CO and HCN gas confirming that no primordial gas is present. The CO and HCN non-detections are consistent with the destruction of volatile-rich Solar System-like comets.

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.

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

While detecting low mass exoplanets at tens of au is beyond current instrumentation, debris discs provide a unique opportunity to study the outer regions of planetary systems. Here we report new ALMA observations of the 80-200 Myr old Solar analogue HD 107146 that reveal the radial structure of its exo-Kuiper belt at wavelengths of 1.1 and 0.86 mm. We find that the planetesimal disc is broad, extending from 40 to 140 au, and it is characterised by a circular gap extending from 60 to 100 au in which the continuum emission drops by about 50%. We also report the non-detection of the CO J=3-2 emission line, confirming that there is not enough gas to affect the dust distribution. To date, HD 107146 is the only gas-poor system showing multiple rings in the distribution of millimeter sized particles. These rings suggest a similar distribution of the planetesimals producing small dust grains that could be explained invoking the presence of one or more perturbing planets. Because the disk appears axisymmetric, such planets should be on circular orbits. By comparing N-body simulations with the observed visibilities we find that to explain the radial extent and depth of the gap, it would be required the presence of multiple low mass planets or a single planet that migrated through the disc. Interior to HD 107146's exo-Kuiper belt we find extended emission with a peak at ~20 au and consistent with the inner warm belt that was previously predicted based on 22$\mu$m excess as in many other systems. This warm belt is the first to be imaged, although unexpectedly suggesting that it is asymmetric. This could be due to a large belt eccentricity or due to clumpy structure produced by resonant trapping with an additional inner planet.

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.

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}$).

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.

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

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.

HR8799 is a benchmark system for direct imaging studies. It hosts two debris belts, which lie internally and externally to four giant planets. This paper considers how the four known planets and a possible fifth planet, interact with the external population of debris through N-body simulations. We find that when only the known planets are included, the inner edge of the outer belt predicted by our simulations is much closer to the outermost planet than recent ALMA observations suggest. We subsequently include a fifth planet in our simulations with a range of masses and semi-major axes, which is external to the outermost known planet. We find that a fifth planet with a mass and semi-major axis of 0.1$\mathrm{M_J}$ and 138au predicts an outer belt that agrees well with ALMA observations, whilst remaining stable for the lifetime of HR8799 and lying below current direct imaging detection thresholds. We also consider whether inward scattering of material from the outer belt can input a significant amount of mass into the inner belt. We find that for the current age of HR8799, only $\sim$1\% of the mass loss rate of the inner disk can be replenished by inward scattering. However we find that the higher rate of inward scattering during the first $\sim$10Myr of HR8799 would be expected to cause warm dust emission at a level similar to that currently observed, which may provide an explanation for such bright emission in other systems at $\sim10$Myr ages.

The young A0V star HR 4796A is host to a bright and narrow ring of dust, thought to originate in collisions between planetesimals within a belt analogous to the Solar System's Edgeworth-Kuiper belt. Here we present high spatial resolution 880$\mu$m continuum images from the Atacama Large Millimeter Array. The 80au radius dust ring is resolved radially with a characteristic width of 10au, consistent with the narrow profile seen in scattered light. Our modelling consistently finds that the disk is also vertically resolved with a similar extent. However, this extent is less than the beam size, and a disk that is dynamically very cold (i.e. vertically thin) provides a better theoretical explanation for the narrow scattered light profile, so we remain cautious about this conclusion. We do not detect $^{12}$CO J=3-2 emission, concluding that unless the disk is dynamically cold the CO+CO$_2$ ice content of the planetesimals is of order a few percent or less. We consider the range of semi-major axes and masses of an interior planet supposed to cause the ring's eccentricity, finding that such a planet should be more massive than Neptune and orbit beyond 40au. Independent of our ALMA observations, we note a conflict between mid-IR pericenter-glow and scattered light imaging interpretations, concluding that models where the spatial dust density and grain size vary around the ring should be explored.

Low-mass stars build a significant fraction of their total mass during short outbursts of enhanced accretion known as FUor and EXor outbursts. FUor objects are characterized by a sudden brightening of ~5 magnitudes at visible wavelengths within one year and remain bright for decades. EXor objects have lower amplitude outbursts on shorter timescales. Here we discuss a 1.3 mm ALMA mini-survey of eight outbursting sources (three FUor, four EXor, and the borderline object V1647 Ori) in the Orion Molecular Cloud. While previous papers in this series discuss the remarkable molecular outflows observed in the three FUor objects and V1647 Ori, here we focus on the continuum data and the differences and similarities between the FUor and EXor populations. We find that FUor discs are significantly more massive (~80-600 M_JUP) than the EXor objects (~0.5-40 M_JUP ). We also report that the EXor sources lack the prominent outflows seen in the FUor population. Even though our sample is small, the large differences in disc masses and outflow activity suggest that the two types of objects represent different evolutionary stages. The FUor sources seem to be rather compact (Rc < 20-40 au) and to have a smaller characteristic radius for a given disc mass when compared to T Tauri stars. V1118 Ori, the only known close binary systemin our sample, is shown to host a disc around each one of the stellar components. The disc around HBC 494 is asymmetric, hinting at a structure in the outer disc or the presence of a second disc.

A decade of surveys has hinted at a possible higher occurrence rate of debris discs in systems hosting low mass planets. This could be due to common favourable forming conditions for rocky planets close in and planetesimals at large radii. In this paper we present the first resolved millimetre study of the debris disc in the 4.6 Gyr old multiplanet system 61 Vir, combining ALMA and JCMT data at 0.86 mm. We fit the data using a parametric disc model, finding that the disc of planetesimals extends from 30 AU to at least 150 AU, with a surface density distribution of millimetre sized grains with a power law slope of 0.1$^{+1.1}_{-0.8}$. We also present a numerical collisional model that can predict the evolution of the surface density of millimetre grains for a given primordial disc, finding that it does not necessarily have the same radial profile as the total mass surface density (as previous studies suggested for the optical depth), with the former being flatter. Finally, we find that if the planetesimal disc was stirred at 150 AU by an additional unseen planet, that planet should be more massive than 10 M$_{\oplus}$ and lie between 10-20 AU. Lower planet masses and semi-major axes down to 4 AU are possible for eccentricities $\gg$ 0.1.

While most of the known debris discs present cold dust at tens of AU, a few young systems exhibit hot dust analogous to the Zodiacal dust. $\eta$ Corvi is particularly interesting as it is old and it has both, with its hot dust significantly exceeding the maximum luminosity of an in-situ collisional cascade. Previous work suggested that this system could be undergoing an event similar to the Late Heavy Bombardment (LHB) soon after or during a dynamical instability. Here we present ALMA observations of $\eta$ Corvi with a resolution of 1."2 (~22au) to study its outer belt. The continuum emission is consistent with an axisymmetric belt, with a mean radius of 152au and radial FWHM of 46au, which is too narrow compared to models of inward scattering of an LHB-like scenario. Instead, the hot dust could be explained as material passed inwards in a rather stable planetary configuration. We also report a 4sigma detection of CO at ~ 20au. CO could be released in situ from icy planetesimals being passed in when crossing the H$_2$O or CO$_2$ ice lines. Finally, we place constraints on hidden planets in the disc. If a planet is sculpting the disc's inner edge, this should be orbiting at 75-100au, with a mass of 3-30 M$_\oplus$ and an eccentricity < 0.08. Such a planet would be able to clear its chaotic zone on a timescale shorter than the age of the system and scatter material inwards from the outer belt to the inner regions, thus feeding the hot dust.

This paper considers the dynamics of scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, depleted Oort Cloud. We use these outcomes to consider which planetary system architectures maximise the observability of specific signatures, given that signatures should be detected first around systems with optimal architectures (if such systems exist in nature). Giant impact debris is most readily detectable for 0.1-10Mearth planets at 1-5au, depending on detection method and spectral type. While A stars have putative giant impact debris at 4-6au consistent with this sweet spot, that of FGK stars is typically <<1au contrary to expectations; an absence of 1-3au giant impact debris could indicate a low frequency of terrestrial planets there. Three principles maximise cometary influx from exo-Kuiper belts: a chain of closely separated planets interior to the belt, none of which is a Jupiter-like ejector; planet masses not increasing strongly with distance (for a net inward torque on comets); ongoing replenishment of comets, possibly by embedded low-mass planets. A high Oort Cloud comet influx requires no ejector and architectures that maximise the Oort Cloud population. Cold debris disks are usually considered classical Kuiper belt analogues. Here we consider the possibility of detecting scattered disk analogues, which could be betrayed by a broad radial profile and lack of small grains, as well as spherical 100-1000au mini-Oort Clouds. Some implications for escaping planets around young stars, detached planets akin to Sedna, and the formation of super-Earths, are also discussed.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

An increasing number of observations have shown that gaseous debris discs are not an exception. However, until now we only knew of cases around A stars. Here we present the first detection of 12CO (2-1) disc emission around an F star, HD 181327, obtained with ALMA observations at 1.3 mm. The continuum and CO emission are resolved into an axisymmetric disc with ring-like morphology. Using a Markov chain Monte Carlo method coupled with radiative transfer calculations we study the dust and CO mass distribution. We find the dust is distributed in a ring with a radius of 86.0 +- 0.4 AU and a radial width of 23.2 +- 1.0 AU. At this frequency the ring radius is smaller than in the optical, revealing grain size segregation expected due to radiation pressure. We also report on the detection of low level continuum emission beyond the main ring out to ~200 AU. We model the CO emission in the non-LTE regime and we find that the CO is co-located with the dust, with a total CO gas mass ranging between 1.2x10^-6 Mearth and 2.9x10^-6 Mearth, depending on the gas kinetic temperature and collisional partners densities. The CO densities and location suggest a secondary origin, i.e. released from icy planetesimals in the ring. We derive a CO+CO2 cometary composition that is consistent with Solar system comets. Due to the low gas densities it is unlikely that the gas is shaping the dust distribution.

Circumstellar asymmetries such as central warps have recently been shown to cast shadows on outer disks. We investigate the hydrodynamical consequences of such variable illumination on the outer regions of a transition disk, and the development of spiral arms. Using 2D simulations, we follow the evolution of a gaseous disk passively heated by the central star, under the periodic forcing of shadows with an opening angle of $\sim$28$^\circ$. With a lower pressure under the shadows, each crossing results in a variable azimuthal acceleration, which in time develops into spiral density waves. Their pitch angles evolve from $\Pi \sim 15^\circ-22^\circ$ at the onset, to $\sim$11$^\circ$-14$^\circ$, over $\sim$65~AU to 150~AU. Self-gravity enhances the density contrast of the spiral waves, as also reported previously for spirals launched by planets. Our control simulations with unshadowed irradiation do not develop structures, except for a different form of spiral waves seen at later times only in the gravitationally unstable control case. Scattered light predictions in the $H$-band show that such illumination spirals should be observable. We suggest that spiral arms in the case-study transition disk HD~142527 could be explained as a result of shadowing from the tilted inner disk.

The identification of on-going planet formation requires the finest angular resolutions and deepest sensitivities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet-disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. These CPDs have eluded unequivocal detection -their identification requires predictions in CPD tracers. In this work, we aim to assess the observability of embedded CPDs with ALMA as features imprinted in the gas kinematics. We use 3D Smooth Particle Hydrodynamic (SPH) simulations of CPDs around 1 and 5 M_Jup planets at large stellocentric radii, in locally isothermal and adiabatic disks. The simulations are then connected with 3D radiative transfer for predictions in CO isotopologues. Observability is assessed by corrupting with realistic long baseline phase noise extracted from the recent HL Tau ALMA data. We find that the presence of a CPD produces distinct signposts: 1) compact emission separated in velocity from the overall circumstellar disk's Keplerian pattern, 2) a strong impact on the velocity pattern when the Doppler shifted line emission sweeps across the CPD location, and 3) a local increase in the velocity dispersion. We test our predictions with a simulation tailored for HD 100546 -which has a reported protoplanet candidate. We find that the CPDs are detectable in all 3 signposts with ALMA Cycle 3 capabilities for both 1 and 5 M_Jup protoplanets, when embedded in an isothermal disk.

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA difference is an opacity or a real dust segregation.

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution.

The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527.

We present an analytical model for light echoes (LEs) coming from circumstellar material (CSM) around Type Ia Supernovae (SNe Ia). Using this model we find two spectral signatures at 4100 {\AA} and 6200 {\AA} that are useful to identify LEs during the Lira law phase (between 35 and 80 days after maximum light) coming from nearby CSM at distances of 0.01-0.25 pc. We analyze a sample of 89 SNe Ia divided in two groups according to their B-V decline rate during the Lira law phase, and search for LEs from CSM interaction in the group of SNe with steeper slopes by comparing their spectra with our LE model. We find that a model with LEs + pure extinction from interstellar material (ISM) fits better the observed spectra than a pure ISM extinction model that is constant in time, but we find that a decreasing extinction alone explains better the observations without the need of LEs, possibly implying dust sublimation due to the radiation from the SN.

Detailed observations of gaps in protoplanetary disks have revealed structures that drive current research on circumstellar disks. One such feature is the two intensity nulls seen along the outer disk of the HD 142527 system, which are particularly well traced in polarized differential imaging. Here we propose that these are shadows cast by the inner disk. The inner and outer disk are thick, in terms of the unit-opacity surface in H-band, so that the shape and orientation of the shadows inform on the three-dimmensional structure of the system. Radiative transfer predictions on a parametric disk model allow us to conclude that the relative inclination between the inner and outer disks is 70+-5 deg. This finding taps the potential of high-contrast imaging of circumstellar disks, and bears consequences on the gas dynamics of gapped disks, as well as on the physical conditions in the shadowed regions.