International Journal of Photoenergy

Volume 2011 (2011), Article ID 316952, 11 pages

http://dx.doi.org/10.1155/2011/316952

## Theoretical Insight into the Spectral Characteristics of Fe(II)-Based Complexes for Dye-Sensitized Solar Cells—Part I: Polypyridyl Ancillary Ligands

^{1}Department of Physics and Materials Science, City University of Hong Kong, Hong Kong^{2}College of Physics Science and Technology, China University of Petroleum, Dongying, Shandong 257061, China

Received 22 May 2011; Revised 7 July 2011; Accepted 25 July 2011

Academic Editor: Gaetano Di Marco

Copyright © 2011 Xiaoqing Lu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The design of light-absorbent dyes with cheaper, safer, and more sustainable materials is one of the key issues for the future development of dye-sensitized solar cells (DSSCs). We report herein a theoretical investigation on a series of polypyridyl Fe(II)-based complexes of FeL_{2}(SCN)_{2}, [FeL_{3}]^{2+}, [FeL′(SCN)_{3}]^{-}, [FeL′_{2}]^{2+}, and FeL′′(SCN)_{2} (L = 2,2′-bipyridyl-4,4′-dicarboxylic acid, L′ = 2,2′,2^{″}-terpyridyl-4,4′,4^{″}-tricarboxylic acid, L^{″} =
4,4^{‴}-dimethyl-2,2′ : 6′,2^{″} :6^{″},2^{‴}-quaterpyridyl-4′,4^{″}-biscarboxylic acid) by density functional theory (DFT) and time-dependent DFT (TD-DFT). Molecular geometries, electronic structures, and optical absorption spectra are predicted in both the gas phase and methyl cyanide (MeCN) solution. Our results show that polypyridyl Fe(II)-based complexes display multitransition characters of Fe → polypyridine metal-to-ligand charge transfer and ligand-to-ligand charge transfer in the range of 350–800 nm. Structural optimizations by choosing different polypyridyl ancillary ligands lead to alterations of the molecular orbital energies, oscillator strength, and spectral response range. Compared with Ru(II) sensitizers, Fe(II)-based complexes show similar characteristics and improving trend of optical absorption spectra along with the introduction of different polypyridyl ancillary ligands.

#### 1. Introduction

Due to the features of low cost and high conversion efficiency, dye-sensitized solar cells (DSSCs) based on organic/inorganic hybrid materials have attracted extensive attention as an alternative to the conventional Si-based solar cells [1–3]. As crucial light-harvesting elements in DSSCs, the photoexcited sensitizers can inject electrons from their excited states into the semiconductor conduction band, and are recharged by an electrolyte to their initial states [4–6]. Early design of sensitizers had focused on transition metal coordinated complexes because of their unique richness of electronic properties for light absorption as well as photochemical and photophysical behaviors, especially for the photoredox activities [7]. These behaviors were originated from the nature of the coordinated complexes, involving the metal-to-ligand charge transfer (MLCT) character and intervening spin-orbit coupling attached to the presence of metal cations with large atomic numbers [7]. Metal coordinated complexes also facilitate the fine tuning of electronic properties by means of the proper choice of both polypyridyl ancillary ligands and metal cations. Up to now, polypyridyl Ru(II)-based complexes are proven to be the most efficient sensitizers employed in DSSCs. Among them, the tetraprotonated Ru(4,4′-dicarboxy-2,2′-bipyridine)_{2}(NCS)_{2} complex (coded as** N3**) [8] and its doubly protonated analog (coded as** N719**) [9] exhibit outstanding properties, that is, broad visible light absorption spectrum, high photo-, electro-, thermal and chemical stability, and strong adsorption on the semiconductor surfaces. Using **N3** as sensitizers, a solar-to-electric energy conversion efficiency of 10% was first attained for DSSCs. Therefore, **N3** is usually considered as a reference dye for DSSCs and is used as a base for designing other sensitizers. At present, the maximum efficiency of DSSCs is 12.1% using Ru(II)-based complex of **C106 **as a sensitizer [10].

However, considering the high cost and limited availability of Ru, novel dyes using cheaper, safer, and more sustainable metal materials are desirable for next-generation DSSCs. Although Fe and Ru are in the same column of the periodic table, Fe is a first row transition metal and is environmental-friendly and inexpensive. It is therefore of an interest to determine whether Fe can be used to replace Ru in sensitizers to exhibit similar photophysical and photochemical properties. Recent studies on Fe(II)-based complexes have been reported [7, 11–20], in which the geometries and absorption spectra were analyzed. However, better understanding of the electronic structures and spectra characteristics of Fe(II)-based complexes for DSSCs is still scarce. In this work, we present a theoretical study on a series of polypyridyl Fe(II)-based complexes of FeL_{2}(SCN)_{2}, [FeL_{3}]^{2+}, [FeL′(SCN)_{3}]^{-}, [FeL′_{2}]^{2+}, and FeL′′(SCN)_{2} (L = 2,2′-bipyridyl-4,4′-dicarboxylic acid, L′= 2,2′,2′′-terpyridyl-4,4′,4′′-tricarboxylic acid, L′′= 4,4′′′-dimethyl-2,2′ : 6′,2′′ : 6′′,2′′′-quaterpyridyl-4′,4′′-biscarboxylic acid) that follow straightforward design strategies applied to Ru(II) sensitizers, by means of DFT/TD-DFT calculations. The molecular geometries, electronic structures, and spectral characteristics are elaborated in both the gas phase and MeCN solution.

