CO2 adsorption on TiO2(110) rutile: insight from dispersion-corrected density functional theory calculations and scanning tunneling microscopy experiments

Adsorption of CO(2) on the rutile(110) surface was investigated using dispersion-corrected density functional theory and scanning tunneling microscopy (STM). On the oxidized surface the CO(2) molecules are found to bind most strongly at the five-fold coordinated Ti sites adopting tilted or flat configurations. The presence of bridging oxygen defects introduces two new adsorption structures, the most stable of which involves CO(2) molecules bound in tilted configurations at the defect sites. Inclusion of dispersion corrections in the density functional theory calculations leads to large increases in the calculated adsorption energies bringing these quantities into good agreement with experimental data. The STM measurements confirm two of the calculated adsorption configurations.     

First principles theoretical study of the hole-assisted conversion of CO to CO2 on the anatase TiO2(101) surface

First principles density functional theory calculations were carried out to investigate the adsorption and oxidation of CO on the positively charged (101) surface of anatase, as well as the desorption of CO(2) from it. We find that the energy gain on adsorption covers the activation energy required for the oxidation, while the energy gain on the latter is sufficient for the desorption of CO(2), leaving an oxygen vacancy behind. Molecular dynamics simulations indicate that the process can be spontaneous at room temperature. The oxidation process described here happens only in the presence of the hole. The possibility of a photocatalytic cycle is discussed assuming electron scavenging by oxygen.     

NO adsorption on MoS(x) clusters: a density functional theory study

The density functional theory (DFT) method has been used to investigate NO probe molecule adsorption on the stoichiometric (Mo(16)S(32)) and nonstoichiometric (Mo(16)S(34) and Mo(16)S(29)) clusters. The calculated adsorption energies indicate that the stoichiometric cluster has stronger NO affinity than the nonstoichiometric surfaces. It is also found that mononitrosyl adsorption is favored at low NO coverage, while dinitrosyl (germinal) and (NO)(2) dimer adsorption at high NO coverage are possible. Strong repulsive interaction has been found for the adsorbed dinitrosyl and (NO)(2) dimer species. In addition, the computed NO stretching frequencies for the mononitrosyl and dinitrosyl species agree well with the experimental data, while those of the dimer species are much lower than the suggested experimental data.     

Supramolecular binding thermodynamics by dispersion-corrected density functional theory

The equilibrium association free enthalpies ΔG(a) for typical supramolecular complexes in solution are calculated by ab initio quantum chemical methods. Ten neutral and three positively charged complexes with experimental ΔG(a) values in the range 0 to -21?kcal?mol(-1) (on average -6?kcal?mol(-1) ) are investigated. The theoretical approach employs a (nondynamic) single-structure model, but computes the various energy terms accurately without any special empirical adjustments. Dispersion corrected density functional theory (DFT-D3) with extended basis sets (triple-ζ and quadruple-ζ quality) is used to determine structures and gas-phase interaction energies (ΔE), the COSMO-RS continuum solvation model (based on DFT data) provides solvation free enthalpies and the remaining ro-vibrational enthalpic/entropic contributions are obtained from harmonic frequency calculations. Low-lying vibrational modes are treated by a free-rotor approximation. The accurate account of London dispersion interactions is mandatory with contributions in the range -5 to -60?kcal?mol(-1) (up to 200?% of ΔE). Inclusion of three-body dispersion effects improves the results considerably. A semilocal (TPSS) and a hybrid density functional (PW6B95) have been tested. Although the ΔG(a) values result as a sum of individually large terms with opposite sign (ΔE vs. solvation and entropy change), the approach provides unprecedented accuracy for ΔG(a) values with errors of only 2?kcal?mol(-1) on average. Relative affinities for different guests inside the same host are always obtained correctly. The procedure is suggested as a predictive tool in supramolecular chemistry and can be applied routinely to semirigid systems with 300-400 atoms. The various contributions to binding and enthalpy-entropy compensations are discussed.     

Density functional study of the interfacial electron transfer pathway for monolayer-adsorbed InN on the TiO(2) anatase (101) surface

