Dramatic increase in the oxygen reduction reaction for platinum cathodes from tuning the solvent dielectric constant

Formic acid stability in different solvents by DFT calculations

CONTEXT

New technologies have been developed toward the use of green energies. The production of formic acid (FA) from carbon dioxide (CO[Formula: see text]) hydrogenation with H[Formula: see text] is a sustainable process for H[Formula: see text] storage. However, the FA adduct stabilization is thermodynamically dependent on the type of solvent and thermodynamic conditions. The results suggest a wide range of dielectric permittivity values between the dimethyl sulfoxide (DMSO) and water solvents to stabilize the FA in the absence of base. The thermodynamics analysis and the infrared and charge density difference results show that the formation of the FA complex with H[Formula: see text]O is temperature dependent and has a major influence on aqueous solvents compared to the FA adduct with amine, in good agreement with the experiment. In these conditions, the stability thermodynamic of the FA molecule may be favorable at non-organic solvents and dielectric permittivity values closer to water. Therefore, a mixture of aqueous solvents with possible ionic composition could be used to increase the thermodynamic stability of H[Formula: see text] storage in CO[Formula: see text] conversion processes.

METHODS

Using the Quantum ESPRESSO package, density functional theory (DFT) calculations were performed with periodic boundary conditions, and the electronic wave functions were expanded in plane waves. For the exchange-correlation functional, we use the vdW-DF functional with the inclusion of van der Waals (vdW) forces. Electron-ion interactions are treated by the projector augmented wave (PAW) method with pseudopotentials available in the PSlibrary repository. The wave functions and the electronic densities were expanded employing accurate cut-off energies of 6.80[Formula: see text]10[Formula: see text] and 5.44[Formula: see text]10[Formula: see text] eV, respectively. The electronic density was computed from the wave functions calculated at the [Formula: see text]-point in the first Brillouin-zone. Each structural optimization was minimized according to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm, with force and energy convergence criteria of 25 meV[Formula: see text]Å[Formula: see text] and 1.36 meV, respectively. The electrostatic solvation effects were performed by the [Formula: see text] package with the Self-Consistent Continuum Solvation (SCCS) approach.

Adsorption Energy in Oxygen Electrocatalysis

Adsorption energy (AE) of reactive intermediate is currently the most important descriptor for electrochemical reactions (e.g., water electrolysis, hydrogen fuel cell, electrochemical nitrogen fixation, electrochemical carbon dioxide reduction, etc.), which can bridge the gap between catalyst's structure and activity. Tracing the history and evolution of AE can help to understand electrocatalysis and design optimal electrocatalysts. Focusing on oxygen electrocatalysis, this review aims to provide a comprehensive introduction on how AE is selected as the activity descriptor, the intrinsic and empirical relationships related to AE, how AE links the structure and electrocatalytic performance, the approaches to obtain AE, the strategies to improve catalytic activity by modulating AE, the extrinsic influences on AE from the environment, and the methods in circumventing linear scaling relations of AE. An outlook is provided at the end with emphasis on possible future investigation related to the obstacles existing between adsorption energy and electrocatalytic performance.

Infusing theory into deep learning for interpretable reactivity prediction

Despite recent advances of data acquisition and algorithms development, machine learning (ML) faces tremendous challenges to being adopted in practical catalyst design, largely due to its limited generalizability and poor explainability. Herein, we develop a theory-infused neural network (TinNet) approach that integrates deep learning algorithms with the well-established d-band theory of chemisorption for reactivity prediction of transition-metal surfaces. With simple adsorbates (e.g., *OH, *O, and *N) at active site ensembles as representative descriptor species, we demonstrate that the TinNet is on par with purely data-driven ML methods in prediction performance while being inherently interpretable. Incorporation of scientific knowledge of physical interactions into learning from data sheds further light on the nature of chemical bonding and opens up new avenues for ML discovery of novel motifs with desired catalytic properties.

An evaluation of solvent effects and ethanol oxidation

Understanding liquid-metal interfaces in catalysis is important, as the liquid can speed up surface reactions, increase the selectivity of products, and open up new favorable reaction pathways. In this work we modeled using density functional theory various steps in ethanol oxidation/decomposition over Rh(111). We considered implicit (continuum), explicit, and hybrid (implicit combined with explicit) solvation approaches, as well as two solvents, water and ethanol. We focused on modeling adsorption steps, as well as C-C/C-H bond scission and C-O bond formation reactions. Implicit solvation had very little effect on adsorption and reaction free energies. However, using the explicit and hybrid models, some free energies changed significantly. Furthermore, ethanol solvent had a more considerable impact than water solvent. We observed that preferred reaction pathways for C-C scission changed depending on the solvation model and solvent choice (ethanol or water). We also applied the bond-additivity solvation method to calculate heats of adsorption. Heats of adsorption and reaction using the bond-additivity model followed the same trends as the other solvation models, but were ∼1.1 eV more endothermic. Our work highlights how different solvation approaches can influence analysis of the oxidation/decomposition of organic surface species.

Selectivity Map for the Late Stages of CO and CO2 Reduction to C2 Species on Copper Electrodes

