Dynamical evolution of two-planet systems and its connection with white dwarf atmospheric pollution

Asteroid material is detected in white dwarfs (WDs) as atmospheric pollution by metals, in the form of gas/dust discs, or in photometric transits. Within the current paradigm, minor bodies need to be scattered, most likely by planets, into highly eccentric orbits where the material gets disrupted by tidal forces and then accreted onto the star. This can occur through a planet-planet scattering process triggered by the stellar mass loss during the post main-sequence evolution of planetary systems. So far, studies of the $N$-body dynamics of this process have used artificial planetary system architectures built ad hoc. In this work, we attempt to go a step further and study the dynamical instability provided by more restrictive systems, that, at the same time allow us an exploration of a wider parameter space: the hundreds of multiple planetary systems found around main-sequence (MS) stars. We find that most of our simulated systems remain stable during the MS, Red and Asymptotic Giant Branch and for several Gyr into the WD phases of the host star. Overall, only $\approx$ 2.3$\%$ of the simulated systems lose a planet on the WD as a result of dynamical instability. If the instabilities take place during the WD phase most of them result in planet ejections with just 5 planetary configurations ending as a collision of a planet with the WD. Finally 3.2$\%$ of the simulated systems experience some form of orbital scattering or orbit crossing that could contribute to the pollution at a sustained rate if planetesimals are present in the same system.Comment: 18 pages, 14 figure

Do instabilities in high-multiplicity systems explain the existence of close-in white dwarf planets?

We investigate the origin of close-in planets and related phenomena orbiting white dwarfs (WDs), which are thought to originate from orbits more distant from the star. We use the planetary architectures of the 75 multiple-planet systems (four, five and six planets) detected orbiting main-sequence stars to build 750 dynamically analogous templates that we evolve to the WD phase. Our exploration of parameter space, although not exhaustive, is guided and restricted by observations and we find that the higher the multiplicity of the planetary system, the more likely it is to have a dynamical instability (losing planets, orbit crossing and scattering), that eventually will send a planet (or small object) through a close periastron passage. Indeed, the fraction of unstable four- to six-planet simulations is comparable to the 25-50$\%$ fraction of WDs having atmospheric pollution. Additionally, the onset of instability in the four- to six-planet configurations peaks in the first Gyr of the WD cooling time, decreasing thereafter. Planetary multiplicity is a natural condition to explain the presence of close-in planets to WDs, without having to invoke the specific architectures of the system or their migration through the von Zeipel-Lidov-Kozai (ZLK) effects from binary companions or their survival through the common envelope phase.Comment: 5 pages, 4 figures, 1 table, accepted to MNRAS Letter

Capture of satellites during planetary encounters A case study of the Neptunian moons Triton and Nereid

Single-binary scattering may lead to an exchange where the single object captures a component of the binary, forming a new binary. This has been well studied in encounters between a star--planet pair and a single star. Here we explore the application of the exchange mechanism to a planet--satellite pair and another planet in the gravitational potential of a central star. As a case study, we focus on encounters between a satellite-bearing object and Neptune. We investigate whether Neptune can capture satellites from that object and if the captured satellites have orbits analogous to the Neptunian moons Triton and Nereid. Using $N$-body simulations, we study the capture probability at different encounter distances. Post-capture, we use a simple analytical argument to estimate how the captured orbits evolve under collisional and tidal effects. We find that the average capture probability reaches $\sim$$10\%$ if Neptune penetrates the donor planet's satellite system. Most moons grabbed by Neptune acquire highly eccentric orbits. Post-capture, around half of those captured, especially those on tight orbits, can be circularised, either by tides only or by collisions+tides, turning into Triton-like objects. Captures further out, on the other hand, stay on wide and eccentric orbits like that of Nereid. Both moon types can be captured in the same encounter and they have wide distributions in orbital inclination. Therefore, Triton naturally has a $\sim$50\% chance of being retrograde. A similar process potentially applies to an exoplanetary system, and our model predicts that exomoons can jump from one planet to another during planetary scattering. Specifically, there should be two distinct populations of captured moons: one on close-in circular orbits and the other on far-out eccentric orbits. The two populations may have highly inclined prograde or retrograde orbits.Comment: A&A in pres

Detailed chemical compositions of the wide binary HD 80606/80607: revised stellar properties and constraints on planet formation

Differences in the elemental abundances of planet hosting stars in binary systems can give important clues and constraints about planet formation and evolution. In this study we performed a high-precision, differential elemental abundance analysis of a wide binary system, HD 80606/80607, based on high-resolution, high signal-to-noise ratio Keck/HIRES spectra. HD 80606 is known to host a four Jupiter mass giant planet while no planet has yet been detected around HD 80607. We determined stellar parameters as well as abundances for 23 elements for these two stars with extremely high precision. Our main results are: (i) we confirmed that the two components share very similar chemical compositions, but HD 80606 is marginally more metal-rich than HD 80607 with an average difference of +0.013 $\pm$ 0.002 dex ($\sigma$ = 0.009 dex) and (ii) there is no obvious trend between abundance differences and condensation temperature. Assuming this binary formed from material with the same chemical composition, it is difficult to understand how giant planet formation could produce the present-day photospheric abundances of the elements we measure. We can not exclude the possibility that HD 80606 might have accreted about 2.5 to 5 $M_{\rm Earth}$ material onto its surface, possibly from a planet destabilised by the known highly-eccentric giant.Comment: 11 pages, 9 figues, accepted for publication in A&