Density functional theory (DFT) in connection with ultrasoft pseudopotential (USP) and generalized gradient spin-polarized approximations (GGSA) is applied to calculate the adsorption energies and structures of monolayer-adsorbed InN on the TiO(2) anatase (101) surface and the corresponding electronic properties, that is, partial density of states (PDOS) for surface and bulk layers of the TiO(2) anatase (101) surface and monolayer-adsorbed InN, to shed light on the possible structural modes for initial photoexcitation within the UV/vis adsorption region followed by fast electron injection through the InN/TiO(2) interface for an InN/TiO(2)-based solar cell design. Our calculated adsorption energies found that the two most probable stable structural modes of monolayer-adsorbed InN on the TiO(2) anatase (101) surface are (1) an end-on structure with an adsorption energy of 2.52 eV through N binding to surface 2-fold coordinated O (O(cn2)), that is, InN-O(cn2), and (2) a side-on structure with an adsorption energy of 3.05 eV through both N binding to surface 5-fold coordinated Ti (Ti(cn5)) and In bridging two surface O(cn2), that is, (O(cn2))(2)-InN-Ti(cn5). Our calculated band gaps for both InN-O(cn2) and (O(cn2))2-InN-Ti(cn5) (including a 1.0-eV correction using a scissor operator) of monolayer-adsorbed InN on the TiO(2) anatase (101) surface are red-shifted to 1.7 eV (730 nm) and 2.3 eV (540 nm), respectively, which are within the UV/vis adsorption region similar to Gratzel's black dye solar cell. Our analyses of calculated PDOS for both surface and bulk layers of the TiO(2) anatase (101) surface and monolayer-adsorbed InN on the TiO(2) anatase (101) surface suggest that the (O(cn2))(2)-InN-Ti(n5) configuration of monolayer-adsorbed InN on the TiO(2) anatase (101) surface would provide a more feasible structural mode for the electron injection through the InN/TiO(2) interface. This is due to the presence of both occupied and unoccupied electronic states for monolayer-adsorbed InN within the band gap TiO(2) anatase (101) surface, which will allow the photoexcitation within the UV/vis adsorption region to take place effectively, and subsequently the photoexcited electronic states will overlap with the unoccupied electronic states around the lowest conduction band of the TiO(2) anatase (101) surface, which will ensure the electron injection through the InN/TiO(2) interface. Finally, another thing worth our attention is our preliminary study of double-layer-adsorbed InN on the TiO(2) anatase (101) surface, that is, (O(cn2))(2)-(InN)(2)-Ti(cn5), with a calculated band gap red-shifted to 2.6 eV (477 nm) and a different overlap of electronic states between double-layer-adsorbed InN and the TiO(2) anatase (101) surface qualitatively indicated that there is an effect of the thickness of adsorbed InN on the TiO(2) anatase (101) surface on both photoexcitation and electron injection processes involved in the photoinduced interfacial electron transfer through InN/TiO(2). A more thorough and comprehensive study of different layers of InN adsorbed in all possible different orientations on the TiO(2) anatase (101) surface is presently in progress.     

Coadsorption properties of CO(2) and H(2)O on TiO(2) rutile (110): A dispersion-corrected DFT study

Adsorption and reactions of CO(2) in the presence of H(2)O and OH species on the TiO(2) rutile (110)-(1×1) surface were investigated using dispersion-corrected density functional theory and scanning tunneling microscopy. The coadsorbed H(2)O (OH) species slightly increase the CO(2) adsorption energies, primarily through formation of hydrogen bonds, and create new binding configurations that are not present on the anhydrous surface. Proton transfer reactions to CO(2) with formation of bicarbonate and carbonic acid species were investigated and found to have barriers in the range 6.1-12.8 kcal∕mol, with reactions involving participation of two or more water molecules or OH groups having lower barriers than reactions involving a single adsorbed water molecule or OH group. The reactions to form the most stable adsorbed formate and bicarbonate species are exothermic relative to the unreacted adsorbed CO(2) and H(2)O (OH) species, with formation of the bicarbonate species being favored. These results are consistent with single crystal measurements which have identified formation of bicarbonate-type species following coadsorption of CO(2) and water on rutile (110).     

Interaction of Pt clusters with the anatase TiO(2)(101) surface: a first principles study

