No Arabic 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 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 ($alphasim0.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.
Previous work concerning planet formation around low-mass stars has often been limited to large planets and individual systems. As current surveys routinely detect planets down to terrestrial size in these systems, a more holistic approach that reflects their diverse architectures is timely. Here, we investigate planet formation around low-mass stars and identify differences in the statistical distribution of planets. We compare the synthetic planet populations to observed exoplanets. We used the Generation III Bern model of planet formation and evolution to calculate synthetic populations varying the central star from solar-like stars to ultra-late M dwarfs. This model includes planetary migration, N-body interactions between embryos, accretion of planetesimals and gas, and long-term contraction and loss of the gaseous atmospheres. We find that temperate, Earth-sized planets are most frequent around early M dwarfs and more rare for solar-type stars and late M dwarfs. The planetary mass distribution does not linearly scale with the disk mass. The reason is the emergence of giant planets for M*>0.5 Msol, which leads to the ejection of smaller planets. For M*>0.3 Msol there is sufficient mass in the majority of systems to form Earth-like planets, leading to a similar amount of Exo-Earths going from M to G dwarfs. In contrast, the number of super-Earths and larger planets increases monotonically with stellar mass. We further identify a regime of disk parameters that reproduces observed M-dwarf systems such as TRAPPIST-1. However, giant planets around late M dwarfs such as GJ 3512b only form when type I migration is substantially reduced. We quantify the stellar mass dependence of multi-planet systems using global simulations of planet formation and evolution. The results compare well to current observational data and predicts trends that can be tested with future observations.
Millimeter observations of CO gas in planetesimal belts show a high detection rate around A stars, but few detections for later type stars. We present the first CO detection in a planetesimal belt around an M star, TWA 7. The optically thin CO (J=3-2) emission is co-located with previously identified dust emission from the belt, and the emission velocity structure is consistent with Keplerian rotation around the central star. The detected CO is not well shielded against photodissociation, and must thus be continuously replenished by gas release from exocomets within the belt. We analyze in detail the process of exocometary gas release and destruction around young M dwarfs and how this process compares to earlier type stars. Taking these differences into account, we find that CO generation through exocometary gas release naturally explains the increasing CO detection rates with stellar luminosity, mostly because the CO production rate from the collisional cascade is directly proportional to stellar luminosity. More luminous stars will therefore on average host more massive (and hence more easily detectable) exocometary CO disks, leading to the higher detection rates observed. The current CO detection rates are consistent with a ubiquitous release of exocometary gas in planetesimal belts, independent of spectral type.
Evidence for exocomets, icy bodies in extrasolar planetary systems, has rapidly increased over the past decade. Volatiles are detected through the gas that exocomets release as they collide and grind down within their natal belts, or as they sublimate once scattered inwards to the regions closest to their host star. Most detections are in young, 10 to a few 100 Myr-old systems that are undergoing the final stages of terrestrial planet formation. This opens the exciting possibility to study exocomets at the epoch of volatile delivery to the inner regions of planetary systems. Detection of molecular and atomic gas in exocometary belts allows us to estimate molecular ice abundances and overall elemental abundances, enabling comparison with the Solar Nebula and Solar System comets. At the same time, observing star-grazing exocomets transiting in front of their star (for planetary systems viewed edge-on) and exozodiacal dust in the systems innermost regions gives unique dynamical insights into the inward scattering process producing delivery to inner rocky planets. The rapid advances of this budding subfield of exoplanetary science will continue in the short term with the upcoming JWST, WFIRST and PLATO missions. In the longer term, the priority should be to explore the full composition of exocomets, including species crucial for delivery and later prebiotic synthesis. Doing so around an increasingly large population of exoplanetary systems is equally important, to enable comparative studies of young exocomets at the epoch of volatile delivery. We identify the proposed LUVOIR and Origins flagship missions as the most promising for a large-scale exploration of exocometary gas, a crucial component of the chemical heritage of young exo-Earths.
The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable $r$-modes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains and unstable $r$-modes. However, based on the resultant distributions alone it is difficult to distinguish between the competing mechanisms.
Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.