What We Can Learn from the Norms of One-Particle Density Matrices, and What We Can’t: Some Results for Interstate Properties in Model Singlet Fission Systems

The Role of Double Excitations in Exciton Dynamics of Multiazobenzenes: Trisazobenzenophane as a Test Case

Molecular exciton dynamics underlie energy and charge transfer processes in organic multichromophoric systems. A particularly interesting class of the latter is multiphotochromic systems made of molecules capable of photochemical transformations. Exciton dynamics in assemblies of photoswitches have been recently investigated using either the molecular exciton model or supermolecular configuration interaction (CI) singles, both approaches being based on a semiempirical Hamiltonian and combined with surface hopping molecular dynamics. Here, we study how inclusion of double excitations in nonadiabatic dynamics simulations affects exciton dynamics of multiazobenzenes, using trisazobenzenophane as an example. We find that both CI singles and CI singles and doubles yield virtually the same time scale of dynamical exciton localization, ∼50 fs for the studied multiazobenzene. However, inclusion of double excitations considerably affects the excited state lifetimes and isomerization quantum yields.

Coupled-cluster treatment of complex open-shell systems: the case of single-molecule magnets

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Accurate Electronic and Optical Properties of Organic Doublet Radicals Using Machine Learned Range-Separated Functionals

Luminescent organic semiconducting doublet-spin radicals are unique and emergent optical materials because their fluorescent quantum yields (Φfl) are not compromised by the spin-flipping intersystem crossing (ISC) into a dark high-spin state. The multiconfigurational nature of these radicals challenges their electronic structure calculations in the framework of single-reference density functional theory (DFT) and introduces room for method improvement. In the present study, we extended our earlier development of ML-ωPBE [J. Phys. Chem. Lett., 2021, 12, 9516-9524], a range-separated hybrid (RSH) exchange-correlation (XC) functional constructed using the stacked ensemble machine learning (SEML) algorithm, from closed-shell organic semiconducting molecules to doublet-spin organic semiconducting radicals. We assessed its performance for a new test set of 64 doublet-spin radicals from five categories while placing all previously compiled 3926 closed-shell molecules in the new training set. Interestingly, ML-ωPBE agrees with the nonempirical OT-ωPBE functional regarding the prediction of the molecule-dependent range-separation parameter (ω), with a small mean absolute error (MAE) of 0.0197 a0-1, but saves the computational cost by 2.46 orders of magnitude. This result demonstrates an outstanding domain adaptation capacity of ML-ωPBE for diverse organic semiconducting species. To further assess the predictive power of ML-ωPBE in experimental observables, we also applied it to evaluate absorption and fluorescence energies (Eabs and Efl) using linear-response time-dependent DFT (TDDFT), and we compared its behavior with nine popular XC functionals. For most radicals, ML-ωPBE reproduces experimental measurements of Eabs and Efl with small MAEs of 0.299 and 0.254 eV, only marginally different from those of OT-ωPBE. Our work illustrates a successful extension of the SEML framework from closed-shell molecules to doublet-spin radicals and will open the venue for calculating optical properties for organic semiconductors using single-reference TDDFT.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Quantification of the Ionic Character of Multiconfigurational Wave Functions: The Qat Diagnostic

The complete active space self-consistent field (CASSCF) method is a cornerstone in modern excited-state quantum chemistry providing the starting point for most common multireference computations. However, CASSCF, when used with a minimal active space, can produce significant errors (>2 eV) even for the excitation energies of simple hydrocarbons if the states of interest possess ionic character. After illustrating this problem in some detail, we present a diagnostic for ionic character, denoted as Q at, that is readily computed from the transition density. A set of 11 molecules is considered to study errors in vertical excitation energies. State-averaged CASSCF obtains a mean absolute error (MAE) of 0.87 eV for the 34 singlet states considered. We highlight a strong correlation between the obtained errors and the Q at diagnostic, illustrating its power to predict problematic cases. Conversely, using multireference configuration interaction with single and double excitations and Pople's size extensivity correction (MR-CISD+P), excellent results are obtained with an MAE of 0.11 eV. Furthermore, correlations with the Q at diagnostic disappear. In summary, we hope that the presented diagnostic will facilitate reliable and user-friendly multireference computations on conjugated organic molecules.

