Density functional theory calculations for the hydrogen evolution reaction in an electrochemical double layer on the Pt(111) electrode

Insight into the catalytic activity of MXenes for hydrogen evolution reaction

MXenes have exhibited great potential as cost-effective electrocatalysts for hydrogen evolution reaction (HER). However, insight into the origin of activity is still missing. Herein, on the basis of a systematical investigation of the HER performance of 20 MXenes (M₂NO₂ and M₂CO₂, M = Sc, Ti, V, Cr, Zr, Nb, Mo, Hf, Ta and W), a Fermi-abundance model is proposed to understand variation of the activity in different MXenes. It is found that the occupied p electronic states of surface O atoms play a decisive role in the HER activity of MXenes. More importantly, Ti₂NO₂ and Nb₂NO₂ are found to be promising HER electrocatalysts with the free energy for hydrogen adsorption close to zero. This work not only provides possible catalysts for HER, the developed Fermi-abundance model but also is applicable to other two-dimensional materials and may serve as a simple descriptor of the intrinsic HER activity.

Effect of hydrophobic cations on the oxygen reduction reaction on single‒crystal platinum electrodes

Highly active catalysts for the oxygen reduction reaction are essential for the widespread and economically viable use of polymer electrolyte fuel cells. Here we report the oxygen reduction reaction activities of single‒crystal platinum electrodes in acidic solutions containing tetraalkylammonium cations with different alkyl chain lengths. The high hydrophobicity of a tetraalkylammonium cation with a longer alkyl chain enhances the oxygen reduction reaction activity. The activity on Pt(111) in the presence of tetra‒n‒hexylammonium cation is eight times as high as that without this cation, which is comparable to the activities on Pt₃Co(111) and Pt₃Ni(111) electrodes. Hydrophobic cations and their hydration shells destabilize the adsorbed hydroxide and adsorbed water. The hydrophobic characteristics of non‒specifically adsorbed cations can prevent the adsorption of poisoning species on the platinum electrode and form a highly efficient interface for the oxygen reduction reaction.

Highly active catalysts for the oxygen reduction reaction are valuable for fuel cells. Here the authors evaluate catalytic activity of single-crystal platinum electrodes in acidic solutions that contain hydrophobic cations, which prevent the adsorption of poisoning species and form an efficient interface.

Hydrogen evolution reaction (HER) on Au@Ag ultrananoclusters as electro-catalysts

We propose and investigate computationally Ag-Au subnanometer clusters as catalysts for the hydrogen evolution reaction (HER). Focusing on Ag12Au, we conduct a complete first-principles study of the HER process on this Au@Ag ultrananocatalyst. After determining the hydrogen-saturated resting state under standard conditions as Ag12AuH11, HER reaction energies and barriers are predicted also including solvent effects using both implicit and explicit models. We find that Ag12Au is a good candidate as a HER catalyst, with good stability and an overall reaction energy barrier of 0.89 eV as an upper bound. We also draw indications for the design of HER subnanometer catalysts.

OH formation and H2 adsorption at the liquid water-Pt(111) interface

The liquid water-Pt(111) interface is studied with constant temperature ab initio molecular dynamics to explore the importance of liquid water dynamics of catalytic reactions such as the oxygen reduction reaction in PEM fuel cells. The structure and energetics of hydroxyls formed at the liquid water-Pt(111) interface are found to be significantly different from those of the hydroxyl formed on a bare Pt(111) surface and the hydroxyl formed on a Pt(111) surface with a static water layer. We identify 1/12 ML *OH, 5/12 ML *OH and 2/3 ML *OH as particularly stable hydroxyl coverages in highly dynamic liquid water environments, which - contrary to static water-hydroxyl models - contain adjacent uncovered Pt sites. Atomic surface oxygen is found to be unstable in the presence of liquid water, in contrast to static atomic level simulations. These results give an improved understanding of hydroxide and surface oxide formation from Pt(111) cyclic voltammetry and allow us to draw detailed connections between the electrostatic potential and the interface structure. The study of hydrogen adsorption at the liquid water-Pt(111) interface finds competitive adsorption between the adsorbed hydrogen atoms and water molecules. This does not adhere with experimental observations, and this indicates that the Pt(111) surface has to be negatively charged for a correct description of the liquid water-Pt(111) interface at potentials where hydrogen adsorption occurs.

