Dear Mohamed,
you are apparently reading here
https://www.diracprogram.org/doc/release-23/tutorials/analysis/projection_analysis.html
but those pages concern projection analysis, not Mulliken population analysis.
You have chosen to use the ZORA Hamiltonian
but please note what is written
Please note that these are only tested for energies and generally do not work for properties !
All the best,
Trond