Multilayer water adsorption on rutile TiO2(110): A first-principles study

Charge Behavior of Terminal Hydroxyl on Rutile TiO2(110)

Titanium dioxide (TiO₂) is of considerable interest as a photocatalyst and a catalyst support. Surface hydroxyl groups (OH) are the most common adsorbates on the TiO₂ surface and are believed to play crucial roles in their applications. Although the characteristics of bridging hydroxyl (OH_br) have been well understood, the adsorption structure and charged states of terminal hydroxyl (OH_t) have not yet been experimentally elucidated at an atomic scale. In this study, we have investigated an isolated OH_t on the rutile TiO₂(110) surface by atomic force microscopy (AFM) and Kelvin probe force microscopy (KPFM). We found that OH_t is in a negatively charged state. The unique characteristic of OH_t is different from that of OH_br and involves the amphoterism and diversity of catalytic reactions of TiO₂.

How the hydroxylation state of the (110)-rutile TiO2 surface governs its electric double layer properties

The hydroxylation state of an oxide surface is a central property of its solid/liquid interface and its corresponding electrical double layer. This study integrated both a reactive force field (ReaxFF) and a non-reactive potential into a hierarchical framework within molecular dynamics (MD) simulations to reveal how the hydroxylation state of the (110)-rutile TiO2 surface affects the electrical double layer properties. The simulation results obtained in the ReaxFF framework have shown that, while water dissociation occurs only at the under-coordinated Ti5c sites on the pristine TiO2 surface, the presence of point defects on the surface facilitates water dissociation at the oxygen vacancy sites, leading to two protonated oxygen bridge atoms for each vacancy site. As a consequence of enhanced water dissociation at the vacancy sites, water dissociation is quenched at the under-coordinated Ti5c sites resulting in two competitive hydroxylation mechanisms on the (110)-TiO2 surface. Using non-reactive MD simulations with hydroxylation states derived from the ReaxFF analysis, we demonstrate that water dissociation at the vacancy sites is a central mechanism governing the structuring of water near the interface. While the structuring of water near the interface is the main contribution to the electric field, water dissociation at the vacancy site enhances the adsorption of the electrolyte ions at the interface. The adsorbed ions lead to an increase of the effective surface charge as well as surface (zeta) potentials which are in the range of experimental observations. Our work provides a hierarchical multiscale simulation approach, covering a series of results with in-depth discussion for atomic/molecular level understanding of water dissociation and its effect on electric double layer properties of TiO2 to advance water splitting.

Unique Structures and Vibrational Spectra of Protic Ionic Liquids Confined in TiO2 Slits: The Role of Interfacial Hydrogen Bonds

The ionic liquid (IL)/titanium dioxide (TiO₂) interface exists in many application systems, such as nanomaterial synthesis, catalysis, and electrochemistry systems. The nanoscale interfacial properties in the above systems are a common issue. However, directly detecting the interfacial properties of nanoconfined ILs by experimental methods is still challenging. To help better learn about the interfacial issue, molecular dynamics simulations have been performed to explore the structures, vibration spectra, and hydrogen bond (HB) properties at the IL/TiO₂ interface. Ethylammonium nitrate (EAN) ILs confined in TiO₂ slit pores with different pore widths were studied. A unique vibrational spectrum appeared for EAN ILs confined in a 0.7 nm TiO₂ slit, and this phenomenon is related to interfacial hydrogen bonds (HBs). An analysis of the HB types indicated that the interfacial NH₃⁺ group of the cations was in an asymmetric HB environment in the 0.7 nm TiO₂ slit, which led to the disappearance of the symmetric N-H stretching mode. In addition, the significant increase in the HB strength between NH₃⁺ groups and the TiO₂ surface slowed down the stretching vibration of the N-H bond, resulting in one peak in the vibrational spectra at a lower frequency. For the first time, our simulation work establishes a molecular-level relationship between the vibrational spectrum and the local HB environment of nanoconfined ILs at the IL/TiO₂ interface, and this relationship is helpful for interface design in related systems.

