Highly efficient photocatalytic reduction of nitrogen into ammonia by single Ru atom catalyst supported by BeO monolayer

Sifei Yun,Beixi Xu,Shujun Li,Wenio Zhu,Shuli Lei,?,Wenyue Guo,?,Ho Ren

aSchool of Materials Science and Engineering,China University of Petroleum(East China),Qingdao 266580,China

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

cInstitute of Mathematics,Free University of Berlin,Arnimallee 6,D-14195 Berlin,Germany

ABSTRACT Ammonia(NH3)is one of the most important building blocks of the chemical industry and a promising sustainable energy carrier.Conventional production of NH3 via the Haber-Bosch process requires high temperature and high pressure,which is energy demanding and suffers safety issues.Photocatalytic nitrogen reduction reaction(NRR)is a green and sustainable route for NH3 production,and has been expected to be an alternative for NH3 production under mild conditions.However,solar-driven N2 activated has appeared as the bottleneck for photocatalytic NRR.In this work,we propose that single Ru atom supported by BeO monolayer is a promising photocatalytic single atom catalyst(SAC)for efficient N2 activation with visible illumination.The high efficiency originates from the enhanced absorption in the visible range,as well as the back-donation mechanism when N2 were adsorbed on the SAC.Our results show that N2 can be efficiently activated by the Ru/BeO SAC and be reduced to NH3 with extremely low limiting potential of ?0.41 V.The NRR process also exhibits dominate selectivity respect to hydrogen evolution.

Keywords:Nitrogen reduction reaction Single atom catalyst Monolayer Ab initio Optical absorption spectra

Ammonia(NH3)is one of the most important building blocks for chemicals,such as fertilizer,explosives,and hydrogen carriers[1,2].Industrial ammonia is produced by Haber-Bosch process under high temperature and high pressure conditions(300–500 °C and 150–300 atm),consuming 1%–2% of the total global energy and releasing millions of tons of CO2[3].Developing an environmentally friendly method under mild conditions is a great challenge to replace Haber-Bosch production of NH3[4].Up to date,solar-driven nitrogen reduction reaction(NRR)appears to be the most promising way due to the superiority with environmentally benign,low-energy consuming and sustainability[5].The great challenge of photocatalytic NRR is effectively activating N2molecule with photogenerated electron transferring to the antibinding orbitals of N2molecule,which requires the photocatalyst chemically adsorbing N2for effective activation of the inert triple bonds[6].N2can be activated by back-donation process.The lonepair electrons of N2would transfer to the empty orbitals of activated sites,while thedorpelectrons of activate sites can be donated into the anti-binding states of N2[7,8].Tremendous experimental efforts have been taken to evaluate the performance of different catalysts,such as metal alloys[9],oxide[10],nitride[11],MOF[6,12],ZIF[13],perovskites[14].Thanks to the development of single atom catalyst(SAC),some SACs have been expected to have superior catalytic performance for NRR[7,15-18].Rational design of NRR SAC is drawing more and more attention recently.

Experimental design of NRR SAC is far from satisfactory and of high-cost.Density functional theory(DFT)calculations have been chosen as a convenient and low-cost approach for catalyst performance evaluation[19–22].More importantly,DFT provides deep insights into reaction mechanism at the atomic level[23–28].A promising NRR SAC has several prerequisites,such as thermodynamic and dynamic stabilities,the selectivity of NRR respect to hydrogen evolution reaction(HER),the ability to activate the inert N≡N triple bond[29],and more importantly,the limiting potential for contract the energy barrier of the rate-limiting step[19].As a general design SAC strategy,substrate materials are decorated by single heteroatom to form SAC,which acts as active site for N2adsorption and protonation with tunable electronic structures[30,31].

Fig.1.(a)Schematic view of single TM anchoring in BeO monolayer.Lime,red and blue violet balls represent Be,O and TM atoms,respectively.(b)Difference between binding energy of single TM atom anchoring on BeO monolayer with a Be vacancy and cohesive energy of isolate TM atoms aggregating metal clusters.(c)Variations of energy and temperature versus the AIMD simulation time for Ru/BeO,and the insets denote the top and side view of Ru/BeO after AIMD simulation last for 5 ps at T=500 K.

