Transient hydroxyl formation from water on oxygen-covered Au(111)

We present evidence for the formation of transient hydroxyls from the reaction of water with atomic oxygen on Au(111) and investigate the effect of adsorbed oxygen on the hydrogen bonding of water. Water is evolved in peaks at 175 and 195 K in temperature programed reaction experiments following adsorption of water on oxygen-covered Au(111). The peak at 175 K is ascribed to sublimation of multilayers of water, whereas the peak at 195 K is associated with oxygen-stabilized water or a water-hydroxyl surface complex. Infrared reflection absorption spectra are consistent with the presence of molecular water over the entire range of coverages studied, indicating that isolated stable hydroxyls are not formed. Isotopic exchange of adsorbed (16)O with H(2)(18)O following adsorption and subsequent temperature programed reaction, however, indicates that transient OH species are formed. The extent of oxygen exchange was considerable--up to 70%. The degree of oxygen exchange depends on the initial coverage of oxygen, the surface temperature when preparing oxygen adatoms, and the H(2)(18)O coverage. The hydroxyls are short-lived, forming and disproportionating multiple times before water desorption during temperature programed reaction. It was also found that chemisorbed oxygen is critical in the formation of hydroxyls and stabilizing water, whereas gold oxide does not contribute to these effects. These results identify transient hydroxyls as species that could play a critical role in oxidative chemical reactions on gold, especially in ambient water vapor. The crystallinity of adsorbed water also depended on the degree of surface ordering and chemical modification based on scanning tunneling microscopy and infrared spectra. These results demonstrate that oxidation of interfaces has a major impact on their interaction with water.     

Adsorption of water on reconstructed rutile TiO2(011)-(2 x 1): Ti=O double bonds and surface reactivity

Recent combined experimental and theoretical studies (Beck et al., Phys. Rev. Lett. 2004, 93, 036104) have provided evidence for Ti=O double-bonded titanyl groups on the reconstructed rutile TiO(2)(011)-(2 x 1) surface. The adsorption of water on the same surface is now investigated to further probe the properties of these groups, as well as to confirm their existence. Ultraviolet photoemission experiments show that water is adsorbed in molecular form at a sample temperature of 110 K. At the same time, the presence of a 3sigma state in the photoemission spectra and work function measurements indicate a significant amount of hydroxyls within the first monolayer of water. At room temperature, scanning tunneling microscopy (STM) suggests that dissociated water is present, and about 30% of the surface active sites are hydroxylated. These findings are well explained by total energy density functional theory calculations and Car-Parrinello molecular dynamics simulations for water adsorption on the titanyl model of TiO(2)(011)-(2 x 1). The theoretical results show that a mixed molecular/dissociative layer is the most stable configuration in the monolayer regime at low temperatures, while complete dissociation takes place at 250 K. The arrangement of the protonated mono-coordinated oxygens in the mixed molecular/dissociated layer is consistent with the observed short-range order of the hydroxyls in the STM images.     

Water-enhanced low-temperature CO oxidation and isotope effects on atomic oxygen-covered Au(111)

Water-oxygen interactions and CO oxidation by water on the oxygen-precovered Au(111) surface were studied by using molecular beam scattering techniques, temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. Water thermally desorbs from the clean Au(111) surface with a peak temperature of approximately 155 K; however, on a surface with preadsorbed atomic oxygen, a second water desorption peak appears at approximately 175 K. DFT calculations suggest that hydroxyl formation and recombination are responsible for this higher temperature desorption feature. TPD spectra support this interpretation by showing oxygen scrambling between water and adsorbed oxygen adatoms upon heating the surface. In further support of these experimental findings, DFT calculations indicate rapid diffusion of surface hydroxyl groups at temperatures as low as 75 K. Regarding the oxidation of carbon monoxide, if a C (16)O beam impinges on a Au(111) surface covered with both atomic oxygen ( (16)O) and isotopically labeled water (H 2 (18)O), both C (16)O (16)O and C (16)O (18)O are produced, even at surface temperatures as low as 77 K. Similar experiments performed by impinging a C (16)O beam on a Au(111) surface covered with isotopic oxygen ( (18)O) and deuterated water (D 2 (16)O) also produce both C (16)O (16)O and C (16)O (18)O but less than that produced by using (16)O and H 2 (18)O. These results unambiguously show the direct involvement and promoting role of water in CO oxidation on oxygen-covered Au(111) at low temperatures. On the basis of our experimental results and DFT calculations, we propose that water dissociates to form hydroxyls (OH and OD), and these hydroxyls react with CO to produce CO 2. Differences in water-oxygen interactions and oxygen scrambling were observed between (18)O/H 2 (16)O and (18)O/D 2 (16)O, the latter producing less scrambling. Similar differences were also observed in water reactivity toward CO oxidation, in which less CO 2 was produced with (16)O/D 2 (16)O than with (16)O/H 2 (16)O. These differences are likely due to primary kinetic isotope effects due to the differences in O-H and O-D bond energies.     

