Next Article in Journal
Facile Synthesis of SnO2 Aerogel/Reduced Graphene Oxide Nanocomposites via in Situ Annealing for the Photocatalytic Degradation of Methyl Orange
Next Article in Special Issue
Reactivity of Atomically Functionalized C-Doped Boron Nitride Nanoribbons and Their Interaction with Organosulfur Compounds
Previous Article in Journal
Novel Magnetic Nanostructured Beads for Cadmium(II) Removal
Previous Article in Special Issue
Computational Investigation of Tuning the Electron-Donating Ability in Metal-Free Organic Dyes Featuring an Azobenzene Spacer for Dye-Sensitized Solar Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure and Electronic Properties of TiO2 Nanoclusters and Dye–Nanocluster Systems Appropriate to Model Hybrid Photovoltaic or Photocatalytic Applications

by
Corneliu I. Oprea
and
Mihai A. Gîrțu
*
Department of Physics and Electronics, Ovidius University of Constanța, 900527 Constanța, Romania
*
Author to whom correspondence should be addressed.
Nanomaterials 2019, 9(3), 357; https://doi.org/10.3390/nano9030357
Submission received: 30 January 2019 / Revised: 18 February 2019 / Accepted: 20 February 2019 / Published: 4 March 2019
(This article belongs to the Special Issue Computational Materials Design for Renewable Energy Applications)

Abstract

:
We report the results of a computational study of TiO2 nanoclusters of various sizes as well as of complex systems with various molecules adsorbed onto the clusters to set the ground for the modeling of charge transfer processes in hybrid organic–inorganic photovoltaics or photocatalytic degradation of pollutants. Despite the large number of existing computational studies of TiO2 clusters and in spite of the higher computing power of the typical available hardware, allowing for calculations of larger systems, there are still studies that use cluster sizes that are too small and not appropriate to address particular problems or certain complex systems relevant in photovoltaic or photocatalytic applications. By means of density functional theory (DFT) calculations, we attempt to find acceptable minimal sizes of the TinO2n+2H4 (n = 14, 24, 34, 44, 54) nanoclusters in correlation with the size of the adsorbed molecule and the rigidity of the backbone of the molecule to model systems and interface processes that occur in hybrid photovoltaics and photocatalysis. We illustrate various adsorption cases with a small rigid molecule based on coumarin, a larger rigid oligomethine cyanine dye with indol groups, and the penicillin V antibiotic having a flexible backbone. We find that the use of the n = 14 cluster to describe adsorption leads to significant distortions of both the cluster and the molecule and to unusual tridentate binding configurations not seen for larger clusters. Moreover, the significantly weaker bonding as well as the differences in the density of states and in the optical spectra suggest that the n = 14 cluster is a poor choice for simulating the materials used in the practical applications envisaged here. As the n = 24 cluster has provided mixed results, we argue that cluster sizes larger than or equal to n = 34 are necessary to provide the reliability required by photovoltaic and photocatalytic applications. Furthermore, the tendency to saturate the key quantities of interest when moving from n = 44 to n = 54 suggests that the largest cluster may bring little improvement at a significantly higher computational cost.

Graphical Abstract

1. Introduction

Nanostructured titanium dioxide is well-known for its numerous and diverse applications [1,2]. Examples of such new uses are in photocatalytic degradation of pollutants [2,3,4], in dye-sensitized solar cells, for energy conversion [5,6], in sensor electronics (for detection of gases, chemicals, or biological materials) [7], in coatings with antibacterial, anticorrosion, antifogging, or self-cleaning properties [8], in drug delivery [9], etc. Due to its abundance, nontoxicity, and high stability under a variety of conditions, titanium dioxide is likely to remain the material of choice for the near future [1].
The physical and chemical properties of TiO2 nanocrystals, dictated to a large extent by their electronic structure, are influenced by the size, shape, organization, and surface properties of the nanoparticles [10,11]. Consequently, understanding the structure-property correlations for TiO2 nanostructured materials is crucial for optimizing their use in practical applications; however, it remains a challenging task even after decades of extensive studies. Density functional theory (DFT) [12,13,14] has, over the past two decades, been the method of choice for most computational studies and quite successful at predicting TiO2 properties and explaining experimental data. Despite some inherent difficulties [14] and limitations [15], DFT methods have been able to treat molecular systems of relatively large sizes, of up to several hundreds of atoms, with moderate computational costs [16,17].
Computational studies have discussed the electronic properties of bulk TiO2 polymorphs [18], the structure and reactivity of anatase surfaces, and the modeling of bare and sensitized TiO2 nanoparticles [16,17]. It was found that rutile is thermodynamically the most stable bulk phase, whereas anatase is very common and stable in nanomaterials [17,19], displaying the highest photocatalytic activity [20,21]. As anatase is the most interesting phase of TiO2 for photovoltaic and photocatalytic applications, it will be the focus of our interest in the following.
In photovoltaic and photocatalytic applications, the key role is played by some electronic properties. For instance, the energy level alignment between the conduction or valence band edges of TiO2 with the ground and excited states of the adsorbed molecule or with the redox level of the electrolyte determine whether a process can take place or not. To investigate the interaction of the individual dye with the TiO2 surface it is natural to use the cluster approach [22], whereas a periodic approach is to be preferred for studying interfaces of periodic crystalline materials, such as perovskites, on TiO2.
Although the TiO2 nanoclusters have been extensively studied theoretically, as reported in recent reviews [16,17], there are still some different perspectives regarding the proper cluster size needed to tackle a particular problem or a certain complex system. For instance, Peng et al. [23] used only a simple TiO6H5 complex to model the anchoring of the dye to the oxide, with the belief that the simple complex would be ‘sufficient for the investigation of geometries, excitations, and electronic structures of the as-investigated dyes before and after binding’. However, they also note that, to study other TiO2 to dyes interactions, their model should be used with precautions [23]. Zhang et al. [24] used the same model, perpetuating the claim that it was suitable for the investigation of organic dyes bound to TiO2. Both studies could only consider bidentate chelate anchoring, in contrast with the experimental and computational data (reviewed in [25]), favoring the bidentate bridging binding of the carboxyl group to the oxide.
Chen et al. [26] chose a Ti5O20H22 model system to describe the binding of several ruthenium sensitizers to the (101) anatase surface, invoking the size limitations needed to adopt a larger basis set. Although they admitted that various models, ranging from 5 to 20 Ti atoms, had been used, showing that the contributions to absorption were not equivalent [26], they went on and used the smallest cluster. There was no discussion of the impact of using one single layer of atoms and of a large number of H atoms to saturate the covalent bonds of O atoms.
Wang et al. [27] used a (TiO2)6 cluster adopted in a previous work, which in fact proposed either a complex with only one Ti atom [24] or a cluster with five Ti atoms [26]. The work influenced Mendizabal et al. [28], who used a cluster with six Ti atoms, Ti6O21H18, cut from the anatase crystal and completed with H atoms to end dangling bonds, and Manzoor et al. [29], who made use of a (TiO2)6 cluster, claiming that it had been ‘extensively employed as a model semiconductor nanoparticle in several studies’. None of these three studies discussed the anchoring, although it has a major role in the charge injection from the dye into the semiconducting oxide.
Analyzing the dependence of the electronic properties of the dye–oxide system on the (TiO2)n nanoclusters’ size, with n = 1, 2, 3, 6, 9, 15, and 38, Sánchez-de-Armas et al. [30] concluded that the minimal TiO4H4 model could reproduce the red shift of the optical spectra upon binding to TiO2. However, they recognized that the electronic structure of the system did change for larger systems, with the orbital on the dye responsible for the excitation moving from below the conduction band edge for the minimal model to well into the conduction band for larger clusters. They concluded that the n = 6 model could provide only a semi-quantitative simulation of all the features in the electronic structure of the complex alizarin–TiO2 system [30]. The same authors argued in subsequent papers that the n = 9 cluster was suitable for larger coumarin-based dyes [31], whereas the n = 15 cluster had the minimal size to provide a complete picture in the case of the smaller catechol molecule binding to the substrate [32].
Early in-depth DFT computational studies of TiO2 clusters were performed by Persson et al. [33,34,35], whose strategy for n equal to 16 and 38 was to remove selected atoms from the (101) surfaces to keep the cluster neutral and stoichiometric. They also performed systematic simulations with the cluster size ranging from n = 16 to 68 [36,37] and from n = 32 to 122 [38,39]. Independently, Jug and co-workers [40] had a different approach, in which they saturated all peripheral oxygen atoms of the clusters with hydrogen atoms and all less than fivefold coordinated titanium atoms with OH groups. They studied (TiO2)n(H2O)m clusters, with n = 33–132 and m = 17–48 [40].
Starting from the original concept of Persson and co-workers [35,36], De Angelis et al. used a cluster with n = 38 [41,42,43] and with n = 82 [44,45]. Both clusters were shown to represent a good trade-off between accuracy and computational convenience and nicely reproduced the main electronic characteristics of TiO2 nanoparticles [41,43], the larger one being used to check the accuracy of the results [46,47]. Larger clusters have also been studied computationally. A study of clusters with n between 58 and 449 and different truncated-bipyramidal shapes chose to saturate the under-coordinated surface atoms by dissociated water molecules [48]. Later, De Angelis and co-workers examined clusters of n = 367 and 411 [49].
In the context of photovoltaic applications [50,51], we showed that, when dealing with small molecules, an n = 24 cluster provides a reasonable compromise between accuracy and computational costs [52,53,54]. In contrast, when studying the photocatalytic activity of TiO2 under UV and visible light [55], we examined the adsorption of various common antibiotics onto TiO2 nanoclusters, and found that the penicillin molecule is strongly distorted when binding to the substrate [56].
In summary, despite the large number of computational studies reported so far for TiO2 clusters, there are still studies that use cluster sizes that are not appropriate to address particular problems or certain complex systems. To find acceptable minimal sizes of the TiO2 nanoclusters, one should answer questions such as: Should the size of the cluster be chosen in correlation with the size of the adsorbed molecule? Does the rigidity of the backbone of the molecule play a role in choosing the appropriate cluster? Additional questions regard the way the anchoring of the adsorbed molecule affects the structure and the optical and charge transfer properties of the system.
In an attempt to find answers to these questions, we report here computational studies of five different TinO2n+2H4 clusters, with n = 14, 24, 34, 44, and 54, to better understand size effects on optimized geometries and charge transfers. Our goal is to explore the structure, the anchoring, electronic spectrum, and optical properties, and the charge injection mechanisms from the adsorbed molecule to such titanium oxide clusters. We choose as adsorbed molecules simple and well-known systems, such as the coumarine-based dye, C343 [52], and penicillin V [56], as well as new molecules, such as OMCD1 [54,57].

2. Materials and Methods

The titanium oxide clusters reported here, regardless of size, were initially cut from the experimental anatase structure [58] to model the (101) surface. In order to ensure the charge neutrality in the presence of under-coordinated Ti atoms [48,59,60], four hydrogen atoms were used to solve the dangling bonds of the oxygen atoms bound to the two Ti atoms at the periphery of the Ti14O30H4 cluster, or to the three-fold coordinated Ti atoms in the corners of Ti24O50H4, Ti34O70H4, Ti44O90H4, and Ti54O110H4 clusters. These structures, for the molecules investigated, C343, OMCD1, and penicillin V, both isolated and adsorbed on the cluster, were energy minimized by DFT [12,13,14] calculations using the hybrid B3LYP exchange-correlation functional [61,62,63] and the split valence 3-21G(d) basis set with polarization functions [64,65].
For the density of electronic states, single-point calculations were performed using the hybrid B3LYP functional with the effective core potentials (ECPs) for Ti, P, S, and I atoms and double-ζ quality basis functions for all atoms via LANL2DZ [66] and accounting for aqueous solvent effects via the conductor-like polarizable continuum model (C-PCM) [67]. The cavity used in the C-PCM calculation was built from spheres centered on heavy nuclei based on the United Atom for Hartree–Fock procedure [68]. The lowest 50 singlet-to-singlet electronic transitions were calculated by Time-Dependent DFT (TD-DFT) [69]. Calculations were performed with the GAUSSIAN09 quantum chemistry package [70], whereas the projection of the density of states and the UV–Vis spectra on different system components were obtained with GaussSum [71].

3. Results

The results of our computational study are reported in three subsections. The first is dedicated to the clusters alone, and the following two to the clusters with adsorbed rigid and flexible molecules, respectively. In each subsection, we start from the optimized geometry and then discuss the electronic spectrum.

3.1. Structure and Properties of TiO2 Nanoclusters

