Exchange Functionals and Basis Sets for Density Functional Theory Studies of Water Splitting on Selected ZnO Nanocluster Catalysts

In this communication, we use density functional theory (DFT) to study the structural (geometry) and electronic properties (vertical detachment energy and electron affinity) of ZnO monomers and dimers that can be used to form ZnO clusters of different sizes, with a view to adapting one or more of them as catalysts or photocatalysts, standing alone or on suitable substrates like graphene, to split water. We also investigate different pairs of exchange functionals and basis sets to optimize their choice in our DFT calculations and to compare the singlet–triplet energy gaps of small ZnO clusters of different sizes to select an optimal cluster size for water splitting. We find that the B3LYP/DGDZVP2 exchange functional/basis set is a reliable combination for use with DFT to calculate the geometry and electronic properties of small ZnO nanoclusters from among several other combinations of exchange functionals and basis sets. Comparisons of the singlet–triplet energy gaps show that the trimer (ZnO)3 has an energy gap of 58.66 k cal/mol. which is approximately equal to the energy of a visible photon at a wavelength of 500 nm, and a HOMO–LUMO gap of 4.4 eV, making it a suitable choice of photocatalyst for the oxidation of water from among six (ZnO)n nanoclusters of monomers, with n ranging from 1 to 6. We used this exchange functional/basis set to study the structural and energetic details of hydration and hydrolysis of water absorbed on the (ZnO)3 nanocatalyst and calculated the corresponding potential energy profiles to identify three sets of singlet–triplet pathways for water splitting. Detailed study of a pathway showed that oxygen is produced after hydrogen, and the rate-determining step is the formation of hydrogen.


