Importance and Nature of Short-Range Excitonic Interactions in Light Harvesting Complexes and Organic Semiconductors

Computational Design of (B)Chl Models: Structural and Chemical Modifications toward Enriched Properties

The functional units of natural photosynthetic systems control the process of converting sunlight into chemical energy. In this article, we explore a series of chemically and structurally modified bacteriochlorophyll and chlorophyll pigments through computational chemistry to evaluate their electronic spectroscopy properties. More specifically, we use multiconfigurational and time-dependent density functional theory methods, along with molecular dynamics simulations, to compute the models' energetics both in an implicit and explicit solvent environment. Structural modifications aimed at reducing the planarity of the macrocycle through alkyl-bridge anchoring reveal the significant role of the curvature in fine-tuning spectral properties, which mimics protein scaffold effects on naturally occurring pigments. Furthermore, chemical substitutions with a carbonyl group show potential for expanding absorption spectra toward the blue region, while incorporating an additional double bond decreases absorption efficiency. These insights lay the groundwork to design novel synthetic pigments, with potential applications in artificial light-harvesting systems and more efficient photovoltaic devices.

Analysis of Site Energies and Excitonic Couplings: The Role of Symmetry and Polarization

An excitonic coupling model is developed based on an equation-of-motion coupled cluster combined with the fragment molecular orbital method. The effects of polarization and excitonic coupling on the splitting of quasi-degenerate levels in systems containing multiple chromophores are elucidated on dimers of formaldehyde, water, formic acid, hydrogen fluoride, and carbon monoxide. It is shown that the level structure is mainly determined by the mutual polarization of chromophores and to a lesser extent by the excitonic coupling. The role of symmetry in excitonic coupling in dimers is discussed. The excitonic coupling between all residues in the photoactive yellow protein (PDB: 2PHY) is analyzed.

Geometry-Independent Ultrafast Energy Transfer in Bioinspired Arrays Containing Electronically Coupled BODIPY Dimers as Energy Donors

In photosynthetic light-harvesting complexes, strong interaction between chromophores enables efficient absorption of solar radiation and has been suggested to enable ultrafast energy funneling to the reaction center. To examine whether similar effects can be realized in synthetic systems, and to determine the mechanisms of energy transfer, we synthesized and characterized a series of bioinspired arrays containing strongly-coupled BODIPY dimers as energy donors and chlorin derivatives as energy acceptors. The BODIPY dimers feature broad absorption in the range of 500-600 nm, complementing the chlorin absorption to provide absorption across the entire visible spectrum. Ultrafast (~10 ps) energy transfer was observed from photoexcited BODIPY dyads to chlorin subunits. Surprisingly, the energy-transfer rate is nearly independent of the position where the BODIPY dimer is attached to the chlorin and of the type of connecting linker. In addition, the energy-transfer rate from BODIPY dimers to chlorin is slower than the corresponding rate in arrays containing BODIPY monomers. The lower rate, corresponding to less efficient through-bond transfer, is most likely due to weaker electronic coupling between the ground state of the chlorin acceptor and the delocalized electronic state of the BODIPY dimer, compared to the localized state of a BODIPY monomer.

Simultaneous Force and Fluorescence Spectroscopy on Single Chains of Polyfluorene: Effect of Intra-Chain Aggregate Coupling

Conjugated polymer chains in compact conformations or in films exhibit spectral features that can be attributed to interactions between individual conjugated segments of the chain, including formation of aggregates or excimers. Here, we use atomic force microscopy (AFM) on single chains of the conjugated polymer polyfluorene (PFO) to control the intersegment interactions by mechanically unfolding the chain. Simultaneously with the force spectroscopy we monitor fluorescence from the single PFO chains using a fluorescence microscope. We found that mechanical stretching of the chain causes disappearance of the green emission band. This observation provides evidence that the green emission originates from an intrachain aggregated state on the self-folded chain, which is decoupled by the stretching. In addition, the stretching upon laser irradiation leads to the appearance of additional features in the force spectra, small force peaks in the initial stages of the unfolding. These features are attributed to a combination of excitonic and van der Waals coupling of a ground-state intrachain aggregate.

Photophysical Properties and Electronic Structure of Hydroporphyrin Dyads Exhibiting Strong Through-Space and Through-Bond Electronic Interactions

