CO adsorption on Cu(1 1 1) and Cu(0 0 1) surfaces: Improving site preference in DFT calculations

CO adsorption on Cu(1 1 1) and Cu(0 0 1) surfaces has been studied within ab initio density functional theory (DFT). The structural, vibrational and thermodynamic properties of the adsorbate–substrate complex have been calculated. Calculations within the generalized gradient approximation (GGA) predict adsorption in the threefold hollow on Cu(1 1 1) and in the bridge-site on Cu(0 0 1), instead of on-top as found experimentally. It is demonstrated that the correct site preference is achieved if the underestimation of the HOMO–LUMO gap of CO characteristic for DFT is corrected by applying a molecular DFT + U approach. The DFT + U approach also produces good agreement with the experimentally measured adsorption energies, while introducing only small changes in the calculated geometrical and vibrational properties further improving agreement with experiment which is fair already at the GGA level.     

Cluster and band structure ab initio calculations on the adsorption of CO on acid sites of the TiO2(110) surface

The interaction of CO with the cationic sites of the TiO₂(110) surface has been investigated with cluster models and band structure calculations. The analysis of Hartree-Fock and correlated wavefunctions has shown that CO adsorbs at a distance of 4.4–4.5 bohr and a binding energy of 0.7–0.8 eV at low coverage. The bond strength is determined by the CO polarization and the CO б-donation to the surface while the electrostatic attraction is almost exactly cancelled by the Pauli repulsion. The adsorption is accompanied by two important measurable features, a considerable blue shift of the CO vibrational frequency and an increase of the CO 5б ionization potential. Both these effects have been analyzed in detail. We found that the CO ω shift is largely due to a combination of the electrostatic Stark effect and of the repulsion occurring when the CO molecule stretches in the presence of the rigid surface. The change in the 1π–5б binding energies, on the other hand, has an entirely electrostatic origin. Neither the ω shift nor the 5б binding energy shift are determined by the б-donation mechanism. Nevertheless, the occurrence of a charge transfer from CO to the empty levels of the Ti centers is well documented by (a) the energy and dipole moment change associated to this mechanism, (b) the expectation value of a projection operator which measures the charge associated with a given orbital, and (c) the CO dynamic dipole moment. The same analyses also rule out the occurrence of a Ti-to-CO back donation.     

First-principles studies of H2S adsorption and dissociation on metal surfaces

Density functional theory calculations were employed to investigate the molecular and dissociative adsorption of H₂S on the closed packed surfaces of a number of important noble metals (Ag(1 1 1), Au(1 1 1) and Cu(1 1 1)) and transition metals (Ir(1 1 1), Ni(1 1 1), Pd(1 1 1) and Pt(1 1 1)). Energy minima corresponding to adsorbed states were identified with H₂S binding preferentially at the top sites. The adsorption of other S moieties (SH and S) was also examined. SH and S were found to prefer bridge sites and hollow sites, respectively. The binding of H₂S and its S-containing dissociated species is stronger on the transition metals. The elementary reactions of abstraction of H from H₂S to form a surface SH intermediate and abstraction of H from SH to form a surface S intermediate as model pathways for the dissociation of H₂S were examined. Our results suggest that H₂S decomposition on the aforementioned transition metal surfaces is more facile, both thermodynamically and kinetically.     

Density functional theory investigation of the structure of SO2 and SO3 on Cu(1 1 1) and Ni(1 1 1)

Density functional theory (DFT) slab calculations, mainly using the generalised gradient approximation, have been used to investigate the minimum energy structures of molecular SO₂ and SO₃ on Cu(1 1 1) and Ni(1 1 1) surfaces. On Ni(1 1 1) the optimal local adsorption structures are in close agreement with experimental results for both molecular species obtained using the X-ray standing wavefield technique, although for adsorbed SO₂ the energetic difference between two alternative lateral positions of the lying-down molecule on the surface is marginally significant. On Cu(1 1 1) the results for adsorbed SO₂, in particular, were sensitive to the DFT functional used in the calculations, but in all cases failed to reproduce the experimentally-established preference for adsorption with the molecular plane perpendicular to the surface. This result is discussed in the context of previously published DFT results for these species adsorbed on Cu(1 0 0). The optimal geometry found for SO₃ on Cu(1 1 1) is similar to that on Ni(1 1 1), providing agreement with experiment regarding the molecular orientation but not the adsorption site.     