The electrochemical CO and CO₂ reduction reactions (CORR and CO₂ RR) using copper catalysts and renewable electricity hold promise as a carbon-neutral route to produce commodity chemicals and fuels. However, the exact mechanisms and structure sensitivity of Cu electrodes toward C₂ products are still under debate. Herein, we investigate ethylene oxide reduction (EOR) as a proxy to the late stages of CORR to ethylene, and the results are compared to those of acetaldehyde reduction to ethanol. Density functional theory (DFT) calculations show that ethylene oxide undergoes ring opening before exclusively reducing to ethylene via *OH formation. Based on generalized coordination numbers (CN), a selectivity map for the late stages of CORR and CO₂ RR shows that sites with moderate coordination (5.9 < CN < 7.5) are efficient for ethylene production, with pristine Cu(100) being more active than defective surfaces such as Cu(311). In contrast, kinks and edges are more active for ethanol production, while (111) terraces are relatively inert.

O2 activation by core-shell Ru13@Pt42 particles in comparison with Pt55 particles: a DFT study

The reaction of O₂ with a Ru₁₃@Pt₄₂ core-shell particle consisting of a Ru₁₃ core and a Pt₄₂ shell was theoretically investigated in comparison with Pt₅₅. The O₂ binding energy with Pt₅₅ is larger than that with Ru₁₃@Pt₄₂, and O-O bond cleavage occurs more easily with a smaller activation barrier (E _a) on Pt₅₅ than on Ru₁₃@Pt₄₂. Protonation to the Pt₄₂ surface followed by one-electron reduction leads to the formation of an H atom on the surface with considerable exothermicity. The H atom reacts with the adsorbed O₂ molecule to afford an OOH species with a larger E _a value on Pt₅₅ than on Ru₁₃@Pt₄₂. An OOH species is also formed by protonation of the adsorbed O₂ molecule, followed by one-electron reduction, with a large exothermicity in both Pt₅₅ and Ru₁₃@Pt₄₂. O-OH bond cleavage occurs with a smaller E _a on Pt₅₅ than on Ru₁₃@Pt₄₂. The lower reactivity of Ru₁₃@Pt₄₂ than that of Pt₅₅ on the O-O and O-OH bond cleavages arises from the presence of lower energy in the d-valence band-top and d-band center in Ru₁₃@Pt₄₂ than in Pt₅₅. The smaller E _a for OOH formation on Ru₁₃@Pt₄₂ than on Pt₅₅ arises from weaker Ru₁₃@Pt₄₂-O₂ and Ru₁₃@Pt₄₂-H bonds than the Pt₅₅-O₂ and Pt₅₅-H bonds, respectively. The low-energy d-valence band-top is responsible for the weak Ru₁₃@Pt₄₂-O and Ru₁₃@Pt₄₂-OH bonds. Thus, the low-energy d-valence band-top and d-band center are important properties of the Ru₁₃@Pt₄₂ particle.

Oxygen Evolution Reaction Kinetic Barriers on Nitrogen-Doped Carbon Nanotubes

We investigate kinetic barriers for the oxygen evolution reaction (OER) on singly and doubly nitrogen-doped single-walled carbon nanotubes (NCNTs) using the climbing image nudged elastic band method with solvent effects represented by a 45-water-molecule droplet. The studied sites were chosen based on a previous study of the same systems utilizing a thermodynamic model which ignored both solvent effects and kinetic barriers. According to that model, the two studied sites, one on a singly nitrogen-doped CNT and the other on a doubly doped CNT, were approximately equally suitable for OER. For the four-step OER process, however, our reaction barrier calculations showed a clear difference in the rate-determining *OOH formation step between the two systems, with barrier heights differing by more than 0.4 eV. Thus, the simple thermodynamic model may alone be insufficient for identifying optimal OER sites. Of the remaining three reaction steps, the two H₂O forming ones were found to be barrierless in all cases. We also performed solvent-free barrier calculations on NCNTs and undoped CNTs. Substantial differences were observed in the energies of the intermediates when the solvent was present. In general, the observed low activation energy barriers for these reactions corroborate both experimental and theoretical findings of the utility of NCNTs for OER catalysis.

How covalence breaks adsorption-energy scaling relations and solvation restores them

It is known that breaking the scaling relations between the adsorption energies of *O, *OH, and *OOH is paramount in catalyzing more efficiently the reduction of O2 in fuel cells and its evolution in electrolyzers. Taking metalloporphyrins as a case study, we evaluate here the adsorption energies of those adsorbates on the metal centers Cr, Mn, Fe, Co, Ni and Cu, using H, F, OH, NH2, CH3, and BH2 as ring ligands. We show that covalence systematically breaks scaling relations under vacuum by strengthening certain M-OOH bonds. However, covalence modifies adsorbate solvation in solution depending on the degree of covalence of the metal-adsorbate bonds. The two effects have similar magnitudes and opposite signs, such that scaling relations are restored in solution. Thus, solvation is a crucial ingredient that must be taken into account in studies aimed at breaking scaling relations in solution. Our findings suggest that the choice of metal and ligand determines the catalytic activity within the limits imposed by scaling relations, whereas the choice of an appropriate solvent can drive such activity beyond those limits.