Electronic interactions between tetrapyrroles are utilized in natural photosynthetic systems to tune the light-harvesting and energy-/charge-transfer processes in these assemblies. Such interactions also can be employed to tailor the electronic properties of tetrapyrrolic dyads and larger arrays for use in materials science and biomedical research. Here, we have utilized static and time-resolved optical spectroscopy to characterize the optical absorption and emission properties of a set of chlorin and bacteriochlorin dyads with varying degrees of through-bond (TB) and through-space (TS) interactions between the constituent macrocycles. The dyads consist of two chlorins or two bacteriochlorins joined by a linker that utilizes a triple-double-triple-bond (enediyne) motif in which the double-bond portion is an ester-substituted ethylene or o-phenylene unit. The photophysical studies are coupled with density functional theory (DFT) calculations to probe the ground-state molecular orbital (MO) characteristics of the dyads and time-dependent DFT calculations (TDDFT) to elucidate excited-state properties. The latter include electronic characteristics of the singlet excited-state manifold and the absorption transitions to these states from the electronic ground state. A comparison of the MO and calculated spectral properties of each dyad with the linker present versus disrupted (by eliminating the double-bond portion) gives insight into the relative contributions of TB versus TS interactions to the electronic properties of the dyads. The results show that the TB and TS contributions are additive (constructively interfere), which is not always the case for molecular dyads. Most of the dyads have shorter lifetimes of the lowest singlet excited state compared to the parent monomer, which derives from increased S₁ → S₀ internal conversion. The enhancement is greater for the dyads in benzonitrile than in toluene. The studies provide insights into the nature of the electronic interactions between the constituents in the tetrapyrrole arrays and how these interactions dictate the spectral properties and excited-state decay characteristics.

Accurate Modeling of Excitonic Coupling in Cyanine Dye Cy3

Accurate modeling of excitonic coupling in molecules is of great importance for inferring the structures and dynamics of coupled systems. Cy3 is a cyanine dye that is widely used in molecular spectroscopy. Its well-separated excitation bands, high sensitivity to the surroundings, and the high energy transfer efficiency make it a perfect choice for excitonic coupling experiments. Many methods have been used to model the excitonic coupling in molecules with varying degrees of accuracy. The atomic transition charge model offers a high-accuracy and cost-effective way to calculating the excitonic coupling. The main focus of this work is to generate high-quality atomic transition charges that can accurately model the Cy3 dye's transition density. The transition density of the excitation of the ground to first excited state is calculated using configuration-interaction singles and time-dependent density functional theory and is benchmarked against the algebraic diagrammatic construction method. Using the transition density we derived the atomic transition charges using two approaches: Mulliken population analysis and charges fitted to the transition electrostatic potential. The quality of the charges is examined, and their ability to accurately calculate the excitonic coupling is assessed via comparison to experimental data of an artificial biscyanine construct. Theoretical comparisons to the supermolecule ab initio couplings and the widely used point-dipole approximation are also made. Results show that using the transition electrostatic potential is a reliable approach for generating the transition atomic charges. A high-quality set of charges, that can be used to model the Cy3 dye dimer excitonic coupling with high-accuracy and a reasonable computational cost, is obtained.

n-Type organic semiconducting polymers: stability limitations, design considerations and applications

This review outlines the design strategies which aim to develop high performing n-type materials in the fields of organic thin film transistors (OTFT), organic electrochemical transistors (OECT) and organic thermoelectrics (OTE). Figures of merit for each application and the limitations in obtaining these are set out, and the challenges with achieving consistent and comparable measurements are addressed. We present a thorough discussion of the limitations of n-type materials, particularly their ambient operational instability, and suggest synthetic methods to overcome these. This instability originates from the oxidation of the negative polaron of the organic semiconductor (OSC) by water and oxygen, the potentials of which commonly fall within the electrochemical window of n-type OSCs, and consequently require a LUMO level deeper than ∼-4 eV for a material with ambient stability. Recent high performing n-type materials are detailed for each application and their design principles are discussed to explain how synthetic modifications can enhance performance. This can be achieved through a number of strategies, including utilising an electron deficient acceptor-acceptor backbone repeat unit motif, introducing electron-withdrawing groups or heteroatoms, rigidification and planarisation of the polymer backbone and through increasing the conjugation length. By studying the fundamental synthetic design principles which have been employed to date, this review highlights a path to the development of promising polymers for n-type OSC applications in the future.

