Electronic structure calculations with GPAW: a real-space implementation of the projector augmented-wave method

Modelling the lattice dynamics in Si(x)Ge(1-x) alloys

The development of simplified models for the simulation of thermodynamic and thermal transport properties in random alloys is of great importance. In this paper we show how a simple second nearest neighbour model can reliably capture the lattice dynamics of Si(x)Ge(1-x) alloys. The model parameters are extracted from DFT-calculated force constant matrices for pure Si, pure Ge and the Si0.5Ge0.5 ordered alloy. We extract the nearest neighbour contributions directly from density functional theory, whereas effective interactions are obtained for the second nearest neighbour contributions. We demonstrate how the thermal properties, including the expansion coefficient, can be reliably reproduced and that the model is transferable to random Si(x)Ge(1-x) alloys.

Interlayer carbon bond formation induced by hydrogen adsorption in few-layer supported graphene

We report on the hydrogen adsorption induced phase transition of a few layer graphene (1 to 4 layers) to a diamondlike structure on Pt(111) based on core level x-ray spectroscopy, temperature programed desorption, infrared spectroscopy, and density functional theory total energy calculations. The surface adsorption of hydrogen induces a hybridization change of carbon from the sp2 to the sp3 bond symmetry, which propagates through the graphene layers, resulting in interlayer carbon bond formation. The structure is stabilized through the termination of interfacial sp3 carbon atoms by the substrate. The structural transformation occurs as a consequence of high adsorption energy.

Layer-resolved study of Mg atom incorporation at the MgO/Ag(001) buried interface

By combining x-ray excited Auger electron diffraction experiments and multiple scattering calculations we reveal a layer-resolved shift for the Mg KL23L23 Auger transition in MgO ultrathin films (4-6 Å) on Ag(001). This resolution is exploited to demonstrate the possibility of controlling Mg atom incorporation at the MgO/Ag(001) interface by exposing the MgO films to a Mg flux. A substantial reduction of the MgO/Ag(001) work function is observed during the exposition phase and reflects both band-offset variations at the interface and band bending effects in the oxide film.

DFT based study of transition metal nano-clusters for electrochemical NH3 production

Theoretical studies of the possibility of producing ammonia electrochemically at ambient temperature and pressure without direct N2 dissociation are presented. Density functional theory calculations were used in combination with the computational standard hydrogen electrode to calculate the free energy profile for the reduction of N2 admolecules and N adatoms on transition metal nanoclusters in contact with an acidic electrolyte. This work has established linear scaling relations for the dissociative reaction intermediates NH, NH2, and NH3. In addition, linear scaling relations for the associative reaction intermediates N2H, N2H2, and N2H3 have been determined. Furthermore, correlations between the adsorption energies of N, N2, and H have been established. These scaling relations and the free energy corrections are used to establish volcanoes describing the onset potential for electrochemical ammonia production and hence describe the potential determining steps for the electrochemical ammonia production. The competing hydrogen evolution reaction has also been analyzed for comparison.

Environmental tight-binding modeling of nickel and cobalt clusters

Tight-binding models derived from density functional theory potentially provide a systematic approach to the development of accurate and transferable models of multicomponent systems. We introduce a systematic methodology for environmental tight binding in which both the overlap and environmental contributions to the electronic structure are included. The parameters of the model are determined directly from ab initio considerations, thus providing a formal conceptual link to density functional approaches. In order to test the validity of the approach, the model is applied to small clusters of Ni and Co, whose electronic structure is largely determined by the interplay of tightly bound d-valent states and the disperse s-states. We numerically illustrate that it is essential to include environmental contributions in the tight-binding approach in order to reliably reproduce the electronic structure of such clusters.

Singular characteristics and unique chemical bond activation mechanisms of photocatalytic reactions on plasmonic nanostructures

The field of heterogeneous photocatalysis has almost exclusively focused on semiconductor photocatalysts. Herein, we show that plasmonic metallic nanostructures represent a new family of photocatalysts. We demonstrate that these photocatalysts exhibit fundamentally different behaviour compared with semiconductors. First, we show that photocatalytic reaction rates on excited plasmonic metallic nanostructures exhibit a super-linear power law dependence on light intensity (rate ∝ intensity(n), with n > 1), at significantly lower intensity than required for super-linear behaviour on extended metal surfaces. We also demonstrate that, in sharp contrast to semiconductor photocatalysts, photocatalytic quantum efficiencies on plasmonic metallic nanostructures increase with light intensity and operating temperature. These unique characteristics of plasmonic metallic nanostructures suggest that this new family of photocatalysts could prove useful for many heterogeneous catalytic processes that cannot be activated using conventional thermal processes on metals or photocatalytic processes on semiconductors.

Graphene coatings: probing the limits of the one atom thick protection layer