■ INTRODUCTION
Metal oxides have been used in chemistry, metallurgy, and engineering over many decades. Depending on the metal, the oxides can be utilized as catalysts, 1,2 sensors, 3,4 medical additives, 5 and in solar cells. 6,7 Scaling down from macro to nano, both the physical and chemical properties of each metal oxide are changed, and nanoscaled metal oxide clusters become a promising material to push the limits of their applications as catalysts. 8−12 Accordingly, there has been a growing interest in the theoretical and experimental study of metal oxide nanoclusters. 13−17 Zinc oxide (ZnO) is such a material with potential applications of interest in many industries including paint, cosmetics, and pharmaceuticals. 18−27 Zinc oxide is a semiconductor with a wide band gap of 3.4 eV at 300 K and has unique magnetic, optical, and electronic properties. 28−30 Moreover, those properties can change considerably when the particle dimension and size change to the nanoscale.
The structural and electronic properties of ZnO clusters have been studied experimentally and theoretically. Bovhyra et al. applied density functional theory (DFT) to investigate the structural and electronic properties of (ZnO) n (n = 34, 60) nanoclusters. 31 Their investigations predicted that the most stable (ZnO) 34 nanoclusters were fullerene-like hollow structures, and the most energetically favorable (ZnO) 60 had a sodalite-type structure made with seven (ZnO) 12 clusters and common quadrangle edges. The energy range of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) for those clusters is 1.93−2.40 eV. A study by Woodley et al. reported on the stable and low energy metastable structures of (ZnO) n clusters, when n = 1−32 using evolutionary algorithm techniques. 32 They found 2D rings for n = 2−7 that were less than 0.8 nm in diameter and spheroidal bubbles up to n = 32 that were less than 1.2 nm in diameter.
Another theoretical study by Szakacs et al. showed the relative thermodynamic stability of two different ZnO nanocluster shapes, namely, needles and plates. The electronic structure calculations of ZnO needle and plate nanoclusters show the (ZnO) n , n = 6, 9, 12, 15, and 18 needles, have higher stability that increased with the number of stacked rings. They also studied Zn 3 O 3 as a successful building block to form larger clusters rather than by the addition of individual or pairs of Zn and O atoms. 33 An extensive theoretical study on the structural and electronic properties of (ZnO) n , n ≤ 168, was carried out by Chen et al. 34 Their major findings include the presence of new magic numbers for (ZnO) n for structures at n = 78, 100, 132, and 168 with the particle stability increasing for structures with single to double and triple layered octahedral cage morphologies. Chen et al. summarized the stability of small ZnO clusters. The lowest energy isomers for (ZnO) n , n = 3−5, were predicted to be ring-like 2-dimensional (2D) structures that had higher stability compared to the open and 3D structures.
To fill the gap in information lacking in previous investigations of anionic ZnO clusters, Gunaratne et al. investigated the relative stabilities and electronic properties of small ZnO clusters using density functional theory (DFT) with both B3LYP and PBE exchange functionals and 6-311++G(3d), aug-cc-pvQZ and LanL2DZ basis sets. 35 They compared their calculations with experimental results, reported the mass spectra of anionic Zn n O m (n = 1−6 and m = 1−7) clusters, and predicted the ground state geometries and stabilities of Zn 3 O m and Zn 3 O m − (m = 3−5). They also determined the adiabatic detachment energies (ADEs) and vertical detachment energies (VDEs) of the anions Zn 3 O 3 − and Zn 3 O 4 − using photoelectron spectroscopy, compared them with their theoretical results, and concluded that the (ZnO) 3 cluster is more stable than its oxygen rich associates. Dheivamalar and Banu 36 very recently reported on absorption mechanism and structural and electronic properties of pyrolle on (ZnO) 3 and (ZnO) 6 using DFT. They showed that pyrolle adsorption significantly changes the electronic properties of the clusters, and they concluded that Zn 3 O 3 nanoclusters are more effective catalysts in solar cells than Zn 6 O 6 .
In this paper, we set the stage for our current and future investigations of selected ZnO nanoclusters as catalysts to split water. We first examined the structural and electronic properties of ZnO monomer and dimer, using 30 different pairs of exchange functionals and basis sets to choose an optimal pair in the theoretical analysis of the structural and electronic properties of ZnO nanoclusters using DFT. Where possible, we compared our results with CCSD(T) calculations using the same basis sets. We calculated the energies of singlet and triplet states of (ZnO) n nanoclusters for n = 1−6 and determined singlet triplet energy gaps with the same functionals and basis sets as a first step toward choosing an optimal cluster size of n = 3, as a potential catalyst to produce H 2 and O 2 gas by splitting water. ZnO is cheaper than TiO 2 , which is widely used as a photocatalyst to split water, but its relative efficiency is unknown. Theoretical calculations at the atomic and molecular level would help us understand the mechanism of water splitting and could lead to better choices and designs of catalysts for water splitting. For example, knowledge of the rate-determining step can be used to design catalysts with a lower activation energy for that step.
The production of H 2 and O 2 from the dissociation or splitting of water is an endothermic reaction that involves two water molecules.
The free energy of formation of H 2 O gas is −54.6 kcal/mol, and the heat of formation is −57.8 ± 0.02 kcal/mol at 298 K temperature, which is comparable to the energy of two photons at 500 nm. 37 This is exploited in photocatalytic water splitting, which captures the energy of photons to split water.
When the photocatalyst is a semiconductor, with the minimal band gap to absorb photons, water splitting occurs indirectly through the photoexcitation of electrons (e − ) from the valence band (VB) to the conduction band (CB) of the catalyst, leaving behind holes (h + ) in the VB.
Charge separation is followed by recombination or transportation to the catalytic surface when water is oxidized by holes to form protons H + and O 2 molecules that are desorbed from the surface.
The protons H + released are reduced by electrons to form H 2 molecules that are also desorbed from the catalytic surface.
In this scenario, hydrogen molecules are formed after oxygen molecules, and the overall efficiency is determined by efficiencies of light harvesting, charge separations of holes and electrons, and the oxidation of water by holes to form oxygen and reduction of hydrogen ions by electrons to form hydrogen. They compete with the recombination of electrons and holes further reducing the efficiency of water splitting. 38 Photochemical water splitting is consequently slower than electrochemical water splitting but has the advantage that it is less expensive and needs no energy other than sunlight.
Here, we used DFT to study the hydration and hydrolysis of two water molecules adsorbed successively on both the singlet and triplet states of (ZnO) 3 nanoclusters that lead to the production of H 2 and O 2 during different stages of water splitting while remaining almost exclusively in these states. The triplet state in DFT is treated as an excited singlet state, and we excluded promotion of electrons from the ground state to excited states. We calculated the potential energy diagrams and elucidated the mechanism of water splitting, which, in this scenario, is not necessarily the same as in photocatalytic water splitting. To the best of our knowledge, detailed theoretical studies of water splitting using ZnO nanoclusters as catalysts have not been reported earlier. Kenmore and Biedermann reported DFT studies of water aggregation and splitting on the ZnO(1010) surface site. 39 There are also accounts of experimental studies of the photocatalytic ability of ZnO nanoparticles. 40 Previously, Fang and Dixon 37 described extensive computational (DFT and CCSD(T)) studies of hydrogen and oxygen production from water splitting on small metal oxide clusters (MO 2 ) n , where n = 1−3 and the metal M is Ti, Zr, or Hf. Nanocatalysts for water splitting have been reviewed recently by Grewe et al. 41 and by Jaffari et al. 42