Fermi Level Determination for Charged Systems via Recursive Density of States Integration

Determining the Fermi level position for a given material is important to understand many of its electronic and chemical properties. Ab initio methods are effective in computing Fermi levels when using charge-neutral supercells. However, in the case where charges are explicitly included, the compensating homogeneous background charge, which is necessary to maintain charge neutrality in periodic models, causes the vacuum potential to be ill-defined - which would otherwise have been a reliable reference potential. Here, we develop a method based on recursively integrating the density of states to determine shifts in the Fermi level upon charging. By introducing incremental charges, one can compute the density of states profile and determine the shift in the Fermi level that corresponds to adding or removing a given increment of charge δq, which allows the evaluation of the Fermi level for any arbitrary charge q. We test this method for a range of materials (graphene, h-BN, C₃N₄, Cu, and MoS₂) and demonstrate that this method can produce a reasonable agreement with models that rely on localized compensating background charges.

Phase Modulation of (1T-2H)-MoSe2/TiC-C Shell/Core Arrays via Nitrogen Doping for Highly Efficient Hydrogen Evolution Reaction

Tailoring molybdenum selenide electrocatalysts with tunable phase and morphology is of great importance for advancement of hydrogen evolution reaction (HER). In this work, phase- and morphology-modulated N-doped MoSe2 /TiC-C shell/core arrays through a facile hydrothermal and postannealing treatment strategy are reported. Highly conductive TiC-C nanorod arrays serve as the backbone for MoSe2 nanosheets to form high-quality MoSe2 /TiC-C shell/core arrays. Impressively, continuous phase modulation of MoSe2 is realized on the MoSe2 /TiC-C arrays. Except for the pure 1T-MoSe2 and 2H-MoSe2 , mixed (1T-2H)-MoSe2 nanosheets are achieved in the N-MoSe2 by N doping and demonstrated by spherical aberration electron microscope. Plausible mechanism of phase transformation and different doping sites of N atom are proposed via theoretical calculation. The much smaller energy barrier, longer HSe bond length, and diminished bandgap endow N-MoSe2 /TiC-C arrays with substantially superior HER performance compared to 1T and 2H phase counterparts. Impressively, the designed N-MoSe2 /TiC-C arrays exhibit a low overpotential of 137 mV at a large current density of 100 mA cm-2 , and a small Tafel slope of 32 mV dec-1 . Our results pave the way to unravel the enhancement mechanism of HER on 2D transition metal dichalcogenides by N doping.

In-situ developed carbon spheres function as promising support for enhanced activity of cobalt oxide in oxygen evolution reaction

Highly active, stable electrocatalyst for oxygen evolution reaction (OER) is sincerely required for the practical application of water splitting to get rid from the sluggish reaction kinetics and the stability issue. Here, Co₃O₄ is studied as OER catalyst and to improve the electrocatalytic activity, carbon is chosen as the conducting support. A simple and cost-effective synthetic route is developed for the synthesis of Co₃O₄ on carbon support following hydrothermal route using various hydrolyzing agents. The heterostructure 'Co₃O₄/C' perform well in OER as a non-precious metal catalyst. The best Co₃O₄/C electrocatalyst can generate 10 and 30 mA/cm² current densities upon application of 1.623 V and 1.678 V vs. RHE whereas, bare Co₃O₄ can generate current density of 10 and 30 mA/cm² upon application of 1.677 and 1.754 V vs. RHE. Carbon in the heterostructure helps to improve the conductivity and at the same time enhances the charge transfer ability which further leads to increase current density and stability to the catalyst. Co₃O₄/C can generate unaltered current density up to 1000 cycles.

