Delocalization effects in singlet fission: Comparing models with two and three interacting molecules

We present surface hopping simulations of singlet fission in 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene (ThBF). In particular, we performed simulations based on quantum mechanics/molecular mechanics (QM/MM) schemes in which either two or three ThBF molecules are inserted in the QM region and embedded in their MM crystal environment. Our aim was to investigate the changes in the photodynamics that are brought about by extending the delocalization of the excited states beyond the minimal model of a dimer. In the simulations based on the trimer model, compared to the dimer-based ones, we observed a faster time evolution of the state populations, with the largest differences associated with both the rise and decay times for the intermediate charge transfer states. Moreover, for the trimer, we predicted a singlet fission quantum yield of similar to 204%, which is larger than both the one extracted for the dimer (similar to 179%) and the theoretical upper limit of 200% for the dimer-based model of singlet fission. Although our study cannot account for the effects of extending the delocalization beyond three molecules, our findings clearly indicate how and why the singlet fission dynamics can be affected.