Mulliken partitioning for two-electron operators

[lelemmen]

Since we already have an intuitive API for Mulliken-partitioning one-electron operators, we could think about implementing two-electron Mulliken partitioning.

This would allow:

• Energy decompositions (i.e. Mulliken-partitioning of the Hamiltonian)
• 'Local' S^2 calculations