Tailoring Bandgap of Perovskite BaTiO₃ by Transition Metals Co-Doping for Visible-Light Photoelectrical Applications: A First-Principles Study

The physical and chemical properties of V-M″ and Nb-M″ (M″ is 3d or 4d transition metal) co-doped BaTiO₃ were studied by first-principles calculation based on density functional theory. Our calculation results show that V-M″ co-doping is more favorable than Nb-M″ co-doping in terms of narrowing the bandgap and increasing the visible-light absorption. In pure BaTiO₃, the bandgap depends on the energy levels of the Ti 3d and O 2p states. The appropriate co-doping can effectively manipulate the bandgap by introducing new energy levels interacting with those of the pure BaTiO₃. The optimal co-doping effect comes from the V-Cr co-doping system, which not only has smaller impurity formation energy, but also significantly reduces the bandgap. Detailed analysis of the density of states, band structure, and charge-density distribution in the doping systems demonstrates the synergistic effect induced by the V and Cr co-doping. The results can provide not only useful insights into the understanding of the bandgap engineering by element doping, but also beneficial guidance to the experimental study of BaTiO₃ for visible-light photoelectrical applications.

New Insight into the Mechanism of the Hydrogen Evolution Reaction on MoP(001) from First Principles

Molybdenum phosphide-based catalysts have recently exhibited excellent catalytic activities for the hydrogen evolution reaction (HER) in wide pH range conditions; the intrinsic reaction mechanism, on the other hand, has not been well established. Herein, by employing the MoP as the prototypical molybdenum phosphide-based catalyst, HER activities in both acid and neutral conditions were studied by conducting periodic density functional theory calculations. Thermodynamic analysis of hydrogen atoms absorbed on both P- and Mo-terminated surfaces were compared, as well as all the reaction energy and activation energy barriers for reactions involved in the HER process. Calculation results revealed that, in an acid condition, the Volmer-Heyrovsky and Volmer-Tafel reaction mechanisms were dominated on the P-terminated and Mo-terminated catalyst surfaces, where Heyrovsky and Volmer reactions were the rate-determining step, respectively. Additionally, water splitting was introduced to the current reaction mechanism and a small reaction activation energy barrier was revealed on the P-terminated surface. Besides, a relevant small activation energy was obtained in the Tafel reaction on the defect of the P-terminated surface in a neutral solution. Theoretical results proved that HER could take place readily on both P- and Mo-terminated catalyst surfaces via different reaction mechanisms in the acid condition from the view of atom scale. More important, computational results uncovered that HER could also occur on the P-terminated surface with the assistance of surface defect in the neutral condition, which sheds new light on the HER mechanism on transition metal phosphite-based catalysts. The doping effect on HER activity was further investigated in theory and calculation results, indicating that catalytic performance could be improved by substitutional doping of the Mo atom with metals such as Mn and W.

First-Principles Approach to Model Electrochemical Reactions: Understanding the Fundamental Mechanisms behind Mg Corrosion

Combining concepts of semiconductor physics and corrosion science, we develop a novel approach that allows us to perform ab initio calculations under controlled potentiostat conditions for electrochemical systems. The proposed approach can be straightforwardly applied in standard density functional theory codes. To demonstrate the performance and the opportunities opened by this approach, we study the chemical reactions that take place during initial corrosion at the water-Mg interface under anodic polarization. Based on this insight, we derive an atomistic model that explains the origin of the anodic hydrogen evolution.

Heterostructures of MXenes and N-doped graphene as highly active bifunctional electrocatalysts

