Search

• Surface electrochemical processes on M@N 4 single-atom catalysts are investigated. • Adsorption of H and OH from electrolyte solution depends on the d -band filling of the metal center. • Mn, Fe, Co, Ru, Rh, and Ir-based SACs are prone to oxidation at anodic potentials. • Stability and selectivity towards oxygen reduction reaction are discussed. • Strategies to identify the oxidation of metal centers are outlined. Single-atom catalysts (SACs) are rapidly developing in various application areas, including electrocatalysis of different reactions, usually taking place under harsh pH/electrode potential conditions. Thus, a full atomic-level understanding of the nature of the active sites under realistic electrochemical conditions is needed, having in mind that the state of SACs active centers could be altered by the adsorption of spectating species. In this contribution, Density Functional Theory is employed to conduct thermodynamic analysis of SACs with metal atoms (Mn, Fe, Co, Ni, Cu, Ru, Rh, Pd, Ag, Ir, Pt, or Au) embedded into N 4 moiety in graphene. Various surface electrochemical processes on such SACs are considered, their Pourbaix plots are constructed, and their activity, selectivity, and stability under operating conditions are discussed. It is demonstrated how adsorption of H, O and OH can cause blockage and restructuring of the active sites and alter the electronic structure. Furthermore, when one deals with metals with lower d -band filling, it is shown that metal center oxidation is preferred over the oxidation of carbon lattice. The effect of the state of the metal center on the reactivity of the carbon lattice is discussed in the case of Fe@N 4 -graphene. Finally, a possible strategy for confirming the changes in the architecture of the SACs' active site by analyzing their vibration spectra is suggested. [Display omitted] [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The core extends from the depth of 2,900 km to the center of the Earth and is composed mainly of an iron-rich alloy with nickel, with 10% of the mass comprised of lighter elements like hydrogen, but the exact composition is uncertain. We present a quantum mechanical first-principles study of the dynamical stability of Fell phases and their phonon densities of states at high pressure. Our free-energy calculations reveal a phonon-driven stabilization of dhcp Fell at low pressures, thus resolving the present contradiction between experimental observations and theoretical predictions. Calculations reveal a complex phase diagram for Fell under pressure with a dhcp → hcp → fcc sequence of structural transitions. [ABSTRACT FROM AUTHOR]

The presence of carbon in the Earth makes the search for high-pressure carbon-containing phases essential for our understanding of mineral compositions of the Earth's mantle. In a recent study Isshiki et al. (2004) demonstrated that magnesite transforms into a new phase at lower mantle pressures. However, the structure of the emerging phase remained unknown. Here we show, by means of first principles calculations, that MgCO3 magnesite can transform into a pyroxene structure at 113 GPa, which further transforms into a CaTiO3-type structure at about 200 GPa. [ABSTRACT FROM AUTHOR]

Understanding the catalytic performance of different materials is of crucial importance for achieving further technological advancements. This especially relates to the behaviors of different classes of catalysts under operating conditions. Here, we analyzed the effects of local coordination of metal centers (Mn, Fe, Co) in graphene-embedded single-atom catalysts (SACs). We started with well-known M@N4-graphene catalysts and systematically replaced nitrogen atoms with oxygen or sulfur atoms to obtain M@OxNy-graphene and M@SxNy-graphene SACs (x + y = 4). We show that local coordination strongly affects the electronic structure and reactivity towards hydrogen and oxygen species. However, stability is even more affected. Using the concept of Pourbaix plots, we show that the replacement of nitrogen atoms in metal coordinating centers with O or S destabilized the SACs towards dissolution, while the metal centers were easily covered by O and OH, acting as additional ligands at high anodic potentials and high pH values. Thus, not only should local coordination be considered in terms of the activity of SACs, but it is also necessary to consider its effects on the speciation of SAC active centers under different potentials and pH conditions. [ABSTRACT FROM AUTHOR]

