High efficient catalysts for hydrogen evolution reaction are crucial for hydrogen production. As potential catalysts for hydrogen evolution reaction, Ni m Mo n (m + n = 5) clusters have been studied in the present work on their structures, stabilities, and chemical reactivity toward H 2 O adsorption, OH bond breaking, and H adsorption with density functional theory based method. Mo is the active site for H 2 O adsorption and the adsorption strengths vary slightly with NiMo ratio. On the other hand, the split H prefers to form a triangle with Ni and Mo stabilized by the synergistic effect of Ni and Mo, i.e. adsorption on Mo and activation (and stabilization) by Ni. Electro-reduction makes the OH bond breaking barrierless on Ni 4 Mo, and this rationalizes the high catalytic efficiency of Ni 4 Mo in hydrogen evolution reaction (in agreement with experiment), and other H 2 O adsorption related reactions.