Currently, the basis field is a map from Atom to BasisFunction objects. Back working on it I found it a little cumbersome and possibly fragile. So I want to add a field Shells to an array of BasisFunctions (possible name shells) and remove the basis field.