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

    原子尺度鈷基氮碳催化劑對析氧反應的構效關系的研究

    2023-02-17 03:55:30吳明亮章燁暉付戰(zhàn)照呂之陽李強王金蘭
    物理化學學報 2023年1期
    關鍵詞:金蘭東南大學工程學院

    吳明亮,章燁暉,付戰(zhàn)照,呂之陽,李強,*,王金蘭,*

    1東南大學物理學院,南京 211189

    2東南大學機械工程學院,江蘇省微納生物醫(yī)療器械設計與制造重點實驗室,南京 211189

    1 Introduction

    With the increase in fossil oil consumption and worldwide attention to the environment, developing sustainable, clean, and renewable energy has become one of the most critical challenges1–3. Storing excess electrical energy in chemical bonds, such as hydrogen production by water splitting driven by electricity, is a promising strategy to alleviate the energy crisis4.In the water splitting process, oxygen evolution reaction (OER)as a half-reaction, involves a multi-step electron transfer process,thermodynamic climbing, and slow kinetic process, which dramatically limits the total water splitting efficiency5,6. At present, nanomaterials based on IrO2/RuO2are generally considered to be the most advanced OER electrocatalysts7. Still and all, their large-scale application is hindered by high cost, low natural reserves, and relatively poor durability8. Thus,exploiting OER electrocatalysts with high-efficiency and low cost is significantly indispensable9–11.

    Recently, transition metal-based (TM, Co, Ni, Fe) M-N-C catalysts have attracted extensive studies due to their low cost,high stability/conductivity, and flexible electronic structure of the active center12–14. Compared with traditional metal oxides bulk, nanoparticles, M-N-C catalyst with the concept of single atom catalysts (SACs) provides the most outstanding atomic efficiency15,16. Moreover, with the help of advanced characterization techniques like X-ray absorption near edge spectroscopy (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy techniques17, the fine structure of M-N-C catalysts can be further explored. The characterization of this kind of structure can bridge the experiment and theory,which is beneficial to further exploring the active origin of the reactive sites. However, the activity of electrocatalytic reactions is highly structurally sensitive, and slight changes in the coordination environment can lead to large changes in the activity18,19. Meanwhile, the dynamic evolution of the structure during the reaction can also lead to changes in reactivity20. This makes it difficult to understand the role played by the active sites to modulate the activity, which further limits the improvement and new design of active catalysts. Among the widely studied metal SACs, the cobalt-based M-N-C catalyst show high oxygen evolution reaction activity. For example, Zhang et al. fabricated an atomically dispersed cobalt- and nitrogen-codoped graphene catalyst with the configurations of Co-N4and Co-N2C2, showing an ultra-low onset overpotential of approximately 210 mV21.Zeng et al. theoretically predicted that the SACs with a Copyrrole-N4center could efficiently catalyze OER compared with pyridine nitrogen coordination22. Lu et al. found that cobalt nanoparticles@N-doped carbon composite exhibited a low onset overpotential of 300 mV for OER, which is similar to the performance of reported single-atom Co-N-C catalysts23.Therefore, the Co-N-C structure is chosen as an example system in this work to study the structure-activity relationship12,24–26.

    In this study, we investigate the Co-N-C catalysts with different configurations (single atom, diatoms, and cluster) to determine the source of OER activity by free energy calculations and polarization curve simulations. The results show that the coordination environment of metal atoms affects the valence state of metal atoms significantly and plays a key role in the reactivity. Furthermore, the effects of long-range spin coupling and short-range weak metal-metal bonding on OER activity are discussed based on the Cu-N4configuration. By replacing the vicinal metal to regulate the weak metal-metal bonding, the extrapolated overpotential of 0.23 V from the OER volcano curve can be obtained, which is comparable to the noble metal catalysts. This study explores the active sites and influencing factors of OER activity through a series of Co-N-C systems,which guide the design of high-efficiency OER electrocatalysts.

    2 Computational methods

    A Co and N co-doped graphene system was constructed using a supercell containing a (6 × 6) graphene sheet, in which the Co doped atoms replaced two adjacent C atoms, and the nitrogen atoms were introduced to replace the carbon atoms, directly connecting with the Co dopant atom. Our calculations in this work are done through the Vienna ab initio Simulation Package(VASP) code27,28. The interaction of ions-electrons is described by the projector augmented wave (PAW) method29. The electron exchange and correlation were approximated by Perdew-Burke-Ernzerhof (PBE)30, and the kinetic energy cutoff of the planewave basis was set to 520 eV. In addition, since the weak interaction cannot be accurately described by the standard PBE functional, the DFT-D3 scheme based Grimme empirical atom pair correction is used to describe the long-range van der Waals interaction31,32. The Brillouin zone was sampled using the Γ-centered k-point grid with the sizes of 5 × 5 × 1 and 7 × 7 × 1 kmeshes for structure relaxation and electronic properties calculation, respectively. A 16 ? (1 ? = 0.1 nm) vacuum along z-direction was employed to avoid the interaction between two periodic units. All atoms were fully relaxed until the energies and residual forces on each atom converged to 10-5eV and 0.02 eV·?-1, respectively, during geometry optimization. The charge transfer between adsorbates and the catalyst was quantitatively described by Bader charge analysis33,34. The climbing image nudged elastic band (CI-NEB) method is used to search for the transition state for the elementary reaction35. With help from the program lobster, the chemical bonding analysis of the metaladsorbate interaction was carried out using the crystal orbital Hamilton population (COHP)36,37.

    Following the kinetic model developed by Hansen et al.38, we simulate the polarization curve of Co-N-C. The O electrochemical oxidation steps are listed by the following equations:

    where the adsorption sites derived from *O2(aq) represent O2in the electrolyte, specific information is given in Supporting Information.

    3 Results and discussion

    3.1 Single atom

    The identification of the active sites, which establishes a definitive correlation between the atomic structures and catalytic performance, plays a key role in the improvement and new design of single atom catalysts (SACs). Recent advances in experimental studies reported that the SACs can be identified as,mostly, an M-N4configuration by sophisticated techniques. For example, Ban et al. used the electron energy loss spectroscopy(EELS) analysis results of C, N, and Co to prove that more Co-N4was induced during the structural differentiation process17.Meanwhile, Wang et al. calculated the formation energy of the Co-N-C structure and found that with the increase of the N dopant concentration, the formation energy decreased, resulting in a higher possibility of observing the most stable Co-N4in the experiment39. However, the current density of the Co-N-C catalyst in the experiments under low overpotential conditions is relatively low, about 5 to 20 mA·cm-2, which means that a small number of active sites can reach limited values and thus the existence of other configurations as the “real” active sites cannot be ruled out21. Therefore, we collect possible configurations for single atom sites based on reported studies, including Co-N3, Co-N1C3, Co-N2C2, Co-N3C1, Co-N3-v and Co-N4structures (Fig.1a)21,39,40. As proposed by prior studies, a four-step OER mechanism was assumed that proceeds through OH*, O*,OOH*, and O2* (the adsorption site denoted with the asterisk)41.The overpotential corresponding to the OER catalyst is an essential parameter for evaluating the catalytic activity, and it is determined by the free energy of the rate-determining step(RDS)12. Several structures are calculated and compared with the existing data, the consistent results validate our computational models and methods21,40. The detailed calculation results for the elemental steps of all the model structures are given in Table S1 (in Supporting Information).

    Fig. 1 Evaluation of catalytic activity of single atoms with different structures. (a) Monoatomic structure of cobalt with different coordination structures, where black, light blue and blue represent C, N and Co atoms, respectively (color online). (b) The relationship between the overpotential of different Co-N-C catalysts and ΔG*OH, where ? is referenced in Ref. 40. (c) Simulated polarization curves of Co-N4,Co-N4-exp, Co-N3C1, Co-N3-v and Co-N1C3, data of Co-N4-exp is referenced in Ref. 21.

    The OH* intermediate, created upon the first electron and proton transfer from water, is expected to initiate the formation of the first chemical bond of O2(O―O), with corresponding free energy differences ΔG*OH. Fig. 1b is obtained according to the overpotential of the Co-N-C catalyst and the ΔG*OH, which was used as a predictive theoretical descriptor of OER activity in previous studies42. In general, the overpotential gradually increases as the bonding strength with *OH increases (ΔG*OHgradually decreases). Such trend corresponds to the increase in the adsorption capacity of the catalyst to certain oxygencontaining intermediates. Based on Sabatier principle, an ideal electrocatalyst needs moderate binding with all the involved intermediates to approach the theoretical value of 1.23 V for each reaction step (4.92 V/4)43. It can be found that Co-N4has the highest ΔG*OHof 0.70 eV, which is closer to the theoretical value compared with other configurations, and its overpotential is the lowest overpotential of 0.37 V among all the investigated configurations. As ΔG*OHdecreases, the overpotentials of Co-N3C1, Co-N3-v, and Co-N1C3are 0.64, 0.87 and 0.90 V,respectively (Fig. S1). The difference in these systems can be essentially explained by the valence state of the central Co atom.Based on Bader charge analysis of Co-N3-v, Co-N3C1, and Co-N4(Table S2), each coordinated N atom gets 0.45e on average,and the C atom gets 0.19e. The Co atom in Co-N3-v loses 0.76e,and the Co atom in Co-N3C1and Co-N4loses 0.80e and 0.89e,respectively, Consistently, it can be found that the adsorption strength of Co-N3C1for oxygen-containing intermediates is slightly stronger than Co-N4. Similarly, for Co-N3-v, the coordination atom number reduction increases the valence state of Co and further increase the adsorption strength of oxygencontaining intermediates on Co-N3-v.

    To further understand the influence of Co-N-C coordination structure on reaction overpotential, four reported typical configurations such as Co-N4, Co-N3-v, Co-N3C1and Co-N1C3are selected for further study, and the overpotentials follow the order of Co-N4< Co-N3C1< Co-N3-v < Co-N1C3. From Fig. S1,it can be found that the rate-determining step (RDS) of Co-N4,RDS is OH* oxidized to *O, and the RDS of other catalysts is the oxidation of O* to OOH*12. From the orbital perspective of*OOH formation, an orbital of σ character (highest occupied molecular orbital of OH-) is close to the M-O π* orbital (lowest unoccupied molecular orbital). The attaction of OH-to MO bond leads to the formation of an O―O σ bond in *OOH and the breakage of one M=O π bonds (Fig. S2). In that case, it is requried that the energy levels of M―O π* orbital and OH-σ should be as close as possible to form a new O―O bond.Therefore, the adsorption of *O can be neither too strong(decrease the energy level of M―O π* orbital ) nor too weak(increase the energy level of M―O π* orbital)44. Among these catalysts, the ΔG*OHand ΔGO*of the Co-N1C3system are the lowest. From the adsorption energy diagram (Fig. S3), it is found that the chemical adsorption of the corresponding intermediates is also strong. This leads to nucleophilic OH-attacking the strongly adsorbed *O to form *OOH with a higher energy barrier and overpotential.

    As mentioned above, the Co-N4structure shows ΔG*OHof 0.70 eV, which is much less than the theoretical value of 1.23 eV, suggesting the activity can be improved by modulating the charge state of the metal center and the binding energies with the oxygen-containing intermediates. For example, He et al.enhanced the intrinsic activity by adding additional graphitic N doping surrounding the Co-N4moieties, the essence of which is to increase the N doping concentration near the metal center,increase the valence state of the metal center, and promote moderate adsorption of reaction intermediates45. It is also worth noting that the overpotentials in the OER are different at the same N concentration. For instance, Zhang et al. calculated the overpotential of three CoN2C2structures, Co-N2-O (1.81 V), Co-N2-F (1.69 V) and Co-N2-S (1.36 V)39. Thus, the chemical environments of the active metal atoms, including coordination number and local composition, are critical to the activity of SACs.

    To give an intuitive demonstration of the OER catalytic performance for Co-N-C, we simulated its polarization curve and compared it with experimental values. At a current density of 10 mA·cm-2, the overpotential of Co-N4is 0.33 V, which was much lower than Co-N3C1(η10= 0.49 V), Co-N3-v (η10= 0.71V)and Co-N1C3(η10= 0.77 V). The activity trend of Co-N-C is consistent with our previous calculation results. It can be observed that the polarization curve of Co-N4structure is closest to the experimental value (η10= 0.39 V)21, thus the dominant SAC sites can be attributed to such a strong coordinated configuration. However, the existence of other configurations,such as Co-N3C1, cannot be ruled out, the gap between the experiment and the theoretical value could be because the active site during the experiment is a mixture of multiple Co singleatom configurations. To sum up, the single-atom structure of Co-N4in Co-N-C has effective coordination, moderate adsorption strength for OER reaction intermediates and the lowest reaction overpotential among all the investigated SAC structures.

    3.2 Diatoms and clusters

    Due to the lack of precise atomic position control in current synthesis methods, and the catalyst is accompanied by structural reconstruction in the electrocatalytic process, it is possible that M atoms migrate and agglomerate into metal clusters during experimental processes46. For diatoms, we design four structures based on reported experimental results, i.e., Co2-type1, Co2-type2, Co2-type3 and Co2-type4 respectively, as shown in Fig. S4. Considering the three-coordination structure of the Co atom in Co2-type1, we compared it with fourcoordination Co-N4and clusters of atomic states to find the nature of OER overpotential change caused by the structure. In addition, as the distance between active sites decreases, there are interactions between metal centers, and the distance between active centers gradually decreases from Co2-type2 (dCo-Co= 4.98 ?) to Co2-type3 (dCo-Co= 4.08 ?) to Co2-type4 (dCo-Co= 2.25 ?).Therefore, we further study the metal-metal interactions by taking these three diatomic structures for comparison.

    Firstly, we further study the oxygen evolution activity of Co-N-C structures containing a metal dimer and four-atom cluster,denoted as Co2-type1 and Co4-type1 (Fig. 2a), respectively. The corresponding OER free energy diagrams are calculated and presented in Fig. 2b. From the free energy diagram, compared with Co-N4, we find that Co2-type1 and Co4-type1 have strong adsorption capacity for the reaction intermediates, leading to considerably large overpotentials of 1.82 V and 0.96 V,respectively, which are much higher than that of Co-N4(0.37 V).It can be noted that free energy of the first step is an energy decreasing process, indicating that the adsorption of *OH can be a thermodynamically favorable process. The high affinity of the Co dimers and clusters towards oxygen intermediates can be attributed to the low oxidation state of the metal atoms. As proved by Bader charge analysis, the Co atom in Co-N4loses 0.89e (0.89e-), which is much more than Co2-type1 (0.67e-) and Co4-type1 (0.10e-). Accordingly, the accumulated charge around the Fermi level of Co2-type1 and Co4-type1 (Fig. S5) can be transferred to oxygen atoms effectively, forming strong chemical bonding between Co atoms and oxygen-containing adsorbates.

    The above results also suggest that the Co atoms in dimers and clusters present high reactivity and can be unstable in the current form during OER processes. Alternatively, the free energy changes indicate that metal sites can be adsorbed and stabilized by *OH and *O with accessible potentials. In that case, the “real”catalytic models can be changed to the structures with adsorbed*OH and *O, as the following OER processes is restricted for pristine Co-N-C catalysts due to large overpotentials. Wang et al. recently discovered through DFT calculation and experimental methods that during the ORR reaction, Co2-type1 would spontaneously form *OH as the modified ligand and bind to atomic metal sites, which has higher overall ORR catalytic activity compared with traditional metal surfaces24.Subsequently, the catalytic performance of Co2-type1 and Co4-type1 after pre-adsorption reaction intermediates (such as O,OH, O2) are evaluated. The results show that the oxygen evolution activity is significantly improved by the preadsorption, the calculated overpotentials drops largely from 1.82 V to 0.48 V, 0.40 V and 0.52 V for *OH, *O and *O2preadsorbed Co2-type1 structure. Similarly, the overpotential of Co4-type1-O (Fig. S6) decreases from 0.96 V to 0.40 V with preadsorbed *O species. Meanwhile, the polarization curves of the catalysts after pre-adsorption were simulated (Fig. 2d). At a current density of 10 mA·cm-2, Co2-type1-O (η10= 0.42 V), Co2-type1-OH (η10= 0.42 V) and Co2- type1-O2(η10= 0.46 V) is similar to Co-N4(η10= 0.33 V). It can be noticed that the overpotentials of the cobalt dimers and clusters are significantly improved after pre-adsorption, where the overpotentials of Co2-type1-O are close to the singly occupied Co-N4structures. In addition, the breakage of O=O and desorption of O2are also considered. The calculation results show that O2adsorption can stabilize on the metal dimer and cluster (Fig. S7 and Table S3),the second O2molecule can desorb (binding energy -0.15 eV)without affecting OER activity or causing catalyst oxidation and blocking metal sites after Co2-type1 pre-adsorbs one O2molecule (Fig. S8).

    Fig. 2 Evaluation of catalytic activity of diatoms and clusters. (a) Structures for Co2-type1 and Co4-type1. (b) Reaction energy diagram of Co-N4, Co2-type1, and Co4-type1. (c) Reaction energy diagram of Co2-type1 after pre-adsorption of OH, O and O2 (Co2-type1-OH,Co2-type1-O, Co2-type1-O2).(d) Simulated polarization curves of Co-N4, Co2-type1-OH, Co2-type1-O and Co2-type1-O2.

    The pre-adsorbed models can also be compared to cobalt(oxygen) hydroxide, which have shown promising OER activity in recent studies47,48. and the high-valence cobalt Co3+(β-CoOOH) is considered to be the active catalytic center of OER process49. The high-valence Co in the (oxygen) hydroxide exhibits the t2g5eg1electronic configuration of the intermediate spin state, and the egorbital directly participates in the σ-bonding between adsorbed species on the surface, the appropriate egorbital filling is beneficial to the oxygen-containing intermediates in the process of oxygen electrocatalysis (such as OH*, O* and OOH*) interact50. Therefore, the catalytic activity of Co-based catalysts is highly depending on the high coordination number and the oxidation states of the metal sites,rather than composition forms (single atom or metal clusters).

    3.3 Indirect metal-metal interactions

    According to the above results, the high coordination is a prerequisite for the OER process and Co-N4presents the moderate adsorption energies for each reaction intermediates.Due to the uneven distribution of reaction sites and different metal doping concentration applied in the experiments, when the distance between active centers decreases, there can be interactions between two adjacent reaction sites, such as the long-range spin coupling effects between metal centers and weak metal-metal bonding51,52. Considering the spin nature of electrons in the OER mechanism, the adsorption of OER reaction intermediates at the metal center can be strongly regulated by the spin state, which has been extensively studied53.However, the effect of inter-site spin coupling is less discussed,thus several reported configurations with different Co-Co distances are considered to study the long-range interactions.Firstly, Co2-type2 and Co2-type3 were established to study the influence of spin coupling between adjacent metal centers on catalytic performance (Fig. 3a), the distances between the two Co atoms are 4.96 ? and 4.08 ?, respectively. As shown in Fig.3a, the energy differences of Co2-type2 and Co2-type3 between ferrimagnetic and anti-ferrimagnetic (energy favorable state) are 60 meV and 156 meV, respectively, indicating that the spin coupling effect becomes more pronounced as the distance decreases. Comparing the change in reaction free energy (Fig.3b), the overpotential of Co2-type2 with the antiferromagnetic state is similar to that of Co-N4, 0.39 V and 0.37 V, respectively.On the contrary, the overpotential of Co2-type3 (ηOER= 0.48 V)is higher than that of Co-N4(ηOER= 0.37 V) and the RDS of Co2-type3 (oxidizes O* to OOH*) is different from Co-N4(oxidizes OH* to O*).

    Fig. 3 Identification of the indirect inter-site interaction. (a) Spin density of Co2-type2 (Ⅰ), Co2-type3 (Ⅱ ), CoFe -type3 (Ⅲ ) and CoNi- type3 (Ⅳ).The yellow or blue color indicates positive and negative spin density. The isosurface is plotted with 0.003 eV·?-3. (b) Reaction energy diagram of Co-N4 (blue), Co2-type2 (red) and Co2-type3 (black). (c) Adsorption energies of oxygen-containing intermediates (*OH, *O, *OOH) on CoFe-type3, Co2-type3 and CoNi-type3, respectively. Color online.

    To further study the effect of spin coupling between active centers on the activity of OER catalyst, we replaced one of the two Co atoms with the strong magnetic Fe and weak magnetic Ni. The structure and spin density are shown in Fig. 3a, and the spin configurations are presented in Fig. S9. Then, the OER process was calculated at the Co site for comparison (Fig. S10).The free energy diagrams of CoFe-type3, Co2-type3, and CoNitype3 show that the overpotential and RDS of OER have no apparent change after the replacement. Compared with the adsorption energy of the reaction intermediates of Co-Fe, Co-Co, Co-Ni in Fig. 3d, there is no significant change, which shows that the spin coupling between active centers has a minor effect on the activity of the reaction centers. To clarify the difference between Co2-type3 and Co-N4, Bader charge analysis is performed (Table S1) and the results show that the charge of the N atom located between the two Co centers is 0.45e-lower than the surrounding N atoms of Co-N4, indicating the local environment is changed. Therefore, the difference in RDS and OER performance between Co2-type3 and Co-N4can be attributed to the influence of the coordination environment,rather than the spin coupling between the two reaction centers.

    As the distance between two Co-N4structures approaches and forms a closely connected structure, i.e., Co2-type4 (Fig. 4a), the short distance can induce metal-metal bonding between two adjacent active centers thus have influence on the catalytic performance. The bonding is expected to be weak as the high coordination nature of the metal sites. As from the electron localization functions (ELF) contour analysis (Fig. 4a, bulk Co is shown in Fig. S11), it can be regarded that there is a weak bond formation between the two Co atoms and the electrons (density)are mainly localized on the surrounding N atoms. The free energy diagram of Co2-type4 is calculated and compared with Co-N4, as seen in Fig. 4b. The overpotential of Co2-type4(yellow) is 0.36 V, which is the same as Co-N4(black), 0.37 V.However, the free energy change of Co2-type4 on each step is different from Co-N4. To further determine the effect of the weak metal-metal bonding, a Co atom is replaced in Co2-type4 with Mn, Fe, Ni, and Cu, and the free energy diagram of the oxygen evolution reaction at the same Co site is calculated (Fig. 4b). The result shows that when the adjacent metal changes, the change of the RDS for the type4 configuration mainly depends on the species O*. Moreover, ΔGO*changes obviously in a linear way and the adsorption strength of O* decreases in the order of Mn >Fe > Co > Ni > Cu (Fig. S12). In addition, we considered the influence of the pre-adsorption intermediates at the ortho metal site on the performance of the Co site, and the results showed thatortho metals of type4 configuration can adsorb oxygencontaining intermediates, but the effect on the OER performance of the Co metal atoms is negligible (Fig. S13).

    Fig. 4 Identification of the weak metal-metal interaction. (a) Structures for Co2-type4 and the electron localization functions (ELF) contour.(b) Reaction energy diagram of Co-N4, CoMn-type4, CoFe-type4, Co2-type4, CoNi-type4, and CoCu-type4. (c) The scaling relationship between ΔG*OOH and ΔG*OH. (d) The scaling relationship between ηOER versus ΔGO* - ΔG*OH. The theoretical volcano plot derived from the data is displayed in black. Color online.

    To further understand the change rule of O binding strength on Co atom, we carried out Crystal Orbital Hamilton Population(COHP) analysis. The projected COHP (pCOHP) and the integrated COHP (ICOHP) for the Co―O bond in O-Co2-type4,O-CoNi-type4, and O-CoCu-type4 systems with consistent adsorption sites are calculated (Fig. S14). In the figure -pCOHP is plotted so that the population on the right of the x axis indicates the bonding state and those on the left indicate the antibonding state. The ICOHP is obtained by integrating pCOHP with Ef,which is a measure of bonding strength. The more negative the ICOHP value is, the stronger the Co―O bond will be. From the COHP analysis, we note that from Co2-type4 to CoCu-type4, the part of the antibonding state above the Fermi level gradually increases, and the antibonding state is less filled, which enhances the adsorption of O*. The ICOHP value of Co2-type4 is -0.76 eV, and its Co―O bond strength is the weakest (Fig. S14).Followed by CoNi-type4 (-0.91 eV) and CoCu-type4 (-0.96 eV), the Co-O bond is strengthened in turn, which is consistent with the above analysis. Therefore, the weak metal-metal bonding can influence the OER activity and the adsorption energy of O* can be regulated to optimize the OER performance,especially the hetero-metals.

    According to previous findings, the activity of OER is strongly related to ΔGO*; thus, it can be used as an adequate descriptor to predict the overall OER reactivity, which involves multiple intermediates. Subsequently, the minimum overpotential can be predicted by establishing the scaling relationship between the adsorption free energies of different adsorption species. As shown in Fig. 4c, the linear relationship between ΔGOOH*and ΔG*OHis ΔGOOH*= 0.69ΔGO*+ 3.32, and the coefficient of determination (R2) is 0.945, indicating that a fine linear relationship is obtained. This correlation shows that the catalytic activity of OER can be described by the free energy of adsorption of one or two intermediates. By plotting the relationship between the oxygen evolution overpotential of different Co-N-C structures and ΔGO*- ΔG*OH, it is found that the binding strength of the adsorbed species, as well as the overpotential, is in a volcanic shape as shown in Fig. 4d. Besides,it is worth mentioning that the RDS of the system on the left side of the x axis 1.47 eV is *O oxidized to *OOH, and the RDS of the system on the right of 1.47 eV is the oxidation of *OH to *O.The reason is that the system on the left has strong adsorption on oxygen-containing intermediates, and the stronger adsorption of OH also activates OH, which reduces the energy barrier in the step of *OH oxidation to *O. However, the strong adsorption of O* on metal sites leads to the nucleophilic ―OH attacking *O,forming *OOH with a high energy barrier (―OH combines with O―M to form O―O―σ bond and break one of the M―O π bonds at the same time), which limits the catalytic process44.The higher ΔGO*value on the right side of 1.47 V on the x axis in Fig. 4d indicates that the adsorption capacity of *OH on the catalyst is weak, which is not conducive to OH activation,resulting in a higher energy barrier for *OH dissociation from H+to form O*. For all the calculated systems, CoNi-type4 lies on the top of the volcano curve and present the best OER performance, consistent with the free energy calculation results for all the elementary steps. Besides, the volcano curve shows that the catalytic activity of OER reaches its peak when ΔGO*is close to 1.47 eV, and the minimum ηOERis 0.23 V, which exceeds the limit overpotential of 0.37 V in the volcano curve of previous OER catalysts, breaking the OOH to OH ratio relationship54,55.

    4 Conclusion

    In this work, a series of Co-N-C catalysts with different configurations (single atom, diatoms, and cluster) have been established to determine the origin of the OER activity. Based on in-depth analyses including overpotential, adsorption energy,PDOS, Bader charge and COHP, the results demonstrate that the high coordination of metal sites with a high valence state is required to ensure the moderate adsorption strength of the reaction intermediates, which makes Co-N4as the best catalyst among all the established models. Metal dimers and clusters need to be further re-constructed by the interaction with catalyst supports (interface engineering) and the adsorption of solvation ions and reaction intermediates. Besides, the binding strength of the reaction intermediates can be also modulated by the indirect change of the coordination environment, which endows the change of the micro-environment of the catalysts to be one of the common ways to promote the reaction activity. Moreover,hetero-atom doping with weak metal-metal interactions is proven to be another way to promote reaction activity. In a combination of these strategies, a volcano curve is established by associating the overpotential of different Co-N-C structures with ΔGO*, and an ultralow overpotential of 0.23 V is obtained,which could break the theoretical optimum value (0.37 V). This study explores the structure-activity relationship on OER through the refined Co-N-C structures, which bridges the gap between the experiments and theoretical calculations and thus provides a theoretical basis for the rational design of highefficiency OER electrocatalysts.

    Acknowledgements: The authors thank the computational resources from the Big Data Center of Southeast University and National Supercomputing Center of Tianjin.

    Supporting Information:available free of charge via the internet at http://www.whxb.pku.edu.cn.

    猜你喜歡
    金蘭東南大學工程學院
    福建工程學院
    福建工程學院
    《東南大學學報(醫(yī)學版)》稿約
    《東南大學學報(醫(yī)學版)》稿約
    《東南大學學報(醫(yī)學版)》稿約
    《東南大學學報(醫(yī)學版)》稿約
    福建工程學院
    福建工程學院
    萬有引力定律
    完美火炬手
    久久精品国产鲁丝片午夜精品 | 亚洲精品一区av在线观看| 18禁黄网站禁片免费观看直播| 长腿黑丝高跟| 毛片女人毛片| 老师上课跳d突然被开到最大视频| 日日摸夜夜添夜夜添av毛片 | 九九爱精品视频在线观看| 国产中年淑女户外野战色| 日本一本二区三区精品| 日韩中文字幕欧美一区二区| 搡女人真爽免费视频火全软件 | 亚洲精品久久国产高清桃花| 久久久久精品国产欧美久久久| 欧美绝顶高潮抽搐喷水| 亚洲成人中文字幕在线播放| 国产黄色小视频在线观看| 国产精品美女特级片免费视频播放器| 国产精品电影一区二区三区| 18禁在线播放成人免费| 中文字幕精品亚洲无线码一区| 人妻丰满熟妇av一区二区三区| 在线观看免费视频日本深夜| 日韩欧美三级三区| 欧美色视频一区免费| 国产综合懂色| 午夜福利在线观看免费完整高清在 | 国产精品亚洲美女久久久| 在线观看av片永久免费下载| 亚洲最大成人av| 久99久视频精品免费| 亚洲三级黄色毛片| 国产精品综合久久久久久久免费| 精品一区二区三区av网在线观看| 成年女人看的毛片在线观看| 亚洲国产欧洲综合997久久,| 99精品久久久久人妻精品| 久久久久久久亚洲中文字幕| 国产精品一区二区性色av| 高清日韩中文字幕在线| 亚洲av免费高清在线观看| 国产又黄又爽又无遮挡在线| 三级男女做爰猛烈吃奶摸视频| avwww免费| 搡老妇女老女人老熟妇| 色av中文字幕| 成年版毛片免费区| 精品一区二区免费观看| 一本一本综合久久| 色综合婷婷激情| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| 国产成人av教育| av专区在线播放| 人妻夜夜爽99麻豆av| 久久中文看片网| 久久久久九九精品影院| 免费观看精品视频网站| 亚洲 国产 在线| 男女边吃奶边做爰视频| 能在线免费观看的黄片| 亚洲av中文字字幕乱码综合| 2021天堂中文幕一二区在线观| 日韩强制内射视频| 成年免费大片在线观看| 国产毛片a区久久久久| www日本黄色视频网| 国产伦人伦偷精品视频| 精品日产1卡2卡| 99久久九九国产精品国产免费| 男人舔女人下体高潮全视频| 久久久久国内视频| 免费观看精品视频网站| 国产精品福利在线免费观看| 日韩欧美免费精品| 精品午夜福利视频在线观看一区| 白带黄色成豆腐渣| 美女黄网站色视频| 亚洲成人久久爱视频| 国产av麻豆久久久久久久| 亚洲精品亚洲一区二区| а√天堂www在线а√下载| 婷婷丁香在线五月| 欧美区成人在线视频| 日韩欧美国产一区二区入口| 2021天堂中文幕一二区在线观| 国产av在哪里看| 亚洲av成人av| 国产精品美女特级片免费视频播放器| 久久久久久久久中文| 最好的美女福利视频网| 亚洲国产精品成人综合色| 成熟少妇高潮喷水视频| 欧美日韩黄片免| 精品久久久久久久久av| 狂野欧美白嫩少妇大欣赏| 国产精品人妻久久久久久| 久久久午夜欧美精品| 深爱激情五月婷婷| 高清在线国产一区| 不卡一级毛片| 一级毛片久久久久久久久女| 久久精品国产亚洲av天美| 免费av不卡在线播放| 一进一出好大好爽视频| 午夜爱爱视频在线播放| 午夜视频国产福利| 欧美bdsm另类| 99在线视频只有这里精品首页| 成人一区二区视频在线观看| 波多野结衣高清无吗| or卡值多少钱| 人妻久久中文字幕网| 麻豆成人av在线观看| 免费高清视频大片| 美女 人体艺术 gogo| 搡老熟女国产l中国老女人| 久久国产乱子免费精品| 亚洲最大成人手机在线| 熟妇人妻久久中文字幕3abv| 免费av毛片视频| 久久国内精品自在自线图片| netflix在线观看网站| 村上凉子中文字幕在线| 国产亚洲av嫩草精品影院| 日日夜夜操网爽| 男人舔奶头视频| 欧美日韩乱码在线| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 国产精品三级大全| 亚洲人成网站在线播| 又紧又爽又黄一区二区| 欧美一区二区国产精品久久精品| 成人亚洲精品av一区二区| 婷婷精品国产亚洲av| 国产精品国产三级国产av玫瑰| 99精品久久久久人妻精品| 尾随美女入室| 搡女人真爽免费视频火全软件 | 最新在线观看一区二区三区| 亚洲中文字幕日韩| 我要搜黄色片| av在线老鸭窝| 日本与韩国留学比较| 国产高清有码在线观看视频| 级片在线观看| 亚洲自拍偷在线| 色精品久久人妻99蜜桃| 国产亚洲精品av在线| 日韩中字成人| 美女cb高潮喷水在线观看| 可以在线观看毛片的网站| 亚洲真实伦在线观看| 日本爱情动作片www.在线观看 | 在线天堂最新版资源| 精品久久久久久久久av| 亚洲精品乱码久久久v下载方式| 成人特级黄色片久久久久久久| 有码 亚洲区| 精品无人区乱码1区二区| 免费看光身美女| 99久久精品热视频| 成年免费大片在线观看| 成年女人看的毛片在线观看| 欧美成人性av电影在线观看| 九色国产91popny在线| ponron亚洲| 亚洲av一区综合| 亚洲图色成人| 嫩草影视91久久| 亚洲七黄色美女视频| 中文在线观看免费www的网站| 日韩人妻高清精品专区| 午夜老司机福利剧场| 国内精品宾馆在线| 国产亚洲精品av在线| 午夜福利欧美成人| 在线免费十八禁| 国产探花极品一区二区| 精品久久久久久久久av| av黄色大香蕉| 日本免费一区二区三区高清不卡| 日本a在线网址| 午夜久久久久精精品| 男人舔奶头视频| 国产精品久久久久久av不卡| 99久久精品一区二区三区| 亚洲精品日韩av片在线观看| 久久午夜福利片| 成人欧美大片| 精华霜和精华液先用哪个| 嫁个100分男人电影在线观看| 看片在线看免费视频| av国产免费在线观看| 男插女下体视频免费在线播放| 中文字幕高清在线视频| 精品乱码久久久久久99久播| 欧美成人免费av一区二区三区| 中文亚洲av片在线观看爽| 亚洲av日韩精品久久久久久密| 一卡2卡三卡四卡精品乱码亚洲| 人妻夜夜爽99麻豆av| а√天堂www在线а√下载| 国产成人a区在线观看| 最好的美女福利视频网| 春色校园在线视频观看| 伊人久久精品亚洲午夜| 丰满乱子伦码专区| 国产高清三级在线| avwww免费| 亚洲人与动物交配视频| 成人国产综合亚洲| 亚洲av成人av| 人妻久久中文字幕网| 一区二区三区高清视频在线| 俺也久久电影网| 日韩欧美在线乱码| 精品一区二区免费观看| 极品教师在线视频| 日韩精品青青久久久久久| 最近最新免费中文字幕在线| 国产伦精品一区二区三区四那| 亚洲精品成人久久久久久| 日韩亚洲欧美综合| 在线观看66精品国产| 中文字幕熟女人妻在线| 精品一区二区三区av网在线观看| 黄色女人牲交| 免费看日本二区| 999久久久精品免费观看国产| 天美传媒精品一区二区| 亚洲午夜理论影院| 99在线人妻在线中文字幕| 人妻久久中文字幕网| 999久久久精品免费观看国产| 国产一区二区激情短视频| 一边摸一边抽搐一进一小说| 亚洲成a人片在线一区二区| 精品久久久噜噜| 中文资源天堂在线| 黄色配什么色好看| 国产白丝娇喘喷水9色精品| 波多野结衣高清无吗| 欧美性猛交╳xxx乱大交人| 午夜爱爱视频在线播放| 人妻丰满熟妇av一区二区三区| 国产高潮美女av| 琪琪午夜伦伦电影理论片6080| 亚洲人成伊人成综合网2020| netflix在线观看网站| av在线观看视频网站免费| 我要搜黄色片| 欧美日韩亚洲国产一区二区在线观看| 欧美不卡视频在线免费观看| 精品国产三级普通话版| 午夜老司机福利剧场| 亚洲在线自拍视频| 22中文网久久字幕| 美女xxoo啪啪120秒动态图| 精品日产1卡2卡| 国产熟女欧美一区二区| 九色成人免费人妻av| 老司机午夜福利在线观看视频| av在线蜜桃| 日韩av在线大香蕉| 亚洲,欧美,日韩| 亚洲av成人av| 亚洲在线观看片| 男女下面进入的视频免费午夜| 欧美xxxx性猛交bbbb| 日本黄色视频三级网站网址| 欧美zozozo另类| 日韩在线高清观看一区二区三区 | 美女高潮的动态| 国产精品综合久久久久久久免费| 日韩亚洲欧美综合| 变态另类成人亚洲欧美熟女| 少妇猛男粗大的猛烈进出视频 | 国产激情偷乱视频一区二区| 亚洲av成人av| 欧美一级a爱片免费观看看| 日韩精品有码人妻一区| 久久草成人影院| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 日本熟妇午夜| 国产精品女同一区二区软件 | 日本a在线网址| 午夜日韩欧美国产| 国产亚洲91精品色在线| 亚洲在线自拍视频| 日韩欧美一区二区三区在线观看| 亚洲久久久久久中文字幕| 午夜影院日韩av| 亚洲国产色片| 成人av一区二区三区在线看| www.www免费av| 久久久久国产精品人妻aⅴ院| 欧美性猛交黑人性爽| 又黄又爽又刺激的免费视频.| 一区福利在线观看| 日本黄色片子视频| 亚洲精品一区av在线观看| 久99久视频精品免费| av.在线天堂| 久久亚洲真实| 色综合色国产| 国产高清有码在线观看视频| 亚洲最大成人中文| 美女xxoo啪啪120秒动态图| 国产激情偷乱视频一区二区| 久久久久久久午夜电影| 成人欧美大片| 男女边吃奶边做爰视频| 免费一级毛片在线播放高清视频| АⅤ资源中文在线天堂| 精品久久久噜噜| 国产私拍福利视频在线观看| 久久人妻av系列| 午夜a级毛片| 特大巨黑吊av在线直播| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 两个人的视频大全免费| 简卡轻食公司| 中文字幕av成人在线电影| 亚洲四区av| 国产精品久久久久久久电影| 色精品久久人妻99蜜桃| 白带黄色成豆腐渣| 嫩草影院入口| 波野结衣二区三区在线| 中国美女看黄片| 少妇高潮的动态图| 久久久久久久精品吃奶| 深爱激情五月婷婷| 91狼人影院| 久久久久久久午夜电影| 18+在线观看网站| 亚洲内射少妇av| 精品久久久久久久久久免费视频| 波多野结衣高清作品| 日韩欧美国产在线观看| 国产精品国产三级国产av玫瑰| 亚洲图色成人| 91麻豆精品激情在线观看国产| 欧美高清性xxxxhd video| 日本色播在线视频| 欧美日韩亚洲国产一区二区在线观看| 久久亚洲精品不卡| 久久久久精品国产欧美久久久| 久久精品国产亚洲av香蕉五月| 99热这里只有精品一区| 国产主播在线观看一区二区| 国产久久久一区二区三区| 国产不卡一卡二| 女的被弄到高潮叫床怎么办 | 老师上课跳d突然被开到最大视频| 久久99热这里只有精品18| 亚洲av电影不卡..在线观看| 久久精品久久久久久噜噜老黄 | 少妇的逼好多水| 人妻丰满熟妇av一区二区三区| 国产高清视频在线观看网站| 成人综合一区亚洲| 国产一区二区亚洲精品在线观看| 亚洲精华国产精华液的使用体验 | 国产av不卡久久| 精品久久久久久久久久免费视频| 精品久久久久久久久av| 在线国产一区二区在线| 亚洲性夜色夜夜综合| 18+在线观看网站| 日韩强制内射视频| 国产免费男女视频| а√天堂www在线а√下载| 天天一区二区日本电影三级| 国产精品久久久久久av不卡| 熟女电影av网| 亚洲国产色片| 日日啪夜夜撸| 成年女人毛片免费观看观看9| 亚洲美女黄片视频| 成人亚洲精品av一区二区| 一区福利在线观看| 老女人水多毛片| 深夜精品福利| 精品久久久久久久人妻蜜臀av| 有码 亚洲区| 波多野结衣巨乳人妻| 噜噜噜噜噜久久久久久91| 韩国av一区二区三区四区| 我要搜黄色片| 亚洲av五月六月丁香网| 一本久久中文字幕| 国产一区二区三区视频了| 1024手机看黄色片| 日韩欧美国产一区二区入口| 国产精品伦人一区二区| 亚洲美女视频黄频| 日韩欧美国产在线观看| 国产国拍精品亚洲av在线观看| 99久久精品一区二区三区| 综合色av麻豆| 欧美人与善性xxx| 日日干狠狠操夜夜爽| 在线国产一区二区在线| 国产亚洲91精品色在线| 日本在线视频免费播放| 日本免费一区二区三区高清不卡| 亚洲精品久久国产高清桃花| 精品一区二区三区视频在线观看免费| 熟女人妻精品中文字幕| 国产三级在线视频| 如何舔出高潮| 国产精品一及| 日日摸夜夜添夜夜添小说| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 一夜夜www| 又黄又爽又免费观看的视频| 在线观看午夜福利视频| 日韩欧美 国产精品| 精品人妻熟女av久视频| 国产精品福利在线免费观看| 亚洲人成网站在线播| 国产精品一及| 麻豆一二三区av精品| 嫩草影院精品99| 日韩国内少妇激情av| 欧美黑人巨大hd| 亚洲成人免费电影在线观看| 麻豆成人午夜福利视频| 亚洲美女黄片视频| 色尼玛亚洲综合影院| 色吧在线观看| 可以在线观看毛片的网站| 嫩草影院精品99| 看免费成人av毛片| 久久久久九九精品影院| 91午夜精品亚洲一区二区三区 | 日韩欧美一区二区三区在线观看| 人妻丰满熟妇av一区二区三区| 91久久精品国产一区二区三区| 搞女人的毛片| 久久天躁狠狠躁夜夜2o2o| 熟女电影av网| 午夜激情欧美在线| 国产一区二区在线观看日韩| 中文在线观看免费www的网站| 亚洲熟妇熟女久久| 国产成人一区二区在线| 亚洲精华国产精华精| 日本黄大片高清| 午夜激情福利司机影院| 在线观看一区二区三区| 国产色爽女视频免费观看| 99在线人妻在线中文字幕| 在线观看舔阴道视频| 成人鲁丝片一二三区免费| 国产精品无大码| 日本在线视频免费播放| 九色成人免费人妻av| 久久久久免费精品人妻一区二区| 欧美黑人巨大hd| 国产精品一区二区性色av| 在线看三级毛片| avwww免费| 精品一区二区三区av网在线观看| 国产精品精品国产色婷婷| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 日韩欧美国产一区二区入口| 国产精品久久电影中文字幕| 深夜精品福利| 亚洲性久久影院| 一个人观看的视频www高清免费观看| 精品一区二区三区视频在线观看免费| 日本撒尿小便嘘嘘汇集6| 精品国产三级普通话版| 成人一区二区视频在线观看| 日本a在线网址| 亚洲成人免费电影在线观看| 免费在线观看成人毛片| 久久这里只有精品中国| 国产欧美日韩精品一区二区| 两个人视频免费观看高清| 久久精品91蜜桃| 在线看三级毛片| 久久这里只有精品中国| 99久久精品一区二区三区| 狠狠狠狠99中文字幕| 国产精品亚洲一级av第二区| 亚洲人与动物交配视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精品一区二区免费欧美| 免费看美女性在线毛片视频| 在线a可以看的网站| 国产精品美女特级片免费视频播放器| 3wmmmm亚洲av在线观看| 成人综合一区亚洲| 男女之事视频高清在线观看| 永久网站在线| 男人狂女人下面高潮的视频| 美女高潮喷水抽搐中文字幕| 伦精品一区二区三区| 午夜福利在线观看吧| 成人毛片a级毛片在线播放| 蜜桃亚洲精品一区二区三区| 丰满人妻一区二区三区视频av| 亚洲精华国产精华液的使用体验 | 最近中文字幕高清免费大全6 | 欧美黑人巨大hd| 国产高清视频在线播放一区| 成人亚洲精品av一区二区| 久久精品影院6| 久久久久久久久久成人| a在线观看视频网站| 成人国产一区最新在线观看| 日本在线视频免费播放| a级毛片a级免费在线| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| 日本 av在线| 欧美最黄视频在线播放免费| 亚洲aⅴ乱码一区二区在线播放| 国产一级毛片七仙女欲春2| 亚洲一区高清亚洲精品| 国产91精品成人一区二区三区| 在线播放无遮挡| 午夜老司机福利剧场| 精品人妻1区二区| 天堂影院成人在线观看| 国产精品女同一区二区软件 | 99热这里只有是精品在线观看| 99热网站在线观看| 看黄色毛片网站| 国模一区二区三区四区视频| 亚洲真实伦在线观看| 国产大屁股一区二区在线视频| 淫秽高清视频在线观看| 免费av不卡在线播放| 很黄的视频免费| 91麻豆精品激情在线观看国产| 联通29元200g的流量卡| 男女之事视频高清在线观看| 1000部很黄的大片| 亚洲avbb在线观看| a级毛片免费高清观看在线播放| 一个人看视频在线观看www免费| 在现免费观看毛片| 五月伊人婷婷丁香| 亚洲三级黄色毛片| 99久久无色码亚洲精品果冻| 日韩欧美国产一区二区入口| 一本精品99久久精品77| 在线观看美女被高潮喷水网站| 不卡一级毛片| 精品99又大又爽又粗少妇毛片 | 少妇高潮的动态图| 最后的刺客免费高清国语| 身体一侧抽搐| xxxwww97欧美| 精品午夜福利在线看| 亚洲无线观看免费| 免费搜索国产男女视频| 日本-黄色视频高清免费观看| 免费在线观看成人毛片| 搡女人真爽免费视频火全软件 | 在线观看66精品国产| 国产视频内射| 在线免费观看的www视频| 老司机深夜福利视频在线观看| 午夜精品久久久久久毛片777| 热99在线观看视频| 91在线观看av| 成年人黄色毛片网站| 成人三级黄色视频| 啦啦啦观看免费观看视频高清| 成人午夜高清在线视频| 婷婷精品国产亚洲av| 真人一进一出gif抽搐免费| av在线老鸭窝| 国产极品精品免费视频能看的| 亚洲av免费高清在线观看| 成人特级黄色片久久久久久久| 日韩亚洲欧美综合| 国产精品久久久久久亚洲av鲁大| 免费看av在线观看网站| 久久99热这里只有精品18| 国产激情偷乱视频一区二区| а√天堂www在线а√下载| 亚洲人成网站在线播放欧美日韩| 岛国在线免费视频观看| 日本-黄色视频高清免费观看| 国产真实乱freesex| 神马国产精品三级电影在线观看| 亚洲人成伊人成综合网2020| 国产精品久久久久久亚洲av鲁大| 久久精品国产99精品国产亚洲性色| 国产 一区精品| 高清在线国产一区| АⅤ资源中文在线天堂| 国产伦一二天堂av在线观看| 久久午夜福利片| 国产爱豆传媒在线观看| 精品99又大又爽又粗少妇毛片 | 97超级碰碰碰精品色视频在线观看| 亚洲精品成人久久久久久| 国产v大片淫在线免费观看| 一区福利在线观看| 久久久久久久亚洲中文字幕| 国产一区二区激情短视频| 精品福利观看|