We have investigated the bonding and adhesion of the TiN(001)/fcc-Fe(111) interface with and without defects using density functional theory. Substitutions in Fe, vacancies and impurities in TiN as well as the influence of vertical stress are studied. The performed bonding analysis shows that the main interaction between Fe and TiN originated from the Fe-N bonding, regardless of the specific set-up. For the ideal TiN(001)/fcc-Fe(111) interface, a complete separation of the interface is preferred to any transfer of iron atoms to the TiN surface. The presence of Ti and N vacancies or alloying elements at the interface further promotes this trend. Application of stress to the system lowers sliding barriers along the Tin/Fe interface. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We present an ab initio density functional theory study of the binding behavior of CO and O2 molecules to two- and three-dimensional isomers of Au13 in order to investigate the potential catalytic activity of this cluster towards low-temperature CO oxidation. First, we scanned the potential energy surface of Au13 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 O2 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 Au13. O2, 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 O2 binding. We did not find a synergetic effect due to the co-adsorption of both molecular species. On the three-dimensional isomer, we find O2 dissociation to be exothermic with an dissociation barrier of 1.44 eV. [ABSTRACT FROM AUTHOR]

Abstract We study the behaviour of the lattice parameter and elastic moduli of Ce 1- x Gd x O 2- x /2 for the random (fluorite-like) and C-type ordered oxygen-vacancy configurations [Žguns et al., PCCP 20 (2018) 11805-11818]. For the fluorite phase, elastic moduli decrease linearly with Gd concentration. For the C-type phase, the bulk, shear and Young moduli are found to be systematically larger and the lattice parameter smaller than those for disordered fluorite phase. Essentially the linear behaviour of the bulk modulus and lattice parameter depending on the degree of the C-type order is found. Our findings explain the experimentally observed elastic moduli of Ce 1- x Gd x O 2- x /2. Graphical abstract Unlabelled Image [ABSTRACT FROM AUTHOR]

Highly effective yet affordable non-noble metal catalysts are a key component for advances in hydrogen generation via electrolysis. The synthesis of catalytic heterostructures containing established Ni in combination with surface NiO, Ni(OH)2, and NiOOH domains gives rise to a synergistic effect between the surface components and is highly beneficial for water splitting and the hydrogen evolution reaction (HER). Herein, the intrinsic catalytic activity of pure Ni and the effect of partial electrochemical oxidation of ultra-smooth magnetron sputter-deposited Ni surfaces are analyzed by combining electrochemical measurements with transmission electron microscopy, selected area electron diffraction, X-ray photoelectron spectroscopy, and atomic force microscopy. The experimental investigations are supplemented by Density Functional Theory and Kinetic Monte Carlo simulations. Kinetic parameters for the HER are evaluated while surface roughening is carefully monitored during different Ni film treatment and operation stages. Surface oxidation results in the dominant formation of Ni(OH)2, practically negligible surface roughening, and 3–5 times increased HER exchange current densities. Higher levels of surface roughening are observed during prolonged cycling to deep negative potentials, while surface oxidation slows down the HER activity losses compared to as-deposited films. Thus, surface oxidation increases the intrinsic HER activity of nickel and is also a viable strategy to improve catalyst durability. [ABSTRACT FROM AUTHOR]

The growing energy crisis put an emphasis on the development of novel efficient energy conversion and storage systems. Here we show that surface modification of cobalt by a fast galvanic displacement with rhodium significantly affects the activity towards hydrogen (HER) and oxygen evolution reactions (OER) in alkaline media. After only 20 s of galvanic displacement, the HER overpotential is reduced by 0.16 V and OER overpotential by 0.06 V. This means that the predicted water splitting voltage is reduced from 2.03 V (clean Co anode and cathode) to 1.81 V at 10 mA cm−2 (Rh-exchanged Co electrode). During the galvanic displacement process, the surface roughness of the Co electrode does not suffer significant changes, which suggests an increase in the intrinsic catalytic activity. Density Functional Theory calculations show that the reactivity of the Rh-modified Co(0001) surface is modified compared to that of the clean Co(0001). In the case of HER, experimentally observed activity improvements are directly correlated to the weakening of the hydrogen-surface bond, confirming the beneficial role of Rh incorporation into the Co surface. [ABSTRACT FROM AUTHOR]

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 H ads acceptor and enable fast H ads recombination. General instructions for the choice of the catalyst|support combination can be used to design new advanced HER catalysts. [ABSTRACT FROM AUTHOR]

