skip to main content


Title: Revisiting trends in the exchange current for hydrogen evolution
Nørskov and collaborators proposed a simple kinetic model to explain the volcano relation for the hydrogen evolution reaction on transition metal surfaces such that j 0 = k 0 f (Δ G H ) where j 0 is the exchange current density, f (Δ G H ) is a function of the hydrogen adsorption free energy Δ G H as computed from density functional theory, and k 0 is a universal rate constant. Herein, focusing on the hydrogen evolution reaction in acidic medium, we revisit the original experimental data and find that the fidelity of this kinetic model can be significantly improved by invoking metal-dependence on k 0 such that the logarithm of k 0 linearly depends on the absolute value of Δ G H . We further confirm this relationship using additional experimental data points obtained from a critical review of the available literature. Our analyses show that the new model decreases the discrepancy between calculated and experimental exchange current density values by up to four orders of magnitude. Furthermore, we show the model can be further improved using machine learning and statistical inference methods that integrate additional material properties.  more » « less
Award ID(s):
1809085
PAR ID:
10417862
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Catalysis Science & Technology
Volume:
11
Issue:
20
ISSN:
2044-4753
Page Range / eLocation ID:
6832 to 6838
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Solid‐state electrocatalysis plays a crucial role in the development of renewable energy to reshape current and future energy needs. However, finding an inexpensive and highly active catalyst to replace precious metals remains a big challenge for this technology. Here, tri‐molybdenum phosphide (Mo3P) is found as a promising nonprecious metal and earth‐abundant candidate with outstanding catalytic properties that can be used for electrocatalytic processes. The catalytic performance of Mo3P nanoparticles is tested in the hydrogen evolution reaction (HER). The results indicate an onset potential of as low as 21 mV, H2formation rate, and exchange current density of 214.7 µmol s−1g−1cat(at only 100 mV overpotential) and 279.07 µA cm−2, respectively, which are among the closest values yet observed to platinum. Combined atomic‐scale characterizations and computational studies confirm that high density of molybdenum (Mo) active sites at the surface with superior intrinsic electronic properties are mainly responsible for the remarkable HER performance. The density functional theory calculation results also confirm that the exceptional performance of Mo3P is due to neutral Gibbs free energy (ΔGH*) of the hydrogen (H) adsorption at above 1/2 monolayer (ML) coverage of the (110) surface, exceeding the performance of existing non‐noble metal catalysts for HER.

     
    more » « less
  2. null (Ed.)
    By coupling a newly developed quantum-electronic-state-selected supersonically cooled vanadium cation (V + ) beam source with a double quadrupole-double octopole (DQDO) ion–molecule reaction apparatus, we have investigated detailed absolute integral cross sections ( σ 's) for the reactions, V + [a 5 D J ( J = 0, 2), a 5 F J ( J = 1, 2), and a 3 F J ( J = 2, 3)] + CH 4 , covering the center-of-mass collision energy range of E cm = 0.1–10.0 eV. Three product channels, VH + + CH 3 , VCH 2 + + H 2 , and VCH 3 + + H, are unambiguously identified based on E cm -threshold measurements. No J -dependences for the σ curves ( σ versus E cm plots) of individual electronic states are discernible, which may indicate that the spin–orbit coupling is weak and has little effect on chemical reactivity. For all three product channels, the maximum σ values for the triplet a 3 F J state [ σ (a 3 F J )] are found to be more than ten times larger than those for the quintet σ (a 5 D J ) and σ (a 5 F J ) states, showing that a reaction mechanism favoring the conservation of total electron spin. Without performing a detailed theoretical study, we have tentatively interpreted that a weak quintet-to-triplet spin crossing is operative for the activation reaction. The σ (a 5 D 0 , a 5 F 1 , and a 3 F 2) measurements for the VH + , VCH 2 + , and VCH 3 + product ion channels along with accounting of the kinetic energy distribution due to the thermal broadening effect for CH 4 have allowed the determination of the 0 K bond dissociation energies: D 0 (V + –H) = 2.02 (0.05) eV, D 0 (V + –CH 2 ) = 3.40 (0.07) eV, and D 0 (V + –CH 3 ) = 2.07 (0.09) eV. Detailed branching ratios of product ion channels for the titled reaction have also been reported. Excellent simulations of the σ curves obtained previously for V + generated by surface ionization at 1800–2200 K can be achieved by the linear combination of the σ (a 5 D J , a 5 F J , and a 3 F J ) curves weighted by the corresponding Boltzmann populations of the electronic states. In addition to serving as a strong validation of the thermal equilibrium assumption for the populations of the V + electronic states in the hot filament ionization source, the agreement between these results also confirmed that the V + (a 5 D J , a 5 F J , and a 3 F J ) states prepared in this experiment are in single spin–orbit states with 100% purity. 
    more » « less
  3. Abstract

    Cold H+produced via charge exchange reactions between ring current ions and exospheric neutral hydrogen constitutes an additional source of cold plasma that further contributes to the plasmasphere and affects the plasma dynamics in the Earth's magnetosphere system; however, its production and associated effects on the plasmasphere dynamics have not been fully assessed and quantified. In this study, we perform numerical simulations mimicking an idealized three‐phase geomagnetic storm to investigate the role of heavy ion composition in the ring current (O+vs. N+) and exospheric neutral hydrogen density in the production of cold H+via charge exchange reactions. It is found that ring current heavy ions produce more than 50% of the total cold H+via charge exchange reactions, and energetic N+is more efficient in producing cold H+via charge exchange reactions than O+. Furthermore, the density structure of the cold H+is highly dependent on the mass of the parent ion; that is, cold H+deriving from charge exchange reactions involving energetic O+with neutral hydrogen, populates the lower L‐shells, while cold H+deriving from charge exchange reactions involving energetic N+with neutral hydrogen populates the higher L‐shells. In addition, the density of cold H+produced via charge exchange reactions involving N+can be peak at values up to one order of magnitude larger than the local plasmaspheric density, suggesting that solely considering the supply of cold plasma from the ionosphere to the plasmasphere can lead to a significant underestimation of plasmasphere density.

     
    more » « less
  4. We generalize Hermite interpolation with error correction, which is the methodology for multiplicity algebraic error correction codes, to Hermite interpolation of a rational function over a field K from function and function derivative values. We present an interpolation algorithm that can locate and correct <= E errors at distinct arguments y in K where at least one of the values or values of a derivative is incorrect. The upper bound E for the number of such y is input. Our algorithm sufficiently oversamples the rational function to guarantee a unique interpolant. We sample (f/g)^(j)(y[i]) for 0 <= j <= L[i], 1 <= i <= n, y[i] distinct, where (f/g)^(j) is the j-th derivative of the rational function f/g, f, g in K[x], GCD(f,g)=1, g <= 0, and where N = (L[1]+1)+...+(L[n]+1) >= C + D + 1 + 2(L[1]+1) + ... + 2(L[E]+1) where C is an upper bound for deg(f) and D an upper bound for deg(g), which are input to our algorithm. The arguments y[i] can be poles, which is truly or falsely indicated by a function value infinity with the corresponding L[i]=0. Our results remain valid for fields K of characteristic >= 1 + max L[i]. Our algorithm has the same asymptotic arithmetic complexity as that for classical Hermite interpolation, namely soft-O(N). For polynomials, that is, g=1, and a uniform derivative profile L[1] = ... = L[n], our algorithm specializes to the univariate multiplicity code decoder that is based on the 1986 Welch-Berlekamp algorithm. 
    more » « less
  5. Understanding H 2 binding and activation is important in the context of designing transition metal catalysts for many processes, including hydrogenation and the interconversion of H 2 with protons and electrons. This work reports the first thermodynamic and kinetic H 2 binding studies for an isostructural series of first-row metal complexes: NiML, where M = Al ( 1 ), Ga ( 2 ), and In ( 3 ), and L = [N( o -(NCH 2 P i Pr 2 )C 6 H 4 ) 3 ] 3− . Thermodynamic free energies (Δ G °) and free energies of activation (Δ G ‡ ) for binding equilibria were obtained via variable-temperature 31 P NMR studies and lineshape analysis. The supporting metal exerts a large influence on the thermodynamic favorability of both H 2 and N 2 binding to Ni, with Δ G ° values for H 2 binding found to span nearly the entire range of previous reports. The non-classical H 2 adduct, (η 2 -H 2 )NiInL ( 3 -H 2 ), was structurally characterized by single-crystal neutron diffraction—the first such study for a Ni(η 2 -H 2 ) complex or any d 10 M(η 2 -H 2 ) complex. UV-Vis studies and TD-DFT calculations identified specific electronic structure perturbations of the supporting metal which poise NiML complexes for small-molecule binding. ETS-NOCV calculations indicate that H 2 binding primarily occurs via H–H σ-donation to the Ni 4p z -based LUMO, which is proposed to become energetically accessible as the Ni(0)→M( iii ) dative interaction increases for the larger M( iii ) ions. Linear free-energy relationships are discussed, with the activation barrier for H 2 binding (Δ G ‡ ) found to decrease proportionally for more thermodynamically favorable equilibria. The Δ G ° values for H 2 and N 2 binding to NiML complexes were also found to be more exergonic for the larger M( iii ) ions. 
    more » « less