Diffusion and reaction pathways of water near fully hydrated TiO2 surfaces from ab initio molecular dynamics

Ab initio molecular dynamics simulations are reported for water-embedded TiO₂ surfaces to determine the diffusive and reactive behavior at full hydration. A three-domain model is developed for six surfaces [rutile (110), (100), and (001), and anatase (101), (100), and (001)] which describes waters as "hard" (irreversibly bound to the surface), "soft" (with reduced mobility but orientation freedom near the surface), or "bulk." The model explains previous experimental data and provides a detailed picture of water diffusion near TiO₂ surfaces. Water reactivity is analyzed with a graph-theoretic approach that reveals a number of reaction pathways on TiO₂ which occur at full hydration, in addition to direct water splitting. Hydronium (H₃O⁺) is identified to be a key intermediate state, which facilitates water dissociation by proton hopping between intact and dissociated waters near the surfaces. These discoveries significantly improve the understanding of nanoscale water dynamics and reactivity at TiO₂ interfaces under ambient conditions.

Modeling of solid–liquid interfaces using scaled charges: rutile (110) surfaces

The first application of the electronic continuum correction model with scaled charges to molecular dynamics simulations of solid–liquid interfaces.

Ion Exchange Thermodynamics at the Rutile-Water Interface: Flow Microcalorimetric Measurements and Surface Complexation Modeling of Na-K-Rb-Cl-NO3 Adsorption

Flow microcalorimetry was used to investigate the energetics associated with Rb⁺, K⁺, Na⁺, Cl^-, and NO₃^- exchange at the rutile-water interface. Heats of exchange reflected differences in bulk hydration/dehydration enthalpies (Na⁺ > K⁺ > Rb⁺, and Cl^- > NO₃^-) such that exchanging Na⁺ or Cl^- from the surface was exothermic, reflecting their greater bulk hydration enthalpies. Exchange heats were measured at pH 2, 3.25, 5.8, and 11 and exhibited considerable differences as well as pH dependence. These trends were rationalized with the aid of a molecularly constrained surface complexation model (SCM) that incorporated the inner-sphere binding observed for the cations on the rutile (110) surface. Explicitly accounting for the inner-sphere binding configuration differences between Rb⁺, K⁺, and Na⁺, as well as accompanying differences in negative surface charge development, resulted in much better agreement with measured exchange ratios than by considering bulk hydration enthalpies alone. The observation that calculated exchange ratios agreed with those measured experimentally lends additional credence to the SCM. Consequently, flow microcalorimetry and surface complexation modeling are a useful complement of techniques for probing the energetics associated with ion exchange and adsorption processes and should also serve to help validate molecular simulations of interfacial energetics.

Modelling heterogeneous interfaces for solar water splitting

The generation of hydrogen from water and sunlight offers a promising approach for producing scalable and sustainable carbon-free energy. The key of a successful solar-to-fuel technology is the design of efficient, long-lasting and low-cost photoelectrochemical cells, which are responsible for absorbing sunlight and driving water splitting reactions. To this end, a detailed understanding and control of heterogeneous interfaces between photoabsorbers, electrolytes and catalysts present in photoelectrochemical cells is essential. Here we review recent progress and open challenges in predicting physicochemical properties of heterogeneous interfaces for solar water splitting applications using first-principles-based approaches, and highlights the key role of these calculations in interpreting increasingly complex experiments.

First Principles Study on the Interaction Mechanisms of Water Molecules on TiO₂ Nanotubes

The adsorption properties of water molecules on TiO₂ nanotubes (TiO₂NT) and the interaction mechanisms between water molecules are studied by first principles calculations. The adsorption preferences of water molecules in molecular or dissociated states on clean and H-terminated TiO₂NT are evaluated. Adsorption of OH clusters on (0, 6) and (9, 0) TiO₂ nanotubes are first studied. The smallest adsorption energies are -1.163 eV and -1.383 eV, respectively, by examining five different adsorption sites on each type of tube. Eight and six adsorption sites were considered for OH adsorbtion on the H terminated (0, 6) and (9, 0) nanotubes. Water molecules are reformed with the smallest adsorption energy of -4.796 eV on the former and of -5.013 eV on the latter nanotube, respectively. For the adsorption of a single water molecule on TiO₂NT, the molecular state shows the strongest adsorption preference with an adsorption energy of -0.660 eV. The adsorption of multiple (two and three) water molecules on TiO₂NT is also studied. The calculated results show that the interactions between water molecules greatly affect their adsorption properties. Competition occurs between the molecular and dissociated states. The electronic structures are calculated to clarify the interaction mechanisms between water molecules and TiO₂NT. The bonding interactions between H from water and oxygen from TiO₂NT may be the reason for the dissociation of water on TiO₂NT.