#### 2. Computational Methods

All calculations were performed using the Gaussian 09 program package [21]. The geometrical optimizations for all Fe(II)-based complexes were carried out using the B3LYP exchange correlation functional in conjunction with the pure 6-31G* and mixed DZVP/3-21G* basis sets, in both the gas phase and MeCN solution. The reason why we chose two basis sets was to eliminate the errors caused by improper selection of basis set. In the pure 6-31G* basis set, all atoms were described in the same basis set. In the mixed DZVP/3-21G* basis set, the DZVP set [22] supplemented with three sets of uncontracted pure angular momentum functions was chosen for Fe, and the 3-21G* basis set for nonmetal atoms. For Ru(II)-based complexes, the B3LYP functional in conjunction with the pure 3-21G* and mixed SDD/6-31G* basis sets was employed in the geometrical optimizations. The SDD set along with the quasirelativistic MWB (28 core electrons) pseudopotential [23, 24] was chosen for Ru and the 6-31G* basis set for the remaining nonmetal atoms. For simplicity, B3LYP/DZVP was used to stand for the B3LYP functional in conjunction with the mixed DZVP/3-21G* basis set and B3LYP/SDD for the B3LYP functional in conjunction with the mixed SDD/6-31G* basis set.

TD-DFT calculations were used to investigate the optical properties of the sensitizers. The lowest vertical excitation energy of the system at the ground state geometry can be accurately and efficiently calculated by TD-DFT [25, 26]. The 70 lowest transitions, up to the energy level of at least 4.50 eV, were taken into account in the calculations of the absorption spectra. The solvent effects were evaluated using the nonequilibrium implementation of the conductor-like polarizable continuum model (C-PCM) [27–30]. In particular, optimized structures and solvation energies were computed by a cavity model [31] coupled to C-PCM [30]. This approach provides results very close to those obtained by the original dielectric model for high dielectric constant solvents but is significantly more effective in geometry optimizations and less prone to numerical errors arising from the small part of the solute electron cloud lying outside the cavity [30].

#### 3. Results and Discussion

In the following sections, we first establish the accuracy for the chosen levels of theory. Then, we present the results of molecular geometries, molecular orbitals and electronic structures, electronic excitations, and absorption spectra of Fe(II)-based complexes, as well as the comparison with those of Ru(II)-based complexes.

##### 3.1. Calibration

It is essential to calibrate the quality of the geometries obtained at different levels because the accuracy of the calculations of electronic structures used to produce reliable absorption spectra depends highly on the geometrical parameters. However, since there is no optimized complex which corresponds exactly to a experimental result, the closest complex [FeL_{3}]^{2+} (L = 2,2′-bipyridyl-4,4′-dicarboxylic acid) is chosen for comparison with the experimental complex, [Fe(bpy)_{3}]^{2+} [11, 13]. Table 1 lists the selected theoretical and experimental parameters of the ground state structures, and the corresponding labeling scheme is given in Figure 1.

The bond length of Fe–N in the gas-phase [FeL_{3}]^{2+} is 1.997 Å calculated at the B3LYP/6-31G* level and 1.992 Å at the B3LYP/DZVP level (see Table 1), in good agreement with the solid-state distances of 1.965–1.967 Å [11, 13]. The bond lengths of C–N and C–C in the N1 containing pyridine ring are in the range of 1.347–1.373 Å and 1.390–1.396 Å, compared to the experimental values of 1.338–1.359 Å and 1.358–1.385 Å, respectively. The bite-angles of the bpy-domains, ∠N1–Fe–N2, is 81.3° obtained at the B3LYP/6-31G* level and 81.5° at the B3LYP/DZVP level, very close to the experimental values of 81.5–81.8° [11, 13]. Moreover, a maximum skewing of 2.5° of the calculated interbpy angles (∠N1–Fe–N3 and ∠N1–Fe–N5) is in good agreement with the available experimental data. In addition, previous calculations using the B3LYP functional with pure DZVP [19], TZVP [13], and 6-311 + G**[13] basis sets are also listed in Table 1 for comparison. The deviations are less than 3% between our results and the experimental and previous calculated values. Considering that our study is performed in the gas phase and does not account for the packing and counter-ion effects that are present in the solid state, our calculated geometries match excellently with experimental results. The good agreements demonstrate the reliability of the chosen model and level of theory.

Considering that TD-DFT excitation energies are dependent on the DFT functionals, benchmark calculations of the absorption spectra using the CAM-B3LYP functional are also provided for Fe(II)-based complex of **1** in both the gas phase and MeCN solution. Seen from Figure S1 (in Supplementary Material available online at doi: 10.1155/2011/316952), these two DFT functionals indeed give consistent absorption spectra. These results justify the reliability of our model and levels of theory. In addition, benchmark calculations containing diffusion functions (e.g., 6-31 + G*) were performed, and the results were very similar to those without the diffusion functions. This is in good agreement with our previous studies on Ru(II)-based complexes [24]. Because the results of Fe(II)-based complexes calculated at the B3LYP/6-31G* and B3LYP/DZVP levels are very close to each other, only the B3LYP/6-31G* results will be used for clarification, unless otherwise stated. In the case of Ru(II)-based complexes, the B3LYP/SDD results will be chosen to compare with Fe(II)-based complexes.