Spin-Orbit Couplings of Open-Shell Systems with Restricted Active Space Configuration Interaction

In this work we perform electronic structure calculations to unravel the origin of spin-orbit couplings (SOCs) in open-shell molecules. For that, we select systems displaying di or polyradical character, e.g., trimethylene, and analyze the changes in the magnitude of SOC constants along molecular distortions of ethylene and in the presence of intermolecular interactions between open and closed-shell moieties in the O₂-C₂H₄ system. Calculations were performed by using nonrelativistic wave functions obtained with the restricted active space configuration interaction (RASCI) method, in conjunction with a recent implementation for the calculation of SOC based on the spin-orbit mean field approximation. Our results demonstrate the suitability of RASCI in the calculation of SOCs of open-shell systems, while providing a deep understanding of the relationship between couplings and the nature of the electronic states. Moreover, we introduce a new definition of the SOC constant for the study of molecular aggregates.

On the Cooperative Origin of Solvent-Enhanced Symmetry-Breaking Charge Transfer in a Covalently Bound Tetracene Dimer Leading to Singlet Fission

Singlet fission in covalently bound acene dimers in solution is driven by the interplay of excitonic and singlet correlated triplet ¹(TT) states with intermediate charge-transfer states, a process which depends sensitively on the solvent environment. We use high-level electronic structure methods to explore this singlet fission process in a linked tetracene dimer, with emphasis on the symmetry-breaking mechanism for the charge-transfer (CT) states induced by low-frequency antisymmetric vibrations and polar/polarizable solvents. A combination of the second-order algebraic diagrammatic construction (ADC(2)) and density functional theory/multireference configuration interaction (DFT/MRCI) methods are employed, along with a state-specific conductor-like screening model (COSMO) solvation model in the former case. This work quantifies, for the first time, an earlier mechanistic proposal [Alvertis et al., J. Am. Chem. Soc. 2019, 141, 17558] according to which solvent-induced symmetry breaking leads to a high-energy CT state which interacts with the correlated triplet state, resulting in singlet fission. An approximate assessment of the nonadiabatic interactions between the different electronic states underscores that the CT states are essential in facilitating the transition from the bright excitonic state to the ¹(TT) state leading to singlet fission. We show that several types of symmetry-breaking inter- and intra-fragment vibrations play a crucial role in a concerted mechanism with the solvent environment and with the symmetric inter-fragment torsion, which tunes the admixture of excitonic and CT states. This offers a new perspective on how solvent-induced symmetry-breaking CT can be understood and how it cooperates with intramolecular mechanisms in singlet fission.

Atomic Group Decomposition of Charge Transfer Excitation Global Indexes

A model for decomposing the Le Bahers, Adamo, and Ciofini Charge Transfer (CT) Excitations global indexes ( J. Chem. Theory Comput. 2011, 7, 2498-2506) into molecular subdomains contributions is presented and a software, DOCTRINE (atomic group Decomposition Of the Charge TRansfer INdExes) for the implementation of this novel model has been coded. Although our method applies to any fuzzy or to any disjoint exhaustive partitioning of the real space, it is here applied using a definition of chemically relevant molecular subdomains based on the Atoms in Molecules Bader basins. This choice has the relevant advantage of associating intra or inter subdomain contributions to rigorously defined quantum objects, yet bearing a clear chemical meaning. Our method allows for a quantitative evaluation of the subdomain contributions to the charge transfer, the charge transfer excitation length and the dipole moment change upon excitation. All these global indexes may be obtained either from the electron density increment or the electron density depletion upon excitation. However, the subdomain contributions obtained from the two distributions generally differ, therefore allowing to distinguish whether the contribution to a given property of a given subdomain is dominated by one of the two distributions or if both are playing a significant role. As a toy system for the first application of our model, a typical [D-π-A, π = conjugated bridge] compound belonging to the merocyanine dyes family is selected, and the first four excited states of this compound in a strongly polar protic solvent and in a weakly polar solvent are thoroughly investigated.

