The analysis becomes more complex for molecules with two heteroatoms as (i) there is a variation in the S0 aromaticity with the relative positions of the two heteroatoms, (ii) with two lone-pairs there are several n,* states since there are two (near-degenerate) lone-pair orbitals in addition to the two (near-degenerate) * orbitals (potentially leading to multiconfigurational character), and (iii) there can be a variation among diheteroaromatics as to which n,* state is lowest in energy. To facilitate, we split the diheteroaromatics in two subgroups: one composed of those with two different heteroatoms ( -14) and one of those with two equal heteroatoms (7 -9 and 15 -19). Additionally, we consider the protonated species of two compounds in the latter subgroup (9H + and 19H + ). Those with two different heteroatoms should, viewed simplistically, have the highest n orbital dominated by the least electronegative element, and their first n,* state may resemble those of the monoheteroaromatics with the same heteroatom. Yet, we will see that this is not necessarily the case. We calculate MCI throughout the group and MICD for selected compounds.