Bulk BeO is known to have a polar wurtzite structure with high thermal conductivity and hardness[32].It was claimed that few-layer BeO two-dimensional(2D)sheets with honeycomb structure can be synthesized through a wet chemistry approach[33].Recently,BeO monolayer is experimentally achievable on Ag(111)thin film by epitaxial growth[34].Theoretical study predicted that the electronic and magnetic properties of BeO monolayer can be tuned by intrinsic vacancy,transition metal substitutions and ribbon width[35–37].The optical properties of BeO monolayer has been reported for the heterostructure with germanene[38].BeO monolayer can also be used as two dimensional field-effect transistors on account of its insulting character[39].Compared with the bulk catalysts,2D catalysts has be reported to possess unique electronic structures and photocatalytic properties[40–42].Inspired by these,we thus wonder whether the BeO monolayer can be used for high efficient NRR catalyst.In this work,we rationally designed an excellent photocatalyst with thermodynamic stability,high activity and selectivity for NRR through single transition metal(TM)atom(TM=Sc,Ti,V,Cr,Mn,Fe,CO,Ni,Cu,Zn,Y,Zr,Nb,Mo,Ru,Rh,Pd,Ag and Cd)supported by BeO monolayer by using spin-polarized density functional theory.

All calculations were performed by using the Viennaab initioSimulation Package(VASP)[43].The projected augmented wave pseudopotential method with a cutoff energy of 500 eV for planewave basis was set to treat the electron-ion interactions[44].All atomic structures were fully relaxed geometry optimization with a converging tolerance of 0.02 eV/?A for residual force and 10?6eV for energy.The van der Waals(vdW)interactions were described by the semi-empirical DFT-D3 method[45].A vacuum space in thez-direction was set to 20 ?A for preventing the interaction between periodic images.All catalysts are modeled by single TM atoms substituting one Be in BeO supercells(4 × 4 × 1).The Brillouin Zone was sampled with the Gamma centered mesh with a 3 × 3 × 1kpoint grid for geometry relaxation and a 5 × 5 × 1k-point grid for electron structure analysis.And the band structures and optical absorption spectra were calculated by HSE06 method.The dynamical stability of SACs were confirmed byab initiomolecular dynamics(AIMD)simulations in canonical ensemble with the Nose-Hoover thermostat.

In order to investigate NRR process,six steps of coupled proton and electron transfer(CPET)are considered and the computational hydrogen electrode(CHE)model is employed in free energy profiles for reduction reaction,where the free energy of H++ e?is replaced by one-half of the chemical potential of the hydrogen molecule[46].The reaction Gibbs free energy of each elementary step is calculated by the following equation:ΔG=ΔE+ΔZPE-TΔS,whereΔEis the reaction energy,ΔZPEis the difference of zero point energy andΔSis the difference of reaction entropy between the product and reactant of each elementary step,respectively.Tis the absolute temperature.The entropy contribution of each gas phase species is obtained from NIST database[47].The optical absorption spectra is calculated based on the complex dielectric function[48].

As shown in Fig.1a,a series of TM atoms(TM=Sc,Ti,V,Cr,Mn,Fe,Co,Ni,Cu,Zn,Y,Zr,Nb,Mo,Ru,Rh,Pd,Ag and Cd)supported by BeO monolayer are screened as candidates for NRR.Binding energy(Eb)and cohesive energy(Ec)are calculated to determine the stability of SACs:

whereETM/BeO,EBeO?,ETMandrepresent the energy of a TM atom anchored on Be vacancy,the energy of BeO monolayer with a Be vacancy,the energy of an isolated single TM atom in vacuum,and the energy of a TM atom in the bulk crystal unit cell,respectively.

The anchoring single TM atom on BeO monolayer is thermodynamically favorable than diffusion and aggregation onceEb–Ec>0.As shown in Fig.1b,all SACs are stable.All calculated values are list in Table S1(Supporting information).Moreover,the thermal stability of Ru/BeO is further evaluated by AIMD simulation.Our results show that Ru/BeO is the best NRR catalyst among these 19 SACs.As shown in Fig.1c,during the AIMD equilibration,both the energy and temperature of Ru/BeO monolayer only oscillate within small ranges and there is no significant structure distortion at an elevated temperature of 500 K.The thermodynamical and dynamical stabilities of the Ru/BeO SAC meet the prerequisite as a promising NRR candidate.

N2adsorption is the first process for NRR and effective activation of N2is vital for further protonation.As shown in Fig.S1(Supporting information),the interaction between N2and BeO monolayer is physical adsorption with adsorption energy around?0.10 eV.Besides,the distance between N and active site(Be or O)varies in the range of 2.93～3.33 ?A,which is larger than covalent N-Be(1.73 ?A)and N–O(1.34 ?A)bond lengths based on covalent radii[49].The bond length of N≡N is 1.11 ?A,which is equal to free N2molecular in our calculation.Therefore,N2can not be effectively activated by pure BeO monolayer and that SACs consist of single TM atoms and BeO monolayer are considered for further study.

