Theoretical investigation of electronic excitation energy transfer in bichromophoric assemblies

Panchromatic Light-Capturing Bis-styryl BODIPY-Perylenediimide Donor-Acceptor Constructs: Occurrence of Sequential Energy Transfer Followed by Electron Transfer

Two wide-band-capturing donor-acceptor conjugates featuring bis-styrylBODIPY and perylenediimide (PDI) have been newly synthesized, and the occurrence of ultrafast excitation transfer from the 1 PDI* to BODIPY, and a subsequent electron transfer from the 1 BODIPY* to PDI have been demonstrated. Optical absorption studies revealed panchromatic light capture but offered no evidence of ground-state interactions between the donor and acceptor entities. Steady-state fluorescence and excitation spectral recordings provided evidence of singlet-singlet energy transfer in these dyads, and quenched fluorescence of bis-styrylBODIPY emission in the dyads suggested additional photo-events. The facile oxidation of bis-styrylBODIPY and facile reduction of PDI, establishing their relative roles of electron donor and acceptor, were borne out by electrochemical studies. The electrostatic potential surfaces of the S1 and S2 states, derived from time-dependent DFT calculations, supported excited charge transfer in these dyads. Spectro-electrochemical studies on one-electron-oxidized and one-electron-reduced dyads and the monomeric precursor compounds were also performed in a thin-layer optical cell under corresponding applied potentials. From this study, both bis-styrylBODIPY⋅+ and PDI⋅- could be spectrally characterizes and were subsequently used in characterizing the electron-transfer products. Finally, pump-probe spectral studies were performed in dichlorobenzene under selective PDI and bis-styrylBODIPY excitation to secure energy and electron-transfer evidence. The measured rate constants for energy transfer, kENT , were in the range of 1011  s-1 , while the electron transfer rate constants, kET , were in the range of 1010  s-1 , thus highlighting their potential use in solar energy harvesting and optoelectronic applications.

Benchmarking Aspects of Ab Initio Fragment Models for Accurate Excimer Potential Energy Surfaces

While Coupled-Cluster methods have been proven to provide an accurate description of excited electronic states, the scaling of the computational costs with the system size limits the degree for which these methods can be applied. In this work different aspects of fragment-based approaches are studied on noncovalently bound molecular complexes with interacting chromophores of the fragments, such as π-stacked nucleobases. The interaction of the fragments is considered at two distinct steps. First, the states localized on the fragments are described in the presence of the other fragment(s); for this we test two approaches. One method is founded on QM/MM principles, only including the electrostatic interaction between the fragments in the electronic structure calculation with Pauli repulsion and dispersion effects added separately. The other model, a Projection-based Embedding (PbE) using the Huzinaga equation, includes both electrostatic and Pauli repulsion and only needs to be augmented by dispersion interactions. In both schemes the extended Effective Fragment Potential (EFP2) method of Gordon et al. was found to provide an adequate correction for the missing terms. In the second step, the interaction of the localized chromophores is modeled for a proper description of the excitonic coupling. Here the inclusion of purely electrostatic contributions appears to be sufficient: it is found that the Coulomb part of the coupling provides accurate splitting of the energies of interacting chromophores that are separated by more than 4 Å.

A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates

We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.

Electronic Couplings and Electrostatic Interactions Behind the Light Absorption of Retinal Proteins

The photo-functional chromophore retinal exhibits a wide variety of optical absorption properties depending on its intermolecular interactions with surrounding proteins and other chromophores. By utilizing these properties, microbial and animal rhodopsins express biological functions such as ion-transport and signal transduction. In this review, we present the molecular mechanisms underlying light absorption in rhodopsins, as revealed by quantum chemical calculations. Here, symmetry-adapted cluster-configuration interaction (SAC-CI), combined quantum mechanical and molecular mechanical (QM/MM), and transition-density-fragment interaction (TDFI) methods are used to describe the electronic structure of the retinal, the surrounding protein environment, and the electronic coupling between chromophores, respectively. These computational approaches provide successful reproductions of experimentally observed absorption and circular dichroism (CD) spectra, as well as insights into the mechanisms of unique optical properties in terms of chromophore-protein electrostatic interactions and chromophore-chromophore electronic couplings. On the basis of the molecular mechanisms revealed in these studies, we also discuss strategies for artificial design of the optical absorption properties of rhodopsins.

