The role of the environment in electronic energy transfer: a molecular modeling perspective

Optical properties and exciton transfer between N-heterocyclic carbene iridium(III) complexes for blue light-emitting diode applications from first principles

N-heterocyclic carbene (NHC) iridium(III) complexes are considered as promising candidates for blue emitters in organic light-emitting diodes. They can play the roles of the emitter as well as of electron and hole transporters in the same emission layer. We investigate optical transitions in such complexes with account of geometry and electronic structure changes upon excitation or charging and exciton transfer between the complexes from first principles. It is shown that excitation of NHC iridium complexes is accompanied by a large reorganization energy ∼0.7 eV and a significant loss in the oscillator strength, which should lead to low exciton diffusion. Calculations with account of spin-orbit coupling reveal a small singlet-triplet splitting ∼0.1 eV, whereas the oscillator strength for triplet excitations is found to be an order of magnitude smaller than for the singlet ones. The contributions of the Förster and Dexter mechanisms are analyzed via the explicit integration of transition densities. It is shown that for typical distances between emitter complexes in the emission layer, the contribution of the Dexter mechanism should be negligible compared to the Förster mechanism. At the same time, the ideal dipole approximation, although giving the correct order of the exciton coupling, fails to reproduce the result taking into account spatial distribution of the transition density. For charged NHC complexes, we find a number of optical transitions close to the emission peak of the blue emitter with high exciton transfer rates that can be responsible for exciton-polaron quenching. The nature of these transitions is analyzed.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0-0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang-Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang-Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0-0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

Emissive and reactive excimers in a covalently-linked supramolecular multi-chromophoric system with a balanced rigid-flexible structure

A novel multi-chromophoric system, triad, in which two styrylbenzoquinoline (SBQ) photochromes are connected by a balanced rigid-flexible linker comprising 2,3-naphthylene framework (a residue of 3-oxy-2-naphthoic acid) and tetramethylene groups, was designed and synthesized to study an excimer formation in the excited state. The ¹H NMR data testified that triad exists in solution as folded conformers with asymmetric parallel-displaced SBQ units. Under light irradiation, in the triad, competitive photoisomerization and [2 + 2] photocycloaddition reactions were observed, both reactions being reversible. The photocycloaddition resulted in a tetrasubstituted cyclobutane. The red-shifted fluorescence spectrum and the appearance of a long-lived component in the triad fluorescence decay indicated formation of an 'emissive' excimer. The photocycloaddition is assumed to occur in a 'reactive' excimer, in which the ethylene groups of the SBQ photochromes are located at a distance sufficient for the formation of the σ-bonds between them. Quantum-chemical density functional theory (DFT) calculations at M06-2X/6-31G* level predicted the existence of the triad conformers with π-stacking interaction of SBQ photochromes, the structure of which is pre-organized for the excimer formation and photocycloaddition. For the first time, both emissive and reactive excimers were experimentally observed in the multi-chromophoric system with two diarylethylene photochromes undergoing [2 + 2] photocycloaddition.

Real-Time Propagation TDDFT and Density Analysis for Exciton Coupling Calculations in Large Systems

Photoactive systems are characterized by their capacity to absorb the energy of light and transform it. Usually, more than one chromophore is involved in the light absorption and excitation transport processes in complex systems. Linear-Response Time-Dependent Density Functional (LR-TDDFT) is commonly used to identify excitation energies and transition properties by solving the well-known Casida’s equation for single molecules. However, in practice, LR-TDDFT presents some disadvantages when dealing with multichromophore systems due to the increasing size of the electron–hole pairwise basis required for accurate evaluation of the absorption spectrum. In this work, we extend our local density decomposition method that enables us to disentangle individual contributions into the absorption spectrum to computation of exciton dynamic properties, such as exciton coupling parameters. We derive an analytical expression for the transition density from Real-Time Propagation TDDFT (P-TDDFT) based on Linear Response theorems. We demonstrate the validity of our method to determine transition dipole moments, transition densities, and exciton coupling for systems of increasing complexity. We start from the isolated benzaldehyde molecule, perform a distance analysis for π-stacked dimers, and finally map the exciton coupling for a 14 benzaldehyde cluster.

Is energy transfer limiting multiphotochromism? answers from ab initio quantifications

Dithienylethenes (DTEs) can be assembled to form supramolecular multiphotochromic systems that are highly functional molecular architectures of potential interest for building complex optoelectronic devices. Yet even simple DTE dimers relying on an organic linker may suffer from a partial photoactivity, i.e., only one of the two switches does isomerise. Experimentally, this limited photochromism has been attributed to an excited state energy transfer (EET) between the two DTEs of the multimer; this EET taking place instead of the desired photoinduced cyclisation of the DTE. However, no clear evidences of this phenomenon have been provided so far. In this work we propose the first rationalisation of this potential parasite photoinduced event using a computational approach based on Time-Dependent Density Functional Theory (TD-DFT) for the calculation of the electronic coupling in DTE dimers. Besides quantifying EET in several systems, we dissect the role of through-bond and through-space mechanisms on this process and clarify their dependence on both the nature and length of the bridge separating the two photochromes. The theoretical data obtained in this framework are in full agreement with the experimental outcomes and pave the way toward a molecular design of coupled, yet fully functionals, DTE-based multiswitches.

