Probing the low-barrier hydrogen bond in hydrogen maleate in the gas phase: a photoelectron spectroscopy and ab initio study

The strength of the low-barrier hydrogen bond in hydrogen maleate in the gas phase was investigated by low-temperature photoelectron spectroscopy and ab initio calculations. Photoelectron spectra of maleic and fumaric acid monoanions (cis-/trans-HO(2)CCH=CHCO(2)(-)) were obtained at low temperatures and at 193 nm photon energy. Vibrational structure was observed for trans-HO(2)CCH=CHCO(2)(-) due to the OCO bending modes; however, cis-HO(2)CCH=CHCO(2)(-) yielded a broad and featureless spectrum. The electron binding energy of cis-HO(2)CCH=CHCO(2)(-) is about 1 eV blue-shifted relative to trans-HO(2)CCH=CHCO(2)(-) due to the formation of intramolecular hydrogen bond in the cis-isomer. Theoretical calculations (CCSD(T)/ aug-cc-pVTZ and B3LYP/aug-cc-pVTZ) were carried out to estimate the strength of the intramolecular hydrogen bond in cis-HO(2)CCH=CHCO(2)(-). Combining experimental and theoretical calculations yields an estimate of 21.5 +/- 2.0 kcal/mol for the intramolecular hydrogen bond strength in hydrogen maleate.     

On the stability of the elusive HO3 radical

The dissociation of HO(3) into OH + O(2) has been studied in a systematic and consistent way using the multireference configuration interaction method. Upon extrapolation of the calculated raw energies to the complete basis set limit and using jointly with a recent realistic estimate of the zero-point vibrational energy, the energy for OO-OH bond-breaking in the trans isomer is predicted to be of D(0) = (2.4 ± 0.1) kcal mol(-1), where the uncertainty reflects only the one inherent to the extrapolation. The average value so obtained falls short of the commonly accepted experimental counterpart by 0.5 kcal mol(-1). Reasons for the deviation are advanced, as well as an estimate of the binding energy for the cis-HO(3) isomer which is predicted to have a somewhat smaller binding energy than trans-HO(3), but likewise the latter dissociates without a barrier to the same products.     

Is HO3 minimum cis or trans? An analytic full-dimensional ab initio isomerization path

The minimum energy path for isomerization of HO(3) has been explored in detail using accurate high-level ab initio methods and techniques for extrapolation to the complete basis set limit. In agreement with other reports, the best estimates from both valence-only and all-electron single-reference methods here utilized predict the minimum of the cis-HO(3) isomer to be deeper than the trans-HO(3) one. They also show that the energy varies by less than 1 kcal mol(-1) or so over the full isomerization path. A similar result is found from valence-only multireference configuration interaction calculations with the size-extensive Davidson correction and a correlation consistent triple-zeta basis, which predict the energy difference between the two isomers to be of only Δ = -0.1 kcal mol(-1). However, single-point multireference calculations carried out at the optimum triple-zeta geometry with basis sets of the correlation consistent family but cardinal numbers up to X = 6 lead upon a dual-level extrapolation to the complete basis set limit of Δ = (0.12 ± 0.05) kcal mol(-1). In turn, extrapolations with the all-electron single-reference coupled-cluster method including the perturbative triples correction yield values of Δ = -0.19 and -0.03 kcal mol(-1) when done from triple-quadruple and quadruple-quintuple zeta pairs with two basis sets of increasing quality, namely cc-cpVXZ and aug-cc-pVXZ. Yet, if added a value of 0.25 kcal mol(-1) that accounts for the effect of triple and perturbative quadruple excitations with the VTZ basis set, one obtains a coupled cluster estimate of Δ = (0.14 ± 0.08) kcal mol(-1). It is then shown for the first time from systematic ab initio calculations that the trans-HO(3) isomer is more stable than the cis one, in agreement with the available experimental evidence. Inclusion of the best reported zero-point energy difference (0.382 kcal mol(-1)) from multireference configuration interaction calculations enhances further the relative stability to ΔE(ZPE) = (0.51 ± 0.08) kcal mol(-1). A scheme is also suggested to model the full-dimensional isomerization potential-energy surface using a quadratic expansion that is parametrically represented by a Fourier analysis in the torsion angle. The method illustrated at the raw and complete basis-set limit coupled-cluster levels can provide a valuable tool for a future analysis of the available (incomplete thus far) experimental rovibrational data.     

Ab initio and analytic intermolecular potentials for Ar-CH(3)OH

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.     

