The highest stability of the substitutional doping emerges from the interplay between chargetransfer, polarization, and electrostatic stabilizing physical effects. Otherwise, the adatom doping is mainly dominated by stabilizing charge transfer effects. Moreover, the magnitude of the polarization stabilizing effects is related to the polar character of the single closed-shell interactions between dopants and phosphorous atoms, while the fulfill of 3d orbitals controls the magnitude of stabilizing charge-transfer energies. Also, it was found that the coordination environment and ligand-field exerted on the dopant atom are related to the magnitude of the stabilizing electrostatics interactions. These results serve as a guideline to understand the stability of doped-phosphorene materials, expanding the understanding of the physics and chemistry of low-dimensional phosphorene-based materials.