Adsorbate-substrate and adsorbate-adsorbate interactions of Na and K adlayers on Al(111)

Theoretical investigation of the electronic structure of the Alq(3)/Mg interface

We have studied the atomic geometries and the electronic properties of the tris-(8-hydroxyquinoline) aluminum (Alq(3))/Mg interface by using density functional theoretical calculations. We have found that the chemical bond is formed between the O atoms of Alq(3) and the substrate Mg atoms, and the stability of the interface structures depends on the number of O-Mg bonds. In the up configurations, where two or three O-Mg chemical bonds are formed and the Alq(3) molecular dipoles are oriented up to the vacuum side, the work function is decreased by as much as 1.1 eV or more. The interface dipole is dominated by the orientation of the molecular dipoles of Alq(3). The interface gap state reported from experiments is ascribed to the highest occupied molecular orbital (HOMO) levels of the down configurations, which may coexist with the dominant up configurations.

How does an external electrical field affect adsorption patterns of thiol and thiolate on the gold substrate?

The responsive behavior of methanethiol and methylthiolate molecules on the Au(111) surface with an applied electrical potential is studied, and it is shown how the sulfur adsorption site, the S-H bond orientation and the interacting energy change with an external electric field strength. The electron charge density corresponding to an electric field minus that obtained in zero field, with zero-field optimal geometry, is calculated to explain the responsive behavior. The interacting energy for the intact methanethiol adsorption is larger than that for the dissociative one, showing that an external electric field cannot make the hydrogen dissociate from the sulfur.

Surface magnetoelectric effect in ferromagnetic metal films

A surface magnetoelectric effect is revealed by density-functional calculations that are applied to ferromagnetic Fe(001), Ni(001), and Co(0001) films in the presence of an external electric field. The effect originates from spin-dependent screening of the electric field which leads to notable changes in the surface magnetization and the surface magnetocrystalline anisotropy. These results are of considerable interest in the area of electrically controlled magnetism and magnetoelectric phenomena.

Theory and simulation in heterogeneous gold catalysis

This critical review covers the application of quantum chemistry to the burgeoning area of the heterogeneous oxidation by Au. We focus on the most established reaction, the oxidation of CO at low temperature. The review begins with an overview of the methods available comparing the treatment of the electron-electron interaction and relativistic effects. The structure of Au particles and their interaction with oxide reviews is then discussed in detail. Calculations of the adsorption and reaction of CO and O2 are then considered and results from isolated and supported Au clusters compared (155 references).

Doping graphene with metal contacts

Making devices with graphene necessarily involves making contacts with metals. We use density functional theory to study how graphene is doped by adsorption on metal substrates and find that weak bonding on Al, Ag, Cu, Au, and Pt, while preserving its unique electronic structure, can still shift the Fermi level with respect to the conical point by approximately 0.5 eV. At equilibrium separations, the crossover from p-type to n-type doping occurs for a metal work function of approximately 5.4 eV, a value much larger than the graphene work function of 4.5 eV. The numerical results for the Fermi level shift in graphene are described very well by a simple analytical model which characterizes the metal solely in terms of its work function, greatly extending their applicability.

The reaction pathways for HSCH3 adsorption on Au(111): a density functional theory study

Density functional theory was used to investigate the reaction pathways for HSCH(3) adsorption on Au(111) at low coverage. A molecular adsorbed state was found with the S atom bond on Top sites (E approximately -0.38 eV) and molecular adsorption is nonactivated. The H-SCH(3) dissociation process is energetically less favorable and becomes slightly exothermic only when surface relaxation is considered (DeltaE approximately -0.2 eV). All the reaction pathways present a sizable activation energy barrier, with the lowest being approximately 0.52 eV (0.41 eV taking into account slab relaxation). In the corresponding saddle point of the potential energy surface, the S atom of the methylthiolate molecule is placed on Top sites and the H near a Bridge site. The high barrier obtained explains the complete absence of reactive methanethiol dissociation found in recent experiments.

Chemisorption of hydrogen molecules on carbon nanotubes: charging effect from first-principles calculations

