We examine extensively the effect of the different crystallization pro
cesses, related to the presence of major and minor chemical species, o
n the binding energy and the cooling time of old white dwarfs. We use
improved equations of state for the solid and the liquid, and crystall
ization diagrams calculated within the modern theory of freezing. We s
how that, in spite of their small abundance, trace elements severely a
lter the cooling process and lengthen the cooling time of the star for
a given luminosity by several gigayears. In particular, Ne-22 is show
n to provide enough gravitational energy at crystallization to sustain
the star at the same luminosity for a time larger than the one due to
the crystallization of C/O itself. These calculations demonstrate the
necessity of including a proper treatment of crystallization in moder
n white dwarf cooling theory. We also consider the effect of an initia
l composition gradient in the distribution of carbon and oxygen throug
hout the star. Finally, we show that a substantial portion of the inte
rior of massive white dwarfs is already in a quantum state in the flui
d phase and that Debye cooling probably occurs prior to crystallizatio
n in these stars.