# Print dipole moment contribution per orbital

Hi,

I know from the Hartree-Fock method, the electric part of dipole moment may be written as
-P_pq x_pq, where P_pq = \sum_i c_ip c_iq is the density matrix and x_pq = <p|x|q>

My question is, is there any method to print the contribution for each orbital, e.g., - c_ip c_iq x_pq? Possibly for natural orbital as well

For example, the following input


geom = """
nocom
noreorient
H 0. 0.  0.00
F  0. 0. 0.92
"""

mol = psi4.geometry(geom)

set {
reference rhf
basis cc-pvdz
}

set freeze_core false

energy('hf')
properties('hf', properties=['dipole'])


  Dipole Moment: [e a0]