The adsorption of Pt(n)() (n = 1-3) clusters on the defect-free anatase TiO(2)(101) surface has been studied using total energy pseudopotential calculations based on density functional theory. The defect-free anatase TiO(2)(101) surface has a stepped structure with a step width of two O-Ti bond distances in the (100) plane along the [10] direction and the edge of the step is formed by 2-fold-coordinated oxygen atoms along the [010] direction. For a single Pt adatom, three adsorption sites were found to be stable. Energetically, the Pt adatom prefers the bridge site formed by 2 2-fold-coordinated oxygen atoms with an adsorption energy of 2.84 eV. Electronic structure analysis showed that the Pt-O bonds formed upon Pt adsorption are covalent. Among six stable Pt(2) adsorption configurations examined, Pt(2) was found to energetically favor the O-O bridge sites on the step edge along [010] with the Pt-Pt bond axis perpendicular to [010]. In these configurations, one of the Pt atoms occupies the same O-O bridge site as for a single Pt adatom and the other one either binds a different 2-fold-coordinated oxygen atom on the upper step or a 5-fold-coordinated Ti atom on the lower terrace. Three triangular and three open Pt(3) structures were determined as minima for Pt(3) adsorption on the surface. Platinum trimers adsorbed in triangular structures are more stable than in open structures. In the most stable configuration, Pt(3) occupies the edge O-O site with the Pt(3) plane being upright and almost perpendicular to the [001] terrace. The preference of Pt(n)() to the coordinately unsaturated 2-fold-coordinated oxygen sites indicates that these sites may serve as nucleation centers for the growth of metal clusters on the oxide surface. The increase in clustering energy with increasing size of the adsorbed Pt clusters indicates that the growth of Pt on this surface will lead to the formation of three-dimensional particles.     

Adsorption-site-dependent electronic structure of catechol on the anatase TiO2(101) surface

The adsorption of catechol (1,2-benzendiol) on the anatase TiO(2)(101) surface was studied with synchrotron-based ultraviolet photoemission spectroscopy (UPS), X-ray photoemission spectroscopy (XPS), and scanning tunneling microscopy (STM). Catechol adsorbs with a unity sticking coefficient and the phenyl ring intact. STM reveals preferred nucleation at step edges and subsurface point defects, followed by 1D growth and the formation of a 2 × 1 superstructure at full coverage. A gap state of ～1 eV above the valence band maximum is observed for dosages in excess of ～0.4 Langmuir, but such a state is absent for lower coverages. The formation of the band gap states thus correlates with the adsorption at regular lattice sites and the onset of self-assembled superstructures.     

Successful a priori modeling of CO adsorption on Pt(111) using periodic hybrid density functional theory

The adsorption of CO on the surface of metals such as Pt(111) is of great interest owing to the industrial importance of the catalytic oxidation of pollutant CO. To date, reliable high-level calculations of this process have not been possible, a situation often referred to as the "CO/Pt(111) puzzle". Standard generalized-gradient-approximation density functional theory approaches fail to capture key details of the binding, such as the location of the adsorption site, while cluster approaches using alternative methods show some but insufficient improvement. Using a new computational methodology combining hybrid density functionals containing non-local Hartree-Fock exchange with periodic imaging plane-wave-based techniques, we demonstrate that key aspects of the adsorption of CO on Pt(111), including the identification of the absorption site and CO frequency change, can now be adequately modeled. The binding is dominated by both CO dative covalent bonding and metal-to-molecule pi back-bonding, effects requiring realistic alignment of both the molecular HOMO and LUMO orbitals with respect to the metal Fermi energy.     

Describing weak interactions of biomolecules with dispersion-corrected density functional theory

Interaction energies of the biomolecules in the JSCH-2005 database are calculated with density functional theory using the exchange-correlation functional BLYP augmented with dispersion-corrected atom-centered potentials (DCACPs). The results are in excellent agreement with extrapolated CCSD(T) complete basis set limit references with unsigned mean errors of less than 1.6 kcal mol(-1). Geometry optimisations all reach stable configurations that are close to the MP2-optimised reference geometries.     