##### 3.2. Molecular Geometries

The optimized molecular structures of the polypyridyl Fe(II)-based complexes of FeL_{2}(SCN)_{2}, [FeL_{3}]^{2+}, [FeL′(SCN)_{3}]^{-}, [FeL′_{2}]^{2+}, and FeL′′(SCN)_{2} are shown in Figure 2, together with the atom numbering scheme. Table 2 summarizes the selected geometrical parameters of the fully optimized Fe(II)-based complexes at the B3LYP/6-31G* level and **N3** at the B3LYP/DZVP level in MeCN solution.

Fe coordination complexes in its Fe(II) form have septet ground state with d^{5}s^{1} configuration, which allows it to bind to six neighboring nitrogen ions to form octahedral complexes. The geometry of **1** with a cis arrangement of the thiocyanate is optimized to be a pseudooctahedral C_{2} symmetry, while the geometry of the trans isomer **1′** exhibits an almost D_{2} symmetry, as shown in Figure 2. For** 1** and **1′**, the bond lengths of Fe–N (thiocyanate) are calculated to be almost the same, while the Fe–N (bpy) bond lengths are ~1.98 Å and ~2.02 Å, respectively. The bite angles of the bpy-domains (∠N1–Fe–N2 and ∠N3–Fe–N4) show a deviation of 1.2° between these two isomers (see Table 2). These differences of both bond lengths and angles correspond to the difference that the trans isomer **1′** is 12.90 kcal mol^{−1} higher in the binding energy than the corresponding cis isomer **1**. To evaluate the protonation/deprotonation effect as appearing between **N3** and its analog **N719**, the H^{+} cations on carboxylic groups are substituted by methyl (**1 **→ **2**). This substitution has slight effect on the geometries, especially for the skeleton parameters (see Table 2). Thiocyanate ligands are usually chosen as the electron-donating groups in the most efficient Ru(II)-based sensitizers. However, from a chemical point of view, the thiocyanate seems to be the weakest part [6]. This situation can be improved by substituting thiocyanate in **1** with a dcbpy ligand (**1 **→ **3**). The interbpy angle in **3 **is distorted somewhat, as indicated by the alteration of ∠N2–Fe–N3 (2.8°). Relatively, distortions are very small for both bond lengths and bite angles of **3** with respect to those of **1**.

The optimized geometry of **4** shows almost a C_{s} symmetric arrangement of the ligands, and the octahedral dispositions of the N atoms around the Fe center are distorted due to the different nature of thiocyanate and tricarboxylic terpyridine (tctpy) ligands. The difference in bond lengths of Fe–N (thiocyanate: 1.946–1.954 Å) and Fe–N (bpy: 1.877–1.991 Å) confirms this point (see Table 2). However, the dihedral angle of ∠N1–N2–N3–N4 (−0.1°) mirrors the fact that the involved N atoms are almost in the same plane, quite different from that of the trans isomer **1′**(22.8°). Substituting tctpy for thiocyanates makes the complex **5** a C_{2} symmetry, in which the geometrical parameters including the bond lengths and angles related to two tctpy ligands are exactly the same with each other.

For the dicarboxylic quaterpyridine (dcqpy) containing Fe(II)-based complex, **6**, the Fe–N (dcqpy) bond lengths are in the range of 1.880–2.063 Å calculated at the 6-31G* level in MeCN solution. A comparison of the optimized parameters shows that complex** 6** has shorter Fe–N (thiocyanate) distances with respect to that in **1′**(1.933 versus 1.944–1.945 Å). This may result from the weaker interaction between Fe and a tetradentate ligand with respect to that between Fe and two bidentate ligands, and thus the stronger interaction between Fe and thiocyanate ligands shortens the Fe–N (thiocyanate) distances. Similar phenomena also appear for the Ru(II)-based complexes when comparing dcqpy containing complex **N886** with dcbpy containing complex **N3 **[32].

##### 3.3. Molecular Orbitals and Electronic Structures

Since the crucial electronic excitations occur from the highest occupied molecular orbitals (HOMOs) to the lowest unoccupied molecular orbitals (LUMOs), it is therefore important to form efficient charge-separated states with the HOMOs localized on the donor subunit and the LUMOs on the acceptor subunit. The molecular orbital energy levels and the HOMO-LUMO gaps of Fe(II)-based complexes are depicted in Figure 3. In view of the distributions of frontier orbital energies, similar molecular orbital shapes would be envisaged for the complexes with the same ancillary ligand. Therefore, **1**, **4**, and **6 **are chosen to exhibit the electron distributions of the frontier molecular orbitals, as shown in Figure 4.