Theoretical characterization of the gas-phase O(3)HO hydrogen-bonded complex.

We report a theoretical study on two gas-phase hydrogen-bonded complexes formed between ozone and hydroxyl radical that have relevance to atmospheric chemistry. This study was carried out by using CASSCF, CASPT2, QCISD, and CCSD(T) theoretical approaches in conjunction with the 6-311+G(2df,2p) and aug-cc-pVTZ basis sets. Both complexes have a planar structure and differ from each other in the orientation of the electronic density of the unpaired electron associated with the HO radical moiety. Our calculations predict their stabilities to be 0.87 and 0.67 kcal mol(-1), respectively, at 0 K and show the importance of anharmonic effects in computing the red shift of the HO stretch originating from the hydrogen-bonding interaction. We also report two transition states involving the movement of the HO moiety on the potential energy surfaces of these hydrogen-bonded complexes.     

On the nature of the unusually long OO bond in HO(3) and HO(4) radicals

The HO(3) and HO(4) polyoxide radicals have attracted some attention due to their potential role in ozone chemistry. Experimentally, the geometrical structure of HO(3) is known whereas that of HO(4) is not. Moreover, the existence of the latter radical has been questioned. Theoretical calculations on the two species have been reported before, showing important structural differences depending on the computational level. Both radicals present an unusually long OO bond (around 1.7-1.8 A) that can be associated with an intricate interaction between HO, or HO(2), with O(2). The nature of such interaction is investigated in detail using large scale ab initio methods (CASSCF, CASPT2, MRCI, QCISD, CCSD(T)) and density functional techniques (B3LYP) in connection with extended basis sets. Stabilization enthalpies at 298 K with respect to HO (or HO(2)) and O(2) have been computed amounting to -3.21 kcal mol(-1) for HO(3) (trans conformation) and 11.33 kcal mol(-1) for HO(4) (cis conformation). The corresponding formation enthalpies are 6.12 and 11.83 kcal mol(-1). The trans conformation of HO(4) is less stable than the cis one by 6.17 kcal mol(-1). Transition states for HO(4) dissociation and for cis/trans conversion are also described.     

CO oxidation catalyzed by a single gold atom: benchmark calculations and the performance of DFT methods

Quantum chemical calculations were carried out on CO oxidation catalyzed by a single gold atom. To investigate the performance of density functional theory (DFT) methods, 42 DFT functionals have been evaluated and compared with high-level wavefunction based methods. It was found that in order to obtain accurate results the functionals used must treat long range interaction well. The double-hybrid mPW2PLYP and B2PLYP functionals are the two functionals with best overall performance. CAM-B3LYP, a long range corrected hybrid GGA functional, also performs well. On the other hand, the popular B3LYP, PW91, and PBE functionals do not show good performance and the performance of the latter two are even at the bottom of the 42 functionals. Our accurate results calculated at the CCSD(T)/aug-cc-pVTZ//mPW2PLYP/aug-cc-pVTZ level of theory indicate that Au atom is a good catalysis for CO oxidation. The reaction follows the following mechanism where CO and O(2) adsorb on Au atom forming an Au(OCOO) intermediate and subsequently O(2) transfer one oxygen atom to CO to form CO(2) and AuO. Then AuO reacts with CO to form another CO(2) to complete the catalytic cycle. The overall energy barrier at 0 K for the first CO oxidation step (Au + CO + O(2)→ AuO + CO(2)) is just 4.8 kcal mol(-1), and that for the second CO oxidation step (AuO + CO → Au + CO(2)) is just 1.6 kcal mol(-1).     

A theoretical study of surface reduction mechanisms of CeO(2)(111) and (110) by H(2).

Reaction mechanisms for the interactions between CeO(2)(111) and (110) surfaces are investigated using periodic density functional theory (DFT) calculations. Both standard DFT and DFT+U calculations to examine the effect of the localization of Ce 4f states on the redox chemistry of H(2)-CeO(2) interactions are described. For mechanistic studies, molecular and dissociative local minima are initially located by placing an H(2) molecule at various active sites of the CeO(2) surfaces. The binding energies of physisorbed species optimized using the DFT and DFT+U methods are very weak. The dissociative adsorption reactions producing hydroxylated surfaces are all exothermic; exothermicities at the DFT level range from 4.1 kcal mol(-1) for the (111) to 26.5 kcal mol(-1) for the (110) surface, while those at the DFT+U level are between 65.0 kcal mol(-1) for the (111) and 81.8 kcal mol(-1) for the (110) surface. Predicted vibrational frequencies of adsorbed OH and H(2)O species on the surfaces are in line with available experimental and theoretical results. Potential energy profiles are constructed by connecting molecularly adsorbed and dissociatively adsorbed intermediates on each CeO(2) surface with tight transition states using the nudged elastic band (NEB) method. It is found that the U correction method plays a significant role in energetics, especially for the intermediates of the exit channels and products that are partially reduced. The surface reduction reaction on CeO(2)(110) is energetically much more favorable. Accordingly, oxygen vacancies are more easily formed on the (110) surface than on the (111) surface.     

