The performance of hybrid density functionals in solid state chemistry:: the case of BaTiO3

The performance of hybrid exchange functionals is investigated within the framework of Density Functional Theory (DFT), in describing the properties of cubic and tetragonal ferroelectric phases of BaTiO3. A variable fraction zeta of non-local HF exchange is mixed with Becke's GGA formulation, and coupled with the LYP correlation functional. Standard DFT functionals (LDA, PW, B3LYP) and uncorrelated HF are also added for comparison. Inclusion of the non-local HF exchange has a noticeable effect on the calculated electronic density, and several observables show clear trends as a function of the mixing parameter zeta; these include not only the band gap, but also lattice parameter and bulk modulus, structural distortions, relative ionic sizes, electronic polarizability and polarization of the tetragonal phase. In this context, it is shown that hybrid functionals yield results that are generally in better agreement with experiment than local (LDA and GGA) DFT formulations. A higher fraction of HF exchange than the 20% optimized in the B3LYP functional for molecular systems is required in the solid state. This is due to the localization effect played by the Madelungfield on the electronic distribution; in this situation, the self-interaction error in local DFT functionals is more pronounced, and is empirically corrected by a high fraction of exact HF exchange component in the hybrid functional.