Controlling HER activity and stability of $γ$- and 6,6,12-Graphyne through engineered B-N doping: DFT and Reactive MD simulations

Juan Gomez Quispe, Matheus Medina, Subhendu Mishra, Douglas S Galvao, Abhishek Singh, Pedro Alves da Silva Autreto

arXiv:2601.15424·cond-mat.mtrl-sci·Published 2026-01-21

Graphynes offer a chemically heterogeneous $sp/sp^{2}$ carbon framework with distinct electronic regimes and site-selective reactivity. Here, Density Functional Theory and Reactive Molecular Dynamics Simulations are combined to evaluate pristine, B-doped, N-doped, and B-N co-doped $γ$-graphyne and 6,6,12-graphyne (meta/ortho/para). $γ$-graphyne is a semiconductor, while 6,6,12-graphyne exhibits an anisotropic Dirac-like semi-metallic dispersion. B/N substitution reconstructs near-$E_F$ states via dopant $π$ hybridization, and B-N pairing stabilizes defects through donor-acceptor compensation, with the ortho substitutions being the most favorable. Hydrogen adsorption remains weak on pristine lattices but becomes locally optimized upon doping, with near thermo-neutral $ΔG_{\mathrm{ads}}$ 'hot spots' predominantly on $sp$-proximate carbon sites adjacent to the dopants. Reactive MD at 300 K further reveals an activity stability trade-off: B-N ortho in $γ$-graphyne sustains controlled hydrogen uptake without catastrophic bond scission, whereas B-N meta/para degrade, and 6,6,12-graphyne is generally more susceptible to over-hydrogenation. These results identify the B-N geometry as a key design variable for graphyne-based HER catalysts, which require both a favorable $ΔG_{\mathrm{ads}}$ and finite-temperature hydrogenation stability.

TopicsAtomistic Modeling of Sulfides and Minerals

Tagsadsorption density-functional-theory hydrogen-evolution molecular-dynamics

arXiv categoriescond-mat.mtrl-sci, cond-mat.other, physics.app-ph, physics.comp-ph

arXiv abstract pagePDF