In order to make clear whether NRR or HER would dominate the reaction,the adsorption Gibbs free energies of H and N2with end-on and side-on pattern on SACs are calculated to evaluate the competitive adsorption between H2and N2.As shown in Fig.2a and Table S2(Supporting information),7 of 19 SACs are ruled out because these SACs neither adsorb H nor N2,namely,Fe/BeO,Ni/BeO,Cu/BeO,Zn/BeO,Pd/BeO,Ag/BeO,and Cd/BeO.H is the most favorable species on 4 SACs,that is Sc/BeO,V/BeO,Cr/BeO and Y/BeO,implying HER will dominate the reaction on these SACs.For the rest 8 SACs,Ti/BeO,Mn/BeO,Co/BeO and Rh/BeO prefer to adsorb N2with end-on pattern,and Zr/BeO,Nb/BeO,Mo/BeO and Ru/BeO prefer to side-on pattern,while both patterns would take place on Ti/BeO,Mn/BeO,Zr/BeO,Mo/BeO,and Ru/BeO.And the selectivity of HER/NRR(Table S3 in Supporting information)can be calculated by assuming that the occupation of SAC follows Boltzmann distribution:

Fig.2.(a)The histogram of adsorption Gibbs free energies of H and N2 with endon and side-on pattern on SACs.(b)The Gibbs free energies of the first and last protonation reaction.The horizontal dashed line denotes the selective criteria of 0.65 eV[50,51].Blue,orange and green bars represent H adsorption,N2 with endon pattern and N2 with side-on pattern,respectively.

whereΔGiis the adsorption Gibbs free energy of H and N2with end-on and side-on adsorption configurations,kis Boltzmann constant,andTis the reaction temperature of 298.15 K.So 8 SACs will be assessed in detail for protonation in the following parts.

The first and last protonation step of these 8 SACs,i.e.,hydrogenation of?N2into?N2H,?NH2into?NH3,is shown in Fig.2b.The first protonation step involves breaking the strong N≡N triple bond of the N2molecule,which can be uphill or downhill step determined by different active sites[21,50].For these 8 SACs,most of the first protonation steps are uphill,except for?N2(e)into?N2H on Mo/BeO.Previous studies show that 0.65 eV is the criteria for initial screening protonation step[50,51].The systems withΔGN2→N2H>0.65 eV are ruled out,including Mn/BeO,Co/BeO,and Rh/BeO,for the end-on adsorption configuration and Mn/BeO,Nb/BeO and Ru/BeO for the side-on adsorption configuration.

In the last protonation step of NRR,the adsorption intermediate of?NH2is much more stable than?NH3because thesp3hybrid orbitals of NH3molecule are fully filled,whilesp3hybrid orbital is half-filled in NH2[52].This step would consume extra energy to facilitate protonation of?NH2.We rule out the systems withΔGN2H→N3H>0.65 eV and then only Ru/BeO system for end-on adsorption is left for further study.

In order to assess the catalytic performance of Ru/BeO for the reduction of N2into NH3in detail,four typical mechanisms with all possible reaction intermediates are taken into consideration,i.e.,distal,alternating,mixed and enzymatic(Fig.S2 in Supporting information)[53].For end-on configuration,N2can be hydrogenated to NH3viathree different ways,distal,alternating and mixed mechanisms.For side-on configuration,N2is protonated by proton-electron pairs alternately and two NH3molecules are released in succession.

Free energy diagrams of Ru/BeO monolayer are illustrated in Fig.3.The free energy change for first hydrogenation of N2is significantly reduced compared to previous studies,with a free energy uphill of 0.41 eV and 0.69 eV for end-on and side-on adsorption,respectively.Previous studies show that the first protonation steps on Ru-based catalysts,such as Ru(0001)surface,Ru1-N3 and Ru1-N4 consume 1.08,0.73 and 0.77 eV energy[54,55].For the subsequent elementary steps through distal pathway in Fig.3a,the intermediate of?NNH could be hydrogenated to?NNH2with slight large free energy of 0.24 eV.The first NH3will be released into surroundings and?N remains on Ru/BeO accompany by the free energy downhill by 1.39 eV.The remaining?N would be protonated in the subsequent steps to form?NH,?NH2and?NH3,with free energy changes of 0.27,?0.47 and 0.18 eV,respectively.And the potential determining step of distal mechanism is the fourth protonation elementary step(?N + H++ e?=?NH).Finally,the second NH3will be released from Ru/BeO with 2.08 eV energy demanding,while this process does not involve hydrogenation.The produced NH3would be dissolved to NH4+in solution releasing energy,which can promote?NH3desorption[20,22,54].