Reactive wetting properties of TiO2 nanoparticles predicted by ab initio molecular dynamics simulations

Small-sized wet TiO2 nanoparticles have been investigated by ab initio molecular dynamics simulations. Chemical and physical adsorption of water on the TiO2-water interface was studied as a function of water content, ranging from dry nanoparticles to wet nanoparticles with monolayer coverage of water. The surface reactivity was shown to be a concave function of water content and driven by surface defects. The local coordination number at the defect was identified as the key factor to decide whether water adsorption proceeds through dissociation or physisorption on the surface. A consistent picture of TiO2 nanoparticle wetting at the microscopic level emerges, which corroborates existing experimental data and gives further insight into the molecular mechanisms behind nanoparticle wetting. These calculations will facilitate the engineering of metal oxide nanoparticles with a controlled catalytic water activity.

Sulfur Dioxide Capture by Heterogeneous Oxidation on Hydroxylated Manganese Dioxide

Here we demonstrate that sulfur dioxide (SO2) is efficiently captured via heterogeneous oxidation into sulfate on the surface of hydroxylated manganese dioxide (MnO2). Lab-scale activity tests in a fluidized bed reactor showed that the removal efficiency for a simulated flue gas containing 5000 mg·Nm(-3) SO2 could reach nearly 100% with a GHSV (gas hourly space velocity) of 10000 h(-1). The mechanism was investigated using a combination of experimental characterizations and theoretical calculations. It was found that formation of surface bound sulfate proceeds via association of SO2 with terminal hydroxyls. Both H2O and O2 are essential for the generation of reactive terminal hydroxyls, and the indirect role of O2 in heterogeneous SO2 oxidation at low temperature was also revealed. We propose that the high reactivity of terminal hydroxyls is attributed to the proper surface configuration of MnO2 to adsorb water with degenerate energies for associative and dissociative states, and maintain rapid proton dynamics. Viability analyses suggest that the desulfurization method that is based on such a direct oxidation reaction at the gas/solid interface represents a promising approach for SO2 capture.

Facet-dependent solar ammonia synthesis of BiOCl nanosheets via a proton-assisted electron transfer pathway

Under the pressure of a fossil fuels shortage and global climate change, solar ammonia synthesis and the need to develop N2 fixation under mild conditions is becoming more urgent need; however, their intrinsic mechanisms still remain unclear. Herein, we demonstrate that the kinetic inertia of N2 can be overcome using oxygen vacancies (OVs) of BiOCl as the catalytic centers to create lower energy molecular steps, which are amendable for the solar light driven N-N triple bond cleavage via a proton-assisted electron transfer pathway. Moreover, the distinct structures of OVs on different BiOCl facets strongly determine the N2 fixation pathways by influencing both the adsorption structure and the activation level of N2. The fixation of terminal end-on bound N2 on the OVs of BiOCl {001} facets follows an asymmetric distal mode by selectively generating NH3, while the reduction of side-on bridging N2 on the OVs of BiOCl {010} facets is more energetically favorable in a symmetric alternating mode to produce N2H4 as the main intermediate.

Titania-catalyzed radiofluorination of tosylated precursors in highly aqueous medium

