Search

Quantum devices depend on addressable elements, which can be modified separately and in their mutual interaction. Self-assembly at surfaces, for example, formation of a porous (metal-) organic network, provides an ideal way to manufacture arrays of identical quantum boxes, arising in this case from the confinement of the electronic (Shockley) surface state within the pores. We show that the electronic quantum box state as well as the interbox coupling can be modified locally to a varying extent by a selective choice of adsorbates, here C60, interacting with the barrier. In view of the wealth of differently acting adsorbates, this approach allows for engineering quantum states in on-surface network architectures.

The increasing demand and high prices of advanced catalysts motivate a constant search for novel active materials with reduced contents of noble metals. The development of thin films and core-shell catalysts seems to be a promising strategy along this path. Using density functional theory we have analyzed a number of surface properties of supported bimetallic thin films with the composition A(3)B (where A = Pt and Pd, and B = Cu, Ag and Au). We focus on the surface segregation, dissolution stability and surface electronic structure. We also address the chemisorption properties of Pd3Au thin films supported by different substrates, by probing the surface reactivity with CO. We find a strong influence of the support in the case of mono- and bilayers, while the surface strain seems to be the predominant factor in determining the surface properties of supported trilayers and thicker films. In particular, we show that the studied properties of the supported trilayers can be predicted from the lattice mismatch between the overlayer and the support. Namely, if the strain dependence of the corresponding quantities for pure strained surfaces is known, the properties of strained supported trilayers can be reliably estimated. The obtained results can be used in the design of novel catalysts and predictions of the surface properties of supported ultrathin catalyst layers.

Vacancies in graphene present sites of altered chemical reactivity and open possibilities to tune graphene properties by defect engineering. The understanding of chemical reactivity of such defects is essential for successful implementation of carbon materials in advanced technologies. We report the results of a systematic DFT study of atomic adsorption on graphene with a single vacancy for the elements of rows 1-6 of the periodic table of elements (PTE), excluding lanthanides. The calculations have been performed using the PBE, long-range dispersion interaction-corrected PBE (PBE+D2 and PBE+D3) and non-local vdW-DF2 functionals. We find that most elements strongly bind to the vacancy, except for the elements of groups 11 and 12, and noble gases, for which the contribution of dispersion interaction to bonding is most significant. The strength of the interaction with the vacancy correlates with the cohesive energy of the elements in their stable phases: the higher the cohesive energy is, the stronger bonding to the vacancy can be expected. As most atoms can be trapped at the SV site we have calculated the potentials of dissolution and found that in most cases the metals adsorbed at the vacancy are more "noble" than they are in their corresponding stable phases.

The incorporation of trace amounts (<0.2%) of Co and Ni noticeably enhanced the catalytic activity of nitrogen-free ordered mesoporous carbon (OMC) towards oxygen reduction reaction (ORR). (Co,Ni)-doped OMCs were characterized by N2-adsorption measurements, X-ray powder diffraction, field emission scanning electron microscopy and Raman spectroscopy methods, and their ORR activity was estimated by voltammetry on rotating disk electrode in acidic and alkaline media. (Co,Ni)-doped OMCs show modest activities in acidic media, while the catalytic activity in alkaline media is rather high. The measured activities are compared to the Pt-based and Pt-free ORR catalysts reported in the literature. The number of electrons consumed per O2in metal-doped OMCs was found to vary between 2 and 4, which is advantageous in comparison to metal-free OMC. Also, the mass activities of metal-doped OMCs were found to be up to 2.5 times higher compared to that of metal-free OMC. We suggest that the ORR activity is governed by a balance between (i) textural properties, which determine the electrochemically accessible surface of the catalyst and which are influenced by the addition of a metal precursor, and (ii) novel active sites formed upon the introduction of metals into the carbon structure. In particular, our Density Functional Theory calculations suggest that Co and Ni atoms embedded into the single vacancies of graphene can activate the O2molecule and contribute to the decomposition of peroxide.

We examine the effects of the dopant type and the dopant distribution on the ion diffusion in ceria doped with rare-earth elements (Pr, Nd, Pm, Sm, Eu, and Gd). Diffusion is simulated by means of a Kinetic Monte Carlo method using input transition rates derived from diffusion barriers calculated in the framework of density functional theory (DFT). Based on diffusion simulations, we discuss the characteristics of the dopants in terms of the diffusion barriers, and study oxygen ion trajectories for different dopants and distributions. Our simulations show a trend of increasing ion diffusivity with increasing atomic number for all distributions.

