Where in the source code is the density built?

Avogadro can plot densities out of molden files. An electron density is the square of two vectors of molecular orbitals. I was wondering where in the source code the electron density is built such that I am able to tinker around with how this is done. Specifically I am looking to replace certain d-orbitals in one of the MO vectors iff in the other MO vector is a p-orbital, etc. I’d imagine there to be a big loop structure somewhere in the source code where it is looped over all GTOs but I didn’t find it. If I have that, recompiling it and running it should be doable.

Thanks for the help!

See, for example in Avo2:

1 Like

This topic was automatically closed 3 days after the last reply. New replies are no longer allowed.