Coverage dependence and hydroperoxyl-mediated pathway of catalytic water formation on Pt (111) surface

Hydrogen oxidation on Pt (111) surface is modeled by density functional theory (DFT). Previous DFT calculations showed too large O2 dissociation barriers, but we find them highly coverage dependent: when the coverage is low, dissociation barriers close to experimental values (approximately 0.3 eV) are obtained. For the whole reaction, a new pathway involving hydroperoxyl (OOH) intermediate is found, with the highest reaction barrier of only approximately 0.4 eV. This may explain the experimental observation of catalytic water formation on Pt (111) surface above the H2O desorption temperature of 170 K, despite that the direct reaction between chemisorbed O and H atoms is a highly activated process with barrier approximately 1 eV as previous calculations showed.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Structural requirements and reaction pathways in dimethyl ether combustion catalyzed by supported Pt clusters

The identity and reversibility of the elementary steps required for catalytic combustion of dimethyl ether (DME) on Pt clusters were determined by combining isotopic and kinetic analyses with density functional theory estimates of reaction energies and activation barriers to probe the lowest energy paths. Reaction rates are limited by C-H bond activation in DME molecules adsorbed on surfaces of Pt clusters containing chemisorbed oxygen atoms at near-saturation coverages. Reaction energies and activation barriers for C-H bond activation in DME to form methoxymethyl and hydroxyl surface intermediates show that this step is more favorable than the activation of C-O bonds to form two methoxides, consistent with measured rates and kinetic isotope effects. This kinetic preference is driven by the greater stability of the CH3OCH2* and OH* intermediates relative to chemisorbed methoxides. Experimental activation barriers on Pt clusters agree with density functional theory (DFT)-derived barriers on oxygen-covered Pt(111). Measured DME turnover rates increased with increasing DME pressure, but decreased as the O2 pressure increased, because vacancies (*) on Pt surfaces nearly saturated with chemisorbed oxygen are required for DME chemisorption. DFT calculations show that although these surface vacancies are required, higher oxygen coverages lead to lower C-H activation barriers, because the basicity of oxygen adatoms increases with coverage and they become more effective in hydrogen abstraction from DME. Water inhibits reaction rates via quasi-equilibrated adsorption on vacancy sites, consistent with DFT results indicating that water binds more strongly than DME on vacancies. These conclusions are consistent with the measured kinetic response of combustion rates to DME, O2, and H2O, with H/D kinetic isotope effects, and with the absence of isotopic scrambling in reactants containing isotopic mixtures of 18O2-16O2 or 12CH3O12CH3-13CH3O13CH3. Turnover rates increased with Pt cluster size, because small clusters, with more coordinatively unsaturated surface atoms, bind oxygen atoms more strongly than larger clusters and exhibit lower steady-state vacancy concentrations and a consequently smaller number of adsorbed DME intermediates involved in kinetically relevant steps. These effects of cluster size and metal-oxygen bond energies on reactivity are ubiquitous in oxidation reactions requiring vacancies on surfaces nearly saturated with intermediates derived from O2.     

Evidence for molecularly chemisorbed oxygen on TiO2 supported gold nanoclusters and Au(111)

