We present the results of an extended theoretical study of the structure, phonon, electronic and optical properties of 2D alkaline-earth metal silicides, germanides and stannides (2D Me2X, where Me=Mg, Ca, Sr, Ba and X=Si, Ge, Sn). The performed analysis has shown the occurrence of the pseudo passivation effect and ionic chemical bonding in these 2D Me2X. In addition, the preformed investigation of their phonon spectra has shown the absence of imaginary frequencies indicating the stability of these 2D structures. The band structure calculations performed using the hybrid functional have revealed that all 2D Me2X are semiconductors with the gap varying from 0.12 to 1.01 eV. Among them Mg- and Ca-based 2D materials are direct band-gap semiconductors with the first direct transition having appreciable oscillator strength. We also propose to consider ternary 2D silicides, germanides and stannides with different Me atoms as a feasible way to modify properties of parent 2D Me2X.
By means of ab initio calculations, we have estimated stability of 2D Me2X (Me = Mg, Ca, Sr, Ba and X = Si, Ge, Sn) in the T and Td phases, which are similar to the ones of 2D transition metal chalcogenides, in addition to their phonon spectra. The T phase is found to be more stable for 2D Ca2X, Sr2X and Ba2X, whereas the Td phase is predicted to be the ground state for 2D Mg2X. We have also discussed that imaginary frequencies in the calculated phonon spectra of 2D Me2X, which appeared in the vicinity of the Gamma point, were not necessarily associated with the dynamic instability.
By means of ab-initio techniques we have investigated changes in the structure and electronic properties of alkaline-earth metal silicide (Ca2Si, Mg2Si and MgCaSi) nanotubes caused by the curvature-induced effects. It is revealed that the curvature-induced effects can: 1) stabilize Mg2Si nanotubes in a phase, which is metastable for the parent 2D Mg2Si; 2) lead to an energy gain as a result of 2D to nanotube structural transformation in the case of ternary MgCaSi nanotubes; 3) modify the band dispersion and band gaps for nanotubes with the diameters less than 30 angstrom. In addition, Mg2Si and MgCaSi nanotubes are found to be direct band-gap (0.5-1.2 eV) materials with appreciable oscillator strength of the first direct transitions.
Using density-functional theory in combination with the exact muffin-tin orbital (EMTO) method and coherent potential approximation, we investigate the alloying effect on the tetragonality of Fe-C solid solution forming the basis of steels. In order to assess the accuracy of our approach, first we perform a detailed study of the performance of the EMTO method for the Fe(16)C(1) binary system by comparing the EMTO results to those obtained using the projector augmented wave method. In the second step, we introduce different substitutional alloying elements (Al, Cr, Co, Ni) into the Fe matrix and study their impact on the structural parameters. We demonstrate that a small amount of Al, Co, and Ni enhances the tetragonal lattice ratio of Fe(16)C(1) whereas Cr leaves the ratio almost unchanged. The obtained trends are correlated with the single-crystal elastic parameters calculated for carbon-free alloys.
We show that the formation of the wetting layer and the experimentally observed continuous shift of the H2O-OH balance toward molecular water at increasing coverage on a TiO2(110) surface can be rationalized on a molecular level. The mechanism is based on the initial formation of stable hydroxyl pairs, a repulsive interaction between these pairs, and an attractive interaction with respect to water molecules. The experimental data are obtained by synchrotron radiation photoelectron spectroscopy and interpreted with the aid of density functional theory calculations and Monte Carlo simulations.
We present an ab initio density functional theory study of the binding behavior of CO and O(2) molecules to two-and three-dimensional isomers of Au(13) in order to investigate the potential catalytic activity of this cluster towards low-temperature CO oxidation. First, we scanned the potential energy surface of Au(13) and studied the effect of spin-orbit coupling on the relative stabilities of the 21 isomers we identified. While spin-orbit coupling increases the stability of the three-dimensional more than the two-dimensional isomers, the ground state structure at 0 K remains planar. Second, we systematically studied the binding of CO and O(2) molecules onto the planar and three-dimensional structures lowest in energy. We find that the isomer dimensionality has little effect on the binding of CO to Au(13). O(2), on the other hand, binds significantly to the three-dimensional isomer only. The simultaneous binding of multiple CO molecules decreases the binding energy per molecule. Still, the CO binding remains stronger than the O(2) binding. We did not find a synergetic effect due to the co-adsorption of both molecular species. On the three-dimensional isomer, we find O(2) dissociation to be exothermic with an dissociation barrier of 1.44 eV.
We performed a systematic density functional (DF) study of the adsorption of copper, silver, and gold adatoms on pristine graphene, especially accounting for van der Waals (vdW) interactions by the vdW-DF and PBE + D2 methods. In particular, we analyze the preferred adsorption site (among top, bridge, and hollow positions) together with the corresponding distortion of the graphene sheet and identify diffusion paths. Both vdW schemes show that the coinage metal atoms do bind to the graphene sheet and that in some cases the buckling of the graphene layer can be significant. Only the results for silver are qualitatively at variance with those obtained with the generalized gradient approximation, which gives no binding in this case. However in all three cases, we observe some quantitative differences between the vdW-DF and PBE + D2 methods. For instance the adsorption energies calculated with the PBE + D2 method are systematically higher than the ones obtained with vdW-DF. Moreover, the equilibrium distances computed with PBE + D2 are shorter than those calculated with the vdW-DF method.
Motivated by the experimentally observed high mobility of gold atoms on graphene and their tendency to form nanometer-sized clusters, we present a density functional theory study of the ground state structures of small gold clusters on graphene, their mobility and clustering. Our detailed analysis of the electronic structures identifies the opportunity to form strong gold-gold bonds and the graphene-mediated interaction of the pre-adsorbed fragments as the driving forces behind gold's tendency to aggregate on graphene. While clusters containing up to three gold atoms have one unambiguous ground state structure, both gas phase isomers of a cluster with four gold atoms can be found on graphene. In the gas phase the diamond-shaped Au-4(D) cluster is the ground state structure, whereas the Y-shaped Au-4(Y) becomes the actual ground state when adsorbed on graphene. As we show, both clusters can be produced on graphene by two distinct clustering processes. We also studied in detail the stepwise formation of a gold dimer out of two pre-adsorbed adatoms, as well as the formation of Au-3. All reactions are exothermic and no further activation barriers, apart from the diffusion barriers, were found. The diffusion barriers of all studied clusters range from 4 to 36 meV only, and are substantially exceeded by the adsorption energies of -0.1 to -0.59 eV. This explains the high mobility of Au1-4 on graphene along the C-C bonds.
In order to explain the experimentally found catalytic characteristics of Au1-4 /MgO (100) we have performed a comprehensive density functional study of these systems and their ability to (co)adsorb CO and O2 molecules. Starting from the carefully determined ground-state structures we have analyzed binding mechanisms, the influence of spin-orbit coupling, and charge redistributions in Au1-4 /MgO+CO (O2). Experimentally Au1,2 /MgO were found to be inactive under a mixed atmosphere. We show that O2 strongly binds to Au1 /MgO that prevents coadsorption. Although a catalytic reaction cycle towards CO oxidation, analogous to the gas phase reaction involving Au 2-, is energetically possible for Au2 /MgO, the cluster will get blocked by a strongly bound CO. On the other hand, the catalytic activity of Au3,4 /MgO could be explained by their ability to coadsorb CO and O2, hence indicating the occurrence of a Langmuir- Hinshelwood-type reaction mechanism for these clusters.
The electronic structure and thermodynamic properties of CeO2 and Ce2O3 have been studied from first principles by the all-electron projector-augmented-wave (PAW) method, as implemented in the ab initio total-energy and molecular-dynamics program VASP (Vienna ab initio simulation package). The local density approximation (LDA)+U formalism has been used to account for the strong on-site Coulomb repulsion among the localized Ce 4f electrons. We discuss how the properties of CeO2 and Ce2O3 are affected by the choice of U as well as the choice of exchange-correlation potential, i.e., the local density approximation or the generalized gradient approximation. Further, reduction of CeO2, leading to formation of Ce2O3 and CeO2-x, and its dependence on U and exchange-correlation potential have been studied in detail. Our results show that by choosing an appropriate U it is possible to consistently describe structural, thermodynamic, and electronic properties of CeO2, Ce2O3, and CeO2-x, which enables modeling of redox processes involving ceria-based materials.
The authors have used density functional theory calculations to investigate how the redox thermodynamics and kinetics of CeO2 are influenced by forming solid solutions with TiO2, ZrO2, HfO2, and ThO2. Reduction is facilitated by dissolving TiO2 (largest improvement), HfO2, or ZrO2 (least improvement), while ThO2 makes reduction slightly more difficult. The migration barrier is much lower in the neighborhood of a Ti (largest decrease), Hf, or Zr (least decrease), while the binding energy of solute ions and vacancies increases in the same sequence. They rationalize the properties of ceria solid solutions in terms of defect cluster relaxations.
We have used density functional theory calculations within the LDA+U formulation to investigate how small amounts of dissolved SiO2, GeO2, SnO2, or PbO2 affect the redox thermodynamics of ceria (CeO2). Compared to pure ceria, reduction is facilitated and the reducibility increases in the sequence of CeO2-SnO2, CeO2-GeO2, and CeO2-SiO2, which correlates with the decrease of the ionic radii of the solutes. For low solute concentrations, there is an inverse relation between high reducibility and the solution energy of tetravalent solutes. CeO2-PbO2 is unique in the sense that the initial reduction occurs by Pb(IV)double right arrow Pb(II) instead of the usual Ce(IV)double right arrow Ce(III) reaction. Among the investigated ceria compounds, CeO2-PbO2 has the lowest reduction energy and rather low solution energy. We have studied how the solution and reduction energies depend on the concentration of Si, Ge, Sn, Pb, Ti, Zr, Hf, and Th solute ions. While the solution energy increases monotonously with concentration, the reduction energy first decreases, as compared to pure ceria (except for Th, which exhibits a small increase), and with further increase of solute concentration, it either remains almost constant (Zr, Hf, and Th) or slightly increases (Ti, Si, Ge, and Sn).
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.
We propose a procedure to accurately describe the structural parameters of an incommensurate phase using ab initio methods by approximating it with a set of analogous commensurate supercells. We apply this approach to obtain the structural parameters of the Sc-II phase, which has recently been identified as a complex incommensurate structure similar to Sr-V. The calculated incommensurate ratio gamma, lattice parameters, and Wyckoff positions of Sc-II are in excellent agreement with the available experimental data. Our results show that gamma increases with pressure up to 60 GPa approaching but never reaching the commensurate value 4/3. Hence calculations do not confirm the prediction made based on the reanalyzing of experimental data. When pressure exceeds 70 GPa, gamma shows a sharp decrease that might be considered as the precursor of a new structural phase transition.
Single layer low-dimensional materials are presently of emerging interest, including in the context of magnetism. In the present report, on-surface supramolecular architecturing was further developed and employed to create surface supported two-dimensional binary spin arrays on atomically clean non-magnetic Au(111). By chemical programming of the modules, different checkerboards were produced combining phthalocyanines containing metals of different oxidation and spin states, diamagnetic zinc, and a metal-free 'spacer'. In an in-depth, spectro-microscopy and theoretical account, we correlate the structure and the magnetic properties of these tunable systems and discuss the emergence of 2D Kondo magnetism from the spin-bearing components and via the physico-chemical bonding to the underlying substrate. The contributions of the individual elements, as well as the role of the electronic surface state in the bottom substrate, are discussed, also looking towards further in-depth investigations.
In pursuit of the ideal photocatalyst, cheap and stable semiconductor TiO(2)is considered to be a good choice if one is able to reduce its band gap and decrease the recombination rate of charge carriers. The approach that offers such improvements for energy conversion applications is the modification of TiO(2)with nitrogen and noble metals. However, the origin of these improvements and possibilities for further design of single-atom catalysts are not always straightforward. To shed light on the atomic-scale picture, we modeled the nitrogen-doped (001) anatase TiO(2)surface as a support for palladium and platinum single-atom deposition. The thermodynamics of various synthesis routes for Pd/Pt deposition and nitrogen doping is considered based on density functional theory (DFT)-calculated energies, highlighting the effect of nitrogen doping on metal dimer formation and metal-support interaction. XPS analysis of the valence band of the modified TiO(2)nanocrystals, and the calculated charge transfer and electronic structure of single-atom catalysts supported on the (001) anatase TiO(2)surface provide an insight into modifications occurring in the valence zone of TiO(2)due to nitrogen doping and Pd/Pt deposition at the surface. DFT results also show that substitutional nitrogen doping significantly increases metal-support interaction, while interstitial nitrogen doping promotes only Pt-support interaction.
The face-centered-cubic (fcc) Lennard-Jones (LJ) model can be considered as a representative model of a simple solid. We investigate the mechanism of melting at the limit of superheating in the fcc LJ solid by means of the procedure recently developed by us [Phys. Rev. B 73, 012201 (2006)]. Insight into the mechanism of melting was gained by studying diffusion and defects in the fcc LJ solid by means of molecular dynamics simulations. We found that the limit of superheating achieved by us is likely to be the highest so far. We also found that the size of the cluster which ignites the melting is very small (down to five to six atoms, depending on the size of the supercell) and closely correlates with the linear size of a supercell when the number of atoms varies between 500 and 13 500.
It is well established that at a pressure of several megabars and low temperature Fe is stable in the hexagonal-close-packed (hcp) phase. However, there are indications that on heating a high-pressure hcp phase of Fe transforms to a less dense (open structure) phase. Two phases have been suggested as candidates for these high-temperature stable phases: namely, body-centered-cubic and body-centered-tetragonal (bct) phases. We performed first-principles molecular dynamics and phonon analysis of the bct Fe phase and demonstrated its dynamical instability. This allows us to dismiss the existence of the bct Fe phase under the high-pressure high-temperature conditions of the Earth's inner core.
Earth's solid-iron inner core has a low rigidity that manifests itself in the anomalously low velocities of shear waves as compared to shear wave velocities measured in iron alloys. Normally, when estimating the elastic properties of a polycrystal, one calculates an average over different orientations of a single crystal. This approach does not take into account the grain boundaries and defects that are likely to be abundant at high temperatures relevant for the inner core conditions. By using molecular dynamics simulations, we show that, if defects are considered, the calculated shear modulus and shear wave velocity decrease dramatically as compared to those estimates obtained from the averaged single-crystal values. Thus, the low shear wave velocity in the inner core is explained.
Earth's solid- iron inner core is elastically anisotropic. Sound waves propagate faster along Earth's spin axis than in the equatorial plane. This anisotropy has previously been explained by a preferred orientation of the iron alloy hexagonal crystals. However, hexagonal iron becomes increasingly isotropic on increasing temperature at pressures of the inner core and is therefore unlikely to cause the anisotropy. An alternative explanation, supported by diamond anvil cell experiments, is that iron adopts a body- centered cubic form in the inner core. We show, by molecular dynamics simulations, that the body- centered cubic iron phase is extremely anisotropic to sound waves despite its high symmetry. Direct simulations of seismic wave propagation reveal an anisotropy of 12%, a value adequate to explain the anisotropy of the inner core.
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.
Hole localization accompanying the formation of a cation vacancy in bulk MgO and CaO and at the (100) MgO/CaO interface is described using the Heyd-Scuseria-Ernzerhof hybrid functionals and DFT + U method. The ground state is found to be the O-1-O-1 bipolaronic configuration both in bulk oxides and at their interfaces. The ground-state magnetic configuration is a triplet, which has an energy only about 1-2 meV lower than that of the singlet state. The one-centered O-2-O-0 bipolaron was found to be metastable with its stability being enhanced at the interfaces compared to that in bulk oxides. Possible transition between different biplaronic configurations at the interface are analyzed for chosen configurations.
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.
Multipolaron solutions were studied in the framework of the Holstein one-dimensional molecular crystal model. The study was performed in the continuous limit where the crystal model maps into the nonlinear Schrödinger equation for which a new periodic dnoidal solution was found for the multipolaron system. In addition, the stability of the multi-polaron solutions was examined, and it was found that cnoidal and dnoidal solutions stabilize in different ranges of the parameter space. Moreover, the model was studied under the influence of nonlocal effects and the polaronic dynamics was described in terms of internal solitonic modes.
The results of scanning tunneling spectroscopy of the electronic states of Au nanoclusters on the graphite surface are presented. The tunneling current is found to be different at different points of a rough-surface nano-cluster. The measured differential current-voltage curve of the clusters is nonmonotonic near the Fermi energy, and the tunneling conductance decreases by almost a factor of two as the cluster volume changes from I to 0.1 mm(3). This decrease can be associated with the change in the density of the electronic states near the Fermi energy. The observed features are qualitatively described within the framework of the mechanism of electron localization in disordered systems.
The transport properties of one-dimensional (1D) systems have been studied theoretically. Contradictory experimental results on molecular transport in quasi-1D systems, such as zeolite structures, when both diffusion transport acceleration and the existence of the diffusion mode with lower particle mobility (single-file diffusion (< x(2)> similar to t(1/2))) have been reported, are consolidated in a consistent model. Transition from the single-file diffusion mode to an Einstein-like diffusion < x(2)> similar to t with diffusion coefficient increasing with the density has been predicted to occur at large observation times.
The properties of monoatomic chains have been studied theoretically by means of statistical mechanics methods. The applied approach can be used to evaluate the interatomic distances and lifetimes of one-dimensional (1D) and quasi-1D systems. In particular, we show that the 1D clusters of gold atoms can exist in two states with different lattice parameters (similar to 3.6 and similar to 2.8 angstrom) that can explain the whole variety of experimental observations on monoatomic gold chains without assuming any wire contamination.
Here, we evaluate three different noble metal co-catalysts (Pd, Pt, and Au) that are present as single atoms (SAs) on the classic benchmark photocatalyst, TiO2. To trap the single atoms on the surface, we introduced controlled surface vacancies (Ti3+-Ov) on anatase TiO2 nanosheets by a thermal reduction treatment. After anchoring identical loadings of single atoms of Pd, Pt, and Au, we measure the photocatalytic H-2 generation rate and compare it to the classic nanoparticle co-catalysts on the nanosheets. While nanoparticles yield the well-established the hydrogen evolution reaction activity sequence (Pt > Pd > Au), for the single atom form, Pd radically outperforms Pt and Au. Based on density functional theory (DFT), we ascribe this unusual photocatalytic co-catalyst sequence to the nature of the charge localization on the noble metal SAs embedded in the TiO2 surface.
Today the search for new energy resources is a crucial topic for materials science. The development of new effective catalysts for the oxygen reduction reaction can significantly improve the performance of fuel cells as well as electrocatalytic hydrogen production. This study presents the scalable synthesis of nitrogen-doped graphene oxide for the oxygen reduction reaction. The combination of an ab initio theoretical investigation of the oxygen reduction reaction (ORR) mechanism and detailed electrochemical characterization allowed the identification of electrocatalytically active nitrogen functionalities. The dominant effect on electrocatalytic activity is the presence of graphitic and pyridinic nitrogen and also N-oxide functionalities. The overpotential of ORR for nitrogen-doped graphene oxide prepared by microwave-assisted synthesis outperformed the metal-doped graphene materials. (C) 2018 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.
We combined the hybrid density functional theory (DFT) calculations and X-ray absorption spectroscopy (XAS) experiments in the study of the local atomic structure around Ir ions in ZnO thin films with different iridium content. This was then used in the first principles analysis of the thermoelectric properties of material. The emphasis has been put on the conditions for a positive Seebeck coefficient and p-type electrical conductivity as the functions of the Fermi level. We studied both computationally and experimentally several possible IrOx polyhedra (complexes) with a different number of surrounding oxygens and Ir oxidation states, including those with the formation of peroxide ions (O-2(2-)). In particular, octahedral coordination of iridium ions was identified by reverse Monte Carlo (RMC) simulations of the Ir L-3-edge EXAFS spectra of ZnO:Ir thin films as the predominant complex, which is supported by the calculated lowest interstitial oxygen incorporation energies. All the calculated IrOx (x = 4, 5, 6) complexes, regardless of Ir the oxidation state, demonstrate potential for p-type conduction if the Fermi level (mu(F)) falls in the range of 0-0.8 eV from the valence band maximum (VBM) and the Ir concentration is high enough (12.5% in the present DFT calculations). Even though the corresponding calculated Seebeck coefficient (S) around 80-89 mu V K-1 slightly exceeds the experimental values, we emphasise the presence of an important plateau in the dependence of S on mu(F) in this range for two complexes with the formation of peroxide ions (O-2(2-)). We predicted also that peroxide ions O-2(2-) are characterized by the calculated phonon frequencies of 810-942 cm(-1) in agreement with our previous Raman experimental results. In this light, we discuss the high sensitivity of calculated S(mu(F)) dependences to the atomic and electronic structure.
The melting curve of hydrogen was computed for pressures up to 200 GPa, using molecular dynamics. The inter- and intramolecular interactions were described by the reactive force field (ReaxFF) model. The model describes the pressure-volume equation of state solid hydrogen in good agreement with experiment up to pressures over 150 GPa, however the corresponding equation of state for liquid deviates considerably from density functional theory calculations. Due to this, the computed melting curve, although shares most of the known features, yields considerably lower melting temperatures compared to extrapolations of the available diamond anvil cell data. This failure of the ReaxFF model, which can reproduce many physical and chemical properties (including chemical reactions in hydrocarbons) of solid hydrogen, hints at an important change in the mechanism of interaction of hydrogen molecules in the liquid state.
We have analyzed by density functional theory calculations the structural and magnetic properties of Fe-Co alloys doped by carbon. In analogy with the formation of martensite in steels we predict that such a structure also forms for Fe-Co alloys in a wide range of concentrations. These alloys are predicted to have a stable tetragonal distortion, which in turn leads to an enhanced magnetocrystalline anisotropy energy of up to 0.75 MJ/m(3) and a saturated magnetization field of 1.9 T.
Due to its unique physical and chemical properties, graphene is being considered as a promising material for energy conversion and storage applications. Introduction of functional groups and dopants on/in graphene is a useful strategy for tuning its properties. In order to fully exploit its potential, atomic-level understanding of its interaction with species of importance for such applications is required. We present a DFT study of the interaction of sodium atoms with epoxy-graphene and analyze how this interaction is affected upon doping with boron and nitrogen. We demonstrate how the dopants, combined with oxygen-containing groups alter the reactivity of graphene towards Na. Dopants act as attractors of epoxy groups, enhancing the sodium adsorption on doped oxygen-functionalized graphene when compared to the case of non-doped epoxy-graphene. Furthermore, by considering thermodynamics of the Na interaction with doped epoxy-graphene it has been concluded that such materials are good candidates for Na storage applications. Therefore, we suggest that controlled oxidation of doped carbon materials could lead to the development of advanced anode materials for rechargeable Na-ion batteries.
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.
For carbon-based materials, in contrast to metal surfaces, a general relationship between strain and reactivity is not yet established, even though there are literature reports on strained graphene. Knowledge of such relationships would be extremely beneficial for understanding the reactivity of graphene-based surfaces and finding optimisation strategies which would make these materials more suitable for targeted applications. Here we investigate the effects of compressive and tensile strain (up to +/- 5%) on the structure, electronic properties and reactivity of pure, N-doped and P-doped graphene, using DFT calculations. We demonstrate the possibility of tuning the topology of the graphene surface by strain, as well as by the choice of the dopant atom. The reactivity of (doped) strained graphene is probed using H and Na as simple adsorbates of great practical importance. Strain can both enhance and weaken H and Na adsorption on (doped) graphene. In case of Na adsorption, a linear relationship is observed between the Na adsorption energy on P-doped graphene and the phosphorus charge. A linear relationship between the Na adsorption energy on flat graphene surfaces and strain is found. Based on the adsorption energies and electrical conductivity, potentially good candidates for hydrogen storage and sodiumion battery electrodes are discussed.
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.
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.
In line with a growing interest in the use of graphene-based materials for energy storage applications and active research in the field of rechargeable metal-ion batteries we have performed a DFT based computational study of alkali metal atoms (Li, Na and K) interaction with an oxidized graphene. The presence of oxygen surface groups (epoxy and hydroxyl) alters the chemisorption properties of graphene. In particular, we observe that the epoxy groups are redox active and enhance the alkali metal adsorption energies by a factor of 2 or more. When an alkali metal atom interacts with hydroxyl-graphene the formation of metal-hydroxide is observed. In addition to a potential boost of metal ion storage capability, oxygen functional groups also prevent the precipitation of the metal phase. By simulating lithiation/de-lithiation process on epoxy-graphenes, it was concluded that the oxidized graphene can undergo structural changes during battery operation. Our results suggest that the content and the type of oxygen surface groups should be carefully tailored to maximize the performance of metal-ion batteries. This is mainly related to the control of the oxidation level in order to provide enough active centers for metal ion storage while preserving sufficient electrical conductivity.
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.
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 N-2-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 O-2 in 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 O-2 molecule and contribute to the decomposition of peroxide.
The conductance of monoatomic gold wires containing 3-7 gold atoms has been obtained from ab initio calculations. The transmission is found to vary significantly depending on the wire stretching and the number of incorporated atoms. Such oscillations are determined by the electronic structure of the one-dimensional (1D) part of the wire between the contacts. Our results indicate that the conductivity of 1D wires can be suppressed without breaking the contact.
Hydrogen evolution reaction (HER) is one of the most important reactions in electrochemistry. This is not only because it is the simplest way to produce high purity hydrogen and the fact that it is the side reaction in many other technologies. HER actually shaped current electrochemistry because it was in focus of active research for so many years (and it still is). The number of catalysts investigated for HER is immense, and it is not possible to overview them all. In fact, it seems that the complexity of the field overcomes the complexity of HER. The aim of this review is to point out some of the latest developments in HER catalysis, current directions and some of the missing links between a single crystal, nanosized supported catalysts and recently emerging, single-atom catalysts for HER.
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.