Oxidative species-induced excitonic transport in tubulin aromatic networks: Potential implications for neurodegenerative disease

Oxidative stress is a pathological hallmark of neurodegenerative tauopathic disorders such as Alzheimer's disease and Parkinson's disease-related dementia, which are characterized by altered forms of the microtubule-associated protein (MAP) tau. MAP tau is a key protein in stabilizing the microtubule architecture that regulates neuron morphology and synaptic strength. When MAP tau is degraded in tauopathic disorders, neuron dysfunction results. The precise role of reactive oxygen species (ROS) in the tauopathic disease process, however, is poorly understood. Classically, mitochondrial dysfunction has been viewed as the major source of oxidative stress and has been shown to precede tau and amyloid pathology in various dementias, but the exact mechanisms are not clear. It is known that the production of ROS by mitochondria can result in ultraweak photon emission (UPE) within cells. While of low intensity, surrounding proteins within the cytosol can still absorb these energetic photons via aromatic amino acids (e.g., tryptophan and tyrosine). One likely absorber of these photons is the microtubule cytoskeleton, as it forms a vast network spanning neurons, is highly co-localized with mitochondria, and shows a high density of aromatic amino acids. Functional microtubule networks may traffic this ROS-generated endogenous photon energy for cellular signaling, or they may serve as dissipaters/conduits of such energy to protect the cell from potentially harmful effects. Experimentally, after in vitro exposure to exogenous photons, microtubules have been shown to reorient and reorganize in a dose-dependent manner with the greatest effect being observed around 280nm, in the tryptophan and tyrosine absorption range. In this paper, recent modeling efforts based on ambient temperature experiment are presented, showing that tubulin polymers can feasibly absorb and channel these photoexcitations via resonance energy transfer, on the order of dendritic length scales and neuronal fine structure. Since microtubule networks are compromised in tauopathic diseases such as Alzheimer's and Parkinson's dementias, patients with these illnesses would be unable to support effective channeling of these photons for signaling or dissipation. Consequent emission surplus due to increased UPE production or decreased ability to absorb and transfer may lead to increased cellular oxidative damage, thus hastening the neurodegenerative process.

Relation between Vibrational Dephasing Time and Energy Gap Fluctuations

Dephasing processes are present in basically all applications in which quantum mechanics plays a role. These applications certainly include excitation energy and charge transfer in biological systems. In a previous study, we have analyzed the vibrational dephasing time as a function of energy gap fluctuation for a large set of molecular simulations. In that investigation, individual molecular subunits were the focus of the calculations. The set of studied molecules included bacteriochlorophylls in Fenna-Matthews-Olson and light-harvesting system 2 complexes as well as bilins in PE545 aggregates. The present work extends this study to entire complexes, including the respective intermolecular couplings. Again, it can be concluded that a universal and inverse proportionality exists between dephasing time and variance of the excitonic energy gap fluctuations, whereas the respective proportionality constants can be rationalized using the energy gap autocorrelation functions. Furthermore, these findings can be extended to the gaps between higher-lying neighboring excitonic states.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Localized operator partitioning method for electronic excitation energies in the time-dependent density functional formalism

We extend the localized operator partitioning method (LOPM) [J. Nagesh, A. F. Izmaylov, and P. Brumer, J. Chem. Phys. 142, 084114 (2015)] to the time-dependent density functional theory framework to partition molecular electronic energies of excited states in a rigorous manner. A molecular fragment is defined as a collection of atoms using Becke's atomic partitioning. A numerically efficient scheme for evaluating the fragment excitation energy is derived employing a resolution of the identity to preserve standard one- and two-electron integrals in the final expressions. The utility of this partitioning approach is demonstrated by examining several excited states of two bichromophoric compounds: 9-((1- naphthyl)- methyl)- anthracene and 4-((2- naphthyl)- methyl)- benzaldehyde. The LOPM is found to provide nontrivial insights into the nature of electronic energy localization that is not accessible using a simple density difference analysis.

Photoprotection and triplet energy transfer in higher plants: the role of electronic and nuclear fluctuations

Photosynthetic organisms employ several photoprotection strategies to avoid damage due to the excess energy in high light conditions.

A simple model for exploring the role of quantum coherence and the environment in excitonic energy transfer