Nucleophilic radiofluorination is an efficient synthetic route to many positron-emission tomography (PET) probes, but removal of water to activate the cyclotron-produced [(18)F]fluoride has to be performed prior to reaction, which significantly increases overall radiolabeling time and causes radioactivity loss. In this report, we demonstrate the possibility of (18)F-radiofluorination in highly aqueous medium. The method utilizes titania nanoparticles, 1:1 (v/v) acetonitrile-thexyl alcohol solvent mixture, and tetra-n-butylammonium bicarbonate as a phase-transfer agent. Efficient radiolabeling is directly performed with aqueous [(18)F]fluoride without the need for a drying/azeotroping step to significantly reduce radiosynthesis time. High radiochemical purity of the target compound is also achieved. The substrate scope of the synthetic strategy is demonstrated with a range of aromatic, aliphatic, and cycloaliphatic tosylated precursors.

Aqueous solutions: state of the art in ab initio molecular dynamics

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.

Theoretical study of absorption of 2,2,6,6-tetramethylpiperidine-1-oxoammonium cation (TEMPO) on TiO₂(110) rutile surface

We present a theoretical study of the adsorption of 2,2,6,6-tetramethylpiperidine-1-oxoammonium cation (TEMPO) onto the TiO₂(110) surface rutile, investigating its bonding nature, electron properties and structural stability. Based on the results obtained with the PBE0/def2-SVP method, natural bond orbital (NBO) analysis suggests a bond order for the O--O bond in complexes 5 and 6, of 0.25 and 0.88, respectively. We also described NBOs for the main interactions of the TiO₂-TEMPO complexes.

A quantum mechanical study of water adsorption on the (110) surfaces of rutile SnO2 and TiO2: investigating the effects of intermolecular interactions using hybrid-exchange density functional theory

Periodic hybrid-exchange DFT is used to explore the structure and properties of SnO₂(110) and TiO₂(110) surfaces in contact with water.

Development of a ReaxFF reactive force field for titanium dioxide/water systems

A new ReaxFF reactive force field has been developed to describe reactions in the Ti-O-H system. The ReaxFF force field parameters have been fitted to a quantum mechanical (QM) training set containing structures and energies related to bond dissociation energies, angle and dihedral distortions, and reactions between water and titanium dioxide, as well as experimental crystal structures, heats of formation, and bulk modulus data. Model configurations for the training set were based on DFT calculations on molecular clusters and periodic systems (both bulk crystals and surfaces). ReaxFF reproduces accurately the QM training set for structures and energetics of small clusters. ReaxFF also describes the relative energetics for rutile, brookite, and anatase. The results of ReaxFF match reasonably well with those of QM for water binding energies, surface energies, and H2O dissociation energy barriers. To validate this ReaxFF description, we have compared its performance against DFT/MD simulations for 1 and 3 monolayers of water interacting with a rutile (110) surface. We found agreement within a 10% error between the DFT/MD and ReaxFF water dissociation levels for both coverages.

Adsorption of arginine-glycine-aspartate tripeptide onto negatively charged rutile (110) mediated by cations: the effect of surface hydroxylation

Classical molecular dynamics (MD) simulations were employed to investigate the adsorption behaviors of arginine-glycine-aspartate (RGD) tripeptide onto the negatively charged hydroxylated/nonhydroxylated rutile (110) surfaces, mediated by biologically important cations (Na+ or Ca2+). The simulation results indicate that the inherent nature of the cation determines its binding strength, thereby regulating the adsorption geometry of the peptide. The sparse hydroxyl groups on the nonhydroxylated rutile diminish the probability of H-bond formation between RGD and the surface, resulting in an early desorption of the peptide even with a mediating Na+ ion. In contrast, the negatively charged aspartate (Asp) side chain is bridged to the negatively charged hydroxylated rutile by an inner-sphere Na+ ion, in coordination with the Asp-rutile hydrogen bonds at the anchoring sites. The inner- and outer-sphere Ca2+ ions are demonstrated to be capable of "trapping" RGD on both hydroxylated and nonhydroxylated rutile, in the absence of hydrogen bonds with the surface. The strongly bound inner-sphere mediating Ca2+ ion exerts a "gluing" effect on the Asp side chain, producing a tightly packed RGD-rutile complex, whereas a less localized distribution density of the outer-sphere mediating Ca2+ ion results in a higher mobility of the Asp side chain. The intramolecular interaction is suggested to facilitate the structural stability of RGD adsorbed on the negative rutile in a "horseshoe" configuration.

