When you say, “read the symmetries of my orbitals,” do you mean “access the orbital energies by symmetry”? If so, wfn.epsilon_a() is indeed the function to call. You just don’t want the to_array function. Use .nph instead.
When you say, “read the symmetries of my orbitals,” do you mean “access the orbital energies by symmetry”? If so, wfn.epsilon_a() is indeed the function to call. You just don’t want the to_array function. Use .nph instead.