Density functional calculations for polymers and clusters - progress and limitations

The applicability of density functional (DF) methods has progressed greatly since the first workshop of this series ten years ago. Applications that show both the successes and the limitations can be found in the fields of: (a) the structures of the isomers of atomic clusters. and (b) the structure of organic molecules and polymers, and their reactions with additional molecules. We shall review some of the results and the lessons to be learned from them. (C) 2001 Elsevier Science B.V. All rights reserved.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

Density functional/Monte Carlo study of ring-opening polymerization

Density functional calculations of the structure, potential energy surface and reactivity for organic systems closely related to bisphenol-A-polycarbonate (BPA-PC) provide the basis for a model describing the ring-opening polymerization of its cyclic oligomers by nucleophilic molecules. Monte Carlo simulations using this model show a strong tendency to polymerize that is increased by increasing density and temperature, and is greater in 3D than in 2D. Entropy in the distribution of inter-particle bonds is the driving force for chain formation. (C) 2002 Elsevier Science B.V. All rights reserved.

Density functional study of reactions of phenoxides with polycarbonate

Density functional calculations with simulated annealing have been used to study the reactions of chains of bisphenol A polycarbonate (BPA-PC) with sodium phenoxide (NaOPh), diphenyl carbonate (DPC), and tetraphenylphosphonium phenoxide (PPh4OPh). These calculations extend our work on the reactions of LiOPh, NaOPh, and phenol with the cyclic tetramer of BPA-PC. We study, in particular, chain growth catalyzed by NaOPh and PPh4OH. The energy barriers for reactions with PPh4OPh are somewhat larger than those involving LiOPh and NaOPh, but they are significantly lower than those involving phenol (HOPh), due in part to the collective rearrangement of phenyl groups in the reacting molecules. We discuss in the Appendix the bonds between alkali metal atoms (Na in the present calculations) and other atoms (here oxygen) that are analogous to the more familiar "hydrogen bonds".

Examining the redox and formate mechanisms for water-gas shift reaction on Au/CeO2 using density functional theory

