Mechanistic Study of the Electrochemical Oxygen Reduction Reaction on Pt(111) Using Density Functional Theory

Hydrogen Coupling on Platinum Using Artificial Neural Network Potentials and DFT

To date, the understanding of reactions at solid-liquid interfaces has proven challenging, mainly because of the inaccessible nature of such systems to current experimental techniques with atomic resolution. This has meant that many important features, including free energy barriers and the atomistic structure of intermediates, remain unknown. To tackle these issues, we construct and utilize a high-dimensional neural network (HDNN) potential for the simulation of hydrogen evolution at the HCl(aq)/Pt(111) interface, taking into consideration the influence of adsorbate-adsorbate, adsorbate-solvent interactions, and ion solvation explicitly. Long time scale MD simulations reveal coadsorbed H_ad/H₂O_ad on the surface. The free energy profiles for the Tafel and Heyrovsky type hydrogen coupling are extracted using umbrella sampling. It is found that the preferential mechanism can change depending on the surface coverage, highlighting the dual mechanistic nature for HER on Pt(111). Our work demonstrates the importance of controlling the solvent-substrate interactions in developing catalysts beyond Pt.

Cross-Sphere Electrode Reaction: The Case of Hydroxyl Desorption during the Oxygen Reduction Reaction on Pt(111) in Alkaline Media

Hydroxide ion is a common electrolyte when electrode reactions take place in alkaline media. In the case of oxygen reduction reaction on Pt(111), we demonstrate by ab initio molecular dynamics calculations that the desorption of hydroxyl (OH*) from the electrode surface to form a solvated OH^- is a cross-sphere process, with the OH* reactant in the inner sphere and the OH^- product directly generated in the aqueous outer sphere. Such a mechanism is distinct from the typical inner sphere and outer sphere reactions. It is dictated by the strong hydrogen bonding interactions between a hydroxide ion and water molecules and is facilitated by proton transfer through solvation layers. It should play a significant role whenever OH* desorption, or its reverse, OH^- adsorption, is involved in an electrochemical reaction.

Electrocatalysis for CO2 conversion: from fundamentals to value-added products

This timely and comprehensive review mainly summarizes advances in heterogeneous electroreduction of CO₂: from fundamentals to value-added products.

Explainable and trustworthy artificial intelligence for correctable modeling in chemical sciences

Data science has primarily focused on big data, but for many physics, chemistry, and engineering applications, data are often small, correlated and, thus, low dimensional, and sourced from both computations and experiments with various levels of noise. Typical statistics and machine learning methods do not work for these cases. Expert knowledge is essential, but a systematic framework for incorporating it into physics-based models under uncertainty is lacking. Here, we develop a mathematical and computational framework for probabilistic artificial intelligence (AI)-based predictive modeling combining data, expert knowledge, multiscale models, and information theory through uncertainty quantification and probabilistic graphical models (PGMs). We apply PGMs to chemistry specifically and develop predictive guarantees for PGMs generally. Our proposed framework, combining AI and uncertainty quantification, provides explainable results leading to correctable and, eventually, trustworthy models. The proposed framework is demonstrated on a microkinetic model of the oxygen reduction reaction.

DFT Calculations of the Adsorption States of O2 on OH/H2O-Covered Pt(111)

The adsorption of O₂ on Pt(111) was studied with Density Functional Theory calculations. Various adsorbed states of O₂ were evaluated on clean and OH/H₂O-covered Pt(111) surfaces at the solid/gas and solid/liquid interfaces. The results reveal that the adsorption of O₂ on OH/H₂O-covered Pt(111) surface starts with the physical adsorption of O₂. Two other adsorption states are reachable from the physisorbed state, the end-on, and bridging chemisorbed O₂. Analysis of the energetics of these adsorption states shows that O₂ physically adsorbed at the OH/H₂O-covered Pt( 111) surface is a high energy state that requires activation to transition to the end-on chemisorbed O₂ state. On the other hand, the end-on chemisorbed state can transition to the bridging chemisorbed state with only a small activation energy when a nearby Pt adsorption site is available. Frequency analysis of the physisorbed, end-on, and bridging adsorption states shows that adsorbed O₂ stretching frequencies are close to 1400, 1300, and 900 cm^-1, respectively.