Spectral Tuning Mechanism of Photosynthetic Light-Harvesting Complex II Revealed by Ab Initio Dimer Exciton Model

Excited states of two kinds of bacteriochlorophyll (BChl) aggregates, B850 and B800, in photosynthetic light-harvesting complex II (LH2) are theoretically investigated by developing and using an extended exciton model considering efficiently evaluated excitonic coupling. Our exciton model based on dimer fragmentation is shown to reproduce the experimental absorption spectrum of LH2 with good accuracy, entailing their different redshifts originating from aggregations of B850 and B800. The systematic analysis has been performed on the spectra by quantitatively decomposing their spectral shift energies into the contributions of various effects: structural distortion, electrostatic, excitonic coupling, and charge-transfer (CT) effects. Our results show that the spectral redshift of B800 is mainly attributed to its electrostatic interaction with the protein environment, while that of B850 arises from the marked effect of the excitonic coupling between BChl units. The interchromophore CT excitation also plays a key role in the spectral redshift of B850. This CT effect can be effectively described using our dimer model. This suited characterization reveals that the pronounced CT effect originates from the characteristics of B850 that has closely spaced BChls as dimers. We highlight the importance of the refinement of the crystal structure with the use of quantum chemical methods for prediction of the spectrum.

Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead?

A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.

Theoretical investigation of non-Förster exciton transfer mechanisms in perylene diimide donor, phenylene bridge, and terrylene diimide acceptor systems

The rates of exciton transfer within dyads of perylene diimide and terrylene diimide connected by oligophenylene bridge units have been shown to deviate significantly from those of Förster's resonance energy transfer theory, according to single molecule spectroscopy experiments. The present work provides a detailed computational and theoretical study investigating the source of such a discrepancy. Electronic spectroscopy data are calculated by time-dependent density functional theory and then compared with experimental results. Electronic couplings between the exciton donor and the acceptor are estimated based on both the transition density cube method and transition dipole approximation. These results confirm that the delocalization of the exciton to the bridge parts contributes to significant enhancement of donor-acceptor electronic coupling. Mechanistic details of exciton transfer are examined by estimating the contributions of the bridge electronic states, vibrational modes of the dyads commonly coupled to both donor and acceptor, inelastic resonance energy transfer mechanism, and dark exciton states. These analyses suggest that the contribution of common vibrational modes serves as the main source of deviation from Förster's spectral overlap expression.

Energy Transfer from Perovskite Nanocrystals to Dye Molecules Does Not Occur by FRET

Single formamidinium lead bromide (FAPbBr₃) perovskite nanocubes, approximately 10 nm in size, have extinction cross sections orders of magnitude larger than single dye molecules and can therefore be used to photoexcite one single dye molecule within their immediate vicinity by means of excitation-energy transfer (EET). The rate of photon emission by the single dye molecule is increased by 2 orders of magnitude under excitation by EET compared to direct excitation at the same laser fluence. Because the dye cannot accommodate biexcitons, NC biexcitons are filtered out by EET, giving rise to up to an order-of-magnitude improvement in the fidelity of photon antibunching. We demonstrate here that, contrary to expectation, energy transfer from the nanocrystal to dye molecules does not depend on the spectral line widths of the donor and acceptor and is therefore not governed by Förster's theory of resonance energy transfer (FRET). Two different cyanine dye acceptors with substantially different spectral overlaps with the nanocrystal donor show a similar light-harvesting capability. Cooling the sample from room temperature to 5 K reduces the average transition line widths 25-fold but has no apparent effect on the number of molecules emitting, i.e., on the spatial density of single dye molecules being photoexcited by single nanocrystals. Narrow zero-phonon lines are identified for both donor and acceptor, with an energetic separation of over 40 times the line width, implying a complete absence of spectral overlap-even though EET is evident. Both donor and acceptor exhibit spectral fluctuations, but no correlation is apparent between the jitter, which controls spectral overlap, and the overall light harvesting. We conclude that the energy transfer process is fundamentally nonresonant, implying effective energy dissipation in the perovskite donor because of strong electron-phonon coupling of the carriers comprising the exciton. The work highlights the importance of performing cryogenic spectroscopy to reveal the underlying mechanisms of energy transfer in complex donor-acceptor systems.