Modulating singlet fission through interchromophoric rotation

Singlet fission (SF) is a spin-allowed, exciton-multiplying phenomenon that can be utilized to improve the efficiency of organic solar cells. It is well-understood that SF is sensitive to the local crystal morphology and an appropriately balanced coupling is essential to facilitate efficient SF. In this study, we show how the interchromophoric rotation selectively modulates the interaction between the monomer frontier molecular orbitals, promoting both fast and exothermal SF. We evaluate the effective electronic coupling for SF (V_SF), the square of which is proportional to the SF rate, and the effective energies of the Frenkel exciton (FE/S₁S₀) and triplet pair exciton (TT) in a terrylene dimer model. Optimal interplanar rotation of the chromophoric moieties in slip-stacked arrangements pulls the effective energy of the TT state below that of the FE state. Consequently, SF is favored over competing pathways such as excimer formation, thereby enhancing the overall triplet yield. This work represents a step towards improvising the molecular design guidelines for SF and understanding the importance of interchromophoric rotation over the conventional slip-stacked arrangements for achieving favorable intermolecular electronic coupling towards efficient SF.

Not dark yet for strong light-matter coupling to accelerate singlet fission dynamics

Polaritons are unique hybrid light-matter states that offer an alternative way to manipulate chemical processes. In this work, we show that singlet fission dynamics can be accelerated under strong light-matter coupling. For superexchange-mediated singlet fission, state mixing speeds up the dynamics in cavities when the lower polariton is close in energy to the multiexcitonic state. This effect is more pronounced in non-conventional singlet fission materials in which the energy gap between the bright singlet exciton and the multiexcitonic state is large ( > 0.1 eV). In this case, the dynamics is dominated by the polaritonic modes and not by the bare-molecule-like dark states, and, additionally, the resonant enhancement due to strong coupling is robust even for energetically broad molecular states. The present results provide a new strategy to expand the range of suitable materials for efficient singlet fission by making use of strong light-matter coupling.

Spin-Orbit Natural Transition Orbitals and Spin-Forbidden Transitions

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T₂-S₁ state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb³⁺ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)₃].

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Pushing the Limits of the Donor-Acceptor Copolymer Strategy for Intramolecular Singlet Fission

Donor-acceptor (D-A) copolymers have shown great potential for intramolecular singlet fission (iSF). Nonetheless, very few design principles exist for optimizing these systems for iSF, with very little knowledge about how to engineer them for this purpose. In recent work, a fundamental trade-off between the main electronic ingredients required for iSF capable D-A coplanar copolymers was revealed. Still, further investigations are needed to understand these limitations and learn how to bypass them. In this work, we propose to induce torsion as an effective way to circumvent the limits of the coplanar approach. We disclose the potential of noncoplanar copolymers with inherently low triplet energies that encompass all the characteristics required for iSF beyond the limiting values associated with fully coplanar systems. Our findings shed some light on the electronic structure aspects of D-A copolymers for iSF and offer a new avenue for the rational design of novel promising candidates.

Theoretical Study on Singlet Fission Dynamics in Slip-Stack-like Pentacene Ring-Shaped Aggregate Models

