Search

Plasmonic materials have optical cross sections that exceed by 10-fold their geometric sizes, making them uniquely suitable to convert light into electrical charges. Harvesting plasmon-generated hot carriers is of interest for the broad fields of photovoltaics and photocatalysis; however, their direct utilization is limited by their ultrafast thermalization in metals. To prolong the lifetime of hot carriers, one can place acceptor materials, such as semiconductors, in direct contact with the plasmonic system. Herein, we report the effect of operating temperature on hot electron generation and transfer to a suitable semiconductor. We found that an increase in the operation temperature improves hot electron harvesting in a plasmonic semiconductor hybrid system, contrasting what is observed on photodriven processes in nonplasmonic systems. The effect appears to be related to an enhancement in hot carrier generation due to phonon coupling. This discovery provides a new strategy for optimization of photodriven energy production and chemical synthesis.

A low-cost template-free solution chemical route to highly porous nanocrystalline sponges of ZnO-EuO1.5 with 0-5 mol % Eu is presented. The process uses Zn- and Eu-acetate-nitrate and triethanolamine as precursors in methanol. After evaporation of the solvent and heating at 200 °C for 3 min, crystalline ZnO:Eu sponges with minor amounts of organic residues were obtained. Heating to 400 °C replaced the organics with carbonate, which in its turn was decomposed at temperatures below 600 °C, forming ZnO:Eu sponges. Samples heated to 200-1000 °C for 3 min were studied with XRD, SEM, TEM, TG, XPS, and IR spectroscopy. The ZnO:Eu crystallite sizes could be tuned from below 10 nm for sponges prepared at 200-500 °C, to over 100 nm range at 900 °C, without sintering of the overall microstructure. XRD showed the presence of hexagonal ZnO:Eu (or at 700-1000 °C, ZnO:Eu and cubic Eu2O3) as the only phases present. The ZnO:Eu had slightly larger unit cell dimensions than the literature value of ZnO for samples obtained at 200-600 °C, while the unit cells of samples obtained at higher temperatures were quite close to the value of undoped ZnO. XPS showed that Eu was mainly in its 3+ state and well-distributed within the sponges but segregated at the ZnO sponge surface upon heating at 700-1000 °C, in accordance with XRD studies showing Eu2O3 formation.

Density functional theory calculations (DFT), coupled with microkinetic modelling, have been used to simulate Temperature Programmed Desorption (TPD) experiments for calcined ceria nanopowders with the aim to gain insight into the chemistry governing their high redox activity. Our simulations consider two main nanoparticle models. One is a perfect ceria octahedron supercharged with adsorbed oxygen molecules turned into superoxide ions, as has previously been used to explain the enhanced oxygen storage capacity (OSC) in nanoceria. The other model is a variant where we have introduced oxygen vacancies under ridge Ce ions, thereby reducing their coordination numbers to five. The results from our microkinetic modelling suggest that including such five-coordinated Ce adsorption sites results in a TPD spectrum that better matches the experimental counterpart in terms of both peak position and width. In addition, this new structural model allows for the co-existence of Ce 3 + ions, superoxide ions and O2 molecules, as seen in experiments in the literature.

Adsorbate induced work function modification of Ni have been investigated by means of first-principles calculations. More specifically, the adsorption of Li, Na, Si, Zr, Pd, Pt, or Sn at various coverages on Ni low-index surface models have been considered. In the case of Sn, a more thorough investigation was performed comparing the adsorption as an overlayer structure with the case of surface alloy formation. Our calculations suggest that the most stable Sn@Ni configuration corresponds to a surface alloy, and here the Ni(100)c(2 × 2)-Sn, Ni(110)c(2 × 2)-Sn, and Ni(111)()R30-Sn surface alloys were found to display similar stability. Concerning the induced work function change, a different behaviour as a function of coverage was observed depending on the nature of the Sn@Ni surface model. Both overlayer adsorption and surface alloying were found to induce a work function decrease already at relatively low coverages ( 0.05 atom Å −2), regardless of the underlying surface orientation. However, while the work function obtained for stable surface alloys was found to monotonously decrease as the coverage increases, the work function for the stable overlayer structures goes through a minimum. For all investigated surface modifications, the change in work function was found to be consistent with the orientation of the charge transfer at the adsorbate–surface interface. The computed data in this work may serve as handles for experimental endeavours aiming to optimize properties of active materials through atomic-scale surface engineering.

