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.     

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.     

Pt interactions with annealed and chemically-etched Nb-doped SrTiO3(0 0 1) surfaces: Metal/oxide surface chemical effects on band bending behavior

XPS and LEED have been used to characterize the interaction of sputter-deposited Pt (maximum coverage <5 ML) with Nb-doped SrTiO₃(0 0 1) surfaces prepared either by annealing in O₂ and then UHV, or by chemical-etching in aqua regia. The annealed surface exhibits an ordered (1 × 1) LEED pattern, with additional diffraction spots and streaks indicating the presence of oxygen vacancies. Increasing Pt coverage results in the decrease of the observed Pt(4f_7/2) binding energy and the uniform shift of the Sr(3d), Ti(2p) and O(1s) levels to smaller binding energies, as expected for Pt cluster growth and surface-to-Pt charge donation on an n-type semiconductor. The etched surface is disordered, and exhibits a hydroxylated surface with a contaminant C film of ∼23 ? average thickness. Pt deposition on the etched surface results in an immediate decrease in the intensity of the OH feature in the O(1s) spectrum, and a uniform shift of the Sr(3d), Ti(2p) and O(1s) levels to larger binding energies with increasing Pt coverage. The observed Pt(4f_7/2) binding energy on the etched surface (∼72 eV) is independent of Pt coverage, and indicates substantial electronic charge donation from the Pt to surface hydroxyl species. The observation of band bending towards higher binding energies upon Pt deposition (behavior normally associated with p-type semiconductors) demonstrates that sub-monolayer quantities of adsorbates can alter metal/oxide interfacial charge transfer and reverse the direction of band bending, with important consequences for Schottky barrier heights and device applications.     

A theoretical study of CH4 dissociation on Pt(1 0 0) surface

Density functional theory calculations have been performed on the adsorption of H and CH₃, and the dissociation of CH₄ on Pt(1 0 0) surface. It was found that H was adsorbed on the top and bridge sites, while CH₃ was adsorbed only on the top site. The coadsorption of methyl and hydrogen which has also been investigated shows that the interaction between the two adsorbates is stabilising. In addition, two distinct pathways were explored, differing by the initial adsorbed state of CH₄. They converge readily to the same transition state corresponding to an activation energy value of 0.53 eV. These results compare favourably with existing data in the literature for Pt(1 1 1) and Pt(1 1 0).     

Adsorption of water on (0 0 1) surface of SrTiO3 and SrZrO3 cubic perovskites: Hybrid HF-DFT LCAO calculations

Hybrid HF-DFT LCAO simulations of (0 0 1) surface properties and water adsorption on cubic SrZrO₃ and SrTiO₃ perovskites are performed using a single slab model framework. Three slab models with the different surface termination including 6-11 atomic planes were used for calculations. The effect of the symmetry reduction and the role of an extra atomic layer basis set have been considered for the bare surface slabs. The optimized structures and water adsorption energies have been calculated for the various types of surface coverage. It is shown that the formation of H-bonds between the water hydrogens and surface oxygens, as well as between the water molecules themselves, controls the structure of the water adsorption layers on the perovskite surfaces. Obtained results indicate that the dissociative type of water adsorption is the energetically more favourable for SrO-terminated zirconate surface than for similar titanate surface giving evidence to the more basic nature of oxygen atoms on SrO-terminated SrZrO₃ surface.     

Theoretical simulation of butane isomers adsorption on a Pt(1 0 0) surface

The adsorption of the two butane isomers on Pt(1 0 0) has been characterised with use of density functional simulations. The adsorption energies corresponding to various adsorption configurations were evaluated in good agreement with experimental values. Limited changes of the molecular structure were evidenced. The C-H bond length increases at a degree depending on the surface-hydrogen distance, while the C-C bond length remains similar to that of the free molecule. The surface on-top Pt sites exert a preferential attraction on the molecule, probably through the interaction with the H atoms. The local density of states curves around H as well as C of the adsorbed molecules show dispersed states below the metal Fermi level indicating a molecule-Pt mixing demonstrating a chemical interaction.     