For **1** in MeCN solution, the three highest HOMOs (from HOMO to HOMO−2) are a set of quasi degenerate orbitals, originating from the combination of Fe t_{2g} orbitals (d_{xy}, d_{xz}, and d_{yz}) with the 2p orbital and lone pairs of the thiocyanate ligands. The following HOMO−3 shows similar orbital compositions from thiocyanate ligands but without the contributions from the Fe t_{2g} orbitals with respect to the first three HOMOs (see Figure 4(a)). Compared with those in **N3 **[33], the first four HOMOs in Fe(II)-based complexes **1** show very close orbital compositions and percentages. Not only that, Fe also shows similar field splitting in **1** to that of Ru in **N3**, mirrored by the molecular orbital energy levels of HOMOs shown in Figure 3. The four lowest LUMOs (from LUMO to LUMO+3) of **1**, lying at least 2.73 eV about the HOMO in the range of 0.64 eV, result mainly from the combinations of the C and N atoms localized on the dcbpy ligands (see Figure 4(a)). Since the carboxylic group of the dcbpy ligands has a strong electron-withdrawing ability, the LUMOs are mainly localized on the dcbpy ligands. The anchoring groups involved in those LUMOs would facilitate the electron injection from the excited state of the sensitizers into the semiconductor. Here, it would be specially mentioned that six HOMOs and LUMOs are presented but less levels are observed in Figure 3 due to the degeneracy. For instance, HOMO−1, HOMO−4, LUMO, and LUMO+2 of **1** correspond to HOMO−2, HOMO−5, LUMO+1, and LUMO+3 in energies; thus only eight levels can be observed in Figure 3. In fact, this situation also appears in all the other Fe(II)-based complexes (see Figure 3).

It can be seen from Figure 3 that the energy levels of **2** are very close to those of **1**, indicating that the protonated/deprotonation effect is slight for Fe(II)-based complexes in MeCN solution. However, the isomerization effect has larger influence on the frontier molecular orbitals, mirrored by the alteration of both the energy levels and energy gap. The most obvious alteration between **1** and **1′** is the degenerate orbitals of HOMOs. The first set for **1 **is composed of the three highest HOMOs (from HOMO to HOMO-2) in the range of 0.22 eV, while that for **1′** is composed of the two highest HOMOs in the range of 0.06 eV. This points out the difference of the field splitting of Fe between **1′** and **1**. Due to the alteration of orbital levels, the HOMO-LUMO gap of **1′** reduces to 2.61 eV with respect to 2.73 eV in **1**. Substitution of the thiocyanates in **1** with a dcbpy ligand (**1**→**3**) has great influence on the energy levels and the HOMO-LUMO gap, as shown in Figure 3. All the molecular orbital energies of **3** are downshifted but the HOMOs are stabilized much more than the LUMOs. Accordingly, the HOMO-LUMO gap of **3** is increased to 3.46 eV.

For **4** in MeCN solution, the HOMO shows mainly contributions from the combination of Fe center and three thiocyanate ligands, while the LUMO shows contributions from the tctpy ligand, as shown in Figure 4(b). Compared with **1**, the six highest HOMOs of **4** are unshifted and lie in a smaller range of 0.74 eV, while the LUMO is close to that of **1** but the six lowest LUMOs lie in a larger range of 2.84 eV. Accordingly, the HOMO-LUMO gap of **4** is calculated to be 2.43 eV. The smaller HOMO-LUMO gap would help decrease the excitation energy, and this is beneficial for the absorption spectrum to shift towards the longer wavelength of the visible region. Substitution of the thiocyanates in **4** with a tctpy ligand (**4**→**5**) shows the same trend on the energy levels as **1**→**3** substitution. As seen from Figure 3, both the HOMO and LUMO levels of **5 **are lower than those of **4**, and the HOMO-LUMO gap of **5** is 0.89 eV smaller than that of **4**. Obviously, both the energy levels and HOMO-LUMO gaps of the thiocyanate-substituted complexes (**3** and **5**) are less favorable than those unsubstituted complexes (**1** and **4**) with respect to the reference dye **N3**.

For **6** in MeCN solution, the HOMO also shows main contributions from the combination of Fe center orbital and two thiocyanate ligands, while the LUMO shows contributions from the dcqpy ligand, as shown in Figure 4(c). Fe shows slightly different field splitting in **6** compared with that in **1**, which leads to slightly different energy levels of the HOMOs, as shown in Figure 3. The great difference in energy levels of LUMOs between **1** and **6** correspond to the great structural difference in the polypyridyl ancillary ligands of dcbpy and dcqpy. The HOMO-LUMO gap of 6 is calculated to be 2.50 eV.

##### 3.4. Electronic Excitations and Absorption Spectra

Figure 5 presents the simulated absorption spectra of Fe(II)-based complexes at the B3LYP/6-31G* level and **N3** at the B3LYP/SDD level in MeCN solution. Table 3 lists excitation energies, oscillator strengths, and relative orbital contributions of the optical transitions between 350 and 800 nm for **1**, **4**, and **6**, in which only oscillator strengths larger than 0.01 and orbital percentages larger than 10% are presented. The gas phase absorption spectra and transitions information of **1′**, **2**,** 3**, and **5 **are provided in Figure S2 and Table S1, respectively. As seen from Figure S2, the band shapes of spectra obtained in the gas phase are qualitatively similar to those in MeCN solution, but the oscillator strengths are stronger for most spectra in MeCN solution at the longer wavelength range. In order to depict the real environmental, the results obtained in MeCN solution will be used for analysis.