• Modifying Ni with Rh improves HER/HOR and ORR/OER kinetics in alkaline media. • Oxidation of the Rh-modified electrode has a negative impact on HER activity. • HER is faster, while HOR is slower on Rh-modified electrodes compared to pt. • OER activity of the Rh-modified electrode is high, while ORR kinetics is slower compared to pt. • Slow Rh-oxide reduction and insulating nature of Rh 2 O 3 are responsible for slower HOR and ORR kinetics. Developing novel electrocatalysts for energy conversion applications is of utmost importance for reaching the energy security of modern society. Here we present a comprehensive investigation of rhodium-modified polycrystalline nickel as an electrocatalyst for hydrogen and oxygen electrode reactions in alkaline media. The surface modification of nickel electrodes was achieved by facile galvanic displacement (up to 30 s) from a highly concentrated acidic Rh3+ solution. The results demonstrate a significant enhancement in the electrocatalytic activity for both hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) on the Rh-modified Ni electrodes, positioning galvanic displacement as a viable approach to engineering advanced electrocatalysts for clean energy applications. On the other hand, the hydrogen oxidation (HOR) and oxygen reduction reaction (ORR) activities of the Rh-modified electrodes are lower compared to polycrystalline platinum. It is suggested that semiconducting Rh 2 O 3 has a detrimental role on the HOR and ORR performance, while the activities of HER and OER, dominantly taking place on metallic Rh and conductive RhO 2 , are very high. This research sheds light on the mechanisms underlying the enhanced electrode kinetics on Rh-modified Ni electrodes and provides insights into the development of efficient and cost-effective electrocatalysts for renewable energy technologies. [Display omitted] [ABSTRACT FROM AUTHOR]

As the global energy crisis continues, efficient hydrogen production is one of the hottest topics these days. In this sense, establishing catalytic trends for hydrogen production is essential for choosing proper H2 generation technology and catalytic material. Volcano plots for hydrogen evolution in acidic media are well-known, while a volcano plot in alkaline media was constructed ten years ago using theoretically calculated hydrogen binding energies. Here, for the first time, we show that the volcano-type relationships are largely maintained in a wide range of pH values, from acidic to neutral and alkaline solutions. We do this using theoretically calculated hydrogen binding energies on clean metallic surfaces and experimentally measured hydrogen evolution overpotentials. When metallic surfaces are exposed to high anodic potentials, hydrogen evolution can be boosted or significantly impeded, depending on the type of metal and the electrolyte in which the reaction occurs. Such effects are discussed here and can be used to properly tailor catalytic materials for hydrogen production via different water electrolysis technologies. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Adsorption of 11 non-metals (H, B, C, N, O, F, Si, P, S, Cl and Br) on defect-free MgO(001) surface was investigate using DFT approach. Adsorption energies were found to be between − 0.56 eV (hydrogen adsorption) and − 2.63 eV (carbon adsorption). Charge transfer from substrate to adsorbate was observed to follow the periodicity in the Periodic Table of Elements, as increases from left to right and decreases from top to bottom. All investigated adsorbates prefer oxygen sites on MgO surface. The analysis of adsorbate-MgO(001) electronic structure suggested that the electronic structure of the O adsorption center and adsorbate atom is molecule-like and there is no strong interaction with MgO electronic bands. Based on the obtained dataset for adsorption energies of selected non-metallic adsorbates (X) the reactivity of MgO towards the bond cleavage in the cases of X–X, H–X and HO–X bonds was discussed. Obtained results point to weak reactivity of MgO(001) towards atomic adsorption and low activity for bond cleavage. However, these results can be used as a starting point for the functionalization of MgO, particularly in the cases where bond cleavage activity and surface-mediated stabilization of dissociation products are desired. [ABSTRACT FROM AUTHOR]