We report a systematic investigation of the charging effect on hydrogen molecule chemisorption on (3, 3), (5, 5), (5, 0), and (8, 0) carbon nanotubes by first-principles calculations. The influence of injected charge on the chemisorption energy barriers is found to be sensitive to the nanotube diameter and chirality. The calculated results also indicate that electron injection is more effective in lowering the energy barrier for armchair carbon nanotubes while hole injection is more effective for zigzag nanotubes. The origin of these interesting trends and systematics can be understood by a close examination of the underlying electronic structure and the electron transfer between the hydrogen molecules and the nanotubes.

On the mechanism of low-temperature water gas shift reaction on copper

Periodic, self-consistent density functional theory (DFT-GGA) calculations are used to investigate the water gas shift reaction (WGSR) mechanism on Cu(111). The thermochemistry and activation energy barriers for all the elementary steps of the commonly accepted redox mechanism, involving complete water activation to atomic oxygen, are presented. Through our calculations, we identify carboxyl, a new reactive intermediate, which plays a central role in WGSR on Cu(111). The thermochemistry and activation energy barriers of the elementary steps of a new reaction path, involving carboxyl, are studied. A detailed DFT-based microkinetic model of experimental reaction rates, accounting for both the previous and the new WGSR mechanism show that, under relevant experimental conditions, (1) the carboxyl-mediated route is the dominant path, and (2) the initial hydrogen abstraction from water is the rate-limiting step. Formate is a stable "spectator" species, formed predominantly through CO2 hydrogenation. In addition, the microkinetic model allows for predictions of (i) surface coverage of intermediates, (ii) WGSR apparent activation energy, and (iii) reaction orders with respect to CO, H2O, CO2, and H2.

Graphite and graphene as perfect spin filters

Based upon the observations (i) that their in-plane lattice constants match almost perfectly and (ii) that their electronic structures overlap in reciprocal space for one spin direction only, we predict perfect spin filtering for interfaces between graphite and (111) fcc or (0001) hcp Ni or Co. The spin filtering is quite insensitive to roughness and disorder. The formation of a chemical bond between graphite and the open d-shell transition metals that might complicate or even prevent spin injection into a single graphene sheet can be simply prevented by dusting Ni or Co with one or a few monolayers of Cu while still preserving the ideal spin-injection property.

Effect of carbon adsorption on the isomer stability of Ir4 clusters

The atomic structure and electronic properties of gas-phase and MgO100-supported iridium tetramers are studied using density functional theory. At variance with experimental data, the most stable Ir4 isomer on MgO100 is the square one, as in the gas phase, and the metastable tetrahedral isomer is highly distorted by interactions with the substrate. In the presence of a single carbon adatom, the most stable structure of Ir4 is tetrahedral for both environments and the structural distortion of the adsorbed cluster is reduced. On MgO100, the binding energy of a C adatom to tetrahedral Ir4 is 1.6 eV larger than that to the square isomer, due to strong interactions between C-2p orbitals and a low-energy unoccupied molecular orbital of tetrahedral Ir4.

Adsorption states of methylthiolate on the Au(111) surface

We have studied the adsorption structures of methylthiolate (MeS) on the Au(111) surface using density functional theory. The most stable adsorption structure of MeS is the MeS-Au-MeS configuration, which has been proposed by Maksymovych and co-workers (2006 Phys. Rev. Lett. 97 146103), followed by the bridge configuration. The calculated work function change and vibrational mode analysis for the two configurations agree quite well with the experimental results. These results strongly support the conclusion that the MeS-Au-MeS and the bridge configurations coexist on the MeS adsorbed Au(111) surface.

Charged fullerenes as high-capacity hydrogen storage media

