MO::Basis::LCAO#shells should return a collection of BasisSet::LCAO::Gaussian::Shells
Move to MO::Basis::LCAO::Gaussian
MO::Basis::LCAO::Gaussian#overlap should return a matrix