Dear Psi4 community members,

I am learning to use Psi4 and am still familiarizing myself with the software. At the moment, I am studying the bonding properties of Fe2O3, as a test. However, all the results are given in terms of Gaussian functions, where, for each shell, we have the coefficients and exponents of the Gaussians. My questions is: How can I project the Gaussian functions onto a basis of atomic orbitals? Can this be done? I would very appreciate any help.