Ultra-massive white dwarfs are powerful tools to study various physical processes in the Asymptotic Giant Branch (AGB), type Ia supernova explosions and the theory of crystallization through white dwarf asteroseismology. Despite the interest in these white dwarfs, there are few evolutionary studies in the literature devoted to them. Here, we present new ultra-massive white dwarf evolutionary sequences that constitute an improvement over previous ones. In these new sequences, we take into account for the first time the process of phase separation expected during the crystallization stage of these white dwarfs, by relying on the most up-to-date phase diagram of dense oxygen/neon mixtures. Realistic chemical profiles resulting from the full computation of progenitor evolution during the semidegenerate carbon burning along the super-AGB phase are also considered in our sequences. Outer boundary conditions for our evolving models are provided by detailed non-gray white dwarf model atmospheres for hydrogen and helium composition. We assessed the impact of all these improvements on the evolutionary properties of ultra-massive white dwarfs, providing up-dated evolutionary sequences for these stars. We conclude that crystallization is expected to affect the majority of the massive white dwarfs observed with effective temperatures below 40 000 K. Moreover, the calculation of the phase separation process induced by crystallization is necessary to accurately determine the cooling age and the mass-radius relation of massive white dwarfs. We also provide colors in the GAIA photometric bands for our H-rich white dwarf evolutionary sequences on the basis of new models atmospheres. Finally, these new white dwarf sequences provide a new theoretical frame to perform asteroseismological studies on the recently detected ultra-massive pulsating white dwarfs.

We created a triaxial stellar system through the cold dissipationless collapse of 100,000 particles whose evolution was followed with a multipolar code. Once an equilibrium system had been obtained, the multipolar expansion was freezed and smoothed in order to get a stationary smooth potential. The resulting model was self-consistent and the orbits and Lyapunov exponents could then be computed for a randomly selected sample of 3472 of the bodies that make up the system. More than half of the orbits (52.7%) turned out to be chaotic. Regular orbits were then classified using the frequency analysis automatic code of Carpintero and Aguilar (1998, MNRAS 298(1), 1-21). We present plots of the distributions of the different kinds of orbits projected on the symmetry planes of the system. We distinguish chaotic orbits with only one non-zero Lyapunov exponent from those with two non-zero exponents and show that their spatial distributions differ, that of the former being more similar to the one of the regular orbits. Most of the regular orbits are boxes and boxlets, but the minor axis tubes play an important role filling in the wasp waists of the boxes and helping to give a lentil shape to the system. We see no problem in building stable triaxial models with substantial amounts of chaotic orbits; the difficulties found by other authors may be due not to a physical cause but to a limitation of Schwarzschild's method.

Context. Recent observations of a large number of DA and DB white dwarfs show evidence of debris disks, which are the remnants of old planetary systems. The infrared excess detected with Spitzer and the lines of heavy elements observed in their atmospheres with high-resolution spectroscopy converge on the idea that planetary material accretes onto these stars. Accretion rates have been derived by several authors with the assumption of a steady state between accretion and gravitational settling. The results are unrealistically different for DA and DB white dwarfs. Aims. When heavy matter is accreted onto stars, it induces an inverse μ-gradient that leads to fingering (thermohaline) convection. The aim of this letter is to study the impact of this specific process on the derived accretion rates in white dwarfs and on the difference between DA and DB. Methods. We solve the diffusion equation for the accreted heavy elements with a time-dependent method. The models we use have been obtained both with the IRAP code, which computes static models, and the La Plata code, which computes evolutionary sequences. Computations with pure gravitational settling are compared with computations that include fingering convection. Results. The most important result is that fingering convection has very important effects on DAs but is inefficient in DBs. When only gravitational settling is taken into account, the time-dependent computations lead to a steady state, as postulated by previous authors. When fingering convection is added, this steady state occurs much later. Conclusions. The surprising difference found in the past for the accretion rates derived for DA and DB white dwarfs disappears. The derived accretion rates for DAs are increased when fingering convection is taken into account, whereas those for DBs are not modified. More precise and developed results will be given in a forthcoming paper