Showing 1–50 of 192 results for author: Wyatt, M C

First ALMA observations of the HD 105211 debris disc: A warm dust component close to a gigayear-old star

Authors: Qiancheng Yang, Qiong Liu, Grant M. Kennedy, Mark C. Wyatt, Sarah Dodson-Robinson, Rachel Akeson, Nenghui Liao

Abstract: Most debris discs consist of a gas-poor, cold dust belt located tens to hundreds of astronomical units away from the host star. Many cold dust belts exhibit distinct structures attributed to the dynamic interaction of planetary systems. Moreover, in a few systems, additional warm components can be found closer to the central star, resembling the asteroid belt or zodiacal dust in our Solar System.… ▽ More Most debris discs consist of a gas-poor, cold dust belt located tens to hundreds of astronomical units away from the host star. Many cold dust belts exhibit distinct structures attributed to the dynamic interaction of planetary systems. Moreover, in a few systems, additional warm components can be found closer to the central star, resembling the asteroid belt or zodiacal dust in our Solar System. In this work, we investigate the structure of the disc surrounding the nearby F2V star HD 105211, which has a warm excess and a potential asymmetry in the cold belt. We applied the CASA pipeline to obtain the ALMA 1.3 mm continuum images. Then we constructed the SED and performed MCMC simulations to fit a model to the ALMA visibility data. To characterise the disc asymmetry, we analysed the ALMA images of two individual observation blocks and compared them to the previous Herschel images. Our modelling reveals that the disc is a narrow ring (23.6+-4.6 au) with low eccentricity positioned at a distance of 133.7+-1.6 au from the central star, which differs from the broad disc (100+-20 au) starting at an inner edge of 87+-2.5 au, inferred from the Herschel images. We found that both observation blocks show excess emission at the stellar position, while OB1 shows an offset between the star and the phase centre, and OB2 shows brightness clumps. We used a two-temperature model to fit the infrared SED and used the ALMA detection to constrain the warm component to a nearly pure blackbody model. The relatively low ratio of actual radius to blackbody radius of the HD105211 debris disc indicates that this system is depleted in small grains, which could indicate that it is dynamically cold. The excess emission from the stellar position suggests that there should be a warm mm-sized dust component close to the star, for which we suggest two possible origins: in situ asteroid belt or comet delivery. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: 12 pages, 7 figures, accepted for publication in A&A

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems V: Do Self-Consistent Atmospheric Models Represent JWST Spectra? A Showcase With VHS 1256 b

Authors: Simon Petrus, Niall Whiteford, Polychronis Patapis, Beth A. Biller, Andrew Skemer, Sasha Hinkley, Genaro Suárez, Anna Lueber, Paulina Palma-Bifani, Jordan M. Stone, Johanna M. Vos, Caroline V. Morley, Pascal Tremblin, Benjamin Charnay, Christiane Helling, Brittany E. Miles, Aarynn L. Carter, Jason J. Wang, Markus Janson, Eileen C. Gonzales, Ben Sutlieff, Kielan K. W. Hoch, Mickaël Bonnefoy, Gaël Chauvin, Olivier Absil , et al. (97 additional authors not shown)

Abstract: 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. W… ▽ More 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. △ Less

Submitted 31 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 32 pages, 16 figures, 6 tables, 2 appendices

The effect of sculpting planets on the steepness of debris-disc inner edges

Authors: Tim D. Pearce, Alexander V. Krivov, Antranik A. Sefilian, Marija R. Jankovic, Torsten Löhne, Tobias Morgner, Mark C. Wyatt, Mark Booth, Sebastian Marino

Abstract: 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 ra… ▽ More 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. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 24 pages, 18 figures, accepted for publication in MNRAS

An ALMA Survey of M-dwarfs in the Beta Pictoris Moving Group with Two New Debris Disc Detections

Authors: Patrick F. Cronin-Coltsmann, Grant M. Kennedy, Quentin Kral, Jean-François Lestrade, Sebastian Marino, Luca Matrà, Mark C. Wyatt

Abstract: 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 lig… ▽ More 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\,$μ$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. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS

The \textit{JWST} Early Release Science Program for Direct Observations of Exoplanetary Systems III: Aperture Masking Interferometric Observations of the star HIP\,65426 at $\boldsymbol{3.8\,\rm{μm}}$

Authors: Shrishmoy Ray, Steph Sallum, Sasha Hinkley, Anand Sivamarakrishnan, Rachel Cooper, Jens Kammerer, Alexandra Z. Greebaum, Deepashri Thatte, Cecilia Lazzoni, Andrei Tokovinin, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan , et al. (98 additional authors not shown)