Geometry dependence of excitonic couplings and the consequences for configuration-space sampling

Excitonic coupling plays a key role for the understanding of excitonic energy transport (EET) in, for example, organic photovoltaics. However, the calculation of realistic systems is often beyond the applicability range of accurate wavefunction methods so that lower-scaling semi-empirical methods are used to model EET events. In the present work, the distance and angle dependence of excitonic couplings of dimers of selected organic molecules are evaluated for the semi-empirical long-range corrected density functional based tight binding (LC-DFTB) method and spin opposite scaled second order approximate coupled cluster singles and doubles (SOS-CC2). While semi-empirically scaled methods can lead to slightly increased deviations for excitation energies, the excitonic couplings and their dependence on the dimer geometry are reproduced. LC-DFTB yields a similar accuracy range as density-functional theory (DFT) employing the ωB97X functional while the computation time is reduced by several orders of magnitude. The dependence of the exchange contributions to the excitonic couplings on the dimer geometry is analyzed assessing the calculation of Coulombic excitonic couplings from monomer local excited states only, which reduces the computational effort significantly. The present work is a necessary first step toward the simulation of excitonic energy transport using semi-empirical methods.

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.

Role of Conical Intersections on the Efficiency of Fluorescent Organic Molecular Crystals

Organic molecular crystals are attractive materials for luminescent applications because of their promised tunability. However, the link between the chemical structure and emissive behavior is poorly understood because of the numerous interconnected factors which are at play in determining radiative and nonradiative behaviors at the solid-state level. In particular, the decay through conical intersection dominates the nonadiabatic regions of the potential energy surface, and thus, their accessibility is a telling indicator of the luminosity of the material. In this study, we investigate the radiative mechanism for five organic molecular crystals which display a solid-state emission, with a focus on the role of conical intersections in their photomechanisms. The objective is to situate the importance of the accessibility of conical intersections with regards to emissive behavior, taking into account other nonradiative decay channels, namely, vibrational decay, and exciton hopping. We begin by giving a brief overview of the structural patterns of the five systems within a larger pool of 13 crystals for a richer comparison. We observe that because of the prevalence of sheet like and herringbone packing in organic molecular crystals, the conformational diversity of crystal dimers is limited. Additionally, similarly spaced dimers have exciton coupling values of a similar order within a 50 meV interval. Next, we focus on three exemplary cases, where we disentangle the role of nonradiative decay mechanisms and show how rotational minimum energy conical intersections in vacuum lead to puckered ones in the crystal, increasing their instability upon crystallization in typical packing motifs. In contrast, molecules with puckered conical intersections in vacuum tend to conserve this trait upon crystallization, and therefore, their quantum yield of fluorescence is determined predominantly by other nonradiative decay mechanisms.

FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method

In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.

fromage: A library for the study of molecular crystal excited states at the aggregate scale

The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in‐house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready‐to‐use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.

On the arrangement of chromophores in light harvesting complexes: chance versus design

We used a homogeneous computational approach to derive the excitonic Hamiltonian for five light harvesting complexes containing only one type of chromophore and compare them in terms of statistical descriptors. We then studied the approximate exciton dynamics for the five complexes introducing a measure, the (averaged and time-dependent) inverse participation ratio, that enables the comparison between very diverse complexes on the same ground. We find that the global dynamics are very similar across the set of systems despite the variety of geometric structures of the complexes. In particular, the dynamics of four out of five light harvesting complexes are barely distinguishable with a small variation from the norm seen only for the Fenna-Matthews-Olson complex. We use the information from the realistic Hamiltonians to build a reduced model system that shows how the global dynamics are ultimately dominated by a single parameter, the degree of localization of the excitonic Hamiltonian eigenstates. Considering the physically plausible range of system parameters, the reduced model explains why the dynamics are so similar across most light harvesting complexes containing a single type of chromophore regardless of the detailed pattern of the inter-chromophore excitonic coupling.

Photoisomerization of Enediynyl Linker Leads to Slipped Cofacial Hydroporphyrin Dyads with Strong Through-Bond and Through-Space Electronic Interactions