Oxides with the cubic fluorite structure, e.g., ceria (CeO2), are known to be good solid electrolytes when they are doped with cations of lower valence than the host cations. The high ionic conductivity of doped ceria makes it an attractive electrolyte for solid oxide fuel cells, whose prospects as an environmentally friendly power source are very promising. In these electrolytes, the current is carried by oxygen ions that are transported by oxygen vacancies, present to compensate for the lower charge of the dopant cations. Ionic conductivity in ceria is closely related to oxygen-vacancy formation and migration properties. A clear physical picture of the connection between the choice of a dopant and the improvement of ionic conductivity in ceria is still lacking. Here we present a quantum-mechanical first-principles study of the influence of different trivalent impurities on these properties. Our results reveal a remarkable correspondence between vacancy properties at the atomic level and the macroscopic ionic conductivity. The key parameters comprise migration barriers for bulk diffusion and vacancy-dopant interactions, represented by association (binding) energies of vacancy-dopant clusters. The interactions can be divided into repulsive elastic and attractive electronic parts. In the optimal electrolyte, these parts should balance. This finding offers a simple and clear way to narrow the search for superior dopants and combinations of dopants. The ideal dopant should have an effective atomic number between 61 (Pm) and 62 (Sm), and we elaborate that combinations of Nd/Sm and Pr/Gd show enhanced ionic conductivity, as compared with that for each element separately. [ABSTRACT FROM AUTHOR]

The interest in single-atom catalysts (SACs) is increasing, as these materials have the ultimate level of catalyst utilization, while novel reactions where SACs are used are constantly being discovered. However, to properly understand SACs and to further improve these materials, it is necessary to consider the nature of active sites under operating conditions. This is particularly important when SACs are used as electrocatalysts due to harsh experimental conditions, including extreme pH values or high anodic and cathodic potential. In this contribution, density functional theory-based thermodynamic modelling is used to address the nature of metal centers in SACs formed by embedding single metal atoms (Ru, Rh, Ir, Ni, Pd, Pt, Cu, Ag, and Au) into graphene monovacancy. Our results suggest that none of these SAC metal centers are clean at any potential or pH in the water thermodynamic stability region. Instead, metal centers are covered with Hads, OHads, or Oads, and in some cases, we observed the restructuring of the metal sites due to oxygen incorporation. Based on these findings, it is suggested that setting up theoretical models for SAC modelling and the interpretation of ex situ characterization results using ultra-high vacuum (UHV) techniques requires special care, as the nature of SAC active sites under operating conditions can significantly diverge from the basic models or the pictures set by the UHV measurements. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract The ability to tune charge storage properties of graphene oxide (GO) is of utmost importance for energy conversion applications. Here we show that the electrochemical reduction of GO is highly sensitive to the cations present in the solution. GO is reduced at more negative potential in alkali metal chloride solutions than in alkaline earth metal chlorides. During the reduction, the capacitance of GO increases from 10 to 70 times. The maximum capacitances of reduced GO are between 65 and 130 F g−1, depending on the electrolyte and the presence of conductive additive. We propose that different interactions of cations with oxygen functional groups of GO during the reduction are responsible for the observed effect. This hypothesis has been confirmed by Density Functional Theory calculations of alkali and alkaline earth metals interactions with epoxy-functionalized graphene sheet. Graphical abstract Image Highlights • Electrochemical reduction of graphene oxide is sensitive to cations. • GO is reduced at more negative potentials in alkali metal chloride solutions. • Capacitance of graphene oxide is maximized when reduced at optimal potential. • Capacitance is tuned between 65 and 130 F g−1. • The effect is due to the interactions of cations with oxygen functional groups. [ABSTRACT FROM AUTHOR]

