A quarter of DA white dwarfs are metal polluted, yet elements heavier than helium sink down through the stellar atmosphere on time-scales of days. Hence, these white dwarfs must be the currently ...accreting material containing heavy elements. Here we consider whether the scattering of comets or asteroids from an outer planetary system, following stellar mass-loss on the asymptotic giant branch, can reproduce these observations. We use N-body simulations to investigate the effects of stellar mass-loss on a simple system consisting of a planetesimal belt whose inner edge is truncated by a planet. Our simulations find that, starting with a planetesimal belt population fitted to the observed main-sequence evolution, sufficient mass is scattered into the inner planetary system to explain the inferred heavy element accretion rates. This assumes that a fraction of the mass scattered into the inner planetary system ends up on star-grazing orbits, is tidally disrupted and is accreted on to the white dwarf. The simulations also reproduce the observed decrease in accretion rate with cooling age and predict accretion rates in old (>1 Gyr) white dwarfs, in line with observations. The efficiency we assumed for material scattered into the inner planetary system to end up on star-grazing orbits is based on a solar-like planetary system, since the simulations show that a single planet is not sufficient. Although the correct level of accretion is reproduced, the simulations predict a higher fraction of accreting white dwarfs than observed. This could indicate that the evolved planetary systems are less efficient in scattering bodies on to star-grazing orbits or that dynamical instabilities post-stellar mass-loss cause rapid planetesimal belt depletion for a significant fraction of systems.
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, multiwavelength 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 time-scales, 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 colour 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 disc orbiting WD 1054–226 may be typical of polluted white dwarfs, and detected only via favourable 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.
ABSTRACT
Planets orbiting intermediate- and low-mass stars are in jeopardy as their stellar hosts evolve to white dwarfs (WDs) because the dynamics of the planetary system changes due to the increase ...of the planet:star mass ratio after stellar mass-loss. In order to understand how the planet multiplicity affects the dynamical stability of post-main sequence (MS) systems, we perform thousands of N-body simulations involving planetary multiplicity as the variable and with a controlled physical and orbital parameter space:equal-mass planets; the same orbital spacing between adjacent planet’s pairs; and orbits with small eccentricities and inclinations. We evolve the host star from the MS to the WD phase following the system dynamics for 10 Gyr. We find that the fraction of dynamically active simulations on the WD phase for two-planet systems is $10.2^{+1.2}_{-1.0}$–$25.2^{+2.5}_{-2.2}$ ${{\rm per\,cent}}$ and increases to $33.6^{+2.3}_{-2.2}$–$74.1^{+3.7}_{-4.6}$ ${{\rm per\,cent }}$ for the six-planet systems, where the ranges cover different ranges of initial orbital separations. Our simulations show that the more planets the system has, the more systems become unstable when the star becomes a WD, regardless of the planet masses and range of separations. Additional results evince that simulations with low-mass planets (1, 10 M⊕) lose at most two planets, have a large fraction of systems undergoing orbit crossing without planet losses, and are dynamically active for Gyr time-scales on the WD’s cooling track. On the other hand, systems with high-mass planets (100, 1000 M⊕) lose up to five planets, preferably by ejections, and become unstable in the first few hundred Myr after the formation of the WD.
Pebble accretion is an efficient mechanism that is able to build up the core of the giant planets within the lifetime of the protoplanetary disc gas-phase. The core grows via this process until the ...protoplanet reaches its pebble isolation mass and starts to accrete gas. During the growth, the protoplanet undergoes a rapid, large-scale, inward migration due to the interactions with the gaseous protoplanetary disc. In this work, we have investigated how this early migration would have affected the minor body populations in our solar system. In particular, we focus on the Jupiter Trojan asteroids (bodies in the coorbital resonance 1:1 with Jupiter, librating around the L4 and L5 Lagrangian points called, respectively, the leading and the trailing swarm) and the Hilda asteroids. We characterised their orbital parameter distributions after the disc dispersal and their formation location and compare them to the same populations produced in a classical in situ growth model. We find that a massive and eccentric Hilda group is captured during the migration from a region between 5 and 8 au and subsequently depleted during the late instability of the giant planets. Our simulations also show that inward migration of the giant planets always produces a Jupiter Trojans’ leading swarm more populated than the trailing one, with a ratio comparable to the current observed Trojan asymmetry ratio. The in situ formation of Jupiter, on the other hand, produces symmetric swarms. The reason for the asymmetry is the relative drift between the migrating planet and the particles in the coorbital resonance. The capture happens during the growth of Jupiter’s core and Trojan asteroids are afterwards carried along during the giant planet’s migration to their final orbits. The asymmetry and eccentricity of the captured Trojans correspond well to observations, but their inclinations are near zero and their total mass is three to four orders of magnitude higher than the current population. Future modelling will be needed to understand whether the dynamical evolution of the Trojans over billions of years will raise the inclinations and deplete the masses to observed values.
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 spectra with high signal-to-noise ratio obtained with Keck/HIRES. HD 80606 is known to host a giant planet with the mass of four Jupiters, but no planet has been detected around HD 80607 so far. We determined stellar parameters as well as abundances for 23 elements for these two stars with extremely high precision. Our main results are that (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 ± 0.002 dex (σ = 0.009 dex); and (ii) there is no obvious trend between abundance differences and condensation temperature. Assuming that 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 cannot exclude the possibility that HD 80606 might have accreted about 2.5 to 5 MEarth material onto its surface, possibly from a planet destabilised by the known highly eccentric giant.
ABSTRACT
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${{\ \rm per\ cent}}$ 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 effects from binary companions or their survival through the common envelope phase.
ABSTRACT
Between 25 and 50 ${{\ \rm per\ cent}}$ of white dwarfs (WD) present atmospheric pollution by metals, mainly by rocky material, which has been detected as gas/dust discs, or in the form of ...photometric transits in some WDs. Planets might be responsible for scattering minor bodies that can reach stargazing orbits, where the tidal forces of the WD can disrupt them and enhance the chances of debris to fall on to the WD surface. The planet–planet scattering process can be triggered by the stellar mass-loss during the post main-sequence (MS) evolution of planetary systems. In this work, we continue the exploration of the dynamical instabilities that can lead to WD pollution. In a previous work, we explored two-planet systems found around MS stars and here we extend the study to three-planet system architectures. We evolved 135 detected three-planet systems orbiting MS stars to the WD phase by scaling their orbital architectures in a way that their dynamical properties are preserved using the N-body integrator package mercury. We find that 100 simulations (8.6 ${{\ \rm per\ cent}}$) are dynamically active (having planet losses, orbit crossing, and scattering) on the WD phase, where low-mass planets (1–100 M⊕) tend to have instabilities in Gyr time-scales, while high-mass planets (>100 M⊕) decrease the dynamical events more rapidly as the WD ages. Besides, 19 simulations (1.6 ${{\ \rm per\ cent}}$) were found to have planets crossing the Roche radius of the WD, where 9 of them had planet–star collisions. Our three-planet simulations have a slight increase in percentage of simulations that may contribute to the WD pollution than the previous study involving two-planet systems and have shown that planet–planet scattering is responsible of sending planets close to the WD, where they may collide directly to the WD, become tidally disrupted or circularize their orbits, hence producing pollution on the WD atmosphere.
We study how close-in systems such as those detected by Kepler are affected by the dynamics of bodies in the outer system. We consider two scenarios: outer systems of giant planets potentially ...unstable to planet-planet scattering and wide binaries that may be capable of driving Kozai or other secular variations of outer planets' eccentricities. Dynamical excitation of planets in the outer system reduces the multiplicity of Kepler-detectable planets in the inner system in ~20-25 per cent of our systems. Accounting for the occurrence rates of wide-orbit planets and binary stars, ≈18 per cent of close-in systems could be destabilized by their outer companions in this way. This provides some contribution to the apparent excess of systems with a single transiting planet compared to multiple; however, it only contributes at most 25 per cent of the excess. The effects of the outer dynamics can generate systems similar to Kepler-56 (two coplanar planets significantly misaligned with the host star) and Kepler-108 (two significantly non-coplanar planets in a binary). We also identify three pathways to the formation of eccentric warm Jupiters resulting from the interaction between outer and inner systems: direct inelastic collision between an eccentric outer and an inner planet; secular eccentricity oscillations that may 'freeze out' when scattering resolves in the outer system; and scattering in the inner system followed by 'uplift', where inner planets are removed by interaction with the outer planets. In these scenarios, the formation of eccentric warm Jupiters is a signature of a past history of violent dynamics among massive planets beyond ~1 au.
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 retain 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.45 g 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 and 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. This work shows that our framework is capable of constraining important properties of planet formation models.
ABSTRACT
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 on to the star. This can occur through a planet–planet scattering process triggered by the stellar mass-loss during the post main-sequence (MS) 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 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 ≈2.3 ${{\ \rm per\ cent}}$ 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 five planetary configurations ending as a collision of a planet with the WD. Finally 3.2 ${{\ \rm per\ cent}}$ 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.