Computational prediction of Mo2@g-C6N6 monolayer as an efficient electrocatalyst for N2 reduction

Jijun Wng,Mengyo Shi,Guolin Yi,Lu Wng,Shuli Lei,Ke Xu,Shujun Li,c,Jinshui Mu,*

a Tianjin Key Laboratory of Structure and Performance for Functional Molecules,College of Chemistry,Tianjin Normal University,Tianjin 300387,China

b Hubei Key Laboratory of Low Dimensional Optoelectronic Materials and Devices,Hubei University of Arts and Science,Xiangyang 441053,China

c Institute of Mathematics,Free University of Berlin,Berlin D-14195,Germany

Keywords:g-C6N6 monolayer Double-atom catalysts Nitrogen reduction reaction Hydrogen evolution reaction Limiting potential Density functional theory

ABSTRACT Electrocatalytic nitrogen reduction reaction (NRR) is an environmentally friendly method for sustainable ammonia synthesis under ambient conditions.Searching for efficient NRR electrocatalysts with high activity and selectivity is currently urgent but remains great challenge.Herein,we systematically investigate the NRR catalytic activities of single and double transition metal atoms (TM=Fe,Co,Ni and Mo) anchored on g-C6N6 monolayers by performing first-principles calculation.Based on the stability,activity,and selectivity analysis,Mo2@g-C6N6 monolayer is screened out as the most promising candidate for NRR.Further exploration of the reaction mechanism demonstrates that the Mo dimer anchored on g-C6N6 can sufficiently activate and efficiently reduce the inert nitrogen molecule to ammonia through a preferred distal pathway with a particularly low limiting potential of-0.06 V.In addition,we find that Mo2@g-C6N6 has excellent NRR selectivity over the competing hydrogen evolution reaction,with the Faradaic efficiency being 100%.Our work not only predicts a kind of ideal NRR electrocatalyst but also encouraging more experimental and theoretical efforts to develop novel double-atom catalysts (DACs) for NRR.

N2fixation for NH3synthesis is of great significance since NH3is one of the main feedstocks for chemical and agricultural industries and also a promising carbon-free energy carrier with high energy density for the future economy [1-5].Compared with the conventional energy-and resource-intensive Haber-Bosch process,electrocatalytic N2reduction reaction (NRR),which hydrogenates N2as well as protons and electrons at ambient conditions,is regarded as a more cost-effective,sustainable,and environmentally benign method for NH3synthesis [6-8].However,as a nonpolar homonuclear diatomic molecule,N2is fairly inert due to the extremely thermodynamically stable covalent triple bond(941 kJ/mol) and thus is difficult to break in the NRR process,resulting in low N2-to-NH3conversion efficiency [9-12].Consequently,the rational design of efficient NRR electrocatalysts for NH3synthesis is highly desirable but remains great challenge.

Today,single-atom catalysts (SACs) comprising isolated atoms dispersed on substrates are of particular interest in the field of electrocatalysis due to their exceptional catalytic performances and maximized atomic utilization [13-15].In particular,many SACs of transition-metal (TM) atoms anchored with two-dimensional(2D) materials,including graphene [16,17],h-BN [18,19],borophene[20,21],phosphorene [22-24],graphyne [25,26],porous carbonnitride nanosheets [27-30],nitrogen-doped carbon [31,32]and so on [33-35],have been theoretically predicted as striking NRR electrocatalysts with ultralow limiting potential.However,as multiple reaction intermediates are typically involved in the NRR process,experimental studies reveal that the single active site is still rather hard to break the scaling relations and enhance both the yield rate of NH3and the Faradaic efficiency (FE) simultaneously.To address this issue,much effort has been devoted to designing doubleatom catalysts (DACs),which possess more flexible active sites to tune the adsorption behavior of NRR intermediates [36-40].For instance,Liet al.[41]synthesized atomically dispersed Fe-Mo dimers anchored on N-doped graphene (FeMo@NG) and observed that it exhibited better NRR catalytic activity than its SAC counterparts(Fe@NG and Mo@NG) due to the combination of geometric,ligand and synergistic effects.Hanet al.[42]reported that the asfabricated Pd-Cu diatom catalyst had higher FE and NH3yield rate than the individual single-atom (Pd or Cu) catalyst.In addition,theoretical screening and designing based on first-principles calculations found that Mo-,Fe-,Co-,Mn-and Ru-based DACs can be served as excellent catalysts for electrochemical N2reduction[43-55].