The limitations of graphene as an effective corrosion-inhibiting coating on metal surfaces, here exemplified by the hex-reconstructed Pt(100) surface, are probed by scanning tunneling microscopy measurements and density functional theory calculations. While exposure of small molecules directly onto the Pt(100) surface will lift the reconstruction, a single graphene layer is observed to act as an effective coating, protecting the reactive surface from O(2) exposure and thus preserving the reconstruction underneath the graphene layer in O(2) pressures as high as 10(-4) mbar. A similar protective effect against CO is observed at CO pressures below 10(-6) mbar. However, at higher pressures CO is observed to intercalate under the graphene coating layer, thus lifting the reconstruction. The limitations of the coating effect are further tested by exposure to hot atomic hydrogen. While the coating can withstand these extreme conditions for a limited amount of time, after substantial exposure, the Pt(100) reconstruction is lifted. Annealing experiments and density functional theory calculations demonstrate that the basal plane of the graphene stays intact and point to a graphene-mediated mechanism for the H-induced lifting of the reconstruction.

Rationale for switching to nonlocal functionals in density functional theory

Density functional theory (DFT) has been steadily improving over the past few decades, becoming the standard tool for electronic structure calculations. The early local functionals (LDA) were eventually replaced by more accurate semilocal functionals (GGA) which are in use today. A major persisting drawback is the lack of the nonlocal correlation which is at the core of dispersive (van der Waals) forces, so that a large and important class of systems remains outside the scope of DFT. The vdW-DF correlation functional of Langreth and Lundqvist, published in 2004, was the first nonlocal functional which could be easily implemented. Beyond expectations, the nonlocal functional has brought significant improvement to systems that were believed not to be sensitive to nonlocal correlations. In this paper, we use the example of graphene nanodomes growing on the Ir(111) surface, where with an increase of the size of the graphene islands the character of the bonding changes from strong chemisorption towards almost pure physisorption. We demonstrate how the seamless character of the vdW-DF functionals makes it possible to treat all regimes self-consistently, proving to be a systematic and consistent improvement of DFT regardless of the nature of bonding. We also discuss the typical surface science example of CO adsorption on (111) surfaces of metals, which shows that the nonlocal correlation may also be crucial for strongly chemisorbed systems. We briefly discuss open questions, in particular the choice of the most appropriate exchange part of the functional. As the vdW-DF begins to appear implemented self-consistently in a number of popular DFT codes, with numerical costs close to the GGA calculations, we draw the attention of the DFT community to the advantages and benefits of the adoption of this new class of functionals.

Desorption of n-alkanes from graphene: a van der Waals density functional study

A recent study of temperature-programmed desorption (TPD) measurements of small linear alkane molecules (n-alkanes, with formula C(N)H(2N+2)) from C(0001) deposited on Pt(111) shows a linear relationship of the desorption energy with increasing n-alkane chain length N. We here present a van der Waals density functional study of the desorption barrier energy of the ten smallest n-alkanes (of carbon chain length N = 1-10) from graphene. We find linear scaling with N, including a non-zero intercept with the energy axis, i.e. an offset at the extrapolation to N = 0. This calculated offset is quantitatively similar to the results of the TPD measurements. From further calculations of the polyethylene polymer we offer a suggestion for the origin of the offset.

Nonequilibrium thermodynamics of interacting tunneling transport: variational grand potential, density functional formulation and nature of steady-state forces

The standard formulation of tunneling transport rests on an open-boundary modeling. There, conserving approximations to nonequilibrium Green function or quantum statistical mechanics provide consistent but computational costly approaches; alternatively, the use of density-dependent ballistic-transport calculations (e.g., Lang 1995 Phys. Rev. B 52 5335), here denoted 'DBT', provides computationally efficient (approximate) atomistic characterizations of the electron behavior but has until now lacked a formal justification. This paper presents an exact, variational nonequilibrium thermodynamic theory for fully interacting tunneling and provides a rigorous foundation for frozen-nuclei DBT calculations as a lowest-order approximation to an exact nonequilibrium thermodynamic density functional evaluation. The theory starts from the complete electron nonequilibrium quantum statistical mechanics and I identify the operator for the nonequilibrium Gibbs free energy which, generally, must be treated as an implicit solution of the fully interacting many-body dynamics. I demonstrate a minimal property of a functional for the nonequilibrium thermodynamic grand potential which thus uniquely identifies the solution as the exact nonequilibrium density matrix. I also show that the uniqueness-of-density proof from a closely related Lippmann-Schwinger collision density functional theory (Hyldgaard 2008 Phys. Rev. B 78 165109) makes it possible to express the variational nonequilibrium thermodynamic description as a single-particle formulation based on universal electron-density functionals; the full nonequilibrium single-particle formulation improves the DBT method, for example, by a more refined account of Gibbs free energy effects. I illustrate a formal evaluation of the zero-temperature thermodynamic grand potential value which I find is closely related to the variation in the scattering phase shifts and hence to Friedel density oscillations. This paper also discusses the difference between the here-presented exact thermodynamic forces and the often-used electrostatic forces. Finally the paper documents an inherent adiabatic nature of the thermodynamic forces and observes that these are suited for a nonequilibrium implementation of the Born-Oppenheimer approximation.