Theoretical study of CO2 activation on Pt(111) induced by coadsorbed K atoms

The adsorption of CO₂ on the clean and potassium-precovered Pt(111) surface has been studied by means of the cluster model approach within the hybrid B3LYP density functional theory-based method. On the clean surface, CO₂ is undistorted and weakly bound. The equilibrium position of this physisorbed species appears at a rather large distance from the surface. However, when coadsorbed K atoms are included in the model, a chemisorbed, bent CO₂ species on top of a surface Pt atom is found. The presence of the coadsorbed K is found to be necessary to promote CO₂ chemisorption and activation, the key step in activating the CO₂ molecule being a direct interaction with the coadsorbate. In addition, the calculated vibrational frequencies for this chemisorbed species are in agreement with available experimental data.     

CO and hydrogen adsorption on Pd(2 1 0)

We have studied the adsorption of CO on Pd(2 1 0) by performing density functional theory (DFT) calculations within the generalized gradient approximation. We find a relatively small corrugation in the CO adsorption energies with the two bridge sites being energetically almost degenerate. CO is furthermore known as a strong poison in heterogeneous catalysis. We have therefore also addressed the coadsorption of CO with atomic hydrogen. There is a significant inhibition of the hydrogen adsorption due to the presence of CO which is analysed in terms of the electronic structure of the adsorbate system.     

The adsorption of H2O on TiO2 and SnO2(110) studied by first-principles calculations

First-principles calculations based on density functional theory and the pseudopotential method have been used to investigate the energetics of H₂O adsorption on the (110) surface of TiO₂ and SnO₂. Full relaxation of all atomic positions is performed on slab systems with periodic boundary conditions, and cases of full and half coverage are studied. Both molecular and dissociative (H₂O→OH⁻+H⁻) adsorption are treated, and allowance is made for relaxation of the adsorbed species to unsymmetrica configurations. It is found that for both TiO₂ and SnO₂ an unsymmetrical dissociated configuration is the most stable. The symmetrical molecularly adsorbed configuration is unstable with respect to lowering of symmetry, and is separated from the fully dissociated configuration by at most a very small energy barrier. The calculated dissociative adsorption energies for TiO₂ and SnO₂ are in reasonable agreement with the results of thermal desorption experiments. Calculated total and local electronic densities of states for dissociatively and molecularly adsorbed configurations are presented, and their relation with experimental UPS spectra is discussed.     

On-surface and sub-surface oxygen on ideal and reconstructed Cu(1 0 0)

In order to understand the first steps of the Cu(1 0 0) oxidation we performed first principles calculations for on-surface and sub-surface oxygen on this surface. According to our calculations, the adsorption energies for all on-surface site oxygen atoms increase, whereas the energies of the sub-surface atoms decrease with the increasing oxygen coverage. At coverage 1 ML and higher on the reconstructed surface, structures including both on- and sub-surface atoms are energetically more favourable than structures consisting only of on-surface adsorbates. On the ideal (1 0 0) surface this change can be perceived at coverage 0.75 ML.     

Comparative DFT study of the adsorption of 1,3-butadiene, 1-butene and 2-cis/trans-butenes on the Pt(1 1 1) and Pd(1 1 1) surfaces