We performed geometry optimization of model clusters with a slight deviation from the TiO2 stoichiometry to avoid the problem of the surface states in the gap. We introduced H atoms or –OH groups to terminate the dangling bonds at the periphery, resulting in the Ti14O30H4, Ti24O50H4, Ti34O70H4, Ti44O90H4, and Ti54O110H4 model clusters (TinO2n+2H4). This approach resulted in compact structures with 4-, 5-, and 6-fold coordinated Ti ions, together with 2- and 3-fold coordinated oxygen atoms [48,52]. Following the geometry optimization, the structure is slightly distorted (see Figure 1) to minimize the surface stresses.
One way to examine the distortion of the lattice is by considering the average value of the Ti–O distance and its corresponding standard deviation, for the various clusters, and to compare them to the experimental values for the bulk [72], as shown in Table 1. The average distance is consistently smaller than the experimental value of 1.950 Å, valid for the bulk oxide, which has two different Ti–O distances for axial and equatorial bonds. Important to note is that the average Ti–O distance increases with the cluster size, the distortion being stronger for smaller clusters, as expected. The lattice distortion is also revealed by the distribution of the Ti–O distances, which widens significantly compared to the bulk. Moreover, the range for values of the Ti–O distances grows with cluster size. For smaller clusters, the surface-to-bulk ratio is larger, as shown below, and the distribution of distances narrower. As the cluster size increases, so does the difference between the surface and internal atoms, as the bonds are more distorted at the surface and less inside the cluster.
The values presented in Table 1 are consistent with the results of DFT calculations for stoichiometric (TiO2)n clusters, with n between 16 and 68, by other groups [36]. We also note that the values reported earlier [52,53] for n = 24, 34, 44 have slight deviations from those presented here due to the different DFT functionals and basis sets used.
Another way to examine the distortion of the lattice is by looking at the length and width of the various clusters and comparing those values to the corresponding distances measured in the bulk, as displayed in Table 1. The length and width of the various clusters were determined between peripheral Ti as well as O atoms.
The width of the clusters was measured as the amount the distance between the peripheral Ti atoms fluctuates from cluster to cluster with respect to the 12.04 Å value in the bulk. However, it can be observed that, despite the fluctuations, the width approaches the bulk value as the cluster size increases, as one would expect. When measuring the width between peripheral O atoms, the values are consistently larger than the one for the bulk, the range of width values being narrower.
In contrast to the width, the length of the clusters, measured between peripheral Ti atoms, is consistently smaller than the corresponding values for the bulk. In the case of the cluster lengths, the tendency is, again, to approach the values for the bulk when the cluster size increases. The lengths measured between peripheral O atoms also follow this tendency.
Examining in parallel the width and the length of the clusters, it can be noticed that the two quantities have opposing trends, particularly when measuring the distances between peripheral O atoms. As the cluster size increases, the width approaches the bulk value from above, whereas the length tends to the bulk from below. This tendency is illustrated in the plot of the relative variation of the cluster width and length with respect to the bulk, displayed in Figure 2.
The dependence on n of the relative cluster lengths shows a tendency to converge towards the bulk value, as expected. The relative variations of the cluster length (which takes values over a larger range when measuring the distance between peripheral O atoms) varies in opposition with the cluster width. The distortion allows for larger distances along one direction and shorter distances along the perpendicular direction. Overall, the distortion of the clusters is extremely large for n = 14, below ±10% for n ≥24, and less than ±5% for n = 54.
Earlier in this section, we stated that the distribution of Ti–O distances widens for larger n, when the surface-to-bulk ratio decreases. To support this statement, we present in Table 2 the numbers of Ti and O atoms of various coordination as well as the surface-to-bulk ratio, determined as the fraction of under-coordinated atoms [36]. The number of bulk (6-fold coordinated) Ti atoms increases slowly with n, as the clusters studied here are intended for modeling interface phenomena in hybrid photovoltaic or photocatalytic applications. For such applications, the most important anatase surface is (101) [21], already shown in Figure 1. Consequently, the clusters modeled here are extended along one direction of that particular surface and not inside, towards the bulk.
The surface-to-bulk ratio is also displayed in Figure 3. It can be seen that the slope of the curve decreases with n and there is a slight tendency of saturation for the particular type of surface cluster of interest here. The small clusters have a large number of under-coordinated atoms, which are likely to affect the electronic and optical properties of the system. Although still far from the bulk, the clusters with n ≥44, displaying the levelling-off of the surface-to-bulk ratio, are likely to be better candidates for describing interface phenomena.
The next step from the analysis of the optimized geometry is to compare the various clusters in terms of energy. For that reason, we present in Figure 4 the binding energy per atom, determined as the difference between the energy of the TinO2n+2H4 cluster and the sum of the calculated energies of the constituents (i.e., of the n Ti atoms, 2n+2 O atoms, and four H atoms), all divided by the total number of atoms (3n + 6).
Regardless of the basis sets used and whether the calculation was performed in vacuum or in water solvent, the binding energy per atom increases with n. This result is expected, because the number of fully coordinated atoms increases with the cluster size and, therefore, the binding becomes stronger for the larger clusters. Important to notice is also the fact that the value of the binding energy per atom is levelling-off at high n, with less than 0.04 eV differences for n ≥ 44.
For each cluster, we report in Table 3 the edge of the valence band, considered at the energy of the highest occupied molecular orbital (HOMO), and the edge of the conduction band, at the lowest unoccupied molecular orbital (LUMO). The geometry relaxation leads to band gaps of 4.38–5.15 eV, values that are sensitive to the DFT functional and basis set used. For that reason, the band gaps presented here have slight deviations from those reported earlier for n = 24, 34, 44 [52,53].
Overall, the values obtained for the band gap are larger than the experimental value of ~3.2 eV for anatase titania [50,73]. The TiO2 band gap has been overestimated by many calculations of model clusters, particularly when using hybrid DFT functionals [17,74]. For instance, for the (TiO2)38 cluster, values of 3.5 eV were obtained from configuration interaction computations using about 3000 configurations [33], whereas a band gap of 3.78 eV was determined using DFT methods [41]. It has to be noticed, however, that the configuration interaction (CI) calculation was for a structure cut from the bulk, with no geometry optimization [33], and that the DFT calculation was performed on a structure optimized by Car–Parrinello molecular dynamics [41], and not by DFT, as in our calculation. Also, the two missing O atoms in the n = 38 model are likely to affect the valence band states, as we shall demonstrate later in this section, possibly leading to a smaller calculated gap.
The values presented in Table 3 are, overall, consistent with the results of DFT calculations for stoichiometric (TiO2)n clusters, with n between 16 and 60, for which the band gap was determined to be between 4.55 and 4.94 eV, with an anomalous value of 4.81 eV for n = 68 [36]. We note that the values reported by us earlier [52,53] for n = 24, 34, 44 have slight deviations from those presented here, due to the different DFT functionals and basis sets used.
Our calculations show that the band gap is smallest, closest to the experimental value, for the largest cluster. The gap does not vary monotonically, likely because of the particular symmetry of the n = 24 cluster. The density of states displayed in Figure 5 reflects the decrease in the band gap with increasing cluster size and reveals the anomaly for n = 24. It can be seen that the peak located at about −8 eV, which consists of the top two states in the valence band (of −7.97 and −8.02 eV, to be precise), is relatively distanced from the adjacent lower energy states, situated at −8.25 eV. This anomaly comes from the stronger electron localization on the two central oxygen atoms, which raises the energy of the orbitals by about 0.25 eV away from the next lower states. For the other clusters, the extra localization does not take place and the four highest states in the valence band are all grouped within a narrow 0.035 eV energy interval.
To better understand the nature of the states in the valence and conduction band we examine more closely the density of the state for the cluster with n = 54, as displayed in Figure 6. The valence band is dominated by the contribution of the O atoms. In contrast, the key contribution in the conduction band comes from the Ti atoms, the participation of the oxygen atoms being small. Worth noting is that the introduction of the hydrogen atoms, which end four dangling bonds, has a minor contribution to the overall density of states but plays an important role in removing from the gap the states localized at the surface.
The density of states curves indicate that the relative weight of the O and Ti atoms in the key molecular orbitals is consistent with the contribution of those atoms for the entire corresponding bands. Specifically, the peak at the valence band edge (the HOMO located just below −8 eV) shows a strong O contribution, whereas the peak at the conduction band edge (the LUMO situated above −4 eV) has a significant Ti weight.
To visually illustrate the nature of the key orbitals discussed above, we display in Figure 7 the electronic density of the highest occupied molecular orbitals and the lowest unoccupied molecular orbitals for three of the five clusters. The two orbitals are especially important as they describe the states at the edges of the valence and conduction bands, respectively.
Upon examination of the HOMOs, we observed the two lobes of the p-type orbital, localized especially on O atoms, and colored red in Figure 7. Moreover, comparing the HOMOs of the three clusters, we note the high electron density strongly localized on the two central oxygen atoms, explaining the higher energy anomaly mentioned above. For the other clusters, the electron density is delocalized, with a significant contribution towards the periphery of the structure.
In contrast, the LUMOs show the four lobes characteristic for d-type orbitals of Ti atoms (colored grey in Figure 7). For these states, the electron density is more substantial in the central part of the structure.
At the end of this section it has to be noted that, although the HOMO–LUMO gap mimics the band gap of the semiconducting oxide, a better approximation for the gap energy, particularly when analyzing the optical properties of the system, is the electronic transition energy [17,46]. The former is determined by DFT calculations at the end of the geometry optimization, whereas the latter is computed by means of TD-DFT methods.
The optical gap is reported in Table 3. The values obtained for the optical gap range from 3.66 to 4.23 eV, systematically smaller than those for the HOMO–LUMO gap and closer to the experimental value of ∼3.4 eV, obtained with optical absorption measurements at low temperatures [75]. Although, overall, the optical gap decreases with increasing cluster size, the variation is not monotonic, the n = 24 anomaly being present. The values presented in Table 3 are, overall, consistent with the results of DFT calculations for stoichiometric (TiO2)n clusters, with n between 16 and 46, for which the TD-DFT gap was between 3.97 and 4.09 eV [36].
Adding to the HOMO energy, obtained by DFT, the energy of the optical transition, determined by TD-DFT, we can estimate the conduction band oxidation potentials (CBOP) (see Table 3). This energy is important for applications, to analyze the charge transfer of the adsorbed molecule to the TiO2 nanocluster [50,51]. The transfer can occur only if there is a proper energy level alignment, such that the excited state of the molecule lies higher than the conduction band’s oxidation potential [17,46].

3.2. Structure and Properties of TiO2 Nanoclusters with Rigid Adsorbed Molecules

In this section, we report results of calculations involving models of TiO2 nanoclusters and two adsorbed molecules of various sizes with a rigid backbone. In each case, we start from the optimized geometry, analyzing the binding to the cluster, and then discuss the electronic spectrum. The anchoring modes of the adsorbed molecules to the TiO2 surface are of crucial importance, the bonding type playing an important role in the electron injection [50,51].

3.2.1. Structure and Properties of TiO2 Nanoclusters with a Small Rigid Adsorbed Molecule