Exocomet signatures around the A-shell star Φ\Phi Leo?

We present an intensive monitoring of high-resolution spectra of the Ca {\sc ii} K line in the A7IV shell star $\Phi$ Leo at very short (minutes, hours), short (night to night), and medium (weeks, months) timescales. The spectra show remarkable variable absorptions on timescales of hours, days, and months. The characteristics of these sporadic events are very similar to most that are observed toward the debris disk host star $\beta$ Pic, which are commonly interpreted as signs of the evaporation of solid, comet-like bodies grazing or falling onto the star. Therefore, our results suggest the presence of solid bodies around $\Phi$ Leo. To our knowledge, with the exception of $\beta$ Pic, our monitoring has the best time resolution at the mentioned timescales for a star with events attributed to exocomets. Assuming the cometary scenario and considering the timescales of our monitoring, our results indicate that $\Phi$ Leo presents the richest environment with comet-like events known to date, second only to $\beta$ Pic.Comment: A&A letters, proof-correcte

ESPRESSO Mass determination of TOI-263b: An extreme inhabitant of the brown dwarf desert

The TESS mission has reported a wealth of new planetary systems around bright and nearby stars amenable for detailed characterization of the planet properties and their atmospheres. However, not all interesting TESS planets orbit around bright host stars. TOI-263b is a validated ultra-short period substellar object in a 0.56-day orbit around a faint (V=18.97) M3.5 dwarf star. The substellar nature of TOI-263b was explored using multi-color photometry, which determined a true radius of 0.87+-0.21 Rj, establishing TOI-263b's nature ranging from an inflated Neptune to a brown dwarf. The orbital period-radius parameter space occupied by TOI-263b is quite unique, which prompted a further characterization of its true nature. Here, we report radial velocity measurements of TOI-263 obtained with 3 VLT units and the ESPRESSO spectrograph to retrieve the mass of TOI-263b. We find that TOI-263b is a brown dwarf with a mass of 61.6+-4.0 Mj. Additionally, the orbital period of the brown dwarf is found to be synchronized with the rotation period of the host star, and the system is found to be relatively active, possibly revealing a star--brown dwarf interaction. All these findings suggest that the system's formation history might be explained via disc fragmentation and later migration to close-in orbits. If the system is found to be unstable, TOI-263 is an excellent target to test the migration mechanisms before the brown dwarf becomes engulfed by its parent star.Comment: Accepted for Publication in Astronomy and Astrophysic

The Kepler-11 system: evolution of the stellar high-energy emission and {initial planetary} atmospheric mass fractions

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.Comment: 8 pages, 3 figure

Close-in planets around giant stars: lack of hot-Jupiters and prevalence of multiplanetary systems

Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that 70.0 ± 6.6% of the planets in this regime are inner components of multiplanetary systems. This value is 4.2σ higher than for main-sequence hosts, which we find to be 42.4 ± 0.1%. The properties of the known planets seem to indicate that the closest-in planets (a< 0.06 au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits 0.06 <a< 0.5 au. We discuss the implications of this emerging observational trend in the context of formation and evolution of hot Jupiters

Debris Disks: Probing Planet Formation

Debris disks are the dust disks found around ~20% of nearby main sequence stars in far-IR surveys. They can be considered as descendants of protoplanetary disks or components of planetary systems, providing valuable information on circumstellar disk evolution and the outcome of planet formation. The debris disk population can be explained by the steady collisional erosion of planetesimal belts; population models constrain where (10-100au) and in what quantity (>1Mearth) planetesimals (>10km in size) typically form in protoplanetary disks. Gas is now seen long into the debris disk phase. Some of this is secondary implying planetesimals have a Solar System comet-like composition, but some systems may retain primordial gas. Ongoing planet formation processes are invoked for some debris disks, such as the continued growth of dwarf planets in an unstirred disk, or the growth of terrestrial planets through giant impacts. Planets imprint structure on debris disks in many ways; images of gaps, clumps, warps, eccentricities and other disk asymmetries, are readily explained by planets at >>5au. Hot dust in the region planets are commonly found (<5au) is seen for a growing number of stars. This dust usually originates in an outer belt (e.g., from exocomets), although an asteroid belt or recent collision is sometimes inferred.Comment: Invited review, accepted for publication in the 'Handbook of Exoplanets', eds. H.J. Deeg and J.A. Belmonte, Springer (2018