Move to MO::Basis::LCAO::Gaussian
MO::Basis::LCAO::Gaussian#overlap should return a matrix