Physical origin of general oscillation of structure, surface energy, and electronic property in rutile TiO2 nanoslab

Titanium oxide (TiO(2)) nanostructures have been attracting consistent focus in the past few years because of their enhanced power in solar-energy conversion. Surface and interface play a crucial role in the determination of thermodynamic stability and electronic structure of TiO(2) nanostructures. The rutile (110) nanoslab (NS) has been used as a common subject to investigate the surface relaxation, defect characters, molecule adsorption, and chemically dynamic reaction of TiO(2) nanostructures. Up to date, a long-time standing issue in TiO(2) NS, i.e., the general oscillation of structure, surface energy and electronic property with changing of NS thickness, has not been clear. We have presented a comprehensive investigation on the relationship between surface and oscillation behavior in the TiO(2) (110) NS by the first-principles calculations accompanied with the wave function analysis. We clearly, for the first time, pointed out that the dipoles and surface states bonding induced by the surface-surface interactions are the physical origin of general oscillations in the TiO(2) (110) NS. Our findings not only have a new insight into the basic interactions between surfaces in TiO(2) nanostructures, but also provide useful information for tuning the photocatalytic and photovoltaic properties by surface design.

Polarization- and Azimuth-Resolved Infrared Spectroscopy of Water on TiO2(110): Anisotropy and the Hydrogen-Bonding Network

We have investigated the structure and dynamics of thin water films adsorbed on TiO2(110) using infrared reflection-absorption spectroscopy (IRAS) and ab initio molecular dynamics. Infrared spectra were obtained for s- and p-polarized light with the plane of incidence parallel to the [001] and [11̅0] azimuths for water coverages ≤ 4 monolayers. The spectra indicate strong anisotropy in the water films. The vibrational densities of states predicted by the ab initio simulations for 1 and 2 monolayer coverages agree well with the observations. The results provide new insight into the structure of water on TiO2(110) and resolve a long-standing puzzle regarding the hydrogen bonding between molecules in the first and second monolayers on this surface. The results also demonstrate the capabilities of polarization- and azimuth-resolved IRAS for investigating the structure and dynamics of adsorbates on dielectric substrates.

Determining the orientation and molecular packing of organic dyes on a TiO2 surface using X-ray reflectometry

The determination of the orientation and molecular density for several porphyrin dyes adsorbed on planar TiO(2) surfaces using X-ray reflectometry (XRR) is reported. Adsorption of nanoscale water layers occurred rapidly upon exposure of freshly prepared TiO(2) surfaces to ambient conditions; however, this was successfully eliminated, resulting in clearly discernible adsorbed dye layers for sensitized surfaces. Adsorbed dye orientations, determined from computations constrained by the measured dye layer thickness, were calculated to have a binding tilt angle of 35°-40°. Combining the XXR data with the orientation models indicates that the porphyrins form densely packed surfaces with an intermolecular spacing of 3-4 Å, consistent with π-π stacking interactions. Changes in the molecular size of probe dyes were reflected in corresponding changes in the measured dye layer thickness, confirming the ability of this technique to resolve small variations in dye layer thickness and consequently adsorption orientation. Application of these results to understanding the behavior of dye-sensitized devices is discussed.

Mapping the surface (hydr)oxo-groups of titanium oxide and its interface with an aqueous solution: the state of the art and a new approach