To better understand the role of the cluster size in studying the binding of a small rigid molecule to the TiO2 nanocluster, we chose first the two smallest clusters (n = 14 and 24) and a small molecule with a rigid backbone, a coumarin-based molecule: 11-oxo-2,3,6,7-tetrahydro-1H,5H,11H-pyrano[2,3-f]pyrido[3,2,1-ij]quinoline-10-carboxylic acid (C343). The optimized geometry of C343 has been reported previously by us and other authors [52,76] and, therefore, here we only state that the structures are in agreement with the ones already presented.
Earlier studies have shown that for the organic dyes bearing a carboxylic acid as the anchoring group, the preferred adsorption mode is bidentate bridging, with one proton transferred to a nearby surface oxygen [25,47,77]. Here we analyze the bidentate bridge anchoring of C343 to the model clusters, with two different binding sites in the case of n = 14 and one binding configuration for n = 24. The optimized geometries are displayed in Figure 8 and some of the most important structural parameters are presented in Table 4.
The two anchoring configurations chosen for the n = 14 cluster connect the oxygen atoms of the adsorbed molecule to two neighboring Ti atoms on the same (101) surface of anatase titania. When the geometry optimization starts with the plane of the C343 molecule roughly transverse to the cluster’s width, the final structure exhibits the typical bidentate bridge binding, but with the adsorbed molecule strongly bent towards the cluster (left column of Figure 8). If the optimization is initiated with the C343 molecule oriented roughly along the cluster’s width, the resulting structure is highly unusual, with three Ti–O bonds (middle column of Figure 8).
Comparing geometry parameters for the isolated cluster and the two cluster-molecule systems (Table 4), we note that both the average Ti–O distances and the standard deviation of those distances increase significantly, particularly for the three-bonds case. The n = 14 cluster is significantly distorted in both cases, as also shown by the cluster sizes. For the two-bonds system, both the width and the length of the cluster decrease, whereas for the three-bonds model the width decreases but the length increases.
Table 4 also presents the Ti–O distances at the binding sites. The values are significantly larger than the average distance of 1.950 Å in the bulk anatase titania, suggesting that the bonds are weak, particularly for the smallest cluster. The n = 14 system with two bonds has highly disproportionate Ti–O distances (differences of 0.15 Å), whereas the three-bonds structure has more even distances (differences to within 0.06 Å).
On the other hand, the two-bonds structure has a smaller torsion angle of the carboxyl group relative to the coumarin plane, compared to the three-bonds system, a fact easily seen in Figure 8. This result suggests that the price paid for more even binding distances is a twisted anchor. Further, examining the dihedral angle between the molecule and the Ti atoms in the cluster, we note that the two-bonds system has the C343 dye highly tilted, almost parallel with the cluster. Of the three models analyzed, the molecule is most upright for the n = 24 cluster.
In contrast to the smallest cluster, the n = 24 system is also less distorted, a fact revealed by all key geometry parameters, starting with the average Ti–O distance and its standard deviation and continuing with the width and the length of the cluster. Moreover, the Ti–O distances at the anchoring sites differ by less than 0.03 Å, suggesting a significantly smaller distortion compared to the n = 14 cluster. Also, the torsion angle of the carboxyl group relative to the coumarin plane is relatively small.
To evaluate the role of the basis set on the structural parameters resulting from DFT calculations for the C343-Ti24O50H4 system we also used, aside from 3-21G(d), larger basis sets, LANL2DZ and a combination (DZ&6-31G) consisting of LANL2DZ on Ti and 6-31G(d,p) [78] on O, C, N, and H atoms. Table 4 shows that the average Ti–O distance increases with larger basis sets, getting closer to the value for the bulk. The difference in the length of the Ti–O bond is less than 0.22%, leading to increases in the width and the length of the cluster of about 2.8%.
Puzzled by the unusual bending, by the triple binding, and by the high distortion of the anchoring distances, we explored the energetics of the three systems. Table 5 presents the adsorption energy, determined as the difference between the calculated energy of the constituents (i.e., the oxide cluster and the C343 molecule) and the energy of the system. Binding in the two-bonds C343-Ti14O30H4 system is consistently the weakest, regardless of the basis set used or whether the calculation was in vacuum or in water solvent. Moreover, in one of the calculations the value is positive, indicating that the system would be better off separated. The three-bonds C343-Ti14O30H4 structure is also relatively weakly bound compared to C343-Ti24O50H4, which has the most stable structure of the three.
Next, we examine the density of states for the three model systems considered in this section and compare them to the corresponding curves for the bare TinO2n+2H4 (n = 14, 24) clusters (Figure 9).
The contribution of the C343 molecule to the density of states is more significant in the valence band and particularly in the gap, where the HOMO of the complex system is located. The information provided by the density of state curves in Figure 9 is completed by the electron densities displayed in Figure 10 for the key molecular orbitals. For instance, for all three systems, the HOMO level clearly has a dominant contribution from the adsorbed molecule, the orbital having a π nature due to the conjugated character of the organic dye. In contrast, the LUMOs have the most significant contribution from the oxide, the nature of the molecular orbital being determined by the d atomic orbitals of Ti. However, neither the HOMOs nor the LUMOs are pure. A small weight from the cluster in the HOMOs and from the molecule in the LUMOs (contributions seen both in the density of state curves and in the electron density plots) make possible optical transitions in the visible region, as we shall see in the electronic spectra.
To examine the influence of the molecule adsorbed on the cluster on the density of states of the cluster, we compared the curves with the ones for the bare corresponding cluster (dotted lines in Figure 9). It can be seen that significant differences occur, particularly for the smallest cluster. For instance, the large peaks in the valence band at about −9.32 and −8.55 eV for the bare n = 14 cluster are much less intense for both types of binding and the higher one moves even higher, to −8.46 eV, for the system with three bonds. Similarly, in the conduction band, the peaks at −2.69 and −2.30 eV for the bare cluster are much less intense in the complex system. Moreover, new peaks occur at −2.96 eV. In contrast, the influence of the adsorbed molecule seems to be less significant for the n = 24 cluster.
By enlarging the basis set in the density of states calculations for the C343–Ti24O50H4 system, we found some systematic changes towards higher energies (see Figure S2 of the Supplementary Materials). The valence band edge moves up by up to 0.28 eV, the HOMO by 0.22 eV, and the conduction band edge by 0.19 eV, leading to a decrease in the band gap of 0.09 eV. We note that by using the LANL2DZ basis set on structures optimized with 3-21G(d) and LANL2DZ functions the effect on the energy of the dye states is minor.
Moving on to the electron density of the key orbitals, particularly of the ones in the conduction band, we note that despite the delocalization of the electron over the Ti atoms of the cluster, some charge is present on the atoms of the anchoring group. This fact can be seen in the LUMOs but is much more obvious for the LUMO+5 orbital of the two-bonds and the LUMO+3 of the three-bonds n = 14 complex system, and for the LUMO+2 of the n = 24 model structure (see Figure S3 of the Supplementary Materials). The presence of charge on the anchor plays an important role in charge transfer properties and has to be well-described by the structure used to model the experimental systems.
To understand the effect of the cluster size and the binding configuration on the optical properties of the system, the electronic spectra of the three structures investigated in this section are displayed in Figure 11. The spectra of the complex system of the C343 molecule with the n = 14 cluster exhibits some differences between the two-bonds and the three-bonds configurations. Most peaks occur together, at 575/565 nm, 508/514 nm, or 411/408 nm, for the two-bonds/three-bonds configurations of the n = 14 complex system. There are also two exceptions. The lowest energy transition is relatively low and occurs at 630 nm (~1.97 eV), only for the system with the two-bonds configuration, whereas a high energy transition, of higher intensity, is located at 369 nm and belongs only to the three-bonds system (involving the presence of charge density on the atoms participating in the extra bond).
When comparing the n = 14 with the n = 24 spectra, we notice the expected similarities in terms of peak positions, particularly between the systems with analogous bidentate bridge anchoring. However, the intensity is not as well-described, particularly for the 569 nm and 465 nm transitions of the n = 24 system, which are HOMO→LUMO+2 and HOMO→LUMO+13, respectively (see Table S1 of the Supplementary Materials).
By enlarging the basis set in the TD-DFT calculations for the C343-Ti24O50H4 system, we found a red shift by up to 12 nm of the high wavelength band. This tendency is consistent with the DFT density of states calculations, where we observed unequal shifts up for the HOMO and the conduction band edge, leading to a smaller energy difference between the two.
At this point, it is worth strengthening a remark made earlier, in the discussion of the density of states, regarding the low intensity of the low-energy transitions. These transitions are made possible by the fact that the HOMOs and the LUMOs are not pure but mix to different degrees small weights from the cluster in the HOMOs and from the molecule in the LUMOs. This mixture is important for the photovoltaic and photocatalytic applications.

3.2.2. Structure and Properties of TiO2 Nanoclusters with a Large Rigid Adsorbed Molecule

Having discussed comparatively the anchoring of a small rigid molecule on two clusters of different sizes, we now move on to a larger molecule, but still rigid, to explore the influence of the size of the adsorbed molecule on the entire system. Again, we start from the optimized geometry, analysing the binding to the cluster, and then discuss the electronic spectrum and optical properties.
We report here results obtained for a new dye, an oligomethine cyanine molecule: 5-carboxy-2-(3-(7-((4-(diphenylphospho) phenyl)ethynyl)-1,1,3-trimethyl-1H-benzo[e] indol-2(3H)-ylidene)prop-1-en-1-yl)-1,3,3-trimethyl-3H-indol-1-ium iodide, OMCD1 [54,57] (see Figure 12). We note that the OMCD1 dye has to be a cationic species in order to ensure the π conjugation over the nitrogen atoms in the cyanine constituents. For the calculation, to ensure the overall neutrality, we placed a negative iodine ion in the vicinity of the cyanine parts, at distances of 5 Å and 5.6 Å to the N atoms. This modeling and computational artifice is particularly appropriate when modeling systems with applications in dye-sensitized solar cells, where iodine is already present in the electrolyte.
Comparing geometry parameters for the bare cluster and the cluster-molecule system (Table 6 and Figure 12), we note that the average Ti–O distances and the standard deviation of those distances are both very close, indicating only a small distortion of the cluster due to the adsorbed molecule. Similarly, the width and the length of the cluster both have a small variation from the isolated cluster, confirming the minor distortion caused by the molecule.
Table 6 also presents the Ti–O distances at the binding sites. The values are larger than the average distance of 1.950 Å in the bulk anatase titania. The two bonding Ti–O distances are relatively close, at about 0.02 Å. The torsion angle of the carboxyl group relative to the indoline ring plane is relatively small and the dihedral angle between the molecule and the Ti atoms on the surface of the cluster is close to a right angle, indicating that the molecule is upright. All parameters indicate that the adsorbed molecule, despite its size, has only a minor effect on the cluster, the distortion being small.
The density of states for the model system consisting of the Ti24O50H4 clusters with the adsorbed OMCD1 molecule as well as the corresponding curve for the bare Ti24O50H4 cluster is displayed in Figure 13. As in the case of the smaller adsorbed molecule, presented in the previous section, the contribution of the OMCD1 molecule to the density of states is more significant in the valence band and particularly in the gap, where the HOMO of the complex system is located. The HOMO level clearly has a dominant contribution from the adsorbed molecule, the orbital having a π nature. In contrast, the LUMO has most significant contribution from the oxide, the nature of the molecular orbital being determined by the d atomic orbitals of Ti. We note that in this case too, the HOMO and the LUMO are not pure, a small contribution being present from the cluster in the HOMO and from the molecule in the LUMO.
We noticed in the density of states of the complex system a peak at about −5.86 eV (HOMO-3 in Figure 13), which is due to three occupied 5p orbitals of the negative iodine ion included in the modeling of the OMVD1 molecule. They are localized just below the HOMO and are not contributing to the regeneration process of the dye in the photovoltaic applications.
To examine the influence of the molecule adsorbed on the cluster on the density of states of the cluster, we included in the plot the curve for the bare cluster. It can be seen that some differences occur, for instance, at about −8.95 and −8.59 eV, in the valence band and at −3.20 eV in the conduction band. We note that the states at the edges of the valence and conduction bands are not influenced significantly by the presence of the adsorbed molecule.
The electron density of some key orbitals is displayed in Figure 14. Electron density plots are shown for HOMO-3 (at −5.86 eV), a molecular orbital with a strong contribution from the 5p atomic orbitals of I-, and HOMO (at −5.54 eV), corresponding to the ground state of OMCD1. From the conduction band, we represent the LUMO (at −3.50 eV), the band edge of the oxide, and two molecular orbitals with some charge density delocalized also on the adsorbed molecule, LUMO+2 (at −3.34 eV) and LUMO+11 (at −2.92 eV). These states can be readily identified in the density of states graph, the HOMO, and the excited states of OMCD1 being able to contribute to light absorption followed by a charge transfer in photovoltaic or photocatalytic applications.

3.3. Structure and Properties of TiO2 Nanoclusters with Flexible Adsorbed Molecules

In the previous section, we presented the results of a calculation involving two molecules with a rigid backbone but of different sizes adsorbed on small clusters (n = 14 and 24). In this section, we investigate the case of molecules with a flexible backbone and study the optimized geometry, the binding to the cluster, and the electronic spectrum. The goal is to find the optimal cluster size to model the complex system.
The molecule of choice in this section is penicillin V (phenoxymethylpenicillin, 3,3-Dimethyl-7-oxo-6-(2-phenoxyacetamido)-4-thia-1-azabicyclo[3.2.0]heptane-2-carboxylic acid) (PV), a member of a group of antibiotics still widely used today in human and veterinary medicine, and present in waste water from the pharmaceutical industry and animal farms. The structures of PV, in neutral as well as in deprotonated forms, are presented in Figure 15. It was reported earlier [56] that the deprotonation of the carboxyl group has some influence on the entire structure in the case of penicillin G but less so for penicillin V. Two angles, which are of interest and will be further used, the angle of the PV relative to the carboxyl group and of the rest of the PV with respect to the phenyl group, have values of 28.5° and 1.1° for the neutral molecule and 23.5° and 1.4° for the deprotonated molecule, respectively.
Given the known adsorption tendencies of molecules onto anatase titania [47,77], we performed DFT calculations to determine the optimized geometry of the complex systems starting from the penicillin V molecule with bidentate bridge anchoring through the carboxyl group, with the proton leaving the antibiotic and binding to an O atom of the oxide. Figure 16 and Table 7 present the optimized geometry and the most important structural parameters for the PV molecule adsorbed on all five clusters.
Examining Figure 16, the first and most important observation is that the PV molecule stands out on the smallest cluster, whereas on the other clusters it is strongly distorted, in different ways, depending on the cluster size. The energy of the total system is minimized by bending and twisting of the molecule with a flexible backbone, such that extra bonds can occur and π–π interactions [79] are made possible.
The extra bond occurs, in the case of n = 24, between the central O atom of the PV backbone and a Ti atom on the oxide, the peripheral phenyl group being roughly transverse to the surface. In the case of the n = 34 cluster, the additional bond involves the O atom connected to the phenyl group and the peripheral ring is roughly parallel to the surface. In the case of n = 44 and 54, there are two extra bonds to surface Ti atoms: one involves the O atom next to the phenyl ring (just as for n = 34) and one engaging the S atom. The tendency already observed for n = 34 continues, with the phenyl ring parallel to the surface, at about 2.95 Å, interacting through π–π interactions with the atoms at the surface of the oxide [79].
Comparing the data presented in Table 1 and Table 7, we observed that the presence of the adsorbed penicillin V molecule on the clusters causes some distortions to the structure. The presence of the PV causes the average Ti–O distance to be smaller by 0.65% for the n = 14 cluster. For the other clusters, the average distance is larger, the relative difference decreasing monotonically in absolute value with increasing n, from 0.43% to 0.11%. The standard deviation of the complex systems is significantly different from the values of the isolated cluster only for n = 14, the much wider range of distances indicating significant distortion of the structure. For the cluster width, the relative differences are less than 1% even for n = 14, whereas for the length the relative differences are larger than 2% for n = 14, 24, and 34 and smaller than 0.5% for n = 44 and 54.
The torsion angle between the PV backbone and the Ti–O bond is about 160° for the smallest cluster and around 130° for all others, consistent with the bending of the molecule observed in Figure 16. The angle between the carboxyl branch and the PV backbone is similar to the value for the isolated PV (28.5°) only for the smallest cluster; for the larger ones, the angle is close to 0°. At the other end of the molecule, the PV backbone suffers some mechanical stress, revealed by the torsion angle of the phenyl ring. For the smallest cluster, the angle is close to the value of the isolated PV (1.1°), whereas for the other clusters the differences are significant, going both ways.
The adsorption energy, determined as the difference between the calculated energies of the entire system and of the separated constituents, is displayed in Figure 17. We note that the binding is weak for the smaller clusters and stronger, with a tendency to saturate, for the larger clusters. This result is consistent with the extra bonds that occur when the cluster surface increases, allowing for structural relaxation.
Next, we examined the density of states for the PV molecule adsorbed on TinO2n+2H4 to examine the role of the size of the cluster. An examination of Figure 18 reveals the valence gap located below roughly -8 eV and the conduction band above −4 eV, as seen in Figure 5 for the isolated clusters. What is new is the presence of new states in the gap that were not present in the case of the bare clusters. Both the density of states calculation and the electron density plots for the key orbitals (Figure 19) indicate that the additional states in the gap represent a contribution of the PV molecule.
In the case of the smallest cluster, when the PV molecule is only slightly distorted, the numerous such states that occur are similar to those of the isolated molecule. As the cluster size increases and the energy is lowered by extra bonding and π–π interactions, the additional states in the gap suffer significant changes in their position.
The electron density of the key molecular orbitals of the complex systems, shown in Figure 19, indicate that for all HOMOs the charge is mostly located on the phenyl group, whereas for all LUMOs the charge is distributed on Ti atoms, with a high contribution from d atomic orbitals.
The states in the gap having an obvious dominant PV contribution are in fact mixed, with some small weight from the oxide. This mixture allows, with a small oscillator strength, for transitions that otherwise would have been forbidden. To understand the effect of the cluster size on the optical properties of the system, the electronic absorption spectra of three complex structures are displayed in Figure 20.
The adsorbed PV molecule introduces states in the gap that move the absorption band from the UV towards the visible range of the spectrum. This behaviour is important for both photovoltaic and photocatalytic applications. A careful examination of Figure 20 reveals the presence of weak absorption bands in the blue/violet region of the spectrum, which correspond to HOMO–LUMO transitions. The positions of the optical bands, calculated with TD-DFT methods, vary significantly with the cluster size, consistent with the DFT calculations that provided the density of states.