We perform DFT calculations to investigate the redox and formate mechanisms of water-gas-shift (WGS) reaction on Au/CeO2 catalysts. In the redox mechanism, we analyze all the key elementary steps and find that the OH cleavage is the key step. Three possible pathways of OH cleavage are calculated: (1) OHad '' + *'--> H-ad' + O-ad"; (2) H-ad' + OHad '' --> H-2(g) + O-ad '' + *'; and (3) OHad" + OHad '' --> 2O(ad '') + H-2(g) (*': the free adsorption sites on the oxides; ad': adsorption on the metal; ad": adsorption on the oxide, respectively). In the formate mechanism, we identify all the possible pathways for the formation and decomposition of surface formates in the WGS reaction. It is found that there is a shortcoming in the redox and formate mechanisms which is related to surface oxygen reproduction. Four possible pathways for producing surface oxygen are studied, and all the barriers of the four pathways are more than 1 eV. Our results suggest that the processes to reproduce surface oxygen in the reaction circle are not kinetically easy. (C) 2008 Elsevier B.V. All rights reserved.

A density functional theory study of CO and atomic oxygen chemisorption on Pt(111)

Ab initio total energy calculations within a density functional theory framework have been performed for CO and atomic oxygen chemisorbed on the Pt(111) surface. Optimised geometries and chemisorption energies for CO and O on four high-symmetry sites, namely the top, bridge, fee hollow and hcp hollow sites, are presented, the coverage in all cases being 0.25 ML. The differences in CO adsorption energies between these sites are found to be small, suggesting that the potential energy surface for CO diffusion across Pt(111) is relatively flat. The 5 sigma and 2 pi molecular orbitals of CO are found to contribute to bonding with the metal. Some mixing of the 4 sigma and 1 pi molecular orbitals with metal states is also observed. For atomic oxygen, the most stable adsorption site is found to be the fee hollow site, followed in decreasing order of stability by the hcp hollow and bridge sites, with the top site being the least stable. The differences in chemisorption energies between sites for oxygen are larger than in the case of CO, suggesting a higher barrier to diffusion for atomic oxygen. The co-adsorption of CO and O has also been investigated. Calculated chemisorption energies for CO on an O/fcc-precovered surface show that of the available chemisorption sites, the top site at the oxygen atom's next-nearest neighbour surface metal atom is the most stable, with the other four sites calculated bring at least 0.29 eV less stable. The trend of CO site stability in the coadsorption system is explained in terms of a 'bonding competition' model. (C) 2000 Elsevier Science B.V. All rights reserved.

Methyl chemisorption on Ni(111) and C-H-M multicentre bonding:a density functional theory study

Density functional theory has been used to study the adsorption of CH3 on Ni(111). CH3 is found to adsorb strongly at all four high symmetry sites of the Ni(111) surface. Calculated adsorption energies of CH3 on the different sites are in the following order: hcp approximate to fcc>bridge>top. The bonding and structures of CH3 on the different sites are analysed in detail. An important factor, namely three-centre bonding between carbon, hydrogen and nickel which contributes to the 'soft' C-H vibrational frequency of CH3 on Ni(111), and may determine the preferred chemisorption site, is stressed. (C) 1999 Elsevier Science B.V. All rights reserved.

A density functional theory study on the interaction between chemisorbed CO and S on Rh(111)

Density functional theory calculations are carried out for Rh(111)-p(2 x 2)-CO, Rh(111)-p(2 x 2)-S, Rh(111)-p(2 x 2)-(S + CO), Rh(111)-p(3 x 3)-CO, Rh(111)-p(3 x 3)-S and Rh(111)-p(3 x 3)-(S + CO), aiming to shed some light on the S poisoning effect. Geometrical structures of these systems are optimized and chemisorption energies are determined. The presence of S does not significantly influence the geometrical structure and chemisorption energy of CO and vice versa, which strongly suggests that the interaction between CO and S on the Rh(111) surface is mainly short-range in nature. The long range electronic effect for the dramatic attenuation of the CO methanation activity by sulfur is likely to be incorrect. It is suggested that an ensemble effect may be dominant in the catalytic deactivation. (C) 1999 Elsevier Science B.V. All rights reserved.

A quantitative determination of reaction mechanisms from density functional theory calculations: Fischer-Tropsch synthesis on flat and stepped cobalt surfaces

We systematically investigated the mechanism of the C-1 + C-1 coupling reactions using density functional theory. The activation energies of C-1 + C-1 coupling and carbon hydrogenation reactions on both flat and stepped surfaces were calculated and analyzed. Moreover, the coverages of adsorbed C-1 species were estimated, and the reaction rates of all possible C-1 + C-1 coupling pathways were quantitatively evaluated. The results suggest that the reactions of CH2 + CH2 and CH3 + C at steps are most likely to be the key C-1 + C-1 coupling steps in FT synthesis on Co catalysts. The reactions of C-2 + C-1 and C-3 + C-1 coupling also were studied; the results demonstrate that in addition to the pathways of RCH + CH2 and RCH2 + C, the coupling of RC + C and RC + CH also may contribute to the chain growth after C-1. (C) 2008 Elsevier Inc. All rights reserved.

Methane transformation to carbon and hydrogen on Pd(100): Pathways and energetics from density functional theory calculations

Density functional theory with gradient corrections has been employed to study the reaction pathways and the reaction energetics for the transformations of CH4 to C and H on a Pd(100) surface. On examination of transition state structures identified in each elementary reaction, a clear relationship between the valencies of the CHx fragments and the locations of the transition states emerges. The higher the valency of the CHx fragment, the higher the coordination number of the CHx with the surface atoms. The calculated reaction energetics are in good agreement with the experiments. In addition, calculation results are also used to illustrate an interesting issue concerning the CH3 stability on Pd surfaces. (C) 2002 American Institute of Physics.

Identifying the distinct features of geometric structures for hole trapping to generate radicals on rutile TiO2(110) in photooxidation using density functional theory calculations with hybrid functional

Using density functional theory calculations with HSE 06 functional, we obtained the structures of spin-polarized radicals on rutile TiO2(110), which is crucial to understand the photooxidation at the atomic level, and further calculate the thermodynamic stabilities of these radicals. By analyzing the results, we identify the structural features for hole trapping in the system, and reveal the mutual effects among the geometric structures, the energy levels of trapped hole states and their hole trapping capacities. Furthermore, the results from HSE 06 functional are compared to those from DFT + U and the stability trend of radicals against the number of slabs is tested. The effect of trapped holes on two important steps of the oxygen evolution reaction, i.e. water dissociation and the oxygen removal, is investigated and discussed.

Some Understanding of Fischer-Tropsch Synthesis from Density Functional Theory Calculations

The combination of density functional theory (DFT) calculations and kinetic analyses is a very useful approach to study surface reactions in heterogeneous catalysis. The present paper reviews some recent work applying this approach to Fischer-Tropsch (FT) synthesis. Emphasis is placed on the following fundamental issues in FT synthesis: (i) reactive sites for both hydrogenation and C-C coupling reactions; (ii) reaction mechanisms including carbene mechanism, CO-insertion mechanism and hydroxyl-carbene mechanism; (iii) selectivity with a focus on CH(4) selectivity, alpha-olefin selectivity and chain growth probability; and (iv) activity.