Working code for multiple core-valence exciations.  This code is not
very well tested yet.  The code for particles is incomplete (it just
needs a few lines in the sort_ee_mos() function).