Ordering of dopants and oxygen vacancies is studied for Gd-doped ceria (x(Gd) <= 0.25) by means of a combined density functional theory (DFT) and cluster expansion approach, where the cluster interactions derived from DFT calculations are further used in Monte Carlo simulations. The methodology is meticulously tested and the stability of the obtained solutions with respect to the volume change, applied exchange-correlation approximation and other modelling parameters is carefully analysed. We study Gd and vacancy ordering in the case of thermodynamic equilibrium and vacancy ordering for quenched Gd configurations. We find that at the thermodynamic equilibrium there exists a transition temperature (T-C) below which phase separation into C-type Gd2O3 and pure CeO2 occurs. The phase separation is observed in the whole studied concentration range and the transition temperature increases with concentration from ca. 600 (x(Gd) = 0.03) to 1000 K (x(Gd) = 0.25). Above T-C the distribution of Gd is random, oxygen vacancies tend to cluster in the coordination shells along < 1, 1/2, 0 > and < 1, 1, 1 >, and the nearest neighbour position is preferred for Gd-vacancy. In the quenched Gd case, where Gd atoms are immobilised below 1500 K, the vacancy ordering is significantly frustrated. In fact, we observe an oxygen freezing transition below temperature T-F approximate to T-C - 350 K, which is close to temperatures at which a change in the conductivity slope is observed experimentally.

Recently, nanostructured TiO2 ("black TiO2") has been discovered to absorb visible light, which makes, it an efficient material for water splitting. Hydrogenization has been proposed to be at the origin of this beneficial electronic structure of black TiO2. Here, we investigate, using ab initio methods, alternative mechanisms related to structure modifications in nanoclusters that could be responsible for absorption in the visible range. To that end, we apply a combination of computational structure prediction using simulated annealing and minima-hopping methods based on density-functional theory to predict low-energy configurations and time-dependent density-functional theory (TDDFT) using a hybrid functional with optimized Hartree Fock content to obtain optical absorption edges., QC 20170614

LaxCa1-xMnO3 (LCMO) has been studied in the framework of density functional theory (DFT) using Hubbard-U correction. We show that the formation of spin polarons of different configurations is possible in the G-type antiferromagnetic phase. We also show that the spin-polaron (SP) solutions are stabilized due to an interplay of magnetic and lattice effects at lower La concentrations and mostly due to the lattice contribution at larger concentrations. Our results indicate that the development of SPs is unfavorable in the C- and A-type antiferromagnetic phases. The theoretically obtained magnetic state diagram is in good agreement with previously reported experimental results., QC 20170705

Understanding the ways graphene can be functionalized is of great importance for many contemporary technologies. Using density functional theory calculations we investigate how vacancy formation and substitutional doping by B, N, P and S affect the oxidizability and reactivity of the graphene basal plane. We find that the presence of these defects enhances the reactivity of graphene. In particular, these sites act as strong attractors for OH groups, suggesting that the oxidation of graphene could start at these sites or that these sites are the most difficult to reduce. Scaling between the OH and H adsorption energies is found on both reduced and oxidized doped graphene surfaces. Using the O-2 molecule as a probe we show that a proper modelling of doped graphene materials has to take into account the presence of oxygen functional groups., QC 20170522

H-2 production via water electrolysis plays an important role in hydrogen economy. Hence, novel cheap electrocatalysts for the hydrogen evolution reaction ( HER) are constantly needed. Here, we describe a simple method for the preparation of composite catalysts for H-2 evolution, consisting in simultaneous reduction of the graphene oxide film, and electrochemical deposition of Ni on its surface. The obtained composites (Ni@rGO), compared to pure electrodeposited Ni, show an improved electrocatalytic activity towards HER in alkaline media. We found that the activity of the Ni@rGO catalysts depends on the surface composition ( Ni vs. C mole ratio) and on the level of structural disorder of the rGO support. We suggest that HER activity is improved via H-ads spillover from the Ni particles to the rGO support, where quick recombination to molecular hydrogen is favored. A deeper insight into such a mechanism of H-2 production was achieved by kinetic Monte-Carlo simulations. These simulations enabled the reproduction of experimentally observed trends under the assumption that the support can act as a Hads acceptor. We expect that the proposed procedure for the production of novel HER catalysts could be generalized and lead to the development of a new generation of HER catalysts by tailoring the catalyst/support interface., QC 20170619