Abstract 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. [ABSTRACT FROM AUTHOR]

Oxide materials offer a wide range of interesting physical and chemical properties. Even more versatile behavior of oxides is seen at the nanoscale, qualifying these materials for a number of applications. In this study we used DFT calculations to investigate the physical and chemical properties of small hexagonal MgO nanotubes of different length. We analyzed the effect of Li, B, C, N, and F doping on the properties of the nanotubes. We find that all dopants favor the edge positions when incorporated into the nanotubes. Doping results in the net magnetization whose value depends on the type of the impurity. Using the CO molecule as a probe, we studied the adsorption properties of pristine and doped MgO nanotubes. Our results show that the dopant sites are also the centers of significantly altered chemical reactivity. While pristine MgO nanotubes adsorb CO weakly, very strong adsorption at the dopant sites (B-, C-, and N-doped nanotubes) or neighboring edge atoms (F- and Li-doped nanotubes) is observed. Our results suggest that impurity engineering in oxide materials can be a promising strategy for the development of novel materials with possible use as selective adsorbents or catalysts. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The understanding of atomic adsorption on graphene is of high importance for many advanced technologies. Here we present a complete database of the atomic adsorption energies for the elements of the Periodic Table up to the atomic number 86 (excluding lanthanides) on pristine graphene. The energies have been calculated using the projector augmented wave (PAW) method with PBE, long-range dispersion interaction corrected PBE (PBE+D2, PBE+D3) as well as non-local vdW-DF2 approach. The inclusion of dispersion interactions leads to an exothermic adsorption for all the investigated elements. Dispersion interactions are found to be of particular importance for the adsorption of low atomic weight earth alkaline metals, coinage and s-metals (11th and 12th groups), high atomic weight p-elements and noble gases. We discuss the observed adsorption trends along the groups and rows of the Periodic Table as well some computational aspects of modelling atomic adsorption on graphene. [ABSTRACT FROM AUTHOR]

The accuracy of the theoretical description of materials properties in the framework of density functional theory (DFT) inherently depends on the exchange‐correlation (XC) functional used in the calculations. Here we investigate the influence of the choice of a XC functional (PBE, RPBE, PW91, and PBE0) on the kinetics of the adsorption, diffusion and dissociation of water on the rutile TiO2(110) surface using a combined Kinetic Monte Carlo (KMC) – DFT approach, where the KMC simulations are based on the barriers for the aforementioned processes calculated with DFT. We also test how the adsorption energy of intact and dissociated water molecules changes when dispersion interactions are included into the calculations. We consider the beginning of the water layer formation varying coverage up to 0.2 monolayer (ML) at temperatures up to 180 K. We demonstrate that the dynamics of the simulated water–titania system is extremely sensitive to the choice of the XC functional. [ABSTRACT FROM AUTHOR]

Considering the increasing interest in the use of graphene-based materials for energy conversion and storage applications, we have performed a DFT study of Na interaction with doped graphene, both in non-oxidized and oxidized forms. Oxidation seems to play the crucial role when it comes to the interaction of doped graphene materials with sodium. The dopants act as attractors of OH groups, making the material prone to oxidation, and therefore altering its affinity towards Na. In some cases, this can result in hydroxide or water formation − an irreversible change lethal for battery performance. Our results suggest that one should carefully control the oxidation level of doped graphene-based materials if they are to be used as sodium battery electrode materials as the optimal oxidation level depends on the dopant type. [ABSTRACT FROM AUTHOR]

