Output orbital resolved partial charges

[awvwgk]

Is your feature request related to a problem? Please describe.

Currently, we provide only shell-resolved and atom-resolved partial charges, we could also compute orbital-resolved partial charges

Describe the solution you'd like

Partition the reference occupation from the shells to the orbitals to compute partial charges on each orbital