4. Discussion

We are now in a position to interpret our results and discuss their relevance. We start by noting that our DFT calculations performed for the different TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54) revealed the existence of some distortions with respect to the bulk, which decrease with the increase of n. As the cluster size increases the surface-to-bulk ratio decreases and, as expected, the properties of the model cluster mimic better and better those of the nanoparticles used in typical photovoltaic and photocatalytic applications (whose average size is about 20 nm [50,51], leading to over 15,000 Ti atoms [80]). Our clusters have much smaller sizes, of less than 2 nm, and were intentionally constructed to model interface phenomena, which involve the anatase titania (101) surface to be well-approximated. That is why, when enlarging our clusters, we did not add layers in depth, which would have increased the number of 6-fold coordinated Ti atoms, and decreased the surface-to-bulk ratio. On the contrary, we enlarged the surface to better describe the interface effects, a fact that proved to be crucial in the proper treatment of flexible adsorbed molecules, such as PV.
Our results for the cluster structures compare well, particularly with those obtained by Persson et al. [33,34,35,36,37,38,39], who also used DFT, and to some extent with those of De Angelis et al. [41,42,43,44,45,46,47], who obtained the geometry by molecular dynamics simulations. More recently, Chen and Dixon [81] used a hybrid generic algorithm to determine the equilibrium structure of ultra-small nanoparticles of TiO2, with n from 1 to 384, and for n up to 80 the geometry optimization being finalized with DFT at the B3LYP/DZVP2 level. They report similar tendencies of the TiO2 nanoclusters to display distortions from the bulk, particularly for clusters of smaller sizes. We note that our n = 14 and 24 clusters have a similar structure with their 14a and 24b clusters; for the other sizes, the numbers of Ti atoms are different, similarities being identified between our n = 34 and their 36b as well as our n = 44 and their 48a [82]. Moreover, consistent with our results, their study demonstrated that the band gap of the ultra-small nanoparticles is strongly affected by the particle shape and surface structure, the gap decreasing with cluster size (scaling with n−1/2).
Our calculations indicated, through the tendency to saturate of (i) the surface-to-bulk ratio, (ii) the binding energy per atom, (iii) the DFT calculated HOMO–LUMO band gap as well as (iv) the TD-DFT determined optical gap, that the cluster with n = 54 brings little improvement compared to the one with n = 44, but increases significantly the computational cost. Our findings are consistent with previous reports, particularly those of Persson and co-workers, who used clusters with n = 38 [35] or 46 [37] and, later, 92 [39], and De Angelis and co-workers, who employed clusters with n = 38 [43] and with n = 82 [44].
Jug and co-workers used free and hydrogen-saturated clusters for the simulation of anatase surfaces by means of semi-empirical methods and found that the influence of the cluster size on the convergence of the adsorption energies in the series n = 33, 56, 72, 132 was insignificant [40]. Our DFT calculations, performed on clusters saturated with four peripheral –OH groups, suggest that these groups play a minor role in the density of states of the clusters and, in any case, not for the states that are crucial for photovoltaic or photocatalytic applications. Studying clusters with n = 29 70, and 121, Vorontsov [82,83] used the self-consistent-charge density-functional tight-binding method (SCC-DFTB) and found, on one hand, that the energy of the HOMO and especially the LUMO substantially depends on the location of the surface hydroxyl groups and, on the other hand, that the band gap is underestimated with respect to the experimental value [50].
The study of rigid adsorbed molecules allowed for a double analysis with respect to the size of the cluster as well as the size of the molecule. The anchoring of the C343 molecule on the n = 14 cluster revealed significant distortions of both the cluster and the molecule, with highly either uneven Ti–O bond lengths for the bidentate bridge anchoring, associated with strong bending of the molecule towards the surface or with unusual tridentate binding configurations. Such unusual behavior is not seen for larger clusters, not even for n = 24, as for such structures the adsorbed molecule tends to stand upright, almost perpendicular to the oxide surface. The adsorption energy comparisons showed significantly weaker bonding for n = 14 with respect to n = 24, the increase in size leading to more stable structures. When contrasting the optical spectra, we noticed some similarities in terms of peak positions, particularly between the systems with analogous bidentate bridge anchoring. However, the intensity was not as well-described, particularly for the higher wavelength transitions occurring in the visible range of greater interest for photovoltaic and photocatalytic applications.
As our results contradict those of several previous studies, which used very small clusters, we argue that the TiO6H5 complex used by other authors [23,24] cannot be ‘sufficient for the investigation of geometries, excitations, and electronic structures of the as-investigated dyes before and after binding’ and should be used with precautions. One obvious reason is that such a complex could only consider bidentate chelate anchoring, in contrast with the vast experimental and computational data, indicating bidentate bridging [25]. Moreover, we warn that small clusters, such as Ti5O20H22 [26], (TiO2)6 [27,29,30], Ti6O21H18 [28] (TiO2)9 [31,84], and even (TiO2)15 [32] or (TiO2)16 [85], may not reproduce sufficiently well the electronic structure and optical spectra of the small molecules bound to the substrate.
When attaching the larger but still rigid OMCD1 molecule onto the n = 24 cluster, we observed relatively small distortions of either the cluster or the molecule. The binding configuration was the expected bidentate bridge, with the adsorbed molecule almost perpendicular to the surface. The density of states showed that, although the large molecule adsorbed on the relatively small cluster affected some of the states in the valence and conduction bands, the key orbitals at the edges of the bands were not altered significantly.
The adsorption of a molecule with a flexible backbone on the different clusters showed important differences, which cannot be ignored. The anchoring of the PV molecule on the n = 14 cluster left the molecule almost unaffected. However, that situation proved to be deceiving, as attaching the molecule onto increasingly larger clusters demonstrated that the energy of the system can be lowered substantially by distorting the PV. The bending and twisting of the backbone allowed for energy minimization due to extra bonds and π–π interactions. The tendency to saturate of the adsorption energy indicated that the benefits brought about by the cluster with n = 54 are not significant compared to the one with n = 44. The same tendency to saturate comes from the analysis of the density of states, obtained by DFT calculations, and of the electronic absorption spectra, determined by TD-DFT methods. Our results compare well with the pioneering work of Persson et al. [36,37] and De Angelis et al. [41,44], as well as with more recent reports [86,87,88,89], which use clusters with an n of at least 38.
Now, to link the results of our calculations, summarized above, to practical applications, such as in hybrid photovoltaics or in photocatalytic degradation of pollutants under visible light, we need to address two key processes: light absorption and charge transfer from the molecule to the oxide. For that reason, it is useful to recall the typical approach used to investigate whether a molecule is a good candidate as a dye to sensitize TiO2 in dye-sensitized solar cells (DSSC) [50,51]. One condition that a sensitizer has to fulfil is the proper anchoring to the oxide, which is considered a precondition for a high rate of charge transfer. The –COOH group at the periphery of a dye has the right positions of the atoms and the right symmetry of the molecular orbital to ensure both a strong mechanical binding but also a high overlap of the p orbitals of O and the d orbitals of Ti, for a sizeable electron transfer probability.
A second criterion is the matching of the absorption spectrum of the dye with the solar irradiation spectrum [51]. Our results showed that the optical spectra, provided by TD-DFT calculations for clusters of various sizes, indicate significant changes in both wavelength and intensity of the electronic transitions in the visible range, where the solar intensity reaches its maximum. Moreover, we found that the spectra were generally insufficiently described by the small clusters, sizes of n = 34 and 44 being sufficient for most practical purposes (given the molecules studied here), whereas n = 54 brought little further improvement.
Another condition for the dye is the charge injection into the oxide. The high rate of charge injection for typical organic sensitizers is achieved if the molecule exhibits the “push→pull” effect [51]. This effect is seen in dyes with an electron-rich donor part at one end and an acceptor part at the opposite side, ending in an anchoring group (such as the carboxyl group), and kept together by a π bridge. The key for charge injection is the alternant single/double bond character of the molecular backbone, which allows for electron delocalization. A push-pull effect occurs when the ground state of the dye has the charge localized on the donor part and the excited state, reached after light absorption, has the charge localized on the acceptor part, right next to the TiO2. The push→pull effect can be studied by analyzing how the electron density is localized in the ground and excited states of the dye adsorbed on the cluster. A proper description of the charge delocalization is crucial for this purpose, which, as our results demonstrated, requires a suitable account of the structure, by an adequate choice of the TiO2 cluster.
The charge injection also depends on the proper alignment of the excited state of the dye with respect to the conduction band edge of the oxide. Additionally, another criterion is the dye regeneration, by returning the electron to the dye after passing through and performing work in the external circuit. In DSSCs, this process is performed by the electrolyte (or, equivalently, by a hole conductor) if the redox level of the electrolyte is situated above the ground state of the dye. The calculations leading to the density of states reveal the position of the ground state and the excited state resulting after light absorption, allowing for the analysis of the energy level alignment. Moreover, the use of TD-DFT calculations providing the optical gaps lead to the conduction band oxidation potential (see Figure S5 of the Supplementary materials), which can be used as an alternate to the conduction band edge in the discussion of the alignment.
Our results show that both the conduction band edge and the ground state of the dye are influenced, to various extents, by the size of the cluster. The differences can be observed from the density of states calculations for the C343-TinO2n+2H4 systems, where they are small but significant, but more so for the PV-TinO2n+2H4 structures, where they become important. Under these circumstances, the influence of the cluster size on the proper description of the energy level alignment has to be taken into account when choosing the cluster size used in the simulation.
As a final point, Illas and co-workers [80] showed that the exchange contribution to the hybrid DFT functional strongly affects the calculated bang gap of TiO2. We performed calculations of the C343-Ti24O50H4 system with the B3LYP functional and various basis sets for determining the optimized geometry, the electronic structure, and the UV–Vis spectra. Our calculations are consistent with the early work of Persson and co-workers [36], emphasizing the need for consistency in computations.

5. Conclusions