Using first-principles calculations within density functional theory, we explore systematically the capacity of charged carbon fullerenes Cn (20 <or= n <or= 82) as hydrogen storage media. We find that the binding strength of molecular hydrogen on either positively or negatively charged fullerenes can be dramatically enhanced to 0.18-0.32 eV, a desirable range for potential room-temperature, near ambient applications. The enhanced binding is delocalized in nature, surrounding the whole surface of a charged fullerene, and is attributed to the polarization of the hydrogen molecules by the high electric field generated near the surface of the charged fullerene. At full hydrogen coverage, these charged fullerenes can gain storage capacities of up to approximately 8.0 wt %. We also find that, contrary to intuitive expectation, fullerenes containing encapsulated metal atoms only exhibit negligible enhancement in the hydrogen binding strength, because the charge donated by the metal atoms is primarily confined inside the fullerene cages. These predictions may prove to be instrumental in searching for a new class of high-capacity hydrogen storage media.

Phase diagram of molecular oxygen adsorption on the (111) platinum surface

Molecular oxygen adsorption on the Pt(111) surface is studied based on ab initio computations and thermodynamics. An O2 adsorption phase diagram is determined. There are two possible chemisorbed molecular states: one at a bridge site and another one at an fcc hollow site. While some population in the bridge sites persists at all coverages, the states coexist through the intermediate coverage phases. The relative coverage of the two species on the surface is determined by the competition between the Pt lattice distortion energy (that results from O2 adsorption) and the O2 repulsion energy. Our results give a reasonable explanation for the seemingly contradictory findings in previous experimental and theoretical work.

Reaction of hydrogen with Ag(111): binding states, minimum energy paths, and kinetics

The interaction of atomic and molecular hydrogen with the Ag(111) surface is studied using periodic density functional total-energy calculations. This paper focuses on the site preference for adsorption, ordered structures, and energy barriers for H diffusion and H recombination. Chemisorbed H atoms are unstable with respect to the H(2) molecule in all adsorption sites below monolayer coverage. The three-hollow sites are energetically the most favorable for H chemisorption. The binding energy of H to the surface decreases slightly up to one monolayer, suggesting a small repulsive H-H interaction on nonadjacent sites. Subsurface and vacancy sites are energetically less favorable for H adsorption than on-top sites. Recombination of chemisorbed H atoms leads to the formation of gas-phase H(2) with no molecular chemisorbed state. Recombination is an exothermic process and occurs on the bridge site with a pronounced energy barrier. This energy barrier is significantly higher than that inferred from experimental temperature-programmed desorption (TPD) studies. However, there is significant permeability of H atoms through the recombination energy barrier at low temperatures, thus increasing the rate constant for H(2) desorption due to quantum tunneling effects, and improving the agreement between experiment and theory.

Density functional theory study of CHx (x=1–3) adsorption on clean and CO precovered Rh(111) surfaces

CH(x) (x=1-3) adsorptions on clean and CO precovered Rh(111) surfaces were studied by density functional theory calculations. It is found that CH(x) (x=1-3) radicals prefer threefold hollow sites on Rh(111) surfaces, and the bond strength between CH(x) and Rh(111) follows the order of CH(3)<CH(2)<CH. A slight attraction between adsorbed CH(x) and CH(x), CH(x) and CO, CO and CO radicals/molecules at coverage of 1/9-1/4 ML is found, and considerable repulsion is built up at coverage higher than 1/4 ML. It is found that CH(x) adsorption results in the reduction of work function due to charge transfer from the adsorbates to the substrate. For CH(3) adsorption, charge accumulation between the substrate and carbon atom is identified, which results in C-H stretch mode softening. For CO precovered surfaces, the mode softening is prevented due to the weakened interaction between CH(x) and substrates from repulsive interaction between adsorbates. Conversely, the overall charge transfer from the CH(x) (x=1-3) to the substrate enhances the charge back donation from the substrate to the empty antibonding states of adsorbed carbon monoxide, which results in the softness of the C-O stretch, respectively. The C 1s surface core-level shifts for CH(x) with and without the presence of CO were calculated, and a negative shift with respect to the carbon in atop CO on Rh(111) is found, and the negative shift follows the order of CO<CH(3)<CH(2)<CH. The results are analyzed in details by difference of charge density and projected density of states.

Phonons in the ordered c(2 × 2) phases of Na and Li on Al(001)

