Zongtang Fang

Computational study of H2 and O2 production from water splitting by small (MO2)n clusters (M = Ti, Zr, Hf).

Coupled cluster [CCSD(T)] theory and density functional theory (DFT) have been used to study the production of H2 and O2 from hydrolysis products generated from H2O addition to (MO2)n (M = Ti, Zr, Hf, n = 1-3) clusters on both the lowest singlet and triplet potential energy surfaces (PESs). H2 production occurs via the formation of an M-H containing intermediate followed by H-H recombination and H2 desorption from M(n)O(2n)(OH)2 and M(n)O(2n+2). The hydrogen transfer reactions to form the M-H bond are the rate determining steps and can be considered to be proton coupled, electron transfer (PCET) reactions with one or two electrons being transferred. Oxygen is produced by breaking two weak M-O bonds in an atomic oxygen saturated metal oxide from an M(n)O(2n)•O2 intermediate. On the triplet PES, the activation energies for the first and second H transfer to the metal are calculated to be ~10 to 50 kcal/mol and ~75 to 90 kcal/mol depending on the size of the clusters and the metal. The barriers on the singlet surface for the first and the second H transfer are predicted to be 110 to 140 kcal/mol, in general larger than the H-O bond dissociation energy. The activation barriers for the step of H-H recombination are 15 to 50 kcal/mol, and the H2 desorption energies are less than 10 kcal/mol on the singlet and triplet PESs. The oxygen desorption energies follow the order Ti < Zr < Hf for the triplets and Ti < Zr ≈ Hf for the singlets. The oxygen desorption energy is approximately independent of the size of the cluster for the same metal. The water splitting reactions prefer to take place on the triplet surface. A low excess potential energy is needed to generate 2H2 and O2 from 2H2O after the endothermicity of the reaction is overcome on the triplet PES.

Prediction of Bond Dissociation Energies/Heats of Formation for Diatomic Transition Metal Compounds: CCSD(T) Works

It was recently reported ( J. Chem. Theory Comput. 2015 , 11 , 2036 - 2052 ) that the coupled cluster singles and doubles with perturbative triples method, CCSD(T), should not be used as a benchmark tool for the prediction of dissociation energies (heats of formation) for the first row transition metal diatomics based on a comparison with the experimental thermodynamic values for a set of 20 diatomics. In the present work the bond dissociation energies as well as the heats of formation for those diatomics have been calculated by the Feller-Peterson-Dixon approach at the CCSD(T)/complete basis set (CBS) level of theory including scalar relativistic corrections and correlation of the outer shell of core electrons in addition to the valence electrons. Revised experimental values for the hydrides are presented that are based on new heterolytic R-H bond dissociation energies, which are needed for analysis of the mass spectrometry experiments. The agreement between the calculated bond dissociation energies and the revised experimental values of the hydrides is good. Good agreement of the calculated bond dissociation energies/heats of formation is also found for most of the chlorides, oxides, and sulfides given the experimental error bars from experiment and those of the transition metal atoms in the gas phase. Thus, reliable results can be achieved by the CCSD(T) method at the CBS limit. The use of PW91 orbitals for the CCSD(T) calculations improves the predictions for some compounds with large T1 diagnostics at the HF-CCSD(T) level. The optimized bond distances and calculated vibrational frequencies for the diatomics also agree well with the available experimental values.

Use of Improved Orbitals for CCSD(T) Calculations for Predicting Heats of Formation of Group IV and Group VI Metal Oxide Monomers and Dimers and UCl6

The prediction of the heats of formation of group IV and group VI metal oxide monomers and dimers with the coupled cluster CCSD(T) method has been improved by using Kohn-Sham density functional theory (DFT) and Brueckner orbitals for the initial wave function. The valence and core-valence contributions to the total atomization energies for the CrO3 monomer and dimer are predicted to be significantly larger than when using the Hartree-Fock (HF) orbitals. The predicted heat of formation of CrO3 with CCSD(T)/PW91 is consistent with previous calculations including high-order corrections beyond CCSD(T) and agrees well with the experiment. The improved heats of formation with the DFT and Brueckner orbitals are due to these orbitals being closer to the actual orbitals. Pure DFT functionals perform slightly better than the hybrid B3LYP functional due to the presence of exact exchange in the hybrid functional. Comparable heats of formation for TiO2 and the second- and the third-row group IV and group VI metal oxides are predicted well using either the DFT PW91 orbitals, Brueckner orbitals, or HF orbitals. The normalized clustering energies for the dimers are consistent with our previous work except for a larger value predicted for Cr2O6. The prediction of the reaction energy for UF6 + 3Cl2 → UCl6 + 3F2 was significantly improved with the use of DFT or Brueckner orbitals as compared to HF orbitals.

