Search

Mounting observational evidence indicates that cold circumstellar gas is present in debris disc systems. This work focuses on various dynamical processes that debris disc gas may undergo. We review five mechanisms that can transport angular momentum and their applications to debris discs. These include molecular viscosity, hydrodynamic turbulence, magnetohydrodynamic (MHD) turbulence, magnetized disc winds, and laminar magnetic stress. We find that molecular viscosity can result in α as high as ≲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 discs the MRI can be excited at the mid-plane, while for high-mass discs it may only be operating at z > 2 − 3 H , 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 vary 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 discs. [ABSTRACT FROM AUTHOR]

Extreme debris discs (EDDs) are bright and warm circumstellar dusty structures around main sequence stars. They may represent the outcome of giant collisions occuring in the terrestrial region between large planetesimals or planetary bodies, and thus provide a rare opportunity to peer into the aftermaths of these events. Here, we report on results of a mini-survey we conducted with the aim to increase the number of known EDDs, investigate the presence of solid-state features around 10  μ m in eight EDDs, and classify them into the silica or silicate dominated groups. We identify four new EDDs and derive their fundamental properties. For these, and for four other previously known discs, we study the spectral energy distribution around 10  μ m by means of VLT/VISIR photometry in three narrow-band filters and conclude that all eight objects likely exhibit solid-state emission features from sub-micron grains. We find that four discs probably belong to the silicate dominated subgroup. Considering the age distribution of the entire EDD sample, we find that their incidence begins to decrease only after 300 Myr, suggesting that the earlier common picture that these objects are related to the formation of rocky planets may not be exclusive, and that other processes may be involved for older objects (≳100 Myr). Because most of the older EDD systems have wide, eccentric companions, we suggest that binarity may play a role in triggering late giant collisions. [ABSTRACT FROM AUTHOR]

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 spectral energy distribution 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 × 10−5 M⊕ 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 ∼0.2...0.6 which is in agreement with a pebble pile scenario for planetesimal growth. [ABSTRACT FROM AUTHOR]

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 (ALMA) 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 μ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}$|  per cent that rises to 23.1 |$^{+8.3}_{-5.5}$|  per cent 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. [ABSTRACT FROM AUTHOR]

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 (<50 Myr) planetary systems, if large (10–100 km) planetesimals are present. [ABSTRACT FROM AUTHOR]

Warps have often been used to explain disc properties, but well-characterized 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-infrared intensity drops along the ring around T-Tauri star DoAr 44 have been interpreted in terms of a central warp. We report new ALMA observations of DoAr 44 in the continuum at 230 and 350 GHz (at ∼10 au), along with a new epoch of Spectro Polarimetric High-contrast Exoplanet REsearch (SPHERE)/Infrared Dual-band Imager and Spectrograph differential polarized 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 ξ = 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 ξ ∼ 30 deg to ξ ∼ 20 deg. The intensity dips probably correspond to temperature decrements, as their mm-spectral index, |$\alpha ^{230 \textrm {GHz}}_{350 \textrm {GHz}} \sim$|  2.0 ± 0.1, is indicative of optically thick emission. The azimuth of the two temperature decrements are leading clockwise relative to the infrared-dips, by η = 14.95 deg and η = 7.92 deg. For a retrograde disc, such shifts are expected from a thermal lag and imply gas surface densities of Σg = 117 ± 10 g cm−2 and Σg = 48 ± 10 g cm−2. A lopsided disc, with contrast ratio f r = 2.4 ± 0.5, is also consistent with the large continuum crescent. [ABSTRACT FROM AUTHOR]

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 structures of 16 highly inclined debris discs observed by ALMA. We find a wide range of vertical aspect ratios, h , ranging from 0.020 ± 0.002 (AU Mic) to 0.20 ± 0.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. [ABSTRACT FROM AUTHOR]

The radial structure of debris discs can encode important information about their dynamical and collisional history. In this paper, we present a three-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 ∼ 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 Atacama Large Millimeter/submillimeter Array observations of seven 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 four out of seven 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 centimetres. The remaining systems have inner edges that are much sharper, which requires maximum planetesimal sizes ≳ 10 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 seven 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. [ABSTRACT FROM AUTHOR]

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 β 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 multiwavelength 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 full width at half-maximum 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. [ABSTRACT FROM AUTHOR]

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 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 parametrized as tanh [(r max  − r)/ l out], where r max  approximates the maximum semimajor axis and l out defines the edge smoothness. If the semimajor axis distribution has sharp edges e rms is roughly 1.2 l out/ r max  or e rms = 0.77 l out/ r max  if semimajor axes have diffused due to self-stirring. This model is fitted to Atacama Large Millimeter/submillimeter Array 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 rms values of 0.121 ± 0.05, |$0.15^{+0.07}_{-0.05}$|⁠ , and 0.10 ± 0.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. [ABSTRACT FROM AUTHOR]

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 five 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 ∼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 |$\sim\! 5{{\ \rm per\ cent}}$| 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. [ABSTRACT FROM AUTHOR]

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 in the range of (2–10) × 10−5 M⊕. 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 signal-to-noise ratio 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. [ABSTRACT FROM AUTHOR]

Despite being >10 Myr, there are ∼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 (C0) that shields CO from photodissociating; i.e. these could be 'shielded secondary discs'. New ALMA observations are presented of HD131835 that detect ∼3 × 10−3 M⊕ of C0, the majority 40–200 au 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 α parameter) and/or the CO production rate is high. Shielding by C0 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 α ∼ 10−3; an inner cavity in C0 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 α = 10−3 to 0.1. The shielded secondary nature of these gas discs can be tested by searching for C0, as well as CN, N2, and CH+, which are also expected to be shielded by C0. [ABSTRACT FROM AUTHOR]

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 (Δap < 20RH, m) low-mass planets are the most efficient at delivering material to exozodi regions ( 5–7 per cent 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), while the exozodi delivery rate stays roughly constant, and that future instruments such as JWST could detect it. In fact for η Corvi, the current Herschel upper limit rules out the scattering scenario by a chain of ≲30 M⊕ planets. Finally, we show that exocomets could be efficient at delivering cometary material to inner planets ( 0.1–1 per cent 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 (∼20RH, m). [ABSTRACT FROM AUTHOR]

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. 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. [ABSTRACT FROM AUTHOR]

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 µm continuum images from the Atacama Large Millimeter Array. The 80 au radius dust ring is resolved radially with a characteristic width of 10 au, consistent with the narrow profile seen in scattered light. Our modelling consistently finds that the disc is also vertically resolved with a similar extent. However, this extent is less than the beam size, and a disc that is dynamically very cold (i.e. vertically thin) provides a better theoretical explanation for the narrowscattered light profile, so we remain cautious about this conclusion. We do not detect 12CO J=3-2 emission, concluding that unless the disc is dynamically cold the CO+CO2 ice content of the planetesimals is of order a few per cent 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 40 au. Independent of our ALMA observations, we note a conflict between mid-IR pericentre-glow and scattered light imaging interpretations, concluding that models where the spatial dust density and grain size vary around the ring should be explored. [ABSTRACT FROM AUTHOR]

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 mag at visible wavelengths within 1 yr and remain bright for decades. EXor objects have lower amplitude outbursts on shorter time-scales. Here we discuss a 1.3 mm Atacama Large Millimeter/submillimeter Array (ALMA) mini-survey of eight outbursting sources (three FUors, four EXors, 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 MJup) than the EXor objects (∼0.5–40 MJup). 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 system in our sample, is shown to host a disc around each one of the stellar components. The disc around HBC 494 is asymmetr [ABSTRACT FROM AUTHOR]