Fig.1.(a) The atomic structure of pristine g-C6N6 monolayer.The rhombus denotes the unit cell,and the brown and light-blue spheres stand for C and N atoms,respectively.(b) The optimized most stable configuration and (c) calculated binding energies of TM@g-C6N6 and TM2@g-C6N6.

Recently,a novel porous graphitic carbon nitride (g-C6N6) is successfully synthesizedviathe solvothermal reaction of C3N3Cl3with Na [56].Interestingly,g-C6N6possesses a high specific surface area,excellent thermal and kinetic stability,and outstanding electric conductivity,making it a remarkable substrate to anchor TM atoms for electrocatalytic reactions [57-59].Inspired by these exciting progressive achievements,in the present work,a series of single and double TM atoms (TM=Fe,Co,Ni and Mo) anchored on g-C6N6are designed and explored as NRR electrocatalysts.Among them,Mo2@g-C6N6is screened out as the most promising candidate for electrochemical NRR due to its high stability,activity,and selectivity.Our work indicates that DACs have great chance to replace SACs as the emerging star of atomic catalysts.

All the spin-polarized first-principles calculations are performed using the density functional theory (DFT),as implemented in Viennaab initioSimulation Package (VASP) [60].The interactions between electrons and ion-cores are described by the projector augmented wave (PAW) method [61]and the exchange-correlation interactions are treated by the generalized gradient approximation (GGA) in the form proposed by Perdew,Burke,and Ernzerhof(PBE) [62].Additionally,Grimme’s semiempirical DFT-D3 method is employed to address the weak van der Waals interactions between the adsorbates and substrate [63].The cut-off energy is set to 500 eV for the plane-wave basis in all of our calculations.Γ-centered k-point meshes of 2×2×1 and 4×4×1 based on Monkhorst-Pack scheme [64]are employed for the geometric and electronic structures calculations,respectively.A vacuum space of 20.0 ˚A along thez-axis is adopted to ensure no appreciable interaction between the image layers under periodic boundary condition.The self-consistent convergence criterion for the total energy and Hellmann-Feynman force are smaller than 10-5eV and 0.02 eV/˚A,respectively.More computational details for the NRR are given in the Supporting information.

The optimized geometrical structure of pristine g-C6N6monolayer is depicted in Fig.1a,which contains six carbon atoms and six nitrogen atoms in the hexagonal unit cell with space group P6/MMM and point group D6H-1.The calculated lattice constants are found to bea=b=7.12 ˚A with the C-C and C-N bond distances of 1.51 and 1.34 ˚A,respectively,in good agreement with previous reports [65-67].In this study,a 2×2×1 g-C6N6supercell is employed as an ideal substrate for anchoring the selected four types of common TM atoms (TM=Fe,Co,Ni and Mo).To find the most favorable anchoring sites,we examine two possible configurations of the single and double TM atoms on the g-C6N6monolayer (labeled as TM@g-C6N6and TM2@g-C6N6,respectively),in which TM atoms bond with two N atoms or three N atoms of g-C6N6substrate.As shown in Fig.1b,we find that the single TM atom bonds strongly with adjacent two N atoms at the cavity edge after full optimization except for the Mo atom,while all the double TM atoms are located within the six-remember nitrogen hole and bond with three N atoms.The distances of TM-TM bonds are 2.13 ˚A for Fe,2.18 ˚A for Co,2.34 ˚A for Ni,and 1.78 ˚A for Mo.The planar structure of the g-C6N6monolayer can be maintained well before and after anchoring these TM atoms.Compared to pristine g-C6N6with a band gap of 1.83 eV (Fig.S1 in Supporting information),these TMn@g-C6N6monolayers exhibit metallic features mainly originated from the hybridized states between TM-d and the adjacent N-p orbitals (Fig.S2 in Supporting information),which will be propitious to the electrocatalytic reactions.