Facetting and (n × 1) reconstructions of SrTiO3(1 1 0) surfaces

(n × 1) reconstructions and facetting of the (1 1 0) polar surface of SrTiO₃ are studied by means of a combination of shell model and density functional calculations. The polarity compensation can be achieved through the formation of {1 0 0} nano-facets, which play a crucial role in the reconstruction process. The behaviors of various possible terminations (Sr, Ti, and O) are analyzed, as well as their atomic structure and energetics. Their stability in different chemical environments is discussed, with respect to previous formulations and experimental results. The Sr-terminated surface tends to expose large facets, while the TiO and O terminations are marginally stabilized or even destabilized by (n × 1) reconstructions, respectively. Trend to facetting results from a subtle competition between the thermodynamic stability of the ideal non stoichiometric (n × 1) surfaces, and huge atomic relaxations that contribute to the lowering of the surface energy differently for each termination.     

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.     

The kaolinite (0 0 1) polar basal plane

Kaolinite is thought to be an efficient ice nucleating agent because kaolinite crystals expose perfect unreconstructed (0 0 1) basal surfaces which provide a suitable template for ice growth. However, we show here with the aid of density functional theory calculations that the unreconstructed basal surface is polar. Various mechanisms to eliminate the macroscopic dipole and in so-doing stabilize the basal surface are considered. The most promising option identified so far involves the adsorption of foreign atoms (for example, Na and Cl atoms) on the perfect basal surface, since this yields a non-metallic surface with a cleavage energy lower than the unreconstructed polar basal surface. A quantitative experimental structure determination of the kaolinite basal surface is called for.     

The small unit cell reconstructions of SrTiO3(1 1 1)

We analyze the basic structural units of simple reconstructions of the (1 1 1) surface of SrTiO₃ using density functional calculations. The prime focus is to answer three questions: what is the most appropriate functional to use; how accurate are the energies; what are the dominant low-energy structures and where do they lie on the surface phase diagram. Using test calculations of representative small molecules we compare conventional PBE-GGA with higher-order methods such as the TPSS meta-GGA and on-site hybrid methods PBE0 and TPSSh, the later being the most accurate. There are large effects due to reduction of the metal d oxygen sp hybridization when using the hybrid methods which are equivalent to a dynamical GGA + U, which leads to rather substantial improvements in the atomization energies of simple calibration molecules, even though the d-electron density for titanium compounds is rather small. By comparing the errors of the different methods we are able to generate an estimate of the theoretical error, which is about 0.25 eV per 1 × 1 unit cell, with changes of 0.5-1.0 eV per 1 × 1 cell with the more accurate method relative to conventional GGA. An analysis of the plausible structures reveals a new low-energy TiO₂-rich configuration with octahedral co-ordination. This structure can act as a template for layers of either TiO or Ti₂O₃, consistent with experimental results. The results also suggest that both the fracture surface and the stoichiometric SrTiO₃(1 1 1) surface should spontaneously disproportionate into SrO and TiO₂ rich domains.     

DFT study on dissociative adsorption of SiH4 and GeH4 on SiGe(1 0 0)-2 × 1 surface