The interaction of 1,3-butadiene, 1-butene and 2-cis/trans-butenes on the Pt(1 1 1) and Pd(1 1 1) surfaces has been studied with density functional theory methods (DFT). The same most stable adsorption modes have been found on both metal surfaces with similar adsorption energies. For 1,3-butadiene the 1,2,3,4-tetra-σ adsorption structure is shown to be the most stable one, in competition with a 1,4-metallacycle-type mode, which is only less stable by 10-12 kJ mol⁻¹. On Pt(1 1 1) these total energy calculations were combined with simulations of the vibrational spectra. This confirms that the 1,2,3,4-tetra-σ adsorption is the most probable adsorption structure, but cannot exclude the 1,4-metallacycle as a minority species. Although similar in type and energy, the adsorption on the Pd(1 1 1) surface shows a markedly different geometry, with a smaller molecular distortion upon adsorption. The most stable adsorption structure for the butene isomers is the di-σ-mode. Similarly to the case of the 1,3-butadiene, the adsorption geometry is closer to the gas phase one on Pd than on Pt, hence explaining the different spectroscopic results, without the previously assumed requirement of a different binding mode. Moreover the present study has shown that the different selectivity observed on Pt(1 1 1) and Pd(1 1 1) for the hydrogenation reaction of butadiene cannot be satisfactory explained by the single comparison of the relative stabilities of 1,3-butadiene and 1-butene on these metals.     

First principle calculations of benzotriazole adsorption onto clean Cu(1 1 1)

The adsorption of benzotriazole (BTAH or C₆N₃H₅) on a Cu(1 1 1) surface is investigated by using first principle density functional theory calculations (VASP). It is found that BTAH can be physisorbed (<0.1 eV) or weakly chemisorbed (∼0.43 eV) onto Cu(1 1 1), and the chemical bond is formed through nitrogen sp² lone pairs. The weak chemisorption can be stabilized by reaction with neighboring protonphilic radicals, like OH⁻. Furthermore, the geometries and associated energies of intermolecular hydrogen bonds between adsorbates on Cu(1 1 1) are also calculated. A model of the first layer of BTAH/BTA⁻ on Cu(1 1 1) surface is developed based on a hydrogen bond network structure.     

Density functional slab model studies of water adsorption on flat and stepped Cu surfaces

Density functional and periodic slab model calculations are performed to study adsorption of water on various Cu surfaces, focusing on monomers and dimers at the planar Cu surfaces and monomers at stepped ones. Single water molecules tend to weakly bind to atop positions with the molecular plane basically parallel to the substrate surface on the planar surfaces or the step plane on the stepped surfaces with negligible structural deformation of water. The experimental adsorption energies of water on the (1 1 1) and (1 0 0) surfaces are about twice as large as the theoretical values of monomerically adsorbed water. This phenomenon is demonstrated to be due to formation of water clusters and/or existence of surface defects. It is revealed that the most favorable hexagonal ring superstructure on Cu(1 1 0) is a four-layer structure, not the commonly accepted bi-layer configuration. We found that the adsorption energy of monomeric water correlates linearly with following quantities, respectively: the bond length and the stretching frequency of the Cu-O bond, the coordination number of the surface Cu atom, the surface work function of the clean surface and the 1b₁ MO energy shift with respect to the value in the gas phase.     

A periodic model for the V2O5-TiO2 (anatase) catalyst. Stability of dimeric species

We present periodic DFT calculations to study the structure of the V₂O₅-TiO₂ (anatase) catalyst. Linear and cyclic dimeric V₂O₅ species represent the active phase. The support TiO₂ (anatase) is represented for the perfect (1 0 0) and (0 0 1) surfaces. The maximum interaction between the active phase and the support is favored, and low coverage is assumed. The most stable models allow the compensation of the surface dangling bonds, and can be understood as a continuation of the bulk anatase structure. The more suitable models for studying reactivity possess uncoordinated atoms available for reactivity, such as terminal oxygen atoms in V₂O₅. Relaxation plays an important role in the adsorption systems, and cannot be discarded when modeling the V₂O₅-anatase catalyst.     