We investigate the singlet fission (SF) dynamics of a slip-stack-like pentacene ring-shaped aggregate model, which is constructed by rotating each pentacene unit around its longitudinal axis in an H-aggregate ring. The aggregate size (N) and rotation angle (α) dependences of SF rates and double triplet (TT) yields are clarified using the quantum master equation method. It is found that there exist optimal ranges of the rotation angle α for each N, yielding efficient SF with high SF rates and TT yields. For example, in an 8-mer model, SF rates at α = 23 and 43° are 18.9 and 38.6 times as high as that at α = 30°, respectively, and the TT yields are as high as 0.871, 0.988, and 0.882 at α = 23, 30, and 43°, respectively. Analysis of the relative relaxation factors shows that the many-to-many relaxation paths from adiabatic Frenkel exciton (FE)-like states to TT-like states are opened by tuning α at relevant aggregate sizes, causing fast and high-TT-yield SF, and efficient SF occurs at α = 40° for medium N (7 ≤ N ≤ 10) or at α = 30° for large N (>10). This mechanism is interpreted by the second-order perturbation theory for electronic couplings. Namely, the inequality in the energies of charge-transfer states [CA and AC states, where the cation (C) and anion (A) are located at two neighboring sites in anticlockwise and clockwise directions, respectively] and the change in the amplitude and sign of the couplings between the FE, CT, and TT states are found to cause quantum superposition of the FE and TT states, which contribute to the high TT yield and SF rate. The present results contribute to a deeper understanding of SF dynamics in ring-shaped aggregates as well as to the development of their new design guidelines.

Short-range DFT energy correction to multiconfigurational wave functions for open-shell systems

Electronic structure methods emerging from the combination of multiconfigurational wave functions and density functional theory (DFT) aim to take advantage of the strengths of the two nearly antagonistic theories. One of the common strategies employed to merge wave function theory (WFT) with DFT relies on the range separation of the Coulomb operator in which DFT functionals take care of the short-distance part, while long-range inter-electronic interactions are evaluated by using the chosen wave function method (WFT-srDFT). In this work, we uncover the limitations of WFT-srDFT in the characterization of open-shell systems. We show that spin polarization effects have a major impact on the (short-range) DFT exchange energy and are of vital importance in order to provide a balanced description between closed and open-shell configurations. We introduce different strategies to account for spin polarization in the short range based on the definition of a spin polarized electron density and with the use of short-range exact exchange. We test the performance of these approaches in the dissociation of the hydrogen molecule, the calculation of energy gaps in spin-triplet atoms and molecular diradicals, and the characterization of low-lying states of the gallium dimer. Our results indicate that the use of short-range DFT correlation in combination with a (full-range) multiconfigurational wave function might be an excellent approach for the study of open-shell molecules and largely improves the performance of WFT and WFT-srDFT.

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. I. Theory and implementation

X-ray absorption creates electron vacancies in the core shell. These highly excited states often relax by Auger decay-an autoionization process in which one valence electron fills the core hole and another valence electron is ejected into the ionization continuum. Despite the important role of Auger processes in many experimental settings, their first-principles modeling is challenging, even for small systems. The difficulty stems from the need to describe many-electron continuum (unbound) states, which cannot be tackled with standard quantum-chemistry methods. We present a novel approach to calculate Auger decay rates by combining Feshbach-Fano resonance theory with the equation-of-motion coupled-cluster single double (EOM-CCSD) framework. We use the core-valence separation scheme to define projectors into the bound (square-integrable) and unbound (continuum) subspaces of the full function space. The continuum many-body decay states are represented by products of an appropriate EOM-CCSD state and a free-electron state, described by a continuum orbital. The Auger rates are expressed in terms of reduced quantities, two-body Dyson amplitudes (objects analogous to the two-particle transition density matrix), contracted with two-electron bound-continuum integrals. Here, we consider two approximate treatments of the free electron: a plane wave and a Coulomb wave with an effective charge, which allow us to evaluate all requisite integrals analytically; however, the theory can be extended to incorporate a more sophisticated description of the continuum orbital.

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit

