A density functional theory study of the surface relaxation and reactivity of Cu2O(100)

Density functional theory has been used to investigate the surface relaxation of Cu2O(100) and the adsorption of NO. The calculations indicate the formation of surface copper dimers on relaxation coupled with a large contraction of the spacing between the first and second layers. Local density of states for atoms in the top three layers shows that the third layer copper atoms have the greatest change in bonding character. Adsorption energies have been calculated for the N-down and O-down adsorption of NO on the Cu2O(100) surface. These indicate that N-down adsorption is favoured and that in this case NO-lattice oxygen interactions dominate the adsorbate structure. (C) 2000 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.

Density-potential mapping in time-dependent density-functional theory

The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).

New insight into mechanisms in water-gas-shift reaction on Au/CeO2(111): A density functional theory and kinetic study

Using density functional theory (DFT) and kinetic analyses, a new carboxyl mechanism for the water-gas-shift reaction (WGSR) on Au/CeO2(111) is proposed. Many elementary steps in the WGSR are studied using an Au cluster supported on CeO2(111). It is found that (i) water can readily dissociate at the interface between Au and CeO2; (ii) CO2 can be produced via two steps: adsorbed CO on the Au cluster reacts with active OH on ceria to form the carboxyl (COOH) species and then COOH reacts with OH to release CO2; and (iii) two adsorbed H atoms recombine to form molecular H-2 on the Au cluster. Our kinetic analyses show that the turnover frequency of the carboxyl mechanism is consistent with the experimental one while the rates of redox and formate mechanisms are much slower than that of carboxyl mechanism. It is suggested that the carboxyl pathway is likely to be responsible for WGSR on Au/CeO2.

Dispersion interactions in room-temperature ionic liquids: Results from a non-empirical density functional

The role of dispersion or van de Waals (VDW) interactions in imidazolium-based room-temperature ionic liquids is studied within the framework of density functional theory, using a recently developed non-empirical functional [M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)], as efficiently implemented in the SIESTA code [G. Roman-Perez and J. M. Soler, Phys. Rev. Lett. 103, 096102 (2009)]. We present results for the equilibrium structure and lattice parameters of several crystalline phases, finding a general improvement with respect to both the local density (LDA) and the generalized gradient approximations (GGA). Similar to other systems characterized by VDW bonding, such as rare gas and benzene dimers as well as solid argon, equilibrium distances and volumes are consistently overestimated by approximate to 7%, compared to -11% within LDA and 11% within GGA. The intramolecular geometries are retained, while the intermolecular distances and orientations are significantly improved relative to LDA and GGA. The quality is superior to that achieved with tailor-made empirical VDW corrections ad hoc [M. G. Del Popolo, C. Pinilla, and P. Ballone, J. Chem. Phys. 126, 144705 (2007)]. We also analyse the performance of an optimized version of this non-empirical functional, where the screening properties of the exchange have been tuned to reproduce high-level quantum chemical calculations [J. Klimes, D. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 074203 (2010)]. The results for solids are even better with volumes and geometries reproduced within 2% of experimental data. We provide some insight into the issue of polymorphism of [bmim][Cl] crystals, and we present results for the geometry and energetics of [bmim][Tf] and [mmim][Cl] neutral and charged clusters, which validate the use of empirical force fields. (C) 2011 American Institute of Physics. [doi:10.1063/1.3652897]

Acetaldehyde Production in the Direct Ethanol Fuel Cell: Mechanistic Elucidation by Density Functional Theory

This study employs density functional theory (DFT) calculations to examine the mechanism by which acetaldehyde is formed on platinum in a typical direct ethanol fuel cell (DEFC). A pathway is found involving the formation of a strongly hydrogen-bonded complex between adsorbed ethanol and the surface hydroxyl (OH) species, followed by the facile alpha-dehydrogenation of ethanol, with spontaneous weakening of the hydrogen bond in favor of adsorbed acetaldehyde and water. This mechanism is found to be comparably viable on both the close-packed surface and the monatomic steps. Comparison of further reactions on these two sites strongly indicates that the steps act as net removers of acetaldehyde from the product stream, while the flat surface acts as a net producer.

Density Functional Theory Study on the Cleavage Mechanism of the Carbonyl Bond in Amides on Flat and Stepped Ru Surfaces: Hydrogen-Induced or Direct C-O Bond Breaking?