Silyl groups are included in a wide range of electrolyte additives to enhance the performance of state-of-the-art Li-ion batteries. A recognized representative thereof is tris-(trimethylsilyl)phosphate (TMSPa) which, along with the similarly structured phosphite, has been at the center of numerous electrolyte studies. Even though the silyl group has already been widely reported to be specifically reactive towards fluorides, herein, a reactivity towards several Lewis bases typically found in Li-ion cells is postulated and investigated with the aim to establish a more simplified and generally applicable reaction mechanism thereof. Both gaseous and electrolyte soluble reactants and products are monitored by combining nuclear magnetic resonance and injection cell-coupled mass spectrometry. Experimental observations are supported by computational models. The results clearly demonstrate that the silyl groups react with water, hydroxide, and methoxide and thereby detach in a stepwise fashion from the central phosphate in TMSPa. Intermolecular interaction between TMSPa and the reactants likely facilitates dissolution and lowers the free energy of reaction. Lewis bases are well known to trigger side reactions involving both the Li-ion electrode and electrolyte. By effectively scavenging these, the silyl group can be explained to lower cell impedance and prolong the lifetime of modern Li-ion batteries.

A computational protocol is developed for efficient studies of partially reduced redox-active oxides using the self-consistent charge density functional tight-binding method. The protocol is demonstrated for ceria, which is a prototypical reducible oxide material. The underlying idea is to achieve a consistent (and harmonized) set of Slater–Koster (SK) tables with connected repulsive potentials that enable switching on and off the in-valence description of the Ce 4 f states without serious loss of accuracy in structure and energetics. The implicit treatment of the Ce 4 f states, with the use of f-in-core SK-tables, is found to lead to a significant decrease in computational time. More importantly, it allows for explicit control of the oxidation states of individual Ce atoms. This makes it possible to “freeze” the electronic configuration, thereby allowing the exploration of the energetics for various meta-stable configurations. We anticipate that the outlined strategy can help to shed light on the interplay between the size, shape, and redox activity for nanoceria and other related materials.

The detailed structure of ZnO doped with 5 at.% (metal) of the large aliovalent Eu3+-ions was investigated using EXAFS to describe the local Eu and Zn coordination. The microstructure, crystalline phases, contents and ZnO unit-cell parameters for the ZnO:5%Eu sponges synthesised at 200 to 900 oC were obtained by XRD, SEM, and TEM analysis. XRD showed peaks solely of h-ZnO for the 600 oC sample, while heating at 700 oC and higher caused phase separation into h-ZnO:Eu and c-Eu2O3. XRD showed a close to zero increase in ZnO unit cell-volume of ca. 0.4 vol%, compared to un-doped ZnO for the non-phase separated, clean oxide made at 600 oC. The Zn EXAFS data showed an almost intact local ZnO structure. The Eu EXAFS showed an unusually low coordination number (CN) of ca. 5 for the 200-600 oC samples, while the CN increased for higher temperatures, in concert with the formation of c-Eu2O3. 23 DFT-generated theoretical ZnO structures containing Eu-clusters built from 1 to 4 Eu3+-Zn2+-vacancy- Eu3+ pairs were compared with the experimental data. The lowest formation energies and ZnO unit-cell volume increase versus un-doped ZnO (0.6-0.7 vol%), were obtained when combining two or four Eu3+-Zn2+-vacancy- Eu3+ pairs into Eu4 and Eu8 clusters showing an average Eu CN of ca. 5. These theoretically determined lowest energy structures were all in good agreement with the experimental results obtained by EXAFS and XRD. Photoluminescence excitation and emission spectra performed on the ZnO:5at%Eu sponges obtained at various temperatures, showed strong quenching of the characteristic Eu3+ transitions for samples obtained at 600 and 800 °C, most likely due to changes in the ZnO defect states, which are crucial for Eu3+ excitation, and due to self-quenching upon Eu clustering and c- Eu2O3 phase separation. Thus, the optical data further supported Eu clustering found by EXAFS, DFT and XRD techniques, corroborating structure-property relationships in these materials. Overall, as far as we can find, the findings reported herein point to a doping structure very different from those previously proposed in the literature. It demonstrates that the semiconductor ZnO can host molecular-sized clusters of metal-oxides, very dissimilar to ZnO.

Poly(ethylene oxide) (PEO) is the most widely used compound as a solid-state (solvent-free) polymer electrolyte for Li batteries, mainly due to its low glass transition temperature (Tg) and ability...