Determining the energy gap between the cis and trans isomers of HO3- using geometry optimization within the anti-Hermitian contracted Schrödinger and coupled cluster methods

The cis and trans isomers of the HO3- anion, which are important in proposed mechanisms for ozonization, are studied computationally. Relative energies, geometries, and normal-mode frequencies are calculated with anti-Hermitian contracted Schr?dinger equation (ACSE) and coupled cluster methods. Both the ACSE method and the coupled cluster method with single and double excitations (CCSD) are applied in a correlation-consistent polarized double-zeta basis set (cc-pVDZ). Using coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), we treat the problem with larger basis sets than those in previous work, including correlation-consistent polarized quadruple-zeta basis sets with (aug-cc-pVQZ) and without (cc-pVQZ) diffuse functions, which permit extrapolation of the cis and trans energies to the complete-basis-set limit. The cis isomer is found to be lower in energy than the trans isomer by -3.5 kcal/mol, which is 50% larger in magnitude than the best previous result of -2.2 kcal/mol. The bond lengths between the O2 and OH fragments of the cis- and trans-HO3 are calculated to be 1.713 and 1.857 A, respectively, where both bond lengths are significantly longer than the 1.464 A O-O bond in hydrogen peroxide. In this paper, we extend the ACSE method [Mazziotti, D. A. J. Chem. Phys. 2007, 126, 184101], which computes the two-electron reduced density matrix directly, to include geometry optimization by a Newton's method with numerical derivatives. Calculation of the cis- and trans-HO3- isomers by the ACSE yields energies, geometries, and frequencies that are closer to those from CCSD(T) than those from CCSD.     

Carbene Rearrangements Unsurpassed: Details of the C(7)H(6) Potential Energy Surface Revealed