Robust and Fragile Quantum Effects in the Transfer Kinetics of Delocalized Excitons between B850 Units of LH2 Complexes

Aggregates of light harvesting 2 (LH2) complexes form the major exciton-relaying domain in the photosynthetic unit of purple bacteria. Application of a generalized master equation to pairs of the B850 units of LH2 complexes, where excitons predominantly reside, provides quantitative information on how the inter-LH2 exciton transfer depends on the distance, relative rotational angle, and the relative energies of the two LH2s. The distance dependence demonstrates significant enhancement of the rate due to quantum delocalization of excitons, the qualitative nature of which remains robust against the disorder. The angle dependence reflects isotropic nature of exciton transfer, which remains similar for the ensemble of disorder. The variation of the rate on relative excitation energies of LH2 exhibits resonance peaks, which, however, is fragile as the disorder becomes significant. Overall, the average transfer times between two LH2s are estimated to be in the range of 4-25 ps for physically plausible inter-LH2 distances.

All-DFTB Approach to the Parametrization of the System-Bath Hamiltonian Describing Exciton-Vibrational Dynamics of Molecular Assemblies

Spectral density functions are central to the simulation of complex many body systems. Their determination requires making approximations not only to the dynamics but also to the underlying electronic structure theory. Here, blending different methods bears the danger of an inconsistent description. To solve this issue we propose an all-DFTB approach to determine spectral densities for the description of Frenkel excitons in molecular assemblies. The protocol is illustrated for a model of a PTCDI crystal, which involves the calculation of monomeric excitation energies and Coulomb couplings between monomer transitions, as well as their spectral distributions due to thermal fluctuations of the nuclei. Using dynamically defined normal modes, a mapping onto the standard harmonic oscillator spectral densities is achieved.

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.

The electronic character of PTCDA thin films in comparison to other perylene-based organic semi-conductors: ab initio-, TD-DFT and semi-empirical computations of the opto-electronic properties of large aggregates

Perylene-based compounds are promising materials for opto-electronic thin film devices but despite intense investigations, important details of their electronic structure are still under debate. For perylene-3,4,9,10-tetracarboxylic dianhydrid (PTCDA), the theoretical models predict a different relative energetic order of Frenkel and Charge Transfer (CT) states. Additionally, while one model indicates strong differences between PTCDA on one hand and other perylene-based compounds on the other, recent ab initio computations indicate electronic properties of all perylene-based compounds to resemble each other. Finally, the models disagree about the amount of mixing between CT and Frenkel states. Definitive answers to these questions are difficult because the approaches use various approximations. Up to date, the ab initio based methods employ rather small model systems and neglect environmental effects. In the present work, we improve our former approach by analyzing the effects of the various simplifications. In more detail, we increase the size of the model systems, include environmental effects and investigate the influence of exciton-phonon couplings on the absorption spectrum. The computations for larger aggregates were performed with the ZINDO/S approach, because benchmark computations show that it provides accurate vertical excitation energies for Frenkel, as well as CT states.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Dependence of FRET efficiency on distance in single donor-acceptor pairs

Possibility to create single donor-acceptor (D-A) pairs by attaching dye molecules to various sites of DNA strands with control of the inter-dye distance R enables one to measure average Förster resonance energy transfer (FRET) efficiency E as a function of R. Triplet states of the dyes influence the dependence E(R) considerably. Two types of FRET efficiency are considered: E = EA and E = ED. The efficiency EA(R) = JA(R)/[JA(R) + JD(R)] depends on the donor and the acceptor average intensities JD(R) and JA(R) measured in D- and A-fluorescence, whereas the efficiency ED(R) = 1 - JD(R)/JD(∞) depends only on the intensity of D-fluorescence, so-called the donor quenching method. The shape of the functions ED (R) and EA (R) depends strongly on whether the dyes have blinking fluorescence. FRET efficiencies ED (R) and EA (R) undergo the influence of many experimental factors and therefore, differ considerably from pure FRET efficiencies ED (s) (R) and EA (s) (R). Pure FRET efficiencies ED,A (s) (R) are calculated with the help of rate equations for D-A pairs, whose molecules have triplet states. It is shown how the calculated efficiencies ED,A (s) (R) can be compared to FRET efficiencies measured with the help of the intensities ID,A(R) corrected by cross talk and background light.

