Exciton–Phonon Interaction Model for Singlet Fission in Prototypical Molecular Crystals

Photo-induced dynamics with continuous and discrete quantum baths

The ultrafast quantum dynamics of photophysical processes in complex molecules is an extremely challenging computational problem with a broad variety of fascinating applications in quantum chemistry and biology. Inspired by recent developments in open quantum systems, we introduce a pure-state unraveled hybrid-bath method that describes a continuous environment via a set of discrete, effective bosonic degrees of freedom using a Markovian embedding. Our method is capable of describing both, a continuous spectral density and sharp peaks embedded into it. Thereby, we overcome the limitations of previous methods, which either capture long-time memory effects using the unitary dynamics of a set of discrete vibrational modes or use memoryless Markovian environments employing a Lindblad or Redfield master equation. We benchmark our method against two paradigmatic problems from quantum chemistry and biology. We demonstrate that compared to unitary descriptions, a significantly smaller number of bosonic modes suffices to describe the excitonic dynamics accurately, yielding a computational speed-up of nearly an order of magnitude. Furthermore, we take into account explicitly the effect of a δ-peak in the spectral density of a light-harvesting complex, demonstrating the strong impact of the long-time memory of the environment on the dynamics.

Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.

Theoretical investigation of distal charge separation in a perylenediimide trimer

