Equation-of-motion spin-flip coupled-cluster method with single and double substitutions (EOM-SF-CCSD) is employed to study how substituents affect the electronic states' ordering in meta-xylylene diradicals. The electronegativity of substituents and the incorporation of a heteroatom are found to have a negligible effect. The effect of charges on energy gaps is much more pronounced, in agreemen...