Off-equilibrium reactivity of boron-enriched metal diboride surfaces in electro-reduction conditions

07 March 2025, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Boron-based material, featuring B-dependent reactivity and diverse phases, is emerging as a promising catalyst system. However, catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron enrichment strategy to tune the surface reactivity of the hexagonal face of MoB2. We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime, with an extensive examination on the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active towards HER via Volmer-Heyrovsky mechanism. This work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address them, along with other realistic complexities, in kinetics simulations.

Keywords

Grand Canonical Ensemble
Global Optimization
Ab initio Thermodynamics
Computational Catalysis
Surface Reconstruction
Adsorbate Coverage
Adsorption Configurations
Electrocatalysis
Hydrogen evolution reaction
Boron
Boride

Supplementary materials

Title
Description
Actions
Title
Supporting Information of Off-equilibrium reactivity of boron-enriched metal diboride surfaces in electroreduction conditions
Description
Slab models for the GCGA, GC-DFT, and AIMD simulations; Pourbaix diagram of boron in aqueous media; Bader charge analyses; Electron localization function analyses; Charge density difference analysis; Quantum theory of atoms in molecules analyses; Global and minima structures from global optimization searches; Progress plot of GCGA searches; Similarity map of sampled H-covered B 1/3 structures from GCGA searches; Surface atom corrugation under varying H coverage; Population of different H binding sites and modes under varying H coverage; Energy profiles of in-plane B rearrangement, H diffusion, and Tafel step from NEB calculations; PES of center B/B-H migration, dimerization, trimerization, and borane formation from 2D PES scans; Key AIMD snapshots showing delocalization of protons in hydronium/water layers and micro-solvation clusters; Fitting of decoupled PES to Morse potentials; Combined adiabatic FES of PCET and its response to varying electrode potential and reaction coordinate; Note on details of the GC global optimizations
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.