The composition of molecular species, in MolecularSpecies.Element.__call__ might be the bottleneck of many operations.

If I understand my code correctly, it first creates the corresponding permutation group, and then constructs the molecular species from it. This may not be ideal, because the latter has to decompose it into atomic species. However, we know that the composition of an atomic species and a molecular species is atomic.