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\, \rm 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.

White dwarfs are compact objects with atmospheres containing mainly light elements, hydrogen or helium. Because of their surface high gravitational field, heavy elements diffuse downwards in a very short timescale compared to the evolutionary timescale, leaving the lightest ones on the top of the envelope. This results in the main classification of white dwarfs as hydrogen rich or helium rich. But many helium rich white dwarfs show also the presence of hydrogen traces in their atmosphere, whose origin is still unsettled. Here we study, by means of full evolutionary calculations, the case for a representative model of the "He-H-Z" white dwarfs, a sub-group of helium rich white dwarfs showing both heavy elements and a large amount of hydrogen in their atmosphere. We find it impossible to explain its hydrogen atmospheric content by the convective mixing of a primordial hydrogen present in the star. We conclude that the most likely explanation is the accretion of hydrogen rich material, presumably water-bearing, coming from a debris disk.

Ultra-massive DA WD stars are expected to harbor ONe cores resulting from the progenitor evolution through the Super-AGB phase. As evolution proceeds during the WD cooling phase, a crystallization process resulting from Coulomb interactions in very dense plasmas is expected to occur, leading to the formation of a highly crystallized core. Pulsating ultra-massive WDs offer a unique opportunity to infer and test the occurrence of crystallization in WD interiors as well as physical processes related with dense plasmas. We aim to assess the adiabatic pulsation properties of ultra-massive DA WD with ONe cores. We studied the pulsation properties of ultra-massive DA WD stars with ONe cores. We employed a new set of ultra-massive WD evolutionary sequences of models with stellar masses in the range 1.10 $\leq M_{star}/M_{sun} \leq$ 1.29 computed by taking into account the complete evolution of the progenitor stars and the WD stage. When crystallization set on in our models, we took into account latent heat release and also the expected changes in the core chemical composition that are due to phase separation according to a phase diagram suitable for O and Ne plasmas. We computed nonradial pulsation g-modes of our sequences of models at the ZZ Ceti phase by taking into account a solid core. We explored the impact of crystallization on their pulsation properties, in particular, the structure of the period spectrum and the distribution of the period spacings. We find that it would be possible, in principle, to discern whether a WD has a nucleus made of CO or a nucleus of ONe by studying the spacing between periods. The features found in the period-spacing diagrams could be used as a seismological tool to discern the core composition of ultra-massive ZZ Ceti stars, something that should be complemented with detailed asteroseismic analysis using the individual observed periods.

We present a set of full evolutionary sequences for white dwarfs with hydrogen-deficient atmospheres. We take into account the evolutionary history of the progenitor stars, all the relevant energy sources involved in the cooling, element diffusion in the very outer layers, and outer boundary conditions provided by new and detailed non-gray white dwarf model atmospheres for pure helium composition. These model atmospheres are based on the most up-to-date physical inputs. Our calculations extend down to very low effective temperatures, of $\sim 2\,500$~K, provide a homogeneous set of evolutionary cooling tracks that are appropriate for mass and age determinations of old hydrogen-deficient white dwarfs, and represent a clear improvement over previous efforts, which were computed using gray atmospheres.

A large fraction of white dwarf stars shows photospheric chemical composition polluted by heavy elements accreted from a debris disk. Such debris disks result from the tidal disruption of rocky planetesimals which had survived to whole stellar evolution from the main sequence to the final white dwarf stage. Determining the accretion rate of this material is an important step towards estimating the mass of the planetesimals and towards understanding the ultimate fate of the planetary systems. The accretion of heavy material with a mean molecular weight, $\mu$, higher than the mean molecular weight of the white dwarf outer layers, induces a double-diffusive instability producing fingering convection and extra-mixing. As a result, the accreted material is diluted deep into the star. We explore the effect of this extra-mixing on the abundance evolution of Mg, O, Ca, Fe and Si in the cases of the two well studied polluted DAZ white dwarfs: GD~133 and G~29-38. We performed numerical simulations of the accretion of material with a chemical composition similar to the bulk Earth one. We considered accretion rates from $10^{4}$~g/s to $10^{10}$~g/s. The double-diffusive instability develops on a very short time scale. The accretion rate needed to reproduce the observed abundances exceeds by more than 2 orders of magnitude the rate estimated by neglecting the fingering convection in the case of GD~133, and by approximately 1.7 dex in the case of G~29-38. Our numerical simulations show that fingering convection is an efficient mechanism to mix the accreted material and that it must be taken into account in the determination of accretion rates.

Fingering (thermohaline) convection has been invoked for several years as a possible extra-mixing which could occur in Red Giant stars due to the modification of the chemical composition induced by nuclear reactions in the hydrogen burning zone. Recent studies show however that this mixing is not sufficient to account for the needed surface abundances. A new prescription for fingering convection, based on 3D numerical simulations has recently been proposed (BGS). The resulting mixing coefficient is larger than the ones previously given in the literature. We compute models using this new coefficient and compare them to previous studies. We use the LPCODE stellar evolution code with the GNA generalized version of the mixing length theory to compute Red Giant models and we introduce fingering convection using the BGS prescription. The results show that, although the fingering zone now reaches the outer dynamical convective zone, the efficiency of the mixing is not enough to account for the observations. The fingering mixing coefficient should be increased by two orders of magnitude for the needed surface abundances to be reached. We confirm that fingering convection cannot be the mixing process needed to account for surface abundances in RGB stars.

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 \emph{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. When heavy matter is accreted onto stars, it induces an inverse $\mu$-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. 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. 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. 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.

We compute full evolutionary sequences of red giant branch stars close to the luminosity bump by including state of the art composition transport prescriptions for the thermohaline mixing regimes. In particular we adopt a self-consistent double-diffusive convection theory, that allows to handle the instabilities that arise when thermal and composition gradients compete against each other, and a very recent empirically motivated and parameter free asymptotic scaling law for thermohaline composition transport. In agreement with previous works, we find that during the red giant stage, a thermohaline instability sets in shortly after the hydrogen burning shell (HBS) encounters the chemical discontinuity left behind by the first dredge-up. We also find that the thermohaline unstable region, initially appearing at the exterior wing of the HBS, is unable to reach the outer convective envelope, with the consequence that no mixing of elements that produces a non-canonical modification of the stellar surface abundances occurs. Also in agreement with previous works, we find that by artificially increasing the mixing efficiency of thermohaline regions it is possible to connect both unstable regions, thus affecting the photospheric composition. However, we find that in order to reproduce the observed abundances of red giant branch stars close to the luminosity bump, thermohaline mixing efficiency has to be artificially increased by about 4 orders of magnitude from that predicted by recent 3D numerical simulations of thermohaline convection close to astrophysical environments. From this we conclude the chemical abundance anomalies of red giant stars cannot be explained on the basis of thermohaline mixing alone.

As a first step towards a comprehensive investigation of stellar motions within globular clusters, we present here the results of a study of stellar orbits in a mildly triaxial globular cluster that follows a circular orbit inside a galaxy. The stellar orbits were classified using the frequency analysis code of Carpintero and Aguilar and, as a check, the Liapunov characteristic exponents were also computed in some cases. The orbit families were obtained using different start spaces. Chaotic orbits turn out to be very common and while, as could be expected, they are particularly abundant in the outer parts of the cluster, they are still significant in the innermost regions. Their relevance for the structure of the cluster is discussed.