Structural and electronic properties of ZrX2)and HfX2 (X=S and Se) from first principles calculations.

Early transition metal dichalcogenides (TMDC), characterized by their quasi-two-dimensional layered structure, have attracted intensive interest due to their versatile chemical and physical properties, but a comprehensive understanding of their structural and electronic properties from a first-principles point of view is still lacking. In this work, four simple TMDC materials, MX(2) (M = Zr and Hf, X = S and Se), are investigated by the Kohn-Sham density functional theory (KS-DFT) with different local or semilocal exchange-correlation (xc) functionals and many-body perturbation theory in the GW approximation. Although the widely used Perdew-Burke-Ernzelhof (PBE) generalized gradient approximation (GGA) xc functional overestimates the interlayer distance dramatically, two newly developed GGA functionals, PBE-for-solids (PBEsol) and Wu-Cohen 2006 (WC06), can reproduce experimental crystal structures of these TMDC materials very well. The GW method, currently the most accurate first-principles approach for electronic band structures of extended systems, gives the fundamental band gaps of all these materials in good agreement with the experimental values obtained from optical absorption. The minimal direct gaps from GW are systematically larger than those measured from thermoreflectance by about 0.1-0.3 eV, implying that excitonic effects may be stronger than previously estimated. The calculated density of states from GW quasi-particle band energies agrees very well with photo-emission spectroscopy data. Ionization potentials of these materials are also computed by combining PBE calculations based on the slab model and GW quasi-particle corrections. The calculated absolute band energies with respect to the vacuum level indicate that that ZrS(2) and HfS(2), although having suitable band gaps for visible light absorption, cannot be used for overall water splitting as a result of mismatch of the conduction band minimum with the redox potential of H(+)/H(2).