In the energy range investigated, there are four bands that clearly appear with the intensity maxima centered at ~362, ~395, ~442, and ~571 nm for **1** (see Figure 5). At ~362 nm, the band originates from the starting orbitals of HOMO-3/HOMO-6 combination to the arriving orbitals of LUMO+1/LUMO+2 with the oscillator strength of 0.032. As seen from Figure 4(a), the HOMO−3 → LUMO+2 transition shows the interligand transition from the thiocyanate ligands to the dcbpy ligands, a typical LLCT character. At ~395 nm, the HOMO−4 → LUMO (61%) transition shows the strongest oscillator strength of 0.070 involving the Fe(II) center to the dcbpy ligand and interligand transitions, a multitransition character of MLCT and LLCT. As shown in Table 3, the band at ~442 nm of **1** is composed of a series of peaks from 415.6 nm to 445.3 nm. Among them, the peak at 445.3 nm includes the transitions of HOMO/HOMO-3 as the starting orbitals and LUMO/LUMO+3 as the arriving orbitals, exhibiting a multitransition character of LLCT and MLCT. At ~571 nm, a broad absorption band is observed with the peaks ranging from 559.4 nm to 633.5 nm. The first vertical excitation occurs at 633.5 nm , originated from the HOMO → LUMO+1 (87%) transition with the multitransition character of MLCT and LLCT.

As shown in Figure 5, the calculated absorption spectrum of the derivative **2** is quantitatively similar to that of **1** but with a slight blueshift. This corresponds to the minor alterations of the energy levels caused by the slight structural changes. It can also be confirmed by comparing the excitations of **2** with those of **1** in Tables S1 and 3. The absorption spectrum of isomer **1′** is quite different from that of **1** due to the great difference in geometry. As seen from Figure 5, three bands are found to center at ~382, ~457, and ~625 nm for **1′** in the energy range investigated. The bands at ~625 nm are composed of two peaks, involving the HOMO−1 → LUMO transition at 630.8 nm and the combinations of HOMO/HOMO−1 → LUMO/LUMO + 9 transition at 603.6 nm (see Table S1). The oscillator strength of the first vertical excitation of **1′** is stronger than that of **1**, but the band of **1′** at the longer wavelength region is not as broad as that of **1**. Substitution of the thiocyanates with a dcbpy ligand (**1**→**3**) shrinks the effective spectral response range sharply, corresponding to the enlarged HOMO-LUMO gaps. Only one band is observed at ~450 nm in the energy range investigated, originated from the combinations of HOMO−1/HOMO−2 → LUMO+2 transitions. Compared with Ru(II)-based complex **N3**, the absorption spectra of** 1** and its derivative **2 **have similar band shape. However, the bands in the range of 400–600 nm are only two-thirds as much oscillator strength as those of **N3**.

For **4** in MeCN solution, the calculated absorption spectrum is shifted towards the near-infrared region, which can be comparable with the region that the black dye can reach [34].Three bands are observed with the intensity maxima centered at ~472, ~609, and ~704 nm in the energy range investigated (see Figure 5). At ~472 nm, the bands are composed of two peaks, originated from HOMO−1 → LUMO+2 (92%) transition at 471.1 nm and HOMO → LUMO+2 (84%) transition at 473.4 nm. The oscillator strengths reach 0.054 and 0.064, respectively. An analysis of both transitions shows the multitransition character of MLCT and LLCT involving the Fe center to the tctpy ligand and interligand transitions. At ~609 nm, the HOMO → LUMO+1 (96%) transition has the oscillator strength of 0.031, exhibiting a multitransition character of MLCT and LLCT. The band at ~704 nm is composed of two peaks at 693.0 and 766.8 nm. The peak at 766.8 nm corresponds to the first vertical excitation, originated from the HOMO−1 → LUMO transition with the oscillator strength of 0.017. This extension towards the longer wavelength region would be greatly beneficial for the sensitizers to absorb the low-energy photons and thus improves the DSSCs spectral response range. The substitution of the thiocyanates with a tctpy ligand (**4**→**5**) renders the absorption spectrum of **5** to be redshifted sharply, similar to the effect of substitution of the thiocyanates with a dcbpy ligand (**1**→**3**). The longest wavelength region of the bands for **5 **is at ~543 nm, composed of two peaks from the transitions of HOMO → LUMO and HOMO → LUMO + 1.

For **6** in the energy range investigated, three separated absorption bands are found to center at ~396, ~552, and ~699 nm (see Figure 5). At ~396 nm, the band is composed of the transition of HOMO → LUMO + 3 . At ~552 nm, the band is mainly composed of transitions in the range of 493.4–561.5 nm, originated from the starting orbitals of the HOMO/HOMO + 1/HOMO+2/HOMO+4 combinations to the arriving orbitals of LUMO/LUMO + 1/LUMO+2/LUMO+8 (see Table 3). An analysis of these transition shows that excitations correspond to electron transfer from the Fe 3d orbitals to the unoccupied orbitals of the dcqpy ligand and from the occupied orbitals of the thiocyanate ligands to the unoccupied orbitals of the dcqpy ligand, that is, a mixture of MLCT and LLCT transitions. At ~699 nm, the band is mainly composed of the transitions of HOMO → LUMO (85%) at 682.4 nm and HOMO−1 → LUMO (93%) at 724.5 nm , showing primary MLCT and LLCT characters (see Figure 4(c)). Compared with **1**, the absorption spectrum of dcqpy containing complex **6** is blueshifted towards the longer wavelength range, similar to the phenomenon that occurs for Ru(II)-based sensitizers when going from **N3** to dcqpy containing complex **N886 **[32].

