As a result of competing physical mechanisms, the atmospheric composition of white dwarfs changes throughout their evolution, a process known as spectral evolution. Because of the ambiguity of their atmospheric compositions and the difficulties inherent to the modeling of their dense atmospheres, no consensus exists regarding the spectral evolution of cool white dwarfs ($T_{\rm eff}<6000$ K). In the previous papers of this series, we presented and observationally validated a new generation of cool white dwarf atmosphere models that include all the necessary constitutive physics to accurately model those objects. Using these new models and a homogeneous sample of 501 cool white dwarfs, we revisit the spectral evolution of cool white dwarfs. Our sample includes all spectroscopically identified white dwarfs cooler than 8300 K for which a parallax is available in Gaia DR2 and photometric observations are available in Pan-STARRS1 and 2MASS. Except for a few cool carbon-polluted objects, our models allow an excellent fit to the spectroscopic and photometric observations of all objects included in our sample. We identify a decrease of the ratio of hydrogen to helium-rich objects between 7500 K and 6250 K, which we interpret as the signature of convective mixing. After this decrease, hydrogen-rich objects become more abundant up to 5000 K. This puzzling increase, reminiscent of the non-DA gap, has yet to be explained. At lower temperatures, below 5000 K, hydrogen-rich white dwarfs become rarer, which rules out the scenario according to which accretion of hydrogen from the interstellar medium dominates the spectral evolution of cool white dwarfs.

A new semi-automatic approach is employed to carry out the abundance analysis of high-resolution spectra of HD41076 and HD148330 obtained recently with the spectropolarimetre ESPaDOnS at the CFHT. This approach allows to prepare in a semi-automatic mode the input data for the modified ZEEMAN2 code and to analyse several hundreds of line profiles in sequence during a single run. It also provides more information on abundance distribution for each chemical element at the deeper atmospheric layers. Our analysis of the Balmer profiles observed in the spectra of HD41076 and HD148330 has resulted in the estimates of their effective temperature, gravity, metallicity and radial velocity. The respective models of stellar atmosphere have been calculated with the code PHOENIX and used to carry out abundance analysis employing the modified ZEEMAN2 code. The analysis shows a deficit of the C, N, F, Mg, Ca, Ti, V, Cu, Y, Mo, Sm and Gd, and overabundance of Cr, Mn, Fe, Co, Ni, Sr, Zr, Ba, Ce, Nd and Dy in the stellar atmosphere of HD41076. In the atmosphere of HD148330, the C, N and Mo appear to be underabundant, while the Ne, Na, Al, Si, P, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, Sr, Y, Zr, Ba, Ce, Pr, Nd, Sm, Eu, Gd and Dy are overabundant. We also have found signatures of vertical abundance stratification of Fe, Ti, Cr and Mn in HD41076, and of Fe, Ti, V, Cr, Mn, Y, Zr, Ce, Nd, Sm and Gd in HD148330.

We present the "Montreal White Dwarf Database" (MWDD), an accessible database with sortable/filterable table and interactive plots that will, when fully completed, allow the community to explore the physical properties of all white dwarfs ever analyzed by the Montreal group, as well as display data and analyses from the literature. We present its current capability and show how it will continuously be updated to instantly reflect improvements made on both the theoretical and observational fronts.

We present the first results of the estimation of gravity and effective temperature for some poorly studied chemically peculiar stars that were recently observed with the spectropolarimeter ESPaDOnS at CFHT. We have analyzed the spectra of HD71030, HD95608 and HD116235 to determine their radial velocity, Vsin(i) and the average abundance of several chemical species. We have also analyzed our results to verify for possible vertical abundance stratification of iron and chromium in these stars.