We have performed density functional theory (DFT) calculations to investigate the reaction mechanism of the cleavage of the carbonyl bond in amides on both flat and stepped Ru surfaces. The simplest amide molecule, N,N-dimethylacetamide (DMA), was used as the exemplar model molecule. Through the calculations, the most stable transition states (TSs) in all the pathways on both flat and stepped Ru surfaces are identified. Comparing the energy profiles of different reaction pathways, we find that a direct cleavage mechanism is always energetically favored as compared with an alternative hydrogen-induced mechanism on either the flat or stepped Ru surface. It is easier for the dissociation process to occur on the stepped surface than on the flat surface. However, as compared with the terrace, the superiority of step sites boosting the C-O bond dissociation is not as evident as that on CO dissociation.

Time-dependent density functional theory study of charge transfer in collisions

We study the charge transfer between colliding ions, atoms, or molecules, within time-dependent density functional theory. Two particular cases are presented, the collision between a proton and a Helium atom, and between a gold atom and a butane molecule. In the first case, proton kinetic energies between 16 keV and 1.2 MeV are considered, with impact parameters between 0.31 and 1.9 angstrom. The partial transfer of charge is monitored with time. The total cross-section is obtained as a function of the proton kinetic energy. In the second case, we analyze one trajectory and discuss spin-dependent charge transfer between the different fragments.

Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Effect of spatial nonlocality on the density functional band gap

We show that for a large class of exchange-correlation functionals the local exchange-correlation potential obtained within an optimized effective potential severely underestimates the band gap. On the other hand, the corresponding nonlocal potential obtained from a generalized Kohn-Sham scheme provides a much better description of the band gap, in good agreement with experiments. These results strongly indicate that a local exchange-correlation potential, however good the exchange-correlation approximation, cannot capture the delicate interplay between correlation effects and spatial localization in the KS band structure, unless the (cumbersome) contribution from the derivative discontinuity of the exchange-correlation energy functional is considered.

Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.

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.

Influence of surface structures, subsurface carbon and hydrogen, and surface alloying on the activity and selectivity of acetylene hydrogenation on Pd surfaces: A density functional theory study

The selective hydrogenation of acetylene to ethylene on several Pd surfaces (Pd(111), Pd(100), Pd(211), and Pd(211)-defect) and Pd surfaces with subsurface species (carbon and hydrogen) as well as a number of Pd-based alloys (Pd-M/Pd(111) and Pd-M/Pd(211) (M = Cu, Ag and Au)) are investigated using density functional theory calculations to understand both the acetylene hydrogenation activity and the selectivity of ethylene formation. All the hydrogenation barriers are calculated, and the reaction rates on these surfaces are obtained using a two-step model. Pd(211) is found to have the highest activity for acetylene hydrogenation while Pd(100) gives rise to the lowest activity. In addition, more open surfaces result in over-hydrogenation to form ethane, while the close-packed surface (Pd(111)) is the most selective. However, we also find that the presence of subsurface carbon and hydrogen significantly changes the reactivity and selectivity of acetylene toward hydrogenation on Pd surfaces. On forming surface alloys of Pd with Cu, Ag and Au, the selectivity for ethylene is also found to be changed. A new energy decomposition method is used to quantitatively analyze the factors in determining the changes in selectivity. These surface modifiers are found to block low coordination unselective sites, leading to a decreased ethane production. (C) 2013 The Authors. Published by Elsevier Inc. All rights reserved.

Selective Hydrogenation of Acetylene over Pd-Boron Catalysts: A Density Functional Theory Study

Boron-modified Pd catalysts have shown excellent performance for the selective hydrogenation of alkynes experimentally. In the current work, we investigated the hydrogenation of acetylene on boron-modified Pd(111) and Pd(211) surfaces, utilizing density functional theory calculations. The activity of acetylene hydrogenation has been studied by estimating the effective barrier of the whole process. The selectivity of ethylene formation is investigated from a comparison between the desorption and the hydrogenation of ethylene as well as comparison between the ethylene and the 1,3-butadiene formation. Formation of subsurface carbon and hydrogen on both boron-modified Pd(111) and Pd(211) surfaces has also been evaluated, since these have been reported to affect both the activity and the selectivity of acetylene hydrogenation to produce ethylene on Pd surfaces. Our results provide some important insights into the Pd B catalysts for selective hydrogenation of acetylene and also for more complex hydrogenation systems, such as stereoselective hydrogenation of longer chain alkynes and selective hydrogenation of vegetable oil.