Structural investigation of ternary PdRuM (M = Pt, Rh, or Ir) nanoparticles using first-principles calculations

We perform first-principles calculations and Monte Carlo sampling to investigate the structures of ternary PdRuM (M = Pt, Rh, or Ir) nanoparticles (NPs) with respect to three different spherical shapes. The morphologies include hexagonal close-packed (hcp), truncated-octahedral (fcc), and icosahedral (Ih, fcc) shapes with 57, 55, and 55 atoms, respectively. The calculations show that the atomic position is dominant in determining the stability of the ternary NPs. For bare ternary NPs, Pd and Ru atoms favor a location on the vertex sites and the core, respectively, which can be understood by the surface energy of the corresponding slab models. For single-crystalline NPs, the binary shell could be either a solid solution or a segregation alloy depending on composition and morphology. However, polycrystalline Ih NPs only form segregated binary shells surrounding the Ru core. Such configurations tend to minimize the surface lattice to gain more energy from the d orbital of the transition metals. In addition to the bare NPs, we study the oxidized ternary NPs. The results show that the Ru atoms penetrate outwards from the core to the surface reducing the oxidation formation energy. Furthermore, oxygen adsorption facilitates Pt, Pd, and Pd penetration into the PdRuPt, PdRuRh, and PdRuIr NPs, respectively. Most of the oxide shells are a solid solution, except for the PdRuRh NP with an Ih shape, which is found to be in a segregation shell. The free energy calculation reveals that the pure hcp NPs are thermodynamically unstable under oxygen-rich conditions. This work clearly demonstrates the structural trends of small ternary NPs and their oxidation, unveiling that the structural trends can be understood by the surface formation energy and the interplay between adsorbent and adsorbing oxygen atoms.

Solvated proton and the origin of the high onset overpotential in the oxygen reduction reaction on Pt(111)

For the hydrogenation of O atoms on Pt(111), protonation can be bypassed by hydrolysis as the electrode potential rises.

Insights into how the aqueous environment influences the kinetics and mechanisms of heterogeneously-catalyzed COH* and CH3OH* dehydrogenation reactions on Pt(111)

Water influences catalytic reactions in multiple ways, including energetic and mechanistic effects. While simulations have provided significant insight into the roles that H2O molecules play in aqueous-phase heterogeneous catalysis, questions still remain as to the extent to which H2O structures influence catalytic mechanisms. Specifically, influences of the configurational variability in the water structures at the catalyst interface are yet to be understood. Configurational variability is challenging to capture, as it requires multiscale approaches. Herein, we apply a multiscale sampling approach to calculate reaction thermodynamics and kinetics for COH* dehydrogenation to CO* and CH3OH* dehydrogenation to CH2OH* on Pt(111) catalysts under liquid H2O. We explore various pathways for these dehydrogenation reactions that could influence the overall mechanism of methanol decomposition by including participation of H2O structures both energetically and mechanistically. We find that the liquid H2O environment significantly influences the mechanism of COH* dehydrogenation to CO* but leaves the mechanism of CH3OH* dehydrogenation to CH2OH* largely unaltered.

Recent advancements in Pt-nanostructure-based electrocatalysts for the oxygen reduction reaction

A comprehensive evaluation of Pt-nanostructure-based electrocatalysts for the oxygen reduction reaction.

Interconversion of hydrated protons at the interface between liquid water and platinum

The free energy barriers for hydrogen transfer at the H₂O/Pt(111) interface calculated usingab initiomolecular dynamics and umbrella sampling.

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.

A density functional theory study on oxygen reduction reaction on nitrogen-doped graphene

Nitrogen (N)-doped carbons reportedly exhibit good electrocatalytic activity for the oxygen reduction reaction (ORR) of fuel cells. This work provides theoretical insights into the ORR mechanism of N-doped graphene by using density functional theory calculations. All possible reaction pathways were investigated, and the transition state of each elementary step was identified. The results showed that OOH reduction was easier than O-OH breaking. OOH reduction followed a direct Eley-Rideal mechanism (the OOH species was in gas phase, but H was chemisorbed on the surface) with a significantly low reaction barrier of 0.09 eV. Pathways for both four-electron and two-electron reductions were possible. The rate-determining step of the two-electron pathway was the reduction of O₂ (formation of OOH), whereas that of the four-electron pathway was the reduction of OH into H₂O. After comparing the barriers of the rate-determining steps of the two pathways, we found that the two-electron pathway was more energetically favored than the four-electron pathway.

