Zinc porphene is a two-dimensional material made of fully fused zinc porphyrins in a tetragonal lattice. It has a fully conju-gated π-system, making it similar to graphene. Zinc porphene has recently been synthesized and shown to be a semiconduc-tor (Nat. Comm., 2023, 14, 6308.). This is in contrast with all previous predictions of its electronic structure, which indi-cated metallic conductivity. We show that the gap-opening in zinc porphene is caused by a Peierls distortion of its unit cell from square to rectangular, thus giving the first account of its electronic structure in agreement with experiment. Accounting for this distortion requires a proper treatment of electron delocalization, which can be done using hybrid functionals with a substantial amount of exact exchange. Such a functional, PBE38, is then applied to predict the properties of many first tran-sition row metalloporphenes, some of which have already been prepared. We find that changing the metal strongly affects the electronic structure of metalloporphenes, resulting in a rich variety of both metallic conductors and semiconductors, which may be of great of interest to molecular electronics and spintronics. Properties of these materials are mostly governed by the extent of the Peierls distortion and the number of electrons in their π system, analogous to changes in aromaticity observed in cyclic conjugated molecules upon oxidation or reduction. These results give an account of how the concept of antiaromaticity can be extended to periodic systems.
The efficiency of solar cells may be improved by using singlet fission (SF), in which one singlet exciton splits into two triplet excitons. SF occurs in molecular crystals. A molecule may crystallize in more than one form, a phenomenon known as polymorphism. Crystal structure may affect SF performance. In the common form of tetracene, SF is experimentally known to be slightly endoergic. A second, metastable polymorph of tetracene has been found to exhibit better SF performance. Here, we conduct inverse design of the crystal packing of tetracene using a genetic algorithm (GA) with a fitness function tailored to simultaneously optimize the SF rate and the lattice energy. The property-based GA successfully generates more structures predicted to have higher SF rates and provides insight into packing motifs associated with improved SF performance. We find a putative polymorph predicted to have superior SF performance to the two forms of tetracene, whose structures have been determined experimentally. The putative structure has a lattice energy within 1.5 kJ/mol of the most stable common form of tetracene.
Singlet fission (k(SF)) and excimer formation (k(EXC)) rate constants along with other photophysical properties of thin solid layers of 1,3-diphenylisobenzofuran and 11 of its fluorinated derivatives have been determined. The molecular properties of these compounds are similar, but their crystal packing varies widely. Most of them undergo singlet fission whereas excitation in others is trapped in excimers. The trend in rate constants k(SF) agrees qualitatively with results of calculations by a simplified version of the frontier orbital model for a molecular pair. The main shortcoming of the model is discussed.
Optimized geometries of macrocycles composed of four Pt(PR3)(2) linkers and four rods terminated with ethynyl or pyridyl moieties have been calculated at the PBEO/Def2-TZVPP//PBEO/Def2SVP level of density functional theory for all combinations of cis and trans configurations at the Pt centers. For uncharged complexes with four p,p'-bis(ethynyl)biphenyl rods and neutral Pt centers, the energy of the planar oval-shaped all-trans isomer lies 17 kcal/mol below that of the puckered square-shaped all-cis isomer. In this case, the electronic strain associated with the cis arrangement at the Pt atoms overrules the mechanical strain associated with the bending of triple bonds. For cationic complexes containing two bipyridyl and two p,p'bis(ethynyl)biphenyl rods, with a +1 charge on each Pt center, a puckered rectangular structure with all-cis Pt centers is found to be 38 kcal/mol lower in energy than the isomer with all four Pt centers trans. The results have been analyzed in terms of three additive factors, referred to as electronic, mechanical, and electrostatic strains.
The crucial role of molecular packing in singlet fission in dimers, aggregates, and solids is summarized and the key roles of Davydov splitting and the size of the electronic matrix element are pointed out. Simple rules and computer algorithms for these effects are described and their shortcomings are listed.
09:30 10:00 1.1-O1 Alexandr Zaykov (Institute of Organic Chemistry and Biochemistry of the CAS), Josef Michl, Zdeněk Havlas, Eric Buchanan, Milena Jovanović Singlet Fission: Chromophores for Exciton Downconversion 10:00 10:30 1.1-O2 Benjamin Daiber (Center for Nanophotonics, AMOLF, Science Park 104, The Netherlands), Koen v.d. Hoven, Joris Y. Bodin, Stefan Luxembourg, Moritz Futscher, Bruno Ehrler Efficiency Potential and Application of Singlet Fission Enhanced Silicon Solar Cells using Different Energy Transfer 10:30 11:00 Coffee Break Exciup 1.2
A procedure is described for unbiased identification of all π-electron chromophore pair geometry choices that locally maximize the rate of conversion of a singlet exciton into a singlet biexciton (triplet pair), using a simplified version of the diabatic frontier orbital model of singlet fission (SF). The resulting approximate optimal geometries provide insight and are expected to represent useful starting points for searches by more advanced methods. The general procedure is illustrated on a pair of ethylenes as the simplest model of a π-electron system, but it is applicable to pairs of much larger molecules, with dozens of non-hydrogen atoms, and not necessarily planar. We first examine the value of |TA|2, the square of the electronic matrix element for SF with initial excitation fully localized on partner A, on a grid of several billion geometries within the six-dimensional space of physically realizable possibilities. Several of the optimized pair geometries are somewhat unexpected, but all are found to follow the qualitative guidance proposed earlier. In the neighborhood of each local maximum of |TA|2, consideration of mixing with charge-transfer configurations and of excitonic interaction between partners A and B determines the SF energy balance and yields squared matrix elements |T*|2 and |T**|2 for the lower and upper excitonic states S* and S**, respectively. Assuming Boltzmann populations of these states, the geometry is further optimized to maximize k, the sum of the SF rates obtained from Marcus theory, and this reorders the suitable geometries substantially. At 87 pair geometries, the |T*|2 and |T**|2 values are compared with those obtained from high-level ab initio nonorthogonal configuration interaction calculations and found to follow the same trend. Finally, the biexciton binding energy at the optimized geometries is calculated. Altogether, 13 significant local maxima of SF rate for a pair of ethylenes are identified in the physically relevant part of space that avoids molecular interpenetration in the hard-sphere approximation. The three best geometries are twist-stacked, slip-stacked, and L-shaped. The maxima occur at the (five-dimensional) surfaces of seven six-dimensional "parent" regions of space centered at physically inaccessible geometries at which the calculated SF rate is very large but the two ethylenes interpenetrate. The results are displayed in interactive graphics. The computer code ("Simple") written for these calculations is flexible in that it permits a choice of performing the search for local maxima in six dimensions on |TA|2, |T*|2, or k. It is available as freeware at https://cloud.uochb.cas.cz/simple .
A simplified version of the frontier orbital model has been applied to pairs of C_2, C_2v, C_s, and C_1 symmetry 1,3-diphenylisobenzofuran rotamers to determine their best packing for fast singlet fission (SF). For each rotamer the square of the electronic matrix element for SF was calculated at 2.2 × 10^9 pair geometries and a few thousand most significant physically accessible local maxima were identified in the six-dimensional space of mutual arrangements. At these pair geometries, SF energy balance was evaluated, relative SF rate constants were approximated using Marcus theory, and the SF rate constant k_SF was maximized by further optimization of the geometry of the molecular pair. The process resulted in 142, 67, 214, and 291 unique geometries for the C_2, C_2v, C_s, and C_1 symmetry molecular pairs, respectively, predicted to be superior to the C_2 symmetrized known crystal pair structure. These optimized pair geometries and their triplet biexciton binding energies are reported as targets for crystal engineering and/or covalent dimer synthesis, and as possible starting points for high-level pair geometry optimizations.
Crystal structures, singlet fission (SF) rate constants, and other photophysical properties are reported for three fluorinated derivatives of 1,3-diphenylisobenzofuran and compared with those of the two crystal forms of the parent. The results place constraints on the notion that the effects of molecular packing on SF rates could be studied separately from effects of chromophore structural changes by examining groups of chromophores related by weakly perturbing substitution if their crystal structures are different. The results further provide experimental evidence that dimer-based models of SF are not sufficiently general and that trimer- and possibly even higher oligomer-based or many-body models need to be formulated.
In search for a qualitative understanding of the effects of molecular packing on singlet fission (SF) rate, a simplified version of the frontier orbital model is described and illustrated on a pair of tetracene molecules. To identify all favorable physically accessible pair geometries, all significant local maxima of the square of the electronic matrix element for SF have been located within the six-dimensional space of possible arrangements of two rigid bodies, using a grid of over 4.7 x 10(8) pair geometries. Those at which the molecules interpenetrate were excluded using a hard-sphere model. The effects of intermolecular interaction on the SF energy balance and thus its rate constant k(SF) were approximated using Marcus theory at each of the maxima using the same simplified version of the frontier orbital model. Starting at these local maxima, the pair geometries were optimized for maximum k(SF) and the 21 best are reported along with their computed Davydov splitting and triplet biexciton binding energies. The optimal pair structures at the resulting maxima follow qualitative rules published previously and further elaborated here.
A simplified version of the frontier orbital model for a noncovalent dimer is used to derive guidelines for dimer geometries that maximize the square of the electronic matrix element for singlet fission. The use of the guidelines requires only the knowledge of the highest occupied and lowest unoccupied orbital of the monomer and the overlaps of the atomic orbitals on partner A with those on partner B.
After a brief review of electronic aspects of singlet fission, we describe a systematic simplification of the frontier orbital (HOMO/LUMO) model of singlet fission and Davydov splitting in a pair of rigid molecules. In both instances, the model includes electron configurations representing local singlet excitation on either chromophore, charge transfer in either direction, and triplet excitation in both chromophores (biexciton). The resulting equations are simple enough to permit complete searches for local extrema of the square of the electronic matrix element and to evaluate the effect of intermolecular interactions on the exoergicity of singlet fission and on the biexciton binding energy in the six-dimensional space of rigid dimer geometries. The procedure is illustrated on results for the six best geometries for dimers of ethylene and of an indigoid heterocycle with 24 carbon, nitrogen, and oxygen atoms.