Nanocrystalline ZnO sponges doped with 5 mol% EuO1.5 are obtained by heating metal–salt complex based precursor pastes at 200–900 °C for 3 min. X-ray diffraction, transmission electron microscopy, and extended X-ray absorption fine structure (EXAFS) show that phase separation into ZnO:Eu and c-Eu2O3 takes place upon heating at 700 °C or higher. The unit cell of the clean oxide made at 600 °C shows only ≈0.4% volume increase versus undoped ZnO, and EXAFS shows a ZnO local structure that is little affected by the Eu-doping and an average Eu3+ ion coordination number of ≈5.2. Comparisons of 23 density functional theory-generated structures having differently sized Eu-oxide clusters embedded in ZnO identify three structures with four or eight Eu atoms as the most energetically favorable. These clusters exhibit the smallest volume increase compared to undoped ZnO and Eu coordination numbers of 5.2–5.5, all in excellent agreement with experimental data. ZnO defect states are crucial for efficient Eu3+ excitation, while c-Eu2O3 phase separation results in loss of the characteristic Eu3+ photoluminescence. The formation of molecule-like Eu-oxide clusters, entrapped in ZnO, proposed here, may help in understanding the nature of the unexpected high doping levels of lanthanide ions in ZnO that occur virtually without significant change in ZnO unit cell dimensions.

The Curvature Constrained Splines (CCS) methodology has been used for fitting repulsive potentials to be used in SCC-DFTB calculations. The benefit of using CCS is that the actual fitting of the repulsive potential is performed through quadratic programming on a convex objective function. This guarantees a unique (for strictly convex) and optimum two-body repulsive potential in a single shot, thereby making the parametrization process robust, and with minimal human effort. Furthermore, the constraints in CCS give the user control to tune the shape of the repulsive potential based on prior knowledge about the system in question. Herein, we developed the method further with new constraints and the capability to handle sparse data. We used the method to generate accurate repulsive potentials for bulk Si polymorphs and demonstrate that for a given Slater-Koster table, which reproduces the experimental band structure for bulk Si in its ground state, we are unable to find one single two-body repulsive potential that can accurately describe the various bulk polymorphs of silicon in our training set. We further demonstrate that to increase transferability, the repulsive potential needs to be adjusted to account for changes in the chemical environment, here expressed in the form of a coordination number. By training a near-sighted Atomistic Neural Network potential, which includes many-body effects but still essentially within the first-neighbor shell, we can obtain full transferability for SCC-DFTB in terms of describing the energetics of different Si polymorphs.

We evaluate the qualitative and quantitative accuracy of various flavors of density functionals with and without accounting for dispersion corrections. Our test system is nickel in the form of bulk, surfaces, and nanoparticles for which we compute structural properties, bulk cohesive energies, surface energies, and work functions and compare to experimental data. We find that the inclusion of any dispersion, either by an a posteriori correction or by a self-consistent treatment by explicitly computing the nonlocal correlation contribution to the total energy, has a significant effect on the calculated properties and improves the quantitative comparison to experiments. Besides the quantitative agreement, we also investigate qualitative features by comparing Wulff shapes of metal nanoparticles as obtained using the different density functionals. We find that all tested functionals predict similar Wulff shapes for nickel nanoparticles but still have some small differences. These results show that the relative energies calculated using the semilocal GGA and meta-GGA functionals, with and without dispersion, are quite similar. Our findings can also be generalized to other systems when rationalized in terms of the computed reduced density gradients. We find that the distribution of reduced density gradients in a material is correlated to the steepness of the exchange enhancement factor and propose that this information can be used as a quantitative guide when it comes to picking the most appropriate density functional for specific target systems as well as when it comes to extrapolating DFT data to predict experiments.

The effect of biaxial strain on the band structure of two-dimensional silicon nanosheets (Si NSs) with (111), (110), and (001) exposed surfaces was investigated by means of density functional theory calculations. For all the considered Si NSs, an indirect-to-direct band gap transition occurs as the lateral dimensions of Si NSs increase; that is, increasing lateral biaxial strain from compressive to tensile always enhances the direct band gap characteristics. Further analysis revealed the mechanism of the transition which is caused by preferential shifts of the conduction band edge at a specific k-point because of their bond characteristics. Our results explain a photoluminescence result of the (111) Si NSs [U. Kim et al., ACS Nano 2011, 5, 2176–2181] in terms of the plausible tensile strain imposed in the unoxidized inner layer by surface oxidation.

