Lorenzo
Cupellini
*,
Samuele
Giannini‡
and
Benedetta
Mennucci
*
Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy. E-mail: Benedetta.Mennucci@unipi.it; lorenzo.cupellini@for.unipi.it
First published on 24th November 2017
Photoinduced electron transfer (ET), hole transfer (HT), charge recombination (CR) and energy transfer (EET) are fundamental mechanisms, which occur in both natural and artificial light harvesting systems. Here, we present a computational strategy which determines ET, HT, CR and EET rates in a consistent way and merges them in a kinetic model to reproduce the net excited state dynamics. The effects of the solvent are included in all steps of the calculations making the present strategy a useful tool for a rational design of charge and energy transfer processes in complex systems. An application to covalently linked zinc and free-base porphyrin–naphthalenediimide dyads is presented. For each of the two systems, ultrafast optical spectroscopy experiments have shown a specific photophysics with different processes taking place simultaneously. The model reveals that such a diversity is mainly due to the different relative stability of the charge-separated state, while the electronic couplings for charge and energy transfer processes are quite similar in the two dyads.
Among the donor–acceptor dyads, great interest has been focused on possible materials for organic photovoltaics (OPVs),16 such as dyads based on 1,4,5,8-naphthalenediimides (NDIs). NDI compounds, also known as naphthalene carbodiimides, have attracted much attention due to their chemical characteristics, such as self-organization, and their tunable optical properties.17 NDI-based systems, bonded with different molecules, have been experimentally studied and their excited-state dynamics have been investigated.10,18–20 Yuschenko and coworkers characterized two NDI-based dyads that exhibit a rich photoinduced dynamics after excitation of either of the chromophoric units.21 These dyads are composed of zinc porphyrins (ZnP) or free-base porphyrins (FbP) as electron donors and a hydroquinoxaline-substituted bromo-NDI (HNDI) as an electron acceptor. The two moieties are connected through an essentially rigid amide bond, which allows for short interchromophore separations without altering the electronic structure of the two units (see Fig. 1).
Despite the similar structure, ZnP–HNDI and FbP–HNDI present very different excited-state photophysics, as revealed using time-resolved fluorescence and femtosecond TA experiments.21 Moreover, several processes, such as excitation energy transfer (EET), electron and hole transfers (ET and HT), as well as charge-recombination (CR), occur on similar time scales, preventing a simple interpretation of the excited-state pathways.
In this work, we established a computational protocol to describe EET, HT, ET, and CR in the two porphyrin–HNDI dyads. This protocol uses Density Functional Theory (DFT) based approaches for the calculation of the needed parameters, namely, the excitation energies, the electronic couplings and the Franck–Condon weighted density of states (FCWD), and uses a kinetic model to quantify the time scales involved in the network of charge and energy transfer processes. Diabatization schemes provide a route for computing electronic couplings taking into account the through-bond mechanism, which is a non-negligible effect in these linked systems as confirmed in this work. Marcus–Levich–Jortner's formalism is compared with the more common Marcus equation, starting from ab initio calculations of reorganization energy and driving force. Our results not only reproduce the experimental findings, but they also allow us to gain a deeper insight into the excited-state dynamics of the two investigated dyads. In particular, we highlight that both ET and HT mechanisms contribute to the charge separation in FbP–HNDI, contrary to what the authors assumed in ref. 21, i.e. charge separation occurs via ET. For the same dyad we were also able to calculate the CR rate, which, being much faster than charge separation, could not be estimated from the experiments. We found that the relative stability of the charge-separated state is of great importance, as it influences the different transfer mechanisms and the actual photophysics of the two NDI-based systems.
This paper is organized as follows. In Section 2.1, we present the computational methodologies used to obtain couplings between states participating in the transfer mechanisms. In Section 2.2 we elucidate the Frank–Condon weighted density factor accounting for the solvent effect on reorganization energy and driving force. Computational details are then reported in Section 2.3. In Section 3.1 we analyze excitation energies and coupling values and assess the influence of the covalent bond in ZnP–HNDI on the EET coupling. The rates obtained for different processes such as EET and charge transfer are reported in Sections 3.2 and 3.3 respectively. These results are used in Section 3.4 to build a kinetic model to study the effective dynamics of such systems and compare with experiments.
(1) |
With regard to the EET process, the electronic coupling between two moieties can be calculated, within a first-order perturbative approximation, directly from the transition densities of the noninteracting moieties. This perturbative approach was developed in the TD-DFT framework, with the assumption that the molecular orbitals can be well separated into two moieties.24 In the following we shall refer to this approach as the Direct Coupling (DC) method. Within this framework the coupling becomes:25
(2) |
The DC approach needs to define separated fragments in order to compute the corresponding transition densities. As in our case the donor and acceptor are covalently bonded, we divided the dyads by cutting the single bond between the phenyl bridge and the amide nitrogen, and saturated the free valences with hydrogen atoms.26,27 This approximation is valid if the localized transitions are well separated and do not interact through the covalent bond, which cannot be excluded for this system. Finally, the DC approach can only be used for EET couplings, since the charge-transfer couplings involve states that are not localized on either the donor or the acceptor.
In order to deal with charge-transfer couplings and through-bond effects, we resort to the fragment difference diabatizaion schemes.28,29 Within this framework, an excited-state calculation on the entire system is performed to yield the adiabatic states, which are a combination of diabatic (i.e. charge- or excitation-localized) electronic states. The transformation to diabatic states is achieved through an additional operator that is diagonal for localized states. The Fragment Charge Difference (FCD)30 was introduced for the calculation of charge transfer couplings: it defines the diabatic states as those where the charge is maximally localized on either the donor or the acceptor. The operator to be used is:
(3) |
In analogy to the FCD method, the Fragment Excitation Difference (FED) method32 defines an operator Δx which is diagonal when the excitation is localized on either fragment. Similar to Δq, Δx is defined as:
(4) |
When degenerate or almost degenerate states are present, we report the effective coupling as the square root of the sum of all couplings involving initial and final states, that is, assuming that the degenerate states can be combined into one state maximally coupled to the final state and one uncoupled state. We assume that the uncoupled state does not participate in the observed excited-state dynamics. For example, for two degenerate initial states i1, i2 and two degenerate final states f1, f2, the effective coupling will be V2eff = V2i1,f1 + V2i1,f2 + V2i2,f1 + V2i2,f2. With this definition, the coupling Veff does not change for rotations of pairs of degenerate states.
(5) |
(6) |
The total reorganization energy, λ = λi + λs, is usually estimated as the average of backward and forward reorganization energies (λBw and λFw,34,35 respectively). Due to the large torsion angle between the two moieties which prevents conjugation, we can safely separate the two fragments and define λBw and λFw as:34
λBw = E(D*@D+) − E(D*@D*) + E(A@A−) − E(A@A) λFw = E(D+@D*) − E(D+@D+) + E(A−@A) − E(A−@A−) | (7) |
−ΔGET,HT = E0–0 − (−ΔGCR) −ΔGCR = e(Eox − Ered) + ΔGS + ΔGB | (8) |
Oxidation and reduction potentials were computed as Eox(D) = E(D+@D+) − E(D@D) and Ered(A) = E(A@A) − E(A−@A−), where D is either ZnP or FbP and A is HNDI. The solvent was considered in equilibrium, as discussed in Section 2.1. Finally, the binding energy ΔGB between the two oppositely charged fragments was calculated as the electrostatic interaction between atomic charges of the D and A, scaled by the static dielectric constant ε of THF or toluene to account for the solvent screening, namely:
(9) |
Moving to the EET processes, the FCWD factor is readily accessible from the experimental donor fluorescence FD() and acceptor absorption AA():
(10) |
Once the electronic coupling and the FCWD factors for all the processes were obtained, the rates are calculated as in eqn (1) and inserted in a set of kinetic equations for the state populations Pi:
(11) |
The two solvents, THF and toluene, included in the reference work,21 have been considered in all calculations by employing the PCM in its IEF-PCM formulation.43 The effect of the solvent has been included both in ground and excited state calculations using either an equilibrium or a non-equilibrium model as discussed in Section 2.1.
We note that the CAM-B3LYP functional correctly describes the energy of the Q bands (see Table S1 in the ESI†), while the error on the Ln transition is quite large (∼0.4 eV). In order to benchmark these values, we tested different functionals and the results are reported in Table S3 of the ESI.† All tested functionals overestimate the Ln energy, with only B3LYP showing a rather good description (error of ∼0.1 eV). However, as the transition dipole magnitudes are uniform across all functionals, we can argue that the transition density of Ln is very robust despite the variability in excitation energies. As our calculations are aimed to estimating accurately electronic couplings, a good evaluation of the transition densities is the most important issue28 and we can use the CAM-B3LYP functional as a suitable method to describe these systems including their charge-transfer states.45
We selected the ZnP–HNDI system to assess the through-bond contribution to EET by comparing the DC and FED approaches.32 The analysis of NTOs (see Fig. S2 in the ESI†) shows that the Q states are slightly delocalized over the phenyl bridge; the cut shown in Fig. 1 was chosen so as to maximally separate the transition densities on the two chromophoric units. We considered the real covalently bonded system, as well as a model system in which the two fragments are artificially separated at a distance of 1.2 Å, cutting the covalent bond and saturating the free valences with hydrogen atoms. The resulting couplings are shown in Fig. 2, in vacuo and in solvent. Both methods predict a screening effect due to the solvent that surrounds the dyad, which decreases the EET couplings.25
When the moieties are bonded, the couplings computed with the FED are substantially larger than those computed with the DC approach. Contrary to the FED calculation, DC treats the molecular fragments separately, therefore it cannot include the effect of a covalent bond between the two moieties. The observed difference between FED and DC values could indeed arise from the presence of a covalent bond, but also from the mutual polarization of the two fragments, which is included only in the FED approach. When the two molecules are set apart at a distance of 1.2 Å, the differences between the two approaches almost completely vanish. The remaining small deviation can be ascribed to mutual polarization effects between the two fragments, whereas the large difference observed in the covalently bonded molecules can be regarded as a through-bond effect.
The through-bond contribution to the total coupling can be quantified as:
(12) |
To gain a better understanding of the impact of through-bond effects on the coupling values in ZnP–HNDI, we analyze this contribution as a function of the dihedral angle between the chlorin ring of the ZnP and the phenyl substituent. This angle is assumed to have a quite low activation energy barrier.21 In order to avoid any spurious effect from the PCM cavity, we used in vacuo calculations for this analysis. The calculations were also performed with the smaller 6-31G+(d) basis set to reduce the computational efforts. As we want to evaluate a qualitative trend, rather than retrieving the quantitative results, those two precautions do not influence our analysis. The results are reported in Fig. 3. When the dihedral angle is more planar, the π orbitals of the phenyl become conjugated with those of the porphyrin. This increases the percentage of through-bond transfer to more than 40%. Instead, when the angle is 90 degrees, the through-bond contribution hits a minimum of 30%. At all values of the dihedral angle the Coulomb coupling remains the largest contribution, but the through-bond contribution is never negligible.
Fig. 3 Coupling (solid lines, cm−1) and through-bond contribution (dashed line, %) with respect to the dihedral angle between the chlorin ring of ZnP–HNDI and the phenyl substituent. The FED coupling is shown (blue) along with the DC (red) and the sole Coulomb contribution (yellow). All quantities refer to in vacuo calculations (see the text). Solid colored lines represent fits of the calculated points through eqn (13). |
Both Coulomb and FED couplings can be fitted very well with an expression depending on the squared cosine of the dihedral angle (see Fig. 3):
V(ϕ) = V0 + βcos2(ϕ) | (13) |
For the Coulomb coupling, this can be understood in terms of the orientational factor between transition dipoles. As the through-bond coupling depends on the extent of π-conjugation, we can expect a similar dependence on the dihedral angle. Indeed, the cos2(ϕ) dependence was found to hold experimentally for electron transfer couplings,46 singlet exciton coupling,47 and triplet transfer coupling.48 The latter case is particularly interesting since in triplet transfer the Coulomb contribution is zero.29 The triplet transfer coupling is thus a measure of the through-bond contribution, which is found to have the same torsional dependence.
Dyad | States | THF | Toluene | ||
---|---|---|---|---|---|
V | τ | V | τ | ||
ZnP–HNDI | Q → Ln | 63 | 1.4 (0.5) | 54 | 1.9 (0.5) |
FbP–HNDI | Ln → Qx | 46 | 1.6 (2.7) | 39 | 2.2 (2.7) |
The EET processes reported in Table 1 are in equilibrium with the relative backward reactions. Considering the experimental spectral overlaps for Ln → Q and Q → Ln, of ZnP–HNDI and FbP–HNDI, respectively, the equilibrium constants are 7.4 and 86.7 in favour of Ln and Qx, respectively, for ZnP–HNDI and FbP–HNDI.
ET | HT | CR | ||||
---|---|---|---|---|---|---|
THF | Toluene | THF | Toluene | THF | Toluene | |
ZnP–HNDI | ||||||
λ i | 0.28 | 0.28 | 0.26 | 0.26 | 0.36 | 0.36 |
λ s | 1.01 | 0.06 | 0.97 | 0.01 | 0.97 | 0.01 |
ΔG | −0.35 | 0.13 | −0.30 | 0.18 | −1.75 | −2.23 |
FbP–HNDI | ||||||
λ i | 0.29 | 0.29 | 0.21 | 0.21 | 0.31 | 0.31 |
λ s | 1.02 | 0.06 | 1.02 | 0.06 | 1.02 | 0.06 |
ΔG | −0.08 | 0.33 | −0.23 | 0.18 | −1.82 | −2.23 |
As expected from the larger polarity, in THF, the solvent represents the main contribution to the reorganization energy, even though the internal reorganization is not negligible for these systems. By contrast, toluene gives only a small contribution. In both systems the charge-separated state in toluene is energetically much higher than the one in THF. The driving force for charge separation in toluene is positive in all cases, making the charge-separated state thermodynamically unreachable, in agreement with experimental findings.21 In addition, we observe that in THF ΔGCR < −λ, meaning that charge recombination occurs in the Marcus inverted region, whereas ET and HT occur in the direct region.
The Marcus rates for ET, HT and CR processes were calculated with eqn (5), while the Marcus–Levich–Jortner results were obtained with eqn (6). In the latter case the effective frequency of the mode coupled with the transfer was estimated from a fitting of the absorption spectra of the HNDI molecule as described in the ESI,† obtaining a frequency of 1210 cm−1.
The electronic couplings are reported in Table 3. In order to retrieve completely localized states from this approach and obtain meaningful couplings, we included more than two states in the diabatization procedure, as described in Section S5 of the ESI.†
Dyad | Process | V | τ (Marcus) | τ (MLJ) | |
---|---|---|---|---|---|
ZnP–HNDI | ET | Q → CSS | 119 | 247 | 98 |
HT | Ln → CSS | 225 | 79 | 32 | |
CR | CSS → GS | 73 | 3.2 | 2.7 | |
FbP–HNDI | ET | Qx → CSS | 82 | 43 × 103 | 15 × 103 |
HT | Ln → CSS | 196 | 279 | 131 | |
CR | CSS → GS | 68 | 5.6 | 4.1 |
In ZnP–HNDI the charge-separated state is experimentally observed,21 rising with a time constant of 18 ps and decaying with a time constant of 7.7 ps. Our computed charge separation rates are close to the experiment but slightly slower, with the HT rate being faster than the ET rate. As the Q → Ln EET is faster than both of these rates, the charge separation starts from the Ln state of HNDI and proceeds through hole transfer. Finally, the calculated charge recombination rate is slightly overestimated with respect to the experiment. Nonetheless, the overall time scales of EET, charge separation and recombination are well reproduced by our calculation.
The charge-separated state of FbP–HNDI was only indirectly observed in transient absorption experiments.21 However, the lifetime of the FbP's Qx state shortens in THF compared to toluene. This was explained by the authors with the appearance of a new decay pathway involving the charge-separated state, which recombines much faster to the ground state.21 The charge-separation time constant was estimated to be around 10 ns by comparing the decay rates in toluene and THF. Our calculated ET time constant of ∼15 ns is in agreement with this experimental finding; however, we observe a much faster rate for the HT process. The Ln and Qx states are in equilibrium (see above), as also seen from transient absorption,21 therefore the HT pathway could contribute to the charge separation. In order to obtain a clear picture of the excited-state dynamics in this dyad, a kinetic model is discussed in the following section.
On a final note, we point out that Marcus eqn (5) always gives smaller rates than the MLJ equation, eqn (6). However, the MLJ rates are no more than three times larger than the Marcus rates. As regards the internal, discrete modes of the molecule, Marcus and MLJ equations describe two opposite temperature limits, that is, the high-temperature limit for the Marcus equation, and the zero-temperature limit for the MLJ equation. The similar order of magnitude given by the two theories indicates that the calculated rates are robust with respect to the choice of the approximation, and also to the frequency of the effective internal mode.
To assign these rates, we follow the population evolution of the two dyads reported in Fig. 4. In particular, looking at the excited-state dynamics of ZnP–HNDI in Fig. 4a, we notice that the Q states, initially populated upon excitation, transfer electronic population to the Ln state of HNDI with a time constant of 1.2 ps. The electronic population of the CSS starts rising after a few picoseconds and decays to the ground state following charge recombination. In the meantime, the electronic population of Ln decays with a time constant of 35 ps. Globally, the excited state dynamics for ZnP–HNDI is considerably fast and occurs on a picosecond time scale. The calculated effective rates are fairly close to the experiment, the largest error being on the charge recombination, whose rate is slightly overestimated by our calculations.
The excited-state dynamics of FbP–HNDI takes place on two different time scales (see Fig. 4b). The faster time scale sees the growth of the Qx state following EET from the initially excited Ln state (1.6 ps). On the slower time scale, the CSS begins to rise, but the ultrafast charge recombination consumes its population. Indeed, throughout the kinetic evolution, the CSS population remains always less than 1% of the population of the Qx/Ln states, explaining why the CSS state could not be detected from transient absorption measurements. At longer times Qx decays with a time constant of 4.4 ns, matching the rising time of the GS and the triplet state. The Qx decay rate contains both ET and HT processes, as well as the ISC. Interestingly, in FbP–HNDI, even though the first EET equilibrium is mainly towards the Qx state, the charge separation also proceeds mainly through HT from the Ln state, contrary to what was assumed in the experimental work.21
On the whole, the photophysics of FbP–HNDI in THF is captured by the kinetic model with the charge-separation decay pathway through both ET and HT. The decay rate of the Qx state only depends on the rates of the elementary ET and HT processes, as well as the Qx/Ln equilibrium constant. As long as the CR rate remains much larger than charge separation rates, the population evolution of the Qx and Ln states does not depend on the CR rate. As a consequence, this rate cannot be estimated from experiment only, unless the system is prepared in the charge-separated state. Based on our calculations, the CR rate is on the same order as in ZnP–HNDI, but slightly lower due to the larger driving force, which puts this process further into the inverted Marcus region. Moreover, as both ET and HT contribute to the charge separation, the estimated charge-separation rate cannot be directly assigned to either of these processes.
The diverse photophysics presented by the two dyads, quantitatively reproduced by the model, can be ascribed mainly to energetic factors. The electronic couplings for charge transfer processes are quite similar in the two dyads (Table 3), as can be expected from the structural similarity. However, the higher energy of the charge-separated state determines a completely different dynamics in FbP–HNDI. This suggests that the photophysics of these dyads can be drastically changed just by tuning the energy of the charge-separated state.
We have investigated the photoinduced energy and charge transfers characterizing the excited state dynamics of the ZnP–HNDI and FbP–HNDI dyads for which previous studies based on ultrafast spectroscopic techniques and global analysis of transient absorption21 showed the simultaneous presence of different charge and electron transfer pathways.
In agreement with the conclusions based on the experimental evidence, we have found that the excited state dynamics of ZnP–HNDI occurs on a picosecond time scale and after an EET between ZnP and HNDI, the HT from the latter fragment is the main cause for the considerable fast charge separation in such a system. In the case of FbP–HNDI, we have confirmed that the charge separation is very slow, and the following (fast) charge recombination immediately consumes the population of the charge-separated state, which cannot be detected experimentally. Notably, despite the EET equilibrium is mainly shifted towards the FbP moiety, we have found that the HT pathway from the HNDI still operates as a fast quenching mechanism for the FbP excited state. Such mechanisms contribute, along with the direct ET, to the charge separation rate. In both these donor–acceptor systems, the energetic position of the CSS state is a key component in determining their photophysics. We can therefore expect that, upon varying the oxidation and reduction potentials of the two moieties and/or the polarity of the solvent, a drastic change in their excited state dynamics will be obtained.
The results reported in this work demonstrate that a complete description of photoinduced transfer processes is possible by using an ab initio description of the single processes which are finally integrated in a kinetic model giving the full evolution of the electronic populations. We stress that our rate calculations only needed experimental absorption and fluorescence spectra, while all the other parameters are computed from quantum mechanical calculations. As a matter of fact, the same strategy can be extended to quantum chemical methods which are more accurate than TDDFT for the simulation of the optical spectra, thereby enhancing its independence of the experiments and the predictive potential.
Footnotes |
† Electronic supplementary information (ESI) available: Structure and fragment definition; natural transition orbitals of the molecules considered in this work; excitation energies and functional dependence; basis set dependence; benchmark of FCD couplings; calculation of charge transfer driving force; details of the fitting of spectra; details of the kinetic modelling. See DOI: 10.1039/c7cp07002k |
‡ Current affiliation: Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK |
This journal is © the Owner Societies 2018 |