Single Nb or W Atom-Embedded BP Monolayers as Highly Selective and Stable Electrocatalysts for Nitrogen Fixation with Low-Onset Potentials.
16 Feb 2021-ACS Applied Materials & Interfaces (American Chemical Society (ACS))-Vol. 13, Iss: 8, pp 10026-10036
TL;DR: In this article, a conductor-like screening model was proposed to screen out eligible NRR electrocatalysts from 29 atomic transition metals embedded on the defective boron phosphide (BP) monolayer with B-monovacancy (M/BP single-atom catalysts, SAC, M = Sc-Zn, Y-Cd, and Hf-Hg).
Abstract: Conversion of dinitrogen (N2) molecules into ammonia through electrochemical methods is a promising alternative to the traditional Haber-Bosch process. However, searching for an eligible electrocatalyst with high stability, low-onset potential, and superior selectivity is still one of the most challenging and attractive topics for the electrochemical N2 reduction reaction (NRR). Here, by means of first-principles calculations and the conductor-like screening model, four comprehensive criteria were proposed to screen out eligible NRR electrocatalysts from 29 atomic transition metals embedded on the defective boron phosphide (BP) monolayer with B-monovacancy (M/BP single-atom catalysts, SAC, M = Sc-Zn, Y-Cd, and Hf-Hg). Consequently, the Nb/BP and W/BP SACs are identified as the promising candidates, on which the N2 molecule can only be activated through the enzymatic pathway with the onset potentials of -0.25 and -0.19 V, and selectivities of 90.5 and 100%, respectively. It is worth noting that the W/BP SAC has the lowest overpotential among the 29 systems investigated. The electronic properties were also calculated in detail to analyze the activity origin. Importantly, the Nb/BP and W/BP SACs possess high thermal stabilities due to that their structures can be retained very well up to 1000 and 700 K, respectively. This work not only provides an efficient and reliable method to screen eligible NRR electrocatalysts but also paves a new way for advancing sustainable ammonia synthesis.
Citations
More filters
TL;DR: In this paper , the authors proposed a new capture-backdonation-recapture mechanism, which skillfully uses the electron capture between boron, the substrate, and other nonmetallic elements to solve the above problems.
Abstract: Facing the increasingly serious energy and environmental crisis, the development of heteronuclear metal-free double-atom catalysts is a potential strategy to realize efficient catalytic nitrogen reduction with low energy consumption and no pollution because it could combine the advantages of flexible active sites in double-atom catalysts while also being pollution-free and have high Faraday efficiency in metal-free catalysts simultaneously. However, according to the existing mechanism, the finite orbits of other nonmetallic atoms, except the boron atom, reduce the possibility of metal-free catalysis and hinder the development of heteronuclear metal-free double-atom catalysts. Herein, we propose a new "capture-backdonation-recapture" mechanism, which skillfully uses the electron capture-backdonation-recapture between boron, the substrate, and other nonmetallic elements to solve the above problems. Based on this mechanism, by means of the first-principle calculations, the material structure, adsorption energy, catalytic activity, and selectivity of 36 catalysts are systematically investigated to evaluate their catalytic performance. B-Si@BP1 and B-Si@BP3 are selected for their good catalytic performance and low limiting potentials of -0.14 and -0.10 V, respectively. Meanwhile, the "capture-backdonation-recapture" mechanism is also verified by analyzing the results of adsorption energy and electron transfer. Our work broadens the ideas and lays the theoretical foundation for the development of heteronuclear metal-free double-atom catalysts in the future.
61 citations
TL;DR: In this paper, a novel strategy to promote NRR through boron-transition-metal (TM) hybrid double-atom catalysts (HDACs) has been proposed.
Abstract: Electrocatalytic nitrogen reduction reaction (NRR) is a promising method for sustainable production of NH3, which provides an alternative to the traditional Haber-Bosch process. However, the poor Faraday efficiency caused by N≡N triple bond activation and competitive hydrogen evolution reaction (HER) have seriously hindered the application of NRR. In this work, a novel strategy to promote NRR through boron-transition-metal (TM) hybrid double-atom catalysts (HDACs) has been proposed. The excellent catalytic activity of HDACs is attributed to a significant difference of valence electron distribution between boron and TMs, which could better activate N≡N bonds and promote the conversion of NH2 to NH3 compared with boron or metal single-atom catalysts and traditional double-atom catalysts (DACs). Hence, by means of DFT computations, the stability, activity, and selectivity of 29 HDACs are systematically investigated to evaluate their catalytic performance. B-Ti@g-CN and B-Ta@g-CN are screened as excellent nitrogen-fixing catalysts with particularly low limiting potentials of 0.13 and 0.11 V for NRR and rather high potentials of 0.54 and 0.82 V for HER, respectively. This work provides a new idea for the rational design of efficient nitrogen-fixing catalysts and could also be widely used in other catalytic reactions.
52 citations
TL;DR: In this paper, a series of transition-metal-doped hexagonal boron phosphide (h-BP) monolayers were comprehensively evaluated using density functional theory, and it was shown that the Vdoped h-BP monolayer stands near the top of the volcano plot with the limiting potential of -0.22 V versus a reversible hydrogen electrode, exhibiting the lowest overpotential among the investigated systems.
Abstract: Electrocatalytic conversion of nitrate (NO3-) into ammonia can not only eliminate harmful pollutant but also provide a green method for a low-temperature ammonia synthesis. The electrochemical NO3- reduction reactions (NO3RRs) of a series of transition-metal-doped hexagonal boron phosphide (h-BP) monolayers were comprehensively evaluated using density functional theory. The V-doped h-BP monolayer was found to stand near the top of the volcano plot with the limiting potential of -0.22 V versus a reversible hydrogen electrode, exhibiting the lowest overpotential among the investigated systems in this work. Besides, the competing hydrogen evolution reaction is significantly suppressed due to the weak adsorption of the H atom. Importantly, the structure of the V-doped h-BP monolayer can be retained very well until 900 K, illustrating the initial indication of high thermal stability and great promise for synthesis. This study not only offers an eligible NO3RR electrocatalyst but also provides an atomic understanding of the behind mechanisms of the NO3RR process.
44 citations
TL;DR: In this article, the application potentials of the tetragonal transition-metal phosphide monolayers (M2P, M = Co, Ni, Ru and Pd) as NO3RR electrocatalysts for NH3 synthesis have been comprehensively and originally studied using density functional theory.
26 citations
TL;DR: In this paper , the potentials of 29 transition-metal atoms anchored on the nitrogen-doped BP monolayer (MN3/BP) as efficient NORR catalysts were systematically examined using first-principles calculations.
21 citations
References
More filters
TL;DR: A simple derivation of a simple GGA is presented, in which all parameters (other than those in LSD) are fundamental constants, and only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked.
Abstract: Generalized gradient approximations (GGA’s) for the exchange-correlation energy improve upon the local spin density (LSD) description of atoms, molecules, and solids. We present a simple derivation of a simple GGA, in which all parameters (other than those in LSD) are fundamental constants. Only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked. Improvements over PW91 include an accurate description of the linear response of the uniform electron gas, correct behavior under uniform scaling, and a smoother potential. [S0031-9007(96)01479-2] PACS numbers: 71.15.Mb, 71.45.Gm Kohn-Sham density functional theory [1,2] is widely used for self-consistent-field electronic structure calculations of the ground-state properties of atoms, molecules, and solids. In this theory, only the exchange-correlation energy EXC › EX 1 EC as a functional of the electron spin densities n"srd and n#srd must be approximated. The most popular functionals have a form appropriate for slowly varying densities: the local spin density (LSD) approximation Z d 3 rn e unif
146,533 citations
TL;DR: A new density functional of the generalized gradient approximation (GGA) type for general chemistry applications termed B97‐D is proposed, based on Becke's power‐series ansatz from 1997, and is explicitly parameterized by including damped atom‐pairwise dispersion corrections of the form C6 · R−6.
Abstract: A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.
23,058 citations
TL;DR: The dynamical steady-state probability density is found in an extended phase space with variables x, p/sub x/, V, epsilon-dot, and zeta, where the x are reduced distances and the two variables epsilus-dot andZeta act as thermodynamic friction coefficients.
Abstract: Nos\'e has modified Newtonian dynamics so as to reproduce both the canonical and the isothermal-isobaric probability densities in the phase space of an N-body system. He did this by scaling time (with s) and distance (with ${V}^{1/D}$ in D dimensions) through Lagrangian equations of motion. The dynamical equations describe the evolution of these two scaling variables and their two conjugate momenta ${p}_{s}$ and ${p}_{v}$. Here we develop a slightly different set of equations, free of time scaling. We find the dynamical steady-state probability density in an extended phase space with variables x, ${p}_{x}$, V, \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{}, and \ensuremath{\zeta}, where the x are reduced distances and the two variables \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{} and \ensuremath{\zeta} act as thermodynamic friction coefficients. We find that these friction coefficients have Gaussian distributions. From the distributions the extent of small-system non-Newtonian behavior can be estimated. We illustrate the dynamical equations by considering their application to the simplest possible case, a one-dimensional classical harmonic oscillator.
17,939 citations
TL;DR: In this article, the authors compared the canonical distribution in both momentum and coordinate space with three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 77, 63 (1983); and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983).
Abstract: Three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 48, 1818 (1982)], and Evans and Morriss [Chem. Phys. 77, 63 (1983)]; and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983)] are examined analytically via calculating the equilibrium distribution functions and comparing them with that of the canonical ensemble. Except for effects due to momentum and angular momentum conservation, method (1) yields the rigorous canonical distribution in both momentum and coordinate space. Method (2) can be made rigorous in coordinate space, and can be derived from method (1) by imposing a specific constraint. Method (3) is not rigorous and gives a deviation of order N−1/2 from the canonical distribution (N the number of particles). The results for the constant temperature–constant pressure ensemble are similar to the canonical ensemble case.
13,921 citations
TL;DR: In this paper, the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids is described and the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets.
Abstract: Recent extensions of the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids are described. Furthermore the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets. The method is applied to calculations of the enthalpy of formation of a set of molecules and solids. We find that the present numerical localized basis sets yield improved results as compared to previous results for the same functionals. Enthalpies for the formation of H, N, O, F, Cl, and C, Si, S atoms from the thermodynamic reference states are calculated at the same level of theory. It is found that the performance in predicting molecular enthalpies of formation is markedly improved for the Perdew–Burke–Ernzerhof [Phys. Rev. Lett. 77, 3865 (1996)] functional.
8,496 citations