It has been suggested in the literature that keto‐to‐enol tautomerization plays a vital role for lignin fragmentation under mild conditions. On the other hand, previous modelling has shown that the adsorbed keto form is more stable than enol on the Pd(111) catalyst. The current density functional theory study of lignin model molecules shows that, in the gas‐phase, keto is more stable than enol, but on the Pd surface, we find enol conformers that are at least as stable as keto. This supports the experimental result that the favourable reaction pathway for lignin depolymerization involves keto‐enol tautomerization. An energy decomposition analysis gives insights concerning the origin of the fine energy balance between the keto and enol forms, where the molecule–surface interaction (−7 eV) and the molecular strain energy (+3 eV) are the main contributors to the adsorption energy., It's a fine balance: Experiments have suggested that keto‐enol tautomerism constitutes one important step in lignin's depolymerization reaction under mild conditions. In the gas‐phase, the keto form is more stable than the enol form, but the presented DFT study shows that, on the surface of a Pd catalyst, the enol conformers can be at least as stable as the keto conformers, supporting experiments.

Thin film structures of water on the CeO2(111) surface for coverages between 0.5 and 2.0 water monolayers have been optimized and analyzed using density functional theory (optPBE-vdW functional). We present a new 1.0 ML structure that is both the lowest in energy published and features a hydrogen-bond network extending the surface in one-dimension, contrary to what has been found in the literature, and contrary to what has been expected due to the large bulk ceria cell dimension. The adsorption energies for the monolayer and multilayered water structures agree well with experimental temperature programmed desorption results from the literature, and we discuss the stability window of CeO2(111) surfaces covered with 0.5-2.0 ML of water.

Cesium and methylammonium bismuth iodides (Cs3Bi2I9 and MA3Bi2I9) are new low-toxic and air stable compounds in the perovskite solar cell family with promising characteristics. Here, the electronic structure and the nature of their optical transitions, dielectric constant, and charge carrier properties are assessed for photovoltaic applications with density functional theory (DFT) calculations and experiments. The calculated direct and indirect band gap values for Cs3Bi2I9 (2.17 and 2.0 eV) and MA3Bi2I9 (2.17 and 1.97 eV) are found to be in good agreement with the experimental optical band gaps (2.2, 2.0 eV and 2.4, 2.1 eV for Cs3Bi2I9 and MA3Bi2I9, respectively) estimated for solution-processed films. There is an error cancelation in the DFT calculated band gap similar to that for lead perovskites. However, fully relativistic DFT calculations indicate that the size of the spin orbit coupling (SOC) error cancelation for bismuth perovskite (0.5 eV) is less than for lead perovskite (1 eV), and other factors...

Density functional theory (DFT) calculations, combined with a constrained minima hopping algorithm (global minimum search while preserving the molecular identity), have been performed to investigate important reaction intermediates for the heterogeneously catalyzed β-O-4′ bond cleavage in lignin derivatives. More specifically, we have studied the adsorption properties of a keto tautomer (1-methoxypropan-2-one) and its enol form on a catalytically active Pd(111) surface. In agreement with experiments, we find that for the gas-phase molecules the keto tautomer is the most stable. Interestingly, the enol tautomer has a higher affinity to the Pd catalyst than the keto form, and becomes the most stable molecular form when adsorbed on the catalyst surface. The global minimum complex found on the metal surface corresponds to an enolate structure formed when the enol tautomer chemisorbs onto the surface and donates its π-electrons from the C═C region to two adjacent palladium atoms. The actual formation of a chem...

We have developed a reactive force-field of the ReaxFF type for stoichiometric ceria (CeO2) and partially reduced ceria (CeO2–x). We describe the parametrization procedure and provide results validating the parameters in terms of their ability to accurately describe the oxygen chemistry of the bulk, extended surfaces, surface steps, and nanoparticles of the material. By comparison with our reference electronic structure method (PBE+U), we find that the stoichiometric bulk and surface systems are well reproduced in terms of bulk modulus, lattice parameters, and surface energies. For the surfaces, step energies on the (111) surface are also well described. Upon reduction, the force-field is able to capture the bulk and surface vacancy formation energies (Evac), and in particular, it reproduces the Evac variation with depth from the (110) and (111) surfaces. The force-field is also able to capture the energy hierarchy of differently shaped stoichiometric nanoparticles (tetrahedra, octahedra, and cubes), and ...

We investigate how the redox properties of ceria are affected by the fraction of Fock exchange in screened HSE06-based hybrid density functionals, and we compare with $\mathrm{PBE}+U$ results, and with experiments when available. We find that using $\ensuremath{\sim}15%$ Fock exchange yields a good compromise with respect to structure, electronic structure, and calculated reduction energies, and represents a significant improvement over the $\mathrm{PBE}+U$ results. We also investigate the possibility to use a computationally cheaper HSE06//$\mathrm{PBE}+U$ protocol consisting of structure optimization with $\mathrm{PBE}+U$, a subsequent lattice parameter rescaling step, and, finally, a single-point full hybrid calculation. We find that such a composite computational protocol works very well and yields results in close agreement with those where HSE06 was used also for the structure optimization.