Benchmark Calculations of Energetic Properties of Groups 4 and 6 Transition Metal Oxide Nanoclusters Including Comparison to Density Functional Theory.

The heats of formation and the normalized clustering energies (NCEs) for the group 4 and group 6 transition metal oxide (TMO) trimers and tetramers have been calculated by the Feller-Peterson-Dixon (FPD) method. The heats of formation predicted by the FPD method do not differ much from those previously derived from the NCEs at the CCSD(T)/aT level except for the CrO3 nanoclusters. New and improved heats of formation for Cr3O9 and Cr4O12 were obtained using PW91 orbitals instead of Hartree-Fock (HF) orbitals. Diffuse functions are necessary to predict accurate heats of formation. The fluoride affinities (FAs) are calculated with the CCSD(T) method. The relative energies (REs) of different isomers, NCEs, electron affinities (EAs), and FAs of (MO2)n (M = Ti, Zr, Hf, n = 1-4) and (MO3)n (M = Cr, Mo, W, n = 1-3) clusters have been benchmarked with 55 exchange-correlation density functional theory (DFT) functionals including both pure and hybrid types. The absolute errors of the DFT results are mostly less than ±10 kcal/mol for the NCEs and the EAs and less than ±15 kcal/mol for the FAs. Hybrid functionals usually perform better than the pure functionals for the REs and NCEs. The performance of the two types of functionals in predicting EAs and FAs is comparable. The B1B95 and PBE1PBE functionals provide reliable energetic properties for most isomers. Long range corrected pure functionals usually give poor FAs. The standard deviation of the absolute error is always close to the mean errors, and the probability distributions of the DFT errors are often not Gaussian (normal). The breadth of the distribution of errors and the maximum probability are dependent on the energy property and the isomer.

Properties of Cerium Hydroxides from Matrix Infrared Spectra and Electronic Structure Calculations.

Reactions of laser ablated cerium atoms with hydrogen peroxide or hydrogen and oxygen mixtures diluted in argon and condensed at 4 K produced the Ce(OH)3 and Ce(OH)2 molecules and Ce(OH)2(+) cation as major products. Additional minor products were identified as the Ce(OH)4, HCeO, and OCeOH molecules. These new species were identified from their matrix infrared spectra with D2O2, D2, and (18)O2 isotopic substitution and correlating observed frequencies with values calculated by density functional theory. We find that the amounts of Ce(OH)3 and of the Ce(OH)2(+) cation increase on UV (λ > 220 nm) photolysis, while Ce(OH)2, Ce(OH)4, and HCeO are photosensitive. The observed major species for Ce are in the +III or +II oxidation state, and the minor product, Ce(OH)4, is in the +IV oxidation state. The calculations for the vibrational frequencies with the B3LYP functional agree well with the experiment. The NBO analysis shows significant backbonding to the metal 4f and 5d orbitals for the closed shell species. Most open shell species have the excess spin in the 4f with paired spin in the 5d due to backbonding. The heats of formation of the observed species were derived from the available data from experiment and the calculated reaction energies. The major products in this study are different from similar reactions for Th where the tetrahydroxide was the major species.

Structures and Properties of the Products of the Reaction of Lanthanide Atoms with H2O: Dominance of the +II Oxidation State

The reactions of lanthanides with H2O have been studied using density functional theory with the B3LYP functional. H2O forms an initial Lewis acid-base complex with the lanthanides exothermically with interaction energies from -2 to -20 kcal/mol. For most of the Ln, formation of HLnOH is more exothermic than formation of H2LnO, HLnO + H, and LnOH + H. The reactions to produce HLnOH are exothermic from -25 to -75 kcal/mol. The formation of LnO + H2 for La and Ce is slightly more exothermic than formation of HLnOH and is less or equally exothermic for the rest of the lanthanides. The Ln in HLnOH and LnOH are in the formal +II and +I oxidation states, respectively. The Ln in H2LnO is mostly in the +III formal oxidation state with either Ln-O(-)/Ln-H(-) or Ln-(H2)(-)/Ln=O(2-) bonding interactions. A few of the H2LnO have the Ln in the +IV or mixed +III/+IV formal oxidation states with Ln=O(2-)/Ln-H(-) bonding interactions. The Ln in HLnO are generally in the +III oxidation state with the exception of Yb in the +II state. The orbital populations calculated within the natural bond orbital (NBO) analysis are consistent with the oxidation states and reaction energies. The more exothermic reactions to produce HLnOH are always associated with more backbonding from the O(H) and H characterized by more population in the 6s and 5d in Ln and the formation of a stronger Ln-O(H) bond. Overall, the calculations are consistent with the experiments in terms of reaction energies and vibrational frequencies.

Extending the Row of Lanthanide Tetrafluorides: A Combined Matrix-Isolation and Quantum-Chemical Study.

