Triplet population dynamics of solution cast films of isolated polymorphs of 6,13‐bis(triisopropylsilylethynyl) pentacene (TIPS‐Pn) provide quantitative experimental evidence that triplet excitation energy transfer is the dominant mechanism for correlated triplet pair (CTP) separation during singlet fission. Variations in CTP separation rates are compared for polymorphs of TIPS‐Pn with their triplet diffusion characteristics that are controlled by their crystal structures. Since triplet energy transfer is a spin‐forbidden process requiring direct wavefunction overlap, simple calculations of electron and hole transfer integrals are used to predict how molecular packing arrangements would influence triplet transfer rates. The transfer integrals reveal how differences in the packing arrangements affect electronic interactions between pairs of TIPS‐Pn molecules, which are correlated with the relative rates of CTP separation in the polymorphs. These findings suggest that relatively simple computations in conjunction with measurements of molecular packing structures may be used as screening tools to predict a priori whether new types of singlet fission sensitizers have the potential to undergo fast separation of CTP states to form multiplied triplets.

Computational materials discovery efforts utilize hundreds or thousands of density functional theory (DFT) calculations to predict material properties. Historically, such efforts have performed calculations at the generalized gradient approximation (GGA) level of theory due to its efficient compromise between accuracy and computational reliability. However, high-throughput calculations at the higher metaGGA level of theory are becoming feasible. The Strongly Constrainted and Appropriately Normed (SCAN) metaGGA functional offers superior accuracy to GGA across much of chemical space, making it appealing as a general-purpose metaGGA functional, but it suffers from numerical instabilities that impede it's use in high-throughput workflows. The recently-developed r2SCAN metaGGA functional promises accuracy similar to SCAN in addition to more robust numerical performance. However, its performance compared to SCAN has yet to be evaluated over a large group of solid materials. In this work, we compared r2SCAN and SCAN predictions for key properties of approximately 6,000 solid materials using a newly-developed high-throughput computational workflow. We find that r2SCAN predicts formation energies more accurately than SCAN and PBEsol for both strongly- and weakly-bound materials and that r2SCAN predicts systematically larger lattice constants than SCAN. We also find that r2SCAN requires modestly fewer computational resources than SCAN and offers significantly more reliable convergence. Thus, our large-scale benchmark confirms that r2SCAN has delivered on its promises of numerical efficiency and accuracy, making it a preferred choice for high-throughput metaGGA calculations.

We use native vibrational modes of the model singlet fission chromophore 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS-Pn) to examine the origins of singlet fission in solution between molecules that are not tethered by a covalent linkage. We use the C—H stretch modes of TIPS side groups of TIPS-Pn to demonstrate that singlet fission does not occur by diffusive encounter of independent molecules in solution. Instead, TIPS-Pn molecules aggregate in solution through their TIPS side groups. This aggregation breaks the symmetry of the TIPS-Pn molecules and enables the formation of triplets to be probed through the formally symmetry forbidden symmetric alkyne stretch mode of the TIPS side groups. The alkyne stretch modes of TIPS-Pn are sensitive to the electronic excited states present during the singlet fission reaction and provide unique signatures of the formation of triplets following the initial separation of triplet pair intermediates. These findings highlight the opportunity to leverage structural information from vibrational modes to better understand intermolecular interactions that lead to singlet fission.

Interaction of oxygen vacancies in BaTiO3 as simulated using the ReaxFF reactive force field.