Ab initio study of 2H-MoS2 using Hay and Wadt effective core pseudo-potentials for modelling the (10(1)over-bar-0) surface structure

Periodic Hartree-Fock and DFT methods were employed to calculate the geometric and electronic properties of bulk 2H-MoS2 and of its catalytically important (1010) edge structure. The core electrons of molybdenum and sulfur were represented by the effective core pseudo-potentials, developed by Hay and Wadt. For the calculations (1010) type surface structures were generated by cutting sections from an 2H-MoS2 crystal, which consists of four, six or eight rows of molybdenum and sulfur atoms. The calculated elastic constants follow the experimental constants and show that 2H-MoS2 is an anisotropic covalent compound held together by weak dispersion interactions between neighbouring S-Mo-S units. The relaxation of the (1010) surface of 2H-MOS2 leads to an inner relaxation of the Mo atoms and the formation of weakly coupled surface S-S species. As a result of the surface relaxation, the electronic charge of the surface states is different to that of the bulk states: empty Mo d states move into the band gap, and the surface becomes a better electron acceptor. In all calculations the model consisting of six rows of sulfur and molybdenum atoms is the most favourable one, because it provides an accurate description of the 2H-MoS2 surface structure and enables calculations at reasonable computation time. The calculations based on pseudo-potentials are in good agreement with all-electron calculations and confirm that, in contrast to the semiconducting bulk, the (1010) 2H-MoS2 surface has metallic properties.