Photoisomerization of 3,4-di(methoxycarbonyl)-enediyne linker in hydroporphyrin (chlorin or bacteriochlorin) dyads leads to thermally stable cis isomers, where macrocycles adopt a slipped cofacial mutual geometry with an edge-to-edge distance of ∼3.6 Å (determined by density functional theory (DFT) calculations). Absorption spectra exhibit a significant splitting of the long-wavelength Q_y band, which indicates a strong electronic coupling with a strength of V = ∼477 cm^-1 that increases to 725 cm^-1 upon metalation of hydroporphyrins. Each dyad features a broad, structureless emission band, with large Stokes shift, which is indicative of excimer formation. DFT calculations for dyads show both strong through-bond electronic coupling and through-space electronic interactions, due to the overlap of π-orbitals. Overall, geometry, electronic structure, strength of electronic interactions, and optical properties of reported dyads closely resemble those observed for photosynthetic special pairs. Dyads reported here represent a novel type of photoactive arrays with various modes of electronic interactions between chromophores. Combining through-bond and through-space coupling appears to be a viable strategy to engineer novel optical and photochemical properties in organic conjugated materials.

Discrete π-Stacks of Perylene Bisimide Dyes within Folda-Dimers: Insight into Long- and Short-Range Exciton Coupling

Four well-defined π-stacks of perylene bisimide (PBI) dyes were obtained in solution by covalent linkage of two chromophores with spacer units of different length and sterical demand. Structural elucidation of the folda-dimers by in-depth nuclear magnetic resonance studies and geometry optimization at the level of density functional theory suggest different, but highly defined molecular arrangements of the two chromophores in the folded state enforced by the various spacer moieties. Remarkably, the dye stacks exhibit considerably different optical properties as investigated by UV/vis absorption and fluorescence spectroscopy, despite only slightly different chromophore arrangements. The distinct absorption properties can be rationalized by an interplay of long- and short-range exciton coupling resulting in optical signatures ranging from conventional H-type to monomer like absorption features with low and appreciably high fluorescence quantum yields, respectively. To the best of our knowledge, we present the first experimental proof of a PBI-based "null-aggregate", in which long- and short-range exciton coupling fully compensate each other, giving rise to monomer-like absorption features for a stack of two PBI chromophores. Hence, our insights pinpoint the importance of charge-transfer mediated short-range coupling that can significantly influence the optical properties of PBI π-stacks.

How fine-tuned for energy transfer is the environmental noise produced by proteins around biological chromophores?

Analysis of intermolecular motions of pigment–protein complexes shows no significant difference in influence of local environment despite different biological functions.

How Inter- and Intramolecular Processes Dictate Aggregation-Induced Emission in Crystals Undergoing Excited-State Proton Transfer

Aggregation-induced emission (AIE) offers a route for the development of luminescent technologies with high quantum efficiencies. Excited-state intramolecular proton transfer (ESIPT) coupled to AIE can produce devices with emission across the visible spectrum. We use a combination of theoretical models to determine the factors that mediate fluorescence in molecular crystals undergoing ESIPT. Using two materials based on 2'-hydroxychalcone as exemplar cases, we analyze how inter- and intramolecular processes determine the emissive properties in the crystal environment. This systematic investigation extends the current interpretation of AIE to polar chromophores with multiple decay pathways. We find that population of nonradiative pathways is dictated by the electronic effects of the substituents and the degree of distortion allowed in the crystal environment. Localization of the electron density is crucial to maximize fluorescence via ESIPT. Our conclusions offer design strategies for the development of luminescent molecular crystals.

Excitonic Coupling Modulated by Mechanical Stimuli

Understanding energy transfer is of vital importance in a diverse range of applications from biological systems to photovoltaics. The ability to tune excitonic coupling in any of these systems, however, is generally limited. In this work, we have simulated a new class of single-molecule spectroscopy in which force microscopy is used to control the excitonic coupling between chromophores. Here we demonstrate that the excitonic coupling can be controlled by mechanical manipulation of the molecule (perylenediimide dimers and terrylenediimide-perylenediimide heterodimers) and can be tuned over a broad range of values (0.02-0.15 eV) that correspond to different regimes of exciton dynamics going from the folded to the elongated structure of the dimer. In all of the systems considered here, the switching from high to low coupling takes place simultaneously with the mechanical deformation detected by a strong increase and subsequent decay of the force. These simulations suggest that single-molecule force spectroscopy can be used to understand and eventually aid the design of excitonic devices.