In this study, we present evidence for the existence of a molecularly chemisorbed oxygen species on a Au/TiO2 model catalyst and a Au(111) single crystal following exposure of these samples to an oxygen plasma-jet molecular beam. We present evidence for the molecularly chemisorbed oxygen species from thermal desorption, collision-induced desorption, and heat of adsorption/reaction-induced desorption measurements. Thermal desorption measurements reveal a peak desorption temperature at approximately 145 K which corresponds to an activation energy for desorption of approximately 0.35 eV.     

Computational characterization of reaction intermediates in the photocycle of the sensory domain of the AppA blue light photoreceptor

The AppA protein with the BLUF (blue light using flavin adenine dinucleotide) domain is a blue light photoreceptor that cycle between dark-adapted and light-induced functional states. We characterized possible reaction intermediates in the photocycle of AppA BLUF. Molecular dynamics (MD), quantum chemical and quantum mechanical-molecular mechanical (QM/MM) calculations were carried out to describe several stable structures of a molecular system modeling the protein. The coordinates of heavy atoms from the crystal structure (PDB code 2IYG) of the protein in the dark state served as starting point for 10 ns MD simulations. Representative MD frames were used in QM(B3LYP/cc-pVDZ)/MM(AMBER) calculations to locate minimum energy configurations of the model system. Vertical electronic excitation energies were estimated for the molecular clusters comprising the quantum subsystems of the QM/MM optimized structures using the SOS-CIS(D) quantum chemistry method. Computational results support the occurrence of photoreaction intermediates that are characterized by spectral absorption bands between those of the dark and light states. They agree with crystal structures of reaction intermediates (PDB code 2IYI) observed in the AppA BLUF domain. Transformations of the Gln63 side chain stimulated by photo-excitation and performed with the assistance of the chromophore and the Met106 side chain are responsible for these intermediates.     

二氧化碳和丙烯直接合成甲基丙烯酸的NiPMo/TiO2催化剂的研究

用溶胶-凝胶法以磷钼酸(MPA)的镍盐溶液水解钛酸四丁酯制备了NiPMo/TiO2催化剂.使用ICP、 XRD、 TG-DTA、 IR、 TPD-MS和微反应技术研究了催化剂的化学组成、热稳定性、化学吸附性质和催化反应性能.杂多钼酸盐与TiO2通过O2-在TiO2表面发生了键合.在623 K下,杂多阴离子仍保持原有的Keggin结构.CO2在Lewis酸位Ni(Ⅱ)和Lewis碱位Ni-O-Mo的桥氧协同作用下生成CO2卧式吸附态Ni(Ⅱ)←O-(CO)←(O--Ni).丙烯有多种吸附态在催化剂上吸附.在563 K、 1 MPa和空速1500 h-1的反应条件下,丙烯的摩尔转化率为3.2%,产物MAA选择性为95%.     

Multilayer adsorption of water at a rutile TiO2)(110) surface: towards a realistic modeling by molecular dynamics

We present a model combining ab initio concepts and molecular dynamics simulations for a more realistic treatment of complex adsorption processes. The energy, distance, and orientation of water molecules adsorbed on stoichiometric and reduced rutile TiO(2)(110) surfaces at 140 K are studied via constant temperature molecular dynamics simulations. From ab initio calculations relaxed atomic geometries for the surface and the most probable adsorption sites were derived. The study comprises (i) large two-dimensional surface supercells, providing a realistically low concentration of surface oxygen defects, and (ii) a water coverage sufficiently large to model the onset of the growth of a bulk phase of water on the surface. By our combined approach the influence of both, the metal oxide surface, below, and the bulk water phase, above, on the water molecules forming the interface between the TiO(2) surface and the water bulk layer is taken into account. The good agreement of calculated adsorption energies with experimental temperature programmed desorption spectra demonstrates the validity of our model.     

CO2氧化丙烷脱氢MoO3-V2O5/SiO2催化剂研究