We investigate the role of quantum coherence in modulating the energy transfer rate between two independent energy donors and a single acceptor participating in an excitonic energy transfer process. The energy transfer rate depends explicitly on the nature of the initial coherent superposition state of the two donors and we connect it to the observed absorption profile of the acceptor and the stimulated emission profile of the energy donors. We consider simple models with mesoscopic environments interacting with the donors and the acceptor and compare the expression we obtained for the energy transfer rate with the results of numerical integration.

A curve-crossing model to rationalize and optimize diarylethene dyads

Extra crossing points play a key role in the photochemistry of diarylethene dyads.

Going from photochromic compounds presenting a single switchable function to multi-addressable photochromic multimers remains an extremely difficult task notably because the interactions of several photochromic units through a linker generally result in a substantial loss of photoactivity. Due to their size and the intrinsic complexity of their electronic structure, coupled photochromes also constitute a fundamental challenge for theoretical chemistry. We present here an effective curve-crossing model that, used in connection with easily accessible ab initio data, allows a first understanding of the difficulty to obtain efficient multiphotochromes. Indeed, we demonstrate that extra crossing points, specific to multiphotochromes, have to be passed to ensure reactivity. In addition, the proposed approach allows the definition of an intuitive tilt criterion that can be used to screen a large number of substitution patterns and hence help in the design of new compounds, an aspect that is also developed here. The compatibility of this tilt criterion with previously proposed static Franck–Condon parameters is discussed as well.

The feasibility of coherent energy transfer in microtubules

It was once purported that biological systems were far too 'warm and wet' to support quantum phenomena mainly owing to thermal effects disrupting quantum coherence. However, recent experimental results and theoretical analyses have shown that thermal energy may assist, rather than disrupt, quantum coherent transport, especially in the 'dry' hydrophobic interiors of biomolecules. Specifically, evidence has been accumulating for the necessary involvement of quantum coherent energy transfer between uniquely arranged chromophores in light harvesting photosynthetic complexes. The 'tubulin' subunit proteins, which comprise microtubules, also possess a distinct architecture of chromophores, namely aromatic amino acids, including tryptophan. The geometry and dipolar properties of these aromatics are similar to those found in photosynthetic units indicating that tubulin may support coherent energy transfer. Tubulin aggregated into microtubule geometric lattices may support such energy transfer, which could be important for biological signalling and communication essential to living processes. Here, we perform a computational investigation of energy transfer between chromophoric amino acids in tubulin via dipole excitations coupled to the surrounding thermal environment. We present the spatial structure and energetic properties of the tryptophan residues in the microtubule constituent protein tubulin. Plausibility arguments for the conditions favouring a quantum mechanism of signal propagation along a microtubule are provided. Overall, we find that coherent energy transfer in tubulin and microtubules is biologically feasible.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. I. Full cumulant expansions and system-bath entanglement

We study the Förster resonant energy transfer rate, absorption and emission spectra in multichromophoric systems. The multichromophoric Förster theory (MCFT) is determined from an overlap integral of generalized matrices related to the donor's emission and acceptor's absorption spectra, which are obtained via a full 2nd-order cumulant expansion technique developed in this work. We calculate the spectra and MCFT rate for both localized and delocalized systems, and calibrate the analytical 2nd-order cumulant expansion with the exact stochastic path integral method. We present three essential findings: (i) The role of the initial entanglement between the donor and its bath is found to be crucial in both the emission spectrum and the MCFT rate. (ii) The absorption spectra obtained by the cumulant expansion method are nearly identical to the exact spectra for both localized and delocalized systems, even when the system-bath coupling is far from the perturbative regime. (iii) For the emission spectra, the cumulant expansion can give reliable results for localized systems, but fail to provide reliable spectra of the high-lying excited states of a delocalized system, when the system-bath coupling is large and the thermal energy is small. This paper also provides a simple golden-rule derivation of the MCFT, reviews existing methods, and motivates further developments in the subsequent papers.

Limits and potentials of quantum chemical methods in modelling photosynthetic antennae

A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.

Real-time UV-visible spectroscopy analysis of purple membrane-polyacrylamide film formation taking into account Fano line shapes and scattering

We theoretically and experimentally analyze the formation of thick Purple Membrane (PM) polyacrylamide (PA) films by means of optical spectroscopy by considering the absorption of bacteriorhodopsin and scattering. We have applied semiclassical quantum mechanical techniques for the calculation of absorption spectra by taking into account the Fano effects on the ground state of bacteriorhodopsin. A model of the formation of PM-polyacrylamide films has been proposed based on the growth of polymeric chains around purple membrane. Experimentally, the temporal evolution of the polymerization process of acrylamide has been studied as function of the pH solution, obtaining a good correspondence to the proposed model. Thus, due to the formation of intermediate bacteriorhodopsin-doped nanogel, by controlling the polymerization process, an alternative methodology for the synthesis of bacteriorhodopsin-doped nanogels can be provided.

