“…In fact, although the present study has focused on a class of quasi-1D polyaromatic systems, they exhibit the feature of localizing the radical spin states in the zigzag boundaries commonly observed in many other extended PAHs, for which the present results can be directly transferred. However, there are other classes of PAHs such as zethrenes or diindeno[b,i]anthracenes 57,58 to name a few, where the analysis of reactive centers is not so clear and where calculations as the ones presented here will provide important support by pinpointing the most chemically reactive centers and tuning their native radical structures through the introduction of donor or acceptor hetero-atomic defects in the carbon skeleton. Thus, the tunability of the diradical character of the acenes demonstrated in the present study through chemical substitutions is a potential route to be explored for designing organic molecules with enhanced nonlinear optical activity 59,60 or with optimized efficiency for singlet fission process in organic photo-voltaic devices 61,62 while retaining the feasibility of chemical synthesis.…”