采用表面改性和等体积浸渍法制备了MoO3-V2O5／SiO2复合氧化物，采用TPR、IR、TPD和微反技术表征了催化剂对CO2和C3H8的吸附性能和CO2部分氧化丙烷脱氢的反应性能．结果表明：随MoO3加入量的增加，V=O还原温度降低，晶格氧更易活化．CO2和丙烷的活化温度降低，丙烷转化率增高，丙烯选择性下降．V=O中晶格氧的活化是CO2氧化丙烷制丙烯反应的关键．     

Acetone and water on TiO2(110): competition for sites

The competitive interaction between acetone and water for surface sites on TiO2(110) was examined using temperature programmed desorption (TPD). Two surface pretreatment methods were employed, one involving vacuum reduction of the surface by annealing at 850 K in ultrahigh vacuum (UHV) and another involving surface oxidation with molecular oxygen. In the former case, the surface possessed about 7% oxygen vacancy sites, and in the latter, reactive oxygen species (adatoms and molecules) were deposited on the surface as a result of oxidative filling of vacancy sites. On the 7% oxygen vacancy surface, excess water displaced all but about 20% of a saturated d6-acetone first layer to physisorbed desorption states, whereas about 40% of the first layer d6-acetone was stabilized on the oxidized surface against displacement by water through a reaction between oxygen and d6-acetone. The displacement of acetone on both surfaces is explained in terms of the relative desorption energies of each molecule on the clean surface and the role of intermolecular repulsions in shifting the respective desorption features to lower temperatures with increasing coverage. Although first layer water desorbs from TiO2(110) at slightly lower temperature (275 K) than submonolayer coverages of d6-acetone (340 K), intermolecular repulsions between d6-acetone molecules shift its leading edge for desorption to 170 K as the first layer is saturated. In contrast, the desorption leading edge for first layer water (with or without coadsorbed d6-acetone) shifted to no lower than 210 K as a function of increasing coverage. This small difference in the onsets for d6-acetone and water desorption resulted in the majority of d6-acetone being compressed into islands by water and displaced from the first layer at a lower temperature than that observed in the absence of coadsorbed water. On the oxidized surface, the species resulting from reaction of d6-acetone and oxygen was not influence by increasing water coverages. This species was stable up to 375 K (well past the first layer water TPD feature) where it decomposed mostly back to d6-acetone and atomic oxygen. These results are discussed in terms of the influence of water in inhibiting acetone photo-oxidation on TiO2 surfaces.     

On the mechanism of low-temperature water gas shift reaction on copper

Periodic, self-consistent density functional theory (DFT-GGA) calculations are used to investigate the water gas shift reaction (WGSR) mechanism on Cu(111). The thermochemistry and activation energy barriers for all the elementary steps of the commonly accepted redox mechanism, involving complete water activation to atomic oxygen, are presented. Through our calculations, we identify carboxyl, a new reactive intermediate, which plays a central role in WGSR on Cu(111). The thermochemistry and activation energy barriers of the elementary steps of a new reaction path, involving carboxyl, are studied. A detailed DFT-based microkinetic model of experimental reaction rates, accounting for both the previous and the new WGSR mechanism show that, under relevant experimental conditions, (1) the carboxyl-mediated route is the dominant path, and (2) the initial hydrogen abstraction from water is the rate-limiting step. Formate is a stable "spectator" species, formed predominantly through CO2 hydrogenation. In addition, the microkinetic model allows for predictions of (i) surface coverage of intermediates, (ii) WGSR apparent activation energy, and (iii) reaction orders with respect to CO, H2O, CO2, and H2.     

Theoretical study of the reaction of carbon monoxide with oxygen molecules in the ground triplet and singlet delta states

Quantum chemical calculations were carried out to study the reaction of carbon monoxide with molecular oxygen in the ground triplet and singlet delta states. Transition states and intermediates that connect the reactants with products of the reaction on the triplet and singlet potential energy surfaces were identified on the base of coupled-cluster method. The values of energy barriers were refined by using compound techniques such as CBS-Q, CBS-QB3, and G3. The calculations showed that there exists an intersection of triplet and singlet potential energy surfaces. This fact leads to the appearance of two channels for the triplet CO+O(2)(X(3)Σ(g)(-)) reaction, which produces atomic oxygen in the ground O((3)P) and excited O((1)D) states. The appropriate rate constants of all reaction paths were estimated on the base of nonvariational transition-state theory. It was found that the singlet reaction rate constant is much greater than the triplet one and that the reaction channel CO+O(2)(a(1)Δ(g)) should be taken into consideration to interpret the experimental data on the oxidation of CO by molecular oxygen.     