The widespread use of ceria-based materials and the need to design suitable strategies to prepare eco-friendly CeO2 supports for effective catalytic screening induced us to extend our computational multiscale protocol to the modeling of the hybrid organic/oxide interface between prototypical fluorinated linear alkane chains (polyethylene-like oligomers) and low-index ceria surfaces. The combination of quantum chemistry calculations and classical reactive molecular dynamics simulations provides a comprehensive picture of the interface and discloses, at the atomic level, the main causes of typical adsorption modes. The data show that at room temperature a moderate percentage of fluorine atoms (around 25%) can enhance the interaction of the organic chains by anchoring strongly pivotal fluorines to the channels of the underneath ceria (100) surface, whereas an excessive content can remarkably reduce this interaction because of the repulsion between fluorine and the negatively charged oxygen of the surface.

We propose that crystalline ionic hydrates constitute a valuable resource for benchmarking theoretical methods for aqueous ionic systems. Many such structures are known from the experimental literature, and they contain a large variety of water–water and ion–water structural motifs. Here we have collected a data set (CRYSTALWATER50) of 50 structurally unique “in-crystal” water molecules, involved in close to 100 nonequivalent O–H···O hydrogen bonds. A dozen well-known DFT functionals were benchmarked with respect to their ability to describe these experimental structures and their OH vibrational frequencies. We find that the PBE, RPBE-D3, and optPBE-vdW methods give the best H-bond distances and that anharmonic OH frequencies generated from B3LYP//optPBE-vdW energy scans outperform the other methods, i.e., here we performed B3LYP energy scans along the OH stretching coordinate while the rest of the structure was kept fixed at the optPBE-vdW-optimized positions.

This paper concerns thin water films and their hydrogen-bond patterns on ionic surfaces. As far as we are aware, this is the first time H-bond correlations for surface water and hydroxide species are presented in the literature while hydrogen-bond relations in the solid state have been scrutinized for at least five decades. Our data set, which was derived using density functional theory, consists of 116 unique surface OH groups–intact water molecules as well as hydroxides–on MgO(001), CaO(001) and NaCl(001), covering the whole range from strong to weak to no H-bonds. The intact surface water molecules are found to always be redshifted with respect to the gas-phase water OH vibrational frequency, whereas the surface hydroxide groups are either redshifted (OsH) or blueshifted (OHf) compared to the gas-phase OH– frequency. The surface H-bond relations are compared with the traditional relations for bulk crystals. We find that the “ν(OH) vs R(H···O)” correlation curve for surface water does not coincide with ...

The interaction between water molecules and small Cu clusters (up to a size of four atoms) adsorbed on the nonpolar ZnO(10 (1) over bar0) surface has been studied using hybrid density functional th ...

This thesis discusses the chemistry and physics of Cu and H2O on ZnO surfaces, based primarily on results from quantum chemical calculations. The underlying context is heterogeneous catalysis, where Cu/ZnO-mixtures are used in the industrial synthesis of methanol and in the water gas shift reaction. Electron transfer between small Cu clusters and ZnO is central to this thesis, as are the design and use of models that can describe realistic and very large-scale ZnO surface structures while still retaining the electronic nature of the system. Method and model enhancements as well as tests and validations constitute a large part of this thesis.The thesis demonstrates that the charges of small Cu clusters, adsorbed on the non-polar ZnO(10-10) surface, depend on whether the Cu clusters contain an even or odd number of atoms, and whether water is present (water can induce electron transfer from Cu to ZnO). On the polar Zn-terminated ZnO(0001) surface, Cu becomes negatively charged, which causes it to attract positively charged subsurface defects and to wet the ZnO(0001) surface at elevated temperatures.When a Cu cluster on a ZnO surface becomes positively charged, this happens because it donates an electron to the ZnO conduction band. Hence, it is necessary to use a method which describes the ZnO band gap correctly, and we show that a hybrid density functional, which includes a fraction of Hartree-Fock exchange, fulfills this requirement. When the ZnO conduction band becomes populated by electrons from Cu, band-filling occurs, which affects the adsorption energy. The band-filling correction is presented as a means to extrapolate the calculated adsorption energy under periodic boundary conditions to the zero coverage (isolated adsorbate, infinite supercell) limit.A part of this thesis concerns the parameterization of the computationally very efficient SCC-DFTB method (density functional based tight binding with self-consistent charges), in a multi-scale modeling approach. Our findings suggest that the SCC-DFTB method satisfactorily describes the interaction between ZnO surfaces and water, as well as the stabilities of different surface reconstructions (such as triangularly and hexagonally shaped pits) at the polar ZnO(0001) and ZnO(000-1) surfaces.