To estimate the structural stability of TMn@g-C6N6,we calculate the binding energies (Eb),as shown in Fig.1c.Overall,all double TM anchored on g-C6N6show larger binding energies compared to the corresponding single TM anchored cases,suggesting that the formed metal-metal bonding can offer positive effects on stabilization of these TM dimers.Moreover,theEbof TMn@g-C6N6are more negative than their respective cohesive energies of bulk metal,signifying strong interactions between the TM atoms and g-C6N6monolayer,and superior thermodynamic stability.In addition,we further check the temperature-dependent stability of TMn@g-C6N6byab initiomolecular dynamics (AIMD) simulations as depicted in Fig.S3 (Supporting information).We observe that these TMn@g-C6N6systems display narrow energy fluctuations and slight structural changes during the simulated time,verifying the excellent thermal stability at room temperature.Thus,it is believed that the as-designed TMn@g-C6N6could serve as potential electrocatalyst with high stability.

Following the screening requirements for NRR catalyst [68,69],we turn to examine the Gibbs free energy changes (△G) of the N2adsorption and first hydrogenation process of NRR to produce*N2H on TMn@g-C6N6for identifying the most qualified candidate.The computed △Gfor the two processes are presented in Figs.2a and b,respectively.Our results show that all the TM@g-C6N6systems are incapable of capturing N2molecule through the side-on configuration because of their positive △G*N2,whereas N2molecule can effectively chemisorb on TM@g-C6N6by means of the end-on configuration.However,the first protonation of adsorbed N2(*N2+H++e-→*N2H) require much high energy inputs,such as 1.34 eV for Fe@g-C6N6,1.33 eV for Co@g-C6N6,1.25 eV for Ni@g-C6N6,and 0.51 eV for Mo@g-C6N6,demonstrating that they are inefficient electrocatalysts for N2reduction.In this regard,the TM@g-C6N6systems are eliminated as eligible NRR catalysts.

With respect to TM@g-C6N6,we find that except for Ni2@g-C6N6,N2molecule can be well captured by TM2@g-C6N6regardless of end-on or side-on configuration,as indicated by the calculated△G*N2ranging from-0.46 eV to-0.81 eV.Further,the △Gvalues of*N2H formation on the Fe2@g-C6N6and Co2@g-C6N6are predicted to be at least 0.48 and 0.72 eV,respectively,marking them as inefficient electrocatalysts for N2reduction.Notably,the first hydrogenation process of NRR to produce*N2H on Mo2@g-C6N6demands only 0.11 eV energy input for the side-on configuration,while it even becomes a spontaneous reaction for the end-on configuration.Based on these screening results,only the Mo2@g-C6N6is of interest in virtue of satisfying all the requirements,and thus will be examined in the following discussion.

Fig.2.Calculated Gibbs free energies (△G) for (a) N2 adsorption and (b) *N2H formation on TM@g-C6N6 and TM2@g-C6N6.

Fig.3.(a) End-on and (b) side-on adsorption configurations for N2 on Mo2@g-C6N6 monolayer.(c,d) The corresponding charge density difference,where the isosurface value is set to be 0.005 e/˚A3 and the yellow and cyan show the positive and negative charges,respectively.The computed projected density of states of Mo2@g-C6N6 with the N2 adsorption through (e) end-on and (f) side-on configurations.

Fig.4.(a) Schematic depiction of three possible reaction mechanisms for N2 electroreduction to NH3 on the Mo2@g-C6N6 monolayer.Optimized geometries of all reaction intermediates along the (b) distal,(c) alternating,and (d) enzymatic pathways.

Figs.3a and b depict the optimized structures of N2adsorbed Mo2@g-C6N6viaend-on and side-on configurations,respectively.Compared with the isolated N2molecule (1.12 ˚A),the N-N bond length is significantly stretched by 0.04 ˚A for end-on adsorption and by 0.09 ˚A for side-on adsorption,signifying the effective activation of the inert triple bond of N2.Bader charge analysis shows that there are about 0.49 and 0.70 electrons transferred from Mo2@g-C6N6to the adsorbed N2in end-on and side-on patterns,respectively,which attributes to the different electronegativities of N (3.04) and Mo (1.86).Moreover,as presented in Figs.3c and d,we can observe the obvious charge accumulation and depletion around the anchored Mo dimer and N2for both adsorption configurations.Interestingly,such a charge distribution ideally accords with the "acceptance-donation" process as proposed by Ling and co-workers [70],that is,the Mo dimer donates electrons into the antibonding orbitals of N2and simultaneously accepts lone-pair electrons from it.From the computed projected density of states given in Figs.3e and f,we find that there are significant electronic coupling between the N-2p orbital of the adsorbed N2molecule and Mo-4d orbital,thus give rise to the anti-bonding states located above the Fermi level and bonding states lied below the Fermi level,further supporting the "acceptance-donation" picture.Moreover,after N2adsorption with side-on configuration,the magnetic moments of the Mo dimer and N2molecule are nearly unchanged and remain 0 μB.However,for N2adsorption with end-on configuration,we observe that their corresponding magnetic moments increase to 0.97 μB and 0.24 μB,respectively.As shown in Fig.S4 (Supporting information),the increased magnetic moments attribute to the asymmetrical distributions of Mo-4d and N-2p orbitals,caused by the couplings of Mo 4dz-N 2py,Mo 4dxz-N 2pzand Mo 4dz2-N 2pz.These above analyses demonstrate the successful activation of N2on Mo2@g-C6N6,which would further facilitate the subsequent reaction of converting N2to NH3molecules.