Quantum-Chemical Studies on Excitation Energy Transfer Processes in BODIPY-Based Donor-Acceptor Systems

BODIPY-based excitation energy transfer (EET) cassettes are experimentally extensively studied and serve as excellent model systems for the investigation of photophysical processes, since they occur in any photosynthetic system and in organic photovoltaics. In the present work, the EET rates in five BODIPY-based EET cassettes in which anthracene serves as the donor have been determined, employing the monomer transition density approach (MTD) and the ideal dipole approximation (IDA). To this end, a new computer program has been devised that calculates the direct and exchange contributions to the excitonic coupling (EC) matrix element from transition density matrices generated by a combined density functional and multireference configuration interaction (DFT/MRCI) calculation for the monomers. EET rates have been calculated according to Fermi's Golden Rule from the EC and the spectral overlap, which was obtained from the calculated vibrationally resolved emission and absorption spectra of donor and acceptor, respectively. We find that the direct contribution to the EC matrix element is dominant in the studied EET cassettes. Furthermore, we show that the contribution of the molecular linker to the EET rate cannot be neglected. In our best fragment model, the molecular linker is attached to the donor moiety. For cassettes in which the transition dipole moments of donor and acceptor are oriented in parallel manner, our results confirm the experimental findings reported by Kim et al. [J. Phys. Chem. A 2006, 110, 20-27]. In cassettes with a perpendicular orientation of the donor and acceptor transition dipole moments, dynamic effects turn out to be important.

Hole-transfer induced energy transfer in perylene diimide dyads with a donor–spacer–acceptor motif

Pump–probe spectroscopy, time resolved fluorescence, chemical variation and quantum chemical calculations reveal an efficient energy transfer mechanism enabled by a bright charge transfer state located on the spacer.

Emergence of Coherence through Variation of Intermolecular Distances in a Series of Molecular Dimers

Quantum coherences between electronically excited molecules are a signature of entanglement and play an important role for energy transport in molecular assemblies. Here we monitor and analyze for a homologous series of molecular dimers embedded in a solid host the emergence of coherence with decreasing intermolecular distance by single-molecule spectroscopy and quantum chemistry. Coherent signatures appear as an enhancement of the purely electronic transitions in the dimers which is reflected by changes of fluorescence spectra and lifetimes. Effects that destroy the coherence are the coupling to the surroundings and to vibrational excitations. Complementary information is provided by excitation spectra from which the electronic coupling strengths were extracted and found to be in good agreement with calculated values. By revealing various signatures of intermolecular coherence, our results pave the way for the rational design of molecular systems with entangled states.

Core hole-electron correlation in coherently coupled molecules

We study the core hole-electron correlation in coherently coupled molecules by energy dispersive near edge x-ray absorption fine-structure spectroscopy. In a transient phase, which exists during the transition between two bulk arrangements, 1,4,5,8-naphthalene-tetracarboxylicacid-dianhydride multilayer films exhibit peculiar changes of the line shape and energy position of the x-ray absorption signal at the C K-edge with respect to the bulk and gas phase spectra. By a comparison to a theoretical model based on a coupling of transition dipoles, which is established for optical absorption, we demonstrate that the observed spectroscopic differences can be explained by an intermolecular delocalized core hole-electron pair. By applying this model we can furthermore quantify the coherence length of the delocalized core exciton.

Energy transfer at the single-molecule level: synthesis of a donor-acceptor dyad from perylene and terrylene diimides

In 2004, we reported single-pair fluorescence resonance energy transfer (spFRET), based on a perylene diimide (PDI) and terrylene diimide (TDI) dyad (1) that was bridged by a rigid substituted para-terphenyl spacer. Since then, several further single-molecule-level investigations on this specific compound have been performed. Herein, we focus on the synthesis of this dyad and the different approaches that can be employed. An optimized reaction pathway was chosen, considering the solubilities, reactivities, and accessibilities of the building blocks for each individual reaction whilst still using established synthetic techniques, including imidization, Suzuki coupling, and cyclization reactions. The key differentiating consideration in this approach to the synthesis of dyad 1 is the introduction of functional groups in a nonsymmetrical manner onto either the perylene diimide or the terrylene diimide by using imidization reactions. Combined with well-defined purification conditions, this modified approach allows dyad 1 to be obtained in reasonable quantities in good yield.