Structural and theoretical basis for ligand exchange on thiolate monolayer protected gold nanoclusters

Ligand exchange reactions are widely used for imparting new functionality on or integrating nanoparticles into devices. Thiolate-for-thiolate ligand exchange in monolayer protected gold nanoclusters has been used for over a decade; however, a firm structural basis of this reaction has been lacking. Herein, we present the first single-crystal X-ray structure of a partially exchanged Au(102)(p-MBA)(40)(p-BBT)(4) (p-MBA = para-mercaptobenzoic acid, p-BBT = para-bromobenzene thiol) with p-BBT as the incoming ligand. The crystal structure shows that 2 of the 22 symmetry-unique p-MBA ligand sites are partially exchanged to p-BBT under the initial fast kinetics in a 5 min timescale exchange reaction. Each of these ligand-binding sites is bonded to a different solvent-exposed Au atom, suggesting an associative mechanism for the initial ligand exchange. Density functional theory calculations modeling both thiol and thiolate incoming ligands postulate a mechanistic pathway for thiol-based ligand exchange. The discrete modification of a small set of ligand binding sites suggests Au(102)(p-MBA)(44) as a powerful platform for surface chemical engineering.

The halogen analogs of thiolated gold nanoclusters

Is it possible to replace all the thiolates in a thiolated gold nanocluster with halogens while still maintaining the geometry and the electronic structure? In this work, we show from density functional theory that such halogen analogs of thiolated gold nanoclusters are highly likely. Using Au(25)X(18)(-) as an example, where X = F, Cl, Br, or I replaces -SR, we find that Au(25)Cl(18)(-) demonstrates a high similarity to Au(25)(SR)(18)(-) by showing Au-Cl distances, Cl-Au-Cl angles, band gap, and frontier orbitals similar to those in Au(25)(SR)(18)(-). DFT-based global minimization also indicates the energetic preference of staple formation for the Au(25)Cl(18)(-) cluster. The similarity between Au(m)(SR)(n) and Au(m)X(n) could be exploited to make viable Au(m)X(n) clusters and to predict structures for Au(m)(SR)(n).

One-pot synthesis and characterization of subnanometre-size benzotriazolate protected copper clusters

A simple one-pot method for the preparation of subnanometre-size benzotriazolate (BTA) protected copper clusters, Cu(n)BTA(m), is reported. The clusters were analyzed by optical and infrared spectroscopy, mass spectrometry and transmission electron microscopy together with computational methods. We suggest a structural motif where the copper core of the Cu(n)BTA(m) clusters is protected by BTA-Cu(i)-BTA units.

Packing defects into ordered structures: strands on TiO2

We have studied vicinal TiO2(110) surfaces by high-resolution scanning tunneling microscopy and density functional theory calculations. On TiO2 surfaces characterized by a high density of <111> steps, scanning tunneling microscopy reveals a high density of oxygen-deficient strandlike adstructures. With the help of density functional theory calculations we develop a complete structural model for the entire strand and demonstrate these adstructures to be more stable than an equivalent amount of bulk defects such as Ti interstitials. We argue that strands can form particularly easy on stepped surfaces because building material is available at step sites. The strands on TiO2(110) represent point defects that are densely packed into ordered adstructures.

Systematic study of Au6 to Au12 gold clusters on MgO(100) F centers using density-functional theory

We present an optimized genetic algorithm used in conjunction with density-functional theory in the search for stable gold clusters and O2 adsorption ensembles in F centers at MgO(100). For Au8 the method recovers known structures and identifies several more stable ones. When O2 adsorption is investigated, the genetic algorithm is used to imitate structural fluxionality, increasing the O2 bond strength by up to 1 eV. Extending the method to Au(6,10,12), strong O2 adsorption configurations are found for all sizes. However, the effect of fluxionality appears to wear off with increasing cluster size.

Promoter effect of BaO on CO oxidation on PdO surfaces

The effect of bulk BaO promoter on CO oxidation activity of palladium oxide phase was studied by density functional calculations. A series of BaO(100) supported Pd(x)O(y) thin layer models were constructed, and energy profiles for CO oxidation on the films were calculated and compared with corresponding profiles for the most stable PdO bulk surfaces PdO(100) and PdO(101). The most stable of the thin films typically exhibit the same PdO(100) and PdO(101) surface planes; the PdO(100) dominates already with double layer thickness. The supporting promoter improves the CO oxidation activity of the Pd(x)O(y) phase via a direct electronic effect and introduced structural strain and corrugation. Changes in CO adsorption strength are reflected in oxidation energy barriers, and the promoting effect of even 0.3 eV can be seen locally. Easier oxygen vacancy formation may partially facilitate the reaction.