The vibrational properties of the Al(001)-c(2 × 2)-Na (Li) ordered phases formed by alkali atoms (Na and Li) on the Al(001) surface at low and room temperatures are presented. The equilibrium structural characteristics, phonon dispersions and polarization of vibrational modes as well as the local density of phonon states are calculated using the embedded-atom method. The obtained structural parameters are in close agreement with experimental data.

A Cu/Pt Near-Surface Alloy for Water−Gas Shift Catalysis

The primary route to hydrogen production from fossil fuels involves the water-gas shift (WGS) reaction, and an improvement in the efficiency of WGS catalysts could therefore lead to a major leap forward in the realization of hydrogen economy. On the basis of a combination of high-resolution scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT) calculations, we suggest the existence of a new thermodynamically stable Cu/Pt near-surface alloy (NSA). Temperature-programmed desorption and DFT reveal that this Cu/Pt NSA binds CO significantly more weakly than does Pt alone, thereby implying a considerable reduction in the potential for CO poisoning of the Cu/Pt NSA surface as compared to that of pure Pt. In addition, DFT calculations show that this Cu/Pt NSA is able to activate H2O easily, which is the rate-determining step for the WGS on several metal surfaces, and, at the same time, to bind the products of that reaction and formate intermediates rather weakly, thus avoiding possible poisoning of the catalyst surface. The Cu/Pt NSA is thus a promising candidate for an improved WGS catalyst.

Electronic excitations of glycine, alanine, and cysteine conformers from first-principles calculations

The electronic and optical properties are studied for three conformers of amino acid molecules using gradient-corrected (spin-) density functional theory within a projector-augmented wave scheme and the supercell method. We investigate single-particle excitations such as ionization energies and electron affinities as well as pair excitations. By comparing eigenvalues resulting from several local and nonlocal energy functionals, the influence of treatment of exchange and correlation is demonstrated. The excitations are described within the Delta-self-consistent field method with an occupation number constraint to obtain excitation energies and Stokes shifts. The results are used to also discuss the optical absorption properties. In contrast to the lowest single- and two-particle excitation energies, remarkable changes are found in absorption spectra in dependence on the conformation of the molecule geometry.

Structural and electronic properties of Al12X+ (X=C, Si, Ge, Sn, and Pb) clusters

Using the first-principles method with the generalized gradient approximation, the authors have studied the structural and electronic properties of Al(12)X(+) (X=C, Si, Ge, Sn, and Pb) clusters in detail. The ground state of Al(12)C(+) is a low symmetry C(s) structure instead of an icosahedron. However, the Si, Ge, Sn, and Pb atom doped cationic clusters favor icosahedral structures. The ground states for Al(12)Si(+) and Al(12)Ge(+) are icosahedra, while the C(5nu) structures optimized from an icosahedron with a vertex capped by a tetravalent atom have the highest binding energy for Al(12)Sn(+) and Al(12)Pb(+) clusters. The I(h) structure and the C(5nu) structure are almost degenerate for Al(12)Ge(+), whose binding energy difference is only 0.03 eV. The electronic properties are altered much by removing an electron from the neutral cluster. The binding strength of a valence electron is enhanced, while the binding energy of the cluster is reduced much. Due to the open electronic shell, the band gaps between the highest occupied molecular orbital and the lowest unoccupied molecular orbital are approximately 0.3 eV for the studied cationic clusters.

First-principles study of K and Cs adsorbed on Pd(111)

The adsorptions of K and Cs on Pd(111) were studied by the density functional calculations within the generalized gradient approximation. The site preference, bonding character, work function, and electron structure of the system were analyzed. For K and Cs adsorption, the hcp hollow site was found to be preferred for all the coverages investigated. The calculated adsorption geometries for (2 x 2) and (square root 3 x square root 3)R30 degrees phases are both in reasonable agreement with the observed results. The decrease of the work function upon the adsorption of K and Cs can be attributed to a dipole moment associated with the polarized adsorbate atom, which is characterized by depletion of the electron charge in the alkali metal layer and a charge accumulation in the interface region. Our results indicate that the bonding of alkali metal with the Pd(111) surface has a mixed ionic and metallic bond character at low coverage and a metallic bond of covalent character at high coverage.