Identifying an O2 supply pathway in CO oxidation on Au/TiO2(110): a density functional theory study on the intrinsic role of water

Au catalysis has been one of the hottest topics in chemistry in the last 10 years or so. How O2 is supplied and what role water plays in CO oxidation are the two challenging issues in the field at the moment. In this study, using density functional theory we show that these two issues are in fact related to each other. The following observations are revealed: (i) water that can dissociate readily into OH groups can facilitate O2 adsorption on TiO2; (ii) the effect of OH group on the O2 adsorption is surprisingly long-ranged; and (iii) O2 can also diffuse along the channel of Ti (5c) atoms on TiO2(110), and this may well be the rate-limiting step for the CO oxidation. We provide direct evidence that O2 is supplied by O2 adsorption on TiO2 in the presence of OH and can diffuse to the interface of Au/TiO2 to participate in CO oxidation. Furthermore, the physical origin of the water effects on Au catalysis has been identified by electronic structure analyses: There is a charge transfer from TiO2 in the presence of OH to O2, and the O2 adsorption energy depends linearly on the O2 charge. These results are of importance to understand water effects in general in heterogeneous catalysis.     

CO2 reforming of CH4 on Ni(111): a density functional theory calculation

CO(2) reforming of CH(4) on Ni(111) was investigated by using density functional theory. On the basis of thermodynamic analyses, the first step is CH(4) sequential dissociation into surface CH (CH(4) --> CH(3) --> CH(2) --> CH) and hydrogen, and CO(2) dissociation into surface CO and O (CO(2) --> CO + O). The second step is CH oxygenation into CHO (CH + O --> CHO), which is more favored than dissociation into C and hydrogen (CH --> C + H). The third step is the dissociation of CHO into surface CO and H (CHO --> CO + H). This can explain the enhanced selectivity toward the formation of CO and H(2) on Ni catalysts. It is found that surface carbon formation by the Bouduard back reaction (2CO = C((ads)) + CO(2)) is more favored than by CH(4) sequential dehydrogenation. The major problem of CO(2) reforming of CH(4) is the very strong CO adsorption on Ni(111), which results in the accumulation of CO on the surface and hinders the subsequent reactions and promotes carbon deposition. Therefore, promoting CO desorption should maintain the reactivity and stability of Ni catalysts. The computed energy barriers of the most favorable elementary reaction identify the CH(4) activation into CH(3) and H as the rate-determining step of CO(2) reforming of CH(4) on Ni(111), in agreement with the isotopic experimental results.     

The reaction pathways for HSCH3 adsorption on Au(111): a density functional theory study

Density functional theory was used to investigate the reaction pathways for HSCH(3) adsorption on Au(111) at low coverage. A molecular adsorbed state was found with the S atom bond on Top sites (E approximately -0.38 eV) and molecular adsorption is nonactivated. The H-SCH(3) dissociation process is energetically less favorable and becomes slightly exothermic only when surface relaxation is considered (DeltaE approximately -0.2 eV). All the reaction pathways present a sizable activation energy barrier, with the lowest being approximately 0.52 eV (0.41 eV taking into account slab relaxation). In the corresponding saddle point of the potential energy surface, the S atom of the methylthiolate molecule is placed on Top sites and the H near a Bridge site. The high barrier obtained explains the complete absence of reactive methanethiol dissociation found in recent experiments.     

In situ ATR FTIR study of dextrin adsorption on anatase TiO2

The adsorption of two dextrin-based polymers, a regular wheat dextrin (TY) and a carboxymethyl-substituted (CM) dextrin, onto an anatase TiO(2) particle film has been studied using in situ attenuated total reflection (ATR) FTIR spectroscopy. Infrared spectra of the polymer solutions and the polymer adsorbed at the anatase surface were acquired for two solution conditions: pH 3 and pH 9; below and above the isoelectric point (IEP) of anatase, respectively. Comparison of the polymer solution spectra and the adsorbed layer spectra highlighted a number of spectral differences that were attributed to involvement of the carboxyl group of CM Dextrin interacting with the anatase surface directly and the adsorption of oxidized dextrin chains in the case of regular dextrin (TY) at high pH. The adsorption/desorption kinetics were determined by monitoring spectral peaks of the pyranose ring of both polymers. Adsorption equilibrium was not established for Dextrin TY for many hours, whereas CM Dextrin reached equilibrium in its adsorption within 60 min. The extent of desorption of Dextrin TY (observed by flowing a background electrolyte dextrin-free solution) was extensive at both pH values, which reflects the poor affinity and binding of the polymer on anatase. In contrast, CM Dextrin underwent almost no desorption, indicating a high affinity between the carboxyl groups of the polymer and the anatase surface.     