For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S0S1) and the multiexcitonic state (1TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc-pVnZ series of basis sets, that both the aug-cc-pVDZ and 6-311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

Unconventional singlet fission materials

Singlet fission (SF) is a photophysical downconversion pathway, in which a singlet excitation transforms into two triplet excited states. As such, it constitutes an exciton multiplication generation process, which is currently at the focal point for future integration into solar energy conversion devices. Beyond this, various other exciting applications were proposed, including quantum cryptography or organic light emitting diodes. Also, the mechanistic understanding evolved rapidly during the last year. Unfortunately, the number of suitable SF-chromophores is still limited. This is per se problematic, considering the wide range of envisaged applicability. With that in mind, we emphasize uncommon SF-scaffolds and outline requirements as well as strategies to expand the chromophore pool of SF-materials.

Energetics and optimal molecular packing for singlet fission in BN-doped perylenes: electronic adiabatic state basis screening

Singlet fission has the potential to increase the efficiency of photovoltaic devices, but the design of suitable chromophores is notoriously difficult. Both the electronic properties of the monomer and the packing motif in the crystal have a big impact on the singlet fission efficiency. Using perylene as an example, it is shown that doping with boron and nitrogen not only helps to align the energy levels but also shifts the stacking position that is optimal for singlet fission. Among all perylene derivatives doped with one or two BN groups, we identify the most suitable isomer for singlet fission with the help of TD-DFT and CASPT2 calculations. The optimal relative disposition of the two monomer units in a cofacially stacked homodimer is explored using two semiempirical models for the singlet fission rate: The first one is the well-known diabatic frontier orbital model, while the second treats singlet fission as a non-adiabatic transition and approximates the rate as the length squared of the non-adiabatic coupling vector between eigenfunctions of the diabatic Hamiltonian.

Overlap-Driven Splitting of Triplet Pairs in Singlet Fission

We analyze correlated-triplet-pair (TT) singlet-fission intermediates toward two-triplet separation (T...T) using spin-state-averaged density matrix renormalization group electronic-structure calculations. Specifically, we compare the triplet-triplet exchange (J) for tetracene dimers, bipentacene, a subunit of the benzodithiophene-thiophene dioxide polymer, and a carotenoid (neurosporene). Exchange-split energy gaps of J and 3J separate a singlet from a triplet and a singlet from a quintet, respectively. We draw two new insights: (a) the canonical tetracene singlet-fission unit cell supports precisely three low-lying TT intermediates with order-of-magnitude differences in J, and (b) the separable TT intermediate in carotenoids emanates from a pair of excitations to the second triplet state. Therefore, unlike with tetracenes, carotenoid fission requires above-gap excitations. In all cases, the distinguishability of the molecular triplets-that is, the extent of orbital overlap-determines the splitting within the spin manifold of TT states. Consequently, J represents a spectroscopic observable that distnguishes the resemblance between TT intermediates and the T...T product.

TheoDORE: A toolbox for a detailed and automated analysis of electronic excited state computations

The advent of ever more powerful excited-state electronic structure methods has led to a tremendous increase in the predictive power of computation, but it has also rendered the analysis of these computations much more challenging and time-consuming. TheoDORE tackles this problem through providing tools for post-processing excited-state computations, which automate repetitive tasks and provide rigorous and reproducible descriptors. Interfaces are available for ten different quantum chemistry codes and a range of excited-state methods implemented therein. This article provides an overview of three popular functionalities within TheoDORE, a fragment-based analysis for assigning state character, the computation of exciton sizes for measuring charge transfer, and the natural transition orbitals used not only for visualization but also for quantifying multiconfigurational character. Using the examples of an organic push-pull chromophore and a transition metal complex, it is shown how these tools can be used for a rigorous and automated assignment of excited-state character. In the case of a conjugated polymer, we venture beyond the limits of the traditional molecular orbital picture to uncover spatial correlation effects using electron-hole correlation plots and conditional densities.

Toward an understanding of electronic excitation energies beyond the molecular orbital picture

Can we gain an intuitive understanding of excitation energies beyond the molecular picture?

Spin-flip methods in quantum chemistry

This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.

Stimulated X-ray Raman and Absorption Spectroscopy of Iron-Sulfur Dimers

Iron-sulfur complexes play an important role in biological processes such as metabolic electron transport. A detailed understanding of the mechanism of long-range electron transfer requires knowledge of the electronic structure of the complexes, which has traditionally been challenging to obtain, either by theory or by experiment, but the situation has begun to change with advances in quantum chemical methods and intense free electron laser light sources. We compute the spectra for stimulated X-ray Raman spectroscopy (SXRS) and absorption spectroscopy of homovalent and mixed-valence [2Fe-2S] complexes, using the ab initio density matrix renormalization group algorithm. The simulated spectra show clear signatures of the theoretically predicted dense low-lying excited states within the d-d manifold. Furthermore, the difference in spectral intensity between the absorption-active and Raman-active states provides a potential mechanism to selectively excite states by a proper tuning of the excitation pump, to access the electronic dynamics within this manifold.

Quantitative El-Sayed Rules for Many-Body Wave Functions from Spinless Transition Density Matrices

One-particle transition density matrices and natural transition orbitals enable quantitative description of electronic transitions and interstate properties involving correlated many-body wave functions within the molecular orbital framework. Here we extend the formalism to the analysis of tensor properties, such as spin-orbit couplings (SOCs), which involve states of different spin projection. By using spinless density matrices and Wigner-Eckart's theorem, the approach allows one to treat the transitions between states with arbitrary spin projections in a uniform way. In addition to a pictorial representation of the transition, the analysis also yields quantitative contributions of hole-particle pairs into the overall many-body matrix elements. In particular, it helps to rationalize the magnitude of computed SOCs in terms of El-Sayed's rules. The capabilities of the new tool are illustrated by the analysis of the equation-of-motion coupled-cluster calculations of two transition metal complexes.

Linker-Dependent Singlet Fission in Tetracene Dimers

Separation of triplet excitons produced by singlet fission is crucial for efficient application of singlet fission materials. While earlier works explored the first step of singlet fission, the formation of the correlated triplet pair state, the focus of recent studies has been on understanding the second step of singlet fission, the formation of independent triplets from the correlated pair state. We present the synthesis and excited-state dynamics of meta- and para-bis(ethynyltetracenyl)benzene dimers that are analogues to the ortho-bis(ethynyltetracenyl)benzene dimer reported by our groups previously. A comparison of the excited-state properties of these dimers allows us to investigate the effects of electronic conjugation and coupling on singlet fission between the ethynyltetracene units within a dimer. In the para isomer, in which the two chromophores are conjugated, the singlet exciton yields the correlated triplet pair state, from which the triplet excitons can decouple via molecular rotations. In contrast, the meta isomer in which the two chromophores are cross-coupled predominantly relaxes via radiative decay. We also report the synthesis and excited-state dynamics of two para dimers with different bridging units joining the ethynyltetracenes. The rate of singlet fission is found to be faster in the dimer with the bridging unit that has orbitals closer in energy to that of the ethynyltetracene chromophores.

Theoretical Modeling of Singlet Fission

Singlet fission is a photophysical reaction in which a singlet excited electronic state splits into two spin-triplet states. Singlet fission was discovered more than 50 years ago, but the interest in this process has gained a lot of momentum in the past decade due to its potential as a way to boost solar cell efficiencies. This review presents and discusses the most recent advances with respect to the theoretical and computational studies on the singlet fission phenomenon. The work revisits important aspects regarding electronic states involved in the process, the evaluation of fission rates and interstate couplings, the study of the excited state dynamics in singlet fission, and the advances in the design and characterization of singlet fission compounds and materials such as molecular dimers, polymers, or extended structures. Finally, the review tries to pinpoint some aspects that need further improvement and proposes future lines of research for theoretical and computational chemists and physicists in order to further push the understanding and applicability of singlet fission.

Theoretical Studies of Singlet Fission: Searching for Materials and Exploring Mechanisms

In this Review article, a survey is given for theoretical studies in the subject of singlet fission. Singlet fission converts one singlet exciton to two triplet excitons. With the doubled number of excitons and the longer lifetime of the triplets, singlet fission provides an avenue to improve the photoelectric conversion efficiency in organic photovoltaic devices. It has been a subject of intense research in the past decade. Theoretical studies play an essential role in understanding singlet fission. This article presents a Review of theoretical studies in singlet fission since 2006, the year when the research interest in this subject was reignited. Both electronic structure and dynamics studies are covered. Electronic structure studies provide guidelines for designing singlet fission chromophores and insights into the couplings between single- and multi-excitonic states. The latter provides fundamental knowledge for engineering interchromophore conformations to enhance the fission efficiency. Dynamics studies reveal the importance of vibronic couplings in singlet fission.

Can TDDFT Describe Excited Electronic States of Naphthol Photoacids? A Closer Look with EOM-CCSD

The ¹L_b and ¹L_a excited states of naphthols are characterized by using time-dependent density functional theory (TDDFT), configuration interaction with singles (CIS), and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) methods. TDDFT fails dramatically at predicting the energy and ordering of the ¹L_a and ¹L_b excited states as observed experimentally, while EOM-CCSD accurately predicts the excited states as characterized by natural transition orbital analysis. The limitations of TDDFT are attributed to the absence of correlation from doubly excited configurations as well as the inconsistent description of excited electronic states of naphthol photoacids revealed by excitation analysis based on the one-electron transition density matrix.

