  It has been shown that the shape of the luminosity function of white dwarfs
can be a powerful tool to check for the possible existence of DFSZ-axions. In
particular, Isern et al. (2008) showed that, if the axion mass is of the order
of a few meV, then the white dwarf luminosity function is sensitive enough to
detect their existence. For axion masses of about $m_a > 5$ meV the axion
emission can be a primary cooling mechanism for the white dwarf and the
feedback of the axion emission into the thermal structure of the white dwarf
needs to be considered. Here we present computations of white dwarf cooling
sequences that take into account the effect of axion emission in a self
consistent way by means of full stellar evolution computations. Then, we study
and discuss the impact of the axion emission in the white dwarf luminosity
function.
