We study the evolution of close binary systems in order to account for the existence of the recently observed binary system containing the most massive millisecond pulsar ever detected, PSR J0740+6620, and its ultra-cool helium white dwarf companion. In order to find a progenitor for this object we compute the evolution of several binary systems composed by a neutron star and a normal donor star employing our stellar code. We assume conservative mass transfer. We also explore the effects of irradiation feedback on the system. We find that irradiated models also provide adequate models for the millisecond pulsar and its companion, so both irradiated and non irradiated systems are good progenitors for PSR J0740+6620. Finally, we obtain a binary system that evolves and accounts for the observational data of the system composed by PSR J0740+6620 (i.e. orbital period, mass, effective temperature and inferred metallicity of the companion, and mass of the neutron star) in a time scale smaller than the age of the Universe. In order to reach an effective temperature as low as observed, the donor star should have an helium envelope as demanded by observations.

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.