By means of ab initio calculations we demonstrate sizable anisotropy in surface energy of the gamma-WO3 (001), (010) and (100) surfaces. The (001) surface has the smallest surface energy followed by the (010) and (100) surfaces. Their surface band structures are characterized by dispersion of bands near the gap region and by band-gap values similar to the ones of the bulk. The role of surface atoms in stabilizing the band gap is revealed. Variations in the position of the Fermi level are traced with respect to the different surface reconstructions., QC 20170524

Simulations of diffusion in solids often produce poor statistics of diffusion events. We present an analytical expression for the statistical error in ion conductivity obtained in such simulations. The error expression is not restricted to any computational method in particular, but valid in the context of simulation of Poisson processes in general. This analytical error expression is verified numerically for the case of Gd-doped ceria by running a large number of kinetic Monte Carlo calculations., QC 20160913

A first-principles nonequilibrium molecular dynamics (NEMD) study employing the color-diffusion algorithm has been conducted to obtain the bulk ionic conductivity and the diffusion constant of gadolinium-doped cerium oxide (GDC) in the 850-1150 K temperature range. Being a slow process, ionic diffusion in solids usually requires simulation times that are prohibitively long for ab initio equilibrium molecular dynamics. The use of the color-diffusion algorithm allowed us to substantially speed up the oxygen-ion diffusion. The key parameters of the method, such as field direction and strength as well as color-charge distribution, have been investigated and their optimized values for the considered system have been determined. The calculated ionic conductivity and diffusion constants are in good agreement with available experimental data., QC 20160205

Rational catalyst design presents one of the main paradigms in the contemporary materials science. Although the electronic structure calculations can be used to search for possible candidates, realistic supported catalysts are difficult to address in this way. In this contribution we use conceptual model of the supported hydrogen evolution reaction (HER) catalyst and investigate possible processes using Kinetic Monte Carlo simulations. In specific, we look at the possibility to boost H-2 production by the H spillover to the support and the tailoring of the catalyst deposit. Different scenarios were considered depending on the nature of the HER rate determining step (RDS) on the catalyst surface and the effects of the rates of elementary processes, catalyst dispersion and morphology are analyzed. Metals with low affinity towards hydrogen should be used as catalyst supports, while H spillover can boost H-2 production if Tafel or Heyrovsky reaction is the RDS on the catalyst surface. However, this can be achieved only if the catalyst dispersion is high, while the support has to act as a Hads acceptor and enable fast Hads recombination. General instructions for the choice of the catalystlsupport combination can be used to design new advanced HER catalysts. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved., QC 20160319

In this contribution we inspect the adsorption of H, OH, Cl and Pt on oxidized graphene using DFT calculations. The introduction of epoxy and hydroxyl groups on the graphene basal plane significantly alters its chemisorption properties, which can be attributed to the deformation of the basal plane and the type and distribution of these groups. We show that a general scaling relation exists between the hydrogen binding energies and the binding energies of other investigated adsorbates, which allows for a simple probing of the reactivity of oxidized graphene with only one adsorbate. The electronic states of carbon atoms located within the 2 eV interval below the Fermi level are found to be responsible for the interaction of the basal plane with the chosen adsorbates. The number of electronic states situated in this energy interval is shown to correlate with hydrogen binding energies., QC 20160324

Graphene synthesized by reduction of graphene oxide, depending on the degree of reduction, retains a certain amount of surface OH groups. Considering the surface OH groups/graphene layer system by means of density functional theory calculations, we evidenced the tendency of OH groups to cluster, resulting in enhanced system stability and no band gap opening. In the oxygen concentration range between 1.8 and 8.47 at%, with the addition of each new OH group, integral binding energy decreases, while differential binding energy shows the boost at even numbers of OH groups. Furthermore, we found that the clustering of OH groups over graphene basal plane plays a crucial role in enhancing the interactions with alkali metals. Namely, if alkali metal atoms interact with individual OH groups only, the interaction leads to an irreversible formation of MOH phase. When alkali atoms interact with clusters containing odd number of OH groups, a reversible transfer of an electron charge from the metal atom to the substrate takes place without OH removal. The strength of the interaction in general increases from Li to K. In an experimental investigation of a graphene sample which dominantly contains OH groups, we have shown that the trend in the specific interaction strength reflects to gravimetric capacitances measured in alkali metal chloride solutions. We propose that the charge stored in OH groups which interact with alkali metal cation and the pi electronic system of the graphene basal plane presents the main part of its pseudocapacitance., QC 20160816