As an alternating pathway,the hydrogenation of?NNH into?NNH2presents an even higher free energy,leading to a much negative limiting potential(?0.69 V)than the distal mechanism(?0.41 V)in Figs.3b and c.Limiting potential(UL)is the required potential to overcome the maximumΔGamong the six elementary reactions defined asUL=-ΔGmax/e[19].Occasionally,N2can be protonated by a mixed pathway that shifts from distal to alternating.It can be found that the potential determining step of mixed pathway in Fig.3c is the fourth protonation elementary step(?NHNH2+ H++ e?=?NH2NH2)and the onset potential is same as distal pathway.With respect to the enzymatic pathway,protonation also occurs alternatively to the two N of N2.As shown in Fig.3d,the potential determining step is the first hydrogenation step,with an uphill free energy of 0.69 eV,then the subsequent elementary steps are exothermic except for the last two steps.To the best of our knowledge,the limiting potential of Ru/BeO(?0.41 V)is comparable to other catalysis that have been reported,such as?0.20 V of B/g-C3N4,?0.31 V of B@g-CN,and ?0.33 V of Ru@g-C3N4[7,20,21].

The photoconversion efficiency of Ru/BeO is very important for highly efficient photocatalytic NRR.The band gap of pure BeO monolayer in this work is 5.4 eV in PBE functional and 6.84 eV for HSE06 functional(Fig.S3 in Supporting information),which is corresponding to previous studies with 4.0 eV to 7.4 eV[39,56–59]and band gap decrease to 1.22 eV and 2.89 eV after Ru doped in BeO monolayer with PBE and HSE06,respectively.As shown in Fig.S4(Supporting information),the DOS and LDOS indicate that the Ru states are introduced once Ru atom substitution,leading to decrease the band gap.And the introduced peaks near the Fermi level(0 eV)with the energy ranges(?0.50,0)and(0.85,1.25)eV,are visualized by LDOS plotted in the Fig.S4b.Pure BeO monolayer leads to less adsorption in the visible light owing to its relatively large band gap.Compared to pure BeO monolayer,the adsorption coefficient in the visible light and infrared(IR)light can be enhanced by Ru/BeO monolayer.As shown in Fig.4a,the main adsorption peak of BeO monolayer is located at～200 nm,meaning that pure BeO adsorb most ultraviolet(UV)light and very limited visible light.After single Ru atom decoration,the adsorption coeffi-cient will be enhanced,while adsorption between visible light and infrared light can be greatly enhanced by at least 1 order of magnitude compared to pure BeO monolayer.Therefore,Ru/BeO monolayer with superior optical absorption property maximize utilization of solar energy to drive nitrogen reduction.

As an excellent photocatalyst for NRR,the band edge positions of catalyst should match the N2/NH3potentials[20,60].Fig.4b displays that the conduction band minimum(CBM)lies above the N2/NH3potential and the valence band maximum(VBM)is lower than the N2/NH3potential.The photogenerated electrons on CBM can be transferred to N2for proton coupled electron process effectively,rather than recombine with the holes on VBM.Thus,the photoexcited electrons and holes would be effectively separated by Ru/BeO,leading to a higher photoconversion efficiency.These results indicate that Ru/BeO a promising photocatalyst for NRR with the merits of suitable visible light absorption and suitable band edge positions.

In summary,based on spin-polarized DFT calculations and CHE model,we rationally designed an efficient photocatalyst Ru/BeO with high thermodynamic stability,high activity and high NRR selectivity.Our results show that Ru/BeO exposes a relatively low limiting potential of ?0.41 V for NRR.In addition,as a promising photocatalyst,Ru/BeO shows visible light adsorption and suitable band edge positions.This work is expected to expand the application of 2D BeO monolayer and motivate more research efforts in photocatalysts.

Fig.3.Free energy diagrams for N2 reducing to NH3 on Ru/BeO via(a)distal,(b)alternating,(c)mixed and(d)enzymatic mechanisms under different applied potentials.Orange shadow represents the limiting rate step.

Fig.4.(a)Optical adsorption spectra of pure BeO and Ru/BeO,which are illustrated by the blue and red lines,respectively.(b)Band edge position of Ru/BeO.EAVS and ENHE denote energy levels relative to absolute vacuum scale(AVS)and normal hydrogen electrode(NHE),respectively.

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

We acknowledge the financial support from the National Natural Science Foundation of China(NSFC,Nos.21773309,21776315),the Fundamental Research Funds for the Central Universities(Nos.19CX05001A,20CX05010A),Hubei University of Arts and Science(Nos.2020kypytd002,2020kypytd003)and Xiangyang Science and Technology Research and Development(No.2020YL09).

Supplementary materials

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