Only the neutral tetrafluorides of Ce, Pr, and Tb as well as the [LnF7 ](3-) anions of Dy and Nd, with the metal in the +IV oxidation state, have been previously reported. We report our attempts to extend the row of neutral lanthanide tetrafluorides through the reaction of laser-ablated metal atoms with fluorine and their stabilization and characterization by matrix-isolation IR spectroscopy. In addition to the above three tetrafluorides, we found two new tetrafluorides, (3) NdF4 and (7) DyF4 , both of which are in the +IV oxidation state, which extends this lanthanide oxidation state to two new metals. Our experimental results are supported by quantum-chemical calculations and the role of the lanthanide oxidation state is discussed for both the LnF4 and [LnF4 ](-) species. Most of the LnF4 species are predicted to be in the +IV oxidation state and all of the [LnF4 ](-) anions are predicted to be in the +III oxidation state. The LnF4 species are predicted to be strong oxidizing agents and the LnF3 species are predicted to be moderate to strong Lewis acids.

1,2-Ethanediol and 1,3-Propanediol Conversions over (MO3)3 (M = Mo, W) Nanoclusters: A Computational Study

The dehydration and dehydrogenation reactions of one and two 1,2-ethanediol and 1,3-propanediol molecules on (MO3)3 (M = Mo, W) nanoclusters have been studied computationally using density functional and coupled cluster (CCSD(T)) theory. The reactions are initiated by the formation of a Lewis acid-base complex with an additional hydrogen bond. Dehydration is the dominant reaction proceeding via a metal bisdiolate. Acetaldehyde, the major product for 1,2-ethanediol, is produced by α-hydrogen transfer from one CH2 group to the other. For 1,3-propanediol, the C-C bond breaking pathways to produce C2H4 and HCH═O simultaneously and proton transfer to generate propylene oxide have comparable barrier energies. The barrier to produce propanal from the propylene oxide complex is less than that for epoxide release from the cluster. On the Mo3O9 cluster, a redox reaction channel for 1,2-ethanediol to break the C-C bond to form two formaldehyde molecules and then to produce C2H4 is slightly more favorable than the formation of acetaldehyde. For W(VI), the energy barrier for the reduction pathway is larger due to the lower reducibility of W3O9. Similar reduction on Mo(VI) for 1,3-propanediol to form propene is not a favorable pathway compared with the other pathways as additional C-H bond breaking is required in addition to breaking a C-C bond. The dehydrogenation and dehydration activation energies for the selected glycols are larger than the reactions of ethanol and 1-propanol on the same clusters. The CCSD(T) method is required because density functional theory with the M06 and B3LYP functionals does not predict quantitative energies on the potential energy surface. The M06 functional performs better than does the B3LYP functional.

Formation and Characterization of Homoleptic Thorium Isocyanide Complexes

Homoleptic thorium isocyanide complexes have been prepared via the reactions of laser-ablated thorium atoms and (CN)2 in a cryogenic matrix, and the structures of the products were characterized by infrared spectroscopy and theoretical calculations. Thorium atoms reacted with (CN)2 under UV irradiation to form the oxidative addition product Th(NC)2, which was calculated to have closed-shell singlet ground state with a bent geometry. Further reaction of Th(NC)2 and (CN)2 resulted in the formation of Th(NC)4, a molecule with a tetrahedral geometry. Minor products such as ThNC and Th(NC)3 were produced upon association reactions of CN with Th and Th(NC)2. Homoleptic thorium cyanide isomers Th(CN)x (x = 1-4) are predicted to be less stable than the corresponding isocyanides. The C-N stretches of thorium cyanides were calculated to be between 2170 and 2230 cm-1 at the B3LYP level, more than 120 cm-1 higher than the N-C stretches of isocyanides and with much weaker intensities. No experimental absorptions appeared where Th(CN)x should be observed.

Energetic Properties and Electronic Structure of [Si,N,S] and [Si,P,S] Isomers

Correlated molecular orbital theory at the coupled cluster CCSD(T) level with augmented correlation consistent basis sets has been used to predict the structure and energetic properties of the isomers of [Si,N,S] and [Si,P,S]. The predicted ground states are linear (2)SNSi and cyclic (2)SPSi. The other two isomers are predicted to be ∼20 to 50 kcal/mol less stable than the ground state. The excess spin is mainly on S for (2)SNSi and on P for (2)SPSi. The calculated total atomization energies with the CBS limits derived from different methods differ by ∼2 kcal/mol. The results provide the best available heats of formation for these species. The bond dissociation energies (BDEs) in (2)SNSi are comparable to those in the corresponding diatomic molecules. For cyclic (2)SPSi, the formation of (4)P + (2)SSi requires less energy than the other bond dissociation processes. The BDEs in the higher energy isomers are substantially smaller than the corresponding diatomic species.