As illustrated in Fig.4a,three possible reaction mechanisms,named distal,alternating,and enzymatic,are considered for the transformation from N2to NH3.Both distal and alternating mechanisms start from the adsorbed N2with end-on configuration,while the enzymatic mechanism is initiated with the side-on adsorbed N2.Whereafter,each mechanism involves six consecutive protonation plus reduction steps,accompanied by the production of two NH3molecules.The optimized geometries of all reaction intermediates along these three pathways on Mo2@g-C6N6monolayer are given in Figs.4b-d,and the predicted Gibbs free energy profiles are illustrated in Figs.5a-c.

As NRR proceeds along the distal or alternating pathway,the first protonation plus reduction step is similar,and found to be exothermic by 0.02 eV with the further elongation of N-N bond(from 1.16 ˚A to 1.25 ˚A).It is noteworthy that this first hydrogenation is normally non-spontaneous for most studied electrocatalysts since the free energy change is thermodynamically uphill.For the following elementary steps through distal pathway,the*N-NH can be spontaneously hydrogenated to*NH2by reacting with the proton-electron pair four times.Specifically,as shown in Fig.5a,the Gibbs reaction free energies of these four hydrogenation steps are-0.10,-0.22,-0.80 and-0.81 eV,respectively.Finally,the process of*NH2+H++e-→*NH3is slightly uphill by 0.06 eV in the Gibbs free energy profile.Overall,the potentiallimiting step (PDS) of the distal mechanism is the sixth protonation plus reduction step owing to its maximum △Gof 0.06 eV.For the case of the alternating mechanism,as given in Fig.5b,the Gibbs free energy changes along the subsequent protonation plus reduction steps are calculated respectively to be-0.02,0.07,-0.34,0.22,-1.88 and 0.06 eV.Thus,the hydrogenation of*NH-NH2into*NH2-NH2is the PDS with a maximal energy demand of 0.22 eV.Regarding of the enzymatic mechanism,the first protonation plus reduction step is slightly uphill energetically,with a △Gvalue of 0.11 eV as demonstrated in Fig.5c.After that,the subsequent elementary reaction steps are all exothermic,except for the formation of*NH-*NH2(*NH-*NH+H++e-→*NH-*NH2) and second NH3(*NH2+H++e-→*NH3).Among the hydrogenation steps of the enzymatic mechanism,the step converting*N-*NH to*NH-*NH2has the largest free energy change (0.28 eV),which can be identified as the PDS.Of note,the desorption of the second NH3from Mo2@g-C6N6monolayer is not a problematic obstacle.As demonstrated in previous studies [55,66],the*NH3species can be released in the form of NH4+under acidic electrochemical conditions.Moreover,the free energy change for the release of the final product NH3is only 0.82 eV,which can be easily overcame in consideration of the released energy (about 2.60 eV,Fig.5) in the proceeding hydrogenation steps.These ensure the rapid removal of the produced NH3under ambient conditions for recovering the catalysts.

Fig.5.Gibbs free energy profiles of N2 reduction on Mo2@g-C6N6 monolayer via the (a) distal,(b) alternating and (c) enzymatic mechanisms at zero and applied potentials.(d) Charge variation of three moieties for the NRR along the preferred distal pathway.

