Would it make sense to move this functionality into the topology class? That way the same code could be used, for instance, by the HDF5File loader when it's asked to load up only a subset of the atom_indices.