Atomic and molecular adsorption on Pt(1 1 1)

The adsorption of several atomic (H, O, N, S, and C) and molecular (N₂, HCN, CO, NO, and NH₃) species and molecular fragments (CN, CNH₂, NH_2, NH, CH₃, CH₂, CH, HNO, NOH, and OH) on the (1 1 1) facet of platinum, an important industrial and fuel cell catalyst, was studied using self-consistent periodic density functional theory (DFT-GGA) calculations at a coverage of 1/4 ML. The best binding site, energy, and position, as well as an estimated diffusion barrier, of each species were determined. The binding strength for all the species can be ordered as follows: N₂ < NH₃ < HCN < NO < CO < CH₃ < OH < NH₂ < H < CN < NH < O < HNO < CH₂ < NOH < CNH₂ < N < S < CH < C. Although the atomic species generally preferred fcc sites, there was no clear trend in site preference by the molecular species or molecular fragments. The vibrational frequencies of all the stable adsorbates in their best and second best adsorption sites were calculated and found to be in good agreement with experimental values reported in the literature. Finally, the decomposition thermochemistry of NOH, HNO, NO, NH₃, N₂, CO, and CH₃ was analyzed.     

A theoretical analysis of adsorption and dissociation of CH3OH on the stoichiometric SnO2(110) surface

A theoretical analysis based on the Hartree–Fock pseudopotential method and a density-functional theory calculation using a hybrid combination of general gradient approximation with pseudopotential procedure has been carried out to study the adsorption and dissociation of methanol on the stoichiometric SnO₂(110) surface. The dependence of the results upon model system and computing method is discussed. An optimization procedure of adsorbate and substrate atom positions on a six-layer slab model has been selected to characterize the corresponding geometric parameters, adsorption energy and charge-transfer processes related with the molecularly adsorbed CH₃OH and dissociative channels to yield methoxy or methyl fragments. In the high-coverage limit (θ=1), we find that dissociation of the methanol molecule via the heterolytic cleavage of the C---O bond is favoured. At lower coverage (θ=1/2), this channel and the molecularly adsorbed methanol present similar adsorption energies.     

CO adsorption on the Pt(1 1 1) surface: a comparison of a gradient corrected functional and a hybrid functional

The adsorption of CO on the Pt(1 1 1) surface in a pattern has been studied with the gradient corrected functional of Perdew and Wang and the B3LYP hybrid functional. A slab which is periodic in two dimensions is used to model the system. The Perdew-Wang functional incorrectly gives the fcc site as the most favorable adsorption site, in accord with a set of previous studies. The B3LYP functional gives the top site as the preferred site. This confirms results from cluster studies where it was suggested that the different splitting, dependent on the functional, between highest occupied and lowest unoccupied molecular orbital, could be the reason for this change of the adsorption site. This is supported by an analysis based on the projected density of states and the Mulliken population.     

First-principles studies of the oxygen adsorption on unreconstructed and reconstructed Ni(1 1 0) surfaces

First-principles calculations have been performed to investigate the adsorption of oxygen on unreconstructed and reconstructed Ni(1 1 0) surfaces. The energetics, structural, electronic and magnetic properties are given in detail. For oxygen adsorption on unreconstructed surface, (n×1)(n=2,3) substrate with oxygen atom on short-bridge site is found to be the most stable adsorption configuration. Whereas energetically most favorable adsorption phase of reconstructed surface is p(n×1) substrate with oxygen atom located at long-bridge site. Our calculations suggest that the surface reconstruction is induced by the oxygen adsorption. We also find there are redistributions of electronic structure and electron transfer from the substrate to adsorbate. Our calculations also indicate surface magnetic moment is enhanced on clean surfaces and oxygen atoms are magnetized weakly after oxygen adsorption. Interestingly, adsorption on unreconstructed surface does not change surface magnetic moment. However, adsorbate leads to reduction of surface magnetic moment in reconstructed system remarkably.     