MXenes with versatile chemistry and superior electrical conductivity are prevalent candidate materials for energy storage and catalysts. Inspired by recent experiments of hybridizing MXenes with carbon materials, here we theoretically design a series of heterostructures of N-doped graphene supported by MXene monolayers as bifunctional electrocatalysts for the oxygen reduction reaction (ORR) and hydrogen evolution reaction (HER). Our first-principles calculations show that the graphitic sheet on V2C and Mo2C MXenes are highly active with an ORR overpotential down to 0.36 V and reaction free energies for the HER approaching zero, both with low kinetic barriers. Such outstanding catalytic activities originate from the electronic coupling between the graphitic sheet and the MXene, and can be correlated with the pz band center of surface carbon atoms and the work function of the heterostructures. Our findings screen a novel form of highly active electrocatalysts by taking advantage of the fast charge transfer kinetics and strong interfacial coupling of MXenes, and illuminate a universal mechanism for modulating the catalytic properties of two-dimensional hybrid materials.

Atomic-Scale Simulation of Electrochemical Processes at Electrode/Water Interfaces under Referenced Bias Potential

Based on constant Fermi-level molecular dynamics and a proper alignment scheme, we perform simulations of the Pt(111)/water interface under variable bias potential referenced to the standard hydrogen electrode (SHE). Our scheme yields a potential of zero charge μ_pzc of ∼0.22 eV relative to the SHE and a double layer capacitance C_dl of ≃19 μF cm^-2, in excellent agreement with experimental measurements. In addition, we study the structural reorganization of the electrical double layer for bias potentials ranging from -0.92 eV to +0.44 eV and find that O_down configurations, which are dominant at potentials above the pzc, reorient to favor H_down configurations as the measured potential becomes negative. Our modeling scheme allows one to not only access atomic-scale processes at metal/water interfaces, but also to quantitatively estimate macroscopic electrochemical quantities.

The stability limits of highly active nitrogen doped carbon ORR nano-catalysts: a mechanistic study of degradation reactions

A new approach in electrode catalysis bearing immense potential for electrochemical technologies is the prospect of carbon-based electrodes. Pristine carbon nanostructures are relatively inert and modifications like nitrogen doping are known for their beneficial effects on the electrochemical activity of carbon nanomaterials in both alkaline and acidic media. However, the long-term stability of these materials, especially in an acidic environment, is rarely mentioned. Here, we evaluate the stability and long-term degradation of nitrogen doped graphene flakes as an oxygen reduction electrocatalyst with theoretical and experimental techniques. We assume that nitrogen dopants in the graphene sheet interact with e- and H+ at the electrode-electrolyte interface, leading to NH3 scission and continuous catalyst deactivation. With Density Functional Theory calculations, NH3 scission pathways of pyridinic, graphitic and pyrrolic nitrogen species were analyzed and compared under different operating conditions which are relevant for low and intermediate temperature fuel cells. The computational results are correlated with electrochemical measurements in solid acid fuel cells in a humidified oxygen environment at 240 °C.

Hydrogen bonding effect between active site and protein environment on catalysis performance in H2-producing [NiFe] hydrogenases

The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol^-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H⁺/e^- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S^- in cysteine residues shows a larger influence on the energy profile than that via the CN^- ligands on Fe.

A theoretical study of the effect of a non-aqueous proton donor on electrochemical ammonia synthesis

Ammonia synthesis is one of the most studied reactions in heterogeneous catalysis. To date, however, electrochemical N₂ reduction in aqueous systems has proven to be extremely difficult, mainly due to the competing hydrogen evolution reaction (HER). Recently, it has been shown that transition metal complexes based on molybdenum can reduce N₂ to ammonia at room temperature and ambient pressure in a non-aqueous system, with a relatively small amount of hydrogen output. We demonstrate that the non-aqueous proton donor they have chosen, 2,6-lutidinium (LutH⁺), is a viable substitute for hydronium in the electrochemical process at a solid surface, since this donor can suppress the HER rate. We also show that the presence of LutH⁺ can selectively stabilize the *NNH intermediate relative to *NH or *NH₂via the formation of hydrogen bonds, indicating that the use of non-aqueous solvents can break the scaling relationship between limiting potential and binding energies.

A kinetic model of water adsorption, clustering and dissociation on the Fe3S4{001} surface

