Here the PBEsol generalized gradient approximation (GGA), TPSS meta-GGA, and HSE06 hybrid functionals have been employed to demonstrate the XC-functional effects on the calculated structural, phonon and thermodynamic properties of rutile and anatase TiO2.  Nevertheless, the thermodynamic properties of rutile and anatase TiO2 estimated by harmonic approximations are in excellent experimental agreement and are effectively invariant to the level of theory employed in the DFT XC-functional.
THERM