Abstract
The metal-insulator transition of VO2 so far has evaded an accurate description by density functional theory. The screened hybrid functional of Heyd, Scuseria and Ernzerhof leads to reasonable solutions for both the low-temperature monoclinic and high-temperature rutile phases only if spin polarization is excluded from the calculations. We explore whether a satisfactory agreement with experiment can be achieved by tuning the fraction of Hartree Fock exchange (alpha) in the density functional. It is found that two branches of locally stable solutions exist for the rutile phase for 12.5% <= alpha <= 20%. One is metallic and has the correct stability as compared to the monoclinic phase, the other is insulating with lower energy than the metallic branch. We discuss these observations based on the V 3d orbital occupations and conclude that alpha = 10% is the best possible choice for spin-polarized VO2 calculations. (c) 2014 Elsevier B.V. All rights reserved.