The interaction of water with catalyst surfaces is a common process which requires investigation. Here, we have employed density functional theory calculations to investigate the adsorption of up to ten water molecules on the {001} surface of greigite (Fe₃S₄), which owing to its redox properties, is of increasing interest as a catalyst, e.g. in electro-catalysis. We have systematically analyzed and characterized the modes of water adsorption on the surface, where we have considered both molecular and dissociative adsorption processes. The calculations show that molecular adsorption is the predominant state on these surfaces, from both a thermodynamic and kinetic point of view. We have explored the molecular dispersion on the surface under different coverages and found that the orientation of the molecule, and therefore the surface dipole, depends on the number of adsorbed molecules. The interactions between the water molecules become stronger with an increasing number of water molecules, following an exponential decay which tends to the interaction energy found in bulk water. We have also shown the evolution of the infra-red signals as a function of water coverage relating to the H-bond networks formed on the surface. Next we have included these results in a classical micro-kinetic model, which introduced the effects of temperature in the simulations, thus helping us to derive the water cluster size on the greigite surface as a function of the initial conditions of pressure, temperature and external potential. The kinetic model concluded that water molecules agglomerate in clusters instead of wetting the surface, which agrees with the low hydrophilicity of Fe₃S₄. Clusters consisting of four water molecules was shown to be the most stable cluster under a wide range of temperatures and external potential.

Modeling the active sites of Co-promoted MoS2 particles by DFT

The atomic-scale structure of the Co-promoted MoS₂ catalyst (CoMoS), used for hydrodesulfurization and as a potential replacement for platinum in the acidic hydrogen evolution reaction has been analyzed by modeling its sites using density functional theory and applying thermochemical corrections to account for different reaction conditions. The equilibrium structures of the edges, basal plane and corners have been found and used to obtain a picture of an ideal CoMoS nanoparticle under hydrodesulfurization and hydrogen evolution reaction conditions. Under hydrodesulfurization conditions small energy differences between structures having an additional or missing sulfur atom relative to the equilibrium structures have been observed for the edges and corners explaining their activity towards hydrodesulfurization at the atomic scale. The lack of these small energy differences at the basal plane explains why it is inert towards hydrodesulfurization. The adsorption free energy of hydrogen was calculated and used as a descriptor for qualifying each site in the context of hydrogen evolution, finding that the corner site should perform better than the edges.

Removal of Cr(iii)/Cr(vi) from wastewater using defective porous boron nitride: a DFT study

In this work, by using DFT calculation, we revealed that p-BN with vacancy defects could efficiently remove both Cr(iii) and Cr(vi) ions from an aqueous solution due to large E_ads values even in the higher Cr(iii)/Cr(vi) coverage.

Hydrogen adsorption trends on Al-doped Ni2P surfaces for optimal catalyst design

Nanoparticles of nickel phosphide are promising materials to replace the currently used rare Pt-group metals at cathode-side electrodes in devices for electrochemical hydrogen production.

MXene nanoribbons as electrocatalysts for the hydrogen evolution reaction with fast kinetics

MXenes, a new class of two-dimensional materials, arouse great interest due to their diverse chemistries, superior electrical conductivity and stability.

A post Gurney quantum mechanical perspective on the electrolysis of water: ion neutralization in solution

