I think it makes more sense to move all of the AO-basis integrals from Orbitals to Integrals and rename these classes to MOIntegrals and AOIntegrals, respectively. Most codes will require one or the other but not both, so I think this is better aligned with the expected usage. The new get_ao_1e_mean_field function can simply take mo_coefficients and mo_slices as arguments.
I think it makes more sense to move all of the AO-basis integrals from
OrbitalstoIntegralsand rename these classes toMOIntegralsandAOIntegrals, respectively. Most codes will require one or the other but not both, so I think this is better aligned with the expected usage. The newget_ao_1e_mean_fieldfunction can simply takemo_coefficientsandmo_slicesas arguments.