The Fenna-Matthews-Olson protein revisited: a fully polarizable (TD)DFT/MM description

We report a combined molecular dynamics and quantum mechanics (QM)/molecular mechanics (MM) analysis of the excitonic properties of the Fenna-Matthews-Olson (FMO) protein by using a polarizable MM model combined with a time-dependent density functional theory description. Overall, our results indicate that structural fluctuations, electrostatic interactions, and short-range quantum effects can significantly modulate the model Hamiltonian parameters (site energies and couplings). We find that the specific interactions with the axial ligand and the hydrogen-bonded residues are responsible for the energy ladder, with their effects being mainly due to electrostatic interactions, but with short-range quantum contributions that are not negligible. In addition, a striking modulation of the screening effects experienced by the BChl pairs, due to the heterogeneous polarizability of the FMO and solvent environment, is observed. Finally, we find that the exciton model gives a reliable description of the delocalized excited states in the complex.

Understanding photosynthetic light-harvesting: a bottom up theoretical approach

We discuss a bottom up approach for modeling photosynthetic light-harvesting. Methods are reviewed for a full structure-based parameterization of the Hamiltonian of pigment-protein complexes (PPCs). These parameters comprise (i) the local transition energies of the pigments in their binding sites in the protein, the site energies; (ii) the couplings between optical transitions of the pigments, the excitonic couplings; and (iii) the spectral density characterizing the dynamic modulation of pigment transition energies and excitonic couplings by protein vibrations. Starting with quantum mechanics perturbation theory, we provide a microscopic foundation for the standard PPC Hamiltonian and relate the expressions obtained for its matrix elements to quantities that can be calculated with classical molecular mechanics/electrostatics approaches including the whole PPC in atomic detail and using charge and transition densities obtained with quantum chemical calculations on the isolated building blocks of the PPC. In the second part of this perspective, the Hamiltonian is utilized to describe the quantum dynamics of excitons. Situations are discussed that differ in the relative strength of excitonic and exciton-vibrational coupling. The predictive power of the approaches is demonstrated in application to different PPCs, and challenges for future work are outlined.

Theory of excitonic couplings in dielectric media : foundation of Poisson-TrEsp method and application to photosystem I trimers

The Poisson-TrEsp method (where TrEsp stands for transition charges from electrostatic potentials) has been successfully applied to calculate excitonic couplings in a variety of pigment-protein complexes. It relies on an isomorphism that allows for relating the excitonic coupling between transition densities in dielectric media to their Coulomb coupling. This isomorphism was derived by Hsu et al. (J. Chem. Phys. 114, 3065, (2001)) using time-dependent density functional response theory. In this article, we provide an alternative and simple derivation by first-order perturbation theory. An application of Poisson-TrEsp to photosystem I trimers reveals that the local field correction/screening factor depends on the mutual orientation of the pigments and on the dielectric boundaries rather than on distance. A mean correction factor of f = 0.69 is determined for this system.

Quantum chemical description of absorption properties and excited-state processes in photosynthetic systems

The theoretical description of the initial steps in photosynthesis has gained increasing importance over the past few years. This is caused by more and more structural data becoming available for light-harvesting complexes and reaction centers which form the basis for atomistic calculations and by the progress made in the development of first-principles methods for excited electronic states of large molecules. In this Review, we discuss the advantages and pitfalls of theoretical methods applicable to photosynthetic pigments. Besides methodological aspects of excited-state electronic-structure methods, studies on chlorophyll-type and carotenoid-like molecules are discussed. We also address the concepts of exciton coupling and excitation-energy transfer (EET) and compare the different theoretical methods for the calculation of EET coupling constants. Applications to photosynthetic light-harvesting complexes and reaction centers based on such models are also analyzed.

Excited-state calculations with TD-DFT: from benchmarks to simulations in complex environments

In this perspective, we present an overview of recent progress on Time-Dependent Density Functional Theory (TD-DFT) with a specific focus on its accuracy and on models able to take into account environmental effects, including complex media. To this end, we first summarise recent benchmarks and define an average TD-DFT accuracy in reproducing excitation energies when a conventional approach is used. Next, coupling of TD-DFT with models able to account for different kinds of interactions between a central chromophore and nearby chemical objects (solvent, organic cage, metal as well as semi-conducting surface) is investigated. Examples of application to excitation properties are presented, allowing to briefly describe several recent computational strategies. In addition, an extension of TD-DFT to describe a phenomenon involving interacting chromophores, e.g. the electronic energy transfer (EET), is presented to illustrate that this methodology can be applied to processes beyond the vertical excitation. This perspective therefore aims to provide to non-specialists a flavour of recent trends in the field of simulations of excited states in "realistic" situations.