The rearrangement of phenylcarbene (1) to 1,2,4,6-cycloheptatetraene (3) has been studied theoretically, using SCF, CASSCF, CASPT2N, DFT (B3LYP), CISD, CCSD, and CCSD(T) methods in conjunction with the 6-31G, 6-311+G, 6-311G(2d,p), cc-pVDZ, and DZd basis sets. Stationary points were characterized by vibrational frequency analyses at CASSCF(8,8)/6-31G and B3LYP/6-31G. Phenylcarbene (1) has a triplet ground state ((3)A") with a singlet-triplet separation (DeltaE(ST)) of 3-5 kcal mol(-)(1). In agreement with experiment, chiral 3 is the lowest lying structure on this part of the C(7)H(6) potential energy surface. Bicyclo[4.1.0]hepta-2,4,6-triene (2) is an intermediate in the rearrangement of 1 into 3, but it is unlikely to be observable experimentally due to a barrier height of only 1-2 kcal mol(-)(1). The enantiomers of 3 interconvert via the (1)A(2) state of cycloheptatrienylidene (4) with an activation energy of 20 kcal mol(-)(1). The "aromatic" (1)A(1) state, previously believed to be the lowest singlet state of 4, is roughly 10 kcal mol(-)(1) higher in energy than the (1)A(2) state, and, in violation of Hund's rule, (3)A(2) is also calculated to lie above (1)A(2) in energy. Thus, even if (3)A(2) were populated, it is likely to undergo rapid intersystem crossing to (1)A(2). We suggest (3)B(1)-4 is the metastable triplet observed by EPR.     

Benchmark calculations on the adiabatic ionization potentials of M-NH(3) (M=Na,Al,Ga,In,Cu,Ag)

The ground states of the M-NH(3) (M=Na,Al,Ga,In,Cu,Ag) complexes and their cations have been studied with density functional theory and coupled cluster [CCSD(T)] methods. The adiabatic ionization potentials (AIPs) of these complexes are calculated, and these are compared to results from high-resolution zero-electron kinetic energy photoelectron spectroscopy. By extrapolating the CCSD(T) energies to the complete basis set (CBS) limit and including the core-valence, scalar relativistic, spin-orbit, and zero-point corrections, the CCSD(T) method is shown to be able to predict the AIPs of these complexes to better than 6 meV or 0.15 kcal/mol. 27 exchange-correlation functionals, including one in the local density approximation, 13 in the generalized gradient approximation (GGA), and 13 with hybrid GGAs, were benchmarked in the calculations of the AIPs. The B1B95, mPW1PW91, B98, B97-1, PBE1PBE, O3LYP, TPSSh, and HCTH93 functionals give an average error of 0.1 eV for all the complexes studied, with the B98 functional alone yielding a maximum error of 0.1 eV. In addition, the calculated metal-ammonia harmonic stretching frequencies with the CCSD(T) method are in excellent agreement with their experimental values, whereas the B3LYP method tends to underestimate these stretching frequencies. The metal-ammonia binding energies were also calculated at the CCSD(T)/CBS level, and are in excellent agreement with the available experimental values considering the error limits, except for Ag-NH(3) and Ag(+)-NH(3), where the calculations predict stronger bond energies than measured by about 4 kcal/mol, just outside the experimental error bars of +/-3 kcal/mol.     

Theoretical study of CCl(4) adsorption and hydrogenation on a Pt (111) surface

The adsorption and hydrogenation of carbon tetrachloride (CCl(4)) on a Pt (111) surface have been investigated using density functional theory (DFT). We have performed calculations on the adsorption energies and structures of CCl(4) on four different adsorption sites of a Pt (111) surface using the full adsorbate geometry optimization method. The results show that the adsorption energy of all of the potential sites is less than -17 kcal/mol, which indicates that CCl(4) is physiosorbed on a Pt (111) surface through van der Waals interactions. The dissociation and hydrogenation pathways were investigated by a transition state search. For the Pt(15), Pt(19), and Pt(25) cluster surfaces, the activation energies of dissociation obtained in this work are 15.69, 16.94, and 16.77 kcal/mol, respectively. The hydrogenation of CCl(3). was studied at the on-top site of the Pt(15) cluster, and the calculated activation energy is 5.06 kcal/mol. The small activation energies indicate that the Pt (111) surface has high catalytic activity for the CCl(4) hydrogenation reaction. In addition, the Hirshfeld population analysis reveals that the charge transfer from the Pt (111) surface to the adsorbates occurs in both the dissociation and hydrogenation pathways.     

Computational study of the release of H2 from ammonia borane dimer (BH3NH3)2 and its ion pair isomers

High-level electronic structure calculations have been used to map out the relevant portions of the potential energy surfaces for the release of H2 from dimers of ammonia borane, BH3NH3 (AB). Using the correlation-consistent aug-cc-pVTZ basis set at the second-order perturbation MP2 level, geometries of stationary points were optimized. Relative energies were computed at these points using coupled-cluster CCSD(T) theory with the correlation-consistent basis sets at least up to the aug-cc-pVTZ level and in some cases extrapolated to the complete basis set limit. The results show that there are a number of possible dimers involving different types of hydrogen-bonded interactions. The most stable gaseous phase (AB)2 dimer results from a head-to-tail cyclic conformation and is stabilized by 14.0 kcal/mol with respect to two AB monomers. (AB)2 can generate one or two H2 molecules via several direct pathways with energy barriers ranging from 44 to 50 kcal/mol. The diammoniate of diborane ion pair isomer, [BH4-][NH3BH2NH3+] (DADB), is 10.6 kcal/mol less stable than (AB)2 and can be formed from two AB monomers by overcoming an energy barrier of approximately 26 kcal/mol. DADB can also be generated from successive additions of two NH3 molecules to B2H6 and from condensation of AB with separated BH3 and NH3 molecules. The pathway for H2 elimination from DADB is characterized by a smaller energy barrier of 20.1 kcal/mol. The alternative ion pair [NH4+][BH3NH2BH3-] is calculated to be 16.4 kcal/mol above (AB)2 and undergoes H2 release with an energy barrier of 17.7 kcal/mol. H2 elimination from both ion pair isomers yields the chain BH3NH2BH2NH3 as product. Our results suggest that the neutral dimer will play a minor role in the release of H2 from ammonia borane, with a dominant role from the ion pairs as observed experimentally in ionic liquids and the solid state.     

Periodic density functional theory study of methane activation over La(2)O(3): activity of O(2-), O(-), O(2)(2-), oxygen point defect,and Sr(2+)-doped surface sites

Results of gradient-corrected periodic density functional theory calculations are reported for hydrogen abstraction from methane at O(s)(2-), O(s)(-), O(2)(s)(2-) point defect, and Sr(2+)-doped surface sites on La(2)O(3)(001). The results show that the anionic O(s)(-) species is the most active surface oxygen site. The overall reaction energy to activate methane at an O(s)(-) site to form a surface hydroxyl group and gas-phase (*)CH(3) radical is 8.2 kcal/mol, with an activation barrier of 10.1 kcal/mol. The binding energy of hydrogen at an site O(s)(-) is -102 kcal/mol. An oxygen site with similar activity can be generated by doping strontium into the oxide by a direct Sr(2+)/La(3+) exchange at the surface. The O(-)-like nature of the surface site is reflected in a calculated hydrogen binding energy of -109.7 kcal/mol. Calculations indicate that surface peroxide (O(2(s))(2-)) sites can be generated by adsorption of O(2) at surface oxygen vacancies, as well as by dissociative adsorption of O(2) across the closed-shell oxide surface of La(2)O(3)(001). The overall reaction energy and apparent activation barrier for the latter pathway are calculated to be only 12.1 and 33.0 kcal/mol, respectively. Irrespective of the route to peroxide formation, the O(2)(s)(2-) intermediate is characterized by a bent orientation with respect to the surface and an O-O bond length of 1.47 A; both attributes are consistent with structural features characteristic of classical peroxides. We found surface peroxide sites to be slightly less favorable for H-abstraction from methane than the O(s)(-) species, with DeltaE(rxn)(CH(4)) = 39.3 kcal/mol, E(act) = 47.3 kcal/mol, and DeltaE(ads)(H) = -71.5 kcal/mol. A possible mechanism for oxidative coupling of methane over La(2)O(3)(001) involving surface peroxides as the active oxygen source is suggested.     

Thermodynamics and kinetics of oxygen-induced segregation of 3d metals in Pt-3d-Pt(111) and Pt-3d-Pt(100) bimetallic structures

The stability of subsurface 3d transition metals (3d represents Ni, Co, Fe, Mn, Cr, V, and Ti) in Pt(111) and Pt(100) was examined in vacuum and with 0.5 ML atomic oxygen by a combined experimental and density functional theory (DFT) approach. DFT was used to predict the trends in the binding energy of oxygen and in the stability of 3d metals to remain in the subsurface layer. DFT calculations predicted that for both (111) and (100) crystal planes the subsurface Pt-3d-Pt configurations were thermodynamically preferred in vacuum and that the surface 3d-Pt-Pt configurations were preferred with the adsorption of 0.5 ML atomic oxygen. Experimentally, the DFT predictions were verified by using Auger electron spectroscopy to monitor the segregation of Ni and Co in Pt-3d-Pt structures on polycrystalline Pt foil, composed of mainly (111) and (100) facets. The activation barrier for the oxygen-induced segregation of Ni was found to be 17+/-1 kcal/mol attributed to the Pt(111) areas and 27+/-1 kcal/mol attributed to the Pt(100) areas of the Pt foil. For Pt-Co-Pt, the activation barrier was found to be 10+/-1 kcal/mol and was attributed to the Pt(111) areas of the Pt foil. The Bronsted-Evans-Polanyi relationship was utilized to predict the activation barriers for segregation of the other Pt-3d-Pt(111) and Pt-3d-Pt(100) systems. These results are further discussed in connection to the activity and stability for cathode bimetallic electrocatalysts for proton exchange membrane fuel cells.     

Reactions of hydrazoic acid on TiO2 nanoparticles: an experimental and computational study

This article reports the results of a computational and experimental study on the reaction of hydrazoic acid, HN3, adsorbed on 15-20 nm TiO2 particle films. Experimentally, FTIR spectra of HN3(a) have been measured by varying HN3 dosage, UV irradiation time and surface annealing temperature. Three sharp peaks, related to v(a)(NNN) of HN3(a) and N3(a) with different configurations in the 2000-2200 cm(-1) region, and a broad band absorption, related to associated and isolated HN(a) and HO(a) adsorptions in the 3000-3800 cm(-1) region, have been detected. Computationally, molecular structures, vibrational frequencies and adsorption energies of possible adsorbates including HN3 and its derivatives, N3, N2, NH, and H, have been predicted by first-principles calculations based on the density functional theory (DFT) and the pseudopotential method. On the basis of the experimental and computational results, the peak appeared at 2075 cm(-1), which increases at a faster rate with HN3 exposure time, is attributed to a stable adsorbate, N3-Ti(a), with the predicted adsorption energy, E(ads) = 13 kcal/mol. The peak at 2118 cm(-1), which survives at the highest surface temperature in the heating experiment, is attributable to the most stable adsorbate, Ti-N2N(H)-O(a) with E(ads) = 36 kcal/mol. The peak at 2170 cm(-1), which vanishes most readily in all of the aforementioned experiments, is related to less stable molecular adsorbates, end-on HN3-Ti(a) with E(ads) = 5 kcal/mol and side-on HN(N2)-Ti(a) with E(ads) = 8 kcal/mol. A potential energy diagram for the formation of various absorbates with their transition states has been established for the HN3/TiO2 system. On the basis of the predicted desorption energies, the four most stable products of the HN3 reaction on TiO2 are H-O(a), 118 kcal/mol; HN-O(a), 85 kcal/mol; Ti-N2N(H)-O(a), 36 kcal/mol; and N3-O(a), 19 kcal/mol.     

The self-assembly mechanism of the Lindqvist anion [W(6)O(19)](2-) in aqueous solution: a density functional theory study

The formation mechanism is always a fundamental and confused issue for polyoxometalate chemistry. Two formation mechanisms (M1 and M2) of the Lindqvist anion [W(6)O(19)](2-) have been adopted to investigate it's self-assembly reaction pathways at a density functional theory (DFT) level. The potential energy surfaces reveal that both the mechanisms are thermodynamically favorable and overall barrierless at room temperature, but M2 is slightly dominant to M1. The formation of the pentanuclear species [W(5)O(16)](2-) and [W(5)O(15)(OH)](-) are recognized as the rate-determining steps in the whole assembly polymerization processes. These two steps involve the highest energy barriers with 30.48 kcal mol(-1) and 28.90 kcal mol(-1), respectively, for M1 and M2. [W(4)O(13)](2-) and [W(4)O(12)(OH)](-) are proved to be the most stable building blocks. In addition, DFT results reveal that the formation of [W(3)O(10)](2-) experiences a lower barrier along the chain channel.     

Computational study on the reactions of H2O2 on TiO2 anatase (101) and rutile (110) surfaces

This study investigates the adsorption and reactions of H(2)O(2) on TiO(2) anatase (101) and rutile (110) surfaces by first-principles calculations based on the density functional theory in conjunction with the projected augmented wave approach, using PW91, PBE, and revPBE functionals. Adsorption mechanisms of H(2)O(2) and its fragments on both surfaces are analyzed. It is found that H(2)O(2) , H(2)O, and HO preferentially adsorb at the Ti(5c) site, meanwhile HOO, O, and H preferentially adsorb at the (O(2c))(Ti(5c)), (Ti(5c))(2), and O(2c) sites, respectively. Potential energy profiles of the adsorption processes on both surfaces have been constructed using the nudged elastic band method. The two restructured surfaces, the 1/3 ML oxygen covered TiO(2) and the hydroxylated TiO(2), are produced with the H(2)O(2) dehydration and deoxidation, respectively. The formation of main products, H(2)O(g) and the 1/3 ML oxygen covered TiO(2) surface, is exothermic by 2.8 and 5.0 kcal/mol, requiring energy barriers of 0.8 and 1.1 kcal/mol on the rutile (110) and anatase (101) surface, respectively. The rate constants for the H(2)O(2) dehydration processes have been predicted to be 6.65 × 10(-27) T(4.38) exp(-0.14 kcal mol(-1)/RT) and 3.18 × 10(-23) T(5.60) exp(-2.92 kcal mol(-1)/RT) respectively, in units of cm(3) molecule(-1) s(-1).     

Adequate representation of charge polarization effects leads to a successful treatment of the CF4 + SiCl4 → CCl4 + SiF4 reaction by density functional theory

Adequate polarization functions reduce the error of density functional theory (DFT) for the heat of reaction for CF(4) + SiCl(4) from ～9-12 kcal mol(-1) to ～2-4 kcal mol(-1), and using an improved density functional further reduces it to ～1 kcal mol(-1). This reaction was previously identified as a stumbling block for DFT, but we show that the problem with the previous calculations was not DFT but rather inadequate basis sets to account for intramolecular charge polarization.