Primary intermediates of oxygen photoevolution reaction on TiO2 (Rutile) particles, revealed by in situ FTIR absorption and photoluminescence measurements

Primary intermediates of oxygen photoevolution (water photooxidation) reaction at the TiO2 (rutile)/aqueous solution interface were investigated by in situ multiple internal reflection infrared (MIRIR) absorption and photoluminescence (PL) measurements. UV irradiation of TiO2 in the presence of 10 mM Fe3+ in the solution caused the appearance of a new peak at 838 cm(-1) and a shoulder at 812 cm(-1). Detailed investigations of the effects of solution pH, the presence of methanol as a hole scavenger, and isotope exchange in water (H2(16)O-->H2(18)O) on the spectra have shown that the 838- and 812-cm(-1) bands can be assigned to the O-O stretching mode of surface TiOOH and TiOOTi, respectively, produced as primary intermediates of the oxygen photoevolution reaction. The results give strong support to our previously proposed mechanism that the oxygen photoevolution is initiated by a nucleophilic attack of a H2O molecule on a photogenerated hole at a surface lattice O site, not by oxidation of surface OH group by the hole. The conclusion is supported by PL measurements. A plausible reaction scheme is proposed for the oxygen photoevolution on TiO2 (rutile) in aqueous solutions of pH less than about 12.     

Ab initio molecular dynamics simulations of the oxygen reduction reaction on a Pt(111) surface in the presence of hydrated hydronium (H3O)(+)(H2O)2: direct or series pathway?

Car-Parrinello molecular dynamics simulations have been performed to investigate the oxygen reduction reaction (ORR) on a Pt(111) surface at 350 K. By progressive loading of (H3O)(+)(H2O)(2,3) + e- into a simulation cell containing a Pt slab and O2 for the first reduction step, and either products or intermediate species for the subsequent reduction steps, the detailed mechanisms of the ORR are well illustrated via monitoring MD trajectories and analyzing Kohn-Sham electronic energies. A proton transfer is found to be involved in the first reduction step; depending on the initial proton-oxygen distance, on the degree of proton hydration, and on the surface charge, such transfer may take place either earlier or later than the O2 chemisorption, in all cases forming an adsorbed end-on complex H-O-O*. Decomposition of H-O-O* takes place with a rather small barrier, after a short lifetime of approximately 0.15 ps, yielding coadsorbed oxygen and hydroxyl (O + HO*). Formation of the one-end adsorbed hydrogen peroxide, HOO*H, is observed via the reduction of H-O-O*, which suggests that the ORR may also proceed via HOO*H, i.e., a series pathway. However, HOO*H readily dissociates homolytically into two coadsorbed hydroxyls (HO* + HO*) rather than forming a dual adsorbed HOOH. Along the direct pathway, the reduction of H-O* + O* yields two possible products, O* + H2O* and HO* + HO*. Of the three intermediates from the second electron-transfer step, HOO*H from the series pathway has the highest energy, followed by O* + H2O* and HO* + HO* from the direct pathway. It is therefore theoretically validated that the O2 reduction on a Pt surface may proceed via a parallel pathway, the direct and series occurring simultaneously, with the direct as the dominant step.     

Modeling water exchange on an aluminum polyoxocation

For the first time, water exchange on a polymeric complex has been modeled using a combination of gas-phase ab initio calculations and molecular dynamics (MD) simulations. The GaO4Al12(OH)24(H2O)12(7+)aq ion (GaAl12) was chosen because high-quality experimental data exist, including an activation enthalpy (+63 +/- 7 kJ/mol) and an activation volume (+3 +/- 1 cm3/mol). We took a two-step approach. First, the local solvent structure and the initial states for reaction were inferred from the molecular dynamics simulations. Second, we used this information to evaluate initial-state structures in the ab initio calculations. The energy differences between the initial and transition states from the ab initio calculations varied from +59 kJ/mol to +53 kJ/mol depending upon details, closely approximating the activation enthalpy.     