An exciton-phonon (ex-ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex-ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore's reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.

Singlet fission is incoherent in pristine orthorhombic single crystals of rubrene: no evidence of triplet-pair emission

Singlet fission (SF) and its inverse, triplet-triplet annihilation (TTA), are promising strategies for enhancing photovoltaic efficiencies. However, detailed descriptions of the processes of SF/TTA are not fully understood, even in the most well-studied systems. Reports of the photophysics of crystalline rubrene, for example, are often inconsistent. Here we attempt to resolve these inconsistencies using time-resolved photoluminescence and transient absorption spectroscopy of 'pristine' rubrene orthorhombic single crystals. We find the reported time-resolved photoluminescence behaviour that hinted at triplet-pair emission is found only at specific sites on the crystals and likely arises from surface defects. Using transient absorption spectroscopy of the same crystals, we also observe no evidence of instantaneous generation of triplet-pair population with ∼100 fs excitation, independent of excitation wavelength (532 nm, 495 nm) or excitation angle. Our results suggest that SF occurs incoherently on a relatively slow (picosecond) timescale in rubrene single crystals, as expected from the original theoretical calculations. We conclude that the sub-100 fs formation of triplet pairs in crystalline rubrene films is likely to be due to static disorder.

Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions

Quantum dynamics (QD) simulation is a powerful tool for interpreting ultrafast spectroscopy experiments and unraveling their microscopic mechanism in out-of-equilibrium excited state behaviors in various chemical, biological, and material systems. Although state-of-the-art numerical QD approaches such as the time-dependent density matrix renormalization group (TD-DMRG) already greatly extended the solvable system size of general linearly coupled exciton-phonon models with up to a few hundred phonon modes, the accurate simulation of larger system sizes or strong system-environment interactions is still computationally highly challenging. Based on quantum information theory (QIT), in this work, we realize that only a small number of effective phonon modes couple to the excitonic system directly regardless of a large or even infinite number of modes in the condensed phase environment. On top of the identified small number of direct effective modes, we propose a hierarchical mapping (HM) approach through performing block Lanczos transformations on the remaining indirect modes, which transforms the Hamiltonian matrix to a nearly block-tridiagonal form and eliminates the long-range interactions. Numerical tests on model spin-boson systems and realistic singlet fission models in a rubrene crystal environment with up to 7000 modes and strong system-environment interactions indicate HM can reduce the system size by 1-2 orders of magnitude and accelerate the calculation by ∼80% without losing accuracy.

Plenty of Room on the Top: Pathways and Spectroscopic Signatures of Singlet Fission from Upper Singlet States

We investigate dynamic signatures of the singlet fission (SF) process triggered by the excitation of a molecular system to an upper singlet state S_N (N > 1) and develop a computational methodology for the simulation of nonlinear spectroscopic signals revealing the S_N → TT₁ SF in real time. We demonstrate that SF can proceed directly from the upper state S_N, bypassing the lowest excited state, S₁. We determine the main S_N → TT₁ reaction pathways and show by computer simulation and spectroscopic measurements that the S_N-initiated SF can be faster and more efficient than the traditionally studied S₁ → TT₁ SF. We claim that the S_N → TT₁ SF offers novel promising opportunities for engineering SF systems and enhancing SF yields.

Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials

The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH₂O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C_2nH_2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.

Dynamics of Excitons in Conjugated Molecules and Organic Semiconductor Systems

The exciton, an excited electron-hole pair bound by Coulomb attraction, plays a key role in photophysics of organic molecules and drives practically important phenomena such as photoinduced mechanical motions of a molecule, photochemical conversions, energy transfer, generation of free charge carriers, etc. Its behavior in extended π-conjugated molecules and disordered organic films is very different and very rich compared with exciton behavior in inorganic semiconductor crystals. Due to the high degree of variability of organic systems themselves, the exciton not only exerts changes on molecules that carry it but undergoes its own changes during all phases of its lifetime, that is, birth, conversion and transport, and decay. The goal of this review is to give a systematic and comprehensive view on exciton behavior in π-conjugated molecules and molecular assemblies at all phases of exciton evolution with emphasis on rates typical for this dynamic picture and various consequences of the above dynamics. To uncover the rich variety of exciton behavior, details of exciton formation, exciton transport, exciton energy conversion, direct and reverse intersystem crossing, and radiative and nonradiative decay are considered in different systems, where these processes lead to or are influenced by static and dynamic disorder, charge distribution symmetry breaking, photoinduced reactions, electron and proton transfer, structural rearrangements, exciton coupling with vibrations and intermediate particles, and exciton dissociation and annihilation as well.

Identifying Hidden Intracell Symmetries in Molecular Crystals and Their Impact for Multiexciton Generation

Organic molecular crystals are appealing for next-generation optoelectronic applications due to their multiexciton generation processes that can increase the efficiency of photovoltaic devices. However, a general understanding of how crystal structures affect these processes is lacking, requiring computationally demanding calculations for each material. Here we present an approach to understand and classify organic crystals and elucidate multiexciton processes. We show that organic crystals that are composed of two sublattices are well-approximated by effective fictitious systems of higher translational symmetry. Within this framework, we derive hidden selection rules in crystal pentacene and predict that the bulk polymorph supports fast Coulomb-mediated singlet fission with a transition rate about 2 orders of magnitude faster than that of the thin-film polymorph, a result confirmed with many-body perturbation theory calculations. Our approach is based on density-functional theory calculations and provides design principles for the experimental and computational discovery of new materials with tailored excitonic properties.

Evaluating the Electronic Structure of Coexisting Excitonic and Multiexcitonic States in Periodic Systems: Significance for Singlet Fission

Singlet fission (SF) in organic molecular solids is an example of a process that is challenging to describe with the most common electronic structure approaches. It involves optically bright singlet excited states delocalized over many molecules, which could be efficiently treated by density functional theory, and multiexcitonic localized states that have to be studied with wavefunction methods, usually with small clusters considering their expensive computational costs. In this work, we propose a methodology to combine multiconfigurational wavefunction calculations with reduced Hamiltonian to investigate the electronic structure of large clusters or fully periodic systems. The method is applied to the prototypical SF materials tetracene and pentacene. The results allow one to study how states of different natures (excitonic, charge-transfer, and multiexcitonic) coexist and are contaminated by their couplings in large or periodic systems. Novel insights are therefore possible. For example, because the excitonic bands are relatively broad with respect to the multiexcitonic states, there are limited regions of the crystal momentum space where the transition between the two is more likely.

Quantum dynamics simulation of intramolecular singlet fission in covalently linked tetracene dimer

In this work, we study singlet fission in tetracene para-dimers, covalently linked by a phenyl group. In contrast to most previous studies, we account for the full quantum dynamics of the combined excitonic and vibrational system. For our simulations, we choose a numerically unbiased representation of the molecule's wave function, enabling us to compare with experiments, exhibiting good agreement. Having access to the full wave function allows us to study in detail the post-quench dynamics of the excitons. Here, one of our main findings is the identification of a time scale t₀ ≈ 35 fs dominated by coherent dynamics. It is within this time scale that the larger fraction of the singlet fission yield is generated. We also report on a reduced number of phononic modes that play a crucial role in the energy transfer between excitonic and vibrational systems. Notably, the oscillation frequency of these modes coincides with the observed electronic coherence time t₀. We extend our investigations by also studying the dependency of the dynamics on the excitonic energy levels that, for instance, can be experimentally tuned by means of the solvent polarity. Here, our findings indicate that the singlet fission yield can be doubled, while the electronic coherence time t₀ is mainly unaffected.

Vacancy control in acene blends links exothermic singlet fission to coherence

The fission of singlet excitons into triplet pairs in organic materials holds great technological promise, but the rational application of this phenomenon is hampered by a lack of understanding of its complex photophysics. Here, we use the controlled introduction of vacancies by means of spacer molecules in tetracene and pentacene thin films as a tuning parameter complementing experimental observables to identify the operating principles of different singlet fission pathways. Time-resolved spectroscopic measurements in combination with microscopic modelling enables us to demonstrate distinct scenarios, resulting from different singlet-to-triplet pair energy alignments. For pentacene, where fission is exothermic, coherent mixing between the photoexcited singlet and triplet-pair states is promoted by vibronic resonances, which drives the fission process with little sensitivity to the vacancy concentration. Such vibronic resonances do not occur for endothermic materials such as tetracene, for which we find fission to be fully incoherent; a process that is shown to slow down with increasing vacancy concentration.

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.

Free-triplet generation with improved efficiency in tetracene oligomers through spatially separated triplet pair states

Singlet fission (SF) can potentially boost the efficiency of solar energy conversion by converting a singlet exciton (S₁) into two free triplets (T₁ + T₁) through an intermediate state of a correlated triplet pair (TT). Although efficient TT generation has been recently realized in many intramolecular SF materials, their potential applications have been hindered by the poor efficiency of TT dissociation. Here we demonstrate that this can be overcome by employing a spatially separated ¹(T…T) state with weak intertriplet coupling in tetracene oligomers with three or more chromophores. By using transient magneto-optical spectroscopic methods, we show that free-triplet generation can be markedly enhanced through the SF pathway that involves the spatially separated ¹(T…T) state rather than the pathway mediated by the spatially adjacent TT state, leading to a marked improvement in free-triplet generation with an efficiency increase from 21% for the dimer to 85% (124%) for the trimer (tetramer).

Bright Frenkel Excitons in Molecular Crystals: A Survey

We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models. Furthermore, this large set of data is used to search for technologically interesting and less common properties. For instance, we considered the largest excitonic bandwidth that is achievable within known molecular crystals and identified materials with strong super-radiant states. Finally, we explored the possibility that strong excitonic coupling can be used to generate emissive states in the near-infrared region in materials formed by molecules with bright visible absorption and we could identify the maximum allowable red shift in this material class. These insights with the associated searchable database provide practical guidelines for designing materials with interesting optical properties.

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.

Singlet Fission Dynamics in Tetracene Single Crystals Probed by Polarization-Dependent Two-Dimensional Electronic Spectroscopy

The exact mechanism of endothermic singlet fission in crystalline polyacene remains to be clarified. It has been elusive whether the excess energy of vibrational hot states and the upper branch of Davydov splitting is important for the energy compensation. Here, we probe the excited-state specified singlet fission dynamics in tetracene single crystals by polarization-dependent two-dimensional electronic spectroscopy (2DES). While a major spectral transfer with a characteristic lifetime of 86 ps is observed to be largely independent of the excitation energy due to formation of the spatially separated triplet pairs (¹(T···T)), the excitation-energy dependent subpicosecond dynamics show marked differences for different states probed, implying the possible involvement of a coherently formed triplet pair state (¹(TT)). Analysis of coherent vibrational modes suggests the coupling to high energy modes may offset the energy difference between singlet and triplet pair states. Moreover, the beating map of the low frequency mode indicates a vibrational hot state violating the aggregation behavior of Davydov exciton, which can be explained as a resonance of the ¹(TT) state. These results suggest that the coherent vibronic mixing between local excitation and triplet pair states is essential for the singlet fission dynamics in molecule aggregates.

FB-REDA: fragment-based decomposition analysis of the reorganization energy for organic semiconductors

We present a fragment-based decomposition analysis tool (FB-REDA) for the reorganisation energy (λ). This tool delivers insights on how to rationally design low-λ organic semiconductors. The contribution of the fragment vibrational modes to the reorganization energy is exploited to identity the individual contributions of the molecular building blocks. The usefulness of the approach is demonstrated by offering three strategies to reduce the reorganization energy of a promising dopant-free hole transport material (TPA1PM, λ = 213 meV). A reduction of nearly 50% (TPD3PM, λ = 108 meV) is achieved. The proposed design principles are likely transferable to other organic semiconductors exploiting common molecular building blocks.