“…They embody the basic physics of such systems including the strong correlation effects in aromatic π-systems 84,89,92 and the self-interaction effects introduced by the nitrogen lone pairs. 51,59 Another advantage of these systems is that they are well-characterized experimentally [93][94][95][96][97][98][99][100][101][102][103][104][105][106][107][108] and well-studied theoretically by high-level wavefunction and Green's function methods. 93,104,[109][110][111][112][113][114][115][116][117][118][119][120][121][122] We calculate the electronic structure of benzene, pyridine, and the diazines using: (i) semi-local and hybrid DFT (ii) G 0 W 0 , (iii) ev-scGW, (iv) scGW 0 , (v) scGW, and (vi) G 0 W 0 combined with the second-order exchange self-energy (2OX), as an attempt to go beyond the GW approximation.…”