A self-consistent DFT + DMFT scheme in the projector augmented wave method: applications to cerium, Ce2O3 and Pu2O3 with the Hubbard I solver and comparison to DFT + U

An implementation of full self-consistency over the electronic density in the DFT + DMFT framework on the basis of a plane wave–projector augmented wave (PAW) DFT code is presented. It allows for an accurate calculation of the total energy in DFT + DMFT within a plane wave approach. In contrast to frameworks based on the maximally localized Wannier function, the method is easily applied to f electron systems, such as cerium, cerium oxide (Ce2O3) and plutonium oxide (Pu2O3). In order to have a correct and physical calculation of the energy terms, we find that the calculation of the self-consistent density is mandatory. The formalism is general and does not depend on the method used to solve the impurity model. Calculations are carried out within the Hubbard I approximation, which is fast to solve, and gives a good description of strongly correlated insulators. We compare the DFT + DMFT and DFT + U solutions, and underline the qualitative differences of their converged densities. We emphasize that in contrast to DFT + U, DFT + DMFT does not break the spin and orbital symmetry. As a consequence, DFT + DMFT implies, on top of a better physical description of correlated metals and insulators, a reduced occurrence of unphysical metastable solutions in correlated insulators in comparison to DFT + U.

The electronic structure of Ge9[Si(SiMe3)3]3-: a superantiatom complex

We report on the electronic structure of Ge(9)[Si(SiMe(3))(3)](3)(-). Systematic density functional theory analysis of the electronic shell structure of the cluster and its derivatives reveals that the Ge(9)[Si(SiMe(3))(3)](3)(-) and its neutral counterpart have electronic shells that can be explained using the superatom model. The ligand-core interaction of these complexes is distinctly different from previously identified gold, gallium, and aluminium superatom complexes, indicating an electron-donating rather than electron-withdrawing ligand. We modify the electron-counting rule for this case and introduce a simple picture for superatom and superantiatom complexes. Discussions comparing shell models, Zintl clusters, the superhalogen Al(13) and superatom complexes to Ge(9)[Si(SiMe(3))(3)](3)(-) are presented.

Robust conductance of dumbbell molecular junctions with fullerene anchoring groups

The conductance of a molecular wire connected to metallic electrodes is known to be sensitive to the atomic structure of the molecule-metal contact. This contact is to a large extent determined by the anchoring group linking the molecular wire to the metal. It has been found experimentally that a dumbbell construction with C(60) molecules acting as anchors yields more well-defined conductances as compared to the widely used thiol anchoring groups. Here, we use density functional theory to investigate the electronic properties of this dumbbell construction. The conductance is found to be stable against variations in the detailed bonding geometry and in good agreement with the experimental value of G=3×10(-4) G(0). Electron tunneling across the molecular bridge occurs via the lowest unoccupied orbitals of C(60) which are pinned close to the Fermi energy due to partial charge transfer. Our findings support the original motivation to achieve conductance values more stable towards changes in the structure of the molecule-metal contact leading to larger reproducibility in experiments.

Towards quantitative accuracy in first-principles transport calculations: The GW method applied to alkane/gold junctions

The calculation of the electronic conductance of nanoscale junctions from first principles is a long-standing problem in the field of charge transport. Here we demonstrate excellent agreement with experiments for the transport properties of the gold/alkanediamine benchmark system when electron-electron interactions are described by the many-body GW approximation. The conductance follows an exponential length dependence: G(n) = G(c) exp(-βn). The main difference from standard density functional theory (DFT) calculations is a significant reduction of the contact conductance, G(c), due to an improved alignment of the molecular energy levels with the metal Fermi energy. The molecular orbitals involved in the tunneling process comprise states delocalized over the carbon backbone and states localized on the amine end groups. We find that dynamic screening effects renormalize the two types of states in qualitatively different ways when the molecule is inserted in the junction. Consequently, the GW transport results cannot be mimicked by DFT calculations employing a simple scissors operator.

Oxidation of Pt(111) under near-ambient conditions

The oxidation of Pt(111) at near-ambient O2 pressures has been followed in situ using x-ray photoelectron spectroscopy (XPS) and ex situ using x-ray absorption spectroscopy (XAS). Polarization-dependent XAS signatures at the O K edge reveal significant temperature- and pressure-dependent changes of the Pt-O interaction. Oxide growth commences via a PtO-like surface oxide that coexists with chemisorbed oxygen, while an ultrathin α-PtO2 trilayer is identified as the precursor to bulk oxidation. These results have important implications for understanding the chemical state of Pt in catalysis.