Theoretical modeling of the diphosphene (HP=PH) ligand