Miftahussurur Hamidi Putra1, Sebastian Seidenath2, Stephan Kupfer2, Stefanie Gräfe2, Axel Groß1,3. 1. Institute of Theoretical Chemistry, Ulm University, 89069, Ulm, Germany. 2. Institute for Physical Chemistry (IPC) and Abbe Center of Photonics, Friedrich Schiller University Jena, Helmholtzweg 4, 07743, Jena. 3. Helmholtz Institute Ulm (HIU), Electrochemical Energy Storage, 89069, Ulm, Germany.
Abstract
Conductive polymers represent a promising alternative to semiconducting oxide electrodes typically used in dye-sensitized cathodes as they more easily allow a tuning of the physicochemical properties. This can then also be very beneficial for using them in light-driven catalysis. In this computational study, we address the coupling of Ru-based photosensitizers to a polymer matrix by combining two different first-principles electronic structure approaches. We use a periodic density functional theory code to properly account for the delocalized nature of the electronic states in the polymer. These ground state investigations are complemented by time-dependent density functional theory simulations to assess the Franck-Condon photophysics of the present photoactive hybrid material based on a molecular model system. Our results are consistent with recent experimental observations and allow to elucidate the light-driven redox chemical processes - eventually leading to charge separation - in the present functional hybrid systems with potential application as photocathode materials.
Conductive polymers represent a promising alternative to semiconducting oxide electrodes typically used in dye-sensitized cathodes as they more easily allow a tuning of the physicochemical properties. This can then also be very beneficial for using them in light-driven catalysis. In this computational study, we address the coupling of Ru-based photosensitizers to a polymer matrix by combining two different first-principles electronic structure approaches. We use a periodic density functional theory code to properly account for the delocalized nature of the electronic states in the polymer. These ground state investigations are complemented by time-dependent density functional theory simulations to assess the Franck-Condon photophysics of the present photoactive hybrid material based on a molecular model system. Our results are consistent with recent experimental observations and allow to elucidate the light-driven redox chemical processes - eventually leading to charge separation - in the present functional hybrid systems with potential application as photocathode materials.
The effective and efficient conversion of solar energy into electrical energy is challenging but might enable a significant reduction of the carbon emission in the world which will even have a beneficial economic impact.
The market of solar cells is still dominated by silicon‐based photovoltaic devices.[
,
,
] However, dye‐sensitized solar cells (DSSC)
have the potential to become an alternative to silicon‐based solar cells, as they are associated with a very low weight, low costs, ease of production, and can be used in a more flexible manner than the conventional silicon‐based technology.[
,
,
,
,
]In general, there are two kinds of DSSCs, i. e. n‐type and p‐type. In n‐type DSSCs, the electrons are injected from the dye to the photoanode (mainly TiO2). In contrast, hole injection occurs in p‐type DSSCs from the dyes to the photocathode made of, e. g., NiO.[
,
,
] Typically, p‐type DSSCs feature an efficiency of about 2.5 %, roughly one order of magnitude lower than the efficiency of n‐type DSSCs (13.1 %).
Interestingly enough, tandem systems combining n‐type with p‐type DSSCs were predicted to achieve efficiencies of up to around 40 %.
Still, there is a strong need to increase the efficiency of p‐type DSSCs to make them more competitive, e. g. by developing new photocathode materials with adjustable electronic properties. In this sense, replacing inorganic by organic cathode materials is highly promising.The fundamental concept as well as the underlying light‐driven processes in p‐type DSSCs are illustrated in Figure 1. Firstly, an electron of the dye is excited by a photon into an excited state. Subsequently, the hole is injected from the dye to the highest occupied molecular orbital (HOMO) of the photocathode and the hot electron is transferred from the dye to the redox or catalytic agent.[
,
] In the operation of the DSSC, undesired electron transfer pathways should be suppressed. Such undesired processes are, for example, the electron transfer from the dye to the photocathode or the direct recombination of the hole and the electron from the excited state of dye to the HOMO of the photocathode.
In addition, the electron may also be transferred from the excited state of the photocathode to the dye, as the energy gap of this semiconductor is typically comparable to the energy gap of the dye.
Figure 1
Schematic presentation of the processes occurring in a p‐DSSC of a photoactive complex coupled to a PPV polymer including the full DSSC circuit. The left panel illustrates the beneficial processes whereas the right panel depicts the parasitic phenomena.
Schematic presentation of the processes occurring in a p‐DSSC of a photoactive complex coupled to a PPV polymer including the full DSSC circuit. The left panel illustrates the beneficial processes whereas the right panel depicts the parasitic phenomena.The concept of p‐DSSCs can also be extended to photochemical molecular devices (PMD) by combining a photosensitizer and a catalyst connected by an electron relay.[
,
,
] In such a PMD, the electron is directly transferred to the catalytic agent, e. g. a hydrogen evolution catalyst, instead to a redox agent, such as I
/I−.[
,
,
,
] Furthermore, apart from just addressing the properties of the DSSCs, hole conducting organic polymers have been put forward as viable alternatives to NiO or other metal oxide materials, as the tuning of the physicochemical properties of polymers is easier
than those of NiO.Thus the integration of molecular light‐driven catalysts into soft matter matrices presents an attractive route for tailored photochemical molecular devices as the embedding can give control over charge transfer processes, photochemical reactivity and degradation resistance.
However, in order to do this in a rational fashion, it is necessary to understand the properties of polymers and its coupling to the photosensitizer. In order to contribute to an improved understanding of this coupling, we have performed quantum chemical calculations on a soft matter embedded Ru(II)‐based photosensitizer. The choice of the considered systems has been motivated by a recent experimental study.
In order to properly capture both the delocalized nature of electronic states of the polymers as well as the local nature of the electronic excitations of the photoactive complex, we combined two complementary approaches: On the one hand, we have carried out periodic density functional theory (DFT) calculations modelling in principle an infinitely long polymer coupled to a periodic array of photoactive transition metal complexes. This allows to capture the delocalized nature of the electronic system of the polymers. On the other hand, in order to assess light‐driven processes such as electronic excitation spectra, we employed time‐dependent density functional theory (TDDFT) onto local models of the complexes.The organization of this paper is as follows: We will first present details of the quantum chemistry methods used in this work. Then, the results of our quantum chemistry calculations will be presented and discussed together with a comparison with the available experimental data. Finally, conclusions are provided.
Computational Details
In the first computational approach used in the study, we modeled the polymer as an infinite chain with a periodic set‐up.
As the polymer, we considered the conductive polymer poly(p‐phenylenevinylene) (PPV) as functionalized PPV derivates are commonly used in polymer solar cells
and have already been functionalized with a Ru(II)‐polypyridyl‐based photosensitizer.
In this experiment,
two photosensitizers with either an electron‐withdrawing ethyl ester group (‐COOEt) or an electron‐donating tert‐butyl group as substituent on their bipyridine ligands were considered. In order to reduce the numerical effort, in the calculations we replaced these substituents by smaller electron‐withdrawing and ‐donating groups, ‐COOH (Ru1) and CH3 (Ru2), respectively, which, however, mimic the electronic effects of the substituents used in the experiment. The structure of the polymer is illustrated in Figure 2a, while the PPV‐embedded photosensitizers (Ru1 and Ru2) – chemically linked by alkyl side chains – are depicted in Figure 2b. In the periodic cell, we include four repeat units of the polymer and replace the remaining C7H15 alkyl chains by a methyl group to reduce the computational effort. Such a substitution is known to only slightly modify the electronic structure of the polymer.
Furthermore, we also included two PF
anions in the periodic unit cell, as indicated in Figure 2b, in order to counter‐balance the charge (
) of the ruthenium complex.
Figure 2
(a) Structure of PPV polymer; alkyl chains are replaced by methyl groups to reduce the computational cost. (b) Unit cell used in the periodic DFT calculations to model the polymer and the Ru(II) dye linked to the polymer. Structure of the Ru(II)‐based photosensitizer linked to the PPV polymer (PPV−Ru) together with two PF6
‐ counter ions for charge balancing. Two different substituents have been used, R=COOH and R=CH3, yielding the hybrid materials PPV−Ru1 and PPV−Ru2, respectively.
(a) Structure of PPV polymer; alkyl chains are replaced by methyl groups to reduce the computational cost. (b) Unit cell used in the periodic DFT calculations to model the polymer and the Ru(II) dye linked to the polymer. Structure of the Ru(II)‐based photosensitizer linked to the PPV polymer (PPV−Ru) together with two PF6
‐ counter ions for charge balancing. Two different substituents have been used, R=COOH and R=CH3, yielding the hybrid materials PPV−Ru1 and PPV−Ru2, respectively.The periodic quantum chemical calculations have been performed employing density functional theory (DFT) as implemented in the FHI‐AIMS code[
,
] that utilizes a basis of numeric atom‐centered orbitals. This allows an efficient parallelization of the periodic DFT calculations, even for hybrid functional calculations[
,
,
] which are typically very time‐consuming for infinite systems. To perform the geometry optimization, we used the Perdew‐Burke‐Ernzerhof (PBE) exchange‐correlation functional within the generalized gradient approximation (GGA)
with a light basis set.
The van der Waals dispersion interaction is taken into account using the Tkatchenko‐Scheffler approach.
Relativistic energy corrections are included within the ZORA approximation.
The PBE functional like all GGA functionals underestimates the energy gap. In order to obtain a reliable electronic structure of the polymer linked to the Ru(II) dye, we performed single points calculations of the relaxed structure obtained from PBE optimization by means of hybrid functionals such as B3LYP and PBE0.[
,
,
]To describe the PPV polymer alone, we used a cell size
Å3 with a
k‐point set for the geometry optimization and a
k‐point set for the determination of the electronic gap. For the combined PPV+Ru system, we extended the unit cell to a
supercell with respect to the polymer repeat units with a cell size of
Å3. Because of the significantly larger unit cell, we used a
k‐point set for the PPV+Ru system for both the geometry optimization as well as for the determination of the electronic structure.In order to estimate the change of the properties of polymer and the dye, we determined the electronic properties of short PPV oligomers, i. e. the monomer, dimer and trimer, and the photosensitizer in various implicit solvents, such as gas phase, toluene, acetone, acetonitrile and water represented through their corresponding dielectric constants. These calculations were performed with the Gaussian 16 package.
We used PBE0 as the exchange‐correlation functional
and the D3 dispersion correction scheme with BJ damping to account for dispersion effects.
The implicit solvent was model employing the SMD method.In the analysis of our results, we use the concept of the hole injection efficiency (HJE)
[
,
,
,
] which corresponds to the difference in the energy levels of the HOMO at the dye and the polymer,and which allows to estimate the charge transfer between polymer and dye. On the other hand, a large negative hole injection energy
corresponds to an indication that hole injection might in fact. be more favorable which can be quantified using the electron injection energy
, which is given by the difference in the energies of the lowest unoccupied molecular orbital (LUMO),Equivalently, a large negative value of
indicates a higher likelihood for an electron transfer process from the dye to the polymer. Consequently, efficient p‐DSSCs are characterized by a small negative or even positive value of
and a large negative value of
.In order to directly access the excited states involved in a charge separation process, we also employed Kohn‐Sham time‐dependent density functional theory (TDDFT). Thus, excitation energies, oscillator strengths, transition dipoles and excited state densities are available for analysis. A major advantage compared to simpler models is the intrinsic treatment of dynamic electron correlation in DFT – an advantage that excited state calculations using TDDFT inherit. In particular, TDDFT struggles with describing Rydberg or charge transfer states, which can be overcome by applying range‐separated or hybrid functional with a high degree of exact HF exchange.[
,
,
] In this case, the PBE0 functional (25 % HF exchange) was chosen as previous joint spectroscopic‐theoretical studies highlighted its capability to unambiguously assess the Franck‐Condon photophysics, i. e. excitation energy, transition dipole moments, electronic characters and excited state gradients, of transition metal complexes.[
,
,
]Still, periodic DFT codes typically do not contain TDDFT implementations to derive oscillator strengths. Therefore, we used local models of the extended PPV−Ru systems based on the periodic structures obtained in the calculations described above in order to account for the computational demands of TDDFT calculations. We took the frozen sequence of four benzene units with one Ru dye from the unit cell used in the periodic DFT calculations depicted in Figure 2b as input for our TDDFT calculations. The C−C bonds broken in this process were capped by hydrogen atoms. The positions of these hydrogens were then optimized at the def2‐SVP/PBE[
,
] level of DFT using the ORCA program suite[
,
] while the remaining atoms remained frozen.TDDFT calculations were done using the NWChem code
with the PBE0 functional to maintain comparability with all other results. The def2‐SVP basis set was used together with its corresponding ECPs while the Tamm‐Dancoff approximation (TDA) was applied.[
,
,
] All calculations were done in the gas phase. The lowest 100 singlet states were calculated, covering an energy range of up to 4.1 eV (302 nm).
Results and Discussion
In order to validate the reliability of our results, we start the presentation of our results with a comparison of the calculated electronic properties of the PPV polymer, such as the work function, the position of HOMO and LUMO levels, and the energy gap, with the corresponding experimental results[
,
,
,
,
] (see Table 1). We have evaluated the electronic properties of the polymer in the structure relaxed with the PBE functional using the different functionals, PBE and the hybrid functionals B3LYP and PBE0. The position of HOMO and LUMO levels are given with respect to the vacuum level. First of all the very well‐known fact is reproduced that GGA functionals such as PBE severely underestimate HOMO‐LUMO gaps[
,
] which can be traced back to unphysical self‐Coulomb repulsion contained in these functionals.[
,
] In contrast, the hybrid functionals, i. e. PBEO and B3LYP, including a certain portion of exact exchange do a much better job in this respect.
Table 1
Energy gap, work function (
), and the positions of the HOMO and LUMO levels of PPV calculated using the PBE, B3LYP and PBE0 functionals and compared with experimental values.
No.
XC Corr.
Energy Gap (eV)
Φ
(eV)
HOMO (eV)
LUMO (eV)
1
PBE
1.0144
3.521
−4.052
−3.038
2
B3LYP
2.017
3.210
−4.650
−2.633
3
PBE0
2.271
3.102
−4.824
−2.554
4
Exp.
2.1–2.4[15, 52, 53, 54]
4.000[55]
−4.90 to −5.41[15, 52, 53, 54]
−2.70 to −3.01[15, 52, 53, 54]
Energy gap, work function (
), and the positions of the HOMO and LUMO levels of PPV calculated using the PBE, B3LYP and PBE0 functionals and compared with experimental values.No.XC Corr.Energy Gap (eV)(eV)HOMO (eV)LUMO (eV)1PBE1.01443.521−4.052−3.0382B3LYP2.0173.210−4.650−2.6333PBE02.2713.102−4.824−2.5544Exp.2.1–2.4[
,
,
,
]4.000−4.90 to −5.41[
,
,
,
]−2.70 to −3.01[
,
,
,
]Among the two hybrid functionals, PBE0 outperforms B3LYP with respect to the position of the HOMO and LUMO levels and the energy gap, as far as the comparison with the experiment is concerned. Interestingly enough, PBE reproduces the measured position of the LUMO level better than the two hybrid functionals with PBE0 exhibiting the largest deviation with respect to the experiment. Still we will use PBE0 in the following in the discussion of the electronic properties of the PPV−Ru systems as it yields the best agreement with the measured electronic levels of the PPV polymer.Now, we consider the electronic levels of the joint system with the Ru dye linked to the PPV polymer. As a first step, we show in Figure 3 the band structures of the polymer alone and the two joint systems calculated with the PBE0 functional. In addition, the total density of states (DOS) and the local density of states (LDOS) projected onto the polymer, the Ru atom, and the phenatroline and bipyridine moiety of the photosensitizer are plotted. Figure 3a illustrates the one‐dimensional band structure of the isolated PPV polymer. The polymer has a direct band gap at the Γ point which widens towards the X point. Note that the variation of the electronic bands with the wavevector k, the dispersion, is a measure of the delocalization, localized states just have a flat dispersion. The changes of the bandstructure upon attaching the Ru complex with the ‐COOH substituent to the polymer (PPV−Ru1) are shown in Figure 3b. In addition to the PPV states, now a high number of dispersionless, i. e., k‐independent states appears. The vanishing dispersion indicates that these are localized states of the Ru complex. In fact, the lowest unoccupied state is also dispersionless and located at the phenatroline and bipyridine moiety of the photosensitizer, as the local density of states plotted in Figures 3c and e indicates. This demonstrates that the LUMO of the combined system is located within the photoactive complex. Furthermore, the plotted DOS also confirms that the HOMO remains to be located at the polymer. The PPV bands hardly change upon the attachment of the Ru complex showing that there is obviously a rather small coupling between the electronic states of the PPV and the Ru complex.
Figure 3
Bandstructures and densities of states of the PPV polymer alone and coupled to the Ru photosensitizer with −COOH (PPV−Ru1) and −CH3 (PPV−Ru2) substituents in the gas phase calculated with the PBE0 functional. a, b, d) Band structure of the isolated PPV polymer, the PPV−Ru1 and the PPV−Ru2 systems, respectively; c, e) density of the states (DOS) of the PPV−Ru1 and the PPV−Ru2 systems, respectively. In addition to the total DOS, we also plot the local DOS of the PPV polymer, of the Ru center and of the phenatroline and bipyridine moiety denoted by Ligand‐Dye. In the presentation of the DOS, a smearing of 50 meV has been applied.
Bandstructures and densities of states of the PPV polymer alone and coupled to the Ru photosensitizer with −COOH (PPV−Ru1) and −CH3 (PPV−Ru2) substituents in the gas phase calculated with the PBE0 functional. a, b, d) Band structure of the isolated PPV polymer, the PPV−Ru1 and the PPV−Ru2 systems, respectively; c, e) density of the states (DOS) of the PPV−Ru1 and the PPV−Ru2 systems, respectively. In addition to the total DOS, we also plot the local DOS of the PPV polymer, of the Ru center and of the phenatroline and bipyridine moiety denoted by Ligand‐Dye. In the presentation of the DOS, a smearing of 50 meV has been applied.Upon changing the substituent to −CH3 (PPV−Ru2), the PPV‐derived band structure hardly changes, as Figure 3d shows, but the positions of the dispersionless states are altered. This becomes even more obvious upon comparing the densities of state in Figure 3e and Figure 3c. The peaks of the LDOS at the Ru atom are at different energetic positions. Furthermore, the LUMO of the Ru2 complex is found at higher energies with respect to the valence band edge of the polymer, overall a considerably modified electronic structure results. In summary, the results presented in Figure 3 show that the electronic properties of the polymer are hardly modified by attaching the Ru complex, irrespective of the nature of the two substituents used. The electronic states of the Ru complex hardly couple to the electronic states of the polymer, but they are sensitive to the nature of the substituent being used.Next, we turn to a local picture of the electronic orbitals which has been obtained by extracting the position‐resolved variation of the HOMO and the LUMO, as illustrated in Figure 4. There, also the corresponding energy levels of the PPV−Ru with the ‐COOH and ‐CH3 substituents are given, now with respect to the vacuum level in order to account for work function effects. In fact, our analysis of the spatial distribution of the electronic orbitals of the Ru photosensitizer linked to the polymer only yields orbitals that are almost entirely located either on the polymer or on the photosensitizer so that the distinction between polymer and photosensitizer orbitals is unique which confirms all the findings we derived from the analysis of the band structure and the density of states. Obviously, the polymer and the attached photosensitizer are not strongly electronically coupled. Thus, the alkyl chain connector considered here apparently acts as a geometric rope immobilizing the photosensitizer rather than as an electronic wire coupling the electronic states of polymer and the photosensitizer.
Figure 4
HOMO and the LUMO levels of the PPV polymer and the Ru photosensitizer with the −COOH (a) and −CH3 (b) substituents, respectively, in the gas phase taken with respect to the electronic vacuum level. The numbers in blue represent the HOMO‐LUMO levels of PPV, those in red color denote the HOMO‐LUMO levels of the photosensitizer. These levels have been identified through a charge analysis of the electronic orbitals which is illustrated by the charge isosurfaces in the ball‐and‐stick presentation of the polymer and the Ru photosensitizer.
HOMO and the LUMO levels of the PPV polymer and the Ru photosensitizer with the −COOH (a) and −CH3 (b) substituents, respectively, in the gas phase taken with respect to the electronic vacuum level. The numbers in blue represent the HOMO‐LUMO levels of PPV, those in red color denote the HOMO‐LUMO levels of the photosensitizer. These levels have been identified through a charge analysis of the electronic orbitals which is illustrated by the charge isosurfaces in the ball‐and‐stick presentation of the polymer and the Ru photosensitizer.As Figure 4 confirms, the energies of the HOMO and the LUMO localized on the polymer with the two different substituents are very similar, and the HOMO‐LUMO gap is almost identical. Hence, this again shows that the substituents do apparently not have a significant influence on the electronic properties of the polymer which is not too surprising considering the fact that they are not directly interacting with the polymer. Furthermore, a closer analysis yields that the HOMO on the polymer has a π‐orbital character whereas the LUMO is of
‐orbital character. Interestingly, the HOMO and LUMO levels localized on the Ru photosensitizer differ by a larger amount, as indicated by the different substitution pattern with the electron‐withdrawing (–COOH) vs. the electron‐donating (−CH3) substituents. The corresponding levels of the Ru1 complex with the −COOH substituent are by about 0.3 eV lower than the Ru2 complex with the −CH3 substituent, also in agreement with the results shown in Figure 3. As Figure 4 illustrates, for both substituents the HOMO of the photosensitizer remains to be localized within the
orbitals at the ruthenium metal center. On the other hand, the spatial distribution of the LUMO depends on the choice of the substituent. For the −COOH substituent, the electron within the LUMO are accumulated at the
‐orbital of the bipyridine moiety, whereas in the −COOH case the LUMO electrons are concentrated in the
orbital of the phenantroline moiety.The calculated hole and electron injection energies
and
for the polymer and the photosensitizer in the gas phase have been listed in Table 2 in order to qualitatively assess the hole and electron transfer. The hole injection energy
with respect to the PPV−Ru1 combination is by about 0.4 eV lower than the one of the PPV−Ru2 combination. Hence we predict that the PPV−Ru1 system is associated with a more efficient hole transfer than PPV−Ru2. On the other hand, PPV−Ru2 has a electron injection energy
that is lower by about 0.35 eV than
of PPV−Ru1. This indicates that PPV−Ru2 will have a faster charge recombination than PPV−Ru1.
Table 2
Position of the HOMO and the LUMO and hole and electron injection efficiencies ▵E
HJE and ▵E
EJE calculated for the PPV−Ru system in gas phase and with
water molecules explicitly included. All experimental values are derived from Ref. [15]. The experimental HOMO and LUMO level are estimated by the corresponding redox levels, and the hole and electron injection efficiencies are derived from the measured cyclic voltammograms. All values are given in eV.
Parameter
Gas Phase
NH2O
=10
NH2O
=30
Experiment[15]
Ru1
Ru2
Ru1
Ru2
Ru1
Ru2
Ru1
Ru2
EPPVHOMO
−5.040‐
−5.083
−4.999
−5.012
−4.921
−4.987
−5.41
−5.41
EPPVHOMO
−2.774
−2.812
−2.729
−2.740
−2.681
−2.709
−3.01
−3.01
EDyeLUMO
−7.222
−6.884
−6.932
−6.612
−6.687
−6.372
−5.78
−5.50
EDyeLUMO
−3.729
−3.415
−3.487
−2.968
−3.162
−2.713
−3.35
−3.03
ΔEHJE
−2.182
−1.801
−1.933
−1.601
−1.766
−1.386
−0.37
−0.09
ΔEEJE
0.955
0.604
0.758
0.228
0.482
0.004
0.34
0.02
Position of the HOMO and the LUMO and hole and electron injection efficiencies ▵E
HJE and ▵E
EJE calculated for the PPV−Ru system in gas phase and with
water molecules explicitly included. All experimental values are derived from Ref. [15]. The experimental HOMO and LUMO level are estimated by the corresponding redox levels, and the hole and electron injection efficiencies are derived from the measured cyclic voltammograms. All values are given in eV.ParameterGas Phase=10=30ExperimentRu1Ru2Ru1Ru2Ru1Ru2Ru1Ru2−5.040‐−5.083−4.999−5.012−4.921−4.987−5.41−5.41−2.774−2.812−2.729−2.740−2.681−2.709−3.01−3.01−7.222−6.884−6.932−6.612−6.687−6.372−5.78−5.50−3.729−3.415−3.487−2.968−3.162−2.713−3.35−3.03−2.182−1.801−1.933−1.601−1.766−1.386−0.37−0.090.9550.6040.7580.2280.4820.0040.340.02These findings are confirmed by the transient absorption spectroscopy measurement
which find that the hole injection process on PPV−Ru1 occurs within
ps which is significantly shorter that the hole injection time
ps for PPV−Ru2. Furthermore, PPV−Ru1 also has a longer lifetime of the charge separated state (
fs) than PPV−Ru2 (
ps).
This comparison with the experiment indicates that the analysis of
and
indeed allows to qualitatively predict the hole injection and charge recombination properties of a conjugated polymer functionalized with Ru photosensitizers.Still, such a photocathode system is usually not operated under gas‐phase conditions but in the presence of an electrolyte.
In fact, it is well‐known that the presence of an electrolyte can significantly modify the electronic properties of photoactive complexes.[
,
] In a first step, we calculated the electronic properties of the isolated PPV oligomers (monomer, dimer, and trimer) and the Ru2 dye separately in the presence of implicit solvents. In the case of the dye, we added two PF
counter ions in order balance the positive charge of the Ru dye. The HOMO, LUMO and the energy gap of the considered systems are plotted as a function of the dielectric constant of the solvent in Figure 5. In general, we find a relatively weak dependence of these electronic properties of the PPV oligomers on the dielectric constant. There is a slight downshift of the energy levels of the HOMO and the LUMO with increasing dielectric constant, but the absolute shifts are below 0.25 eV. On the other hand, introducing implicit solvents modifies the HOMO and LUMO levels of the Ru2 dye complex significantly. Both the HOMO and the LUMO are shifted up by up to 1 eV with increasing dielectric constant, but in such a way that the HOMO‐LUMO gap after an initial sharp rise then decreases again. Thus the electronic levels of the dye are obviously more sensitive with respect to the presence of a solvent than the electronic levels of the polymer.
Figure 5
Energy levels of the HOMO (▴), the LUMO (▾), and the energy gap (•) of isolated PPV oligomers and the Ru2 dye as a function of the dielectric constant of an implicit solvent they are embedded in.
Energy levels of the HOMO (▴), the LUMO (▾), and the energy gap (•) of isolated PPV oligomers and the Ru2 dye as a function of the dielectric constant of an implicit solvent they are embedded in.As the FHI‐AIMS code used in this study does not allow to include implicit solvents in its implementation for periodic systems, we tried to estimate the effect of the presence of a solvent by adding
and 30 explicit water molecules to the simulation cell and relaxing them. Note that by this procedure the statistical nature of liquid water is not properly accounted for, however, performing a statistical average over many possible water structures would be much too time‐consuming. The calculated energy levels and the hole and electron injection efficiences without and with water molecules are presented in Table 2. We also compare the calculated values with the experiment.
and
are estimated from the experiment by analyzing the cyclic voltammograms in Ref. [15], and the HOMO and LUMO levels are estimated from the reported oxidation and reduction levels, as done before.[
,
,
]According to Table 2, the addition of the water molecules hardly changes the local HOMO and LUMO levels on the PPV polymer which is in nice agreement with the calculations including implicit solvents presented in Figure 5a. However, the HOMO and LUMO levels are shifted up by about 0.6 eV upon adding 30 water molecules to the simulation cell, again in good agreement with the implicit solvent results (Figure 5b). Consequently, also the values of the hole and electron injection efficiency
and
become modified, both are reduced by about 0.4 eV upon adding 30 explicit water molecules.As far as the comparison between calculated and measured values listed in Table 2 is concerned, a qualitative and almost semi‐quantitative agreement is observed. This is rather satisfying given first the uncertainties of the derivation of the numbers from the experimental results. Second, the regular periodic structure of the PPV−Ru system illustrated in Figure 2b certainly represents an idealization with respect to the periodicity of the attached photosensitizer, the assumption of a straight polymer, the exact length of the linker between polymer and photosensitizer and the representation of the aqueous environment by a finite number of explicit water molecules in just one structure. Still, our DFT calculations yield a HOMO‐LUMO gap of the Ru photosensitizer that is by roughly 1 eV larger than the experimentally derived value which has also been found in similar systems.[
,
]Time‐dependent density functional theory, on the other hand, is well established to access and predict the light‐driven processes in transition metal complexes[
,
] and soft matter embedded light‐harvesting units.
Through TDDFT, individual excited states and their properties such as oscillator strength and electron density become accessible, which allows fundamental insights into the mechanisms of light‐driven processes such as electron or hole injection. Further, orbital‐based Kohn‐Sham TDDFT has the advantage of providing easy means of visualization in the familiar orbital picture or with charge density difference plots. It is an excellent method to couple with the experiment: Actual measurements help verify calculated data, while the calculation can access properties that are difficult to ascertain for the experimentalist, such as the presence of optically inaccessible “dark” states, such as the desired charge‐transfer states in the present hybrid material. Therefore, as a complementary approach, we address the PPV−Ru complexes also by orbital‐based TDDFT calculations using local models for the complexes.A first comparison between TDDFT calculations and experiment is given by the comparison of the TDDFT spectra with experimentally measured ones. Figure 6 shows calculated and experimental spectra for PPV−Ru1 and PPV−Ru2 side by side. The spectra are in good agreement, especially considering the simplified geometries, omitted solvent and computational concessions towards the basis size. Table 3 lists the most prominent excitations in detail, while their electronic nature is visualized by charge density difference (CDD) plots, see Figure 6.
Figure 6
Visualization of calculated UV/Vis absorption spectra (solid) and excitations for PPV−Ru1 (a) and PPV−Ru2 (b). Experimental spectra (dashed) are adapted from Ref. [15] and have been scaled arbitrarily. Calculated transitions have been broadened to a full width at half maximum of 0.5 eV. Charge density difference plots (isovalue=0.002) for some important transitions are included. Holes are depicted in orange, electrons in blue. Calculated data was obtained at the PBE0/def2‐SVP level of TDDFT. Energies and oscillator strengths can be found in Table 3.
Table 3
Selected excitations contributing to UV/Vis absorption spectra of PPV−Ru1 and PPV−Ru2 as well as (dipole‐forbidden) excitations involved in charge separation. Experimental values (measured in CHCl3) are adapted from ref. [15]. Calculated data obtained at the PBE0/def2‐SVP level of TDDFT.
Compound
Excitation
Type
fosc
Ecalc/eV
Eexp/eV
λcalc
/nm
λexp
/nm
PPV−Ru1
S1
PPV→Dye
0.00
1.45
–
855
–
S2
PPV→Dye
0.00
1.57
–
789
–
S3
PPV→Dye
0.00
1.82
–
680
–
S4
PPV→Dye
0.00
1.86
–
666
–
S5
PPV→Dye
0.00
1.93
–
642
–
S6
PPV→Dye
0.00
2.05
–
604
–
S7
PPV→Dye
0.00
2.25
–
550
–
S8
PPV→Dye
0.00
2.30
–
540
–
S9
PPV→Dye
0.00
2.31
–
537
–
S10
PPV→Dye
0.00
2.34
–
530
–
S12
PPV→Dye
0.00
2.42
–
512
–
S13
π→π*
3.20
2.47
2.56
502
485
S25
MLCT
0.17
2.81
2.60
441
476
S35
MLCT
0.20
3.04
2.83
408
438
S68
π→π*
0.41
3.67
4.01
338
309
PPV−Ru2
S1
PPV→Dye
0.00
1.62
–
766
–
S2
PPV→Dye
0.00
1.82
–
681
–
S3
PPV→Dye
0.00
2.11
–
588
–
S4
PPV→Dye
0.00
2.16
–
573
–
S5
PPV→Dye
0.00
2.21
–
562
–
S6
PPV→Dye
0.00
2.31
–
537
–
S8
π→π*
3.27
2.48
2.58
500
481
S16
MLCT
0.09
2.82
2.71
440
458
S23
MLCT
0.12
3.01
2.88
412
430
S51
π→π*
0.63
3.68
3.79
337
327
Visualization of calculated UV/Vis absorption spectra (solid) and excitations for PPV−Ru1 (a) and PPV−Ru2 (b). Experimental spectra (dashed) are adapted from Ref. [15] and have been scaled arbitrarily. Calculated transitions have been broadened to a full width at half maximum of 0.5 eV. Charge density difference plots (isovalue=0.002) for some important transitions are included. Holes are depicted in orange, electrons in blue. Calculated data was obtained at the PBE0/def2‐SVP level of TDDFT. Energies and oscillator strengths can be found in Table 3.Selected excitations contributing to UV/Vis absorption spectra of PPV−Ru1 and PPV−Ru2 as well as (dipole‐forbidden) excitations involved in charge separation. Experimental values (measured in CHCl3) are adapted from ref. [15]. Calculated data obtained at the PBE0/def2‐SVP level of TDDFT.CompoundExcitationTypef
oscE
calc/eVE
exp/eV/nm/nmPPV−Ru1S1PPV→Dye0.001.45–855–S2PPV→Dye0.001.57–789–S3PPV→Dye0.001.82–680–S4PPV→Dye0.001.86–666–S5PPV→Dye0.001.93–642–S6PPV→Dye0.002.05–604–S7PPV→Dye0.002.25–550–S8PPV→Dye0.002.30–540–S9PPV→Dye0.002.31–537–S10PPV→Dye0.002.34–530–S12PPV→Dye0.002.42–512–S133.202.472.56502485S25MLCT0.172.812.60441476S35MLCT0.203.042.83408438S680.413.674.01338309PPV−Ru2S1PPV→Dye0.001.62–766–S2PPV→Dye0.001.82–681–S3PPV→Dye0.002.11–588–S4PPV→Dye0.002.16–573–S5PPV→Dye0.002.21–562–S6PPV→Dye0.002.31–537–S83.272.482.58500481S16MLCT0.092.822.71440458S23MLCT0.123.012.88412430S510.633.683.79337327Both experimental and calculated spectra are dominated by a peak at ca. 2.5 eV, which stems from a
transition localized on the PPV polymer. This excitation energy is in rather good agreement with the band gap of the extended polymer states derived from the periodic DFT calculations (see Figure 3). The intensity of this peak is possibly overestimated by the calculation due to the linear arrangement of the phenylenevinylene chain in the model system which amplifies the transition dipoles in the aromatic plane. The Ru(II) dye is hardly involved in this transition at all, as already identified in the analysis of the band structure, and it is also observed in the measured spectrum of pure PPV‐related compounds.PPV−Ru1 as well as PPV−Ru2 feature a variety of moderately intense MLCT bands in the energetic region around 3 eV – overlapping with the bright
transition in PPV. Therefore, a precise assignment of the predicted MLCT transitions with respect to the experimental data is hampered. Yet, these MLCT transitions are likely to be of fundamental importance for the charge separation process. Since these transitions are completely localized on the ruthenium photosensitizer, some insight can instead be gained from experimental measurements of the unfunctionalized Ru1 and Ru2 complexes. Compared to those values, calculated results are systematically blue‐shifted by 0.1 eV to 0.2 eV which can likely be traced to the influence of the CHCl3 solvent used in the experiment.Furthermore, an additional
absorption feature is predicted around 4 eV for both molecules, outside the visible spectrum. Therefore, and not surprisingly, the absorption spectra of PPV−Ru1 and PPV−Ru2 are governed by the separate photophysical properties of the polymer and the Ru‐based photosensitizer, while any significant light‐driven interactions between the polymer and the photosensitizer – i. e. strongly dipole‐allowed CT states – are not observed. However, dipole‐forbidden CT transitions would have negligible intensity in an experimental measurement and are optically inaccessible from the ground state and only experimentally detectable based on transient absorption techniques. Of particular interest is a series of dark charge transfer excitations present in both hybrid systems at slightly differing energies, starting from 1.45 eV in PPV−Ru1 and from 1.62 eV in PPV−Ru2 which are characterized by very low oscillator strengths. All of these states feature a distinct charge‐separated nature, where an electron has been transferred from the PPV backbone to the ligand sphere of the respective ruthenium complex.The TDDFT results provide a conclusive picture of possible charge separation processes. As all charge‐separated states in the low energy region are dark, initial excitation must occur into a higher excited state such as via
(PPV−Ru1: S8; PPV−Ru2: S13) or 1MLCT (PPV−Ru1: S25, S35; PPV−Ru2: S16, S23) excitation in the visible region between 2.5 eV and 3 eV. From these states, the molecule can relax into a charge‐separated state through internal conversion and/or intersystem crossing.However, the pronounced excited state lifetimes of 3MLCT states – rapidly accessible via intersystem crossing – suggests that the most important pathway to charge separation occurs through MLCT transitions. Experiments yield lifetimes of the charge‐separated states of (
ps) for PPV−Ru1 and (
ps) for PPV−Ru2.
At these timescales, it is most likely that the charge‐separated states have triplet character, which delays the spin‐forbidden transition back towards the singlet ground state.Both PPV−Ru1 and PPV−Ru2 show certain similarities but also some characteristic differences in the TDDFT calculations. First, excitations centered on the aromatic polymer chain are hardly influenced by the presence of linker and photosensitizer, which is consistent with other findings both here (see, e. g., Figure 3) and elsewhere.
Energies of the MLCT bands in the two compounds are also similar, in spite of the fact that oscillator strengths are higher for PPV−Ru1. The main differences lie in the low‐lying dark states with charge transfer character. The S1 state is lower in PPV−Ru1 by ca. 0.2 eV than in PPV−Ru2, and similar findings hold for other states of this type. This trend is most easily attributed to the electron‐withdrawing effect of the −COOH substituents on the photosensitizer. As a result, PPV−Ru1 features a higher number of low energy charge‐separated states which are accessible upon internal conversion. Although the time scale of charge separation cannot be directly deduced from the density of states, we still speculate that this finding might also allow to explain why charge separation occurs faster in PPV−Ru1 (
ps) than in PPV−Ru2 (
ps).As highlighted above, the TDDFT predictions are in excellent agreement with both the results of the ground‐state periodic DFT calculations and the experiment data. In particular, the ability to assess the overshadowed MLCT transitions as well as to elucidate dark charge‐separated transitions is a valuable tool in predicting and analyzing the photophysics and subsequent excited state relaxation pathways in the present dye‐sensitized polymers.
Conclusion
Using both periodic density functional theory (DFT) calculations and local time‐dependent DFT calculations, we have addressed the coupling of photoactive Ru(II)‐based complexes with two different substituents (−CH3 vs. −COOH) to a polymer matrix. Conducting polymers represent a promising alternative to semiconducting oxide electrodes typically used in dye‐sensitized cathodes. Whereas the electronic properties of the Ru dye depend on the choice of the substituents, those of the polymer are insensitive to the presence of the attached dye. The electronic states of the polymer and the dye do not couple, hence the calculated photoinduduced electronic excitations, which agree well with experimental results, remain localized either on the polymer or the dye. However, we identify possible routes towards photoinduced excitations that would lead to a charge transfer from the polymer to the dye which might then for example be used to induce catalytic reactions at suitable metal centers attached to the dye.
Conflict of interest
The authors declare no conflict of interest.As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.Supporting InformationClick here for additional data file.
Authors: Stephan Kupfer; Julien Guthmuller; Maria Wächtler; Sebastian Losse; Sven Rau; Benjamin Dietzek; Jürgen Popp; Leticia González Journal: Phys Chem Chem Phys Date: 2011-07-26 Impact factor: 3.676
Authors: Ruri A Wahyuono; Bianca Seidler; Sebastian Bold; Andrea Dellith; Jan Dellith; Johannes Ahner; Pascal Wintergerst; Grace Lowe; Martin D Hager; Maria Wächtler; Carsten Streb; Ulrich S Schubert; Sven Rau; Benjamin Dietzek Journal: Sci Rep Date: 2021-02-02 Impact factor: 4.379