Intermolecular Packing Effects on Singlet Fission in Oligorylene Dimers

Using the density functional theory method, the crystalline packing effect on the singlet fission (SF) rate of oligorylenes, some of which are found to exhibit SF in crystal forms, is revealed by evaluating the effective electronic coupling (|V_eff|), the square of which is proportional to the SF rate. The |V_eff| values for terrylene and quaterrylene dimer models are investigated for a variety of slip-stacked forms. It is found that these values show similar dependences on the intermolecular packing as a function of lateral and longitudinal displacements of monomer frameworks, and that they are maximized in several configurations of one monomer slipped from another along the longitudinal axis. The present estimation method of the SF rate is also found to qualitatively explain the experimental SF rate difference between terrylene derivatives with different packing forms. Furthermore, by analyzing the effect of electronic couplings on the adiabatic electronic states related to SF, we predict several favorable molecular packings leading to a fast SF with a high triplet yield.

Role of the Dark 2Ag State in Donor-Acceptor Copolymers as a Pathway for Singlet Fission: A DMRG Study

The mechanism of intramolecular singlet fission in donor-acceptor-type copolymers, especially the role of the dark 2A_g state, is not so clear. In this Letter, the electronic structure of the benzodithiophene (B)-thiophene-1,1-dioxide (TDO) copolymer is calculated by density matrix renormalization group theory with the Pariser-Parr-Pople model. We find that the dark 2A_g state is the lowest singlet excited state and is nearly degenerate with the 1B_u state. So, a fast internal conversion from 1B_u to 2A_g state is highly possible. The 2A_g state has a strong triplet pair character, localized on two neighboring acceptor units, which indicates that it is an intermediate state for the intramolecular singlet fission process. With the increase of the donor-acceptor push-pull strength in our model, this triplet pair character of the 2A_g state becomes more prominent, and meanwhile the binding energy of this coupled triplet pair state decreases, which favors the separation into two uncoupled triplet states. We propose a model in which the competition between the singlet fission process and the nonradiative decay process from the 2A_g state would determine the final quantum yield.