##### 3.5. Fe(II)-Based Complexes versus Ru(II)-Based Complexes

Benefiting from the introduction of different polypyridyl ancillary ligands, the performance of Ru(II)-based complexes has been improved in some respects. For example, when going from dcbpy containing **N3** to tctpy containing black dye, the spectral response range is highly extended and thus the efficiency of DSSCs is enhanced [34]. How do these ligands affect the performance of Fe(II)-based complexes and to what extent can it be confirmed by comparing these complexes with different polypyridyl ancillary ligands?. The key factors for confirmation include energy levels, absorption bands, and oscillator strengths. The energy levels of a sensitizer would determine the electron transfer processes which have great influence on the performance of the DSSCs [35]. As shown in Figure 3*, ***1**,** 4**, and **6** provide suitable controls of the energetic and redox properties, in which the energy levels of the HOMO and LUMO match well with the redox potential of and the conduction band of TiO_{2}. The HOMO-LUMO gap is in sequence of **4** (2.43 eV) < **6** (2.50 eV) < **1 **(2.73 eV), that is, the HOMO-LUMO gap decreases by 0.30 and 0.23 eV when going from **1** to **4** and **6**. It means that **4** and **6 **have values of HOMO-LUMO gap much closer to that of **N3** with respect to the complex **1**. In addition, as seen from Table 3, the lowest vertical excitation energies are also in the sequence of **4** (766.8 nm) < **6** (724.5 nm) < **1 **(633.5 nm), corresponding to the sequence of the HOMO-LUMO gaps. Both facts indicate that the absorption bands are tuned to longer wavelengths when going from **1** to **4** and **6**, thus enhancing its spectral response to the near-infrared light.

On the other hand, the improving trend of Fe(II)-based complexes caused by using the different polypyridyl ancillary ligands can also be shown by comparing with the corresponding Ru(II)-based complexes. The selected absorption spectral properties of Fe(II)- and Ru(II)-based complexes are listed in Table 4. The maximum absorption bands of Fe(II)- and Ru(II)-based complexes show similar trend along with the introduction of different polypyridyl ancillary ligands. As shown in Table 4, the maximum absorption band of **1** centered at 571 nm is redshifted with respect to the corresponding values of **4** (704 nm) and **6 **(699 nm), in good agreement with the alterations from **N3 **(534 nm) to black dye (605 nm) and **N886** (637 nm). Furthermore, the relative oscillator strength of the maximum absorption band, , is calculated to be 2.42 at the B3LYP/6-31G* level. This value is in line with the experimental value of 1.90 obtained by comparing the molar absorption coefficient of **N3** at 534 nm (1.42 × 10^{4} M^{−1} cm^{−1}) [8] with that of black dye at 605 nm (0.748 × 10^{4} M^{−1} cm^{−1}) [34]. This means that the Fe(II)- and Ru(II)- based complexes show similar relation between the structural optimizations and the spectral properties. At present, although the efficiencies using Fe(II)-based complexes are not comparable with those using the Ru(II) sensitizers, the greater abundance of Fe can offer considerable opportunities for reduction in cost. Considering that Fe(II)-based complexes exhibit optical properties and trend similar to Ru(II) sensitizers and still do not undergo the barrage of investigations, it can be envisaged that Fe(II)-based complexes would be competent sensitizers for DSSCs.

#### 4. Conclusions

We have performed DFT/TD-DFT calculations to investigate the molecular geometries, electronic structures, and spectral characteristics of Fe(II)-based complexes in both the gas phase and MeCN solution. The good agreement between the calculated results with the experimental observations ensures the reliability of the chosen theoretical methods. Our results show that the first three HOMOs originate from the combination of Fe 3d orbitals with the 2p orbital and lone pairs of the thiocyanate ligands, while the first four LUMOs are composed of the antibonding combinations of the C and N 2p orbitals. The protonated/deprotonation effect is slight, while the isomerization effect and the solvent effect have great influence on the spectra characteristics of Fe(II)-based complexes. The spectra show ligand-to-ligand charge transfer (LLCT) character at shorter wavelength region, whereas the spectra display multitransition characters of Fe → polypyridine metal-to-ligand CT (MLCT) and LLCT transitions at the longer wavelength of the visible region. Structural optimizations by choosing different polypyridyl ancillary ligands lead to alterations of the molecular orbital energies, oscillator strength, and spectral response range. Compared with Ru(II) sensitizers, Fe(II)-based complexes show similar characteristics and improving trend of optical absorption spectra along with the introduction of different polypyridyl ancillary ligands.

#### Acknowledgments

This work was supported by a Grant from City University of Hong Kong (Project no. 7002657) and Centre for Applied Computing and Interactive Media (ACIM). The authors would also like to thank National Natural Science Foundation of China (10979077, 21003158) and the Fundamental Research Funds for the Central Universities (09CX05002A) for the financial support.

#### References

