Singlet fission (SF) offers opportunities for wavelength-selective processing of solar photons with an end goal of achieving higher efficiency inexpensive photovoltaic or solar-fuels-producing devices. In order to evaluate new molecular design strategies and for theoretical exploration of dynamics, it is important to put in place tools for efficient calculation of the electronic coupling between single-exciton reactant and multiexciton product states. For maximum utility, the couplings should be calculated at multiple nuclear geometries (rather than assumed constant everywhere, i.e., the Condon approximation) and we must be able to evaluate couplings for covalently linked multichromophore systems. With these requirements in mind, here we discuss the simplest methodology possible for rapid calculation of diabatic one-electron coupling matrix elements-based on Boys localization and rediagonalization of molecular orbitals. We focus on a covalent species called BT1 that juxtaposes two tetracene units in a partially cofacial geometry via a norbornyl bridge. In BT1, at the equilibrium C2v structure, the "nonhorizontal" couplings between HOMOs and LUMOs (t(HL) and t(LH)) vanish by symmetry. We then explore the impact of molecular vibrations through the calculation of t(AB) coupling gradients along 183 normal modes of motion. Rules are established for the types of motions (irreducible representations in the C2v point group) that turn on tHL and tLH values as well as for the patterns that emerge in constructive versus destructive interference of pathways to the SF product. For the best modes, calculated electronic coupling magnitudes for SF (at root-mean-squared deviation in position at 298 K), are within a factor of 2 of that seen for noncovalent tetracene dimers relevant to the molecular crystal. An overall "effective" electronic coupling is also given, based on the Stuchebrukhov formalism for non-Condon electron transfer rates.
Find this issue from the American Chemical Society