“…Consequently, traditional TD-DFT methods for studying excited states miserably fail in describing singlet-triplet inversion, as was noted early on. 2,46,47 To address ST-inversion, high-quality but computationally demanding methods are typically used, spanning from ADC(2), [2][3][4]7,8 CC2, 2,3,8,9 EOM-CCSD, 2,4,5,48,49 Mukherjee's multireference coupled-cluster with noniterative triple excitations (Mk-MRCCSD(T)), 50 CC3, 51 Green's functions methods, 10 CIS(D), 3,6 SA-CASSCF, 3,8,50 and NEVPT2. 3,4,8,50 We adopt a different strategy and conduct an extensive study of various triangular boron carbon nitrides using the Pariser-Parr-Pople (PPP) model Hamiltonian to account for correlated electrons in π-conjugated molecules.…”