- B. O'Regan and M. Grätzel, “A low-cost, high-efficiency solar cell based on dye-sensitized colloidal TiO
_{2}films,”*Nature*, vol. 353, no. 6346, pp. 737–740, 1991. View at Google Scholar · View at Scopus - M. Grätzel, “Photoelectrochemical cells,”
*Nature*, vol. 414, no. 6861, pp. 338–344, 2001. View at Publisher · View at Google Scholar · View at Scopus - M. Grätzel, “Solar cells to dye for,”
*Nature*, vol. 421, no. 6923, pp. 586–587, 2003. View at Publisher · View at Google Scholar · View at Scopus - G. Boschloo and A. Hagfeldt, “Characteristics of the iodide/triiodide redox mediator in dye-sensitized solar cells,”
*Accounts of Chemical Research*, vol. 42, no. 11, pp. 1819–18, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Grätzel, “Recent advances in sensitized mesoscopic solar cells,”
*Accounts of Chemical Research*, vol. 42, no. 11, pp. 1788–1798, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Hagfeldt, G. Boschloo, L. Sun, L. Kloo, and H. Pettersson, “Dye-sensitized solar cells,”
*Chemical Reviews*, vol. 110, no. 11, pp. 6595–6663, 2010. View at Publisher · View at Google Scholar · View at Scopus - P. P. Lainé, S. Campagna, and F. Loiseau, “Conformationally gated photoinduced processes within photosensitizer-acceptor dyads based on ruthenium(II) and osmium(II) polypyridyl complexes with an appended pyridinium group,”
*Coordination Chemistry Reviews*, vol. 252, no. 23-24, pp. 2552–2571, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. K. Nazeeruddin, A. Kay, I. Rodicio, R. Humphry-Baker, E. Müller, and P. Liska, “Conversion of light to electricity by cis-X
_{2}bis(2,2'-bipyridyl-4,4'-dicarboxylate)ruthenium(II) charge-transfer sensitizers (X = Cl^{-}, Br^{-}, I^{-}, CN^{-}, and SCN^{-}) on nanocrystalline TiO_{2}electrodes,”*Journal of the American Chemical Society*, vol. 115, no. 14, pp. 6382–6390, 1993. View at Google Scholar · View at Scopus - M. K. Nazeeruddin, F. De Angelis, S. Fantacci et al., “Combined experimental and DFT-TDDFT computational study of photoelectrochemical cell ruthenium sensitizers,”
*Journal of the American Chemical Society*, vol. 127, no. 48, pp. 16835–16847, 2005. View at Publisher · View at Google Scholar · View at Scopus - Q. J. Yu, Y. H. Wang, Z. Yi et al., “High-efficiency dye-sensitized solar cells: The influence of lithium ions on exciton dissociation, charge recombination, and surface states,”
*ACS Nano*, vol. 4, no. 10, pp. 6032–6038, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. Dick and B. A. Gregg, “Crystal structure of tris(2,2'-bipyridine)iron(II) bis(hexafluorophosphate), (C
_{10}H_{8}N_{2})_{3}Fe(PF_{6})_{2},”*New Crystal Structures*, vol. 213, no. 2, p. 356, 1998. View at Google Scholar · View at Scopus - S. Ferrere and B. A. Gregg, “Photosensitization of TiO
_{2}by [Fe^{II}(2,2'-bipyridine-4,4'-dicarboxylic acid)_{2}(CN)_{2}]: Band selective electron injection from ultra-short-lived excited states,”*Journal of the American Chemical Society*, vol. 120, no. 4, pp. 843–844, 1998. View at Publisher · View at Google Scholar · View at Scopus - L. M. Lawson Daku, A. Vargas, A. Hauser, A. Fouqueau, and M. E. Casida, “Assessment of density functionals for the high-spin/low-spin energy difference in the low-spin iron(II) tris(2,2'-bipyridine) complex,”
*ChemPhysChem*, vol. 6, no. 7, pp. 1393–1410, 2005. View at Publisher · View at Google Scholar · View at Scopus - S. Ferrere, “New photosensitizers based upon [Fe(L)
_{2}(CN)_{2}] and [Fe(L)_{3}] (L = substituted 2,2'-bipyridine): Yields for the photosensitization of TiO_{2}and effects on the band selectivity,”*Chemistry of Materials*, vol. 12, no. 4, pp. 1083–1089, 2000. View at Publisher · View at Google Scholar · View at Scopus - S. Ferrere, “New photosensitizers based upon [Fe
^{II}(L)_{2}(CN)_{2}] and [Fe^{II}L_{3}], where L is substituted 2,2'-bipyridine,”*Inorganica Chimica Acta*, vol. 329, no. 1, pp. 79–92, 2002. View at Publisher · View at Google Scholar · View at Scopus - F. De Angelis, A. Tilocca, and A. Selloni, “Time-dependent DFT study of [Fe(CN)
_{6}]^{4-}sensitization of TiO_{2}nanoparticles,”*Journal of the American Chemical Society*, vol. 126, no. 46, pp. 15024–15025, 2004. View at Publisher · View at Google Scholar · View at Scopus - W. Gawelda, A. Cannizzo, V.-T. Pham, F. Van Mourik, C. Bressler, and M. Chergui, “Ultrafast nonadiabatic dynamics of [Fe
^{II}(bpy)_{3}]^{2+}in solution,”*Journal of the American Chemical Society*, vol. 129, no. 26, pp. 8199–8206, 2007. View at Publisher · View at Google Scholar · View at Scopus - A. Govindasamy, C. Lv, H. Tsuboi et al., “Theoretical investigation of the photophysical properties of black dye sensitizer [(H
_{3}-tctpy)M(NCS)_{3}] (M = Fe, Ru, Os) in dye sensitized solar cells,”*Japanese Journal of Applied Physics*, vol. 46, no. 4 B, pp. 2655–2660, 2007. View at Publisher · View at Google Scholar · View at Scopus - B. D. Alexander, T. J. Dines, and R. W. Longhurst, “DFT calculations of the structures and vibrational spectra of the [Fe(bpy)
_{3}]^{2+}and [Ru(bpy)_{3}]^{2+}complexes,”*Chemical Physics*, vol. 352, no. 1–3, pp. 19–27, 2008. View at Publisher · View at Google Scholar · View at Scopus - P. V. Bernhardt, G. K. Boschloo, F. Bozoglian, A. Hagfeldt, M. Martínez, and B. Sienra, “Tailoring mixed-valence Co
^{III}/Fe^{II}complexes for their potential use as sensitizers in dye sensitized solar cells,”*New Journal of Chemistry*, vol. 32, no. 43, pp. 705–711, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. J. Frisch, G. W. Trucks, H. B. Schlegel et al.,
*Gaussian 09*, Revision A.1, Gaussian, Wallingford, Conn, USA, 2009. - S. Chiodo, N. Russo, and E. Sicilia, “Newly developed basis sets for density functional calculations,”
*Journal of Computational Chemistry*, vol. 26, no. 2, pp. 175–183, 2005. View at Publisher · View at Google Scholar · View at Scopus - D. Andrae, U. Häußermann, M. Dolg, H. Stoll, and H. Preuß, “Energy-adjusted ab initio pseudopotentials for the second and third row transition elements,”
*Theoretica Chimica Acta*, vol. 77, no. 2, pp. 123–141, 1990. View at Publisher · View at Google Scholar - X. Lu, S. Wei, C. -M. Lawrence Wu, W. Guo, and L. Zhao, “Theoretical characterization of ruthenium complexes containing functionalized bithiophene ligands for dye-sensitized solar cells,”
*Journal of Organometallic Chemistry*, vol. 696, no. 8, pp. 1632–1639, 2011. View at Publisher · View at Google Scholar - R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, “An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules,”
*Journal of Chemical Physics*, vol. 109, no. 19, pp. 8218–8224, 1998. View at Publisher · View at Google Scholar · View at Scopus - G. Scalmani, M. J. Frisch, B. Mennucci, J. Tomasi, R. Cammi, and V. Barone, “Geometries and properties of excited states in the gas phase and in solution: theory and application of a time-dependent density functional theory polarizable continuum model,”
*Journal of Chemical Physics*, vol. 124, no. 9, Article ID 094107, pp. 1–15, 2006. View at Publisher · View at Google Scholar · View at Scopus - M. Cossi and V. Barone, “Time-dependent density functional theory for molecules in liquid solutions,”
*Journal of Chemical Physics*, vol. 115, no. 10, pp. 4708–4717, 2001. View at Publisher · View at Google Scholar · View at Scopus - S. Miertš, E. Scrocco, and J. Tomasi, “Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects,”
*Chemical Physics*, vol. 55, no. 1, pp. 117–129, 1981. View at Google Scholar · View at Scopus - M. Cossi, V. Barone, R. Cammi, and J. Tomasi, “Ab initio study of solvated molecules: a new implementation of the polarizable continuum model,”
*Chemical Physics Letters*, vol. 255, no. 4–6, pp. 327–335, 1996. View at Google Scholar · View at Scopus - V. Barone and M. Cossi, “Quantum calculation of molecular energies and energy gradients in solution by a conductor solvent model,”
*Journal of Physical Chemistry A*, vol. 102, no. 11, pp. 1995–2001, 1998. View at Google Scholar · View at Scopus - V. Barone, M. Cossi, and J. Tomasi, “A new definition of cavities for the computation of solvation free energies by the polarizable continuum model,”
*Journal of Chemical Physics*, vol. 107, no. 8, pp. 3210–3221, 1997. View at Google Scholar · View at Scopus - C. Barolo, M. K. Nazeeruddin, S. Fantacci et al., “Synthesis, characterization, and DFT-TDDFT computational study of a ruthenium complex containing a functionalized tetradentate ligand,”
*Inorganic Chemistry*, vol. 44, no. 20, pp. 6852–6864, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Fantacci, F. De Angelis, and A. Selloni, “Absorption spectrum and solvatochromism of the [Ru(4,4′-COOH-2,2′-bpy)
_{2}(NCS)_{2}] molecular dye by time dependent density functional theory,”*Journal of the American Chemical Society*, vol. 125, no. 14, pp. 4381–4387, 2003. View at Publisher · View at Google Scholar · View at Scopus - M. K. Nazeeruddin, P. Péchy, T. Renouard et al., “Engineering of efficient panchromatic sensitizers for nanocrystalline TiO
_{2}-based solar cells,”*Journal of the American Chemical Society*, vol. 123, no. 8, pp. 1613–1624, 2001. View at Publisher · View at Google Scholar · View at Scopus - R. Argazzi, N. Y. Murakami Iha, H. Zabri, F. Odobel, and C. A. Bignozzi, “Design of molecular dyes for application in photoelectrochemical and electrochromic devices based on nanocrystalline metal oxide semiconductors,”
*Coordination Chemistry Reviews*, vol. 248, no. 13-14, pp. 1299–1316, 2004. View at Publisher · View at Google Scholar · View at Scopus