Two-dimensional Ir cluster lattice on a graphene moiré on Ir(111)

Lattices of Ir clusters have been grown by vapor phase deposition on graphene moirés on Ir(111). The clusters are highly ordered, and spatially and thermally stable below 500 K. Their narrow size distribution is tunable from 4 to about 130 atoms. A model for cluster binding to the graphene is presented based on scanning tunneling microscopy and density functional theory. The proposed binding mechanism suggests that similar cluster lattices might be grown of materials other than Ir.

Surface Dipoles and Work Functions of Alkylthiolates and Fluorinated Alkylthiolates on Au(111)

We study the dipole formation at the surface formed by -CH(3) and -CF(3) terminated short-chain alkylthiolate monolayers on Au(111). In particular, we monitor the change in work function upon chemisorption using density functional theory calculations. We separate the surface dipole into two contributions, resulting from the gold-adsorbate interaction and the intrinsic dipole of the adsorbate layer, respectively. The two contributions turn out to be approximately additive. Adsorbate dipoles are defined by calculating dipole densities of free-standing molecular monolayers. The gold-adsorbate interaction is, to a good degree, determined by the Au-S bond only. This bond is nearly apolar and its contribution to the surface dipole is relatively small. The surface dipole of the self-assembled monolayer is then dominated by the intrinsic dipole of the thiolate molecules. Alkylthiolates increase the work function of Au(111), whereas fluorinated alkylthiolates decrease it.

Structures of dense glycine and alanine adlayers on chiral Cu(3,1,17) surfaces

Density Functional Theory calculations have been used to predict the structures of dense glycine and alanine adlayers on Cu(3,1,17)(S). Facets of this chiral Cu surface result from adsorbate-induced surface reconstruction when glycine or alanine are adsorbed and annealed on Cu(100). We have calculated the surface energy changes associated with this surface reconstruction. Our results allow the enantiospecificity of this reconstruction following adsorption of enantiopure or racemic alanine on Cu(100) to be discussed. The overall stability of glycine and alanine adlayers on Cu(3,1,17)(S) arises from an interplay between the formation of chemical bonds with the Cu surface, deformations in the adsorbed molecules during adsorption, and intermolecular hydrogen bonds within the adlayer; none of these factors individually dominates.

Predicting Catalysis: Understanding Ammonia Synthesis from First-Principles Calculations

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.

Density Functional Theory Study of Ligand Binding on CdSe (0001), (0001̄), and (112̄0) Single Crystal Relaxed and Reconstructed Surfaces: Implications for Nanocrystalline Growth

To gain a better understanding of the influence of ligand-surface interactions on nanocrystalline growth, periodic density functional theory calculations were employed in the study of the binding of organic ligands on the relaxed nonpolar (1120) and polar Se terminated (0001) surfaces and the relaxed and vacancy and adatom reconstructed Cd terminated (0001) surface. We examined chemisorption properties of phosphine, amine, phosphine oxide, carboxylic acid, and phosphinic acid model ligands, including preferred binding sites and geometries, vibrational frequencies, and binding energetics, and compared findings to intrinsic growth via addition of CdSe molecules or Cd and Se atoms. Our results indicate that binding of the ligands is preferred in the electron-poor 1-fold sites on all surfaces, with secondary coordination of the acidic ligands through the hydroxyl hydrogen to the electron-rich surface sites. In general ligand adsorption directly obstructs binding sites for growth species on the (1120) surface and only indirectly on the two polar surfaces. The order of binding affinities on the (1120) and (0001) surfaces is PH(3) < OPH(3) approximately HCOOH < NH(3) < OPH(2)OH and that on the (0001) surface is OPH(3) approximately HCOOH < OPH(2)OH < NH(3) < PH(3). Our findings corroborate the experimental observation that incorporation of the nonbulky phosphinic acid-type ligands with high affinity and high selectivity for both the (1120) and (0001) surfaces strongly enhances unidirectional growth on the (0001) surface, while incorporation of either bulky ligands or ligands with moderate affinity does not. Higher affinity of all traditionally used ligands for the (1120) surface compared to the (0001) surface also suggests that new ligands should be engineered to achieve the synthesis of novel shapes that require preferential growth on the (1120) surface.