Hydrocarbon adsorption in an aqueous environment: A computational study of alkyls on Cu(111)

Hydrocarbon chains are important intermediates in various aqueous-phase surface processes, such as CO2 electroreduction, aqueous Fischer-Tropsch synthesis, and aqueous phase reforming of biomass-derived molecules. Further, the interaction between water and adsorbed hydrocarbons represents a difficult case for modern computational methods. Here, we explore various methods for calculating the energetics of this interaction within the framework of density functional theory and explore trade-offs between the use of low water coverages, molecular dynamics approaches, and minima hopping for identification of low energy structures. An effective methodology for simulating low temperature processes is provided by using a unit cell in which the vacuum space is filled with water, employing the minima hopping algorithm to search for low-lying minima, and including dispersion (van der Waals) interactions. Using this methodology, we show that a high coverage of adsorbed alkyls is destabilized by the presence of water, while a low coverage of alkyls is stabilized. Solvation has a small effect on the energetics of hydrocarbon chain growth, generally decreasing its favorability at low temperatures. We studied higher temperatures by running molecular dynamics simulations starting at the minima found by the minima hopping algorithm and found that increased temperatures facilitate chain growth. The self-consistent continuum solvation method effectively describes the alkyl-water interaction and is in general agreement with the explicit solvation results in most cases, but care should be taken at high alkyl coverage.

Adsorption of asphaltenes on the calcite (10.4) surface by first-principles calculations

While the asphaltene-resin dimer leads to π–π stacking outward from surface, sulphide group favors this nanoaggregate adsorption on calcite that is energetically similar to the asphaltene monomer adsorbed through less steric and dielectric effects.

Catalytic activity of Pt38 in the oxygen reduction reaction from first-principles simulations

Mechanism of OH_ads/Pt₃₈ diffusion via transient hydronium species in first-principles molecular dynamics simulations.

Influence of nitrogen precursors on the structure, composition, and oxygen reduction reaction performance of dual heteroatom doped carbon nanohorns

The different N-sources used in the synthesis would cause different ways the doped N atoms interact with B or P dopants, and result in the enhanced or declined electrocatalytic activity of dual-doped carbon nanohorns in oxygen reduction reaction.

The atomistic origin of the extraordinary oxygen reduction activity of Pt3Ni7 fuel cell catalysts

Recently Debe et al. reported that Pt3Ni7 leads to extraordinary Oxygen Reduction Reaction (ORR) activity. However, several reports show that hardly any Ni remains in the layers of the catalysts close to the surface ("Pt-skin effect"). This paradox that Ni is essential to the high catalytic activity with the peak ORR activity at Pt3Ni7 while little or no Ni remains close to the surface is explained here using large-scale first-principles-based simulations. We make the radical assumption that processing Pt-Ni catalysts under ORR conditions would leach out all Ni accessible to the solvent. To simulate this process we use the ReaxFF reactive force field, starting with random alloy particles ranging from 50% Ni to 90% Ni and containing up to ∼300 000 atoms, deleting the Ni atoms, and equilibrating the resulting structures. We find that the Pt3Ni7 case and a final particle radius around 7.5 nm lead to internal voids in communication with the exterior, doubling the external surface footprint, in fair agreement with experiment. Then we examine the surface character of these nanoporous systems and find that a prominent feature in the surface of the de-alloyed particles is a rhombic structure involving 4 surface atoms which is crystalline-like but under-coordinated. Using density-functional theory, we calculate the energy barriers of ORR steps on Pt nanoporous catalysts, focusing on the Oad-hydration reaction (Oad + H2Oad → OHad + OHad) but including the barriers of O2 dissociation (O2ad → Oad + Oad) and water formation (OHad + Had → H2Oad). We find that the reaction barrier for the Oad-hydration rate-determining-step is reduced significantly on the de-alloyed surface sites compared to Pt(111). Moreover we find that these active sites are prevalent on the surface of particles de-alloyed from a Pt-Ni 30 : 70 initial composition. These simulations explain the peak in surface reactivity at Pt3Ni7, and provide a rational guide to use for further optimization of improved catalytic and nanoporous materials.

Optimizing the oxygen evolution reaction for electrochemical water oxidation by tuning solvent properties

Electrochemical water-based energy cycles provide a most promising alternative to fossil-fuel sources of energy. However, current electrocatalysts are not adequate (high overpotential, lack of selectivity toward O2 production, catalyst degradation). We propose here mechanistic guidelines for experimental examination of modified catalysts based on the dependence of kinetic rates on the solvent dielectric constant. To illustrate the procedure we consider the fcc(111) platinum surface and show that the individual steps for the oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) change systematically with the polarizability of the medium. Thus changing this environmental variable can be used to tune the rate determining steps and the barriers, providing a means for screening and validating new systems to optimize the rate determining steps for the ORR and OER reaction pathways.

A computationally designed titanium-mediated amination of allylic alcohols for the synthesis of secondary allylamines

A direct amination on allylic alcohols under mild conditions was enlightened by computational investigations and implemented in secondary allylamines synthesis.