The adsorption and dissociation of a formic acid molecule (HCOOH) on a partially reduced rutile TiO2-x (110) surface and the subsequent transformations of the adsorbed fragments are studied via qua ...

We have calculated the stabilities of some reactive oxygen species (ROS) in stoichiometric bulk ceria and at the low-index (111) and (110)-surfaces, both in vacuum and in the presence of additional O-2 molecules. We find that the formation of intrinsic ROS, here oxygen superoxides (O-2(-)) and peroxides (O-2(2-)), is always endothermic at vacuum conditions and that the superoxide formation always leads to a higher formation energy than the peroxide formation. In the presence of additional O-2 molecules, intrinsic peroxide formation becomes exothermic at the (110)-surface in conjunction with the formation of extrinsic superoxide ions from adsorbed O-2 molecules. This coexistence of intrinsic and extrinsic ROS species is anticipated to be stable at low temperatures, and can be important for understanding the ROS chemistry for nanoceria used in low-temperature applications. [GRAPHICS] Oxygen-assisted reduction of the ceria(110) surface.

Using hybrid density functional theory, we investigate structural and electronic properties of small Cu-n clusters (with n = 7. For n = 5, both neutral planar and positively charged polyhedral conf ...

We investigate the effects of anion doping with fluorine impurities on the chemistry of the CeO

In this work, a range of van der Waals type density functionals are applied to the H2O/NaCl(001) and H2O/MgO(001) interface systems to explore the effect of an explicit dispersion treatment. The functionals we use are the self-consistent vdW functionals vdW-DF, vdW-DF2, optPBE-vdW, optB88-vdW, optB86b-vdW, and vdW-DF-cx, as well as the dispersion-corrected PBE-TS and PBE-D2 methods; they are all compared with the standard PBE functional. For both NaCl(001) and MgO(001), we find that the dispersion-flavoured functionals stabilize the water-surface interface by approximately 20%-40% compared to the PBE results. For NaCl(001), where the water molecules remain intact for all overlayers, the dominant contribution to the adsorption energy from "density functional theory dispersion" stems from the water-surface interactions rather than the water-water interactions. The optPBE-vdW and vdW-DF-cx functionals yield adsorption energies in good agreement with available experimental values for both NaCl and MgO. To probe the strengths of the perturbations of the adsorbed water molecules, we also calculated water dipole moments and found an increase up to 85% for water at the MgO(001) surface and 70% at the NaCl(001) surface, compared to the gas-phase dipole moment.

The role of hydrogen and other defects in the stabilization of polar oxide interfaces is a matter of significant fundamental and practical interest. Using experimental (scanning tunneling microscopy, x-ray photoelectron spectroscopy) and theoretical (density functional theory) surface science techniques, we find that the polar Zn-terminated ZnO(0001) surface becomes excessively Zn deficient during high-temperature annealing (780 K) in ultrahigh vacuum (UHV). The Zn vacancies align themselves into rows parallel to the $[10\overline{1}0]$ direction, and the remaining surface Zn ions alternately occupy wurtzite (hcp) and zinc-blende (fcc) lattice positions, giving a characteristic ``striped'' $c(\sqrt{12}\ifmmode\times\else\texttimes\fi{}\sqrt{12})\mathrm{R}30$\ifmmode^\circ\else\textdegree\fi{} surface morphology with three types of rows: wurtzite Zn, zinc-blende Zn, and Zn vacancies. Interstitial H plays a central role in such a reconstruction, as it helps to compensate the excessive Zn deficiency. We propose a model in which hydrogen occupies positions in half of the vacancy rows to form hydroxide ions that can participate in hydrogen bonds in the O subsurface layer as a result of the mixed wurtzite/zinc-blende stacking.

The electronic and structural properties of the Ge/GeO2 interface are addressed through a density functional simulation scheme which includes the use of hybrid functionals for achieving accurate band gaps, band offsets and defect levels. The present work discusses the germanium dangling bond levels, the thermodynamics of GeOx, the stability of the oxygen vacancy across Ge/HfO2 interfaces, the atomic structure of GeOx, electron and hole trapping in GeOx, and the band alignment at the Ge/GeO2 interface. (C) 2011 Elsevier B.V. All rights reserved.

An atomistic model of substoichiometric germanium oxide is generated through ab initio molecular dynamics. The resulting structure shows a predominance of threefold coordinated Ge and O atoms. We also generate substoichiometric models through bond-switching Monte-Carlo simulations, which preserve the fourfold Ge and the twofold O coordinations. These differing structures are energetically competitive. Alignment of their electron densities of states to that of GeO"2 reveals that the band-gap reduction is similar for both structures, mainly occurring through a shift of the valence band edge.

Hybrid density functional calculations applied to defect charge transition levels are explored in the attempt to overcome the band-gap problem of semilocal density functionals. Charge transition levels of a large set of point defects calculated with semilocal and hybrid density functionals are found to correspond closely when aligned with respect to the average electrostatic potential. This strongly suggests that the defect levels defined in this way are already accurately described at these levels of theory. In particular, this then also applies to the energy separation between different defect levels, which is directly accessible experimentally. At variance, within the same alignment scheme, the band edges obtained with hybrid functionals are found to undergo significant shifts with respect to those obtained with semilocal functionals. While these shifts systematically give larger band gaps, the agreement with experiment is not always satisfactory when a fixed fraction of exact exchange is admixed. This describes a current limitation of hybrid functional schemes. In the attempt of identifying a viable theoretical description within the class of one-parameter hybrid functionals based on bare exchange, we explore the validity of the empirical procedure which consists in tuning the fraction of nonlocal exchange to a value which gives a theoretical band gap reproducing the experimental one. Comparisons with experiment for band offsets and specific defect levels record very encouraging results. Despite its inherent limitations, such an empirical scheme based on hybrid functionals represents a definite improvement with respect to semilocal functionals.

The electronic and structural properties at the Ge/GeO2 interface are addressed through a density-functional simulation scheme, which includes the use of hybrid functionals for achieving defect levels and band gaps comparable to experiment. The present work discusses the determination of the germanium dangling bond levels, the band alignment in a model structure of the Ge/GeO2 interface, the thermodynamics of GeOx, the stability of the oxygen vacancy across Ge/GeOx/HfO2 stacks, the structural and electronic properties of GeOx, and its ability to trap electrons.

We introduce a theoretical scheme to study defect energy levels and band alignments at semiconductor-oxide interfaces. The scheme relies on hybrid functionals to overcome the band gap underestimation typically achieved with semilocal density functionals. For atomically localized defects, the more accurate description achieved with, hybrid functionals does not lead to significant shifts of the charge transition levels, provided these levels are referred to a common reference potential. This result effectively decouples the shifts of the band edges with respect to the defect levels. We also show that relative shifts of conduction and valence band edges as determined by exact nonlocal exchange. lead to band offsets in excellent agreement with experimental values for several semiconductor-oxide interfaces. The proposed scheme is illustrated through a series of applications, including the dangling bond defects in silicon and germanium, the charge state of the O-2 molecule during silicon oxidation, and the oxygen vacancy in Si-SiO2 HfO2 stacks. (C) 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

We calculate defect levels of dangling bonds in germanium using hybrid density functionals. To validate our approach, we first consider the dangling bond in silicon finding two well-separated defect levels, in excellent correspondence with their experimental location. Application of our scheme to the dangling bond in germanium then yields two very close defect levels lying just above the valence band, consistent with the location of the experimentally determined charge neutrality level. The small correlation energy and the proximity to the valence band edge provide an explanation for the absence of clear defect signatures in electrical and electron spin resonance experiments. © 2008 Elsevier Ltd. All rights reserved.

Models of amorphous (HfO"2)"x (SiO"2)"1"-"x for varying hafnium content x are generated by ab initio molecular dynamics. The structural properties are analyzed in terms of pair distribution functions from which typical bond lengths and coordination numbers are extracted. Electronic band gaps are calculated with a hybrid density functional and are found to decrease in a nonlinear way with increasing x, in accord with experimental observations. Static dielectric constants are evaluated through the application of a finite electric field. The calculated values show a linear dependence on x, supporting recent experimental data.

Experimental characterization is facing severe difficulties in elucidating the nature of defects at interfaces between silicon and amorphous high-@k oxides. This motivates recourse to first principle simulation approaches based on density functional theory. A major issue for these approaches is the choice of the model structure. Indeed, the most trivial choice of parent crystalline structures may fail in describing the amorphous nature of the oxide or the specific bonding arrangements occurring at the interface. A second critical issue relates to the determination of energy levels with respect to band edges of the semiconductor and the oxide, which often represent the only contact with experimental data. An improved description of the electronic structure is required going beyond standard density functional approaches that severely underestimate electronic band gaps. Here, these issues are illustrated by a comparative study of the oxygen vacancy in monoclinic HfO"2, in amorphous HfO"2, and at the Si-HfO"2 interface.

We provide an explanation for the experimental finding of a dramatically enhanced low-temperature oxygen storage capacity for small ceria nanoparticles. At low temperature, small octahedral ceria nanoparticles will be understoichiometric at both oxidizing and reducing conditions without showing explicit oxygen vacancies. Instead, rather than becoming stoichiometric at oxidizing conditions, such particles are stabilized through oxygen adsorption forming superoxo (O2(-)) ions and become in this way supercharged with oxygen. The supercharging effect is size-dependent and largest for small nanoparticles where it gives a direct increase in the oxygen storage capacity and simultaneously provides a source of active oxygen species at low temperatures.

Density functional calculations are reported for NO2 adsorption on BaO(1 0 0) films supported on Pt. The metal is predicted to largely stabilize NO2 adsorption. The phenomenon is most pronounced for a two layer BaO film, but clearly present for films as thick as five oxide layers. From the different contributions to the adsorption energy, the stabilization is found to originate from an increased oxide–metal bonding. The results stress the complex nature of molecular-surface interactions in multi-component systems, and have implications for the understanding of surface processes in heterogeneous catalysis.

The nature of NO2 and nitrite/nitrate pairs adsorbed on BaO has been studied within the density functional theory. Size effects are investigated by comparing results for the BaO molecule, BaO clusters [(BaO)x, x = 4, 6, 9, 12] and BaO(1 0 0). The adsorption energies show weak dependence on size, which is a manifestation of the fast size convergence of the BaO electronic structure and the local character of the NO2 adsorption bond. Nitrite/nitrate pair formation is associated with a significant energy gain. For BaO(1 0 0), the stability of the pair is insensitive to the nitrite–nitrate separation, a finding that demonstrates a surface mediated non-local mechanism of molecular pair formation on oxide surfaces. The results have implications for the understanding of NOx storage and reduction catalysts.

We review our understanding of the NO2 interaction with BaO. The presented picture has evolved from calculations using the density functional theory, reactor measurements and vibrational spectroscopy of surface species during NO2 storage on a BaO/Al2O3 catalyst and BaO powder. The DFT calculations predict nitrite formation on the BaO(1 0 0) surface upon NO2 adsorption. A particular stable adsorption configuration is a nitrate-nitrite pair geometry, with the nitrate bonded to a barium cation and the nitrite involving a surface oxygen anion. This configuration is supported by vibrational spectroscopy. Only nitrite formation is observed for NO exposure to BaO powder at low temperatures, whereas signatures at similar to250degreesC reveal the formation of surface nitrates. Continuous NO2 adsorption at this temperature results in surface Ba(NO3)(2), as well as bulk Ba(NO3)(2), formation.

Structural and energetic surface properties of the alkaline earth metal oxides MgO, CaO, SrO, and BaO are investigated within the density functional theory. In particular, structural distortions (relaxation and rumpling) and surface energies are studied for the (100) and (110) surfaces. A different sign of the rumpling is calculated for MgO compared to the other oxides. This difference is carefully examined and discussed in terms of differences in electronic properties, as well as nearest and next-nearest neighbor interactions. The stability of the results is investigated with respect to model approximations, such as the finite number of atomic layers used to describe semi-infinite surfaces, and choice of exchange-correlation energy functional.

The density functional theory has been applied to study the interaction of platinum (atoms, dimers, and layers with equal number of platinum atoms as oxygen surface atoms) with the (100) surfaces of MgO and BaO. The systems are characterized with respect to preferred geometrical configurations and electronic structure. Platinum is found to interact considerably stronger (∼1.5 eV) with BaO than with MgO. The mapping of the potential energy surface suggests that monomer diffusion at room temperature is unlikely on both studied supports. Dimer diffusion, however, should be possible to observe at room temperature on MgO(100). For layer adsorption, a marked difference is predicted between MgO and BaO. Whereas the Pt/MgO interaction is weak, the Pt/BaO interaction is close to that of atoms and dimers. Substrate induced modifications of chemical properties are also investigated by studying variations in CO adsorption energy. CO is predicted to bind more strongly to supported atoms and dimers than to the correspo...