Sm-doped ceria is a prospective electrolyte material for intermediate-temperature solid-oxide fuel cells (IT-SOFC). Equilibrium ab initio molecular dynamics (AIMD) studies of oxygen ion diffusion in this material are currently impractical due to the rareness of diffusive events on the accessible timescale. To overcome this issue we have performed ab initio non-equilibrium molecular dynamics calculations of Sm-doped ceria using the color-diffusion algorithm. Applying an external force field we have been able to increase the frequency of diffusive events over the simulation time, while keeping the physical mechanism of diffusion intact. We have investigated the temperature dependence of the maximum strength of the applied external field that could be used while maintaining the response of the system in a linear regime. This allows one to obtain the diffusivity at zero field. The bulk ionic conductivity has been calculated and found to match the experimental data well. We have also compared the description of the diffusion process by our method to previous findings and show that the migration mechanism and site preference of oxygen vacancies with respect to the Sm dopants is well reproduced. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

• Fast galvanic displacement of Ni with Rh significantly improves hydrogen evolution activity. • After 30 s of exchange, smooth Ni disk shows higher HER activity than polycrystalline Pt. • The same exchange protocols are applied for commercial Ni foam and electrodeposited Ni. • After 10 min of exchange modified Ni foam shows HER overpotential of −0.07 V at −10 mA cm−2. • Surface modification of Ni with Rh has minor impact on the price of the catalyst compared to energy savings. Energy-efficient hydrogen production is one of the key factors for advancing hydrogen-based economy. Alkaline water electrolysis is the main route for the production of high-purity hydrogen, but further improvements of hydrogen evolution reaction (HER) catalysts are still needed. Industrial alkaline electrolysis relies on Ni-based catalysts, and here we describe a drastic improvement of HER activity of Ni in alkaline media using several model catalysts for HER, obtained upon nickel surface modification in the aqueous solution of rhodium salts, where a spontaneous deposition of rhodium takes place, based on the chemical displacement reaction 3Ni + 2Rh3+ = 3Ni2+ + 2Rh. In the case of smooth Ni-poly electrodes, HER activity surpasses the activity of Pt-poly after just 30 s of exchange with Rh. SEM analysis showed that Rh is uniformly distributed, and that surface roughness changes are lower than 10%, which is in agreement with the electrochemical measurements. Furthermore, XPS analysis has shown effective incorporation of Rh in the surface, while DFT calculations suggest that hydrogen binding is significantly weakened on the Rh-modified Ni surfaces. Such tuning of the hydrogen binding energy is seen as the main factor governing HER activity improvements. The same galvanic displacement protocols were employed for nickel foam electrodes and electrodeposited Ni on Ti mesh. In both cases, somewhat longer Rh exchange times are needed to obtain superior activities than for the smooth Ni surface, but within 10 min. HER overpotentials corresponding to −10 mA cm−2 for nickel foam and electrodeposited Ni electrodes, after modification with Rh, amounted to only −0.07 and −0.09 V, respectively. Thus, it is suggested that a fast spontaneous displacement of Ni with Rh could effectively boost HER in alkaline media with minor cost penalties with regards to energy saving in the electrolysis process. [Display omitted] [ABSTRACT FROM AUTHOR]