77Se solid-state NMR of As2Se3, As4Se4 and As4Se3 crystals: a combined experimental and computational study

(77)Se NMR parameters for three prototypical crystalline compounds (As2Se3, As4Se4 and As4Se3) have been determined from solid-state NMR spectra in the framework of an investigation concerning AsxSe(1-x) glass structure understanding. Density functional NMR calculations using the gauge including projector augmented wave methodology have been performed on X-ray and optimized crystal structures for a set of selenium-based crystals. These theoretical results have been combined with the experimental data in order to achieve a precise assignment of the spectral lines. This work and the high sensitivity of solid-state NMR to local order show that the structure of As4Se3 should be reinvestigated using state-of-the-art diffraction techniques. Calculations performed on several molecules derived from the crystal structures have demonstrated the limited effect of interlayer or intermolecular interactions on the isotropic chemical shifts. These interactions are therefore not responsible for the unexpected large chemical shift difference observed between these three systems that could mostly be attributed to the presence of short rings.

Prediction of standard interfacial electron-transfer rate constants for the Ru(NH3)6(3+/2+) couple through omega-hydroxyalkanethiol self-assembled monolayers on gold electrodes

Two relatively simple approaches are developed and used to calculate (predict) the standard interfacial electron-transfer (ET) rate constants (k degrees) of the Ru(NH3)6(3+/2+) couple dissolved in aqueous electrolyte solutions in contact with Au electrodes coated with self-assembled monolayers (SAMs) composed of HS(CH2)nOH as functions of both n and temperature. These approaches are suggested by the conclusion reached by Smalley et al. (J. Electroanal. Chem. 2006, 589, 1-6) that the interfacial ET rate of a solution-dissolved redox couple in contact with a SAM is, within 1 order of magnitude, the same as the (normalized) interfacial ET rate of a similar attached (as a constituent of a similar SAM) couple. The calculations, therefore, employ the measured electronic coupling of the attached (to Au electrodes through alkanethiolate bridges) -PyRu(NH3)5(3+/2+) couple. The two approaches also both include dynamic solvent effects on the ET kinetics and the influence of electronic coupling on the activation barrier for the ET reaction. At T=298 K and n=3, 11, and 14, the predicted rate constants are in very good agreement with the existing measurements of k degrees. However, for n<3 at 298 K, the predicted rate constants are extremely large (i.e., >4.5 cm s(-1)) and do not tend toward a limiting value. Additionally, even if the electronic coupling between a Au electrode and a Ru(NH3)6(3+/2+) moiety located at the surface of the SAM is >0.1 eV, the calculated standard rate constant is not directly proportional to the inverse of the longitudinal dielectric time of the solvent. A primary reason for both the absence of a limiting value for the predicted k degrees's at 298 K and the attenuated influence of dynamic solvent effects is the activation energy barrier suppression caused by large values of the electronic coupling.

Electrocatalysis of oxygen reduction and small alcohol oxidation in alkaline media

We present here a critical review of several technologically important electrocatalytic systems operating in alkaline electrolytes. These include the oxygen reduction reaction (ORR) occurring on catalysts containing Pt, Pd, Ir, Ru, or Ag, the methanol oxidation reaction (MOR) occurring on Pt-containing catalysts, and the ethanol oxidation reaction (EOR) occurring on Ni-Co-Fe alloy catalysts. Each of these catalytic systems is relevant to alkaline fuel cell (AFC) technology, while the ORR systems are also relevant to chlor-alkali electrolysis and metal-air batteries. The use of alkaline media presents advantages both in electrocatalytic activity and in materials stability and corrosion. Therefore, prospects for the continued development of alkaline electrocatalytic systems, including alkaline fuel cells, seem very promising.