甲醛在CeO2(111)表面吸附的密度泛函理论研究

采用基于第一性原理的密度泛函理论和周期平板模型, 研究了甲醛在以桥氧为端面的CeO2(111)稳定表面上的吸附行为. 通过对不同覆盖度, 不同吸附位的甲醛吸附构型、吸附能及电子态密度的分析发现, 甲醛在CeO2(111)表面存在化学吸附与物理吸附两种情况. 化学吸附结构中甲醛的碳、氧原子分别与表面的氧、铈原子发生相互作用, 形成CH2O2物种; 吸附能随着覆盖度的增加而减小. 与自由甲醛分子相比, 物理吸附的甲醛构型变化不大, 其吸附能较小. 利用CNEB(climbing nudged elastic band)方法计算了甲醛在CeO2(111)表面的初步解离反应活化能(约1.71 eV), 远高于甲醛脱附能垒, 这与甲醛在清洁CeO2(111)表面程序升温脱附实验中产物主要为甲醛的结果相一致.     

Oxygen vacancy promoting catalytic dehydration of formic acid on TiO2(110) by in situ scanning tunneling microscopic observation

The catalytic dehydration reaction processes of formic acid on a TiO2(110) surface at 350 K have been studied to visualize reaction intermediates and their dynamic behaviors by scanning tunneling microscopy. Three types of configurations of adsorbed formates on the surface were identified by their shapes and positions in STM images. Successive STM observations revealed transformations among the three configurations, i.e., bridge formate on a 5-fold coordinated Ti4+ row, bridge formate on an oxygen vacancy site with an oxygen atom of formate and on a 5-fold coordinated Ti4+ ion and with the other formate oxygen atom, and a monodentate formate on an oxygen vacancy site with an oxygen atom of formate. The decomposition of the monodentate formate to carbon monoxide and hydroxyl was also imaged, which is a rate-determining step in the catalytic dehydration of formic acid. Combined with first-principle DFT calculations, the overall reaction processes of the catalytic dehydration of formic acid on the surface have been elucidated. Oxygen vacancies on the surface that can be produced by dehydration of two hydroxyls in situ under the catalytic reaction conditions are essential for the reaction.     

Interaction between water molecules and zinc sulfide nanoparticles studied by temperature-programmed desorption and molecular dynamics simulations

We have investigated the bonding of water molecules to the surfaces of ZnS nanoparticles (approximately 2-3 nm sphalerite) using temperature-programmed desorption (TPD). The activation energy for water desorption was derived as a function of the surface coverage through kinetic modeling of the experimental TPD curves. The binding energy of water equals the activation energy of desorption if it is assumed that the activation energy for adsorption is nearly zero. Molecular dynamics (MD) simulations of water adsorption on 3 and 5 nm sphalerite nanoparticles provided insights into the adsorption process and water binding at the atomic level. Water binds with the ZnS nanoparticle surface mainly via formation of Zn-O bonds. As compared with bulk ZnS crystals, ZnS nanoparticles can adsorb more water molecules per unit surface area due to the greatly increased curvature, which increases the distance between adjacent adsorbed molecules. Results from both TPD and MD show that the water binding energy increases with decreasing the water surface coverage. We attribute the increase in binding energy with decreasing surface water coverage to the increasing degree of surface under-coordination as removal of water molecules proceeds. MD also suggests that the water binding energy increases with decreasing particle size due to the further distance and hence lower interaction between adsorbed water molecules on highly curved smaller particle surfaces. Results also show that the binding energy, and thus the strength of interaction of water, is highest in isolated nanoparticles, lower in nanoparticle aggregates, and lowest in bulk crystals. Given that water binding is driven by surface energy reduction, we attribute the decreased binding energy for aggregated as compared to isolated particles to the decrease in surface energy that occurs as the result of inter-particle interactions.