DFT and in situ EXAFS investigation of gold/ceria-zirconia low-temperature water gas shift catalysts: Identification of the nature of the active form of gold

A combined experimental and theoretical investigation of the nature of the active form of gold in oxide-supported gold catalysts for the water gas shift reaction has been performed. In situ extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) experiments have shown that in the fresh catalysts the gold is in the form of highly dispersed gold ions. However, under water gas shift reaction conditions, even at temperatures as low as 100 degrees C, the evidence from EXAFS and XANES is only 14 consistent with rapid, and essentially complete, reduction of the gold to form metallic clusters containing about 50 atoms. The presence of Au-Ce distances in the EXAFS spectra, and the fact that about 15% of the gold atoms can be reoxidized after exposure to air at 150 degrees C, is indicative of a close interaction between a fraction (ca. 15%) of the gold atoms and the oxide support. Density functional theory (DFT) calculations are entirely consistent with this model and suggest that an important aspect of the active and stable form of gold under water gas shift reaction conditions is the location of a partially oxidized gold (Audelta+) species at a cerium cation vacancy in the surface of the oxide support. It is found that even with a low loading gold catalysts (0.2%) the fraction of ionic gold under water gas shift conditions is below the limit of detection by XANES (<5%). It is concluded that under water gas shift reaction conditions the active form of gold comprises small metallic gold clusters in intimate contact with the oxide support.

DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

A DFT study of the chain growth probability in Fischer-Tropsch synthesis

The chain growth probability (alpha value) is one of the most significant parameters in Fischer-Tropsch (FT) synthesis. To gain insight into the chain growth probability, we systematically studied the hydrogenation and C-C coupling reactions with different chain lengths on the stepped Co(0001) surface using density functional theory calculations. Our findings elucidate the relationship between the barriers of these elementary reactions and the chain length. Moreover, we derived a general expression of the chain growth probability and investigated the behavior of the alpha value observed experimentally. The high methane yield results from the lower chain growth rate for C-1 + C-1 coupling compared with the other coupling reactions. After C-1, the deviation of product distribution in FT synthesis from the Anderson-Schulz-Flory distribution is due to the chain length-dependent paraffin/olefin ratio. (C) 2008 Elsevier Inc. All rights reserved.

Chain growth mechanism in Fischer-Tropsch synthesis: A DFT study of C-C coupling over Ru, Fe, Rh, and Re surfaces

A quantitative approach is used to understand the chain growth mechanism in FT synthesis on the Ru, Fe, Rh, and Re surfaces. The C-C coupling reactions are extensively calculated on the stepped metal surfaces. Combining the coupling barriers and reactant stabilities, we investigate the reaction rates of all possible C, + C-1 coupling pathways on the metal surfaces. It is found that (i) all the transition-state structures are similar on these surfaces, while some coupling barriers are very different; (ii) the dominant chain growth pathways on these surfaces are different: C + CH and CH + CH on Rh and Ru surfaces, C + CH3 on Fe surface, and C + CH on Re surface. The common features of the major coupling reactions together with those on the Co surface are also discussed.

Reactivity of the 4d transition metals toward N hydrogenation and NH dissociation: A DFT-based HSAB analysis

The density functional theory (DFT) based hard-soft acid-base (HSAB) reactivity indices, including the electrophilicity index, have been successfully applied to many areas of molecular chemistry. In this work we test the applicability of such an approach to fundamental surface chemistry. We have considered, as prototypical surface reactions, both the hydrogenation of atomic nitrogen and the dissociative adsorption of the NH molecular radical. By use of a DFT methodology, the minimum energy reaction pathways, and corresponding reaction barriers, of the above reactions over Zr(001), Nb(110), Mo(110), Tc(001), Ru(001), Rh(111), and Pd(111) have been determined. By consideration of the chemical potential and chemical hardness of the surface metal atoms, and the principle of electronegativity equalization, it is found that the charge transferred to the NH radical during the process of dissociative adsorption correlates very well with that determined by Mulliken population analysis. Furthermore, it is found that the stability of the NH/surface transition state complex relates directly to this charge transfer and that the trend in transition state stability predicted by a HSAB; treatment correlates very strongly with that determined by DFT calculations. With regards to N hydrogenation, we find that during the course of the reaction, H loses cohesion to the surface, as it must migrate from a 3-fold hollow site to either a bridge or top site, to react with N. Partial density of states (PDOS) and Mulliken population analysis reveal that this loss of bonding is accompanied by charge transfer from H to the surface metal atoms. Moreover, by simple modeling, we show that the reaction barriers are directly proportional to this mandatory charge transfer. Indeed, it is found that the reaction barriers correlate very well with the electrophilicity index of the metal atoms.