Density functional theory study of formic acid adsorption on anatase TiO2(001): geometries, energetics, and effects of coverage, hydration, and reconstruction

We present density functional theory calculations and first-principles molecular dynamics simulations of formic acid adsorption on anatase TiO(2)(001), the minority surface exposed by anatase TiO(2) nanoparticles. A wide range of factors that may affect formic acid adsorption, such as coverage, surface hydration, and reconstruction, are considered. It is found that (i) formic acid dissociates spontaneously on unreconstructed clean TiO(2)(001)-1 x 1, as well as on the highly reactive ridge of the reconstructed TiO(2)(001)-1 x 4 surface; (ii) on both the 1 x 1 and 1 x 4 surfaces, various configurations of dissociated formic acid exist with adsorption energies of about 1.5 eV, which very weakly depend on the coverage; (iii) bidentate adsorption configurations, in which the formate moiety binds to the surface through two Ti-O bonds, are energetically more favored than monodentate ones; (iv) partial hydration of TiO(2)(001)-1 x 1 tends to favor the bidentate chelating configuration with respect to the bridging one but has otherwise little effect on the adsorption energetics; and (v) physical adsorption of formic acid on fully hydrated TiO(2)(001)-1 x 1 is also fairly strong. Comparison of the present results for formic acid adsorption with those for water and methanol under similar conditions provides valuable insights to the understanding of recent experimental results concerning the coadsorption of these molecules.     

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.     

Theoretical investigation on RuO2 nanoclusters adsorbed on TiO2 rutile (110) and anatase (101) surfaces

Minimizing the energy of an $N$ -electron system as a functional of a two-electron reduced density matrix (2-RDM), constrained by necessary $N$ -representability conditions (conditions for the 2-RDM to represent an ensemble $N$ -electron quantum system), yields a rigorous lower bound to the ground-state energy in contrast to variational wave function methods. We characterize the performance of two sets of approximate constraints, (2,2)-positivity (DQG) and approximate (2,3)-positivity (DQGT) conditions, at capturing correlation in one-dimensional and quasi-one-dimensional (ladder) Hubbard models. We find that, while both the DQG and DQGT conditions capture both the weak and strong correlation limits, the more stringent DQGT conditions improve the ground-state energies, the natural occupation numbers, the pair correlation function, the effective hopping, and the connected (cumulant) part of the 2-RDM. We observe that the DQGT conditions are effective at capturing strong electron correlation effects in both one- and quasi-one-dimensional lattices for both half filling and less-than-half filling.     

Growth and organization of an organic molecular monolayer on TiO2: catechol on anatase (101)

Anatase TiO(2) is a widely used photocatalytic material, and catechol (1,2-benzendiol) is a model organic sensitizer for dye-sensitized solar cells. The growth and the organization of a catecholate monolayer on the anatase (101) surface were investigated with scanning tunneling microscopy and density functional theory calculations. Isolated molecules adsorb preferentially at steps. On anatase terraces, monodentate ('D1') and bidentate ('D2') conformations are both present in the dilute limit, and frequent interconversions can take place between these two species. A D1 catechol is mobile at room temperature and can explore the most favorable surface adsorption sites, whereas D2 is essentially immobile. When a D1 molecule arrives in proximity of another adsorbed catechol in an adjacent row, it is energetically convenient for them to pair up in nearest-neighbor positions taking a D2-D2 or D2-D1 configuration. This intermolecular interaction, which is largely substrate mediated, causes the formation of one-dimensional catecholate islands that can change in shape but are stable to break-up. The change between D1 and D2 conformations drives both the dynamics and the energetics of this model system and is possibly of importance in the functionalization of dye-sensitized solar cells.     

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.