We reported DFT calculations of several TiO2 clusters and complex systems consisting of the cluster with an adsorbed molecule to study the influence of the cluster size on the electrical and optical properties of materials/interfaces, having in mind hybrid photovoltaic and photocatalytic applications.
We performed calculations for different TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), which revealed the existence of some distortions with respect to the bulk, which are larger for smaller clusters. The tendency of the cluster with n = 54 to saturate the surface-to-bulk ratio, the binding energy per atom, the electronic (DFT) band gap, as well as the optical (TD-DFT) gap brings about little improvement compared to the one with n = 44, but increases significantly the computational cost.
The study of the complex systems, consisting of the cluster with the adsorbed molecule, demonstrated for the n = 14 case significant distortions of both the cluster and the molecule and unusual tridentate binding configurations. Such behaviour is not seen for larger clusters, not even for n = 24. The adsorption energy comparisons showed for n = 14 significantly weaker bonding, differences in the density of states (particularly the positions of the key orbitals and their relative intensity), and in the optical spectra (both in terms of peak wavelength and intensity). Under these circumstances, we conclude that the n = 14 cluster is a poor choice for simulating the materials used in the practical applications envisaged here.
The use of the n = 24 cluster has provided mixed results. First, the inner symmetry of the cluster structure causes stronger electron localization of the electron density at the centre of the cluster, causing the HOMO to move higher in energy and split away from the immediately lower energy states. This anomaly leads to a smaller DFT band gap and TD-DFT optical gap. The cluster showed some distortions but described qualitatively well the binding of the rigid molecules, even for those of a larger size. For flexible molecules, the n = 24 cluster indicated correctly the bending and twisting of the adsorbed molecule but, due to its small size, was not able to describe the correct geometry and the π–π interactions. Consequently, we conclude that the n = 24 cluster can be used well for qualitative calculations, particularly when the adsorbed molecule is rigid, or to obtain some initial information to set up further more elaborate, computational experiments. However, if accuracy is important or the molecule is flexible, the structure is unlikely to provide the precision required by such applications as hybrid photovoltaics or photocatalytic degradation of pollutants under visible light.
When the molecule is adsorbed on the larger clusters, starting from n = 34, the structure and the electronic and optical properties are described in a more reliable way.
We reviewed the important criteria used in computational studies to explain or predict the suitability of various dyes as TiO2 sensitizers in DSSCs or the likelihood of the degradation of organic pollutants under visible light irradiation using anatase titania as a photocatalyst. We emphasized the influence of the cluster size on the proper description of the anchoring to the oxide, ruling out the n = 14 cluster for all systems and using with high caution the n = 24 cluster when molecules are flexible.
The analysis of the matching between the absorption spectrum of the dye with the solar irradiation spectrum relies on quality UV–Vis simulated spectra. Our results showed significant changes in both wavelength and intensity of the electronic transitions in the visible range, the spectra being qualitatively described by the n = 24 cluster, and more accurate results being obtained for larger systems.
The charge injection into the oxide is influenced by the anchoring, by the push→pull effects determined by the charge distribution of the ground and excited states of the dye, and by the energy level alignment between the excited state of the dye and the conduction band edge of the oxide. Hence, the accuracy of the geometry optimization, of the density of states, and of the optical spectrum calculations are very important. Such accuracy is obtained with larger clusters.
The dye regeneration, by passing the electron from the electrolyte back to the dye, can happen only if the redox level of the electrolyte lies higher than the ground state of the dye. A proper approach to check this criterion requires a reliable calculation of the density of states, which, again, is high for larger clusters.
In summary, our results lead us to the conclusion that cluster sizes larger than or equal to n = 34 are necessary to provide the reliability required by photovoltaic and photocatalytic applications. The tendency to saturation of the key quantities of interest when moving from n = 44 to n = 54 suggests that the largest cluster may bring little improvement at a significantly higher computational cost. Probably, for many molecules of practical interest, sizes of n = 34 and particularly n = 44 may be an appropriate choice.
Finally, an important word of caution is mandatory. Consistency in calculations may be an even more important rule in the simulation of practical systems, as the functionals and the basis sets used in the DFT and TD-DFT computations may have an even stronger influence on the final results than the size of the TiO2 cluster.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-4991/9/3/357/s1, Figure S1: a Ti54O110H4 nanocluster illustrating how the widths and the lengths reported in Table 1 were measured, Figure S2: Density of states of the C343-Ti24O50H4 system calculated by DFT using the B3LYP functional and three different basis sets, Figure S3: Electron density of key orbitals C343-TinO2n+2H4 (n = 14, 24), Table S1: Energy, wavelength, oscillator strength and composition of the most intense optical transition for C343-TinO2n+2H4 clusters (n = 14, 24), Figure S4: Structure of PV-Ti14O30H4 to illustrate the geometrical parameters mentioned in Table 7, Figure S5: Energy level alignment for TinO2n+2H4 (n = 14, 24, 34, 44, 54).

Author Contributions

Conceptualization, C.I.O. and M.A.G.; methodology, C.I.O.; validation, M.A.G.; formal analysis, C.I.O.; investigation, C.I.O.; resources, M.A.G.; data curation, C.I.O.; writing (original draft preparation), C.I.O.; writing (review and editing), M.A.G.; visualization, C.I.O.; supervision, M.A.G.; project administration, M.A.G.

Funding

This research was funded, in part, by SNSF and UEFISCDI, under the Romanian-Swiss Research Programme, through the grant RSRP #IZERO-142144/1—PN-II-ID-RSRP-1/2012 and by JINR (Dubna, Russia) –Romania grant programme, project no. 91, topic number 04-4-1121-2015/2020, JINR 322/21.05.2018.

Acknowledgments