Electron fluxes crossing the interface between a metallic conductor and an aqueous environment are important in many fields; hydrogen production, environmental scanning tunnelling microscopy, scanning electrochemical microscopy being some of them. Gurney (Gurney 1931 Proc. R. Soc. Lond.134, 137 (doi:10.1098/rspa.1931.0187)) provided in 1931 a scheme for tunnelling during electrolysis and outlined conditions for it to occur. We measure the low-voltage current flows between gold electrodes in pure water and use the time-dependent behaviour at voltage switch-on and switch-off to evaluate the relative contribution to the steady current arising from tunnelling of electrons between the electrodes and ions in solution and from the neutralization of ions adsorbed onto the electrode surface. We ascribe the larger current contribution to quantum tunnelling of electrons to and from ions in solution near the electrodes. We refine Gurney's barrier scheme to include solvated electron states and quantify energy differences using updated information. We show that Gurney's conditions would prevent the current flow at low voltages we observe but outline how the ideas of Marcus (Marcus 1956 J. Chem. Phys.24, 966-978 (doi:10.1063/1.1742723)) concerning solvation fluctuations enable the condition to be relaxed. We derive an average barrier tunnelling model and a multiple pathways tunnelling model and compare predictions with measurements of the steady-state current-voltage relation. The tunnelling barrier was found to be wide and low in agreement with other experimental studies. Applications as a biosensing mechanism are discussed that exploit the fast tunnelling pathways along molecules in solution.

Nature of Highly Active Electrocatalytic Sites for the Hydrogen Evolution Reaction at Pt Electrodes in Acidic Media

The hydrogen evolution reaction (HER) is one of the two processes in electrolytic water splitting. Known for more than two centuries, the HER still receives great attention in fundamental and applied science in view of its apparent simplicity (only two electrons are transferred), fast kinetics in acidic media, and promising technological applications in electrolyzers. However, the exact nature of active catalytic sites for this reaction is often uncertain, especially at nonuniform metal electrodes. Identification of such centers is important, as the HER will probably be central in future energy provision schemes, and it is simultaneously a convenient model reaction to study structure-composition-activity relations in catalysis. In this work, using simple coordination-activity considerations, we outline the location and geometric configuration of the active sites at various model Pt single-crystal electrodes. We show that when the coordination of such surface sites is optimized and their density at the surface is maximized, the experimental-specific HER activities are among the highest reported in the literature for pure platinum with a well-defined surface structure under similar conditions.

Hydrogen adsorption on doped MoS2 nanostructures

Electrochemical devices for efficient production of hydrogen as energy carrier rely still largely on rare platinum group metal catalysts. Chemically and structurally modified metal dichalcogenide MoS₂ is a promising substitute for these critical raw materials at the cathode side where the hydrogen evolution reaction takes place. For precise understanding of structure and hydrogen adsorption characteristics in chemically modified MoS₂ nanostructures, we perform comprehensive density functional theory calculations on transition metal (Fe, Co, Ni, Cu) doping at the experimentally relevant MoS₂ surfaces at substitutional Mo-sites. Clear benefits of doping the basal plane are found, whereas at the Mo- and S-edges complex modifications at the whole edge are observed. New insight into doping-enhanced activity is obtained and guidance is given for further experiments. We study a machine learning model to facilitate the screening of suitable structures and find a promising level of prediction accuracy with minimal structural input.

Highly Active Tungsten Oxide Nanoplate Electrocatalysts for the Hydrogen Evolution Reaction in Acidic and Near Neutral Electrolytes

An efficient, cost-effective, and earth-abundant catalyst that could drive the production of hydrogen from water without or with little external energy is the ultimate goal toward hydrogen economy. Herein, nanoplates of tungsten oxide and its hydrates (WO₃·H₂O) as promising electrocatalysts for the hydrogen evolution reaction (HER) are reported. The square-shaped and stacked WO₃·H₂O nanoplates are synthesized at room temperature under air in ethanol only, making it as a promising green synthesis strategy. The repeated electrochemical cyclic voltammetry cycles modified the surface of WO₃·H₂O nanoplates to WO₃ as confirmed by X-ray photoelectron and Auger spectroscopy, which leads to an improved HER activity. Hydrogen evolution is further achieved from distilled water (pH 5.67) producing 1 mA cm^-2 at an overpotential of 15 mV versus the reversible hydrogen electrode. Moreover, WO₃·H₂O and WO₃ nanoplates demonstrate excellent durability in acidic and neutral media, which is highly desirable for practical application. Improved hydrogen evolution by WO₃(200) when compared to that by Pt(111) is further substantiated by the density functional theory calculations.