• Electrochemical reduction of graphene oxide proceeds with activation energy below 30 kJ mol−1. • Reduction of graphene oxide films is faster in KCl than in LiCl-containing electrolyte. • Lateral conductivity through the film is restored rapidly, at deep negative potentials. • Transversal conductivity is restored locally, through the growth of conductive islands which coalesce. • Isolated OH groups are easy to reduce, while cluster oxygen functional groups are more stable. Graphene oxide finds applications in different fields of science, including energy conversion. Electrochemical reduction of graphene oxide (GO) significantly improves its conductivity. However, the kinetics of this process depends on the solvent, supporting electrolyte, pH, and numerous other factors. Most studies report the macroscopic views and ex-situ properties of reduced GO. To expand the knowledge about GO reduction, in this study, we used cyclic voltammetry (CV), simultaneous 2 points and 4 points resistance measurement (s24), conductive atomic force microscopy (AFM), and theoretical calculations. Using CV, we demonstrated that the choice of supporting electrolyte (KCl or LiCl) influences the potential range in which electrochemical GO reduction occurs. The activation energy of this process was estimated to be below 30 kJ mol‒1 in both electrolytes, being significantly lower than that required for thermal reduction of GO. Simultaneous in situ s24 resistance measurements suggest that GO films reach a highly conductive state at deep negative potentials, with an abrupt, irreversible switch from non-conductive to the conductive state. However, conductive AFM presents a more exact picture of this process: the reduction of GO films starts locally while the formed conductive islands grow during the reduction. This mechanism was confirmed by theoretical calculations indicating that the reduction starts on isolated oxygen-functional groups over the GO basal plane, while clustered OH groups are more difficult to reduce. The presented results can help in tailoring reduced GO for a particular electrochemical application by precisely controlling the reduction degree and percentage of the conductive area of the reduced GO films. [Display omitted] [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we show that electrochemical quasi-steady state potential curves of disordered tungsten, titanium, cerium and iridium oxides exhibit fine structure in striking agreement with the density of electronic states, as obtained from ab initio calculations for the crystalline counterparts. The ability to probe the electronic structure by our electrochemical technique is restricted to disordered non-metallic materials. It is highly probable that the probed electronic states are required to be localized, in order for the technique to give a good picture of their density. The electrochemical density of states is often smaller than the computed one due to kinetic effects, i.e. very slow relaxations of the charge carriers. By combining electrochemical and electrical conductivity measurements for the case of tungsten oxide, we have verified the localized character of the probed states and estimated their localization radius. [Copyright &y& Elsevier]

Using x-ray photoelectron spectroscopy of the oxygen 1s core level, the ratio between intact (D2O) and dissociated (OD) water in the hydrated stoichiometric TiO2(110) surface is determined at varying coverage and temperature. In the submonolayer regime, both the D2O∶OD ratio and the core-level binding energy of D2O (ΔBE) decrease with temperature. The observed variations in ΔBE are shown with density functional theory to be governed crucially and solely by the local hydrogen bonding environment, revealing a generally applicable classification and details about adsorption motifs. [ABSTRACT FROM AUTHOR]

• Modelling of oxygen reduction in Li-O 2 and Na-O 2 cells is analysed. • Graphene-based materials are considered as cathode catalysts. • Modelling approach is revisited by combining different levels of theory. • General guidelines for identification of new catalysts are provided. Understanding the reactions in M-O 2 cells (M = Li or Na) is of great importance for further advancement of this promising technology. Computational modelling can be helpful along this way, but an adequate approach is needed to model such complex systems. We propose a new scheme for modelling processes in M-O 2 cells, where reference energies are obtained from high-level theory, CCSD(T), while the interactions of reaction intermediates with catalyst surfaces are extracted from computationally less expensive DFT. The approach is demonstrated for the case of graphene-based surfaces as model catalysts in Li-O 2 and Na-O 2 cells using the minimum viable mechanism. B-doped graphene was identified as the best catalyst amongst considered surfaces, while pristine graphene performs poorly. Moreover, we show that the inclusion of dispersion corrections for DFT has a significant impact on calculated discharge and charge potentials and suggests that long-range dispersion interactions should always be considered when graphene-based materials are modelled as electrocatalysts. Finally, we offer general guidelines for designing new ORR catalysts for M-O 2 cells in terms of the optimization of the interactions of catalyst surface with reaction intermediates. Image, graphical abstract [ABSTRACT FROM AUTHOR]

• Topology of graphene is tuned by combination of strain and choice of the dopant atom. • Reactivity is probed with H and Na - simple adsorbates of great practical importance. • Strain can both enhance and weaken H and Na adsorption on (doped) graphene. • Linear relationship between Na adsorption energy on P-doped graphene and P charge. • Good candidates for hydrogen storage and sodium-ion battery electrodes are discussed. 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 sodium-ion battery electrodes are discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]