Density Functional Theory (UB3LYP/6-31111G(d,p)) calculations of the affinity of the pentaaqua nickel(II) complex for a set of phosphoryl [O@P(H)(CH 3 )(PhR)], imino [HN@C(CH 3 )(PhR)], thiocarbonyl [S@C(CH 3 )(PhR)] and carbonyl [O@C(CH 3 )(PhR)] ligands were performed, where R@NH 2 , OCH 3 , OH, CH 3 , H, Cl, CN, and NO 2 is a substituent at the para-position of a phenyl ring.The affinity of the pentaaqua nickel(II) complex for these ligands was analized and quantified in terms of interaction enthalpy (DH), Gibbs free energy (DG 298 ), geometric and electronic parameters of the resultant octahedral complexes. The DH and DG 298 results show that the ligand coordination strength increases in the following order: carbonyl < thiocarbonyl < imino < phosphoryl. This coordination strength order is also observed in the analysis of the metal-ligand distances and charges on the ligand atom that interacts with the Ni(II) cation. The electronic character of the substituent R is the main parameter that affects the strength of the metal-ligand coordination. Ligands containing electron-donating groups (NH 2 , OCH 3 , OH) have more exothermic DH and DG 298 than ligands with electron-withdrawing groups (Cl, CN, NO 2 ). The metal-ligand interaction decomposed by means of the energy decomposition analysis (EDA) method shows that the electronic character of the ligand modulates all the components of the metal-ligand interaction. The absolute softness of the free ligands is correlated with the covalent contribution to the instantaneous interaction energy calculated using the EDA method.