Doping of simple oxide materials can give rise to new exciting physical and chemical properties and open new perspectives for a variety of possible applications. Here we use density functional theory calculations to investigate the B-, C- and N-doped MgO(001) surfaces. We have found that the investigated dopants induce magnetization of the system amounting to 3, 2 and 1 mu(B) for B, C and N, respectively. The dopants are found to be in the X2- state and tend to segregate to the surface. These impurity sites also present the centers of altered chemical reactivity. We probe the chemisorption properties of the doped MgO(001) surfaces with the CO molecule and atomic O. The adsorption of CO is much stronger on B- and C-doped MgO(001) compared to pure MgO(001) as the impurity sites serve as potent electron donors. The situation is similar to the case of atomic oxygen, for which we find the adsorption energy of -8.78 eV on B-doped MgO(001). The surface reactivity changes locally around the dopant atom, which is mainly restricted to its first coordination shell. The presented results suggest doped MgO as a versatile multifunctional material with possible use as an adsorbent or a catalyst.

Electron localization and polaron mobility in oxygen-deficient as well as Li-doped monoclinic tungsten trioxide have been studied in the adiabatic limit in the framework of density functional theory. We show that small polarons formed in the presence of oxygen vacancy prefer the bipolaronicW(5+)-W5+ configuration, whereas the W6+-W4+ configuration is found to be metastable. Our calculations suggest that bipolarons are tightly bound by the vacancy and therefore largely immobile. On the contrary, polarons formed as a result of Li intercalation can be mobile; the activation energy for polaron jumping in this case varies between 98 and 124 meV depending on the crystallographic direction. The formation of W5+-W5+ bipolarons in Li-WO3 is possible. When situated along [001] the bipolaronic configuration is 8 meV lower in energy than two separate W5+ polarons.

We have performed a numerical study of the process of electron localization in reduced ceria. Our results show that different localized charge distributions can be attained in a bulk system by varying the lattice parameter. We demonstrate that the effect of electron localization is mainly determined by lattice relaxation and an accurate account for the effects of electronic correlation is necessary to achieve localized charge distribution.

A large number of computational studies have been devoted to the investigation of monometallic clusters supported by MgO. However, in practice, catalysis shows that multicomponent catalytic systems often win in catalytic performance over single component systems. In this study, the geometrical and electronic structure, stability and chemisorption properties of M1M2 metal dimers (M1, M2 = Ru, Rh, Pd, Ir, Pt) supported by defect free MgO(001) have been investigated in the framework of density functional theory. The oxygen sites of MgO(001) are the preferred adsorption sites for all the studied clusters, the majority of them adsorbing parallel to the surface with metal atoms attached to two surface oxygen atoms. The energetics of M1M2 + MgO(001) formation shows that the adsorption complexes are stable and benefit from metal-oxygen and metal-metal interaction. The chemisorption properties of Pd and Pt atoms in PdM2 and PtM2 dimers are studied using CO as a probe molecule. A linear relationship between the CO chemisorption and the d-band center position of the reacting atom in the dimer is observed, extending the d-band center model to the case of highly under-coordinated metal atoms supported by a non-conductive material.

The changes in structural and electronic properties, occurring in one monolayer of MoS2 at different concentrations of oxygen atoms doping and vacancies are investigated by means of ab initio computer simulation. The substitution of sulphur atoms by oxygen ones reduces the band gap for high concentrations only, transforming direct-gap semiconductor into an indirect one, whereas a smaller concentration of oxygen practically does not influence the gap. The presence of sulphur vacancies strongly reduces the band gap, leading to bands overlapping at high concentration and appearance of new bands at the gap region, which are determined by Mo 4d states with the mixture of S 3p states, at low concentrations., QC 20150522

Molybdenum supported thin films of ScN, MgO and NaF with a Cu adatom have been studied in the framework of density functional theory. We have observed a charge transfer from the metal/film interface to the Cu atom and investigated its relation to surface and interface deformations. We find that a weak interaction between the metal and the film is a promising prerequisite for adatom charging. The detailed study of Cu/NaF/Mo and NaF/Mo indicates that the distortion of the NaF film caused by the Cu adsorption has essentially anharmonic character, as it is coupled to a strong charge redistribution in the system., QC 20160218