SiH₄ and GeH₄ dissociative adsorptions on a buckled SiGe(1 0 0)-2 × 1 surface have been analyzed using density functional theory (DFT) at the B3LYP level. The Ge alloying in the Si(1 0 0)-2 × 1 surface affects the dimer buckling and its surface reactivity. Systematic Ge influences on the reaction energetics are found in SiH₄ and GeH₄ reactions with four dimers of Si^∗-Si, Ge^∗-Si, Ge^∗-Ge, and Si^∗-Ge (∗ denotes the protruded atom). On a half H-covered surface, the energy barriers for silane and germane adsorption are higher than those on the pristine surface. The energy barrier for silane adsorption is higher than the corresponding barrier for germane adsorption. Rate constants are also calculated using the transition-state theory. We conclude that the SiGe surface reactivity in adsorption reaction depends on the Ge presence in dimer form. If the surface Ge is present in form of Ge^∗-Ge, the surface reactivity decreases as the Ge^∗-Ge content increases. If the surface Ge prefers to be in form of Ge^∗-Si at low Ge contents, the surface reactivity increases first, then decreases at high Ge surface contents when Ge^∗-Ge prevails. The calculated rate constant ratio of GeH₄ adsorption on Si^∗-Si over Ge^∗-Ge at 650 °C is 2.1, which agrees with the experimental ratio of GeH₄ adsorption probability on Si(1 0 0) over Si(1 0 0) covered by one monolayer Ge. The experimental ratio is 1.7 measured through supersonic molecular beam techniques. This consistency between calculation and experimental results supports that one monolayer of Ge on Si(1 0 0) exists in form of Ge^∗-Ge dimer.     

CO methanation on Ni(1 1 1) and modified Ni3Al(1 1 1) surfaces: A first-principle study

The adsorption energies of intermediates in CO methanation on the modified Ni₃Al(1 1 1) surface and the Ni(1 1 1) surface are calculated using density functional theory. A microkinetic analysis based on the calculated adsorption energies is performed to explain the different kinetics of CO methanation catalyzed by Ni₃Al and Ni powders. The electronic structures of different atoms on the modified Ni₃Al surface are also presented, and correlate well with the adsorption energies and geometries.     

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.     

Density functional study of hypophosphite adsorption on Ni (1 1 1) and Cu (1 1 1) surfaces

Surface structures and electronic properties of hypophosphite, H₂PO₂⁻, molecularly adsorbed on Ni(1 1 1) and Cu(1 1 1) surfaces are investigated in this work by density functional theory at B3LYP/6-31++g(d, p) level. We employ a four-metal-atom cluster as the simplified model for the surface and have fully optimized the geometry and orientation of H₂PO₂⁻ on the metal cluster. Six stable orientations have been discovered on both Ni (1 1 1) and Cu (1 1 1) surfaces. The most stable orientation of H₂PO₂⁻ was found to have its two oxygen atoms interact the surface with two PO bonds pointing downward. Results of the Mulliken population analysis showed that the back donation from 3d orbitals of the transition metal substrate to the unfilled 3d orbital of the phosphorus atom in H₂PO₂⁻ and 4s orbital's acceptance of electron donation from one lone pair of the oxygen atom in H₂PO₂⁻ play very important roles in the H₂PO₂⁻ adsorption on the transition metals. The averaged electron configuration of Ni in Ni₄ cluster is 4s^0.634p^0.023d^9.35 and that of Cu in Cu₄ cluster is 4s^1.004p^0.033d^9.97. Because of this subtle difference of electron configuration, the adsorption energy is larger on the Ni surface than on the Cu surface. The amount of charge transfers due to above two donations is larger from H₂PO₂⁻ to the Ni surface than to the Cu surface, leading to a more positively charged P atom in Ni_nH₂PO₂⁻ than in Cu_nH₂PO₂⁻. These results indicate that the phosphorus atom in Ni_nH₂PO₂⁻ complex is easier to be attacked by a nucleophile such as OH⁻ and subsequent oxidation of H₂PO₂⁻ can take place more favorably on Ni substrate than on Cu substrate.     

TiO2-rich reconstructions of SrTiO3(0 0 1): a theoretical study of structural patterns

We have recently reported structure solutions for the (2 × 1) and c(4 × 2) reconstructions of SrTiO₃(0 0 1) based on high-resolution electron microscopy, direct methods analysis of diffraction data and density functional theory. Both reconstructions were found to be TiO₂-rich and feature a single overlayer of TiO₂ stoichiometry on top of a bulk-like TiO₂ layer. Qualitatively, the two reconstruction geometries differ in the cation sub-lattice of the overlayer only, where Ti atoms occupy half of the fivefold cation sites. In the present work we use density functional theory to generate a number of variations of this structural motif in search of patterns of stability. We find a reliable predictor for the reconstruction energy in the ability of oxygen atoms to relax vertically out of the overlayer plane to minimize non-bonded oxygen-oxygen repulsions. Out-of-plane relaxation of oxygen atoms in turn is modulated by the number and relative position of coordinating Ti atoms, which yields simple empirical rules as to how cations are distributed in low energy reconstructions.     

