We present the results of DFT calculations of free-energy profiles along the reaction pathways starting from 50% coverage of MoS2 M-edge and 100% coverage of S-edge by sulfur, and leading to 37% coverage, i.e. creation of anionic vacancies, upon reduction by dihydrogen and production of H2S. Significant entropic and enthalpic corrections to electronic energies are deduced from the sets of normal modes vibration frequencies computed for all stationary and transition states. On that basis, we revisit and discuss the surface phase diagrams for M- and S-edges as a function of temperature. H-2 partial pressure and H2S/H-2 molar ratio, with respect to ranges of conditions relevant to industrial hydrotreating operations. We show that in such conditions, anionic vacancies on the M-edge, and surface SH groups on the M- and S-edges, may coexist at equilibrium. Moderate activation barriers connect stationary states along all paths explored. (C) 2011 Elsevier Inc. All rights reserved.