Density functional theory (DFT) calculations combined with surface thermodynamic arguments and the Gibbs-Curie-Wulff equilibrium morphology formalism have been employed to explore the effect of the reaction conditions, temperature (T), and gas-phase partial pressures (p(H2) and p(H2S)) on the stability of nickel sulfide (Ni3S2) surfaces. Furthermore, the strength and nature of chemical bonds for selected Ni3S2 surface cuts were investigated with the quantum theory of atoms in molecules methodology. A particular analysis of the electrostatic potential within this theoretical framework is performed to study the potential activity of nickel sulfide nanoparticles as hydrodesulfurization (HDS) catalysts. The calculated thermodynamic surface stabilities and the resulting equilibrium morphology model suggest that unsupported Ni3S, nanoparticles mainly expose (11 (1) over bar) and (11 (1) over bar) type surface faces in HDS conditions. Analysis of the electrostatic potential mapped onto a selected electron density isocontour (0.001 au) on those expose surface reveals a poor potential reactivity toward electron-donating reagents (i.e., low Lewis acidity). Consequently, a very low attraction between coordinatively unsaturated active sites (Lewis sites) exposed at the catalytic particles and the S atoms coming from reagent polluting molecules does inactive these kinds of particles for HDS.