In this article the "titanium oxide/electrolyte solution" interface is studied by taking in advantage the recent developments in the field of Surface and Interface Chemistry relevant to this oxide. Ab-initio calculations were performed in the frame of the DFT theory for estimating the charge of the titanium and oxygen atoms exposed on the anatase (1 0 1), (1 0 0), (0 0 1), (1 0 3)(f) and rutile (1 1 0) crystal faces. These orientations have smaller surface energy with respect to other ones and thus it is more probable to be the real terminations of the anatase and rutile nanocrystallites in the titania polycrystalline powders. Potentiometric titrations for obtaining "fine structured" titration curves as well as microelectrophoresis and streaming potential measurements have been performed. On the basis of ab-initio calculations, and taking into account the relative contribution of each crystal face to the whole surface of the nanocrystals involved in the titania aggregates of a suspension, the three most probable surface ionization models have been derived. These models and the Music model are then tested in conjunction with the "Stern-Gouy-Chapman" and "Basic Stern" electrostatic models. The finally selected surface ionization model (model A) in combination with each one of the two electrostatic models describes very well the protonation/deprotonation behavior of titania. The description is also very good if this model is combined with the Three Plane (TP) model. The application of the "A/(TP)" model allowed mapping the surface (hydr)oxo-groups [TiO(H) and Ti(2)O(H)] of titania exposed in aqueous solutions. At pH>pzc almost all terminal oxygens [TiO] are non-protonated whereas even at low pH values the non-protonated terminal oxygens predominate. The acid-base behavior of the bridging oxygens [Ti(2)O] is different. Thus, even at pH=10 the greater portion of them is protonated. The application of the "A/TP" model in conjunction with potentiometric titrations, microelectrophoresis and streaming potential experiments allowed mapping the "titania/electrolyte solution" interface. It was found that the first (second) charged plane is located on the oxygen atoms of the first (second) water overlayer at a distance of 1.7 (3.4) A from the surface. The region between the surface and the second plane is the compact layer. The region between the second plane and the shear plane is the stagnant diffuse part of the interface, with an ionic strength dependent width, ranging from 20 (0.01 M) up to 4 A (0.3 M). The region between the shear plane and the bulk solution is the mobile diffuse part, with an ionic strength dependent width, ranging from 10 (0.01 M) up to 2 A (0.3 M). At I>0.017 M the mean concentration of the counter ions is higher in the stagnant than in the mobile part of the diffuse layer. For a given I, removal of pH from pzc brings about an increase of the mean concentration in the interfacial region and a displacement of the counter ions from the mobile to the stagnant part of the diffuse layer. The mean concentration of the counter ions in the compact layer is generally lower than the corresponding ones in the stagnant and mobile diffuse layers. The mobility of the counter ions in the stagnant layer decreases as pH draws away from pzc or ionic strength increases.

Hydrogen bonding, H-D exchange, and molecular mobility in thin water films on TiO2(110)

Hydrogen bonding, H-D exchange, and molecular mixing in water films on TiO2(110) have been studied using water electron-stimulated desorption. For T<70 K, films with different water isotopes adsorbed on the Ti4+ and bridging oxygen rows can be prepared. For T>70 K, H-D and molecular exchange between these "layers" occur with a distribution of activation energies and is surprisingly efficient. The results demonstrate that all the water molecules that are directly bound to TiO2(110) -- i.e., in the first and second monolayers -- are also hydrogen bonded to each other.

Reaction intermediates in the photoreduction of oxygen molecules at the (101) TiO2 (anatase) surface

The structural, electronic, and vibrational properties of intermediates of the O(2) photoreduction at the (101) TiO(2) (anatase) surface have been investigated by performing ab initio density functional calculations. In detail, a recently proposed approach has been used where molecules on the surface are treated like surface defects. Thus, by applying theoretical methods generally used in the physics of semiconductors, we successfully estimate the location and donor/acceptor character of the electronic levels induced by an adsorbed molecule in the TiO(2) energy gap, both crucial for the surface-molecule charge-transfer processes, and investigate the formation and the properties of charged intermediates. The present approach permits a view of the O(2) photoreduction process through several facets, which elucidates the molecule-surface charge-transfer conditions and reveals the key role played by charged intermediates. A comparison of present results with those of a highly sensitive IR (infrared) spectroscopy study of intermediates of the O(2) photoreduction leads to a deeper understanding of this process and to revised vibrational-line assignments and reaction paths.

First-principles calculations of hydrogen diffusion on rutile TiO2(110) surfaces