Abstract: We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{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λ/D$ f… ▽ More We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{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λ/D$ for an interferometer), which are inaccessible with the classical inner working angles of the \textit{JWST} coronagraphs. When combined with \textit{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 contrast of $Δm_{F380M}{\sim }7.8$\,mag relative to the host star at a separation of ${\sim}0.07\arcsec$ but detect no additional companions interior to the known companion HIP\,65426\,b. Our observations thus rule out companions more massive than $10{-}12\,\rm{M\textsubscript{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 \textit{JWST} is sensitive to planetary mass companions orbiting at the water frost line, even for more distant stars at $\sim$100\,pc. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening essentially unexplored parameter space. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 15 pages, 9 figures, submitted to ApJ Letters

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems IV: NIRISS Aperture Masking Interferometry Performance and Lessons Learned

Authors: Steph Sallum, Shrishmoy Ray, Jens Kammerer, Anand Sivaramakrishnan, Rachel Cooper, Alexandra Z. Greebaum, Deepashri Thatte, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Sasha Hinkley, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan, Sarah K. Betti, Anthony Boccaletti , et al. (98 additional authors not shown)

Abstract: 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… ▽ More 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 λ/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. △ Less

Submitted 11 March, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

The Eccentric Kozai-Lidov Mechanism as the Cause of Exocomet Transits of KIC 8462852

Authors: Steven D. Young, Mark C. Wyatt

Abstract: KIC 8462852 is a star in the Kepler field that exhibits almost unique behaviour. The deep, irregular and aperiodic dips in its light curve have been interpreted as the breakup of a large exocomet on a highly eccentric orbit whose post-disruption material obscures the star. It is hypothesised that a nearby M-dwarf, recently confirmed to be bound to the system, could be exciting planetesimals in a s… ▽ More KIC 8462852 is a star in the Kepler field that exhibits almost unique behaviour. The deep, irregular and aperiodic dips in its light curve have been interpreted as the breakup of a large exocomet on a highly eccentric orbit whose post-disruption material obscures the star. It is hypothesised that a nearby M-dwarf, recently confirmed to be bound to the system, could be exciting planetesimals in a source belt to high eccentricities if its orbit is highly misaligned with the belt: an effect known as the 'Eccentric Kozai-Lidov Mechanism'. To quantify how often this effect is expected to occur, this paper presents a Monte Carlo model of wide binary stars with embedded, misaligned planetesimal belts. These belts collisionally erode over time until they are excited to high eccentricities on secular timescales by a companion star if its orbit is sufficiently misaligned. The large planetesimals then produce an observable dimming signature in the light curve for a set period of time which may or may not overlap with similar events. The model finds that, for dimming events that persist for 100 yr, the most likely companion stars are located at $10^2 - 10^4$ au, the most likely belts are at $10^2-10^3$ au and the system age is most likely to be $10^2 - 10^3$ Myr. However, the probability of observing one or more stars exhibiting this phenomenon in the Kepler field is $1.3 \times 10^{-3}$, such that it is unlikely this mechanism is driving the observations of KIC 8462852. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 23 pages, 24 figures and 2 tables. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Secondary gas in debris discs released following the decay of long-lived radioactive nuclides, catastrophic or resurfacing collisions

Authors: Amy Bonsor, Mark C. Wyatt, Sebastian Marino, Björn J. R. Davidsson, Quentin Kral, Philippe Thebault

Abstract: 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 dec… ▽ More 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. △ Less

Submitted 22 September, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: Accepted by MNRAS, 16 pages. Updated following review comments

Deprojecting and constraining the vertical thickness of exoKuiper belts

Authors: James Terrill, Sebastian Marino, Richard A. Booth, Yinuo Han, Jeff Jennings, Mark C. Wyatt

Abstract: 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 brig… ▽ More 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. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS. 17 pages. 16 figures

Formation of Gaps in Self-gravitating Debris Disks by Secular Resonance in a Single-planet System. II. Towards a Self-consistent Model

Authors: Antranik A. Sefilian, Roman R. Rafikov, Mark C. Wyatt

Abstract: High-resolution observations of several debris disks reveal structures such as gaps and spirals, suggestive of gravitational perturbations induced by underlying planets. Most existing studies of planet--debris disk interactions ignore the gravity of the disk, treating it as a reservoir of massless planetesimals. In this paper, we continue our investigation into the long-term interaction between a… ▽ More High-resolution observations of several debris disks reveal structures such as gaps and spirals, suggestive of gravitational perturbations induced by underlying planets. Most existing studies of planet--debris disk interactions ignore the gravity of the disk, treating it as a reservoir of massless planetesimals. In this paper, we continue our investigation into the long-term interaction between a single eccentric planet and an external, massive debris disk. Building upon our previous work, here we consider not only the axisymmetric component of the disk's gravitational potential, but also the non-axisymmetric torque that the disk exerts on the planet (ignoring for now only the non-axisymmetric component of the disk self-gravity). To this goal, we develop and test a semi-analytic `$N$-ring' framework that is based on a generalized (softened) version of the classical Laplace--Lagrange secular theory. Using this tool, we demonstrate that even when the disk is less massive than the planet, not only can a secular resonance be established within the disk that leads to the formation of a wide gap, but that the very same resonance also damps the planetary eccentricity $e_p$ via a process known as resonant friction. The resulting gap is initially non-axisymmetric (akin to those observed in HD 92945 and HD 206893), but evolves to become more axisymmetric (similar to that in HD 107146) as $e_p(t)\rightarrow0$ with time. We also develop analytic understanding of these findings, finding good quantitative agreement with the outcomes of the $N$-ring calculations. Our results may be used to infer both the dynamical masses of (gapped) debris disks and the dynamical history of the planets interior to them, as we exemplify for HD 206893. △ Less

Submitted 14 July, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: Accepted version (ApJ): 35 Pages (including 8 pages of Appendices), 15 Figures, 1 Table, 3 Animations (see Ancillary files). A copy of the N-ring code is made publicly available at http://doi.org/10.6084/m9.figshare.22598377

Inner edges of planetesimal belts: collisionally eroded or truncated?

Authors: Amaia Imaz Blanco, Sebastian Marino, Luca Matrà, Mark Booth, John Carpenter, Virginie Faramaz, Thomas Henning, A. Meredith Hughes, Grant M. Kennedy, Sebastián Pérez, Luca Ricci, Mark C. Wyatt

Abstract: 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… ▽ More 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. △ Less

Submitted 16 May, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS, 21 pages, 11 figures

Stirred but not shaken: a multi-wavelength view of HD 16743's debris disc

Authors: Jonathan P. Marshall, Julien Milli, Elodie Choquet, Carlos del Burgo, Grant M. Kennedy, Francisca Kemper, Mark C. Wyatt, Quentin Kral, Remi Soummer

Abstract: Planetesimals -- asteroids and comets -- are the building blocks of planets in protoplanetary discs and the source of dust, ice and gas in debris discs. Along with planets they comprise the left-over material after star formation that constitutes a planetary system. Planets influence the dynamics of planetesimals, sculpting the orbits of debris belts to produce asymmetries or gaps. We can constrai… ▽ More Planetesimals -- asteroids and comets -- are the building blocks of planets in protoplanetary discs and the source of dust, ice and gas in debris discs. Along with planets they comprise the left-over material after star formation that constitutes a planetary system. Planets influence the dynamics of planetesimals, sculpting the orbits of debris belts to produce asymmetries or gaps. We can constrain the architecture of planetary systems, and infer the presence of unseen planetary companions, by high spatial resolution imaging of debris discs. HD~16743 is a relatively young F-type star that hosts a bright edge-on debris disc. Based on far-infrared \textit{Herschel} observations its disc was thought to be stirred by a planetary companion. Here we present the first spatially resolved observations at near-infrared and millimetre wavelengths with \textit{HST} and ALMA, revealing the disc to be highly inclined at $87\fdg3~^{+1\fdg9}_{-2\fdg5}$ with a radial extent of 157.7$^{+2.6}_{-1.5}$~au and a FWHM of 79.4$^{+8.1}_{-7.8}$~au ($ΔR/R = 0.5$). The vertical scale height of the disc is $0.13~\pm~0.02$, significantly greater than typically assumed unstirred value of 0.05, and could be indicative of stirring of the dust-producing planetesimals within the disc by bodies at least a few times the mass of Pluto up to 18.3~$M_{\oplus}$ in the single object limit. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 12 pages, 7 figures, 4 tables, accepted for publication in MNRAS

The clumpy structure of $ε$ Eridani's debris disc revisited by ALMA

Authors: Mark Booth, Tim D. Pearce, Alexander V. Krivov, Mark C. Wyatt, William R. F. Dent, Antonio S. Hales, Jean-François Lestrade, Fernando Cruz-Sáenz de Miera, Virginie C. Faramaz, Torsten Löhne, Miguel Chavez-Dagostino

Abstract: $ε… ▽ More $ε$ Eridani is the closest star to our Sun known to host a debris disc. Prior observations in the (sub-)millimetre regime have potentially detected clumpy structure in the disc and attributed this to interactions with an (as yet) undetected planet. However, the prior observations were unable to distinguish between structure in the disc and background confusion. Here we present the first ALMA image of the entire disc, which has a resolution of 1.6"$\times$1.2". We clearly detect the star, the main belt and two point sources. The resolution and sensitivity of this data allow us to clearly distinguish background galaxies (that show up as point sources) from the disc emission. We show that the two point sources are consistent with background galaxies. After taking account of these, we find that resolved residuals are still present in the main belt, including two clumps with a $>3σ$ significance -- one to the east of the star and the other to the northwest. We perform $n$-body simulations to demonstrate that a migrating planet can form structures similar to those observed by trapping planetesimals in resonances. We find that the observed features can be reproduced by a migrating planet trapping planetesimals in the 2:1 mean motion resonance and the symmetry of the most prominent clumps means that the planet should have a position angle of either ${\sim10^\circ}$ or ${\sim190^\circ}$. Observations over multiple epochs are necessary to test whether the observed features rotate around the star. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Has the dust clump in the debris disk of Beta Pictoris moved?

Authors: Yinuo Han, Mark C. Wyatt, William R. F. Dent

Abstract: The edge-on debris disk of the nearby young star Beta Pictoris shows an unusual brightness asymmetry in the form of a clump. The clump has been detected in both the mid-IR and CO and its origin has so far remained uncertain. Here we present new mid-IR observations of Beta Pic to track any motion of the dust clump. Together with previous observations, the data span a period of 12 years. We measured… ▽ More The edge-on debris disk of the nearby young star Beta Pictoris shows an unusual brightness asymmetry in the form of a clump. The clump has been detected in both the mid-IR and CO and its origin has so far remained uncertain. Here we present new mid-IR observations of Beta Pic to track any motion of the dust clump. Together with previous observations, the data span a period of 12 years. We measured any projected displacement of the dust clump over the 12-yr period to be $0.2^{+1.3}_{-1.4}$ au away from the star based on the median and 1$σ$ uncertainty, and constrain this displacement to be <11 au at the 3$σ$ level. This implies that the observed motion is incompatible with Keplerian motion at the 2.8$σ$ level. It has been posited that a planet migrating outwards may trap planetesimals into a 2:1 resonance, resulting in the observed clump at pericentre of their orbits that trails the planet. The observed motion is also incompatible with such resonant motion at the 2.6$σ$ level. While Keplerian motion and resonant motion is still possible, the data suggest that the dust clump is more likely stationary. Such a stationary dust clump could originate from the collision or tidal disruption of a planet-sized body, or from secular perturbations due to a planet that create regions with enhanced densities in the disk. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 14 pages, 8 figures, 5 tables. Published in MNRAS

Journal ref: MNRAS 519(3), 3257-3270 (2023)

High resolution ALMA and HST imaging of $κ$CrB: a broad debris disc around a post-main sequence star with low-mass companions

Authors: J. B. Lovell, M. C. Wyatt, P. Kalas, G. M. Kennedy, S. Marino, A. Bonsor, Z. Penoyre, B. J. Fulton, N. Pawellek

Abstract: $κ$CrB is a ${\sim}2.5\,$Gyr old K1 sub-giant star, with an eccentric exo-Jupiter at ${\sim}2.8\,$au and a debris disc at tens of au. We present ALMA Band 6 ($1.3\,$mm) and HST scattered light ($0.6\,μ$m) images, demonstrating $κ$CrB's broad debris disc, covering an extent $50{-}180\,$au in the millimetre (peaking at $110\,$au), and $51{-}280\,$au in scattered light (peaking at $73\,… ▽ More $κ$CrB is a ${\sim}2.5\,$Gyr old K1 sub-giant star, with an eccentric exo-Jupiter at ${\sim}2.8\,$au and a debris disc at tens of au. We present ALMA Band 6 ($1.3\,$mm) and HST scattered light ($0.6\,μ$m) images, demonstrating $κ$CrB's broad debris disc, covering an extent $50{-}180\,$au in the millimetre (peaking at $110\,$au), and $51{-}280\,$au in scattered light (peaking at $73\,$au). By modelling the millimetre emission, we estimate the dust mass as ${\sim}0.016\,M{\oplus}$, and constrain lower-limit planetesimal sizes as $D_{\rm{max}}{>}1\,$km and the planetesimal belt mass as $M_{\rm{disc}}{>}1\,M_{\oplus}$. We constrain the properties of an outer body causing a linear trend in 17 years of radial velocity data to have a semi-major axis $8{-}66\,$au and a mass $(0.4{-}120)\,M_{\rm{Jup}}$. There is a large inner cavity seen in the millimetre emission, which we show is consistent with carving by such an outer massive companion with a string of lower mass planets. Our scattered light modelling shows that the dust must have a high anisotropic scattering factor ($g{\sim}0.8{-}0.9$) but an inclination ($i{\sim}30{-}40\,$degrees) that is inferred to be significantly lower than the $i{\sim}61\,$degrees millimetre inclination. The origin of such a discrepancy is unclear, but could be caused by a misalignment in the micron and millimetre sized dust. We place an upper limit on the CO gas mass of $M_{\rm{CO}}{<}(4.2{-}13){\times}10^{-7}\,M_{\oplus}$, and show this to be consistent with levels expected from planetesimal collisions, or from CO-ice sublimation as $κ$CrB begins its giant branch ascent. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS. 21 pages, 14 figures

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems II: A 1 to 20 Micron Spectrum of the Planetary-Mass Companion VHS 1256-1257 b

Authors: Brittany E. Miles, Beth A. Biller, Polychronis Patapis, Kadin Worthen, Emily Rickman, Kielan K. W. Hoch, Andrew Skemer, Marshall D. Perrin, Niall Whiteford, Christine H. Chen, B. Sargent, Sagnick Mukherjee, Caroline V. Morley, Sarah E. Moran, Mickael Bonnefoy, Simon Petrus, Aarynn L. Carter, Elodie Choquet, Sasha Hinkley, Kimberly Ward-Duong, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Laurent Pueyo, Shrishmoy Ray, Karl R. Stapelfeldt , et al. (79 additional authors not shown)

Abstract: 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… ▽ More 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 $μ$m to 20 $μ$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. △ Less

Submitted 4 July, 2024; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 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

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems I: High Contrast Imaging of the Exoplanet HIP 65426 b from 2-16 $μ$m

Authors: Aarynn L. Carter, Sasha Hinkley, Jens Kammerer, Andrew Skemer, Beth A. Biller, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Simon Petrus, Jordan M. Stone, Kimberly Ward-Duong, Jason J. Wang, Julien H. Girard, Dean C. Hines, Marshall D. Perrin, Laurent Pueyo, William O. Balmer, Mariangela Bonavita, Mickael Bonnefoy, Gael Chauvin, Elodie Choquet, Valentin Christiaens, Camilla Danielski, Grant M. Kennedy, Elisabeth C. Matthews, Brittany E. Miles , et al. (86 additional authors not shown)

Abstract: 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 exo… ▽ More 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 $μ$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$σ$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $μ$m and MIRI at 11.3 $μ$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 $μ$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. △ Less

Submitted 3 May, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

Authors: Devin Sullivan, David Wilner, Luca Matra, Mark C. Wyatt, Sean M. Andrews, Meredith A. MacGregor, Brenda Matthews

Abstract: Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust… ▽ More Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.5 arcsec (200 au) and a median rms noise of 54 mircoJy/beam, which corresponds to a fractional luminosity $L_{dust}/L_{star} \sim 10^{-4}$ for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations.The ALMA non-detections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (HII 1132 and HD 22680) show 24 micron excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, accepted for publication in AJ

Vertical evolution of exocometary gas: I. How vertical diffusion shortens the CO lifetime

Authors: S. Marino, G. Cataldi, M. R. Jankovic, L. Matrà, M. C. Wyatt

Abstract: Bright debris discs can contain large amounts of CO gas. This gas was thought to be a protoplanetary remnant until it was recently shown that it could be released in collisions of volatile-rich solids. As CO is released, interstellar UV radiation photodissociates CO producing CI, which can shield CO allowing a large CO mass to accumulate. However, this picture was challenged because CI is ineffici… ▽ More Bright debris discs can contain large amounts of CO gas. This gas was thought to be a protoplanetary remnant until it was recently shown that it could be released in collisions of volatile-rich solids. As CO is released, interstellar UV radiation photodissociates CO producing CI, which can shield CO allowing a large CO mass to accumulate. However, this picture was challenged because CI is inefficient at shielding if CO and CI are vertically mixed. Here, we study for the first time the vertical evolution of gas to determine how vertical mixing affects the efficiency of shielding by CI. We present a 1D model that accounts for gas release, photodissociation, ionisation, viscous evolution, and vertical mixing due to turbulent diffusion. We find that if the gas surface density is high and the vertical diffusion weak ($α_{\rm v}/α<[H/r]^2$) CO photodissociates high above the midplane, forming an optically thick CI layer that shields the CO underneath. Conversely, if diffusion is strong ($α_{\rm v}/α>[H/r]^2$) CI and CO become well mixed, shortening the CO lifetime. Moreover, diffusion could also limit the amount of dust settling. High-resolution ALMA observations could resolve the vertical distribution of CO and CI, and thus constrain vertical mixing and the efficiency of CI shielding. We also find that the CO and CI scale heights may not be good probes of the mean molecular weight, and thus composition, of the gas. Finally, we show that if mixing is strong the CO lifetime might not be long enough for CO to spread interior to the planetesimal belt where gas is produced. △ Less

Submitted 27 June, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS (18 pages, 14 figures)

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

Authors: Charles L. H. Hull, Haifeng Yang, Paulo C. Cortés, William R. F. Dent, Quentin Kral, Zhi-Yun Li, Valentin J. M. Le Gouellec, A. Meredith Hughes, Julien Milli, Richard Teague, Mark C. Wyatt

Abstract: We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polariza… ▽ More We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $χ$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $β$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 22 pages, 13 figures, 2 tables. Accepted for publication in the Astrophysical Journal. Materials accessible in the online version of the ApJ article include the FITS files used to make the ALMA images

ALMA's view of the M-dwarf GSC 07396-00759's edge-on debris disc: AU Mic's coeval twin

Authors: Patrick F. Cronin-Coltsmann, Grant M. Kennedy, Christian Adam, Quentin Kral, Jean-François Lestrade, Sebastian Marino, Luca Matrà, Simon J. Murphy, Johan Olofsson, Mark C. Wyatt

Abstract: 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 sys… ▽ More 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. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 14 pages, 11 figures, 1 table, Accepted for publication in MNRAS

Rave: A non-parametric method for recovering the surface brightness and height profiles of edge-on debris disks

Authors: Yinuo Han, Mark C. Wyatt, Luca Matra

Abstract: Extrasolar analogues of the Solar System's Kuiper belt offer unique constraints on outer planetary system architecture. Radial features such as the sharpness of disk edges and substructures such as gaps may be indicative of embedded planets within a disk. Vertically, the height of a disk can constrain the mass of embedded bodies. Edge-on debris disks offer a unique opportunity to simultaneously ac… ▽ More Extrasolar analogues of the Solar System's Kuiper belt offer unique constraints on outer planetary system architecture. Radial features such as the sharpness of disk edges and substructures such as gaps may be indicative of embedded planets within a disk. Vertically, the height of a disk can constrain the mass of embedded bodies. Edge-on debris disks offer a unique opportunity to simultaneously access the radial and vertical distribution of material, however recovering either distribution in an unbiased way is challenging. In this study, we present a non-parametric method to recover the surface brightness profile (face-on surface brightness as a function of radius) and height profile (scale height as a function of radius) of azimuthally symmetric, edge-on debris disks. The method is primarily designed for observations at thermal emission wavelengths, but is also applicable to scattered light observations under the assumption of isotropic scattering. By removing assumptions on underlying functional forms, this algorithm provides more realistic constraints on disk structures. We also apply this technique to ALMA observations of the AU Mic debris disk and derive a surface brightness profile consistent with estimates from parametric approaches, but with a more realistic range of possible models that is independent of parametrisation assumptions. Our results are consistent with a uniform scale height of 0.8 au, but a scale height that increases linearly with radius is also possible. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

Comet fragmentation as a source of the zodiacal cloud

Authors: Jessica K. Rigley, Mark C. Wyatt

Abstract: Models of the zodiacal cloud's thermal emission and sporadic meteoroids suggest Jupiter-family comets (JFCs) as the dominant source of interplanetary dust. However, comet sublimation is insufficient to sustain the quantity of dust presently in the inner solar system, suggesting that spontaneous disruptions of JFCs may supply the zodiacal cloud. We present a model for the dust produced in comet fra… ▽ More Models of the zodiacal cloud's thermal emission and sporadic meteoroids suggest Jupiter-family comets (JFCs) as the dominant source of interplanetary dust. However, comet sublimation is insufficient to sustain the quantity of dust presently in the inner solar system, suggesting that spontaneous disruptions of JFCs may supply the zodiacal cloud. We present a model for the dust produced in comet fragmentations and its evolution. Using results from dynamical simulations, the model follows individual comets drawn from a size distribution as they evolve and undergo recurrent splitting events. The resulting dust is followed with a kinetic model which accounts for the effects of collisional evolution, Poynting-Robertson drag, and radiation pressure. This allows to model the evolution of both the size distribution and radial profile of dust, and we demonstrate the importance of including collisions (both as a source and sink of dust) in zodiacal cloud models. With physically-motivated free parameters this model provides a good fit to zodiacal cloud observables, supporting comet fragmentation as the plausibly dominant dust source. The model implies that dust in the present zodiacal cloud likely originated primarily from disruptions of $\sim$ 50 km comets, since larger comets are ejected before losing all their mass. Thus much of the dust seen today was likely deposited as larger grains $\sim$ 0.1 Myr in the past. The model also finds the dust level to vary stochastically; e.g., every $\sim$ 50 Myr large (>100 km) comets with long dynamical lifetimes inside Jupiter cause dust spikes with order of magnitude increases in zodiacal light brightness lasting $\sim$ 1 Myr. If exozodiacal dust is cometary in origin, our model suggests it should be similarly variable. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 25 pages, 28 figures, 2 tables. Accepted for publication in MNRAS

Carbon monoxide gas produced by a giant impact in the inner region of a young system

Authors: Tajana Schneiderman, Luca Matrà, Alan P. Jackson, Grant M. Kennedy, Quentin Kral, Sebastián Marino, Karin I. Öberg, Kate Y. L. Su, David J. Wilner, Mark C. Wyatt

Abstract: 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 impact… ▽ More 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. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 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

Journal ref: Nature 598, 425-428 (2021)

Debris Disks

Authors: Mark C. Wyatt

Abstract: Debris disk is a catch-all term that can be used to refer to any component of a planetary system which is not an actual planet. In the Solar System this refers to the asteroids and comets in the Asteroid and Kuiper belts as well as the dust and gas derived from them, such as the zodiacal cloud. Studying the structure of extrasolar debris disks provides unique constraints on the underlying planetar… ▽ More Debris disk is a catch-all term that can be used to refer to any component of a planetary system which is not an actual planet. In the Solar System this refers to the asteroids and comets in the Asteroid and Kuiper belts as well as the dust and gas derived from them, such as the zodiacal cloud. Studying the structure of extrasolar debris disks provides unique constraints on the underlying planetary system and on the processes of planet formation and protoplanetary disk evolution. Debris disks also have important implications for processes affecting the planets, such as impact events that may strip an atmosphere or deliver volatiles. Their presence or absence also has consequences both positive and negative for the detectability of exoplanets in the system. Here I describe five of the most important ways in which debris disks will contribute to our understanding of extrasolar planetary systems in the coming years. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: To appear as a book chapter in "ExoFrontiers: Big questions in exoplanetary science", Ed. N Madhusudhan (Bristol: IOP Publishing Ltd) AAS-IOP ebooks https://iopscience.iop.org/bookListInfo/aas-iop-astronomy

Relentless and Complex Transits from a Planetesimal Debris Disk

Authors: J. Farihi, J. J. Hermes, T. R. Marsh, A. J. Mustill, M. C. Wyatt, J. A. Guidry, T. G. Wilson, S. Redfield, P. Izquierdo, O. Toloza, B. T. Gänsicke, A. Aungwerojwit, V. S. Dhillon, A. Swan

Abstract: This article reports quasi-continuous transiting events towards WD 1054-226 at d=36.2 pc and V=16.0 mag, based on simultaneous, high-cadence, multi-wavelength imaging photometry using ULTRACAM over 18 nights from 2019 to 2020 March. The predominant period is 25.02 h, and corresponds to a circular orbit with blackbody Teq = 323 K, where a planetary surface can nominally support liquid water. The li… ▽ More This article reports quasi-continuous transiting events towards WD 1054-226 at d=36.2 pc and V=16.0 mag, based on simultaneous, high-cadence, multi-wavelength imaging photometry using ULTRACAM over 18 nights from 2019 to 2020 March. The predominant period is 25.02 h, and corresponds to a circular orbit with blackbody Teq = 323 K, where a planetary surface can nominally support liquid water. The light curves reveal remarkable night-to-night similarity, with changes on longer timescales, and lack any transit-free segments of unocculted starlight. The most pronounced dimming components occur every 23.1 min -- exactly the 65th harmonic of the fundamental period -- with depths of up to several per cent, and no evident color dependence. Myriad additional harmonics are present, as well as at least two transiting features with independent periods. High-resolution optical spectra are consistent with stable, photospheric absorption by multiple, refractory metal species, with no indication of circumstellar gas. Spitzer observations demonstrate a lack of detectable dust emission, suggesting that the otherwise hidden circumstellar disk orbiting WD 1054-226 may be typical of polluted white dwarfs, and only detected via favorable geometry. Future observations are required to constrain the orbital eccentricity, but even if periastron is near the Roche limit, sublimation cannot drive mass loss in refractory parent bodies, and collisional disintegration is necessary for dust production. △ Less

Submitted 5 January, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 22 pages, 2 tables, and 15 figures including the appendix, accepted to MNRAS

High resolution ALMA and HST images of q$^1$ Eri: an asymmetric debris disc with an eccentric Jupiter

Authors: J. B. Lovell, S. Marino, M. C. Wyatt, G. M. Kennedy, M. A. MacGregor, K. Stapelfeldt, B. Dent, J. Krist, L. Matrà, Q. Kral, O. Panić, T. D. Pearce, D. Wilner

Abstract: We present \textit{ALMA} 1.3 mm and 0.86 mm observations of the nearby (17.34 pc) F9V star q1 Eri (HD 10647, HR 506). This system, with age ${\sim}1.4$ Gyr, hosts a ${\sim}2$ au radial velocity planet and a debris disc with the highest fractional luminosity of the closest 300 FGK type stars. The \textit{ALMA} images, with resolution ${\sim}0.5''$, reveal a broad (34{-}134 au) belt of millimeter em… ▽ More We present \textit{ALMA} 1.3 mm and 0.86 mm observations of the nearby (17.34 pc) F9V star q1 Eri (HD 10647, HR 506). This system, with age ${\sim}1.4$ Gyr, hosts a ${\sim}2$ au radial velocity planet and a debris disc with the highest fractional luminosity of the closest 300 FGK type stars. The \textit{ALMA} images, with resolution ${\sim}0.5''$, reveal a broad (34{-}134 au) belt of millimeter emission inclined by $76.7{\pm}1.0$ degrees with maximum brightness at $81.6{\pm}0.5$ au. The images reveal an asymmetry, with higher flux near the southwest ansa, which is also closer to the star. Scattered light observed with the Hubble Space Telescope is also asymmetric, being more radially extended to the northeast. We fit the millimeter emission with parametric models and place constraints on the disc morphology, radius, width, dust mass, and scale height. We find the southwest ansa asymmetry is best fitted by an extended clump on the inner edge of the disc, consistent with perturbations from a planet with mass $8 M_{\oplus} {-} 11 M_{\rm Jup}$ at ${\sim}60$ au that may have migrated outwards, similar to Neptune in our Solar System. If the measured vertical aspect ratio of $h{=}0.04{\pm}0.01$ is due to dynamical interactions in the disc, then this requires perturbers with sizes ${>}1200$ km. We find tentative evidence for an 0.86 mm excess within 10 au, $70{\pm}22\, μ$Jy, that may be due to an inner planetesimal belt. We find no evidence for CO gas, but set an upper bound on the CO gas mass of $4{\times}10^{-6}$ M$_{\oplus}$ ($3\,σ$), consistent with cometary abundances in the Solar System. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS. Paper: 21 pages, appendix: 4 pages. 16 figures

ALMA imaging of the M-dwarf Fomalhaut C's debris disc

Authors: Patrick F. Cronin-Coltsmann, Grant M. Kennedy, Paul Kalas, Julien Milli, Cathie J. Clarke, Gaspard Duchêne, Jane Greaves, Samantha M. Lawler, Jean-François Lestrade, Brenda C. Matthews, Andrew Shannon, Mark C. Wyatt

Abstract: Fomalhaut C (LP 876-10) is a low mass M4V star in the intriguing Fomalhaut triple system and, like Fomalhaut A, possesses a debris disc. It is one of very few nearby M-dwarfs known to host a debris disc and of these has by far the lowest stellar mass. We present new resolved observations of the debris disc around Fomalhaut C with the Atacama Large Millimetre Array which allow us to model its prope… ▽ More Fomalhaut C (LP 876-10) is a low mass M4V star in the intriguing Fomalhaut triple system and, like Fomalhaut A, possesses a debris disc. It is one of very few nearby M-dwarfs known to host a debris disc and of these has by far the lowest stellar mass. We present new resolved observations of the debris disc around Fomalhaut C with the Atacama Large Millimetre Array which allow us to model its properties and investigate the system's unique history. The ring has a radius of 26 au and a narrow full width at half maximum of at most 4.2 au. We find a 3$σ$ upper limit on the eccentricity of 0.14, neither confirming nor ruling out previous dynamic interactions with Fomalhaut A that could have affected Fomalhaut C's disc. We detect no $^{12}$CO J=3-2 emission in the system and do not detect the disc in scattered light with HST/STIS or VLT/SPHERE. We find the original Herschel detection to be consistent with our ALMA model's radial size. We place the disc in the context of the wider debris disc population and find that its radius is as expected from previous disc radius-host luminosity trends. Higher signal-to-noise observations of the system would be required to further constrain the disc properties and provide further insight to the history of the Fomalhaut triple system as a whole. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 15 pages, 10 figures, 2 tables, Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 504, Issue 3, July 2021, Pages 4497-4510

A molecular wind blows out of the Kuiper belt

Authors: Quentin Kral, J. E. Pringle, Aurélie Guilbert-Lepoutre, Luca Matrà, Julianne I. Moses, Emmanuel Lellouch, Mark C. Wyatt, Nicolas Biver, Dominique Bockelée-Morvan, Amy Bonsor, Franck Le Petit, G. Randall Gladstone

Abstract: Gas has been detected in many exoplanetary systems ($>$10 Myr), thought to be released in the destruction of volatile-rich planetesimals orbiting in exo-Kuiper belts. In this letter, we aim to explore whether gas is also expected in the Kuiper belt (KB) in our Solar System. To quantify the gas release in our Solar System, we use models for gas release that have been applied to extrasolar planetary… ▽ More Gas has been detected in many exoplanetary systems ($>$10 Myr), thought to be released in the destruction of volatile-rich planetesimals orbiting in exo-Kuiper belts. In this letter, we aim to explore whether gas is also expected in the Kuiper belt (KB) in our Solar System. To quantify the gas release in our Solar System, we use models for gas release that have been applied to extrasolar planetary systems, as well as a physical model that accounts for gas released due to the progressive internal warming of large planetesimals. We find that only bodies larger than about 4 km can still contain CO ice after 4.6 Gyr of evolution. This finding may provide a clue as to why Jupiter-family comets, thought to originate in the Kuiper belt, are deficient in CO compared to Oort-clouds comets. We predict that gas is still produced in the KB right now at a rate of $2 \times 10^{-8}$ M$_\oplus$/Myr for CO and orders of magnitude more when the Sun was younger. Once released, the gas is quickly pushed out by the Solar wind. Therefore, we predict a gas wind in our Solar System starting at the KB location and extending far beyond with regards to the heliosphere with a current total CO mass of $\sim 2 \times 10^{-12}$ M$_\oplus$. We also predict the existence of a slightly more massive atomic gas wind made of carbon and oxygen (neutral and ionized) with a mass of $\sim 10^{-11}$ M$_\oplus$. We predict that gas is currently present in our Solar System beyond the Kuiper belt and that although it cannot be detected with current instrumentation, it could be observed in the future with an in situ mission using an instrument similar to Alice on New Horizons with larger detectors. Our model of gas release due to slow heating may also work for exoplanetary systems and provide the first real physical mechanism for the gas observations. △ Less

Submitted 9 September, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: accepted for publication as a Letter to the editor in A&A; abstract shortened; 15 pages, 5 figures, 4 tables

Journal ref: A&A 653, L11 (2021)

The HOSTS survey: evidence for an extended dust disk and constraints on the presence of giant planets in the Habitable Zone of $β$ Leo

Authors: D. Defrère, P. M. Hinz, G. M. Kennedy, J. Stone, J. Rigley, S. Ertel, A. Gaspar, V. P. Bailey, W. F. Hoffmann, B. Mennesson, R. Millan-Gabet, W. C. Danchi, O. Absil, P. Arbo, C. Beichman, M. Bonavita, G. Brusa, G. Bryden, E. C. Downey, S. Esposito, P. Grenz, C. Haniff, J. M. Hill, J. M. Leisenring, J. R. Males , et al. (16 additional authors not shown)

Abstract: The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Tel… ▽ More The young (50-400 Myr) A3V star $β$ Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot ($\sim$1600$^\circ$K), warm ($\sim$600$^\circ$K), and cold ($\sim$120$^\circ$K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47\%$\pm$0.050\% within the central 1.5 au, rising to 0.81\%$\pm$0.026\% within 4.5 au, outside the habitable zone of $β$~Leo. This dust level is 50 $\pm$ 10 times greater than in the solar system's zodiacal cloud. Poynting-Robertson drag on the cold dust detected by Spitzer and Herschel under-predicts the dust present in the habitable zone of $β$~Leo, suggesting an additional delivery mechanism (e.g.,~comets) or an additional belt at $\sim$5.5 au. A model of these dust components is provided which implies the absence of planets more than a few Saturn masses between $\sim$5 au and the outer belt at $\sim$40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8~$μ$m wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5 au to 50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comments: 11 pages, 9 figures, accepted for publication in Astronomical Journal

Large Interferometer For Exoplanets (LIFE): I. Improved exoplanet detection yield estimates for a large mid-infrared space-interferometer mission

Authors: S. P. Quanz, M. Ottiger, E. Fontanet, J. Kammerer, F. Menti, F. Dannert, A. Gheorghe, O. Absil, V. S. Airapetian, E. Alei, R. Allart, D. Angerhausen, S. Blumenthal, L. A. Buchhave, J. Cabrera, Ó. Carrión-González, G. Chauvin, W. C. Danchi, C. Dandumont, D. Defrère, C. Dorn, D. Ehrenreich, S. Ertel, M. Fridlund, A. García Muñoz , et al. (46 additional authors not shown)

Abstract: One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measur… ▽ More One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measures the thermal emission of exoplanets. For this, we have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo simulations of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Within a 2.5-year initial search phase, an interferometer consisting of four 2 m apertures with a total instrument throughput of 5% covering a wavelength range between 4 and 18.5 $μ$m could detect up to ~550 exoplanets with radii between 0.5 and 6 R$_\oplus$ with an integrated SNR$\ge$7. At least ~160 of the detected exoplanets have radii $\le$1.5 R$_\oplus$. Depending on the observing scenario, ~25-45 rocky exoplanets (objects with radii between 0.5 and 1.5 $_{\oplus}$) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With an aperture size of 3.5 m, the total number of detections can increase to up to ~770, including ~60-80 rocky, eHZ planets. With 1 m aperture size, the maximum detection yield is ~315 exoplanets, including $\le$20 rocky, eHZ planets. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. (abridged) △ Less

Submitted 20 April, 2022; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: Accepted for publication by A&A - some typos corrected and affiliations updated; 14 pages main text (incl. 14 figures); first paper in the LIFE paper series; papers II (arXiv:2203.00471) and III (arXiv:2112.02054) are also available

Journal ref: A&A 664, A21 (2022)

Rapid CO gas dispersal from NO Lup's class III circumstellar disc

Authors: J. B. Lovell, G. M. Kennedy, S. Marino, M. C. Wyatt, M. Ansdell, M. Kama, C. F. Manara, L. Matrà, G. Rosotti, M. Tazzari, L. Testi, J. P. Williams

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

Submitted 26 November, 2020; originally announced November 2020.

Comments: 6 pages, 4 figures. Accepted for publication in MNRAS following peer review (this is a pre-copyedited, author-produced PDF version)

Formation of Gaps in Self-Gravitating Debris Disks by Secular Resonance in a Single-planet System I: A Simplified Model

Authors: Antranik A. Sefilian, Roman R. Rafikov, Mark C. Wyatt

Abstract: Spatially resolved images of debris disks frequently reveal complex morphologies such as gaps, spirals, and warps. Most existing models for explaining such morphologies focus on the role of massive perturbers (i.e. planets, stellar companions), ignoring the gravitational effects of the disk itself. Here we investigate the secular interaction between an eccentric planet and a massive, external debr… ▽ More Spatially resolved images of debris disks frequently reveal complex morphologies such as gaps, spirals, and warps. Most existing models for explaining such morphologies focus on the role of massive perturbers (i.e. planets, stellar companions), ignoring the gravitational effects of the disk itself. Here we investigate the secular interaction between an eccentric planet and a massive, external debris disk using a simple analytical model. Our framework accounts for both the gravitational coupling between the disk and the planet, as well as the disk self-gravity -- with the limitation that it ignores the non-axisymmetric component of the disk (self-)gravity. We find generally that even when the disk is less massive than the planet, the system may feature secular resonances within the disk (contrary to what may be naively expected), where planetesimal eccentricities get significantly excited. Given this outcome we propose that double-ringed debris disks, such as those around HD 107146 and HD 92945, could be the result of secular resonances with a yet-undetected planet interior to the disk. We characterize the dependence of the properties of the secular resonances (i.e. locations, timescales, and widths) on the planet and disk parameters, finding that the mechanism is robust provided the disk is massive enough. As an example, we apply our results to HD 107146 and find that this mechanism readily produces $\sim 20$ au wide non-axisymmetric gaps. Our results may be used to set constraints on the total mass of double-ringed debris disks. We demonstrate this for HD 206893, for which we infer a disk mass of $\approx 170$ Earth masses by considering perturbations from the known brown dwarf companion. △ Less

Submitted 8 January, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ: 25 Pages, 13 Figures, 1 Table, 3 Animations (see Ancillary files)

ALMA Survey of Lupus Class III Stars: Early Planetesimal Belt Formation and Rapid Disk Dispersal

Authors: J. B. Lovell, M. C. Wyatt, M. Ansdell, M. Kama, G. M. Kennedy, C. F. Manara, S. Marino, L. Matrà, G. Rosotti, M. Tazzari, L. Testi, J. P. Williams

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

Submitted 23 October, 2020; originally announced October 2020.

Comments: 30 pages, 20 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Evolution of the Earth's Atmosphere during Late Veneer Accretion

Authors: Catriona A. Sinclair, Mark C. Wyatt, Alessandro Morbidelli, David Nesvorny

Abstract: Recent advances in our understanding of the dynamical history of the Solar system have altered the inferred bombardment history of the Earth during accretion of the Late Veneer, after the Moon-forming impact. We investigate how the bombardment by planetesimals left-over from the terrestrial planet region after terrestrial planet formation, as well as asteroids and comets, affects the evolution of… ▽ More Recent advances in our understanding of the dynamical history of the Solar system have altered the inferred bombardment history of the Earth during accretion of the Late Veneer, after the Moon-forming impact. We investigate how the bombardment by planetesimals left-over from the terrestrial planet region after terrestrial planet formation, as well as asteroids and comets, affects the evolution of Earth's early atmosphere. We develop a new statistical code of stochastic bombardment for atmosphere evolution, combining prescriptions for atmosphere loss and volatile delivery derived from hydrodynamic simulations and theory with results from dynamical modelling of realistic populations of impactors. We find that for an initially Earth-like atmosphere impacts cause moderate atmospheric erosion with stochastic delivery of large asteroids giving substantial growth ($\times 10$) in a few $\%$ of cases. The exact change in atmosphere mass is inherently stochastic and dependent on the dynamics of the left-over planetesimals. We also consider the dependence on unknowns including the impactor volatile content, finding that the atmosphere is typically completely stripped by especially dry left-over planetesimals ($<0.02 ~ \%$ volatiles). Remarkably, for a wide range of initial atmosphere masses and compositions, the atmosphere converges towards similar final masses and compositions, i.e. initially low mass atmospheres grow whereas massive atmospheres deplete. While the final properties are sensitive to the assumed impactor properties, the resulting atmosphere mass is close to that of current Earth. The exception to this is that a large initial atmosphere cannot be eroded to the current mass unless the atmosphere was initially primordial in composition. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 30 pages, accepted for publication in MNRAS

Mutual inclinations between giant planets and their debris discs in HD 113337 and HD 38529

Authors: Jerry W. Xuan, Grant M. Kennedy, Mark C. Wyatt, Ben Yelverton

Abstract: HD 113337 and HD 38529 host pairs of giant planets, a debris disc, and wide M-type stellar companions. We measure the disc orientation with resolved images from Herschel and constrain the three-dimensional orbits of the outer planets with Gaia DR2 and Hipparcos astrometry. Resolved disc modelling leaves degeneracy in the disc orientation, so we derive four separate planet-disc mutual inclination (… ▽ More HD 113337 and HD 38529 host pairs of giant planets, a debris disc, and wide M-type stellar companions. We measure the disc orientation with resolved images from Herschel and constrain the three-dimensional orbits of the outer planets with Gaia DR2 and Hipparcos astrometry. Resolved disc modelling leaves degeneracy in the disc orientation, so we derive four separate planet-disc mutual inclination ($ΔI$) solutions. The most aligned solutions give $ΔI=17-32°$ for HD 113337 and $ΔI=21-45°$ for HD 38529 (both 1$σ$). In both systems, there is a small probability (<0.3 per cent) that the planet and disc are nearly aligned ($ΔI < 3°$). The stellar and planetary companions cause the orbits of disc material to precess about a plane defined by the forced inclination. We determine this as well as the precession time-scale to interpret the mutual inclination results. We find that the debris discs in both systems could be warped via joint influences of the outer planet and stellar companion, potentially explaining the observed misalignments. However, this requires HD 113337 to be old (0.8-1.7 Gyr), whereas if young (14-21 Myr), the observed misalignment in HD 113337 could be inherited from the protoplanetary disc phase. For both systems, the inclination of the stellar spin axis is consistent with the disc and outer planet inclinations, which instead supports system-wide alignment or near alignment. High-resolution observations of the discs and improved constraints on the planetary orbits would provide firmer conclusions about the (mis)alignment status. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS

Solution to the debris disc mass problem: planetesimals are born small?

Authors: Alexander V. Krivov, Mark C. Wyatt

Abstract: Debris belts on the periphery of planetary systems, encompassing the region occupied by planetary orbits, are massive analogues of the Solar system's Kuiper belt. They are detected by thermal emission of dust released in collisions amongst directly unobservable larger bodies that carry most of the debris disc mass. We estimate the total mass of the discs by extrapolating up the mass of emitting du… ▽ More Debris belts on the periphery of planetary systems, encompassing the region occupied by planetary orbits, are massive analogues of the Solar system's Kuiper belt. They are detected by thermal emission of dust released in collisions amongst directly unobservable larger bodies that carry most of the debris disc mass. We estimate the total mass of the discs by extrapolating up the mass of emitting dust with the help of collisional cascade models. The resulting mass of bright debris discs appears to be unrealistically large, exceeding the mass of solids available in the systems at the preceding protoplanetary stage. We discuss this "mass problem" in detail and investigate possible solutions to it. These include uncertainties in the dust opacity and planetesimal strength, variation of the bulk density with size, steepening of the size distribution by damping processes, the role of the unknown "collisional age" of the discs, and dust production in recent giant impacts. While we cannot rule out the possibility that a combination of these might help, we argue that the easiest solution would be to assume that planetesimals in systems with bright debris discs were "born small", with sizes in the kilometre range, especially at large distances from the stars. This conclusion would necessitate revisions to the existing planetesimal formation models, and may have a range of implications for planet formation. We also discuss potential tests to constrain the largest planetesimal sizes and debris disc masses. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review; 20 pages, 11 figures. The version of record is available online at https://doi.org/10.1093/mnras/staa2385

Evidence for a high mutual inclination between the cold Jupiter and transiting super Earth orbiting $π$ Men

Authors: Jerry W. Xuan, Mark C. Wyatt

Abstract: $π$ Men hosts a transiting super Earth ($P\approx6.27$ d, $m\approx4.82$ $M_{\oplus}$, $R\approx2.04$ $R_{\oplus}$) discovered by TESS and a cold Jupiter ($P\approx2093$ d, $m \sin I\approx10.02$ $M_{\rm{Jup}}$, $e\approx0.64… ▽ More $π$ Men hosts a transiting super Earth ($P\approx6.27$ d, $m\approx4.82$ $M_{\oplus}$, $R\approx2.04$ $R_{\oplus}$) discovered by TESS and a cold Jupiter ($P\approx2093$ d, $m \sin I\approx10.02$ $M_{\rm{Jup}}$, $e\approx0.64$) discovered from radial velocity. We use Gaia DR2 and Hipparcos astrometry to derive the star's velocity caused by the orbiting planets and constrain the cold Jupiter's sky-projected inclination ($I_b=41-65^{\circ}$). From this we derive the mutual inclination ($ΔI$) between the two planets, and find that $49^{\circ}< ΔI < 131^{\circ}$ (1$σ$), and $28^{\circ} < ΔI < 152^{\circ}$ (2$σ$). We examine the dynamics of the system using $N$-body simulations, and find that potentially large oscillations in the super Earth's eccentricity and inclination are suppressed by general relativistic precession. However, nodal precession of the inner orbit around the invariable plane causes the super Earth to only transit between 7-22 per cent of the time, and to usually be observed as misaligned with the stellar spin axis. We repeat our analysis for HAT-P-11, finding a large $ΔI$ between its close-in Neptune and cold Jupiter and similar dynamics. $π$ Men and HAT-P-11 are prime examples of systems where dynamically hot outer planets excite their inner planets, with the effects of increasing planet eccentricities, planet-star misalignments, and potentially reducing the transit multiplicity. Formation of such systems likely involves scattering between multiple giant planets or misaligned protoplanetary discs. Future imaging of the faint debris disc in $π$ Men and precise constraints on the stellar spin orientation would provide strong tests for these formation scenarios. △ Less

Submitted 6 August, 2020; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: Corrected Fig 7 and Sec 6.1.1 (HAT-P-11 also weak planet-star coupling). Published as MNRAS 497, 2096-2118 (2020)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 497, Issue 2, September 2020, Pages 2096-2118

Dust size and spatial distributions in debris discs: predictions for exozodiacal dust dragged in from an exo-Kuiper belt

Authors: Jessica K. Rigley, Mark C. Wyatt

Abstract: The SEDs of some nearby stars show mid-infrared excesses from warm habitable zone dust, known as exozodiacal dust. This dust may originate in collisions in a planetesimal belt before being dragged inwards. This paper presents an analytical model for the size distribution of particles at different radial locations in such a scenario, considering evolution due to destructive collisions and Poynting-… ▽ More The SEDs of some nearby stars show mid-infrared excesses from warm habitable zone dust, known as exozodiacal dust. This dust may originate in collisions in a planetesimal belt before being dragged inwards. This paper presents an analytical model for the size distribution of particles at different radial locations in such a scenario, considering evolution due to destructive collisions and Poynting-Robertson (P-R) drag. Results from more accurate but computationally expensive numerical simulations of this process are used to validate the model and fit its free parameters. The model predicts 11 $μ$m excesses ($R_{11}$) for discs with a range of dust masses and planetesimal belt radii using realistic grain properties. We show that P-R drag should produce exozodiacal dust levels detectable with the Large Binocular Telescope Interferometer (LBTI) ($R_{11} > 0.1\%$) in systems with known outer belts; non-detection may indicate dust depletion, e.g. by an intervening planet. We also find that LBTI could detect exozodiacal dust dragged in from a belt too faint to detect at far-infrared wavelengths, with fractional luminosity $f\sim 10^{-7}$ and radius $\sim 10-80$ au. Application to systems observed with LBTI shows that P-R drag can likely explain most (5/9) of the exozodiacal dust detections in systems with known outer belts; two systems ($β$ Uma and $η$ Corvi) with bright exozodi may be due to exocomets. We suggest that the three systems with exozodiacal dust detections but no known belt may have cold planetesimal belts too faint to be detectable in the far-infrared. Even systems without outer belt detections could have exozodiacal dust levels $R_{11} > 0.04\%$ which are problematic for exo-Earth imaging. △ Less

Submitted 2 July, 2020; originally announced July 2020.

Comments: 25 pages, 20 figures, 2 tables. Accepted for publication in MNRAS

Dust Populations in the Iconic Vega Planetary System Resolved by ALMA

Authors: Luca Matrà, William R. F. Dent, David J. Wilner, Sebastián Marino, Mark C. Wyatt, Jonathan P. Marshall, Kate Y. L. Su, Miguel Chavez, Antonio Hales, A. Meredith Hughes, Jane S. Greaves, Stuartt A. Corder

Abstract: 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 s… ▽ More 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. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 18 pages, 3 figures, Accepted for publication in ApJ

The HOSTS survey for exozodiacal dust: Observational results from the complete survey

Authors: Steve Ertel, Denis Defrère, Philip M. Hinz, Bertrand Mennesson, Grant M. Kennedy, William C. Danchi, Christopher Gelino, John M. Hill, William F. Hoffmann, Johan Mazoyer, George Rieke, Andrew Shannon, Karl Stapelfeldt, Eckhart Spalding, Jordan M. Stone, Amali Vaz, Alycia J. Weinberger, Phil Willems, Olivier Absil, Paul Arbo, Vanessa P. Bailey, Charles Beichman, Geoffrey Bryden, Elwood C. Downey, Olivier Durney , et al. (21 additional authors not shown)

Abstract: The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model p… ▽ More The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates. △ Less

Submitted 6 March, 2020; originally announced March 2020.

Comments: accepted for publication in AJ

Population synthesis of exocometary gas around A stars

Authors: S. Marino, M. Flock, Th. Henning, Q. Kral, L. Matrà, M. C. Wyatt

Abstract: 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 prod… ▽ More 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 ($α\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. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS, 23 pages, 13 figures

Are Inner Disc Misalignments Common? ALMA Reveals an Isotropic Outer Disc Inclination Distribution for Young Dipper Stars

Authors: M. Ansdell, E. Gaidos, C. Hedges, M. Tazzari, A. L. Kraus, M. C. Wyatt, G. M. Kennedy, J. P. Williams, A. W. Mann, I. Angelo, G. Duchene, E. E. Mamajek, J. Carpenter, T. L. Esplin, A. C. Rizzuto

Abstract: Dippers are a common class of young variable star exhibiting day-long dimmings with depths of up to several tens of percent. A standard explanation is that dippers host nearly edge-on (70 deg) protoplanetary discs that allow close-in (< 1 au) dust lifted slightly out of the midplane to partially occult the star. The identification of a face-on dipper disc and growing evidence of inner disc misalig… ▽ More Dippers are a common class of young variable star exhibiting day-long dimmings with depths of up to several tens of percent. A standard explanation is that dippers host nearly edge-on (70 deg) protoplanetary discs that allow close-in (< 1 au) dust lifted slightly out of the midplane to partially occult the star. The identification of a face-on dipper disc and growing evidence of inner disc misalignments brings this scenario into question. Thus we uniformly (re)derive the inclinations of 24 dipper discs resolved with (sub-)mm interferometry from ALMA. We find that dipper disc inclinations are consistent with an isotropic distribution over 0-75 deg, above which the occurrence rate declines (likely an observational selection effect due to optically thick disc midplanes blocking their host stars). These findings indicate that the dipper phenomenon is unrelated to the outer (>10 au) disc resolved by ALMA and that inner disc misalignments may be common during the protoplanetary phase. More than one mechanism may contribute to the dipper phenomenon, including accretion-driven warps and "broken" discs caused by inclined (sub-)stellar or planetary companions. △ Less

Submitted 19 December, 2019; v1 submitted 3 December, 2019; originally announced December 2019.

Comments: Accepted to MNRAS (19 pages, 7 figures, 2 tables)

Modelling the distributions of white dwarf atmospheric pollution: a low Mg abundance for accreted planetesimals?

Authors: Samuel G. D. Turner, Mark C. Wyatt

Abstract: The accretion of planetesimals onto white dwarf atmospheres allows determination of the composition of this polluting material. This composition is usually inferred from observed pollution levels by assuming it originated from a single body. This paper instead uses a stochastic model wherein polluting planetesimals are chosen randomly from a mass distribution, finding that the single body assumpti… ▽ More The accretion of planetesimals onto white dwarf atmospheres allows determination of the composition of this polluting material. This composition is usually inferred from observed pollution levels by assuming it originated from a single body. This paper instead uses a stochastic model wherein polluting planetesimals are chosen randomly from a mass distribution, finding that the single body assumption is invalid in >20% of cases. Planetesimal compositions are modelled assuming parent bodies that differentiated into core, mantle and crust components. Atmospheric levels of Ca, Mg and Fe in the model are compared to a sample of 230 DZ white dwarfs for which such pollution is measured. A good fit is obtained when each planetesimal has its core, mantle and crust fractions chosen independently from logit-normal distributions which lead to average mass fractions of $f_\text{Cru}=0.15$, $f_\text{Man}=0.49$ and $f_\text{Cor}=0.36$. However, achieving this fit requires a factor 4 depletion of Mg relative to stellar material. This depletion is unlikely to originate in planetesimal formation processes, but might occur from heating while the star is on the giant branch. Alternatively the accreted material has stellar abundance, and either the inferred low Mg abundance was caused by an incorrect assumption that Mg sinks slower than Ca and Fe, or there are unmodelled biases in the observed sample. Finally, the model makes predictions for the timescale on which the observed pollutant composition varies, which should be the longer of the sinking and disc timescales, implying variability on decadal timescales for DA white dwarfs. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: 20 pages, 19 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 491, February 2020, Pages 4672-4689

Susceptibility of planetary atmospheres to mass loss and growth by planetesimal impacts: the impact shoreline

Authors: M. C. Wyatt, Q. Kral, C. A. Sinclair

Abstract: This paper considers how planetesimal impacts affect planetary atmospheres. Atmosphere evolution depends on the ratio of gain from volatiles to loss from atmosphere stripping f_v; for constant bombardment, atmospheres with f_v<1 are destroyed in finite time, but grow linearly with time for f_v>1. An impact outcome prescription is used to characterise how f_v depends on planetesimal impact velociti… ▽ More This paper considers how planetesimal impacts affect planetary atmospheres. Atmosphere evolution depends on the ratio of gain from volatiles to loss from atmosphere stripping f_v; for constant bombardment, atmospheres with f_v<1 are destroyed in finite time, but grow linearly with time for f_v>1. An impact outcome prescription is used to characterise how f_v depends on planetesimal impact velocities, size distribution and composition. Planets that are low mass and/or close to the star have atmospheres that deplete in impacts, while high mass and/or distant planets grow secondary atmospheres. Dividing these outcomes is an fv=1 impact shoreline analogous to Zahnle & Catling's cosmic shoreline. The impact shoreline's location depends on assumed impacting planetesimal properties, so conclusions for the atmospheric evolution of a planet like Earth with f_v~1 are only as strong as those assumptions. Application to the exoplanet population shows the gap in the planet radius distribution at ~1.5R_earth is coincident with the impact shoreline, which has a similar dependence on orbital period and stellar mass to the observed gap. Given sufficient bombardment, planets below the gap would be expected to lose their atmospheres, while those above could have atmospheres enhanced in volatiles. The level of atmosphere alteration depends on the total bombardment a planet experiences, and so on the system's (usually unknown) other planets and planetesimals, though massive distant planets would have low accretion efficiency. Habitable zone planets around lower luminosity stars are more susceptible to atmosphere stripping, disfavouring M stars as hosts of life-bearing planets if Earth-like bombardment is conducive to the development of life. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted for publication by MNRAS

Extrasolar Kuiper Belts

Authors: Mark C. Wyatt

Abstract: Extrasolar debris disks are the dust disks found around nearby main sequence stars arising from the break-up of asteroids and comets orbiting the stars. Far-IR surveys (e.g., with Herschel) showed that ~20% of stars host detectable dust levels. While dust temperatures suggest a location at 10s of au comparable with our Kuiper belt, orders of magnitude more dust is required implying a planetesimal… ▽ More Extrasolar debris disks are the dust disks found around nearby main sequence stars arising from the break-up of asteroids and comets orbiting the stars. Far-IR surveys (e.g., with Herschel) showed that ~20% of stars host detectable dust levels. While dust temperatures suggest a location at 10s of au comparable with our Kuiper belt, orders of magnitude more dust is required implying a planetesimal population more comparable with the primordial Kuiper belt. High resolution imaging (e.g., with ALMA) has mapped the nearest and brightest disks, providing evidence for structures shaped by an underlying planetary system. Some of these are analogous to structures in our own Kuiper belt (e.g., the hot and cold classical, resonant, scattered disk and cometary populations), while others have no Solar System counterpart. CO gas is seen in some debris disks, and inferred to originate in the destruction of planetesimals with a similar volatile-rich composition to Solar System comets. This chapter reviews our understanding of extrasolar Kuiper belts and of how our own Kuiper belt compares with those of neighbouring stars. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: Chapter to be published in the book "The Trans-Neptunian Solar System" (D. Prialnik, M.A. Barucci and L. Young, eds.). 25 pages, 8 figures

Searching for a cometary belt around Trappist-1 with ALMA

Authors: S. Marino, M. C. Wyatt, G. M. Kennedy, M. Kama, L. Matrà, A. H. M. J. Triaud, Th. Henning

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

Submitted 27 January, 2020; v1 submitted 19 September, 2019; originally announced September 2019.

Comments: accepted for publication in MNRAS, 7 pages, 5 figures

A statistically significant lack of debris discs in medium separation binary systems

Authors: Ben Yelverton, Grant M. Kennedy, Kate Y. L. Su, Mark C. Wyatt

Abstract: We compile a sample of 341 binary and multiple star systems with the aim of searching for and characterising Kuiper belt-like debris discs. The sample is assembled by combining several smaller samples studied in previously published work with targets from two unpublished Herschel surveys. We find that 38 systems show excess emission at 70 or 100 $μ$m suggestive of a debris disc. While nine of the… ▽ More We compile a sample of 341 binary and multiple star systems with the aim of searching for and characterising Kuiper belt-like debris discs. The sample is assembled by combining several smaller samples studied in previously published work with targets from two unpublished Herschel surveys. We find that 38 systems show excess emission at 70 or 100 $μ$m suggestive of a debris disc. While nine of the discs appear to be unstable to perturbations from their host binary based on a simple analysis of their inferred radii, we argue that the evidence for genuine instability is not strong, primarily because of uncertainty in the true disc radii, uncertainty in the boundaries of the unstable regions, and orbital projection effects. The binary separation distributions of the disc-bearing and disc-free systems are different at a confidence level of $99.4\%$, indicating that binary separation strongly influences the presence of detectable levels of debris. No discs are detected for separations between $\sim$25 and 135 au; this is likely a result of binaries whose separations are comparable with typical disc radii clearing out their primordial circumstellar or circumbinary material via dynamical perturbations. The disc detection rate is $19^{+5}_{-3}\%$ for binaries wider than 135 au, similar to published results for single stars. Only $8^{+2}_{-1}\%$ of systems with separations below 25 au host a detectable disc, which may suggest that planetesimal formation is inhibited in binaries closer than a few tens of au, similar to the conclusions of studies of known planet-hosting binaries. △ Less

Submitted 11 July, 2019; v1 submitted 10 July, 2019; originally announced July 2019.

Comments: MNRAS, in press; updated publication year

The REASONS Survey: Resolved Millimeter Observations of a Large Debris Disk Around the Nearby F Star HD 170773

Authors: Aldo G. Sepulveda, Luca Matra, Grant M. Kennedy, Carlos del Burgo, Karin I. Oberg, David J. Wilner, Sebastian Marino, Mark Booth, John M. Carpenter, Claire L. Davies, William R. F. Dent, Steve Ertel, Jean-Francois Lestrade, Jonathan P. Marshall, Julien Milli, Mark C. Wyatt, Meredith A. MacGregor, Brenda C. Matthews

Abstract: 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… ▽ More 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. △ Less

Submitted 20 June, 2019; originally announced June 2019.

Comments: 14 pages, 6 figures, accepted to ApJ

Journal ref: The Astrophysical Journal, 881, 84 (2019)

Oort cloud asteroids: Collisional evolution, the Nice Model, and the Grand Tack

Authors: Andrew Shannon, Alan P. Jackson, Mark C. Wyatt

Abstract: If the Solar system had a history of planet migration, the signature of that migration may be imprinted on the populations of asteroids and comets that were scattered in the planets' wake. Here, we consider the dynamical and collisional evolution of inner Solar system asteroids which join the Oort cloud. We compare the Oort cloud asteroid populations produced by migration scenarios based on the `N… ▽ More If the Solar system had a history of planet migration, the signature of that migration may be imprinted on the populations of asteroids and comets that were scattered in the planets' wake. Here, we consider the dynamical and collisional evolution of inner Solar system asteroids which join the Oort cloud. We compare the Oort cloud asteroid populations produced by migration scenarios based on the `Nice' and `Grand Tack' scenarios, as well as a null hypothesis where the planets have not migrated, to the detection of one such object, C/2014 S3 (PANSTARRS). Our simulations find that the discovery of C/2014 S3 (PANSTARRS) only has a greater than one percent chance of occurring if the Oort cloud asteroids evolved on to Oort cloud orbits when the Solar system was not more than about one million years old, as this early transfer to the Oort cloud is necessary to keep the amount of collisional evolution low. We argue this only occurs when a giant (greater than thirty Earth masses) planet orbits at 1 ~ 2 au, and thus our results strongly favour a `Grand Tack'-like migration having occurred early in the Solar system's history. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: Resubmitted to MNRAS after revision from referee's report