■ COMPUTATIONAL METHOD
We chose six different exchange functionals B3LYP, B3PW91, PBE, BVP86, BP86, and M06 exchange functionals and 30 different basis sets shown in Figure 1 in our DFT calculations of the structural and electronic properties of zinc oxide nanoclusters as potential catalysts.
We optimized the ZnO monomer and dimer structures and energies using DFT and the CCSD(T) theory with same basis sets. We also analyzed bond lengths, bond angles, and Mulliken charges on atoms that characterize the structural properties and the electron affinity (EA), vertical detachment energy (VDE), and energy representing the electronic properties of each optimized structure. Additionally we compared the HOMO− LUMO gap, the difference between singlet and triplet energies and calculated the infrared spectra of (ZnO) n clusters, with n = 1−6, to evaluate candidates with suitable photocatalytic properties. We compared our calculated values with experimental results 35 to identify the best basis set candidates. We studied the densities of state (DOS), the HOMO and LUMO profiles of (ZnO) n , n = 1−6, and (ZnO) n when one water molecule absorbed on each cluster. The DOS plots were calculated using the GaussSum 3.0 Program. 43 We also calculated the infrared (IR) spectra of each of the six (ZnO) n , n = 1−6 nanoclusters and the same (ZnO) n nanoclusters with one water molecule absorbed on each, compared them with previous work where possible, and verified the absence of negative frequencies, confirming their stability. All transition states of the reaction pathways were calculated using the Berny approach 44 and further studied using the intrinsic reaction coordinates method 45,46 using the selected basis set.
All geometry optimization as well as structural and electronic property calculations were carried out using the Gaussian 16 software package 47 and GaussView 6 graphical interface. 48 ■ RESULTS AND DISCUSSION Comparison of Structural Properties. Geometry optimization is the process of finding the positions of nuclei for which the potential energy is minimized. We optimized the ZnO monomer and dimer and compared their structural properties using DFT and CCSD(T) methods with 30 different basis sets. Figure 2 illustrates the optimized monomer (a) and dimer (b) structures of ZnO using DFT with the B3LYP/DGDZVP2 combination of exchange functional and basis set. We selected the bond lengths (2Zn−1O for both monomer and dimer as shown in the Figure 2), dimer bond angles (2Zn−1O−3Zn and 1O−2Zn−4O), and Mulliken charges on each atom labeled in the Supporting Information (SI) Figures S1 and S2 as characteristic structural properties.
The monomer and dimer bond lengths, calculated for 30 different basis sets and 6 different exchange functionals in DFT and the corresponding results from CCSD(T) theory using the same basis sets, are shown in Figures 3 and 4 with the details numbers recorded in Tables S17 and S18, respectively. Other results including the optimized energies of the neutral and anion monomer and dimer are tabulated in the SI.    We compared the Zn−O bond lengths of the monomer and dimer with experimental values (the bond length of Zn−O is 1.719 and 1.787 Å for its anion). 35 The bond lengths for the monomer, calculated using B3LYP functional in DFT, vary widely from 1.609−1.818 Å. However, the 6-311+G(d), 6-311+ +G(d), 6-311++G(d,p), and DGDZVP2 basis sets predict bond lengths in good agreement with the experimental value of 1.719 Å to within ±0.01 Å. Monomer bond lengths calculated with the B3LYP functional are comparably closer to experiment and to within the first decimal using the CCSD(T) method. The Zn−O bond lengths in the dimer are longer than the corresponding bond length in the monomer by about 0.13 Å in both the DFT and CCSD(T) theory. The metal−metal (Zn−Zn) bond distance in the dimer, optimized at the B3LYP/DGDZVP2 level, is 2.373 Å.
We next examined the bond angles, 2Zn−1O−3Zn and 1O− 2Zn−4O, obtained from DFT using the B3LYP functional and the 30 different basis sets. They are displayed in Table S16 of the SI. The 2Zn−1O−3Zn angles range between 73.17°and 79.59°, while the 1O−2Zn−4O angles are larger and between 100.06°a nd 106.86°. Overall they are in excellent agreement with CCSD theory calculations to within 0.6°for the same basis sets as displayed in Figure 1.
A theoretical study by Ugalde et.al on the global minimum for small ZnO clusters concludes that the 1O−2Zn−4O bond angle is 103.7°, and the Zn−O bond length is 1.892 Å for the ZnO dimer using DFT with the B3LYP functional and a relativistic compact effective core potential and the shared-exponent basis set of Stevens, Krauss, Basch, and Jasien (SKBJ). 49 Their theoretical calculations are in fair agreement with our DFT results (102.70°and 1.900 Å) using the same exchange functional and the DGDZVP2 basis set. Our CCSD(T) calculation with the same basis set predicts 103.3°for the dimer bond angle and 1.917 Å for the bond length. These differences are small and may reflect corresponding deficiencies in our choice of unmixed basis sets and exchange functional. In the remainder of this paper, we used B3LYP/DGDZVP2 exchange functional and basis set.
Comparison of Electronic Properties. Focusing on the electronic properties of the ZnO nanoclusters, we compared the experimental 50,51 electron affinities (EAs) and vertical detachment energies (VDEs) of both ZnO monomer and dimer with our DFT and CCSD(T) calculations for the same basis sets. The EAs and VDEs were calculated from the optimized energies using eqs 6 and 7. 52 = The EA of a molecule is the difference in energy between the neutral molecule and the negative ion of the same molecule in their ground states. Positive EA values indicate thermodynamically stable negative ions. Our results are shown in Figure 5 with the details in Table S19 Table  S19, show that the CCSD(T) theory with the 6-31+G, 6-31+    +G, 6-31+G(3df,2p), and aug-cc-PVTZ basis sets predict EA that agree to within 0.06 eV of the experimental value. Our DFT calculations using the B3LYP functional and DGDZVP2 or CEP-121G basis sets are also within the same error with respect to the experiment. The B3PW91 functional with the LANL2DZ basis set, the PBE functional with the cc-PVQZ and aug-cc-PVDZ basis sets, and the M06 functional with the cc-PVQC basis sets are also in accord with the experimental results to within the same error. The EA values of the ZnO dimer shown in Figure 6 and Table S20 are Table S22). We next investigated the choice of a suitable ZnO cluster size as a catalyst for water splitting. We chose six different sizes of ZnO nanoclusters, (ZnO) n , n = 1−6, whose equilibrium structures obtained using DFT and the exchange functional/ basis set (B3LYP/DGDZVP2) are displayed in Figure 9. The Zn−Zn bonds in the ring structures n > 2 are omitted for clarity.
We identified a suitable catalyst to assist in H 2 and O 2 production by splitting water. The catalyst should be able to act alone in the absence of light and as a photocatalyst in the presence of light.
The role of a photocatalyst is to absorb the energy of visible light and elevate the ground state electrons of a molecule to an excited state that enables the reaction to occur faster through a lowering of the activation barrier. Studying the excited state of the desired reaction is thus important. But, excited open shell singlet states are difficult or nearly impossible to investigate using current DFT methods. Since the energy difference between the excited singlet state and the triplet state is comparatively small, we considered the triplet state as an excited singlet state with spin one.
The optimized energies of the singlet and the first excited triplet state of (ZnO) n , n = 1−6, were calculated using the selected methods. The details are shown in the SI (Tables S11 and S12). The energy differences of the singlet and the first excited triplet state in kcal/mol are displayed in Table 1.
On examining the calculated energy differences between the singlet and triplet states of small (ZnO) n clusters, we observed that they are smaller than 57.8 kcal/mol when n = 1 and 2 and greater than 57.8 kcal/mol for n = 4 and 5. Clusters with n = 3 and 6 have energy gaps closer to the energy of two visible photons that can photocatalyze the splitting of two water molecules to produce H 2 and O 2 . These results should be useful in designing further improvements to the performance of catalysts using appropriate substrates like graphene and graphene oxide. We chose the n = 3 as the preferred cluster size for a ZnO nanocatalyst with the desired properties that match the band gap, while taking into account the lower computational costs of using smaller clusters of ZnO monomers as catalysts.
The energy gap between HOMO and LUMO is an important parameter that identifies cluster stability and possible use as a photocatalyst. We calculated the energy gap (E g ) from eq 8 and DFT, where the E HOMO is the energy of the HOMO and E LUMO is the energy of LUMO. We compared the HOMO−LUMO gap of the six ZnO nanoclusters in Table 2 using the B3LYP functional and the four 6-311++G, DGDZVP, DGDZVP2, and aug-cc-PVDZ basis sets. Our results for the HOMO−LUMO gap of desired ZnO clusters using the B3LYP/DGDZVP2 combination of exchange functional and basis set in DFT are shown in Figure 10, and other relevant data are in the SI.
Cluster sizes 4 and 5 have relatively higher HOMO−LUMO gaps (4.66 and 4.79 eV, respectively using the B3LYP/ DGDZVP2 method) indicating that they may be relatively less efficient than the other neutral clusters. (ZnO) 3 and (ZnO) 6 have lower HOMO−LUMO energy gaps (4.42 and 3.76 eV using the same method) compared to the fourth and fifth clusters, possibly indicating that the n = 3 and 6 (ZnO) n clusters may have a comparably higher reactivity. When considering the photocatalytic water splitting reaction to produce H 2 , the minimum photon energy required for the overall reaction is 1.23 eV, which indicates that all the small ZnO nanoclusters studied are possible candidates to catalyze the photosplitting of water from DFT using the B3LYP exchange functional and four different basis sets.   The differences between the DFT results for the same nanocluster and four different basis sets using DFT are small. However, the absolute values could be in error, since DFT is known to underestimate the band gap. Figure 11 shows the DOS plot along with the HOMO, LUMO profile of (ZnO) 3    than on Zn atoms but the acceptor LUMO orbitals are dispersed over all the atoms of the cluster.
H 2 O Adsorption on the (ZnO) 3 Cluster. The charges on the atoms of nanocluster and water molecule steer the water molecule to the Zn adsorption site on the nanocluster to form a bond. To help clarify the steps, we used DFT at the B3LYP/ DGDZVP2 level to calculate the Mulliken charges on each atom of an isolated water molecule and the atoms of the (ZnO) 3 trimer before and after the adsorption of a water molecule. The charges are displayed in Table 3 and show why the postively charged zinc atom is the preferred adsorption site for the negatively charged oxygen atom of water. The atomic charge distribution also plays an important role in explaining not only bond formation but also bond breaking in the hydration and hydrolysis of a water molecule, leading eventually to water splitting. In the (ZnO) 3 Figure 12 shows the DOS plot and the HOMO, LUMO profile of (ZnO) 3 with H 2 O, using the B3LYP/DGDZVP2 level of theory. The DOS shows the HOMO and LUMO energies as −6.97 and −2.32 eV, respectively, and energy gap calculated using eq 3 between HOMO and LUMO is 4.42 eV This yields the Fermi level as −4.64 eV. When we introduce a water molecule to the (ZnO) 3 cluster, it tends to make bonds with Zn and O atoms in the cluster (as shown in Figure 4). The HOMO−LUMO gap increased when the water molecule absorbed on (ZnO) 3 cluster from 4.42 to 4.65 eV. We also calculated the vibrational frequencies in order to check the stability of the optimized geometries. Infrared spectra for (ZnO) n , n = 1−3 and (ZnO) 3 , when one water molecule absorbed on it are shown in Figure 13. The IR spectra for (ZnO) 4 , (ZnO) 5 , (ZnO) 6 , and all the relevant frequencies for  There were no imaginary vibration frequencies in our study of any of the clusters. The monomer ZnO has one wide peak at 708 cm −1 . The dimer ZnO has three major peaks at 213, 441, and 566 cm −1 frequencies. (ZnO) 3 has a major peak at 648 cm −1 and a minor peak at 221 cm −1 , (ZnO) 4 has a major peak at 735 cm −1 and two minor peaks at 176 and 225 cm −1 , (ZnO) 5 has one major peak at 801 cm −1 and three minor peaks at 174, 221, and 518 cm −1 , and (ZnO) 6 has two major peaks at 385 and 581 cm −1 and a minor peak at 491 cm −1 frequencies. When the water molecule absorbed on (ZnO) 3 , the IR peaks shifted toward higher frequency values.
The optimized energies of the singlet and the first excited triplet state of (ZnO) n , n = 1−6, are calculated using the selected methods. The details are shown in the SI. The energy differences of the singlet and the first excited triplet state are displayed in Table 1, and the detailed energies are in Tables S11 and S12.
Hydration and Hydrolysis Reaction Pathways of (ZnO) 3 . First Hydrolysis. The hydrolysis of a water molecule on the (ZnO) 3 catalyst occurs in two stages. In the first stage, a water molecule is adsorbed with its oxygen bonded to any one of three zinc atoms of the (ZnO) 3 nanocatalyst followed by the molecule tipping over to enable proton transfer of a H atom of the same water molecule to an oxygen atom adjacent to the zinc atom already bonded to the water molecule. Dissociation of the absorbed water molecule leaves a hydrogen atom and a OH group bonded to adjacent O and Zn (sites 4 and 5) on the ring of the (ZnO) 3 cluster. The pathways for adsorption, hydration, and hydrolysis are similar in the singlet and triplet states of the catalyst, but the  Second Hydrolysis. There are two vacant Zn sites (a and b) at which a second water molecule can be absorbed on the already hydrated (ZnO) 3 nanocluster. Site a has only one vacant ring oxygen next to it that can accept a H atom from the dissociation of the second water molecule, while site b has two. We characterized the two pathways for hydration and hydrolysis from site b as b1 and b2 and the hydration and hydrolysis from Zn site a as pathway a. Overall, there are six pathways (three sets of singlet and triplet pathways) for the hydrolysis of a second water molecule.    Table 4.    hydrolysis of the second water molecule. Hydrolysis of two water molecules occurs readily on the catalytic surface in the singlet and triplet states in the absence of light. It is not a ratedetermining step in the splitting of water Following our study of the hydrolysis of two water molecules on the (ZnO) 3 nanocatalyst, we turn our attention to further details of the mechanism of water splitting to produce H 2 and O 2 formation and release. Since pathways a and b2 lead to the same hydrolyzed products, we need to consider only the twos sets a and b1 pathways beyond this point to understand the mechanism of formation of H 2 and O 2 during water splitting.
Light is not explicitly considered in our calculations of the pathways to split water using a (ZnO) 3 nanocatalyst, but in the presence of light of appropriate wavelength, one or more photons corresponding to wavelengths of 500 nm can, in principle, promote transitions from the single to triplet states of the nanocatalyst at the beginning or at any other stage of the reaction depending on the energy gap between the states.
We calculated the potential energy surfaces for the formation and release of H 2 molecules and O 2 molecule along pathway b1 after the hydrolysis of two water molecules on the surface of ZnO nanoclusters. The activated and intermediate states as well as the reaction pathways were determined using intrinsic reaction coordinate (IRC) methods. 45,46 Figure 21 shows the potential energy surface for formation and release of the first H 2 molecule along this pathway during water splitting, and the corresponding relative energies are recorded in Table 5.
The formation of the first hydrogen molecule from the hydrolyzed product of pathway b1 requires prior rearrangement of the hydrolyzed ring to bring a ring oxygen close to the oxygen of Zn−O−H group, followed by Zn−H and O−O bonds formed within the cluster after surmounting two barriers with high activation energies of 120.8 and 49.9 kcl/mol, which brings the hydrogen atoms on the −O−H group close to another H atom bonded to a ring oxygen. This is a favorable arrangement that enables bond formation between hydrogen atoms to produce a hydrogen molecule. Formation of the first hydrogen molecule is clearly a rate-determining step for water splitting. Figure 22 displays the potential energy surface for the release of a second hydrogen molecule ion followed by the formation and desorption of an oxygen molecule in the triplet state to complete the splitting of two water molecules and the regeneration of the (ZnO) 3 nanocatalyst in its original state. The barrier for the formation of the second hydrogen molecule is lower than for the first hydrogen molecule, but it is still high in contrast to the relatively low barrier for oxygen formation, which follows the release of hydrogen. The energetic details are recorded in Table 6. The last steps involve two curve crossings, which may require mixed states for a more accurate representation. In the final steps, the oxygen molecule released is in its most stable triplet state and the hydrogen molecules are in their singlet states. The overall catalyzed reaction to split water using a (ZnO) 3 nanocatalyst is   Figure 22 show that the energy change for the splitting of water represented in eq 11 is 106.8 kcal/mol. The literature value of the enthalpy change is 115.6 kcal/mol at 298 K.

■ CONCLUSION
We investigated the structural and electronic properties of ZnO monomer and dimer and the vibrational frequencies of both to  choose suitable basis sets and exchange functionals to study the catalytic effects of ZnO nanoclusters on water splitting. By comparing our calculated values with experiments, we found that the B3LYP/DGDZVP2 combination of exchange functional and basis set would be a reliable and optimal combination to use in calculating the properties of ZnO nanoclusters using DFT with a relatively low computational time compared to other methods (e.g., CCSD(T)). We also studied the singlet− triplet energy difference and HOMO−LUMO gap of (ZnO) n , with n = 1−6, to select the best cluster sizes of (ZnO) n as photocatalysts. From the singlet−triplet energy comparisons, both the n = 3 and n = 6 nanocluster systems are possible photocatalysts for H 2 and O 2 production in a water splitting reaction. We selected (ZnO) 3 rather than (ZnO) 6 as a suitable nanocatalyst to study water splitting, considering the complexity of the reaction and reduction in computational time. Successive hydrolysis of two water molecules on the catalytic surface occur readily with low activation barriers. Water splitting occurs after hydrolysis with the release of hydrogen gas before oxygen. The rate-determining step is controlled by high barriers for hydrogen formation.
■ ASSOCIATED CONTENT