Ti+(CO2)2Ar和Ti+(CO2)n(n=3-7)络合物离子的红外光解离光谱 (cited: 1)

利用脉冲激光溅射-超声分子束载带方法制备了气相Ti⁺(CO₂)₂Ar和Ti⁺(CO₂)_n(n=3-7)络合物离子．采用红外光解离光谱研究了这些选定的质量离子的振动光谱． 对于每一种络合物离子, 在CO伸缩振动频率范围都观察到了振动峰,表明这些离子具有插入的OTi⁺CO(CO₂)_n-1结构． 对于n≦5的OTi⁺CO(CO₂)_n-1离子,其CO振动和CO₂的反对称伸缩振动频率都比自由的CO和CO₂的频率要高,表明CO和CO₂配体与中心金属离子之间主要是静电相互作用．实验结果还表明TiO⁺可以直接络合五个配体(1个CO和4个CO₂分子)．对于n=2络合物体系,除了插入的OTi⁺CO(CO₂)结构以外,还观察到了具有弯曲结构的OCO-Ti⁺-OCO异构体的存在     

Theoretical study of the chemisorption of CO on bimetallic RhCu surfaces and nanoparticles

The CO interaction with bimetallic RhCu surface models representing several compositions has been studied by first principles density functional theory calculations. The analysis of the bare bimetallic clusters Rh(4s) and Cu(3s) core-level binding energies indicates that is not possible to extract information about the oxidation state of the alloy components. The present calculations predict that CO does always sit on top sites, the influence of the alloy composition on the equilibrium geometry and vibrational frequency of CO chemisorbed at a given Rh or Cu site being very small. However, there is a large difference in the structural properties corresponding to CO chemisorption above either Rh or Cu. Therefore, the absolute value of the vibrational frequency of chemisorbed CO does not permit to extract any information about the alloy composition but afford to assign the chemisorption site. Finally, the CO adsorption energy does not follow a monotonic trend with composition. The use of the Constrained Space Orbital Variation analysis permits one to firmly establish that the difference in adsorption energy for different compositions cannot be explained through differences in the σ-donation and π-backdonation mechanisms.     

DFT study of Pt adsorption on low index SrTiO3 surfaces: SrTiO3(1 0 0), SrTiO3(1 1 1) and SrTiO3(1 1 0)

Density Functional Theory has been used to determine the energetically preferred structures of submonolayer, monolayer, and multilayer Pt films on both ideal terminations of SrTiO₃(1 0 0), SrTiO₃(1 1 1), and SrTiO₃(1 1 0). The strength of the resulting metal/metal oxide interfaces was characterized by the adsorption energy of the film and the film’s work of separation. The two polar surfaces, SrTiO₃(1 1 1) and SrTiO₃(1 1 0), form significantly stronger interfaces than the non-polar SrTiO₃(1 0 0) surface. Approximate criteria were applied to predict the growth mode of Pt on each surface.     

Quantum states of hydrogen (H, D, T) atoms on Cu(1 0 0) and (1 1 0) surfaces

We investigate the quantum mechanical behavior of adsorbed hydrogen (H, D, T) on Cu(1 0 0) and (1 1 0) surfaces. We construct potential energy surfaces (PESs) for the motion of the hydrogen H atom on Cu(1 0 0) and (1 1 0) surfaces within the framework of density functional theory. The potential energy takes a minimum value on the hollow site of Cu(1 0 0) and on the short bridge site of Cu(1 1 0). Moreover, we calculate the quantum states of hydrogen atom motion on these calculated PESs. The ground state wave function of the hydrogen atom motion is strongly localized around the hollow site on the Cu(1 0 0) surface. On the other hand, the ground state wave function of the hydrogen atom motion on Cu(1 1 0) is distributed from the short bridge site to two neighboring pseudo-threefold sites. We finally show isotope effects on the quantum states of the motion of hydrogen on both surfaces.