The temperature dependence of the adsorption of NO on Pt{211}: A RAIRS and DFT investigation

RAIRS experiments have been performed to investigate the adsorption of NO on Pt{211}. Results show that adsorption is complex and strongly temperature dependent. At 307 K, three bands are seen at saturation with frequencies of 1801, 1609, and 1576 cm(-1). However, at 120 K only two bands, at 1688 and 1620 cm(-1), are observed. To help with the assignment of these vibrational bands, DFT calculations were also performed. The calculations show that a bridged NO species, bonded to the step edge, is the most stable species on the surface and gives rise to the band observed at 1610-1620 cm(-1). The calculations also suggest that the temperature dependence of NO adsorption on Ptf{211} can be assigned to NO dissociation which occurs at room temperature but not at 120 K. In particular, the RAIRS band observed at 1801 cm(-1), which is observed on adsorption at 307 K but not at 120 K, is tentatively assigned to the formation of an O-NO complex. This species forms when a NO molecule bonds on top of an O atom, which results from the dissociation of NO on the Pt{211} surface at room temperature.

Resonance Raman and DFT studies of tetra-tert-butyl porphine: Assignment of strongly enhanced distortion modes in a ruffled porphyrin

The free-base form of tetra-tert-butyl porphine (TtBP), which has extremely bulky meso substituents, is severely distorted from planarity, with a ruffling angle of 65.5degrees. The resonance Raman spectrum of TtBP (lambda(ex) = 457.9 nm) and its d(2), d(8), and d(10) isotopomers have been recorded, and while the spectra show high-frequency bands similar to those observed for planar meso-substituted porphyrins, there are several additional intense bands in the low-frequency region. Density functional calculations at the B3-LYP/6-31G(d) level were carried out for all four isotopomers, and calculated frequencies were scaled using a single factor of 0.98. The single factor scaling approach was validated on free base porphine where the RMS error was found to be 14.9 cm(-1). All the assigned bands in the high-frequency (> 1000 cm(-1)) region of TtBP were found to be due to vibrations similar in character to the in-plane skeletal modes of conventional planar porphyrins. In the low-frequency region, two of the bands, assigned as nu(8) (ca. 330 cm(-1)) and nu(16) (ca. 540 cm(-1)), are also found in planar porphyrins such as tetra-phenyl porphine (TPP) and tetra-iso-propyl porphine (IPP). Of the remaining three very strong bands, the lowest frequency band was assigned as gamma(12) (pyr swivel, obsd 415 cm(-1), calcd 407 cm(-1) in do). The next band, observed at 589 cm-1 in the do compound (calcd 583 cm(-1)), was assigned as a mode whose composition is a mixture of modes that were previously labeled gamma(13) (gamma(CmCaHmCa)) andy gamma(11) (pyr fold(asym)) in NiOEP. The final strong band, observed at 744 cm(-1) (calcd 746 cm(-1)), was assigned to a mode whose composition is again a mixture of gamma(11) and gamma(13), although here it is gamma(11) rather than gamma(13) which predominates. These bands have characters and positions similar to those of three of the four porphyrin ring-based, weak bands that have previously been observed for NiTPP. In addition there are several weaker bands in the TtBP spectra that are also