Density functional calculations are performed to study the H-atom diffusion on titanium dioxide (110) surface in the cases of water-molecule dissociation and splitting of the adjacent hydroxyl OH pair. It is shown that, when a water molecule is adsorbed at a surface oxygen-vacancy site, a fragment H atom of the water molecule tends to diffuse toward the nearest-neighboring bridging-oxygen sites by using a straight-line or relay-point path. As the result, a pair of surface hydroxyl OH is formed on the same oxygen row. In a thermal process, on the other hand, such OH pair favorably splits only by using a relay-point path, i.e., by transferring one H atom from a bridging-oxygen site to a next-neighboring one along the same oxygen row by way of another in-plane oxygen site. We found that the latter splitting reaction is activated around room temperature.

Computer simulations of surfactant monolayers at solid walls

A series of molecular dynamics simulations to study the structure of a surfactant monolayer near real surfaces was carried out. A comparison of two different surfaces, TiO2 and SiO2, with the same monolayer was performed. Moreover, each surface was modeled by two different approaches, the first model considers the complete structure of a TiO2 (or SiO2) wall, whereas the second model is a continuous solid wall with an effective potential. Both wall models give essentially the same monolayer configuration suggesting that the explicit form of the wall does not play a relevant role to study the structure of surfactant monolayers close to planar surfaces.

Adsorption bond length for H2O on TiO2(110): a key parameter for theoretical understanding

Scanned-energy mode photoelectron diffraction results show the adsorption site of molecular water on TiO2(110) to be atop under-coordinated surface Ti atoms, confirming the results of total energy calculations and STM imaging. However, the Ti-O(water) bond length is 2.21 +/- 0.02 A, much longer than Ti-O bond lengths in strongly chemisorbed species on this surface, but significantly shorter than found in most total energy calculations. The need for theory to describe this weak bond effectively may be a key factor in the controversial problem of understanding this important surface reaction system.

Reactions and clustering of water with silica surface

The interaction between silica surface and water is an important topic in geophysics and materials science, yet little is known about the reaction process. In this study we use first-principles molecular dynamics to simulate the hydrolysis process of silica surface using large cluster models. We find that a single water molecule is stable near the surface but can easily dissociate at three-coordinated silicon atom defect sites in the presence of other water molecules. These extra molecules provide a mechanism for hydrogen transfer from the original water molecule, hence catalyzing the reaction. The two-coordinated silicon atom is inert to the water molecule, and water clusters up to pentamer could be stably adsorbed at this site at room temperature.

Vertical and lateral order in adsorbed water layers on anatase TiO2(101)

The structure and energetics of thin water overlayers on the (101) surface of TiO(2)-anatase have been studied through first-principles molecular dynamics simulations at T = 160 K. At one monolayer coverage, H(2)O molecules are adsorbed at the 5-fold Ti sites (Ti(5c)), forming an ordered crystal-like 2D layer with no significant water-water interactions. For an adsorbed bilayer, H(2)O molecules at both Ti(5c) and bridging oxygen (O(2c)) sites form a partially ordered structure, where the water oxygens occupy regular sites but the orientation of the molecules is disordered; in addition, stress-relieving defects are usually present. When a third layer is adsorbed, very limited parallel and perpendicular order is observed above the first bilayer. The calculated energetics of multilayer adsorption is in good agreement with recent temperature-programmed desorption data.

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.

A density functional theory study of the coadsorption of water and oxygen on TiO2(110)

The behavior of adsorbed water on oxides is of fundamental interest in many areas. Despite considerable attention received recently, our understanding of water chemistry is still short of needs and expectations, particularly on the topic of the coadsorption of water and other species. In this study we carry out density functional theory calculations to investigate the coadsorption of water and oxygen on the TiO(2)(110) surface. We show that oxygen exerts profound influences on the water adsorption, altering the mechanism of water dissociation. On the one hand, the possible dissociation route along [-110] is prohibited due to the weakening of the H bond between water and the lattice bridging oxygen in the presence of the coadsorbed oxygen, and on the other hand the coadsorbed oxygen induces dissociation along [001]. These results lead to a consistent interpretation of experiments. Furthermore, several possible final states and the related formation mechanisms are discussed in detail.