In electrocatalysis,the limiting potential (UL),defined as the lowest applied electrode potential to overcome the maximum positive △Gon the considered pathway,is the most commonly used descriptor for chemical reactivity,which can be calculated by UL=-△Gmax/e [71,72].A less negativeULindicates a faster electrochemical process at a given potential.As discussed above,the limiting potential with respect to the standard hydrogen electrode are calculated to be-0.06 V for distal,-0.22 V for alternating,and-0.28 V for enzymatic.Consequently,we predict that the NRR on Mo2@g-C6N6will prefer to proceedviathe distal mechanism in virtue of its lowest negative value.Remarkably,such a limiting potential is substantially less negative than the equilibrium potential of overall NRR (-0.16 V) as well as other recently reported DACs for NH3synthesis,implying that the applied potential even as low as the equilibrium potential can adequately drive the NRR.Therefore,Mo2@g-C6N6could be expected as a promising candidate electrocatalyst for N2reduction to NH3.

It should be pointed out that conventional DFT method often fails to describe the systems with strongly correlated d-electrons.In order to further evaluate the reliability of the above results,we employ the DFT+U method to re-examine Gibbs free energy profile of N2reduction on Mo2@g-C6N6monolayer through the preferred distal pathway.As shown in Fig.S5 (Supporting information),we can see that the Hubbard U value has little influence on the free energy diagram,and both of the limiting potential and PDS with the Hubbard U correction are in line with the PBE results,suggesting that the PBE calculated results are acceptable in this work.

To gain insight into the superior catalytic activity of Mo2@g-C6N6,we further employ Bader charge analysis to investigate the charge variation of various intermediates along the favorable distal pathway.Herein,we divide each intermediate into three groups,including the g-C6N6,anchored Mo dimer,and adsorbed NxHyspecies.As depicted in Fig.5d,the N2molecule after chemisorption on the Mo dimer accumulates about 0.49 electrons,which is considered to be propitious for the hydrogeneration of*N2into*N2H.In the following NRR steps,it is observed that there are obvious charge fluctuations in both of the g-C6N6and NxHygroups,while the charge value of Mo dimer remains almost unchanged.That is to say,the charge variation of the NxHyspecies mainly associates with that of the g-C6N6.The Mo dimer serves as not only the active site but also electron transmitter between the g-C6N6and NxHyspecies during the entire NRR process.

As one of the prerequisites of an excellent NRR electrocatalyst with high FE,it is crucial to suppress the hydrogen evolution reaction (HER),which is the major competing reaction during the NRR.Therefore,we compute the adsorption free energies of*H (△G*H)at the Mo dimer site of Mo2@g-C6N6monolayer.The △G*His found to be 0.46 eV,which is more positive than the △G*N2(-0.71 eV),suggesting that the active site would be primarily covered by*N2.Moreover,the free energy barrier for HER (0.46 eV) is considerable larger than the PDS barrier for NRR (0.06 eV).In this respect,the FE of Mo2@g-C6N6is estimated to be approximately 100% at room temperature according to the Boltzmann distribution,indicating the high selectivity for NRR.

In conclusion,by performing comprehensive DFT computations,we explore the potential of single and double transition metal atoms (TM=Fe,Co,Ni and Mo) anchored g-C6N6monolayers for electrochemical N2reduction to NH3.Among these systems,Mo2@g-C6N6monolayer is identified as the most promising catalyst candidate toward NRR with the help of prescreening criteria.Our results demonstrate that N2molecule is effectively captured and activated by the Mo dimer anchored on g-C6N6through the"acceptance-donation" process.The subsequent N2reduction reaction on Mo2@g-C6N6monolayer proceed dominantlyviathe distal mechanism with an extremely low limiting potential of-0.06 V.Importantly,the competitive HER can be well inhibited on the Mo2@ g-C6N6monolayer due to a much larger △G*H,ensuring substantial selectivity (～100%) toward NH3synthesis.Hence,the Mo2@g-C6N6monolayer proposed in this work has great potential applications for NRR,highlighting the importance of DACs for NRR.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was supported by the Science &Technology Development Fund of Tianjin Education Commission for Higher Education (No.2020KJ008),by the Natural Science Foundation of Tianjin(No.18JCQNJC76000),and by the College Students’Innovation and Entrepreneurship Training Program of Tianjin (No.202110065112),Science and Technology Research Project of Hubei Provincial Department of Education (No.D20212603),Hubei University of Arts and Science (Nos.2020kypytd002,XK2021024) and China Scholarship Council.

Supplementary materials

Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2021.12.040.