The authors acknowledge Professor Reda M. AbdelAal from Suez University, Egypt for the useful discussions regarding the OMCD1 molecule.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Chen, X.; Selloni, A. Introduction: Titanium Dioxide (TiO2) Nanomaterials. Chem. Rev. 2014, 114, 9281–9282. [Google Scholar] [CrossRef] [PubMed]
  2. Fujishima, A.; Honda, K. Electrochemical Photolysis of Water at a Semiconductor Electrode. Nature 1972, 238, 37–38. [Google Scholar] [CrossRef] [PubMed]
  3. Schneider, J.; Matsuoka, M.; Takeuchi, M.; Zhang, J.; Horiuchi, Y.; Anpo, M.; Bahnemann, D.W. Understanding TiO2 Photocatalysis: Mechanisms and Materials. Chem. Rev. 2014, 114, 9919–9986. [Google Scholar] [CrossRef] [PubMed]
  4. Ma, Y.; Wang, X.; Jia, Y.; Chen, X.; Han, H.; Li, C. Titanium Dioxide-Based Nanomaterials for Photocatalytic Fuel Generations. Chem. Rev. 2014, 114, 9987–10043. [Google Scholar] [CrossRef] [PubMed]
  5. Bai, Y.; Mora-Seró, I.; De Angelis, F.; Bisquert, J.; Wang, P. Titanium Dioxide Nanomaterials for Photovoltaic Applications. Chem. Rev. 2014, 114, 10095–10130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kapilashrami, M.; Zhang, Y.; Liu, Y.-S.; Hagfeldt, A.; Guo, J. Probing the Optical Property and Electronic Structure of TiO2 Nanomaterials for Renewable Energy Applications. Chem. Rev. 2014, 114, 9662–9707. [Google Scholar] [CrossRef] [PubMed]
  7. Bai, J.; Zhou, B. Titanium Dioxide Nanomaterials for Sensor Applications. Chem. Rev. 2014, 114, 10131–10176. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, G.; Yang, H.G.; Pan, J.; Yang, Y.Q.; Lu, G.Q.; Cheng, H.-M. Titanium Dioxide Crystals with Tailored Facets. Chem. Rev. 2014, 114, 9559–9612. [Google Scholar] [CrossRef] [PubMed]
  9. Rajh, T.; Dimitrijevic, N.M.; Bissonnette, M.; Koritarov, T.; Konda, V. Titanium Dioxide in the Service of the Biomedical Revolution. Chem. Rev. 2014, 114, 10177–10216. [Google Scholar] [CrossRef] [PubMed]
  10. Cargnello, M.; Gordon, T.R.; Murray, C.B. Solution-Phase Synthesis of Titanium Dioxide Nanoparticles and Nanocrystals. Chem. Rev. 2014, 114, 9319–9345. [Google Scholar] [CrossRef] [PubMed]
  11. Fattakhova-Rohlfing, D.; Zaleska, A.; Bein, T. Three-Dimensional Titanium Dioxide Nanomaterials. Chem. Rev. 2014, 114, 9487–9558. [Google Scholar] [CrossRef] [PubMed]
  12. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef]
  13. Geerlings, P.; De Proft, F.; Langenaeker, W. Conceptual Density Functional Theory. Chem. Rev. 2003, 103, 1793. [Google Scholar] [CrossRef] [PubMed]
  14. Cohen, A.J.; Mori-Sánchez, P.; Yang, W. Challenges for Density Functional Theory. Chem. Rev. 2011, 112, 289. [Google Scholar] [CrossRef] [PubMed]
  15. Onida, G.; Reining, L.; Rubio, A. Electronic excitations: Density-functional versus many-body Green’s-function approaches. Rev. Mod. Phys. 2002, 74, 601. [Google Scholar] [CrossRef]
  16. Catlow, C.R.A.; Bromley, S.T.; Hamad, S.; Mora-Fonz, M.; Sokol, A.A.; Woodley, S.M. Modelling nano-clusters and nucleation. Phys. Chem. Chem. Phys. 2010, 12, 786–811. [Google Scholar] [CrossRef] [PubMed]
  17. De Angelis, F.; Di Valentin, C.; Fantacci, S.; Vittadini, A.; Selloni, A. Theoretical Studies on Anatase and Less Common TiO2 Phases: Bulk, Surfaces, and Nanomaterials. Chem. Rev. 2014, 114, 9708–9753. [Google Scholar] [CrossRef] [PubMed]
  18. Carp, O.; Huisman, C.L.; Reller, A. Photoinduced reactivity of titanium dioxide. Prog. Solid State Chem. 2004, 32, 33–177. [Google Scholar] [CrossRef]
  19. Conesa, J.C. The Relevance of Dispersion Interactions for the Stability of Oxide Phases. J. Phys. Chem. C 2010, 114, 22718–22726. [Google Scholar] [CrossRef]
  20. Kavan, L.; Grätzel, M.; Gilbert, S.E.; Klemenz, C.; Scheel, H.J. Electrochemical and Photoelectrochemical Investigation of Single-Crystal Anatase. J. Am. Chem. Soc. 1996, 118, 6716–6723. [Google Scholar] [CrossRef]
  21. Selloni, A. Anatase shows its reactive side. Nature Mater. 2008, 7, 613–615. [Google Scholar] [CrossRef] [PubMed]
  22. De Angelis, F. Modeling Materials and Processes in Hybrid/Organic Photovoltaics: From Dye-Sensitized to Perovskite Solar Cells. Acc. Chem. Res. 2014, 47, 3349–3360. [Google Scholar] [CrossRef] [PubMed]
  23. Peng, B.; Yang, S.; Li, L.; Cheng, F.; Chen, J. A density functional theory and time-dependent density functional theory investigation on the anchor comparison of triarylamine-based dyes. J. Chem. Phys. 2010, 132, 034305–034309. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, J.; Li, H.B.; Sun, S.L.; Geng, Y.; Wu, Y.; Su, Z.M. Density functional theory characterization and design of high-performance diarylamine-fluorene dyes with different pi spacers for dye-sensitized solar cells. J. Mater. Chem. 2012, 22, 568. [Google Scholar] [CrossRef]
  25. Pastore, M.; De Angelis, F. Computational modelling of TiO2 surfaces sensitized by organic dyes with different anchoring groups: Adsorption modes, electronic structure and implication for electron injection/recombination. Phys. Chem. Chem. Phys. 2012, 14, 920–928. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, J.; Bai, F.-Q.; Wang, J.; Hao, L.; Xie, Z.-F.; Pan, Q.-J.; Zhang, H.-X. Theoretical studies on spectroscopic properties of ruthenium sensitizers absorbed to TiO2 film surface with connection mode for DSSC. Dyes Pigm. 2012, 94, 459–468. [Google Scholar] [CrossRef]
  27. Wang, J.; Li, H.; Ma, N.N.; Yan, L.K.; Su, Z.M. Theoretical studies on organoimido-substituted hexamolybdates dyes for dye-sensitized solar cells (DSSC). Dyes Pigm. 2013, 99, 440. [Google Scholar] [CrossRef]
  28. Mendizabal, F.; Lopéz, A.; Arratia-Pérez, R.; Zapata-Torres, G. Interaction of LD14 and TiO2 in dye-sensitized solar-cells (DSSC): A density functional theory study. Comput. Theor. Chem. 2015, 1070, 117–125. [Google Scholar] [CrossRef]
  29. Manzoor, T.; Asmi, S.; Niaz, S.; Pandith, A.H. Computational studies on optoelectronic and charge transfer properties of some perylene-based donor-p-acceptor systems for dye sensitized solar cell applications. Int. J. Quantum Chem. 2017, 117, e25332. [Google Scholar] [CrossRef]
  30. Sanchez-de-Armas, R.; Oviedo, J.; San Miguel, M.A.; Ordejon, P.; Pruneda, M.; Sanz, J.F. Real-Time TD-DFT Simulations in Dye Sensitized Solar Cells: The Electronic Absorption Spectrum of Alizarin Supported on TiO2 Nanoclusters. J. Chem. Theory Comput. 2010, 6, 2856. [Google Scholar] [CrossRef] [PubMed]
  31. Sanchez-de-Armas, R.; San-Miguel, M.A.; Oviedo, J.; Marquez, A.; Sanz, J.F. Direct vs Indirect Mechanisms for Electron Injection in Dye-Sensitized Solar Cells. Phys. Chem. Chem. Phys. 2012, 14, 225–233. [Google Scholar] [CrossRef] [PubMed]
  32. Sanchez-de-Armas, R.; San-Miguel, M.A.; Oviedo, J.; Marquez, A.; Sanz, J.F. Electronic structure and optical spectra of catechol on TiO2 nanoparticles from real time TD-DFT simulations. Phys. Chem. Chem. Phys. 2011, 13, 1506–1514. [Google Scholar] [CrossRef] [PubMed]
  33. Persson, P.; Bergstrom, R.; Lunell, S. Quantum Chemical Study of Photoinjection Processes in Dye-Sensitized TiO2 Nanoparticles. J. Phys. Chem. B 2000, 104, 10348–10351. [Google Scholar] [CrossRef]
  34. Persson, P.; Gebhardt, J.C.M.; Lunell, S. The Smallest Possible Nanocrystals of Semiionic Oxides. J. Phys. Chem. B 2003, 107, 3336–3339. [Google Scholar] [CrossRef]
  35. Persson, P.; Lundqvist, M.J. Calculated Structural and Electronic Interactions of the Ruthenium Dye N3 with a Titanium Dioxide Nanocrystal. J. Phys. Chem. B 2005, 109, 11918–11924. [Google Scholar] [CrossRef] [PubMed]
  36. Lundqvist, M.J.; Nilsing, M.; Persson, P.; Lunell, S. DFT Study of Bare and Dye-Sensitized TiO2 Clusters and Nanocrystals. Int. J. Quantum Chem. 2006, 106, 3214–3234. [Google Scholar] [CrossRef]
  37. Lundqvist, M.J.; Nilsing, M.; Lunell, S.; Akermark, B.; Persson, P. Spacer and Anchor Effects on the Electronic Coupling in Ruthenium-bis-Terpyridine Dye-Sensitized TiO2 Nanocrystals Studied by DFT. J. Phys. Chem. B 2006, 110, 20513–20525. [Google Scholar] [CrossRef] [PubMed]
  38. Fredin, L.A.; Wärnmark, K.; Sundström, V.; Persson, P. Molecular and Interfacial Calculations of Iron(II) Light Harvesters. ChemSusChem 2016, 9, 667–675. [Google Scholar] [CrossRef] [PubMed]
  39. Gałynska, M.; Persson, P. Emerging Polymorphism in Nanostructured TiO2: Quantum Chemical Comparison of Anatase, Rutile, and Brookite Clusters. Int. J. Quantum Chem. 2013, 113, 2611–2620. [Google Scholar] [CrossRef]
  40. Homann, T.; Bredow, T.; Jug, K. Adsorption of small molecules on the anatase (100) surface. Surf. Sci. 2004, 555, 135–144. [Google Scholar] [CrossRef]
  41. De Angelis, F.; Tilocca, A.; Selloni, A. Time-Dependent DFT Study of [Fe(CN)6]4- Sensitization of TiO2 Nanoparticles. J. Am. Chem. Soc. 2004, 126, 15024–15025. [Google Scholar] [CrossRef] [PubMed]
  42. Dualeh, A.; De Angelis, F.; Fantacci, S.; Moehl, T.; Yi, C.; Kessler, F.; Baranoff, E.; Nazeeruddin, M.K.; Gräzel, M. Influence of Donor Groups of Organic D-π-A Dyes on Open-Circuit Voltage in Solid-State Dye-Sensitized Solar Cells. J. Phys. Chem. C 2012, 116, 1572–1578. [Google Scholar] [CrossRef]
  43. Ronca, E.; Pastore, M.; Belpassi, L.; Tarantelli, F.; De Angelis, F. Influence of the dye molecular structure on the TiO2 conduction band in dye-sensitized solar cells: Disentangling charge transfer and electrostatic effects. Energy Environ. Sci. 2013, 6, 183–193. [Google Scholar] [CrossRef]
  44. De Angelis, F.; Fantacci, S.; Mosconi, E.; Nazeeruddin, M.K.; Grätzel, M. First-Principles Modeling of the Adsorption Geometry and Electronic Structure of Ru(II) Dyes on Extended TiO2 Substrates for Dye-Sensitized Solar Cell Applications. J. Phys. Chem. C 2010, 114, 6054–6061. [Google Scholar] [CrossRef]
  45. De Angelis, F.; Fantacci, S.; Mosconi, E.; Nazeeruddin, M.K.; Grätzel, M. Absorption Spectra and Excited State Energy Levels of the N719 Dye on TiO2 in Dye-Sensitized Solar Cell Models. J. Phys. Chem. C 2011, 115, 8825–8831. [Google Scholar] [CrossRef]
  46. Lasser, L.; Ronca, E.; Pastore, M.; De Angelis, F.; Cornil, J.; Lazzaroni, R.; Beljonne, D. Energy Level Alignment at Titanium Oxide−Dye Interfaces: Implications for Electron Injection and Light Harvesting. J. Phys. Chem. C 2015, 119, 9899–9909. [Google Scholar] [CrossRef]
  47. Pastore, M.; De Angelis, F. Intermolecular Interactions in Dye-Sensitized Solar Cells: A Computational Modeling Perspective. J. Phys. Chem. Lett. 2013, 4, 956–974. [Google Scholar] [CrossRef] [PubMed]
  48. Li, Y.-F.; Liu, Z.-P. Particle Size, Shape and Activity for Photocatalysis on Titania Anatase Nanoparticles in Aqueous Surroundings. J. Am. Chem. Soc. 2011, 133, 15743–15752. [Google Scholar] [CrossRef] [PubMed]
  49. Nunzi, F.; Mosconi, E.; Storchi, L.; Ronca, E.; Selloni, A.; Grätzel, M.; De Angelis, F. Inherent electronic trap states in TiO2 nanocrystals: Effect of saturation and sintering. Energy Environ. Sci. 2013, 6, 1221–1229. [Google Scholar] [CrossRef]
  50. Grätzel, M. Photoelectrochemical cells. Nature 2001, 414, 338–344. [Google Scholar] [CrossRef] [PubMed]
  51. Hagfeldt, A.; Boschloo, G.; Sun, L.; Kloo, L.; Pettersson, H. Dye-sensitized solar cells. Chem. Rev. 2010, 110, 6595–6663. [Google Scholar] [CrossRef] [PubMed]
  52. Oprea, C.I.; Panait, P.; Cimpoesu, F.; Ferbinteanu, M.; Gîrţu, M.A. Density Functional Theory (DFT) Study of Coumarin-based Dyes Adsorbed on TiO2 Nanoclusters—Applications to Dye-Sensitized Solar Cells. Materials 2013, 6, 2372–2392. [Google Scholar] [CrossRef] [PubMed]
  53. Oprea, C.I.; Panait, P.; Lungu, J.; Stamate, D.; Dumbravă, A.; Cimpoesu, F.; Gîrţu, M.A. DFT study of binding and electron transfer from a metal-free dye with carboxyl, hydroxyl and sulfonic anchors to a titanium dioxide nanocluster. Int. J. Photoenergy 2013, 2013, 893850. [Google Scholar] [CrossRef]
  54. Oprea, C.I.; Panait, P.; AbdelAal, R.M.; Gîrțu, M.A. DFT Calculations of Structure and Optical Properties in Wide Band-Gap Semiconductor Clusters for Dye-Sensitized Solar Cells. In Proceedings of the IEEE International Semiconductor Conference (CAS), Sinaia, Romania, 9–11 October 2018; pp. 17–26. [Google Scholar]
  55. Oprea, C.I.; Panait, P.; Gîrțu, M.A. DFT study of binding and electron transfer from colorless aromatic pollutants to a TiO2 nanocluster: Application to photocatalytic degradation under visible light irradiation. Beilstein J. Nanotechnol. 2014, 5, 1016–1030. [Google Scholar]
  56. Oprea, C.I.; Petcu, L.C.; Gîrțu, M.A. DFT Study of Binding and Electron Transfer from Penicillin to a TiO2 Nanocluster: Applications to Photocatalytic Degradation. In Proceedings of the IEEE E-Health and Bioengineering Conference (EHB), Iași, Romania, 19–21 November 2015; pp. 1–4. [Google Scholar]
  57. Panait, P.; Oprea, C.I.; Essam, Z.M.; Lungu, J.; Georgescu, A.; Dumbravă, A.; AbdelAal, R.M.; Gîrțu, M.A. Photoexcitation processes in dye-sensitized solar cells based on oligomethine cyanine dyes—A combined experimental and theoretical study. J. Photochem. Photobiol. A-Chem. 2019. under review. [Google Scholar]
  58. Horn, M.; Schwerdtfeger, C.F.; Meagher, E.P. Refinement of the structure of anatase at several temperatures. Z. Krist. 1972, 136, 273–281. [Google Scholar] [CrossRef]
  59. Qu, Z.-W.; Kroes, G.-J. Theoretical Study of Stable, Defect-Free (TiO2)n Nanoparticles with n = 10–16. J. Phys. Chem. C 2007, 111, 16808–16817. [Google Scholar] [CrossRef]
  60. Iacomino, A.; Cantele, G.; Ninno, D.; Marri, I.; Ossicini, S. Structural, electronic, and surface properties of anatase TiO2 nanocrystals from first principles. Phys. Rev. B 2008, 78, 075405. [Google Scholar] [CrossRef]
  61. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  62. Lee, C.; Yang, W.; Parr, R.G. Development of the colle-salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  63. Becke, A.D. Density-functional thermo-chemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  64. Pietro, W.J.; Francl, M.M.; Hehre, W.J.; Defrees, D.J.; Pople, J.A.; Binkley, J.S. Self-Consistent Molecular Orbital Methods. 24. Supplemented small split-valence basis-sets for 2nd-row elements. J. Am. Chem. Soc. 1982, 104, 5039–5048. [Google Scholar] [CrossRef]
  65. Dobbs, K.D.; Hehre, W.J. Molecular-orbital theory of the properties of inorganic and organometallic compounds. 5. Extended basis-sets for 1st-row transition-metals. J. Comp. Chem. 1987, 8, 861–879. [Google Scholar] [CrossRef]
  66. Hay, P.J.; Wadt, W.R. Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985, 82, 299–311. [Google Scholar] [CrossRef]
  67. Barone, V.; Cossi, M. Quantum calculation of molecular energies and energy gradients in solution by a conductor solvent model. J. Phys. Chem. A 1998, 102, 1995–2001. [Google Scholar] [CrossRef]
  68. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105, 2999–3093. [Google Scholar] [CrossRef] [PubMed]
  69. Casida, M.E.; Jamorski, C.; Casida, K.C.; Salahub, D.R. Molecular excitation energies to high-lying bound states from time-dependent density-functional response theory: Characterization and correlation of the time dependent local density approximation ionization threshold. J. Chem. Phys. 1998, 108, 4439–4449. [Google Scholar] [CrossRef]
  70. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. GAUSSIAN 09, Revision A.02. Available online: http://gaussian.com/g09citation/ (accessed on 22 June 2018).
  71. O’Boyle, N.M.; Tenderholt, A.L.; Langner, K.M. Cclib: A library for package-independent computational chemistry algorithms. J. Comp. Chem. 2008, 29, 839–845. [Google Scholar]
  72. Lazzeri, M.; Vittadini, A.; Selloni, A. Structure and energetics of stoichiometric TiO2 anatase surfaces. Phys. Rev. B 2011, 63, 155409. [Google Scholar] [CrossRef]
  73. Scanlon, D.O.; Dunnill, C.W.; Buckeridge, J.; Shevlin, S.A.; Logsdail, A.J.; Woodley, S.M.; Catlow, C.R.A.; Powell, J.M.; Palgrave, R.G.; Parkin, I.P.; et al. Band alignment of rutile and anatase TiO2. Nat. Mater. 2013, 12, 798–801. [Google Scholar] [CrossRef] [PubMed]
  74. Zhang, Y.F.; Lin, W.; Li, Y.; Ding, K.N.; Li, J.Q. A Theoretical Study on the Electronic Structures of TiO2: Effect of Hartree−Fock Exchange. J. Phys. Chem. B 2005, 109, 19270–19277. [Google Scholar] [CrossRef] [PubMed]
  75. Tang, H.; Lévy, F.; Berger, H.; Schmid, P.E. Urbach tail of anatase TiO2. Phys. Rev. B 1995, 52, 7771–7774. [Google Scholar] [CrossRef]
  76. Hara, K.; Sato, T.; Katoh, R.; Furube, A.; Ohga, Y.; Shinpo, A.; Suga, S.; Sayama, K.; Sugihara, H.; Arakawa, H. Molecular design of coumarin dyes for efficient dye-sensitized solar cells. J. Phys. Chem. B 2003, 107, 597–606. [Google Scholar] [CrossRef]
  77. Mosconi, E.; Selloni, A.; De Angelis, F. Solvent Effects on the Adsorption Geometry and Electronic Structure of Dye-Sensitized TiO2: A First-Principles Investigation. J. Phys. Chem. C 2012, 116, 5932–5940. [Google Scholar] [CrossRef]
  78. Francl, M.M.; Pietro, W.J.; Hehre, W.J.; Binkley, J.S.; DeFrees, D.J.; Pople, J.A.; Gordon, M.S. Self-Consistent Molecular Orbital Methods. 23. A polarization-type basis set for 2nd-row elements. J. Chem. Phys. 1982, 77, 3654–3665. [Google Scholar] [CrossRef]
  79. Schottel, B.L.; Chifotides, H.T.; Dunbar, K.R. Anion-π interactions. Chem. Soc. Rev. 2008, 37, 68–83. [Google Scholar] [CrossRef] [PubMed]
  80. Ko, K.C.; Lee, J.Y.; Illas, F. Modeling realistic titania nanoparticles. In Frontiers of Nanoscience, Computational Modelling of Nanoparticles; Bromley, S.T., Woodley, S.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2018; Volume 12, pp. 205–238. [Google Scholar]
  81. Chen, M.; Dixon, D.A. Modeling the formation of TiO2 ultra-small nanoparticles. Nanoscale 2017, 9, 7143–7162. [Google Scholar] [CrossRef] [PubMed]
  82. Vorontsov, A.V. Effect of the Structure of Small Anatase Nanoparticles on the Localization of Photogenerated Charge Carriers. Kinet. Catal. 2017, 58, 688–694. [Google Scholar] [CrossRef]
  83. Vorontsov, A.V.; Smirniotis, P.G. Size and surface groups effects in decahedral anatase nanoparticles for photocatalytic applications. J. Photochem. Photobiol. Chem. 2018, 363, 51–60. [Google Scholar] [CrossRef]
  84. Xu, Z.; Li, Y.; Zhang, W.; Yuan, S.; Hao, L.; Xu, T.; Lu, X. DFT/TD-DFT study of novel T shaped phenothiazine-based organic dyes for dye-sensitized solar cells applications. Spectrochim. Acta. A Mol. Biomol. Spectrosc. 2019, 212, 272–280. [Google Scholar] [CrossRef] [PubMed]
  85. Lu, T.-F.; Li, W.; Chen, J.; Tang, J.; Bai, F.-Q.; Zhang, H.-X. Promising pyridinium ylide based anchors towards high-efficiency dyes for dye-sensitized solar cells applications: Insights from theoretical investigations. Electrochim. Acta 2018, 283, 1798–1805. [Google Scholar] [CrossRef]
  86. Feng, S.; Sun, Z.-Z.; Li, Q. How to screen a promising anchoring group from heterocyclic components in dye sensitized solar cell:A theoretical investigation. Electrochim. Acta 2019, 296, 545–554. [Google Scholar] [CrossRef]
  87. Zhu, H.-C.; Zhang, J.; Wang, Y.-L. Adsorption orientation effects of porphyrin dyes on the performance of DSSC: Comparison of benzoic acid and tropolone anchoring groups binding onto the TiO2 anatase (101) surface. Appl. Surf. Sci. 2018, 433, 1137–1147. [Google Scholar] [CrossRef]
  88. Zhu, H.-C.; Li, C.-F.; Fu, Z.-H.; Wei, S.-S.; Zhu, X.-F.; Zhang, J. Increasing the open-circuit voltage and adsorption stability of squaraine dye binding onto the TiO2 anatase (1 0 1) surface via heterocyclic anchoring groups used for DSSC. Appl. Surf. Sci. 2018, 455, 1095–1105. [Google Scholar] [CrossRef]
  89. Zhang, J.; Zhu, H.-C.; Zhong, R.-L.; Wang, L.; Su, Z.-M. Promising heterocyclic anchoring groups with superior adsorption stability and improved IPCE for high-efficiency noncarboxyl dye sensitized solar cells: A theoretical study. Org. Electron. 2018, 54, 104–113. [Google Scholar] [CrossRef]
