• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    非金屬電催化劑環(huán)境條件下氮還原反應(yīng)的研究進(jìn)展

    2021-07-10 10:10:28許桐馬奔原梁杰岳魯超劉倩李廷帥趙海濤羅永嵐盧思宇孫旭平
    物理化學(xué)學(xué)報 2021年7期
    關(guān)鍵詞:劉倩

    許桐,馬奔原,梁杰,岳魯超,劉倩,李廷帥,趙海濤,羅永嵐,盧思宇,孫旭平,*

    1 電子科技大學(xué),基礎(chǔ)與前沿研究院,成都 610054

    2 鄭州大學(xué),化學(xué)學(xué)院與綠色催化研究中心,鄭州 450001

    1 Introduction

    NH3plays an important role in modern society as essential building block in manufacturing fertilizer, aqueous ammonia,plastic, explosive and dyeetc1,2. It also receives attention as green alternative fuel because of its carbon-free nature, large hydrogen capacity (17.6% (w)), high energy density (4.3 kW·h·kg?1) and easy transportation (?33 °C, liquid)3. Currently,industrial NH3production mainly relies on Haber-Bosch process at high pressure and temperature using N2and H2as feeding gas4.However, this process accounts for ~2% global energy consumption and ~1% CO2emission, which significantly stimulates the exploration of alternative approaches.Encouragingly, electrochemical N2reduction reaction (NRR)offers a promising route towards sustainable and environmentally-benign ammonia synthesis at ambient conditions, and its applied potential is compatible with intermittent solar, wind and other renewable energies5–8.However, it needs efficient electrocatalysts to drive N2-to-NH3conversion because of the extremely inert N≡N bond (945 kJ·mol?1bond energy).

    To date, great efforts have been put into exploring highperformance catalysts towards high efficiency and selectivity.Generally, noble-metal catalysts perform efficiently for the NRR, but the scarcity and high cost limit their large-scale applications9–12. Therefore, much attention has focused on earthabundant transition-metal (TM) catalysts which can activate N2by an acceptance-donation process, using empty or unoccupied orbitals to accept the lone-pair electrons of N2meanwhile donating the abundantd-orbital electrons into the antibonding orbitals of N213–35. However, these catalysts may release metal ions, leading to environmental pollution36,37. Nevertheless, most of these TM electrocatalysts show moderate proton adsorption free energy and thed-orbital electrons of TM may also favor the formation of TM―H bonds, facilitating the hydrogen evolution reaction (HER) during the electrocatalytic reaction. Recent years, it is a hot pot to explore metal-free catalysts (MFCs)38,39.MFCs mainly include carbon-based catalysts (CBCs) and some boron-based and phosphorus-based catalysts. Usually, CBCs have porous structure and high surface area, which are favorable to expose more active sites and provide rich accessible channels for mass/electron transfer. Furthermore, heteroatoms doping(such as N, S, O, B, P and F), defect engineering, edge engineering and heterojunction can induce more intrinsic active sites or construct favorable catalytic environment by redistributing charge density and tailoring electronic structure38–41.More importantly, the graphitic carbon of CBCs can greatly enhance the electron conductivity which is necessary to supply enough reactive electrons to active site and antibonding orbital of N241,42. The unsaturated carbon atoms (spandsp2) may favor heteroatom doping by accepting foreign electrons in unoccupied orbitals43,44. Finally, CBCs are generally made up of earthabundant elements (C, B, N, O, S) with good durability in acidic environment, taking the advantages of low cost and stability38,45Moreover, the Lewis acid sites of most metal-free compounds could also accept the lone-pair electron of N2and adsorb N2molecules by forming nonmetal―N bond, further widening their potential in electrocatalytic NRR. However, compared with metal-based catalyst, the occupied orbitals of metal-free catalysts can only form covalent bond or conjugatedπbond,hindering the electron donation from electrocatalyst to N2and molecular activation.

    This review summarizes recent advances in designing and developing MFCs towards NRR under ambient conditions in aqueous media, including compounds, heteroatom-doping, and organic polymers. We also discuss the strategies of boosting NRR performances and outlook the development perspectives of MFCs.

    2 Fundamentals of electrocatalytic N2 reduction

    The inertness of the N2molecule partly originates from the high bond energy (941 kJ·mol?1) of the triple bond. On account of the negative electron affinity, low proton affinity and nonpolarity render N2be Lewis base. The large energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of N2is not benefit for electron transfer processes. There are two reactive mechanisms responsible for the N2-to-NH3conversion, including dissociative and associative pathways41,46. The dissociative pathway mainly concerns the Haber-Bosch process, by high-energy cleaving N≡N to obtain isolated N atoms which then react with H atoms to form NH3molecules. The associative pathway involves the biological and electrochemical N2fixation, for which hydrogenation process occurs but keeping the two N atoms of N2bonding with each other before the first NH3is generated.According to N2adsorption pattern and hydrogenation order, the associative mechanism can be further classified into distal,alternating and enzymatic pathways. In distal pathway,hydrogenation firstly occurs on the N atom away from catalyst surface to generate one NH3, and then another N atom of N2bonded on catalyst surface. As to alternating pathway,hydrogenation occurs on the two N atoms alternatively. The enzymatic pathway also has an alternatively hydrogenation process but simultaneously keeping the two N atoms bonded to catalyst surface before the first NH3desorption. For the associative mechanism, the catalyst is vitally important in favoring N2adsorption, activation, hydrogenation and NH3desorption.

    3 Carbon-based catalysts

    3.1 Heteroatom doping

    3.1.1 N atom doping

    Fig. 1 (a) Schematic of NPC preparation. (b) Contents of pyridinic,pyrrolic and graphitic N in NPCs. (c) NH3 yields of NPC-850 and NPC-1 at ?0.9 V. Reproduced with permission from Ref. 47. Copyright 2018 American Chemical Society. (d) NRR reaction pathway on pyridinic N3 moiety site. Reproduced with permission from Ref. 48.Copyright 2018 Elsevier.

    Recent studies show that N doping plays a significant role in the NRR performances of CBCs. Density functional theory(DFT) calculations confirm that pyridinic N and pyrrolic N are the key active sites for NRR, and graphitic N is unfavorable for N2adsorption. As shown in Fig. 1a, N-doped porous carbons(NPC) were synthesizedviathe pyrolysis of zeolite imidazolate framework (ZIF) at 750, 850, and 950 °C (denoted as NPC-750,NPC-850, and NPC-950, respectively)47. The obtained catalysts not only preserve ZIF’s features of high surface area and porous structure, but also have abundant active sites enriched by thermally induced defects and charge polarization. Among them,NPC-750 catalyst has the smallest surface area and the highest defect content. In 0.05 mol·L?1H2SO4, it achieves the highest FE of 1.42% and a NH3yield of 1.4 mmol·g?1·h?1at -0.9 V. In this regard, the highest NRR performances of NPC-750 mainly benefit from its intrinsic catalytic activity. XPS results show the highest N content of NPC-750, especially pyridinic (6.2% (x,atomic fraction)) and pyrrolic (5.3% (x)) N. For example as shown in Fig. 1b and 1c, NPC-1 sample of higher content of pyridinic N (2.8% (x)) than NPC-850 (2.1% (x)) has better NRR performance even though its total N content is slightly lower than NPC-850 (5.2% (x)vs. 5.5% (x)). To verify NH3not reduced from the N atoms of catalyst, NPC-750 catalyst is performed under Ar for NRR. The resulting electrolyte shows no NH3,confirming this hypothesis. Wu and co-workers48also reported N-doped nanoporous carbons (C-ZIF)viaheating ZIF-8 at 700–1100 °C for different hours. As temperature and time increasing,N doping content reduces with more disordered carbon structures induced. In 0.1 mol·L?1KOH, the sample thermally treated at 1100 °C for 1 h exhibits a maximum NH3yield of 3.4 ×10?6mol·cm?2·h?1with a Faradic efficiency (FE) of 10.2% at?0.3 V. DFT results reveal that the energetically favorable active sites are thermally induced pyridinic N3moieties embedded in graphitic layer. As shown in Fig. 1d, it consists of one vacancy linked with two pyridinic N and one H linked pyridinic N. The vacancy may result from the removal of one pyridinic N caused by high temperature thermal treatment. It suggests that the vacancy can effectively adsorb and reduce N2molecule with the assistance of pyridinic N. In Song's work49, the catalyst was synthesized from C2H2and NH3gas by a plasma method. The catalyst surface is composed of sharp spikes (~50–80 nm length,~1 nm width), which can concentrate the electric field at tips,promoting N2polarization and electrons injection into the antibonding orbitals of N2. In 0.25 mol·L?1LiClO4, the catalyst achieves a high NH3yield of 97.18 μg?h?1?cm?2with a FE of 11.56% at ?1.19 V. In Wang's study50, N-doped nanoporous graphitic carbon membrane (NCM) catalyze the NRR in 0.1 mol·L?1HCl, attaining a NH3yield of 0.08 g?m?2?h?1and a FE of 5.2%. After functionalized by Au nanoparticles, the catalytic efficiency can be greatly improved by attaining the value of 0.36 g·m?2·h?1and 22%.

    Fig. 2 (a) Schematic of the synthesis process of bamboo shoots-derived NC. (b) HRTEM image of NC-800. Reproduced with permission from Ref. 51. Copyright 2019 Multidisciplinary Digital Publishing Institute.(c) SEM image of NCF-900. (d) Schematic of the synthesis process of cicada slough-derived NCF. Reproduced with permission from Ref. 54.Copyright 2018 The Royal Society of Chemistry. (e) Schematic of the synthesis process of alfalfa-derived NPC. (f) SEM image of NPC-500.Reproduced with permission from Ref. 52. Copyright 2019 American Chemical Society.

    Biomass-derived catalysts also have promising NRR application perspectives because of the extensively accessible,renewable and inexpensive features. As shown in Fig. 2a,b, Liet al.51synthesized N-doped porous carbon (NC) from bamboo shootsviadrying, grinding into powder, hydrothermal carbonization and 600–1000 °C pyrolysis. The obtained NC-800 catalyst (800 °C pyrolysis) shows a graphitic carbon structure with total N content of 3.05% (w, mass fraction). XPS spectra reveal the ratios of pyridinic, pyrrolic and graphitic N of 25.7%,6.7% and 67.6%, respectively. Among them, pyridinic N and pyrrolic N are considered as active sites to activate *N252. The large electronegativity of N (3.04) can effectively induce charge density redistribution and tailor electronic structure53.Meanwhile, the large specific surface area and the hierarchically micro-meso-macroporous structure are favorable for active site exposure and mass transport. In 0.1 mol·L?1HCl, the NC-800 catalyst can promote N2-to-NH3conversion, attaining a NH3yield of 16.3 μg·h?1·mg?1and a FE of 27.5% at ?0.35 V. In Yang’s study54, N-doped hierarchical porous carbon foams(NCF, Fig. 2c,d) are obtained by thermal treatment the ballmilled mixture of cicada slough and ZnCl2, following by acid washing away the Zn2+ions. The ZnCl2plays a role of chemical activating agent and pore-fabricating agent, which can fabricate architecture with micropores (~1.92 nm) and mesopores (~25.3 nm). The NCF-900 catalyst (900 °C thermal treatment) presents a large Brunauer-Emmet-Teller (BET) surface area of 1547.13 m2·g?1. These features are favorable for electrolyte entering and shortened transport distance. Furthermore, the NCF-900 catalyst has high content of pyridinic N (1.96% (x)) and abundant defects as confirmed by XPS analysis and Raman spectra, respectively.In 0.1 mol·L?1HCl, it reaches a NH3yield of 15.7 μg·h?1·mg?1and a FE of 1.45% at ?0.2 V. In Zhao’s study52, N-doped porous graphitic carbon (NPGC, Fig. 2e,f) was synthesized from alfalfa with CaCO3and K2C2O4. The NPGC-500 catalyst (500 °C heat treatment) has 6.35% (x) N content and 3D porous structure composed of micro-, meso- and macropores. The temperatureprogrammed desorption (TPD) measurements confirm the strongest N2adsorption capacity of NPGC-500, which is ascribed to the synergy effect between large surface area and high contents of pyridinic and pyrrolic N. In 0.005 mol·L?1H2SO4, it achieves a NH3yield of 1.31 mmol·h?1·g?1and a FE of 9.98% at ?0.4 V. Furthermore, it has long-term durability though15N isotopic labeling experiments and XPS analysis confirm some pyridinic N escaping from catalyst surface. DFT calculations further reveal that the pyridinic N can be reduced into NH3during NRR, with the generated N vacancies as new active sites. It suggests that both N and C vacancies contribute to the remained NRR catalytic activity.

    3.1.2 O atom doping

    Oxygen has a larger electronegativity of 3.44 than that of carbon of 2.55. Oxygen doping leads to positively charged carbon site by forming polarized C―O bond. This Lewis acid site is beneficial to adsorb N2of weak Lewis base41. Our group developed oxygen-doped graphene (O―G, Fig. 3a)viathe pyrolysis of sodium gluconate with Na2CO3at 950 °C for 10 min under Ar gas flow, followed by acid washing to remove remnant Na2CO355. XPS (Fig. 3d) and FTIR results reveal that O doping into C nanosheets mainly adopt C―O, C=O and O―C=O functional groups. As shown in Fig. 3e, this catalyst achieves a NH3yield of 21.3 μg·h?1·mg?1at ?0.55 V and a FE of 12.6% at?0.45 V in 0.1 mol·L?1HCl. The NRR performances show negligible decrease after 5 cycling tests and 22 h long-term electrolysis, indicating its strong electrochemical stability. As shown in Fig. 3f, three oxygen-doping models are designed to study the NRR performancesviaDFT method. In C―O model,oxygen atom replaces one carbon atom in graphene. Mulliken charge analysis shows that oxygen atom has a negative charge of?0.562 ?e?versusnearby carbon atom of 0.252 ?e?. This electronic structure should not only provide dispersed active sites for NRR by stabilizing N2on positively charged carbon atom, but also suppress HER by prohibiting direct H+adsorption41. DFT results accord well with this assumption. The positively charged carbon atom can effectively adsorb N2, and NRR follows alternative pathway with a low free energy barrier of 0.77 eV from *NN to*NNH. In C=O or O―C=O model, C?OH or ―COOH functional group deviates from graphene plane to form a hole defect. The unsaturated carbon atoms around the hole defect adsorb N2using side-on configuration. DFT calculations suggest that both C=O and O―C=O groups lead to the decrease of free energy barrier of 0.34 eV.

    Fig. 3 (a) TEM image of O-G. Reproduced with permission from Ref. 55. Copyright 2019 The Royal Society of Chemistry. (b) TEM image of O-CN. Reproduced with permission from Ref. 46. Copyright 2019 Wiley-VCH. (c) TEM image of O-KFCNTs. Reproduced with permission from Ref. 39. Copyright 2019 The Royal Society of Chemistry. (d) XPS spectra of O-G in C 1s and O 1s regions. (e) NH3 yields and FEs of O-G at different potentials. (f) Optimized C―O, C=O and O―C=O models and corresponding geometric structures of *NNH. Gray, red, blue and white balls represent the C, O, N and H atoms, respectively. Color online. Reproduced with permission from Ref. 55.Copyright 2019 The Royal Society of Chemistry.

    Oxygen-doped porous carbon nanosheets were also synthesizedviaheating sodium citrate at 750 °C for 1 h under Ar gas flow, followed by acid washing to remove residue solid38.Beside of C―O, C=O and O―C=O functional groups, the catalyst has defect-rich amorphous phase and hierarchically porous structure with an average pore size of 13 nm. As a result,it exhibits a NH3yield of 18.03 μg·h?1·mg?1at ?0.55 V and a FE of 12.6% at ?0.45 V in 0.1 mol·L?1HCl. Furthermore, we developed oxygen-doped carbon materialsviadirectly carbonization of biomass materials under inert gas39,46. The oxygen-doped amorphous nanosheets derived from tannin (OCNs, Fig. 3b) shows NRR performances of a NH3yield of 20.15μg·h?1·mg?1and a FE of 4.97% at ?0.6 V in 0.1 mol·L?1HCl.The catalyst derived from kapok fibers (O-KFCNTs, Fig. 3c) has homogeneous oxygen doping, graphitic carbon and hollow microtube morphology, attaining a NH3yield of 25.12μg·h?1·mg?1and a FE of 9.1%. In addition, the two catalysts have excellent selectivity without N2H4detection in electrolytes and good stability after control experiments.

    3.1.3 S atom doping

    Sulfur is important component element of nitrogenase for biological N2fixation and S doping introduces more defects56.Our group developed S-doped carbon nanosphere (S-CNS, Fig.4a) as a NRR catalystviaa hydrothermal reaction of glucose and benzyl disulfide followed by annealing in Ar atmosphere57. XPS results reveal that S doping into CNS mainly adopts ―C―S―C― and ―C―SOx(x= 2, 3, 4) functional groups. Based on Raman spectra shown in Fig. 4e and f, the calculatedID/IGvalue of 1.16 of S-CNS is larger than the 0.91 of CNS. BecauseIDrepresents the peak intensity of defects whileIGcorresponds to thesp2carbon atom vibrate in plane, it concludes that S doping introduces more defects which may induce abundant active sites to promote N2adsorption and reduction. Compared to pristine CNS, S-CNS offers enhanced N2adsorption capability (Fig. 4d).In 0.1 mol·L?1Na2SO4, S-CNS (19.07 μg·h?1·mg?1) shows higher NH3yield than CNS (3.76 μg·h?1·mg?1).

    This strategy also works effectively for graphene. Our group suggests that S-doped graphene achieves a high NH3yield of 27.3 μg·h?1·mg?1and a FE of 11.5% in 0.1 mol·L?1HCl58. DFT results reveal that the carbon atoms near sulfur atoms are active sites in both S-replace-C model (Fig. 4g) and S-hole model (Fig.4h). The two models of active sites can effectively promote N2-to-NH3conversion, and NH3desorption occurs spontaneously nearly without free energy barrier, suggesting the easy release of active sites for next catalytic cycle. Following our studies, Wang developed S-doped three-dimensional graphene (Fig. 4c) as a metal-free NRR electrocatalyst, capable of attaining a large NH3yield of 38.81 μg·h?1·mg?1with a FE of 7.72% in 0.05 mol·L?1H2SO459. Of note, sulfur dots–graphene nanohybrid (Fig. 4b)also performs efficiently for ambient N2-to-NH3conversion with a NH3yield of 28.56 μg·h?1·mg?1and a FE of 7.07% in 0.5 mol·L?1LiClO460.

    Fig. 4 (a) TEM image of S-CNS. Reproduced with permission from Ref. 57. Copyright 2018 Wiley-VCH. (b) TEM image of sulfur dots-graphene.Reproduced with permission from Ref. 60. Copyright 2019 The Royal Society of Chemistry. (c) TEM image of sulfur-doped 3D-graphene.Reproduced with permission from Ref. 59. Copyright 2020 The Royal Society of Chemistry. (d) N2-TPD curves of S-CNS and CNS. Raman spectra of (e) S-CNS and (f) CNS. Reproduced with permission from Ref. 57. Copyright 2018 Wiley-VCH. Free energy profiles for NRR on sulfurdoped graphene based on (g) model 1 and (h) model 2. Gray and yellow spheres denote C and S atoms, respectively. Related charge density differences of adsorption configurations are also shown. Yellow and blue denote the charge accumulation and depletion, respectively.Reproduced with permission from Ref. 58. Copyright 2019 The Royal Society of Chemistry. Color online.

    3.1.4 B atom doping

    In a pioneering work of Légaréet al., it was found that the boron atom of borylene molecular can fix N2effectively61and the dopant of boron introduces more defects. Previous literatures also suggested that B atom with electron-deficient nature prohibits the binding of Lewis acid H+under acidic conditions62–64. In Yu’s work64, B-doped graphene (BG) was synthesizedviathermal reduction of H3BO3and graphene oxide(GO) at 900 °C in H2/Ar gas. The catalysts with different B doping contents are obtained by controlling the mass ratio of H3BO3to graphene oxide, labeled as BG-1 (5 : 1, Fig. 5a), BG-2 (1 : 10), G (undoped) and BOG (5 : 1, with high oxygen content). XPS results reveal that B-substitution-C (BC3, Fig. 5b)is the main B doping form with BG-1 sample having the highest content, while defect form (B4C) and nanosheet edge defect form(BC2O, BCO2) take only small percentage. As shown in Fig. 5c,the calculatedID/IGvalues of BG-1, BOG, BG-2 and G are 1.3,1.3, 1.1 and 0.91, respectively. Because of the larger electronegativity of carbon (2.55) than boron (2.04), B-doping position has electron-deficient environment, which is beneficial to adsorb Lewis base N2and suppress HER. TPD curves confirm the gradually enhanced N2adsorption capability with increasing B doping content (Fig. 5d). The NRR performances are in line with B doping content, with BG-1 showing the highest NH3yield rate of 9.8 μg·h?1·cm?2and a FE of 10.8% at ?0.5 V in 0.05 mol·L?1H2SO4(Fig. 5e,f).

    Fig. 5 (a) TEM image of BG-1. Reproduced with permission from Ref. 64. Copyright 2018 Elsevier. (b) Schematic of the sp3 orbitals of BC3 for binding N2 (left). LUMO (blue) and HOMO (red) of undoped G (G) and BG (right). The position of a single doped boron atom was labeled (right).Reproduced with permission from Ref. 45. Copyright 2019 American Chemical Society. (c) Raman spectra, (d) N2-TPD curves, (e) NH3 yields and(f) FEs of BG-1, BOG, BG-2 and G. Reproduced with permission from Ref. 45. Copyright 2018 Elsevier.

    B atom has three outermost electrons with an electronic configuration of 2s22p1. As B atom is anchored into carbon materials, one ofselectrons will be excited into aporbital meanwhile with thesandporbitals hybridized by the form ofsp,sp2orsp3. As shown in Fig. 5b, bothspandsp2have unoccupied orbitals and onepempty orbital, which can drive the acceptance-donation process with N2. Andpempty orbital accepts the lone-pair electrons of N2meanwhile unoccupied orbitals donates electrons into the antibonding orbital of N≡N.Liuet al.studied boron-anchored graphyne (GY) for NRRviaDFT method45. The B adsorption on GY configuration shows a poor NRR performance, because the adsorbed B atom has similar electronic structure to free B atom with twopempty orbitals to accept electrons from N2, possessing strong interaction with the N-related intermediates to desorption, e.g.,?6 eV adsorption energy of *NH2. In contrast, the B substitution acetylenic C configuration hassp-hybrid B atom, only with onepempty orbital freely accepting electrons from N2. Its NRR profile shows a free energy barrier only of 0.44 eV at the hydrogenation step from *N to *NH. In fact,sp2- andsp3-hybrid B atoms also have NRR performances, e.g., the above discussed B-doped graphene and the B-decorated graphitic-carbon nitride(B/g-C3N4)64,65. Furthermore, Sun’s group highlights that chemical environment around boron site also plays a key role in NRR process because the abundant electrons in substrate may increase the possibility of electron injection and following reduction63.

    3.1.5 P atom doping

    P has a smaller electronegativity of 2.19 than 2.55 of C,thereby P doping into carbon materials provides an effective strategy to boost NRR performances. P atoms can be successfully incorporated into graphene by P―O and P―C bonds. Our group synthesized P-doped graphene (PG, Fig. 6a)viaannealing of the mixture of GO and triphenyphosphine at 1000 °C in Ar atmosphere66. The calculatedID/IGvalue of 1.07 is higher than 1.02 of pristine graphene. By DFT calculations,the NRR catalytic activities of P―O and P―C binds are compared. The results show that P―O bind has not NRR activity while the P site of P―C bind can effectively promote N2-to-NH3conversion. Density of states (DOS, Fig. 6b) shows that P doping introduces new electron state around Fermi energy, which has a high activity to interact with surface molecule. More importantly, the free energy of *H on P active site is as high as 0.34 eV. The charge variation analysis shows that P atom acts as a role of electron transmitter between nearby three C atoms and adsorbed NxHy. As a result, the P-doped graphene can catalyze NRR along a preferable distal pathway with the first hydrogenation as the rate-limiting step (1.21 eV). As shown in Fig. 6c, it attains a NH3yield rate of 32.33 μg·h?1·mg?1and a FE of 20.82% at ?0.65 V in 0.5 mol·L?1LiClO4.

    3.1.6 F atom doping

    Fig. 6 (a) TEM image of PG. (b) Density of states of P atom in PG and C atom in graphene. (c) NH3 yields and FEs of PG/CP. Reproduced with permission from Ref. 66. Copyright 2020 The Royal Society of Chemistry. (d) HRTEM image of d-FG. (e) Top view of d-FG. (f) NH3 yields and FEs of d-PG/CP with different fluorination time at ?0.7 V. Reproduced with permission from Ref. 69. Copyright 2019 The Royal Society of Chemistry.

    Previous literatures suggested that the strong incorporation of fluorine atoms into graphene might break somesp2carbon bonds and induce defects, thereby providing effective active sites67,68.F atoms can be successfully incorporated into graphene by C?F and C =C―F bonds. Our group synthesized defect-rich fluorographene nanosheet (d-FG, Fig. 6d)viaa hydrothermal reaction of graphene and HF at 180 °C69. The hydrothermal fluorination time can change F doping content which has a significant effect on the catalytic ability of d-FG (Fig. 6f).However, high F content is unfavorable because the large amount F may breaksp2conjugated system and reduce electrical conductivity. The d-FG with 3.12% (x) F doping content presents a defect-rich structure, as confirmed by the highID/IGof Raman spectrum and the pronounced signal (g= 1.918) of electron spin resonance (ESR) spectrum. DFT calculations reveal that the two carbon atoms near F-induced pore defects are active sites (Fig.6e), which adsorb and activate N2effectively by a side-on model with N≡N elongated from 0.112 to 0.126 nm. The first hydrogenation from *N2to *N2H even shows downhill free energy change. The preferable mixed pathway for NRR shows a maximal free energy barrier of 0.37 eV from *N*NH to*N*NH2. In 0.1 mol·L?1Na2SO4, it exhibits good NRR performance and stability, with a NH3yield of 9.3 μg·h?1·mg?1and a FE of 4.2% at ?0.7 V. Liuet al.also developed a F-doped 3D porous carbon framework. Zirconium chloride and terephthalic acid were used to the synthesis of octahedral carbon template, with polytetrafluoroethylene as F source70. Because of the different electronegativity between F (3.98) and C (2.55), F bonding to C atom creates Lewis acid site to facilitate the repulsive interaction between Lewis acid site and H+,suppressing HER. In 0.05 mol·L?1H2SO4, it exhibits a NH3yield of 197.7 μg·h?1·mg?1at ?0.2 V and a FE of 54.8% at ?0.3 V.

    3.1.7 Co-doping

    For electrochemical NRR, N2molecule activation is usually difficult and the first hydrogenation step from *N2to *N2H generally shows uphill free energy (ΔGN2-N2H) as a rate-limiting step in overall pathway71,72. In Chen’s study73, carbon-based catalysts were developed, including graphene (G), N-doped graphene (NG), B-doped graphene (BG), graphene with separated B and N doping (BNG-S), graphene with B?N bonds(BNG-B), BCN with B―N bonds (pyrolyzed from the mixture of PEG-2000, boric acid and urea). It suggests that pristine G and mono-doping G are ineffective for N2molecule activation with positive ΔGN2-N2Hof 2.01, 1.99 and 2.04 eV for G, NG and BG catalysts, respectively. However, BNG-B can activate the N2molecule with ΔGN2-N2Hreduced to 1.72 eV. The incorporation of B―N bonds is the active trigger and the edge C atoms near B―N bonds are the active sites for NRR. The BNG-B and BCN catalysts with abundant B―N bonds can greatly promote the NRR performances with obvious increasement of NH3yields and FEs as shown in Fig. 7a. Moreover, BNG-B and BCN offer higher HER overpotentials versus pristine G and mono-doping G catalysts in Fig. 7b, indicating suppressed HER. DOS calculations confirm that the synergistic effect of B―N bonds raises the highest electron state of catalyst towards Fermi energy level, which triggers the beneficial electron exchange between*N2H and active sites. B―N bonds also play a key role in the NRR performances of B,N co-doped porous carbon nanofiber(B/N-CNF, Fig. 7c), which was fabricatedviaelectrospinning of the mixture of boric acid and polyacrylonitrile (PAN) followed by 250 °C air stabilization and 900 °C NH3annealing74. The obtained nanofiber has a diameter of ~100 nm, consisting of amorphous and graphitic carbons (Fig. 7d,e). This B/N-CNF catalyst possesses abundant defects and vacancies as confirmed by the high ratio ofID/IGof Raman spectrum. Compared with mono-doping CNT, the synergistic effect of B and N of B/NCNF catalyst favors robust NRR partial current density, small charge transfer resistance of electrochemical impedance spectroscopy (EIS), lower Tafel slope and higher electrochemical double layer capacitance (Cdl). It is believed that B atoms contribute to N2adsorption and charge transfer while N atoms enhance electrical conductivity. In 0.1 mol·L?1KOH, the B/N-CNF catalyst attains a NH3yield of 32.5 μg·h?1·mg?1at?0.7 V and a FE of 13.2% at ?0.5 V.

    Fig. 7 (a) NH3 yields, FEs and (b) HER performances of G, BG, NG, BNG-S, BNG-B and BCN. Reproduced with permission from Ref. 73.Copyright 2019 Wiley-VCH. (c) NH3 yields of B/N-CNF, B-CNF and N-CNF. (d) FESEM and (e) HRTEM images of B/N-CNF.Reproduced with permission from Ref. 74. Copyright 2019 The Royal Society of Chemistry.

    The large electronegativity of N (3.04) greatly influences the electron density of adjacent C atoms and leads to beneficial Lewis acid sites for N2adsorption. N, P co-doping induces charge redistribution and creates active sites. N, P co-doped porous carbon (NPPC) also harvests NRR performances with a NH3yield of 0.97 μg·h?1·mg?1and a FE of 4.2% at ?0.2 V in 0.1 mol·L?1HCl75. The NPPC catalyst was synthesizedviaannealing of polyaniline aerogels at 900 °C in Ar atmosphere.TEM results reveal a foam architecture structure with massive graphitic carbons existing in it, indicating beneficial high conductivity. The resulting NPPC catalyst promotes NRR to follow an associative pathway as confirmed byin situFourier transform infrared spectroscopy (FTIR). In Tian’s study, N,S codoped graphene (NSG) was synthesizedviamicrowave treatment (2450 MHZ, 800 W) of thiourea and GO76. It suggests that N, S co-doping is beneficial to the NRR of graphene by promoting N2adsorption. In 0.1 mol·L?1HCl, it achieves a NH3yield of 7.7 μg·h?1·mg?1and a FE of 5.8% at ?0.6 V.

    3.2 Organic polymers

    Conductive polymers such as polyaniline (PAN), polypyrrole(PPY) and polythiophene (PTH) also attract attention in the applications of electrocatalysts77,78. In our study79, PAN was electrochemically deposited on carbon paper (PAN/CP). The existence of pyridinic and pyrrolic N play an important role in the electrocatalytic activity of N-doped carbon frameworks75,80.In 0.1 mol·L?1HCl, the PAN/CP catalyst harvests a NH3yield of 5.45 × 10?11mol·s?1·cm?2at ?0.7 V and a FE of 3.76% at ?0.6 V.In our another study81, perylene-3,4,9,10-tetracarboxylic acid(PTCA) nanorods, which have extendedπ–πelectronic interaction and high electron mobility82,83, were decorated on reduced graphene oxide (rGO) (PTCA-rGO, Fig. 8a).Differential charge density shows obvious charge transfer between PTCA and rGO. Compared with PTCA and rGO, the synergy role of PTCA and rGO greatly promotes N2activation and reduction with reduced free energy barriers. The N2-to-NH3conversion follows both distal and alternating pathways because of the negligible difference of maximum free energy barriers. As shown in Fig. 8b, it attains a high NH3yield of 24.7 μg·h?1·mg?1with FE of 6.9% at ?0.5 V in 0.1 mol·L?1HCl. In Chen’s work84,Li+ions were incorporated into poly(N-ethyl-benzene-1,2,4,5-tetracarboxylic diimide) (PEBCD) by immersing in 0.5 mol·L?1Li2SO4solution. The formed O-Li+plays the role of active site,which greatly retards HER. (Fig. 8e). As shown in Fig. 8d, the PEBCD/C (deposited on carbon cloth, Fig. 8c) catalyst has high HER Tafel slopes of 968.5 and 356.3 mV·dec?1in 0.5 mol·L?1Li2SO4and 0.5 mol·L?1H2SO4, respectively. Benefited from the role of Li+incorporation, PEBCD/C can promote the N2-to-NH3conversion by an alternating pathway, attaining a NH3yield of 1.58 μg·h?1·cm?2and a FE of 2.85% at ?0.5 V.

    Fig. 8 (a) TEM image of PTCA-rGO. (b) NH3 yields of different electrodes. Reproduced with permission from Ref. 81. Copyright 2019 The Royal Society of Chemistry. (c) SEM image of PEBCD. (d) HER Tafel plots of PEBCD/C in 0.5 mol·L?1 Li2SO4 and 0.5 mol·L?1 H2SO4. (e) Schematic of Li+association with O sites of PEBCD. Reproduced with permission from Ref. 84. Copyright 2017 American Chemical Society.

    3.3 Carbon nitride

    Lvet al.42synthesized polymeric carbon nitride (PCN)viaannealing polycondensation melamine at 620 °C in Ar atmosphere. The PCN-NV4 sample (4 h annealing, Fig. 9a) has more nitrogen vacancies (NVs) and smaller transfer resistance as confirmed by electron paramagnetic resonance (EPR, Fig. 9b)spectra and Nyquist plots, respectively. DFT calculations reveal that N2molecule can be activated with N≡N elongated to 0.126 nm after adsorption, and the NRR follows a preferable alternating pathway. Because of its 2D nature, abundant NVs and high conductivity, PCN-NV4 promotes the N2-to-NH3conversion with a NH3yield of 8.09 μg·h?1·mg?1and a FE of 11.59% at ?0.2 V in 0.1 mol·L?1HCl. Many theoretical works were also done on carbon nitrides for NRR, including C2N and C9N4monolayers85–87. In Ji’s work85, B-decorated C2N monolayer was designed as catalyst for NRR. DFT calculations reveal that boron atom can be adsorbed steadily on C2N monolayer by linking to the two N atoms of hole edge (Fig. 9c).The B atom with positive charge and magnetic moment plays a role of active site to activate N2molecule through the acceptance-donation mechanism. As shown in Fig. 9d,enzymatic model is the favorable NRR pathway with an energy barrier of 0.15 eV. In Zhang’s work87, B atom was steadily decorated on the bridge site of C9N4monolayer, by linking to the two N atoms of hole edge (Fig. 9e). DOS calculations indicate the metallic nature of B-C9N4monolayer, which is favorable for promoting charge transfer during NRR (Fig. 9f). The B active site can effectively activate N2molecule and catalyze NRR by following a preferable alternating pathway, with the first hydrogenation as the rate-limiting step.

    3.4 Defect engineering

    Fig. 9 (a) TEM image of PCN-NV4. (b) EPR spectra of PCN and PCN-NV4. Reproduced with permission from Ref. 42. Copyright 2018 Wiley-VCH.(c) B decorated C2N monolayer. (d) Free energy profiles for NRR on B decorated C2N monolayer along enzymatic pathway. Reproduced with permission from Ref. 85. Copyright 2019 The Royal Society of Chemistry. (e) B-C9N4 monolayer. (f) Charge density difference of B-C9N4 monolayer with adsorbed N2 via end-on pattern. Red and green denote the electrons accumulation and depletion, respectively. Reproduced with permission from Ref. 87. Copyright 2020 The Royal Society of Chemistry. Color online.

    Hereinbefore, we summarized the N, O, S, B, P, F atom dopant and co-dopants. Pyridinic N and pyrrolic N can increase the electron-donating property of the catalyst, further indicating the tendency of the electron migration from the carbon atoms to the antibonding orbitals of adsorbed N2. The doping of oxygen atoms could transform hybridization state of the carbon-based materials fromsp2intosp3. And synergistic effect in co-dopants could induce an asymmetrical charge distribution on carbonbased materials, thereby enhancing the electron-donor property of the carbon atoms for weakening the N―N bond and facilitating the N2activation. Besides, defect engineering also is an effective strategy to boost the NRR performances of catalysts.In fact, defects break theπconjugation and play the role of active sites.

    Fig. 10 (a) TEM images of pristine CC and CC-450. (b) NH3 yields and FEs of different electrodes. Reproduced with permission from Ref. 88.Copyright 2018 The Royal Society of Chemistry. (c) Schematic of the synthesis process of layered defective graphene. (d) Raman spectra of PG,DG-700, DG-800 and DG-900. (e) N2 adsorption energies of different catalysts with side-on configuration. (f) TEM image of DG-800. Reproduced with permission from Ref. 89. Copyright 2020 The Royal Society of Chemistry. TEM images of (g) CNT and (h) O-CNT. Reproduced with permission from Ref. 90. Copyright 2019 The Royal Society of Chemistry.

    In Li’s study88, carbon cloth with rich defects was obtained by air annealing at 450 °C (CC-450). Raman spectra of CC-45 show a significantly improvedID/IGof 1.25versusthe values of 1.04,1.15 and 1.21 of pristine CC, CC-250 and CC-350, respectively.As observed in Fig. 10a. HRTEM images confirm the abundant defects of CC-450. Moreover, its hierarchical porous structure contributes to abundant edge defects and better N2adsorption capacity. EIS analysis also shows smaller charge transfer impedance. In NRR tests, CC-450 presents a more positive onset potential in N2-staturated electrolyte than pristine CC, CC-250 and CC-350, indicating its favorable NRR kinetics. As shown in Fig. 10b, it offers the best electrocatalytic activity for NRR, with a NH3yield of 2.59 × 10?10mol·cm?2·s?1and a FE of 6.92% at?0.3 V in 0.1 mol·L?1Na2SO4and 0.02 mol·L?1H2SO4.Graphene is an ideal substrate to load heteroatoms by breaking its integrity ofπconjugation. In Du's study89, bulk graphite oxide is annealed in molten salts of LiCl/KCl (mass ration = 45 :55) at 800 °C to exfoliate graphene by salt permeability and peeling off (DG-800, Fig. 10c). As shown in Fig. 10f, they obtain few-layer ultrathin nanosheets (~3 nm) with defective and porous structures (1–5 nm) which expose more edges and active sites. The abundant defects are identified by broaden XRD peaks, highID/IGof Raman spectra (Fig. 10d) and highsp3/sp2ratio of XPS spectra. Compared with pristine graphene (PG),DG-800 can effectively adsorb N2molecules and catalyze NRR because of point defects and edge defects. As shown in Fig. 10e,the calculated N2adsorption energy shows more negative value at edge defect. In 0.01 mol·L?1H2SO4, DG-800 obtains a NH3yield of 4.31 μg·h?1·mg?1and a FE of 8.51% at ?0.4 V. Chemical oxidation is also an effective way to produce oxygen-containing groups and defects in carbon framework. In our previous work90, carbon nanotube (CNT, Fig. 10g) was simply kept in HNO3at 80 °C for 12–48 h to obtain oxidized CNT (O-CNT,Fig. 10i). The higherID/IGof O-CNT compared to CNT confirms its defective structure. XPS spectra identify the existence of C―OH, ―COOH, C=O and C―O groups. DFT calculations reveal that C?O group is the active site which has the most negative N2-adsorption energy of ?0.32 Vversus?0.06, ?0.08 and ?0.08 of C―OH, ―COOH, C=O and C―O, respectively.When tested in 0.1 mol·L?1LiClO4, the O―CNT catalyst exhibits a high NH3yield of 32.33 μg·h?1·mg?1and a FE of 12.50%.

    3.5 Surface modification engineering

    For the defects engineering mentioned before, carbon materials suffer from high-temperature treatment or strong acid for materials preparation. In fact, there is a surface modification method which can be operated in mild conditions. In our previous work91, rGO was surface modified with oxygen-rich tannic acid (TA)viasimple ultrasonic mixture under ambient conditions. As shown in Fig. 11a,b, TA can greatly disperse rGO from bulk suspension to ultrafine emulsion. TEM images and Raman spectra show no obvious change of intrinsic properties.XPS results (Fig. 11d,e) reveal enhanced O=C peak intensity which is attributed to the ―COOH group of TA attached to rGO. It is believed that theπ–πstacking interaction betweenπrich TA and rGO results in intimate contact of the oxygen groups of TA and rGO, which favors the effective manipulation of the electronic properties of rGO. In 0.5 mol·L?1LiClO4, TA-rGO has better NRR performances, with a NH3yield of 17.02 μg·h?1·mg?1and a FE of 4.83% at ?0.75 V higher than the values of 3μg·h?1·mg?1and 1% of rGO. Furthermore, Fig. 11f,g show the coupled TA-rGO catalyst has good cycling stability and longterm durability.

    Fig. 11 Photos of (a) rGO and (b) TA-rGO in water. (c) TEM image of TA-rGO. C 1s XPS spectra of (d) TA-rGO and (e) rGO. (f) NH3 yields and FEs of TA-rGO at ?0.75 V in N2/Ar-saturated electrolyte. (g) Chronoamperometry curve of TA-rGO/CP at ?0.75 V. Reproduced with permission from Ref. 91. Copyright 2019 American Chemical Society.

    4 Boron-based catalysts

    4.1 Boron atom doping

    Boron has three valence electrons with an electronic configuration of 2s22p1, thereby it has electron-deficient nature with two emptyporbitals to accept foreign electrons. The boron monolayer, composed of triangular and hexagonal motifs, is energetically stable. There are many structures with different hexagon hole densities of 0, 1/4, 1/5, 1/6, 1/7, 1/8 and 1/992–94.Among them, α sheet (1/9) andβ12sheet (1/6) are the stable structures95,96. In Liu's work92, NRR onαsheet andβ12sheet were theoretically investigated. The boron sites with 5 or 6 coordination numbers forαsheet, 4 or 5 coordination numbers forβ12sheet, present catalytic activities for NRR performances.Bader charge analysis shows that N2adsorption is achieved through electron transfer from boron to N2(emptyπ* state) but requiring additional energy input, which is ascribed to the electron-deficient nature of boron. However, the metal substrate(Ag or Au) can facilitate electron injection into boron sheet and reduce the N2adsorption energy barrier. In our previous work97,2D boron nanosheets (BNS, Fig. 12a) were exfoliated from bulk boron powders in isopropyl alcohol using ultrasonic cell pulverizer. The obtained BNS have average thickness of ~3 nm and large BET surface area of 16.38 m2·g?1. Its 2D nature can favor the exposure of active sites. DOS calculations indicate the favorable conductor nature of BNS versus the semiconductor nature of bulk boron. This is in accord with the smaller charge transfer resistance of EIS. It achieves a NH3yield 13.22μg·h?1·mg?1and a FE of 4.04% at ?0.8 V in 0.1 mol·L?1H2SO4(Fig. 12b). By DFT method, 1B, 2B and 3B are considered as N2adsorption sites (Fig. 12c). The calculated results show that 3B atom is the preferable N2active site where an obvious hybridization occurs between *N2and the 2porbitals of B atom.The *N2adsorption free energy at 3B site is identified to be?0.66 eV which is more negative than the *H adsorption free energy of 0.23 eV, which suggests excellent NRR selectivity.Moreover, 1B and 2B sites present strong *H adsorption with the free energies of ?1.67 and ?1.87 eV, indicating that 1B and 2B sites may be deactivated by H atoms. Interestingly, 3B site after the H deactivation of 1B and 2B sites (3B(H), Fig. 12d) shows better N2activation and reduction, with a more negative *N2adsorption free energy of ?0.69 eV. The NRR follows distal pathway with the second NH3desorption as the rate-limiting step. Considering the detected O―H bond of BNS by XPS, the*OH adsorption is also investigated. The results show that the 3B site prefers to be saturated by *OH. However, after this *OH deactivation, the adjacent B atom (Fig. 12e) presents good NRR performances by following alternating pathway.

    4.2 Boron-based compounds

    Fig. 12 (a) TEM image of BNS. (b) NH3 yields and FEs of BNS/CP.Optimal structures and charge density differences of N2 adsorption on(c) 3B site, (d) 3B(H) site and (e) oxidized B (104) surface. Pink, blue,red and white balls represent the B, N, O and H atoms, respectively.Reproduced with permission from Ref. 97. Copyright 2019 American Chemical Society. Color online.

    B4C has focused much attention as electrode materials or catalyst substrates for batteries and fuel cells due to its high mechanical strength, good conductivity and electrochemical stability98–101. In fact, it also has good catalytic activity for NRR.In our work6, B4C bulk was exfoliated in ethanol to form nanosheets using ultrasonic cell pulverize. When tested in 0.1 mol·L?1HCl at ?0.75 V, B4C nanosheets offer good electrochemical stability and achieve a high FE of 15.95% and a NH3yield of 26.57 μg·h?1·mg?1. DFT calculations reveal that*N2adsorption and the first hydrogenation proceed spontaneously at B active site with released free energies,indicating a favorable N2activation. As shown in Fig. 13a, the hydrogenation step of *NH2?*NH2→ *NH2?*NH3is the ratelimiting step with a free energy barrier of 0.34 eV.

    Fig. 13 (a) Free energy profiles of NRR on B4C (110) surface. Reproduced with permission from Ref. 6. Copyright 2018 Springer Nature.(b) Free energy profiles of NRR on the zigzag edge of h-BNNS. (c) Nyquist plots of h-BNNS/CP and bulk h-BN/CP. Reproduced with permission from Ref. 104. Copyright 2019 Springer. (d–g) Structural configurations of N2 bonding to end-on B@ZZBN, side-on B@ZZBN,end-on B@ACBN and side-on B@ACBN. Pink, blue, gray and white balls represent the B, N, C and H atoms, respectively.Reproduced with permission from Ref. 106. Copyright 2019 The Royal Society of Chemistry.

    Hexagonal boron nitride (h-BN) has advantages of high mechanical strength, good thermal conductivity and chemical stability102,103. In our work, h-BN nanosheets (h-BNNS) were exfoliated form bulk h-BN in ethanol using ultrasonic cell pulverize for an anticipation of NRR104. AFM high profiles confirm the layered structure of h-BNNS with thickness of ~1.3 nm. EIS shows a smaller charge transfer resistanceversusbulk h-BN (Fig. 13c). In 0.1 mol·L?1HCl, the h-BNNS catalyst catalyzes NRR by following an enzymatic pathway, attaining a FE of 4.7% and a NH3yield of 22.4 μg·h?1·mg?1at ?0.75 V. The catalytic activity is ascribed to the unsaturated boron atoms at the edges. DFT calculations reveal that N2can be readily adsorbed to the two boron atoms at zigzag edge by forming a 5-membered ring. The unsaturated boron atoms effectively activate and reduce *N2, with the first hydrogenation proceeding spontaneously and the hydrogenation step of *NH*NH →*NH2*NH as the rate-limiting step (Fig. 13b). The catalytic mechanism and kinetics of boron-terminated zigzag (B-ZZ)edges were also confirmed by Leeet al.viatheoretical calculations105. It is believed that the boron atoms at B-ZZ edges act as the active sites for NRR because the chemical environment enables the flexible conversion between B (sp2) and B (sp3)during NRR. In Mao’s study106, BN edge was investigated for NRRviaDFT method. At the nitrogen-terminated zigzag edge of pristine BN, the calculated *N2adsorption energy is as high as ?3.02 eV, which also means difficult NH3desorption (?1.45 eV). Hence, single B atom is decorated at BN zigzag or armchair edge to redistribute the charge density around the active site. The decorated B atom prefers to bond with the two edge N atoms to form a 4-membered ring for zigzag edge (B@ZZBN, Fig. 13d and e), while it bonds at the hollow site to generate 5-membered ring for armchair edge (B@ACBN, Fig. 13f,g). The decorated B atoms act as active sites, and end-on model is favorable for *N2adsorption for both the edge cases.Viacalculating the distal and alternating pathways for both B@ZZBN and B@ACBN, the NRR performances are compared. The results show the distal pathway at B@ACBN is the most favorable way for NRR, with an energy barrier of 0.29 eV from *NH to *NH2. Moreover, the*N2adsorption and NH3desorption free energies for B@ACBN are greatly re-optimized to ?1.97 and 0.85 eV respectively.

    5 Phosphorus-based catalysts

    5.1 P atom doping

    Orthorhombic black phosphorus has a layered structure composed of the monolayers with puckered honeycomb-like top view. Because of the week van der Waals force, black phosphorus can be easily exfoliated into monolayer or few-layer structure, which has extensive applications in electronics, fieldeffect transistor and gas senor107–110. In Zhang’s study111, black phosphorus was exfoliated in ice-bath sonicator to form fewlayer nanosheets (FL-BP NSs). AFM analysis confirms the layered morphology with a thickness of ~4.1 nm (5–7 layers).Benefited from the more exposed active sites, sufficient electrons of 2D nature and weak HER, the FL-BP NSs catalyst has good activity and selectivity for NRR. In 0.01 mol·L?1HCl,it achieves a FE of 5.07% at ?0.6 V and a NH3yield of 33.37μg·h?1·mg?1at ?0.7 V. DFT calculations reveal that the P atoms at zigzag and diff-zigzag edges are active sites, which catalyze NRR by following a preferable alternating pathway.

    Fig. 14 (a) TEM image of BP nanoparticles. (b) NH3 yields and FEs of BP/CP at each potential. Densities of active sites on (c)elementary B and (d) BP. Reproduced with permission from Ref. 115.Copyright 2019 The Royal Society of Chemistry. (e) Free energy profiles of NRR on BP (111) surface. Green and purple balls represent the B and P atoms, respectively. Reproduced with permission from Ref. 113. Copyright 2019 The Royal Society of Chemistry. Color online.

    5.2 Phosphorus compounds

    Because of good stability, high mechanical hardness, modest band gap (2.0 eV) and high dielectric constant, cubic BP attracts much attention as catalyst112,113. In Chen’s study113, NRR performances on the (100), (110) and (111) surfaces of BP were investigatedviaDFT method. The (111) surface has the highest catalytic activity towards NRR by following an energetically enzymatic pathway. As shown in Fig. 14e, the hydrogenation from *NH2to *NH3is the rate-limiting step with an energy barrier only of 0.12 eV, suggesting the robust intrinsic catalytic activity of BP. Moreover, the adsorption energy of *N2H (?3.38 eV) is more negative than that of *H (?1.69 eV) and *H2O(?0.95 eV), indicating favorable N2activation and sluggish HER. The adsorbed *H2O prefers to react with N2molecule with the help of H+and electron to form *N2H (*H2O + N2+ H++e?→ *N2H + H2O) with negative free energy, excluding the H2O poisoning effect and suggesting the favorable role of ionic liquid which can dissolve more N2molecules113,114. Our recent work provided a strong experimental support for BP-enabled N2fixation (Fig. 14a,b)115. In 0.1 mol·L?1HCl, BP nanoparticle affords a high NH3yield of 26.42 μg·h?1·mg?1and a high FE of 12.7% at ?0.6 V. It also has good cycling stability and long-term durability. As shown in Fig. 14c,d, BP is more favorable for the exposure of active sites than the elementary B catalyst.Furthermore, the P atoms can indirectly weaken the N≡N triple bond and strengthen the B―N bond by transferring sufficient electrons to B atoms and then to N atoms. Obviously, elemental phosphorus and phosphorus compounds show promising potential as eligible catalysts for NRR.

    6 Conclusion and outlook

    This review provides an overview of recent progress in designing and developing MFCs towards sustainable and environmentally-benign NH3synthesis at ambient conditions.As can be seen, MFCs show benign NRR electrocatalytic performances (Table 1). Although, extensive attempts and progress have been made in previous works towards MFCs electrocatalysts, compared with TM-based electrocatalysts,MFCs lack the combination of unoccupied and occupied valence orbitals, hindering the electron transfer from catalysts to N2and thus showing poor ability in nitrogen activation. To overcome this limitation, asymmetrical charge distribution on the surface active site of MFCs is effective for promoting the charge exchange between catalysts and N2. Additionally, deeper insights into catalysts design and significant improvement of NRR performances are necessary. In the future studies, the following aspects should be considered:

    (I) Theoretical prediction. The possible active sites,adsorption energies, possible catalytic pathways are generally predicted. These predictions can greatly cut down experimental cost. A high-throughput screening is strongly recommended toquickly select out the effective catalysts. The free energy barriers at the hydrogenation steps of *N2→ *N2H and *NH2→ *NH3are set upper limitation to exclude the sluggish catalysts because these two steps usually have the most positive free energy changes among the whole NRR process for most catalysts.

    Table 1 NRR performances of Metal-free electrocatalysts.

    (II) Promotion of the intrinsic activities of catalysts. An ideal catalyst should enable robust N2adsorption, activation, smooth hydrogenation and NH3desorption. The catalyst kinetics can be modulated by redistributing charge density and tailor electronic structureviaintroducing heteroatoms, vacancies and edges unsaturated atoms,etc. For example, heteroatom of large electronegativity can induce the positive charge of adjacent atoms, which as Lewis acid sites are favorable to adsorb the N2of weak Lewis base, meanwhile these Lewis acid sites can effectively suppress HER by prohibiting binding with the Lewis acid H+41,64. In addition, catalysts after modification of heteroatoms or vacancies,etc. may generate spin magnetic mo ment or raise the highest electron state towards Fermi energy,which may trigger the beneficial electron exchange between *N2and active sites as well as the enhanced electron transfer73. The unsaturated coordination sites are also attractive to N2molecule because these sites mean unpaired electrons which are prone to pair with the foreign electrons.

    (III) Exposing more active sites. Catalysts with hierarchical pore structure and high surface area can expose rich active sites and provide accessible paths for mass/electron transfer. In addition, metal-organic frameworks (MOFs) derived metal-free nanomaterials should be paid more attention as an effective way to boost NRR performances due to their good N2adsorption,adjustable morphologies and tunable coordination atoms.

    (IV) Developments of data reproducibility and accuracy. The amount of ammonia produced is usually so small that it cannot be precisely attributed to NRR. Even though the contemporary methods of analysis are highly sensitive, quantification of the amounts of ammonia at levels that hardly exceed the level of the ubiquitous nitrogen-containing contaminants. Firstly, there are many possible sources of NOxcontaminants in the process of NRR experiments, such as NOxcompounds, NOximpurities and gaseous NOxcompounds that are typically present in the nitrogen gas stream. Secondly, electrocatalysts, electrolytes,solvents, ambient air and any other component of the experimental preparation that is nitrogen-based or was prepared using N-containing precursor can be a source of contamination.

    Therefore, given the contaminations above, it is necessary to quantify the concentrations of NOxcompounds in the N2gas supply before NRR experiments. Among other options, alkaline solutions is a good choice for removal of NOx. In the presence of alkali, HNO2and HNO3, can be neutralized to NO2?and NO3?, impeding the reformation of gaseous NOx. Thus, to remove NOxcontaminations from the cell, electrodes and other labware, washing with alkaline solutions is an efficient method.Chemicals, electrolytes and catalytic materials contaminated with nitrogenous species can be purified through recrystallisation or annealing at an appropriate temperature,while subsequent storage should be under vacuum or clean Ar atmosphere.

    猜你喜歡
    劉倩
    古琴藝術(shù)
    作品(2)
    廣告大觀(2019年8期)2019-09-10 23:26:04
    Rethinking Emotional Branding: Challenges, Risks and Unintended Consequences of Emotional Branding
    商情(2018年10期)2018-03-29 07:14:58
    佳作精賞2
    佳作精賞3
    劉倩、王楠作品
    佳作精賞
    Jamais trop tard
    跟我學(xué)英語
    Un enfant ou deux?
    99国产精品一区二区三区| 麻豆国产97在线/欧美| 亚洲av成人一区二区三| 不卡一级毛片| 国产成人欧美在线观看| 日韩欧美免费精品| 嫩草影视91久久| 亚洲国产日韩欧美精品在线观看 | 每晚都被弄得嗷嗷叫到高潮| 又黄又爽又免费观看的视频| 黄色片一级片一级黄色片| 亚洲熟妇中文字幕五十中出| 三级男女做爰猛烈吃奶摸视频| 欧美一级a爱片免费观看看| 国产成人精品无人区| 日本精品一区二区三区蜜桃| 免费在线观看影片大全网站| 日韩欧美国产在线观看| 国产亚洲av嫩草精品影院| 国产一区二区在线av高清观看| 999久久久国产精品视频| 国产1区2区3区精品| 成年免费大片在线观看| 岛国在线免费视频观看| 午夜免费观看网址| 999久久久精品免费观看国产| 亚洲av成人一区二区三| 国产精品 欧美亚洲| 操出白浆在线播放| 俄罗斯特黄特色一大片| 少妇人妻一区二区三区视频| 黄色日韩在线| 天堂网av新在线| 精品一区二区三区四区五区乱码| 少妇裸体淫交视频免费看高清| 美女免费视频网站| 男人舔奶头视频| 曰老女人黄片| 亚洲 欧美一区二区三区| 热99re8久久精品国产| 亚洲七黄色美女视频| 久久久久久久久久黄片| 此物有八面人人有两片| 成人特级黄色片久久久久久久| 91九色精品人成在线观看| 男女那种视频在线观看| 久久草成人影院| 国产亚洲精品久久久com| 亚洲欧美一区二区三区黑人| 国产三级黄色录像| 宅男免费午夜| 久久国产乱子伦精品免费另类| 在线a可以看的网站| 久久精品国产综合久久久| www.熟女人妻精品国产| 国产单亲对白刺激| aaaaa片日本免费| 青草久久国产| 亚洲黑人精品在线| 中文字幕熟女人妻在线| 香蕉av资源在线| av欧美777| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲专区中文字幕在线| 午夜视频精品福利| 女同久久另类99精品国产91| 日韩国内少妇激情av| 99精品在免费线老司机午夜| 欧美成人免费av一区二区三区| 久久性视频一级片| 18禁美女被吸乳视频| 九色成人免费人妻av| 男人舔女人下体高潮全视频| 国产av在哪里看| 日本成人三级电影网站| 久久精品亚洲精品国产色婷小说| 最新美女视频免费是黄的| 久久久久国产一级毛片高清牌| 亚洲精品久久国产高清桃花| 欧美成人性av电影在线观看| 久久精品影院6| 亚洲国产看品久久| 久久久国产欧美日韩av| 一进一出抽搐gif免费好疼| 色老头精品视频在线观看| 中文资源天堂在线| 小蜜桃在线观看免费完整版高清| 亚洲在线观看片| 亚洲国产精品久久男人天堂| 欧美极品一区二区三区四区| 99久久精品国产亚洲精品| 中国美女看黄片| 亚洲七黄色美女视频| 中文字幕熟女人妻在线| 亚洲中文av在线| 色吧在线观看| 国产91精品成人一区二区三区| 在线视频色国产色| 国产精品免费一区二区三区在线| 午夜福利18| 欧美又色又爽又黄视频| ponron亚洲| 国产一区二区三区视频了| 黑人巨大精品欧美一区二区mp4| 亚洲色图av天堂| 又黄又爽又免费观看的视频| 国产亚洲精品综合一区在线观看| 久久久精品大字幕| 黄色丝袜av网址大全| 国产视频内射| 免费在线观看成人毛片| 欧美日韩一级在线毛片| 国产成人影院久久av| 成人午夜高清在线视频| 我要搜黄色片| 91av网站免费观看| 亚洲欧美日韩卡通动漫| 久久久精品欧美日韩精品| 哪里可以看免费的av片| 在线视频色国产色| 婷婷精品国产亚洲av在线| 久久久久久人人人人人| 老司机在亚洲福利影院| 在线十欧美十亚洲十日本专区| 国产免费男女视频| 久久精品亚洲精品国产色婷小说| 色综合婷婷激情| 少妇的逼水好多| 看片在线看免费视频| 久久热在线av| 久久久久国产一级毛片高清牌| 两人在一起打扑克的视频| 国内精品久久久久久久电影| 国产欧美日韩一区二区精品| 在线观看午夜福利视频| av视频在线观看入口| 性色av乱码一区二区三区2| 嫩草影院精品99| 亚洲狠狠婷婷综合久久图片| 床上黄色一级片| 一二三四社区在线视频社区8| 中文字幕人妻丝袜一区二区| 每晚都被弄得嗷嗷叫到高潮| 99久久综合精品五月天人人| 中文字幕高清在线视频| 亚洲一区二区三区色噜噜| 欧美日韩一级在线毛片| 久久久久久大精品| 国产精品一区二区三区四区久久| 美女黄网站色视频| 精品乱码久久久久久99久播| 免费电影在线观看免费观看| 日韩欧美三级三区| 美女高潮喷水抽搐中文字幕| 久久精品国产清高在天天线| 国产精品,欧美在线| 亚洲av美国av| 99精品欧美一区二区三区四区| 99久久精品热视频| 在线免费观看的www视频| 国产伦精品一区二区三区视频9 | 亚洲成a人片在线一区二区| av片东京热男人的天堂| 久久精品国产99精品国产亚洲性色| 国产成人一区二区三区免费视频网站| 亚洲人成电影免费在线| 国产成人一区二区三区免费视频网站| 亚洲国产色片| 亚洲人与动物交配视频| 午夜激情福利司机影院| 在线看三级毛片| 色视频www国产| 精品99又大又爽又粗少妇毛片 | 国产淫片久久久久久久久 | 国产又黄又爽又无遮挡在线| 最新美女视频免费是黄的| 欧美最黄视频在线播放免费| 成人欧美大片| 亚洲欧美日韩高清在线视频| 欧美一级a爱片免费观看看| 欧美性猛交黑人性爽| 欧美zozozo另类| 91麻豆av在线| 日本黄大片高清| 亚洲九九香蕉| 久久国产精品人妻蜜桃| 欧美+亚洲+日韩+国产| 国产精品自产拍在线观看55亚洲| 精品久久久久久成人av| 999久久久精品免费观看国产| 极品教师在线免费播放| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| 成在线人永久免费视频| 精品久久久久久久久久久久久| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频| 91在线观看av| 老汉色av国产亚洲站长工具| 亚洲国产中文字幕在线视频| 欧美日本亚洲视频在线播放| 十八禁人妻一区二区| 一级毛片高清免费大全| 岛国视频午夜一区免费看| 最新在线观看一区二区三区| 久久久久久人人人人人| 十八禁人妻一区二区| 国产熟女xx| 久久精品国产清高在天天线| 人妻久久中文字幕网| 成人av一区二区三区在线看| 久久精品91无色码中文字幕| 久久久久免费精品人妻一区二区| 九色成人免费人妻av| 欧美日韩亚洲国产一区二区在线观看| 日韩精品青青久久久久久| 久久精品91无色码中文字幕| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 国产 一区 欧美 日韩| av天堂在线播放| 欧美成人一区二区免费高清观看 | 国产黄色小视频在线观看| 九色国产91popny在线| 亚洲成人精品中文字幕电影| 99精品久久久久人妻精品| 亚洲午夜精品一区,二区,三区| 国产黄a三级三级三级人| 成熟少妇高潮喷水视频| 日本与韩国留学比较| 国产精华一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久午夜电影| 淫妇啪啪啪对白视频| 中文字幕人成人乱码亚洲影| 国产成人aa在线观看| 国产精品1区2区在线观看.| 琪琪午夜伦伦电影理论片6080| 欧美中文综合在线视频| 97人妻精品一区二区三区麻豆| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 亚洲精品美女久久av网站| 蜜桃久久精品国产亚洲av| 国产三级中文精品| 99精品久久久久人妻精品| 国产黄片美女视频| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 美女黄网站色视频| 亚洲第一电影网av| 午夜a级毛片| 亚洲人与动物交配视频| 麻豆久久精品国产亚洲av| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 亚洲性夜色夜夜综合| 国产成人av激情在线播放| 精华霜和精华液先用哪个| 熟女人妻精品中文字幕| 亚洲 欧美 日韩 在线 免费| 18美女黄网站色大片免费观看| 久久久久免费精品人妻一区二区| 亚洲无线在线观看| 一级毛片女人18水好多| 亚洲精品色激情综合| 中出人妻视频一区二区| 亚洲乱码一区二区免费版| 久久九九热精品免费| 成年版毛片免费区| 国产成人精品无人区| 亚洲国产日韩欧美精品在线观看 | 欧美色视频一区免费| 曰老女人黄片| 岛国在线免费视频观看| 老司机深夜福利视频在线观看| 国产一区二区在线av高清观看| 精品久久久久久久久久久久久| 色在线成人网| 桃色一区二区三区在线观看| 婷婷精品国产亚洲av在线| 男人和女人高潮做爰伦理| 国产精品98久久久久久宅男小说| 久久中文字幕一级| 人人妻人人看人人澡| 日本精品一区二区三区蜜桃| 男女那种视频在线观看| 18禁国产床啪视频网站| 久久久久精品国产欧美久久久| 美女大奶头视频| 美女 人体艺术 gogo| 国产精品久久久久久久电影 | www.精华液| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩高清专用| 精品一区二区三区视频在线观看免费| 又大又爽又粗| 又紧又爽又黄一区二区| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 免费人成视频x8x8入口观看| 亚洲午夜精品一区,二区,三区| 国内毛片毛片毛片毛片毛片| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 真实男女啪啪啪动态图| 香蕉av资源在线| 色av中文字幕| 亚洲精品一区av在线观看| 久9热在线精品视频| 高清毛片免费观看视频网站| 热99re8久久精品国产| 身体一侧抽搐| 亚洲成av人片免费观看| 国产精品久久久久久久电影 | aaaaa片日本免费| 国产精品综合久久久久久久免费| 别揉我奶头~嗯~啊~动态视频| 熟女电影av网| 十八禁人妻一区二区| 婷婷精品国产亚洲av| 一个人免费在线观看电影 | 夜夜爽天天搞| 国产精品1区2区在线观看.| av视频在线观看入口| 国产精品99久久久久久久久| 性色av乱码一区二区三区2| 在线国产一区二区在线| 90打野战视频偷拍视频| 久久精品国产综合久久久| 琪琪午夜伦伦电影理论片6080| 一级毛片女人18水好多| 国产成人福利小说| 国产av麻豆久久久久久久| 波多野结衣高清作品| 黄色丝袜av网址大全| 热99在线观看视频| 伦理电影免费视频| 亚洲片人在线观看| 丰满人妻熟妇乱又伦精品不卡| 男女下面进入的视频免费午夜| 99久国产av精品| 午夜福利18| 两人在一起打扑克的视频| 精品电影一区二区在线| 午夜激情欧美在线| 人人妻,人人澡人人爽秒播| bbb黄色大片| 日韩免费av在线播放| 最新在线观看一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲人成伊人成综合网2020| 最好的美女福利视频网| 少妇的丰满在线观看| 日韩成人在线观看一区二区三区| 性欧美人与动物交配| 国产三级在线视频| 狂野欧美激情性xxxx| 国产男靠女视频免费网站| 三级毛片av免费| 亚洲 国产 在线| 久99久视频精品免费| 悠悠久久av| 黄片大片在线免费观看| 亚洲一区二区三区色噜噜| 久久九九热精品免费| 午夜福利在线观看免费完整高清在 | 亚洲av成人精品一区久久| 国产私拍福利视频在线观看| 91老司机精品| 亚洲精品粉嫩美女一区| 亚洲av片天天在线观看| 欧美日韩乱码在线| 激情在线观看视频在线高清| 日本 欧美在线| 婷婷精品国产亚洲av| 精品久久久久久久末码| 精品国产亚洲在线| 久久久国产成人免费| 成人av一区二区三区在线看| 久久这里只有精品19| 又粗又爽又猛毛片免费看| 热99re8久久精品国产| 波多野结衣高清无吗| 亚洲欧美激情综合另类| 一级a爱片免费观看的视频| 亚洲欧美精品综合久久99| 国产又色又爽无遮挡免费看| 亚洲第一电影网av| 五月玫瑰六月丁香| 一级毛片高清免费大全| 日韩欧美免费精品| 国产精品一区二区三区四区免费观看 | 免费看十八禁软件| 成年女人看的毛片在线观看| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 夜夜夜夜夜久久久久| 国产精品久久久av美女十八| 国产精品 欧美亚洲| 国产精品一及| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 欧美乱码精品一区二区三区| 亚洲欧美一区二区三区黑人| 久久国产精品影院| 麻豆国产av国片精品| 日韩免费av在线播放| x7x7x7水蜜桃| 18禁黄网站禁片免费观看直播| 亚洲avbb在线观看| 久久精品国产清高在天天线| 啦啦啦韩国在线观看视频| 国产一区二区在线av高清观看| 色综合婷婷激情| 91av网站免费观看| 亚洲欧美精品综合一区二区三区| 两性夫妻黄色片| 久久久久久九九精品二区国产| 午夜激情福利司机影院| 色在线成人网| 色尼玛亚洲综合影院| 宅男免费午夜| 看片在线看免费视频| 欧美性猛交╳xxx乱大交人| 亚洲一区二区三区不卡视频| 99国产精品一区二区蜜桃av| 日韩中文字幕欧美一区二区| bbb黄色大片| 午夜福利18| 校园春色视频在线观看| 国产私拍福利视频在线观看| 丁香六月欧美| 日本免费a在线| 老熟妇乱子伦视频在线观看| 久久性视频一级片| 嫩草影视91久久| 女人被狂操c到高潮| 国产不卡一卡二| 日韩欧美国产在线观看| 麻豆国产97在线/欧美| 欧美日本视频| 天天躁日日操中文字幕| 1000部很黄的大片| 日本黄色片子视频| 久久99热这里只有精品18| 国产 一区 欧美 日韩| 国产成人福利小说| 老司机深夜福利视频在线观看| 久久久久久人人人人人| 精品欧美国产一区二区三| av黄色大香蕉| 国产三级在线视频| 欧美国产日韩亚洲一区| 老汉色av国产亚洲站长工具| 国产午夜精品论理片| 夜夜爽天天搞| 亚洲成人久久性| xxxwww97欧美| 黄色女人牲交| 99久久精品一区二区三区| 免费观看人在逋| 网址你懂的国产日韩在线| 国产精华一区二区三区| 亚洲精品在线观看二区| 欧美中文综合在线视频| 长腿黑丝高跟| 国产免费男女视频| 俺也久久电影网| 18禁观看日本| 免费看日本二区| 丁香欧美五月| 99热这里只有是精品50| 天天添夜夜摸| 一级黄色大片毛片| 亚洲人与动物交配视频| 色噜噜av男人的天堂激情| 真实男女啪啪啪动态图| 久久久精品欧美日韩精品| 手机成人av网站| 给我免费播放毛片高清在线观看| 女人高潮潮喷娇喘18禁视频| 成人18禁在线播放| 久久香蕉国产精品| 757午夜福利合集在线观看| 十八禁人妻一区二区| 男人舔女人的私密视频| 精品日产1卡2卡| 美女高潮喷水抽搐中文字幕| 搡老熟女国产l中国老女人| 99久久无色码亚洲精品果冻| 国产高潮美女av| 97超级碰碰碰精品色视频在线观看| 久久亚洲精品不卡| 一级毛片精品| 亚洲精品中文字幕一二三四区| 亚洲中文字幕日韩| av片东京热男人的天堂| 久久精品国产清高在天天线| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放| 亚洲中文字幕一区二区三区有码在线看 | 国产毛片a区久久久久| 久久九九热精品免费| 最好的美女福利视频网| 亚洲天堂国产精品一区在线| 亚洲av电影在线进入| 色老头精品视频在线观看| tocl精华| 亚洲成人中文字幕在线播放| 久久久精品大字幕| 国产黄色小视频在线观看| 全区人妻精品视频| 美女高潮的动态| 国产91精品成人一区二区三区| 日韩欧美在线乱码| 欧美成人一区二区免费高清观看 | 国产精品久久电影中文字幕| 久久香蕉国产精品| 手机成人av网站| 在线观看日韩欧美| 99re在线观看精品视频| 成人一区二区视频在线观看| 亚洲欧美一区二区三区黑人| av在线蜜桃| 婷婷丁香在线五月| 久久久水蜜桃国产精品网| 欧美绝顶高潮抽搐喷水| 亚洲五月婷婷丁香| 日韩欧美精品v在线| 欧美大码av| av在线天堂中文字幕| 国产亚洲精品久久久久久毛片| 日本撒尿小便嘘嘘汇集6| 日韩欧美一区二区三区在线观看| 国产亚洲精品综合一区在线观看| 亚洲精品久久国产高清桃花| 国产精品久久久久久亚洲av鲁大| av黄色大香蕉| 国产精品av久久久久免费| 小说图片视频综合网站| 97超级碰碰碰精品色视频在线观看| 亚洲激情在线av| 国产熟女xx| 日本一二三区视频观看| 欧美午夜高清在线| 香蕉久久夜色| 久久久色成人| 国产又色又爽无遮挡免费看| 99久久精品一区二区三区| 男女视频在线观看网站免费| 日韩人妻高清精品专区| 日本在线视频免费播放| or卡值多少钱| 不卡av一区二区三区| 国产淫片久久久久久久久 | 激情在线观看视频在线高清| 一进一出抽搐动态| 日本 av在线| 一级毛片精品| 久久天躁狠狠躁夜夜2o2o| 日本五十路高清| 两个人看的免费小视频| 一个人免费在线观看电影 | 国产v大片淫在线免费观看| 亚洲精品美女久久av网站| 老熟妇乱子伦视频在线观看| 无人区码免费观看不卡| 他把我摸到了高潮在线观看| 国产午夜精品论理片| 一本一本综合久久| 男女下面进入的视频免费午夜| 精品久久久久久久人妻蜜臀av| 18美女黄网站色大片免费观看| 国产av在哪里看| 精品久久久久久久久久免费视频| av视频在线观看入口| 一区二区三区高清视频在线| 精品久久久久久成人av| 中文字幕熟女人妻在线| 国产极品精品免费视频能看的| 精品久久蜜臀av无| 美女黄网站色视频| 757午夜福利合集在线观看| 国产美女午夜福利| 久久精品国产清高在天天线| 特级一级黄色大片| 欧美日韩一级在线毛片| 99久久精品热视频| 久久精品亚洲精品国产色婷小说| 麻豆成人午夜福利视频| 久久99热这里只有精品18| 一本综合久久免费| 国产真实乱freesex| 国产精品久久久人人做人人爽| 欧美日韩国产亚洲二区| 国产淫片久久久久久久久 | 欧美日韩综合久久久久久 | 国产精品综合久久久久久久免费| www.熟女人妻精品国产| 制服丝袜大香蕉在线| 午夜福利高清视频| 欧美大码av| 亚洲av熟女| 又粗又爽又猛毛片免费看| 男女那种视频在线观看| 免费人成视频x8x8入口观看| 亚洲真实伦在线观看| 亚洲av成人精品一区久久| 亚洲黑人精品在线| 欧美又色又爽又黄视频| 午夜精品一区二区三区免费看| 亚洲国产欧美一区二区综合| 久久久国产欧美日韩av| 看免费av毛片|