Catriona H. McDonald1,2, Dimitri Veras1,2
1 Centre for Exoplanets and Habitability, University of Warwick, Coventry, CV4 7AL, UK
2 Department of Physics, University of Warwick, Coventry, CV4 7AL, UK
White dwarfs provide a unique opportunity to probe exoplanetary chemistry and post-main-sequence planetary system architectures. Between a quarter and half of all white dwarfs observed show evidence of heavy elements closely aligned with terrestrial planet composition in their atmospheres (Koester et al., 2014, Bonsor et al., 2020). The extreme densities of white dwarfs cause these elements to gravitationally settle on fast timescales and disappear from view. So, the material we see must have been recently accreted from a disrupting planetesimal.
Direct observations of two disrupting exo-asteroids have bolstered this idea (Vanderburg et al., 2015, Vanderbosch et al., 2019). The processes which lead to planetesimals breaking up around white dwarfs are still poorly understood and this work aims to further our understanding of this area.
Previous theoretical work has used spherical approximations. However, observations of Solar System asteroids show that they come in a variety of different shapes and sizes (as can be seen in Figure 1) but generally they can be well described by an ellipsoidal model.
Figure 2: Interactive plots showing the six different ellipsoidal shape models used in this work.
Here we move away from the typical spherical shape model that has been used in previous theoretical work and adopt an ellipsoidal shape.
Ellipsoids can be described by two aspect ratios b and c, which relate the shape's different semi-axes (a, b and c in size order) and give an idea of the degree of elongation of the shape.
b = b/a
c = c/a
For example, a typical asteroid could have b = 0.8, whereas the interstellar asteroid 'Oumuamua may have b = 0.1 (Meech 2017)!
We define an ellipsoid as oblate when b=1 and prolate when b=c. Interactive examples of all of the shape models used in this work can be seen to the left. The different models can be selected using the drop down menu.
Following a procedure outlined in Brown et al., 2017 for quasi-spherical planetesimals, we developed an analytical model which considers ellipsoidal asteroids approaching a white dwarf on extremely eccentric (e ∼ 1) orbits. The possible disruption is divided into three outcomes:
The process comes down to identifying two variables for any combination of white dwarf and asteroid properties.
Once these two parameters have been found, they can be run through the following to identify its ultimate destruction.
Figure 3: Logical flowchart depicting how to identify how an asteroid will disrupt in our model.
An asteroid belt similar to the Solar System's main belt could survive a star's giant branch evolutionary phase and provide an ample reservoir to pollute white dwarfs. Here we construct a simplified main belt to investigate the expected debris distribution if every asteroid in the belt was randomly perturbed towards the white dwarf. The main belt is constructed with the following properties:
We find that utilising the triaxial model makes a significant difference to the results produced using a spherical model. The distance from the white dwarf where destruction occurs can change by orders of magnitude, and in some cases the destruction mode itself can change.
For five different white dwarf temperatures, and hence cooling ages, the main belt analogue is put through the analytical framework and the distance from the white dwarf and mode of disruption are plotted.
Figure 4: results legend, showing the parameters included. Colour indicates the shape model used, the marker shape indicates how the asteroid will disrupt and the marker fill indicates the material the asteroid is made of. Dashed lines indicate the shape model is extreme.
Figure 5: these figures show how an asteroid's semi-major axis changes as it approaches a white dwarf. The legend is as indicated in Figure 4. The sizes of some Solar System asteroids are shown on the right hand axis.
We find that all three modes of destruction occur at every cooling age. As the white dwarf cools and its radiative power decreases, sublimation plays less of a role and direct impacts are more likely. Fragmentation occurs at all times when the asteroid is above the limiting size for its material.
For rocky asteroids, different shape models do not affect the disruption outcome. Whereas snowy asteroids with extreme shape models are more susceptible to sublimation at larger sizes than standard shape models.
Triaxial shape models provide a more accurate model for asteroids than a spherical model. We can place constraints on where the planetary debris from a disrupted asteroid will orbit. This simple analytical formulism can be used in other studies to predict the destruction regime of specific asteroids.
Further work on this project intends to include the effect of rotation and to follow the evolution of the fragmentation products.
We invite feedback and opportunities for collaboration in furthering this work.