Figure 1. Top and lateral views of the optimized structures of the Ti14O30H4, Ti24O50H4, Ti34O70H4, Ti44O90H4, and Ti54O110H4 nanoclusters (from left to right, respectively) modeling the anatase titania (101) surface. The geometry optimization was performed by density functional theory (DFT) at the B3LYP/3-21G(d) level. Atom colors are grey for Ti, red for O, and light grey for H. In the top view, the widths and the lengths reported in Table 1 are measured horizontally and vertically, respectively, as shown in Figure S1 of the Supplementary Materials.
Figure 1. Top and lateral views of the optimized structures of the Ti14O30H4, Ti24O50H4, Ti34O70H4, Ti44O90H4, and Ti54O110H4 nanoclusters (from left to right, respectively) modeling the anatase titania (101) surface. The geometry optimization was performed by density functional theory (DFT) at the B3LYP/3-21G(d) level. Atom colors are grey for Ti, red for O, and light grey for H. In the top view, the widths and the lengths reported in Table 1 are measured horizontally and vertically, respectively, as shown in Figure S1 of the Supplementary Materials.
Nanomaterials 09 00357 g001
Figure 2. The relative cluster width, ( w n w bulk ) / w bulk , and length, ( l n l bulk ) / l bulk , of TinO2n+2H4, nanoclusters (where n = 14, 24, 34, 44, 54) with respect to the corresponding value for the bulk anatase TiO2 [72]. The distances are measured between peripheral Ti (grey) or O (red) atoms.
Figure 2. The relative cluster width, ( w n w bulk ) / w bulk , and length, ( l n l bulk ) / l bulk , of TinO2n+2H4, nanoclusters (where n = 14, 24, 34, 44, 54) with respect to the corresponding value for the bulk anatase TiO2 [72]. The distances are measured between peripheral Ti (grey) or O (red) atoms.
Nanomaterials 09 00357 g002
Figure 3. The surface-to-bulk ratio (defined as the fraction of under-coordinated atoms) dependence on the size of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), geometry optimized by DFT/B3LYP/3-21G(d) in vacuum.
Figure 3. The surface-to-bulk ratio (defined as the fraction of under-coordinated atoms) dependence on the size of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), geometry optimized by DFT/B3LYP/3-21G(d) in vacuum.
Nanomaterials 09 00357 g003
Figure 4. Dependence of the binding energy per atom, [ E n n E Ti ( 2 n + 2 ) E O 4 E H ] / ( 3 n + 6 ) , on the size of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated by B3LYP functional and 3-21G(d) or LANL2DZ basis sets in vacuum or in water solvent.
Figure 4. Dependence of the binding energy per atom, [ E n n E Ti ( 2 n + 2 ) E O 4 E H ] / ( 3 n + 6 ) , on the size of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated by B3LYP functional and 3-21G(d) or LANL2DZ basis sets in vacuum or in water solvent.
Nanomaterials 09 00357 g004
Figure 5. The density of states (in arbitrary units) of the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated at the DFT/B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with a full width at half maximum (FWHM) of 0.1 eV.
Figure 5. The density of states (in arbitrary units) of the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated at the DFT/B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with a full width at half maximum (FWHM) of 0.1 eV.
Nanomaterials 09 00357 g005
Figure 6. Total density of states (DoS) of the Ti54O110H4 cluster, calculated at the DFT/B3LYP/LANL2DZ level in water, showing, separately, the contribution of Ti atoms, O atoms, and –OH groups. Energy levels were convoluted with Gaussian distributions with an FWHM of 0.1 eV.
Figure 6. Total density of states (DoS) of the Ti54O110H4 cluster, calculated at the DFT/B3LYP/LANL2DZ level in water, showing, separately, the contribution of Ti atoms, O atoms, and –OH groups. Energy levels were convoluted with Gaussian distributions with an FWHM of 0.1 eV.
Nanomaterials 09 00357 g006
Figure 7. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of TinO2n+2H4 clusters (n = 14, 24, 54, from left to right, respectively), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first two rows display the top and a perspective view of the LUMOs, whereas the bottom two rows show the top and a perspective view of the HOMOs. Atom colors: Ti, grey; O, red; and H, light grey.
Figure 7. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of TinO2n+2H4 clusters (n = 14, 24, 54, from left to right, respectively), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first two rows display the top and a perspective view of the LUMOs, whereas the bottom two rows show the top and a perspective view of the HOMOs. Atom colors: Ti, grey; O, red; and H, light grey.
Nanomaterials 09 00357 g007
Figure 8. Top (first row), side (second row), and front (third row) views of the optimized structures of the C343 dye molecule adsorbed on the Ti14O30H4 cluster (left and middle column) and on the Ti24O50H4 cluster (right column). Calculations performed by DFT at the B3LYP/3-21G(d) level in vacuum. The side and front views also show the three neighbouring Ti atoms to which the molecule can bind. Atom colors are: C, grey; O, red; N, blue; and H, light grey.
Figure 8. Top (first row), side (second row), and front (third row) views of the optimized structures of the C343 dye molecule adsorbed on the Ti14O30H4 cluster (left and middle column) and on the Ti24O50H4 cluster (right column). Calculations performed by DFT at the B3LYP/3-21G(d) level in vacuum. The side and front views also show the three neighbouring Ti atoms to which the molecule can bind. Atom colors are: C, grey; O, red; N, blue; and H, light grey.
Nanomaterials 09 00357 g008
Figure 9. The density of states (DoS) of the systems consisting of the C343 molecule adsorbed on the Ti14O30H4 cluster via (a) two or (b) three O–Ti bonds and on the Ti24O50H4 cluster (c). Dotted curves represent the density of states of the corresponding bare TinO2n+2H4 clusters (n = 14, 24). Calculations were performed by DFT at the B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with of 0.1 eV full width at half maximum.
Figure 9. The density of states (DoS) of the systems consisting of the C343 molecule adsorbed on the Ti14O30H4 cluster via (a) two or (b) three O–Ti bonds and on the Ti24O50H4 cluster (c). Dotted curves represent the density of states of the corresponding bare TinO2n+2H4 clusters (n = 14, 24). Calculations were performed by DFT at the B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with of 0.1 eV full width at half maximum.
Nanomaterials 09 00357 g009
Figure 10. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the systems consisting of the C343 molecule adsorbed on the Ti14O30H4 cluster via two (left column) or three (middle column) O–Ti bonds and on the Ti24O50H4 cluster (right column), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first row displays the LUMOs, whereas the bottom row shows the HOMOs. Atom colors: Ti, grey; O, red; and H, light grey.
Figure 10. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the systems consisting of the C343 molecule adsorbed on the Ti14O30H4 cluster via two (left column) or three (middle column) O–Ti bonds and on the Ti24O50H4 cluster (right column), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first row displays the LUMOs, whereas the bottom row shows the HOMOs. Atom colors: Ti, grey; O, red; and H, light grey.
Nanomaterials 09 00357 g010
Figure 11. Simulated UV–Vis absorption spectra of the systems consisting of the C343 dye molecule adsorbed on the TinO2n+2H4 clusters, calculated in water solvent by TD-DFT at the B3LYP/LANL2DZ level for n = 14. For n = 24, the calculations were performed with the LANL2DZ basis set on the geometry optimized with 3-21G(d) or LANL2DZ and a combination of LANL2DZ for Ti atoms and 6-31G(d,p) for on O, C, N, and H atoms (DZ & 6-31G) on the structure determined with the same functions. The spectral lines were convoluted with Gaussian distributions of 1000 cm−1 full width at half maximum.
Figure 11. Simulated UV–Vis absorption spectra of the systems consisting of the C343 dye molecule adsorbed on the TinO2n+2H4 clusters, calculated in water solvent by TD-DFT at the B3LYP/LANL2DZ level for n = 14. For n = 24, the calculations were performed with the LANL2DZ basis set on the geometry optimized with 3-21G(d) or LANL2DZ and a combination of LANL2DZ for Ti atoms and 6-31G(d,p) for on O, C, N, and H atoms (DZ & 6-31G) on the structure determined with the same functions. The spectral lines were convoluted with Gaussian distributions of 1000 cm−1 full width at half maximum.
Nanomaterials 09 00357 g011
Figure 12. Optimized structures of the OMCD1 molecule (a) and the system OMCD1–Ti24O50H4, side view (b) and front view (c), calculated by DFT at the B3LYP/3-21G(d) level. Atom colors are: C, grey; O, red; N, blue; P, orange; I, purple; and H, light grey.
Figure 12. Optimized structures of the OMCD1 molecule (a) and the system OMCD1–Ti24O50H4, side view (b) and front view (c), calculated by DFT at the B3LYP/3-21G(d) level. Atom colors are: C, grey; O, red; N, blue; P, orange; I, purple; and H, light grey.
Nanomaterials 09 00357 g012
Figure 13. Density of states of the OMCD1–Ti24O50H4 system, compared to the bare Ti24O50H4 cluster, calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The contributions of the various atoms are: Ti, blue; O, red; –OH groups, grey; dye molecules, black; iodine anion, green. Energy levels were convoluted with Gaussian distributions with a full width at half maximum of 0.1 eV.
Figure 13. Density of states of the OMCD1–Ti24O50H4 system, compared to the bare Ti24O50H4 cluster, calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The contributions of the various atoms are: Ti, blue; O, red; –OH groups, grey; dye molecules, black; iodine anion, green. Energy levels were convoluted with Gaussian distributions with a full width at half maximum of 0.1 eV.
Nanomaterials 09 00357 g013
Figure 14. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the system consisting of the OMCD1 molecule adsorbed onto the Ti24O50H4 cluster, calculated by DFT at the B3LYP/LANL2DZ level in water solvent. Atom colors: T, light grey; C, grey; O, red; N, blue; P, orange; I, purple; and H, small light grey.
Figure 14. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the system consisting of the OMCD1 molecule adsorbed onto the Ti24O50H4 cluster, calculated by DFT at the B3LYP/LANL2DZ level in water solvent. Atom colors: T, light grey; C, grey; O, red; N, blue; P, orange; I, purple; and H, small light grey.
Nanomaterials 09 00357 g014
Figure 15. The optimized structure of (a) neutral and (b) deprotonated penicillin V calculated by DFT at the B3LYP/3-21G(d) level. Atom colors: C, grey; O, red; N, blue; S, yellow; and H, light grey.
Figure 15. The optimized structure of (a) neutral and (b) deprotonated penicillin V calculated by DFT at the B3LYP/3-21G(d) level. Atom colors: C, grey; O, red; N, blue; S, yellow; and H, light grey.
Nanomaterials 09 00357 g015
Figure 16. Top view (first column), side view (second column), and perspective view (third column) of the optimized structures of the Penicillin V molecule adsorbed on TinO2n+2H4 clusters, with n = 14, 24, 34, 44, 54, on the first, second, third, fourth, and fifth rows, respectively. Calculations performed by DFT at the B3LYP/3-21G(d) level in vacuum. Atom colors are: C, grey; O, red; N, blue; S, yellow; and H, light grey.
Figure 16. Top view (first column), side view (second column), and perspective view (third column) of the optimized structures of the Penicillin V molecule adsorbed on TinO2n+2H4 clusters, with n = 14, 24, 34, 44, 54, on the first, second, third, fourth, and fifth rows, respectively. Calculations performed by DFT at the B3LYP/3-21G(d) level in vacuum. Atom colors are: C, grey; O, red; N, blue; S, yellow; and H, light grey.
Nanomaterials 09 00357 g016
Figure 17. Adsorption energy of PV on the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated by DFT at the B3LYP/3-21G(d) level in vacuum and in water solvent.
Figure 17. Adsorption energy of PV on the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54), calculated by DFT at the B3LYP/3-21G(d) level in vacuum and in water solvent.
Nanomaterials 09 00357 g017
Figure 18. Density of states (DoS) of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54). Calculations were performed by DFT at the the B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with a 0.1 eV full width at half maximum.
Figure 18. Density of states (DoS) of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54). Calculations were performed by DFT at the the B3LYP/LANL2DZ level in water solvent. Energy levels were convoluted with Gaussian distributions with a 0.1 eV full width at half maximum.
Nanomaterials 09 00357 g018
Figure 19. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 cluster (n = 14, left; n = 24, middle; and n = 34, right column), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first row displays the LUMOs, whereas the bottom row shows the HOMOs. Atom colors: Ti, light grey; C, grey; O, red; N, blue; S, Yellow; and H, very light grey.
Figure 19. Isodensity surfaces (0.03 e/bohr3) of the key molecular orbitals of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 cluster (n = 14, left; n = 24, middle; and n = 34, right column), calculated by DFT at the B3LYP/LANL2DZ level in water solvent. The first row displays the LUMOs, whereas the bottom row shows the HOMOs. Atom colors: Ti, light grey; C, grey; O, red; N, blue; S, Yellow; and H, very light grey.
Nanomaterials 09 00357 g019
Figure 20. Simulated UV–Vis absorption spectra of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 clusters (n = 14, 24, 34, 44), calculated by TD-DFT at the B3LYP/LANL2DZ level in water solvent. The inset shows details of a higher wavelengths range. The spectral lines were convoluted with Gaussian distributions of 1000 cm−1 full width at half maximum.
Figure 20. Simulated UV–Vis absorption spectra of the systems consisting of the PV molecule adsorbed on the TinO2n+2H4 clusters (n = 14, 24, 34, 44), calculated by TD-DFT at the B3LYP/LANL2DZ level in water solvent. The inset shows details of a higher wavelengths range. The spectral lines were convoluted with Gaussian distributions of 1000 cm−1 full width at half maximum.
Nanomaterials 09 00357 g020
Table 1. Structural parameters describing the optimized geometry of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54) performed by DFT at the B3LYP/3-21G(d) level in vacuum. The average Ti–O distance and its standard deviation, and the cluster width and length, defined by the positions of the peripheral Ti or O atoms and determined by calculations are compared to the experimental values of the bulk TiO2 [72]. All quantities are expressed in Å.
Table 1. Structural parameters describing the optimized geometry of TinO2n+2H4 clusters (n = 14, 24, 34, 44, 54) performed by DFT at the B3LYP/3-21G(d) level in vacuum. The average Ti–O distance and its standard deviation, and the cluster width and length, defined by the positions of the peripheral Ti or O atoms and determined by calculations are compared to the experimental values of the bulk TiO2 [72]. All quantities are expressed in Å.
ParameterTi14O30H4Ti24O50H4Ti34O70H4Ti44O90H4Ti54O110H4TiO2
r(Ti–O)1.8341.8761.8891.8961.9011.950
σr(Ti–O)0.0740.1000.1000.1030.1050.022
Width (Ti…Ti)12.4512.1811.8711.9412.0012.04
Width (O…O)-13.2813.3613.2813.2912.83
Length (Ti…Ti)3.557.2010.6714.4718.05
Length (O…O)4.948.1612.0515.6318.76
Length bulk (Ti…Ti)3.787.5911.3715.1418.92
Table 2. The number of Ti or O atoms m-fold coordinated (mx) to O or Ti atoms, respectively, and the surface-to-bulk ratio (s/b, determined as the fraction of under-coordinated atoms [36]).
Table 2. The number of Ti or O atoms m-fold coordinated (mx) to O or Ti atoms, respectively, and the surface-to-bulk ratio (s/b, determined as the fraction of under-coordinated atoms [36]).
ParameterTi14O30H4Ti24O50H4Ti34O70H4Ti44O90H4Ti54O110H4TiO2 (Bulk)
Ti, 6x02468all Ti atoms
Ti, 5x6142230380
Ti, 4x888880
O, 3x618304254all O atoms
O, 2x20283644520
O, 1x444440
s/b (%)86.473.067.364.262.20
Table 3. Valence (EVB) and conduction (ECB) band edges, the band gap (ECBEVB), the optical gap, and the conduction band oxidation potentials (CBOP) (in eV), calculated for clusters of various sizes, using DFT and Time-Dependent DFT (TD-DFT), the B3LYP functional, LANL2DZ basis sets, and water solvent.
Table 3. Valence (EVB) and conduction (ECB) band edges, the band gap (ECBEVB), the optical gap, and the conduction band oxidation potentials (CBOP) (in eV), calculated for clusters of various sizes, using DFT and Time-Dependent DFT (TD-DFT), the B3LYP functional, LANL2DZ basis sets, and water solvent.
ClusterEVB
HOMO
ECB
LUMO
ECBEVB
HL Gap
E0-0
Optical Gap
CBOP
EVB + E0-0
Ti14O30H4−8.52−3.375.154.23−4.29
Ti24O50H4−7.97−3.534.443.66−4.31
Ti34O70H4−8.21−3.704.523.81−4.41
Ti44O90H4−8.17−3.774.403.75−4.42
Ti54O110H4−8.17−3.794.383.75−4.42
HOMO, highest occupied molecular orbital; LUMO, lowest unoccupied molecular orbital.
Table 4. The structural parameters of the model system consisting of the TinO2n+2H4 clusters (n = 14, 24) with an adsorbed C343 molecule, obtained after geometry optimization by DFT with the B3LYP functional, in vacuum. The basis sets used were 3-21G(d) for n = 14 and 3-21G(d), LANL2DZ, and a combination of LANL2DZ on Ti and 6-31G(d,p) on O, C, N, and H atoms (DZ&6-31G) for n = 24. Geometrical parameters describing the adsorption modes: bond lengths (in Å), torsion angle of the carboxyl group relative to the coumarin plane, and dihedral angle defined by the coumarin plane and the plane of the three neighbour Ti atoms (both in deg).
Table 4. The structural parameters of the model system consisting of the TinO2n+2H4 clusters (n = 14, 24) with an adsorbed C343 molecule, obtained after geometry optimization by DFT with the B3LYP functional, in vacuum. The basis sets used were 3-21G(d) for n = 14 and 3-21G(d), LANL2DZ, and a combination of LANL2DZ on Ti and 6-31G(d,p) on O, C, N, and H atoms (DZ&6-31G) for n = 24. Geometrical parameters describing the adsorption modes: bond lengths (in Å), torsion angle of the carboxyl group relative to the coumarin plane, and dihedral angle defined by the coumarin plane and the plane of the three neighbour Ti atoms (both in deg).
ParameterC343 on Ti14O30H4 (Two Bonds)C343 on Ti14O30H4 (Three Bonds) C343 on Ti24O50H4
3-21G(d)LANL2DZDZ&6-31G
r(Ti–O)1.8521.8801.8871.8911.891
σr(Ti–O)0.0900.1190.1080.1210.116
Width (Ti…Ti)12.4112.0711.9812.4412.23
Width (O…O)--13.1313.4113.49
Length (Ti…Ti)3.483.757.257.267.30
Length (O…O)5.095.648.288.368.35
r(O1–Ti)2.1282.0681.9721.9511.956
r(O2–Ti)1.9762.0391.9992.0532.058
r(O3–Ti)-2.096---
<(O–C–C–C)2.031.8−4.30.2−1.3
<(C343;3Ti)8.352.281.788.885.8
combination of LANL2DZ on Ti and 6-31G(d,p) on O, C, N, and H atoms.
Table 5. The adsorption energy of the C343-TinO2n+2H4 systems (n = 14, 24), defined as the difference between the calculated energy of the bound system and the sum of the energies for the separated constituents (isolated Ti2nOn+2H4 and C343), (in eV), calculated by the DFT/B3LYP method with 3-21G(d) or LANL2DZ basis sets, in vacuum or in water solvent.
Table 5. The adsorption energy of the C343-TinO2n+2H4 systems (n = 14, 24), defined as the difference between the calculated energy of the bound system and the sum of the energies for the separated constituents (isolated Ti2nOn+2H4 and C343), (in eV), calculated by the DFT/B3LYP method with 3-21G(d) or LANL2DZ basis sets, in vacuum or in water solvent.
MethodC343-Ti14O30H4, Two-BondsC343-Ti14O30H4, Three-BondsC343-Ti24O50H4
B3LYP/3-21G(d) vacuum−2.473−3.075−3.254
B3LYP/3-21G(d) solvent−1.961−2.574−2.835
B3LYP/LANL2DZ solvent0.429−0.242−1.188
Table 6. The structural parameters of the model system consisting of the Ti24O50H4 clusters with an adsorbed OMCD1 molecule, and of the bare cluster, obtained after geometry optimization by DFT/B3LYP/3-21G(d) in vacuum. Geometrical parameters describing the adsorption modes: bond lengths (in Å), the torsion angle of the carboxyl group relative to the indoline plane, and the dihedral angle defined by the indoline plane and the plane of the three neighbor Ti atoms on the cluster surface (both in deg).
Table 6. The structural parameters of the model system consisting of the Ti24O50H4 clusters with an adsorbed OMCD1 molecule, and of the bare cluster, obtained after geometry optimization by DFT/B3LYP/3-21G(d) in vacuum. Geometrical parameters describing the adsorption modes: bond lengths (in Å), the torsion angle of the carboxyl group relative to the indoline plane, and the dihedral angle defined by the indoline plane and the plane of the three neighbor Ti atoms on the cluster surface (both in deg).
ParameterOMCD1 on Ti24O50H4Ti24O50H4
r(Ti–O)1.8801.876
σr(Ti–O)0.1020.100
Width (Ti…Ti)12.2012.18
Width (O…O)13.2213.28
Length (Ti…Ti)7.217.20
Length (O…O)8.158.16
r(O1–Ti)2.012-
r(O2–Ti)1.989-
<(O–C–C–C)−4.2-
<(Phenyl;3Ti)87.2-
Table 7. The structural parameters of the model system consisting of the TinO2n+2H4 clusters (n = 14, 24) with an adsorbed penicillin V molecule, obtained after geometry optimization by DFT/B3LYP/3-21G(d) in vacuum. Geometrical parameters describing the adsorption configurations: bond lengths (in Å) and torsion angles (in deg.) of the PV relative to the Ti–O bond, to the carboxyl group, and of the rest of the PV with respect to the phenyl group 1.
Table 7. The structural parameters of the model system consisting of the TinO2n+2H4 clusters (n = 14, 24) with an adsorbed penicillin V molecule, obtained after geometry optimization by DFT/B3LYP/3-21G(d) in vacuum. Geometrical parameters describing the adsorption configurations: bond lengths (in Å) and torsion angles (in deg.) of the PV relative to the Ti–O bond, to the carboxyl group, and of the rest of the PV with respect to the phenyl group 1.
ParameterPV on Ti14O30H4PV on Ti24O50H4PV on Ti34O70H4PV on Ti44O90H4PV on Ti54O110H4
r(Ti–O)1.8221.8841.8961.9001.903
σr(Ti–O)0.1800.1040.1040.1030.107
Width (Ti…Ti)12.3412.2111.8711.9512.01
Width (O…O) 13.2713.3413.3513.32
Length (Ti…Ti)3.637.4010.9014.5318.11
Length (O…O)4.818.8211.6415.3018.80
r(O1–Ti)2.0402.0732.0772.0412.044
r(O2–Ti)2.0352.1022.1002.1192.110
<(Ti–O–C–C)159.4129.1129.7130.0130.4
<(O–C–C–N)25.70.81.0−0.8−0.5
<(C–C–O–C)0.6−30.479.947.442.3

Share and Cite

MDPI and ACS Style

Oprea, C.I.; Gîrțu, M.A. Structure and Electronic Properties of TiO2 Nanoclusters and Dye–Nanocluster Systems Appropriate to Model Hybrid Photovoltaic or Photocatalytic Applications. Nanomaterials 2019, 9, 357. https://doi.org/10.3390/nano9030357

AMA Style

Oprea CI, Gîrțu MA. Structure and Electronic Properties of TiO2 Nanoclusters and Dye–Nanocluster Systems Appropriate to Model Hybrid Photovoltaic or Photocatalytic Applications. Nanomaterials. 2019; 9(3):357. https://doi.org/10.3390/nano9030357

Chicago/Turabian Style

Oprea, Corneliu I., and Mihai A. Gîrțu. 2019. "Structure and Electronic Properties of TiO2 Nanoclusters and Dye–Nanocluster Systems Appropriate to Model Hybrid Photovoltaic or Photocatalytic Applications" Nanomaterials 9, no. 3: 357. https://doi.org/10.3390/nano9030357

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop