Search

A comprehensive understanding of the electrode-electrolyte interface in energy conversion systems remains challenging due to the complex and multifaceted nature of interfacial processes. This complexity hinders the development of more efficient electrocatalysts. In this work, we propose a hybrid approach to the theoretical description of the OER process on nickel-iron-based oxyhydroxides ($\gamma$-Ni$_{1-x}$Fe$_x$OOH) electrodes in alkaline media as a model system. Multiple reaction pathways represented by the single- and dual-site mechanisms were investigated by taking into account the realistic structure of the catalyst, the doping, and the solvation effects using a simple and computationally feasible strategy. Accounting for the variable solvation effects considerably affects the predicted overpotential in a roughly linear relationship between overpotential and dielectric constant. By incorporating quantum chemical simulations with kinetic modeling, we demonstrate that tuning the local solvation environment can significantly enhance the OER activity, opening new routine ways for elucidation of the emerging issues of OER processes on transition metal oxide surfaces and design of cost-effective, efficient electrocatalytic systems.

The catalytic activities of high-spin small Fe(III) oxides have been investigated for efficient hydrogen production through ammonia decomposition, using the Artificial Force Induced Reaction (AFIR) method within the framework of density functional theory (DFT) with the B3LYP hybrid exchange-correlation functional. Our results reveal that the adsorption free energy of NH$_3$ on (Fe$_2$O$_3$)$_n$ ($n=1-4$) decreases with increasing cluster size up to $n=3$, followed by a slight increase at $n=4$. The strongest NH$_3$ adsorption energy, 33.68 kcal/mol, was found for Fe$_2$O$_3$, where NH$_3$ interacts with a two-coordinated Fe site, forming an Fe-N bond with a length of 2.11 \AA. A comparative analysis of NH$_3$ decomposition and H$_2$ formation on various Fe(III) oxide sizes identifies the rate-determining steps for each reaction. We found that the rate-determining step for the full NH$_3$ decomposition on (Fe$_2$O$_3$)$_n$ ($n=1-4$) is size-dependent, with the NH$^{*}$ $\rightleftharpoons$ N$^{*}$ + 3H$^{*}$ reaction acting as the limiting step for $n=1-3$. Additionally, our findings indicate that H$_2$ formation is favored following the partial decomposition of NH$_3$ on Fe(III) oxides.

Uncertainty control and scalability to large datasets are the two main issues for the deployment of Gaussian process (GP) models within the autonomous machine learning-based prediction pipelines in material science and chemistry. One way to address both of these issues is by introducing the latent inducing point variables and choosing the right approximation for the marginal log-likelihood objective. Here, we empirically show that variational learning of the inducing points in a molecular descriptor space improves the prediction of energies and atomic forces on two molecular dynamics datasets. First, we show that variational GPs can learn to represent the configurations of the molecules of different types that were not present within the initialization set of configurations. We provide a comparison of alternative log-likelihood training objectives and variational distributions. Among several evaluated approximate marginal log-likelihood objectives, we show that predictive log-likelihood provides excellent uncertainty estimates at the slight expense of predictive quality. Furthermore, we extend our study to a large molecular crystal system, showing that variational GP models perform well for predicting atomic forces by efficiently learning a sparse representation of the dataset., Comment: 17 pages, 6 figures

Due to its exceptional lithium storage capacity silicon is considered as a promising candidate for anode material in lithium-ion batteries (LIBs). In the present work we demonstrate that methods of the soft X-ray emission spectroscopy (SXES) can be used as a powerful tool for the comprehensive analysis of the electronic and structural properties of lithium silicides Li$_{x}$Si forming in LIB's anode upon Si lithiation. On the basis of density functional theory (DFT) and molecular dynamics (MD) simulations it is shown that coordination of Si atoms in Li$_{x}$Si decreases with increase in Li concentration both for the crystalline and amorphous phases. In amorphous a-Li$_{x}$Si alloys Si tends to cluster forming Si-Si covalent bonds even at the high lithium concentration. It is demonstrated that the Si-L$_{2,3}$ emission bands of the crystalline and amorphous Li$_{x}$Si alloys show different spectral dependencies reflecting the process of disintegration of Si-Si network into Si clusters and chains of the different sizes upon Si lithiation. The Si-L$_{2,3}$ emission band of Li$_{x}$Si alloys become narrower and shifts towards higher energies with an increase in Li concentration. The shape of the emission band depends on the relative contribution of the X-ray radiation from the Si atoms having different coordination. This feature of the Si-L$_{2,3}$ spectra of Li$_{x}$Si alloys can be used for the detailed analysis of the Si lithiation process and LIB's anode structure identification.

The four-electron oxygen reduction reaction on Pt catalyst in alkaline solution undergoes proton transfer via tunneling mechanism. The hydrogen/deuterium kinetic isotopic rate constant ratio (kH/kD ) = 32 in a low overpotential region, indicating the importance of the quantum-proton-tunneling at the rate-determining step (RDS). However, kH/kD goes down to 3 in a high overpotential region, suggesting the classical proton-transfer (PT) scheme. Therefore, there is a quantum-to-classical transition of PT process as a function of potential, which is confirmed by theoretical study., Comment: 12 pages, 3 figures, Supplemental Material is available at the web page of PRL

The specific role played by small gold nanoparticles supported on the rutile TiO2(110) surface in the processes of adsorption and dissociation of H2 is discussed. It is demonstrated that the molecular and dissociative adsorption of H2 on Au_n clusters containing n = 1, 2, 8 and 20 atoms depends on cluster size, geometry structure, cluster flexibility and the interaction with the support material. Rutile TiO2(110) support energetically promotes H2 dissociation on gold clusters. It is demonstrated that the active sites towards H2 dissociation are located at corners and edges on the surface of the gold nanoparticle in the vicinity of the support. The low coordinated oxygen atoms on the TiO2(110) surface play a crucial role for H2 dissociation. Therefore the catalytic activity of a gold nanoparticle supported on the rutile TiO2(110) surface is proportional to the length of the perimeter interface between the nanoparticle and the support., Comment: 26 pages, 8 figures, in press Faraday Discussions (2011)

The adsorption of ethylene molecule on neutral, anionic and cationic gold clusters consisting of up to 10 atoms has been investigated using density-functional theory. It is demonstrated that C2H4 can be adsorbed on small gold clusters in two different configurations, corresponding to the pi- and di-sigma-bonded species. Adsorption in the pi-bonded mode dominates over the di-sigma mode over all considered cluster sizes n, with the exception of the neutral C2H4-Au5 system. A striking difference is found in the size-dependence of the adsorption energy of C2H4 bonded to the neutral gold clusters in the pi and di-sigma configurations. The important role of the electronic shell effects in the di-sigma mode of ethylene adsorption on neutral gold clusters is demonstrated. It is shown that the interaction of C2H4 with small gold clusters strongly depends on their charge. The typical shift in the vibrational frequencies of C2H4 adsorbed in the pi- and the di-sigma configurations gives a guidance to experimentally distinguish between the two modes of adsorption., Comment: 30 pages, 10 figures

We demonstrate that the addition of a single carbon impurity leads to significant changes in the thermodynamic properties of Ni clusters consisting of more than a hundred atoms. The magnitude of the change induced is dependent upon the parameters of the Ni-C interaction. Hence, thermodynamic properties of Ni clusters can be effectively tuned by the addition of an impurity of a particular type. We also show that the presence of a carbon impurity considerably changes the mobility and diffusion of atoms in the Ni cluster at temperatures close to its melting point. The calculated diffusion coefficients of the carbon impurity in the Ni cluster can be used for a reliable estimate of the growth rate of carbon nanotubes., Comment: 27 pages, 13 figures

The optimized structure and electronic properties of neutral and singly charged magnesium clusters have been investigated using ab initio theoretical methods based on density-functional theory and systematic post-Hartree-Fock many-body perturbation theory accounting for all electrons in the system. We have systematically calculated the optimized geometries of neutral and singly charged magnesium clusters consisting of up to 21 atoms, electronic shell closures, binding energies per atom, ionization potentials and the gap between the highest occupied and the lowest unoccupied molecular orbitals. We have investigated the transition to the hcp structure and metallic evolution of the magnesium clusters, as well as the stability of linear chains and rings of magnesium atoms. The results obtained are compared with the available experimental data and the results of other theoretical works., Comment: 30 pages, 10 figures, 3 tables

This work is devoted to the elucidation the applicability of jellium model to the description of alkali cluster properties on the basis of comparison the jellium model results with those derived from experiment and within ab initio theoretical framework. On the basis of the Hartree-Fock and local-density approximation deformed jellium model we have calculated the binding energies per atom, ionization potentials, deformation parameters and the optimized values of the Wigner-Seitz radii for neutral and singly charged sodium clusters with the number of atoms N <= 20. These characteristics are compared with the results derived from the ab initio all-electron simulations of cluster electronic and ionic structure based on the density functional theory as well as on the post Hartree-Fock perturbation theory on many-electron correlation interaction. The comparison performed demonstrates the great role of cluster shape deformations in the formation cluster properties and the quite reasonable level of applicability of the deformed jellium model., Comment: 32 pages, 10 figures

Fission of doubly charged metal clusters is studied using the open-shell two-center deformed jellium Hartree-Fock model and Local Density Approximation. Results of calculations of the electronic structure and fission barriers for the symmetric and asymmetric channels associated with the following processes Na_{10}^{2+} --> Na_{7}^{+} + Na_{3}^{+}, Na_{18}^{2+} --> Na_{15}^{+} + Na_{3}^{+} and Na_{18}^{2+} --> 2 Na_{9}^{+} are presented. The role of the exact exchange and many-body correlation effects in metal clusters fission is analysed. It is demonstrated that the influence of many-electron correlation effects on the height of the fission barrier is more profound if the barrier arises nearby or beyond the scission point. The importance of cluster deformation effects in the fission process is elucidated with the use of the overlapping-spheroids shape parametrization allowing one an independent variation of deformations in the parent and daughter clusters., Comment: 21 pages, 7 figures. Submitted to Phys.Rev.A