Blinking fluorescence of single donor-acceptor pairs: important role of "dark'' states in resonance energy transfer via singlet levels

The influence of triplet levels on Förster resonance energy transfer via singlet levels in donor-acceptor (D-A) pairs is studied. Four types of D-A pair are considered: (i) two-level donor and two-level acceptor, (ii) three-level donor and two-level acceptor, (iii) two-level donor and three-level acceptor, and (iv) three-level donor and three-level acceptor. If singlet-triplet transitions in a three-level acceptor molecule are ineffective, the energy transfer efficiency E=I_{A}/(I_{A}+I_{D}), where I_{D} and I_{A} are the average intensities of donor and acceptor fluorescence, can be described by the simple theoretical equation E(F)=FT_{D}/(1+FT_{D}). Here F is the rate of energy transfer, and T_{D} is the donor fluorescence lifetime. In accordance with the last equation, 100% of the donor electronic energy can be transferred to an acceptor molecule at FT_{D}≫1. However, if singlet-triplet transitions in a three-level acceptor molecule are effective, the energy transfer efficiency is described by another theoretical equation, E(F)=F[over ¯](F)T_{D}/[1+F[over ¯](F)T_{D}]. Here F[over ¯](F) is a function of F depending on singlet-triplet transitions in both donor and acceptor molecules. Expressions for the functions F[over ¯](F) are derived. In this case the energy transfer efficiency will be far from 100% even at FT_{D}≫1. The character of the intensity fluctuations of donor and acceptor fluorescence indicates which of the two equations for E(F) should be used to find the value of the rate F. Therefore, random time instants of photon emission in both donor and acceptor fluorescence are calculated by the Monte Carlo method for all four types of D-A pair. Theoretical expressions for start-stop correlators (waiting time distributions) in donor and acceptor fluorescence are derived. The probabilities w_{N}^{D}(t) and w_{N}^{A}(t) of finding N photons of donor and acceptor fluorescence in the time interval t are calculated for various values of the energy transfer rate F and for all four types of D-A pair. Comparison of the calculated D and A fluorescence trajectories with those measured by Weiss and co-workers proves the important role of triplet levels in energy transfer via singlet levels.

The S1/S2 exciton interaction in 2-pyridone·6-methyl-2-pyridone: Davydov splitting, vibronic coupling, and vibronic quenching

The excitonic splitting between the S(1) and S(2) electronic states of the doubly hydrogen-bonded dimer 2-pyridone[middle dot]6-methyl-2-pyridone (2PY·6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C(2h) symmetric (2-pyridone)(2) homodimer, in which the S(1) ← S(0) transition is symmetry-forbidden but the S(2) ← S(0) transition is allowed, the symmetry-breaking by the additional methyl group in 2PY·6M2PY leads to the appearance of both the S(1) and S(2) origins, which are separated by Δ(exp) = 154 cm(-1). When combined with the separation of the S(1) ← S(0) excitations of 6M2PY and 2PY, which is δ = 102 cm(-1), one obtains an S(1)/S(2) exciton coupling matrix element of V(AB, el) = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S(1)/S(2) ← S(0) spectrum of 2PY·6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a(') vibration and to the intermolecular σ(') stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a(')) = 0.15 and C(σ(')) = 0.05, which places this dimer in the weak-coupling limit. However, the S(1)/S(2) state exciton splittings Δ(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra- and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY·6M2PY the electronic splitting is quenched by a factor of ~30 (i.e., the vibronic quenching factor is Γ(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine·6-methyl-2-pyridone (2HP·6M2PY) are also observed and assigned.

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

The key role of the environment in electronic energy transfer has been underscored in recent experimental and theoretical studies. In this perspective, we provide an overview of novel quantum-mechanical methodologies aimed at describing environment effects in energy transfers. The techniques described include continuum dielectric and atomistic descriptions of the surroundings. We discuss the advantages and limitations of each technique, as well as the main insights that have emerged from their application to solvated dyads and photosynthetic pigment-protein complexes. We finally highlight the aspects that still need to be solved in order to provide a full theoretical route to the study of energy transfer phenomena in complex environments.