Density functional theory study of water dissociation in a double water bilayer with or without coadsorption of CO on Pt(111)

Using density functional theory, we investigate the structure of the double water bilayer with or without coadsorption of CO on Pt(111). The double water bilayer consists of two bilayers. Each bilayer is buckled with every second water molecule being closer to the surface than every other water molecule. CO is found to adsorb most strongly when substituting in the first bilayer, the water molecule closest to the surface. Dissociation of H2O in the water bilayer (with or without CO) is further considered. A great number of pathways for the dissociation are studied. These include homolytic pathways where both dissociation products end up adsorbed on the Pt surface and heterolytic pathways where only the OH is adsorbed, while a proton is transferred to the water adlayers. We find that the heterolytic dissociation pathways are energetically more favorable than the homolytic ones, yet they are all rather endothermic. The most favorable pathways found have reaction energies of 0.60 and 0.52 eV without and with CO present. The corresponding activation energies are 0.99 and 0.53 eV, respectively.

The energetics and structure of rutile TiO(2)(110)

Density functional theory and a pseudopotential plane wave method are applied to study electronic and structural properties of the defect-free TiO(2)(110) surface. The variations of the surface energy, work function, and atomic displacements are examined for partially and fully relaxed slabs modelling the rutile (110) surface, and consisting of up to 33 atomic layers. Relatively small relaxations of atomic positions in the outermost layers have a strong influence on the calculated surface energies and work functions. The effect of nonequivalence of the odd-even layer terminations is explored. A simple method is proposed which allows one to estimate accurate surface energies for relaxed systems from calculations for partially relaxed slabs.

Density Functional Theory and DFT+U Study of Transition Metal Porphines Adsorbed on Au(111) Surfaces and Effects of Applied Electric Fields

We apply density functional theory (DFT) and the DFT+U technique to study the adsorption of transition metal porphine molecules on atomistically flat Au(111) surfaces. DFT calculations using the Perdew-Burke-Ernzerhof exchange correlation functional correctly predict the palladium porphine (PdP) low-spin ground state. PdP is found to adsorb preferentially on gold in a flat geometry, not in an edgewise geometry, in qualitative agreement with experiments on substituted porphyrins. It exhibits no covalent bonding to Au(111), and the binding energy is a small fraction of an electronvolt. The DFT+U technique, parametrized to B3LYP-predicted spin state ordering of the Mn d-electrons, is found to be crucial for reproducing the correct magnetic moment and geometry of the isolated manganese porphine (MnP) molecule. Adsorption of Mn(II)P on Au(111) substantially alters the Mn ion spin state. Its interaction with the gold substrate is stronger and more site-specific than that of PdP. The binding can be partially reversed by applying an electric potential, which leads to significant changes in the electronic and magnetic properties of adsorbed MnP and approximately 0.1 A changes in the Mn-nitrogen distances within the porphine macrocycle. We conjecture that this DFT+U approach may be a useful general method for modeling first-row transition metal ion complexes in a condensed-matter setting.

First-principles studies of chiral step reconstructions of Cu(100) by adsorbed glycine and alanine

Adsorption of amino acids on Cu(100) is known experimentally to induce surface reconstructions featuring intrinsically chiral Cu(3,1,17) facets, but no information about the geometry of the molecules on these chiral facets is available. We present density-functional theory calculations for the structure of glycine and alanine at moderate coverages on Cu(3,1,17). As might be expected, molecules prefer to bind at the step edges on this surface rather than on the surface's (100)-oriented terraces. The adsorption of enantiopure alanine on Cu(3,1,17) is predicted to be weakly enantiospecific, with S-alanine being more stable on Cu(3,1,17)(S) than R-alanine. By comparing the surface energies of Cu(100) and Cu(3,1,17) in the presence of adsorbed glycine or alanine, our calculations provide insight into the driving force for chiral reconstructions of Cu(100) by amino acids.