First-principles study on field evaporation of surface atoms from W(0 1 1) and Mo(0 1 1) surfaces

The simulations of field-evaporation processes for surface atoms on W(0 1 1) and Mo(0 1 1) surfaces are implemented using first-principles calculations based on the real-space finite-difference method. The threshold values of the external electric field for evaporation of the surface atoms, which are ∼6 V/Å for tungsten and ∼5 V/Å for molybdenum, are in agreement with the experimental results. While the threshold value of the electric field and the local-field enhancement around the evaporating atoms agree with those expected from the conclusion of the previous study using structureless jellium, the induced charge around the surface atom has a significant difference from that obtained by the jellium model.     

Density functional study of oxygen adsorption on the Mg3Nd (0 0 1) surface

The adsorption of oxygen atoms on Mg₃Nd (0 0 1) surface was studied based on density function theory (DFT), in which the exchange-correlation potential was chosen as the generalized gradient approximation (GGA) in the Perdew and Wang (PW91). The most preferred adsorption position was at the top-hollow site. Upon the optimization on top-hollow site with different coverage, it was found that the adsorption energy decreased with oxygen coverage. The density of states analysis showed that obvious charge transfer took place between O atom and the nearest Nd atom and chemical bond formed between O atom and the nearest Nd atom after O adsorption. The result of surface energy as a function of chemical potential change of oxygen indicated the clean Mg₃Nd (0 0 1) surface was easy to adsorb oxygen and form 1.00 ML surface.     

Effects of N adsorption on the structural and electronic properties of SrTiO3(0 0 1) surface

The atomic configurations, bonding characteristics, and electronic structures of the N-adsorbed (directly and substitutionally) SrTiO₃(0 0 1) surface are studied by using first-principles method based on the density functional theory. From the analysis of the energetics and density of states, it is found that the stability of the directly adsorbed N depends on the relative position of N atom to the surface. To better understand the effects of the substitutionally adsorbed N on the surface, as an example, the behavior of Au atoms adsorbed on the N-substituted surface is discussed in detail. There is clearly a synergy effect between the substitution of N for O_s(I) and the adsorption of Au atoms on the SrTiO₃(0 0 1) surface.     

First-principles study of adsorption and diffusion on Ge/Si(0 0 1)-(2 × 8) and Ge/Si(1 0 5)-(1 × 2) surfaces

Using first-principles total-energy calculations, we have investigated the adsorption and diffusion of Si and Ge adatoms on Ge/Si(0 0 1)-(2 × 8) and Ge/Si(1 0 5)-(1 × 2) surfaces. The dimer vacancy lines on Ge/Si(0 0 1)-(2 × 8) and the alternate S_A and rebonded S_B steps on Ge/Si(1 0 5)-(1 × 2) are found to strongly influence the adatom kinetics. On Ge/Si(0 0 1)-(2 × 8) surface, the fast diffusion path is found to be along the dimer vacancy line (DVL), reversing the diffusion anisotropy on Si(0 0 1). Also, there exists a repulsion between the adatom and the DVL, which is expected to increase the adatom density and hence island nucleation rate in between the DVLs. On Ge/Si(1 0 5)-(1 × 2) surface, the overall diffusion barrier of Si(Ge) along direction is relative fast with a barrier of ∼0.83(0.61) eV, despite of the large surface undulation. This indicates that the adatoms can rapidly diffuse up and down the (1 0 5)-faceted Ge hut island. The diffusion is also almost isotropic along [0 1 0] and directions.