Davide
Moia
*a,
Valérie
Vaissier
a,
Ismael
López-Duarte
b,
Tomás
Torres
bc,
Mohammad K.
Nazeeruddin
d,
Brian C.
O'Regan
e,
Jenny
Nelson
a and
Piers R. F.
Barnes
*a
aDepartment of Physics, Imperial College London, SW7 2AZ, UK. E-mail: davide.moia11@imperial.ac.uk; piers.barnes@imperial.ac.uk
bDepartamento de Química Orgánica, Facultad de Ciencias, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
cIMDEA-Nanociencia, c/Faraday, 9, Campus de Cantoblanco, E-28049 Madrid, Spain
dLaboratory of Photonics and Interfaces, Department of Chemistry and Chemical Engineering, Swiss Federal Institute of Technology, Station 6, CH-1015 Lausanne, Switzerland
eDepartment of Chemistry, Imperial College London, SW7 2AZ, UK
First published on 18th October 2013
We measured the rate of hole hopping between dye molecules on titanium dioxide nanocrystals using cyclic voltammetry. Dyes commonly used in the field of dye sensitized solar cells exhibited efficient intermolecular charge transport, showing apparent diffusion coefficient values between 10−8 up to over 10−7 cm2 s−1 at room temperature. From temperature dependent measurements, we observed that hole transport across dye monolayers is a thermally activated process with Arrhenius activation energies between about 170 and 370 meV depending on the dye. Analysis of the data in terms of non-adiabatic Marcus theory of charge transfer enabled the estimation of the reorganization energy (740–1540 meV) and of an effective electronic coupling for the different systems. The measured reorganization energies show reasonable agreement with values obtained from density functional theory based calculations, validating our computational approach. Finally, we interpret the experimental and calculated data with reference to the chemical structure of the dyes and to the packing of the dyes on the surface of the TiO2 and suggest that delocalization of the HOMO and rigidity of the conjugated molecular structure result respectively in lower outer and inner sphere reorganization energies.
The first demonstration of lateral conduction in a monolayer of dye molecules in a DSSC configuration was shown by Bonhôte et al.12 and was followed by several pieces of work where either hole or electron hopping was demonstrated for several classes of molecules.10,13–20 In these studies, electrochemical, spectroelectrochemical and transient absorption anisotropy measurements were used to assess the diffusion of the charges across the dye sensitized surface, and the observations were discussed in relation to the chemical structure of the molecules and the kinetics of the charge transfer step.
Non-adiabatic Marcus theory is commonly used to express the charge transfer rate, Γ, between a weakly coupled electron donor and electron acceptor as a function of the electronic coupling, J, the reorganization energy, λ, and the difference in free energy between the equilibrium states of products and reactants, ΔG0 = Gproducts − Greactants:21
(1a) |
In the particular case of a self exchange reaction, such as hole hopping between two identical dye molecules anchored to a surface, we can set ΔG0 = 0, if we assume negligible energetic disorder in the system (eqn (1b)). Then
(1b) |
The reorganization energy λ is a measure of the energetic cost of rearranging the molecules' conformation and the surrounding medium upon the charge transfer step. Electrochemistry as well as transient absorption spectroscopy and Raman spectroscopy have been employed to experimentally determine λ for non-adiabatic charge transfer processes referring to eqn (1a) and (1b).22–26 In these studies, heterogeneous charge transfer between an electrode and molecules either in solution or anchored to the electrode surface and homogeneous charge transfer reactions between molecules in solution have been considered. To our knowledge, no measurement of λ for charge transfer between dyes anchored to a surface has previously been reported.
In the case of transport between dyes anchored to a surface, the rate of intermolecular charge transfer between dyes (Γ) is expected to be related to the apparent diffusion coefficient of holes hopping across the surface of a dye sensitized film (Dapp). Dapp can be obtained from electrochemical measurements such as cyclic voltammetry of the films (see Fig. 1a–d).16 Thus, by examining the dependence of Dapp on temperature for a monolayer of dyes, we can use eqn (1b) to derive estimates for λ for thermally activated diffusion of holes across monolayers (Fig. 1e and f). The values of λ can then be compared to those that we obtain using a computational method recently proposed by our group.27
Fig. 1 (a and b) Cyclic voltammetry measurement of dye sensitized titanium dioxide mesoporous films on FTO in an inert electrolyte. When the positive voltage ramp (VWE) is applied to the sample (working electrode of the electrochemical cell), holes injected from the FTO electrode in the HOMO level of the dyes travel across the dye sensitized surface. (c) Fraction of oxidized dyes in the film (concentration of holes, ch, divided by the total concentration of dyes, c0) as a function of position x in the film. The curves are calculated at different points of the cyclic voltammogram using the model described in the Section 2.1 of the ESI.†ch varies throughout the forward (solid lines) and the reverse (dashed lines) voltage scans following Fick's diffusion equation. The boundary conditions are set by the potential VWE at the FTO electrode (x = 0) and the finite thickness of the film. (d) The calculated current density j flowing across the dyes during the cyclic voltammogram is related to the gradient of ch at the FTO/dyed-TiO2 interface. (e and f) We describe the individual charge transfer process with non-adiabatic Marcus theory and estimate the value of the reorganization energy λ which refers to the energetic cost of rearranging the dye molecules' conformation and the surrounding electrolyte. |
The method that we present can be applied to any redox active molecule, whose oxidized state is stable within the timescale of the measurement. The molecule also needs to be able to chemically anchor to a surface. Herein, we examine dye molecules, widely applied in DSSCs, that present a diverse range of chemical structures (Fig. 2) sensitizing films of TiO2 nanoparticles. TiO2 is a convenient scaffold for the formation of dye monolayers and does not participate in the hole hopping process. Note that other inert insulating materials can be used.9,12 The ruthenium bipirydine complexes Z907, N719 were chosen since they represent reference dyes in the field of DSSCs. The dye N820 was also considered to compare the effect of removing the long alkyl chains from the structure of Z907. Organic indoline dyes D131 and D149 have been successfully employed as sensitizers in DSSCs, especially in thinner film devices, on account of their high extinction coefficient. Phthalocyanines (Pc) represent another important family of organic semiconductors, used as sensitizers but also as absorbers and hole transport materials in organic electronic devices.28–33 A series of zinc (TT-1) and ruthenium (TT-35, A2 and A5) phthalocyanines were investigated. TT-1 (ref. 34) and TT-35 (ref. 35) complexes share the same peripheral substitution pattern based on the presence of tert-butyl groups and the carboxy anchoring group directly attached to the Pc ring. In the case of ruthenium phthalocyanines A2 and A5,36 the anchoring group is on one of the axial pyridine ligands and is therefore orthogonal to the Pc plane. The other axial pyridine ligand bears an electron-donor group, namely triphenylamine (A2) and methoxy (A5). The dye zinc protoporphyrin(IX) (ZnPP) was also tested since interest in porphyrin sensitizers has been increasing after high efficiency values for liquid DSSCs were achieved using a molecule from this class.37 Some of these molecules have been already studied as efficient hole transporters in dye sensitized systems.16,17,20,38
The reorganization energy of charge transfer between dyes was calculated using a standard method39 which was modified to incorporate the effect of a solvent.27 First, the geometry of both the hole donor (oxidized dye) and acceptor (unoxidized dye) were optimized with acetonitrile as the surrounding medium using DFT coupled with a polarizable continuum model (PCM). The cavity used in the PCM is made of interlocking spheres centered on each atom, offering a specific mapping of each molecule.40 In the PCM used to calculate these equilibrium geometries, the solvent was characterized by its static dielectric constant which we modified to account for the presence of electrolyte ions. The energy of the molecule in its surrounding medium was calculated under four conditions. The first pair of calculations give the equilibrium energy of the system: the energy of the unoxidized and oxidized states with their surrounding medium in the equilibrium geometries, as described above. Then the energies of the oxidized state in the unoxidized equilibrium geometry and unoxidized state in the oxidized equilibrium geometry were calculated. This second pair of calculations represent the energy of the system under non-equilibrium conditions immediately after the moment of charge transfer. For these non-equilibrium conditions the solvent in the PCM model was described using its optical dielectric constant. The reorganization energy is then calculated by subtracting the energy of the equilibrium states from the non-equilibrium states. This easily accessible method, which uses the Gaussian09 software package,41 accounts for the effects of the surrounding electrolyte (solvent and ions) on the reorganization energy. By repeating the calculations in a vacuum, we can discriminate between the inner sphere reorganization energy (λi), due to the rearrangement of the dyes upon hole transfer, from the outer sphere component (λo), related to the electrolyte.
We consider cases where the penetration depth of the holes into the film is sufficiently short relative to the film thickness at the time corresponding to the oxidation peak. If this condition is valid, the apparent diffusion coefficient of the holes, Dapp [cm2 s−1], can be determined from the peak current density of the cyclic voltammogram, jp [A cm−2], using the following relationship:
(2) |
The Arrhenius plots in Fig. 4a–c show the values of Dapp for different dyes plotted as a function of temperature. They also show lines corresponding to the least square fits of the expression of a thermally activated process to the experimental data. The data quality is sufficient to discern differences in the values of the activation energy (Eact) between some of the systems (Table 1). Fig. 4d shows the values of reorganization energy estimated from these data as described in the following paragraphs where it will be discussed further.
Fig. 4 Temperature dependent values of the hole apparent diffusion coefficient (Dapp) extracted from cyclic voltammetry performed on dye sensitized TiO2 films on FTO in 0.1 M TBAP in acetonitrile. For each dye, we show the fitting of the expression of a thermally activated process Dapp = D0exp(−Eact/kBT) to the data, which enables the estimation of the activation energy Eact related to the charge transfer between dyes. Three classes of dyes were examined: ruthenium bipyridine complexes (a), indoline dyes (b) and phthalocyanines (c). In (d) we show the values of the reorganization energy λmeas that we obtained by fitting the data using eqn (1b) and (3) plotted against Dapp at room temperature for all the dyes. Dashed lines represent the relation between λ and Dapp for different values of the effective electronic coupling Jeff and for a value of the mean intermolecular distance R of 1.35 nm (average over all the dyes considered). |
Dye | c 0 [1019 cm−3] | D app at 300 K [10−8 cm2 s−1] | E act [meV] | J eff [meV]a | λ meas [meV] | λ calc [meV] |
---|---|---|---|---|---|---|
a ×/ Expresses the geometric error. | ||||||
Z907 | 9.1 ± 1.4 | 2.0 ± 0.5 | 274 ± 19 | 6.0 ×/1.5 | 1143 ± 76 | 976 |
N820 | 9.6 ± 1.5 | 3.5 ± 0.9 | 262 ± 25 | 6.4 ×/1.7 | 1095 ± 101 | 979 |
N719 | 6.6 ± 1.2 | — | — | — | — | — |
D131 | 8.8 ± 1.4 | 2.8 ± 0.8 | 373 ± 17 | 51 ×/1.5 | 1538 ± 67 | 1062 |
D149 | 6.6 ± 1.0 | 27 ± 8 | 244 ± 22 | 10 ×/1.6 | 1023 ± 88 | 914 |
ZnPP | 4.2 ± 0.7 | — | — | — | — | — |
TT-1 | 13 ± 2 | 0.39 ± 0.12 | 174 ± 24 | 0.40 ×/1.7 | 742 ± 97 | 710 |
TT-35 | 4.4 ± 0.7 | 3.9 ± 1.1 | 245 ± 18 | 3.2 ×/1.5 | 1026 ± 70 | 680 |
A2 | 7.0 ± 1.1 | 1.7 ± 0.5 | 297 ± 21 | 7.6 ×/1.6 | 1234 ± 83 | 666 |
A5 | 2.7 ± 0.4 | 1.3 ± 0.4 | 281 ± 24 | 3.1 ×/1.7 | 1170 ± 95 | 669 |
Assuming that the temperature dependence of Dapp has the same functional form as the expression of the charge transfer rate shown in eqn (1b), then we are able to determine the reorganization energy describing the charge transfer reaction between the dye molecules anchored on the TiO2 surface. Moreover, considering the simple case of an isotropic arrangement of the dyes on the surface of the TiO2 with constant intermolecular distance R, the relationship describing Dapp can be expressed as:27
(3) |
The results of the fits of eqn (1b) for the values of the effective electronic coupling, Jeff, and of the reorganization energy, λmeas, are shown in Table 1. The values of λmeas are also displayed in Fig. 4d, where they are plotted against the corresponding dye's apparent diffusion coefficient at room temperature. Lines representing the relation between the two quantities according to eqn (1b) and (3) are shown for different values of the effective electronic coupling using the average value of R over the dyes considered of 1.35 nm. Table 1 also shows the estimation of the activation energy, Eact (determined from the Arrhenius equation), the apparent diffusion coefficient at room temperature (Dapp) and the measured dye loading (c0) used for the estimation of Jeff. The measured values of reorganization energy (λmeas) are compared to values calculated via quantum chemical modelling (λcalc) using the method described in ref. 27 for the case of 0.1 molar concentration of ions in acetonitrile as surrounding dielectric medium.
In Fig. 5 we show this comparison between theory and experiments and we display values of λcalc as the sum of the inner sphere reorganization energy components (λi) and the outer sphere ones (λo). The dyes N719 and ZnPP showed highly irreversible behaviour. For these molecules, we did not apply the analysis in terms of hole diffusion. Further details can be found in Section 1.2 of the ESI.†
In the next section we discuss the results in Table 1 and describe the behavior of the dyes in terms of hole diffusion in relation to their chemical structure. We consider the HOMO delocalization, the rigidity and conjugation of the dyes, the presence of specific side groups or ligands, the anchoring mechanism and the configuration of the monolayer to interpret the calculated and measured reorganization energies and effective electronic couplings.
Regarding the calculated inner sphere reorganization energy, we report almost identical values for the ruthenium dyes Z907 and N820. This suggests that the alkyl chains play a negligible role in the rearrangement of the molecular conformation of Z907 upon charge transfer. For the indoline dyes, D131 shows higher λi than D149, possibly due to the reduced size of the conjugated structure. Finally, we observe that the phthalocyanine dyes show very low values of λi. We ascribe this result to the rigidity of the phthalocyanine ring, where the HOMO is mainly located (almost all atoms are braced by bonds to neighboring atoms, see Fig. 2), which constrains the unoxidized and the oxidized conformations to be similar. The planarity of the structure ensures high level of conjugation across the ring that provides easier accommodation for the hole.
Our analysis of the thermally activated charge transfer process with the non-adiabatic Marcus equation is based on reactions between redox sites of identical energy (see eqn (1b)). In reality, we expect that variations in the TiO2 surface, in the configuration of the dyes on the surface and in the dielectric and electrostatic environment result in variations in site energy (energetic disorder). This introduces additional temperature dependence in the charge diffusion process and a distribution of values for the activation energy. If these effects are significant, they could lead to an overestimate of λ in our experimental analysis. This could explain the mismatch between our experiments and calculations.
The close packing of dyes on the TiO2 surface is expected to increase the coupling between molecules and result in faster charge transport. The enhancement of hole diffusion for indoline dyes due to J-aggregation has already been discussed by Fattori et al.17 Moreover, it has been shown that the preferential adsorption mechanism of D149 on the TiO2 surface occurs with lower packing and reduced interaction than for another indoline dye, D102, due to steric effects.44 In our experiments, D131 shows higher packing than D149 (see values of the dye loading c0 in Table 1), and probably adopts out of plane configurations, similar to what has been presented for the indoline dye D102. The difference in packing results in different effective electronic coupling, which is higher in this case for D131 than D149 (see Jeff in Table 1). The faster diffusion that we measure for D149 at room temperature compared to D131 must therefore be due to the difference in reorganization energy.
The phthalocyanine dye TT-1 adsorbs onto the TiO2 with very high packing density. Its slow hole diffusion at room temperature can be ascribed to the presence of the tert-butyl groups or to unfavorable relative orientations of neighboring dyes. The ruthenium phthalocyanines show significantly lower dye loading than TT-1 due to the presence of axial ligands, in the case of TT-35, and to the different anchoring scheme, for A2 and A5, (c0 is particularly low for A5, which we are not able to explain on the basis of the dye's chemical structure). For these dyes, we record apparent diffusion coefficients at room temperature above 10−8 cm2 s−1. All the phthalocyanine dyes investigated here are particularly interesting given the high degree of reversibility we observe when performing multiple CV measurements on the same sample.
The corresponding value of mobility that can be calculated from the measurement of Dapp at room temperature for the dyes investigated in this study ranges between 10−7 and 10−5 cm2 V−1 s−1. These values are lower than the mobilities recorded for solution processable amorphous organic semiconductors (up to 10−2 cm2 V−1 s−1). This may be because of the different dimensionality of a dye monolayer on a mesoporous structure compared to the bulk of a thin film. Mobility greater than 10−2 cm2 V−1 s−1 was reported using solution processable oligothiophenes in self assembled monolayer field effect transistors (SAMFETs), suggesting a potential improvement for the case of dyes that are able to order when sensitizing a surface.11,45 Organic semiconductor monolayers can occupy only about 10% volume fraction of a mesoporous film, and can revolutionize the provision/extraction of charge carriers in the bulk of a film. 2D transport across a semiconductor (Perovskite) sensitized nanocrystal was recently observed by Lee et al. and enabled the fabrication of high efficiency solid state solar cells.46,47
The method that we present has potential application beyond the field of DSSCs for the characterization of charge transport between semiconducting molecules. Traditionally, it has been difficult to resolve the contributions from energetic disorder and reorganization energy from observed temperature dependence of transport in these materials.48 Indeed, there is, as yet, no straightforward technique to determine λ experimentally. The experimental and theoretical methods used here could be applied to organic semiconducting moieties when attached to the surface of a supporting scaffold and used to study the reorganization energy in different dielectric environments. In the surface attached configuration we may expect a lower degree of disorder than in a three-dimensional solid film. However, as mentioned previously in this discussion, variations in the surface and in the configuration of the molecules could still make a significant contribution to the temperature dependence of hole diffusion.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c3sc52359d |
This journal is © The Royal Society of Chemistry 2014 |