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

    Highly efficient photocatalytic conversion of CO2 into CH4 over Cu single atom promoted heterojunction: The effect of uplifted d-band center

    2024-04-05 02:28:46HaoZhangQianLiuZhuruiShen
    Chinese Chemical Letters 2024年2期

    Hao Zhang,Qian Liu,Zhurui Shen

    School of Materials Science and Engineering,Nankai University,Tianjin 300350,China

    Keywords: Cu single atom Heterojunction CO2 photoreduction d-band center DFT

    ABSTRACT Photoreduction of CO2 to solar fuels has caused great interest,but suffers from low catalytic efficiency and poor selectivity.Herein,we designed a S-scheme heterojunction (Cu-TiO2/WO3) with Cu single atom to significantly boost the photoreduction of CO2.Notably,the developed Cu-TiO2/WO3 achieved the solardriven conversion of CO2 to CH4 with an evolution rate of 98.69 μmol g-1 h-1,and the electron selectivity of CH4 reached 88.5%.The yield was much higher than those of pristine WO3,TiO2/WO3 and Cu-TiO2 samples.Experimental and theoretical analysis suggested that the S-scheme heterojunction accelerated charge migration and inhibited the recombination of electron-hole pairs.Importantly,the charge separation effect of the heterojunction meliorated the position of the d-band.The uplifted d-band centers of Cu and Ti on Cu-TiO2/WO3 not only improved the electron interaction between Cu single atoms and substrate-TiO2,accelerated the adsorption and activation of CO2 on the active sites of Cu single atom,but also optimized the Gibbs free energies of CH4 formation pathway,leading to excellent selectivity toward CH4.This work provides new insights into the design of photocatalyst systems with high photocatalytic performance.

    Excessive emissions of carbon dioxide (CO2) caused by the everincreasing consumption of fossil fuels have resulted in many environmental problems [1-8].The conversion of photocatalytic CO2into value-added chemical products (CH4,CO,C2H4,etc.) has attracted significant interest because it can utilize solar light and H2O,that are abundant in nature,as energy and proton sources,respectively,without secondary pollution to the air environment.It is acknowledged that the prerequisite of improving the reduction performance of CO2photocatalysis is the establishment of suitable and high-performance photocatalysts.TiO2is one of the most widely utilized semiconductor materials in catalytic reactions owing to its environmental benignity,cost-effectiveness,and high physical and chemical stability [1,5,6,9].More importantly,the conduction band (CB) of TiO2is relatively negative,which would improve the reduction ability of photogenerated electrons.However,the limited charge separation efficiency and optical absorption hinder the reduction performance of the CO2photocatalysisviaTiO2.

    In the strategies of regulating activity of CO2photoreduction,constructing heterojunction with suitable semiconductors,photocatalysts with excellent performance can be obtained,which would facilitate charge separation by trapping photoelectrons [1,2,10-14].Tungsten oxide (WO3),has attracted much attention with relatively small bandgap (2.2-2.6 eV),low in toxicity as well as stable in acidic and oxidative conditions [10,14].Importantly,WO3has the characteristic of band matching with TiO2[11,13],coupling the two semiconductors to build a heterojunction,which improves the charge separation by trapping the photoelectrons and extending the energy range of photoexcitation toward the visible range,seems to be an efficient strategy to promote the performance of CO2photoreduction.Generally,various heterostructures have been studied for photoreduction of CO2[1,10-14].However,compare with conventional fuel production technologies,the effi-ciency of solar to fuel conversion is still very moderate.Thus,advanced strategies are required to design and fabricate a challenging heterojunction photocatalyst to improve efficiency of CO2conversion.

    Metal co-catalysts have also been widely applied in many photochemical reactions because of their high activity and selectivity [15-24].In particular,when tuned to the atomic scale,metal catalysts can maximize atomic utilization and accelerate the surface catalytic process;therefore,a number of atomically dispersed metals on semiconductors have emerged as effective platforms for establishing high-performance photocatalysts in recent years[4,6,8,19].Compared with noble metals (Au,Rh,and Pt),the earthabundant metal Cu is regarded as one of the most attractive metal co-catalysts for producing hydrocarbons from CO2due to the diversity of reduction products [4,8,20-22].For example,Lietal.[4] developed a photoinduction strategy to achieve the formation of Cu single atoms (SA) on a UiO-66-NH2support.Notably,the introduction of Cu SAs not only significantly enhanced transfer and separation efficiency,but also acted as active centers to facilitate the conversion of CO2to CHO*and CO*intermediates,leading to excellent selectivity toward methanol and ethanol.Likewise,wangetal.[7] designed single-atom Cu modified polymeric carbon nitride (CN),in which the C-Cu-N2single-atom catalytic site activated CO2molecules and reduced the energy barrier toward photocatalytic CO2reduction,resulting simultaneous production of CO,CH4and CH3OH.Inspired by these efforts,the improvement of CO2reduction activity on heterojunction could be achievedvialoading Cu SA on its surface.In fact,there are several reports concerning single atomic multi-heterojunctions on photoreduction of CO2.However,the mechanism of CO2conversion is not well elaborated,especially the effect of charge separation on the d-band of metal single atom.It is known that the interaction between metal d-band and CO2molecular orbitals is directly reflected in the adsorption energy,the position of d-band center also determined by the interact with CO2molecular orbitals of intermediates in catalytic behavior of photocatalyst [25-28].Therefore,the d-band center has proven to be an extremely useful electronic descriptor to explain photocatalytic CO2reduction selectivity.

    Hence,in this work,Cu-TiO2/WO3heterojunction was synthesized by a simple chemisorption strategy to investigate the activity and selectivity of CO2photocatalytic reaction.In the process,Cu single atoms were firstly loaded on the surface of TiO2(Cu-TiO2)viareduction method,and then coupled with WO3to synthesize a S-scheme heterojunction photocatalyst to ensure efficient electron utilization and chemical transformation.Under simulated solar light irradiation,the Cu-TiO2/WO3photocatalyst displayed an obviously enhanced photoreduction activity and selectivity toward CH4.Using first-principles density functional theory (DFT) calculations,we revealed the decisive influence of charge separation effect of heterojunction on interaction between Cu single atom and TiO2substrate,and the d-band position.

    The morphologies and microstructures of as-prepared samples were shown in Fig.1 and Fig.S1 (Supporting information).In the high-resolution transmission electron microscope (HRTEM) images (Fig.1a),bare TiO2was uniform and the lattice fringes of 0.33 nm correspond to the (101) planes of anatase-TiO2[1,6].And in the aberration-corrected high-angle annular dark-field scanning transmission microscopy (AC HAADF-STEM) images,isolated bright spots were observed on Cu-TiO2sample,which indicated that Cu were homogeneously diffused on the surface of TiO2nanoparticles in the form of atoms rather than clusters.After Cu-TiO2was contacted with WO3,a clear boundary line (red) was observed (Fig.1c),and the lattice fringes with distances of 0.33 nm and 0.37 nm were attributed to the (101) crystal plane of TiO2and (200) crystal plane of WO3,respectively [1,6].The result indicated that the Cu-TiO2/WO3heterojunction was successfully obtained.Thereafter,according to the energy-dispersive X-ray (EDX) mapping shown in Figs.1d and e,Cu-TiO2/WO3displayed homogeneous distribution of Ti,O,Cu and W elements.Besides,according to the result of inductive coupled plasma emission spectra (ICP) shown in Table S1(Supporting information),the loading amount of Cu single atom in Cu-TiO2/WO3was just slightly loss,which may be caused by cleaning during in the process of preparation.

    Fig.1.Images of (HR)TEM over TiO2 (a),AC HAADF-STEM over Cu-TiO2 (b),HRTEM over Cu-TiO2/WO3 (c).EDX elemental mapping over Cu-TiO2 (d) and Cu-TiO2/WO3(e).

    To explore the coordination environment and chemical states of Cu atoms,X-ray absorption fine structure (XAFS) was measured on Cu-TiO2/WO3.As shown in Fig.2a,the peak of Cu K-edge Xray absorption near edge structure (XANES) located at 8997 eV of Cu-TiO2/WO3was similar to that of CuO,indicating the existence of Cu single atom with a positive 2 valence.Subsequently,k3-weighted Cu K-edge Fourier transform (FT) curve and wavelet transform (WT) analysis of Cu-TiO2/WO3were shown in Fig.2b and Fig.S2 (Supporting information),except for a main scattering peak at 1.47 ?A attributed to the Cu-O coordination,the peak at 2.24 ?A corresponding to Cu-Cu interaction in Cu-TiO2/WO3was not detected,confirming the single atom dispersion of Cu.According to the fitting results of XANES shown in Table S2 (Supporting information),the Cu single atoms were coordinated with 4.0 ± 0.26 O atoms of TiO2with an average bond length of 1.922 ± 0.006 ?A,which was slight shorter than Ti-O bond (1.94 ?A) due to the lattice distortion and/or contraction in Cu-TiO2/WO3.Consequently,the tetrahedral oxygen-coordinated Cu-O4sites were anchored at the Ti sites,schematically shown in Fig.S3 (Supporting information).

    Fig.2.Cu K-edge XANES (a) and Fourier transform of k3-weighted Cu K-edge EXAFS(b) of Cu-TiO2/WO3 and the reference Cu foil,Cu2O,and CuO.XRD patterns (c) of samples and Cu 2p of XPS (d) over Cu-TiO2/WO3.

    Fig.2c displayed the X-ray diffraction (XRD) patterns of samples.The spectra of TiO2agreed with simulated anatase structure (JCPDS No.21-1272),and the introduction of Cu SA did not change TiO2crystallinity.For Cu-TiO2/WO3heterojunction,besides the peaks of TiO2,the peaks assigned to WO3(yellow line) were observed [1,14],indicating the successful construction of heterojunction.Next,as shown in Fig.S4a (Supporting information),Xray photoelectron spectrum (XPS) of Cu-TiO2/WO3verified the existence of Ti,O,Cu,and W.In Fig.S4b (Supporting information),the Ti 2p spectra of TiO2exhibited Ti 2p1/2and Ti 2p3/2,respectively.A positive (0.21 eV) shift was observed in Cu-TiO2,indicating a decreasing of electron density on Ti owing to the electron enrichment of Cu sites.For Cu-TiO2/WO3,a larger shift was happened (0.45 eV),suggesting that more electrons on Ti were transferred to Cu sites in the heterojunction.In Fig.2d,the spectra of Cu 2p was divided into 2p1/2and 2p3/2on Cu-TiO2[8,21],respectively.While,the both peaks shifted to the lower binding energies (0.28 eV) in Cu-TiO2/WO3,indicating the electron density of Cu sites was increased.A similar tendency with Ti 2p appeared in W 3d XPS shown in Fig.S4c (Supporting information),which indicated that the electrons of W site of heterojunction also decreased.The results suggested that Cu-TiO2/WO3heterojunction accelerated the charge migration behavior.

    Ultraviolet-visible diffuse reflectance spectroscopy (UV-vis DRS)was applied to survey the optical property of samples.In Fig.S5a (Supporting information),comparing with Cu-TiO2,the range of optical absorption over Cu-TiO2/WO3further broaden due to the visible adsorption of WO3,which was beneficial to CO2photoreduction.To reveal band structure of the heterojunction,the bandgap and Mott-Schottky plots were conducted and shown in Figs.S5b-e and Table S3 (Supporting information).According to the results,the valence bands (VB) of TiO2and WO3were 2.44,and 2.53 V,their conduct bands (CB) were -0.61 and 0.02 V [11,13],respectively.The movement of electrons in heterojunction photocatalyst can be described by the work functions (W) of semiconductors.W of TiO2and WO3were 4.69 and 5.05 eV [11,13].When Cu-TiO2/WO3heterojunction was formed,the free electrons would transfer from TiO2to WO3easily until the Fermi levels (Ef) tend to equilibrate at the interface.Concurrently,an internal electric field(IEF) was formed.When the heterojunction was exposed to light,photoelectrons were excited and transferred from VB to CB of WO3and TiO2,respectively.Under the effect of IEF,the photoelectrons on CB of WO3easily migrated and combined with the holes on VB of TiO2.This charge transfer behavior of Cu-TiO2/WO3heterojunction was matched with the S-scheme mechanism [11,13],which was consistent with the XPS results.The illustration of charge transfer under simulated sunlight for efficient photocatalysis over Cu-TiO2/WO3was shown in Fig.S6 (Supporting information).In this regard,efficient charge separation would be achieved.Finally,more photoinduced electrons would participate in the photocatalysis process.Comparing with E0(CO2/CH4=-0.17 Vvs.NHE,pH 0),Cu-TiO2/WO3has the ability of photocatalytic reduction of CO2to CH4.

    To evaluate photocatalytic property of Cu-TiO2/WO3,CO2photoreduction experiment was performed under simulated sunlight irradiation with water as hole sacrificial reagents.To optimize the photoreduction performance,TiO2loaded with different amounts of Cu single atom and the heterojunction with different mass ratios of Cu-TiO2and WO3were also prepared.CO2photoreduction experiments of all samples were shown in Fig.3a and Fig.S7 (Supporting information),both of pure TiO2and WO3converted CO2into CO,with formation rates of 3.28 and 1.08 μmol g-1h-1,respectively.After Cu SA introduced,CH4become the main product of CO2photoreduction,and a small amount of CO and C2H4were also observed.Moreover,comparing to 0.5 and 2.0 mL,1.0 mL of Cu2+solution was the optimal volume and Cu-TiO2exhibited the supreme photocatalytic performance with formation rates of 32.29 μmol g-1h-1for CH4.The results suggested that the Cu single atom center may be a highly active site for CO2to CH4.In addition,too small amount of Cu single atom would lead to insufficient active sites for photoreduction,and excessive ones would result in aggregation,which also affected the activity of CO2photoreduction.For Cu-TiO2/WO3photocatalyst,it found that CH4was still the main product in CO2conversion.And the evolution rate of CH4reached 98.69 μmol g-1h-1,which was much higher than those of pristine WO3,TiO2/WO3and Cu-TiO2samples.The results indicated that the formation rate of CH4were significantly improved due to the charge separation effect of Cu-TiO2/WO3heterojunction.

    Fig.3.Products from CO2 photoreduction reaction (a),labeling test of 13CO2 on Cu-TiO2/WO3 (b),in-situ FTIR spectroscopy measurement over Cu-TiO2/WO3 (c).

    To further verify the role played by Cu single atom and the effect of charge migration on Cu-TiO2/WO3heterojunction,TiO2/WO3sample was prepared in the absence of Cu single atom.In Fig.3a,it found that TiO2/WO3only could convert CO2into CO,with a rate of 6.13 μmol g-1h-1,and no CH4was detected,suggesting that Cu site may be the key active center for the conversion of CO2to CH4on Cu-TiO2/WO3.Next,Cu-TiO2and WO3were mixed by simple mechanical stirring without interfacial chemical contact.This sample was named as Cu-TiO2+WO3and its photoreduction result displayed in Fig.S7b (Supporting information) was just close to that of Cu-TiO2,far lower than that of Cu-TiO2/WO3,indicating that the construction of Cu-TiO2/WO3significantly improved the activity of Cu active sites due to the charge migration and separation,thereby accelerating the conversion of CO2to CH4.In Table S4 (Supporting information),the selectivity of CH4was also calculated.For Cu-TiO2/WO3,the electron selectivity toward CH4was reached 88.5%,which were 15.1% higher than that of Cu-TiO2,indicating that the charge separation effect of the heterojunction promoted the pathway reaction of CO2to CH4.The above results indicated that the construction of Cu-TiO2/WO3photocatalyst not only enhanced the activity of CO2photoreduction,but also increased the selectivity toward CH4.

    To further confirm the photoreduction of CO2,13C isotope labeling and blank control tests were performed.As shown in Fig.3b,the overwhelming peaks ofm/z=17 (13CH4),29 (29CO) and 30(13C2H4) rooted in13CO2sources indicated that all reduction products were originated from the13CO2rather than other carbonaceous impurities.Blank control experiments were shown in Fig.S8 (Supporting information),where no reduction product was detected in N2or in the absence of photocatalysts,light irradiation,and H2O,indicating that CO2gas,photocatalyst,light irradiation,and H2O were the prerequisites for the photoreduction of CO2.Subsequently,the stability of the heterojunction was investigated.As shown in Fig.S9 (Supporting information),the XRD structure of Cu-TiO2/WO3did not change before and after photoreduction reaction.The result proved that Cu-TiO2/WO3has a relatively stable structure.

    Next,we studied charge-transfer properties to investigate the reason of the enhanced photocatalytic activity on Cu-TiO2/WO3.In Fig.S10a (Supporting information),the transient photocurrent of heterojunction was enhanced significantly compared with those of other samples,implying that more effective charge separation over Cu-TiO2/WO3.And in Fig.S10b (Supporting information) Cu-TiO2/WO3showed a marginally lowest PL signal,indicating more effective retardation of the charge recombination [14].Electrochemical impedance spectroscopy (EIS) can also reflects the charge-transfer behavior.As shown in Fig.S10c (Supporting information),Cu-TiO2/WO3has the lowest charge-transfer resistance of all samples,suggesting that more efficient transmission of charge have been achieved in Cu-TiO2/WO3.According to the results,it demonstrated that the S-scheme heterojunction could promote the migration and separation of surface photogenerated charges,thus improving photocatalytic performance of CO2,which provided a clear evidence for electron transfer within several components of Cu-TiO2/WO3.Time-resolved photoluminescence (TRPL) profiles of samples were measured to further reveal the efficient separation of photogenerated carriers.At the excitation wavelengths (EW) of 425 nm,the fitted decay curves disclose the lifetime (τ) and percentage (Rel.%) of charge carriers.As shown in Fig.S10d and Table S5 (Supporting information),Cu-TiO2/WO3has a lowest percentage ofτ1-carriers than other samples,indicating that the charge recombination was further inhibited due to the charge migration[14].The average carrier lifetime of Cu-TiO2/WO3(τa=8.86 ns)was higher than that Cu-TiO2(4.92 ns),which was attributed to efficient transfer and separation of charges over Cu-TiO2/WO3heterojunction.Therefore,it was not surprising that the Cu-TiO2/WO3composite sample exhibit improved photocatalytic CO2reduction performance.

    To understand the reaction mechanism of CO2photoreduction over Cu-TiO2/WO3,insituFTIR measurements were performed to investigate the intermediates in the process of CO2conversion.As shown in Fig.3c,the peaks located at 1075,1145 and 1341 cm-1were detected and ascribed to the CHO*,CH3O*,and COOH*intermediates [29-32],respectively,which were the crucial groups for CO2conversion to CH4.A new peak at 1203 cm-1was observed and ascribed to COCO*group [33],a crucial C-C coupling intermediate of C2H4formation.Additionally,the peak at 1659 cm-1corresponded to the C=O stretching belonging to the CO*intermediate [34].TheinsituFTIR was agreed well with the results of CO2photoreduction test on Cu-TiO2/WO3.Based on the results,it speculated that Cu sites converted the adsorbed CO2to COOH*,CO*,CHO*and CH3O*intermediates with the participation of protons and electrons,which in turn to produce CH4.For weaker adsorbed groups [34],such as CO*,could be desorbed from the Cu active sites to generate CO gas.Meanwhile,partial CO*would also couple with the neighbor CO*to form COCO*intermediates,and finally form C2H4with the assistance of multiple electrons and protons.

    Next,DFT calculation was applied to examined the mechanism of CO2photoreduction on Cu-TiO2/WO3photocatalyst.The model structure of the heterojunction was built by TiO2(101) and WO3(200) planes by replacing the Ti with Cu atoms.To further confirm the active site of the heterojunction,the adsorption energies of CO2on different metal sites were calculated and shown in Table S6 (Supporting information).According to the results,it proved that the Cu single atom site was the optimal active center in Cu-TiO2and Cu-TiO2/WO3photocatalysts for CO2to CH4.Next,Gibbs free energies (ΔG) were calculated in the premise of Cu single atom as active site,and shown in Figs.4a and b and Table S7 (Supporting information).TheΔGof photoreduction of CO2to COOH*,CO*,CHO*,CH3O*and CH4formation were displayed in Fig.4a.Along with the hydrogen proton (H+) produced by water hydrolysis,the C bond of CO2*was interacted with a pair of proton and electron to form COOH*.According to the calculated results,it found that the whole process of CO2-to-CH4conversion on the Cu-TiO2/WO3was more favorable in thermodynamics behavior than that on the Cu-TiO2surface.And the formation of COOH*intermediate was the rate-limiting step for both photocatalysts.For Cu-TiO2/WO3,an energy of 0.54 eV was needed to complete the step of CO2→COOH*.By contrast,ΔGof 0.79 eV (ΔGCOOH* -ΔGCO2)was essential in endothermic reaction of Cu-TiO2,suggesting the structure of Cu-TiO2/WO3could efficiently lower the formation energy of key COOH*intermediate.From the results,it suggested that more electrons were transferred to Cu active sites from TiO2due to the charge separation effect of the heterojunction,which helped to optimize the rate-limiting step and promote the conversion of CO2-to-CH4.

    Fig.4.Free-energy diagrams of CO2 photoreduction over Cu-TiO2 and Cu-TiO2/WO3 calculated from DFT (a),structure models of CO2 to CH4 pathway over Cu-TiO2/WO3(b),the partial density of states (PDOS) of Cu-TiO2 (c) and Cu-TiO2/WO3 (d) and the dashed lines represent the positions of d-band centers.

    Typically,the adsorption and activation of CO2is also critical for photocatalytic reduction [29].Thus,the molecule structure of CO2before and after adsorption on the catalyst surface was calculated.In Table S8 (Supporting information),comparing with Cu-TiO2,the CO2adsorbed on Cu-TiO2/WO3has a smaller C-O-C angle,and longer C=O bonds,resulting in more curved configuration of adsorbed CO2and easier cleavage of C=O bonds [35,36].Furthermore,in Table S9 (Supporting information),the adsorption energy of CO2on heterojunction was lower than Cu-TiO2,indicating CO2adsorption and activation were easier to occurred on its surface.

    Electronic properties of the catalyst were further carried out and analyzed to figure out the underlying factors governing the significant activity of the heterojunction,Firstly,the feature of CO2adsorption was further evaluated through the partial density of states (PDOS) on Cu active site,and the results were shown in Fig.S11 (Supporting information).When the O atom of CO2molecule was adsorbed on the Cu-TiO2/WO3surface,the peaks shifted to lower energy than that of Cu-TiO2,suggesting that the orbital coupling of CO2molecule with Cu active site of Cu-TiO2/WO3was more prone to carry out [37].Next,the d-band center of Cu site was further explored to elaborate the adsorption behavior.In Fig.4c and Table S10 (Supporting information),it suggested that the d-band centers of Cu in two catalysts depended on its local coordination environment,and was calculated to be -1.49 (Cu-TiO2/WO3) and -1.77 eV (Cu-TiO2).Further analysis found that a more remarkable overlap between Cu 3d in Cu-TiO2/WO3and O 2p of CO2at the identical position.Apparently,more substantial orbital and electron interactions of Cu sites and O were present than that of Cu in Cu-TiO2,which uplifted the d-band center of Cu-TiO2/WO3closer to the Fermi level.Such changes in the d-band centers of Cu were precisely correlated with the adsorption energy of key COOH*intermediates.In addition,the d-band of Ti was also calculated to elaborate the electron interaction between the Cu SA and the substrate TiO2in the heterojunction.As shown in Fig.4d,comparing with Cu-TiO2,the d-band of Ti become narrowed in Cu-TiO2/WO3,and the center of the d-band was also uplifted and closer to the fermi level,indicating more stronger interaction between the Cu SA sites and TiO2substrate.

    The S-scheme heterojunction of coupling TiO2supported single Cu atoms with WO3provided an efficient method for charge separation and migration.The results showed that Cu single atom played a key role as active centers in CO2photoreduction over Cu-TiO2/WO3photocatalyst,and the heterojunction significantly improved the mobility efficiency of surface photogenerated charge based on the unique advantages of S-scheme.Moreover,the charge separation effect of the heterojunction further promoted the adsorption and activation of CO2on Cu active sites.Importantly,the uplifted d-band centers of Cu and Ti of Cu-TiO2/WO3photocatalyst enhanced the electron interaction between Cu single atom and TiO2substrate,and optimized the Gibbs free energies of CH4formation.The specific pattern of combining single-atom and heterojunction fully exploits the highly efficient charge separation to open up a new channel for the design of efficient photocatalyst systems.

    Declaration of competing interest

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

    Acknowledgments

    This work was supported by the grants from the National Natural Science Foundation of China (Nos.21872102 and 22172080) and the Tianjin “Project+Team” innovation team,2020.

    Supplementary materials

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

    久久天躁狠狠躁夜夜2o2o| 久久久国产精品麻豆| 久久久久久九九精品二区国产| 精品国内亚洲2022精品成人| 亚洲男人的天堂狠狠| 久久久久九九精品影院| 九色成人免费人妻av| 亚洲欧美日韩无卡精品| 少妇熟女aⅴ在线视频| 伦理电影免费视频| 亚洲欧美精品综合久久99| xxxwww97欧美| 青草久久国产| 黄色视频,在线免费观看| 亚洲精品美女久久久久99蜜臀| 亚洲成人中文字幕在线播放| 国产亚洲av高清不卡| 一个人看的www免费观看视频| 亚洲自偷自拍图片 自拍| 国产av一区在线观看免费| 午夜精品久久久久久毛片777| 欧美性猛交黑人性爽| 精品国产三级普通话版| 精品不卡国产一区二区三区| 国产成人av激情在线播放| 美女高潮的动态| 亚洲精品久久国产高清桃花| 男女床上黄色一级片免费看| av天堂在线播放| av福利片在线观看| 国产aⅴ精品一区二区三区波| 国产精品美女特级片免费视频播放器 | 午夜福利在线在线| 欧美不卡视频在线免费观看| 亚洲午夜理论影院| 黄色丝袜av网址大全| 久久精品国产99精品国产亚洲性色| 欧美成狂野欧美在线观看| 两个人视频免费观看高清| 午夜亚洲福利在线播放| 精品乱码久久久久久99久播| 亚洲欧美日韩无卡精品| 久久香蕉精品热| 成人高潮视频无遮挡免费网站| 精品久久久久久久久久久久久| 欧美日韩一级在线毛片| 亚洲精品中文字幕一二三四区| 日韩欧美三级三区| 久久精品91蜜桃| 一级毛片女人18水好多| 一进一出抽搐动态| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久毛片微露脸| 久久精品综合一区二区三区| 精品久久久久久,| 久久久久亚洲av毛片大全| 亚洲av成人不卡在线观看播放网| 亚洲色图av天堂| 亚洲中文日韩欧美视频| 亚洲欧美日韩高清在线视频| 国产激情久久老熟女| 欧美色欧美亚洲另类二区| 久久久久国内视频| 性色av乱码一区二区三区2| 天堂动漫精品| 少妇丰满av| 伊人久久大香线蕉亚洲五| 亚洲 欧美 日韩 在线 免费| 搞女人的毛片| 看片在线看免费视频| 高清在线国产一区| 欧美中文日本在线观看视频| 亚洲国产色片| 久久精品91无色码中文字幕| 久久久成人免费电影| 伊人久久大香线蕉亚洲五| 亚洲熟女毛片儿| 日韩精品中文字幕看吧| 9191精品国产免费久久| 国产 一区 欧美 日韩| 久久草成人影院| 在线观看免费午夜福利视频| 性色avwww在线观看| 国产成人精品无人区| or卡值多少钱| 露出奶头的视频| 精品国产乱码久久久久久男人| 亚洲精品国产精品久久久不卡| 亚洲精品粉嫩美女一区| 久久久久久久久久黄片| 黄色女人牲交| 亚洲欧美日韩无卡精品| av中文乱码字幕在线| 日韩免费av在线播放| av天堂中文字幕网| 欧美乱色亚洲激情| 久久久久久久久中文| 免费观看精品视频网站| 成人一区二区视频在线观看| 制服人妻中文乱码| 99久久成人亚洲精品观看| 亚洲av电影不卡..在线观看| 亚洲欧美激情综合另类| 精品一区二区三区视频在线 | 欧美黑人欧美精品刺激| 午夜a级毛片| 国产一区在线观看成人免费| 国产高清激情床上av| 一二三四社区在线视频社区8| 亚洲狠狠婷婷综合久久图片| 免费看十八禁软件| 成年女人毛片免费观看观看9| 亚洲人成网站高清观看| 99热精品在线国产| 精品国内亚洲2022精品成人| 亚洲午夜精品一区,二区,三区| 久久久久国产精品人妻aⅴ院| 天天一区二区日本电影三级| 一个人看视频在线观看www免费 | 精品一区二区三区av网在线观看| 久久伊人香网站| 日本黄色片子视频| 亚洲人成伊人成综合网2020| 国产日本99.免费观看| 国产日本99.免费观看| 在线免费观看不下载黄p国产 | 欧美乱妇无乱码| 国产成人影院久久av| 久久性视频一级片| 亚洲欧美日韩高清专用| 男人和女人高潮做爰伦理| 亚洲成人久久爱视频| 级片在线观看| 欧美日韩亚洲国产一区二区在线观看| 又爽又黄无遮挡网站| 国产视频内射| 最新在线观看一区二区三区| 久久久久性生活片| 美女 人体艺术 gogo| 国产成人av教育| 精品熟女少妇八av免费久了| 女警被强在线播放| 久99久视频精品免费| 在线观看日韩欧美| 好男人在线观看高清免费视频| 毛片女人毛片| 欧美日韩福利视频一区二区| 麻豆成人午夜福利视频| 麻豆成人午夜福利视频| 国产美女午夜福利| 亚洲成a人片在线一区二区| 国产69精品久久久久777片 | 成人特级黄色片久久久久久久| 精品一区二区三区视频在线观看免费| 亚洲精品在线美女| 亚洲av成人不卡在线观看播放网| 亚洲av成人不卡在线观看播放网| 久久精品夜夜夜夜夜久久蜜豆| 少妇丰满av| 亚洲av免费在线观看| 91av网站免费观看| 亚洲国产精品成人综合色| 日日干狠狠操夜夜爽| 欧美激情在线99| 精品久久久久久久末码| 国产成人影院久久av| 男人的好看免费观看在线视频| 黄色成人免费大全| 99国产精品一区二区三区| 搞女人的毛片| 老司机在亚洲福利影院| 国产精品亚洲一级av第二区| 精品一区二区三区视频在线观看免费| 久久久久久人人人人人| 亚洲真实伦在线观看| 国产野战对白在线观看| 亚洲真实伦在线观看| 国产不卡一卡二| 国产伦精品一区二区三区视频9 | 欧美日韩瑟瑟在线播放| 757午夜福利合集在线观看| 夜夜看夜夜爽夜夜摸| 1000部很黄的大片| 欧美av亚洲av综合av国产av| 九九久久精品国产亚洲av麻豆 | 极品教师在线免费播放| 性色av乱码一区二区三区2| 精品久久久久久久末码| 精品无人区乱码1区二区| 校园春色视频在线观看| 国产精品一区二区免费欧美| 99精品在免费线老司机午夜| 搡老熟女国产l中国老女人| 婷婷亚洲欧美| 国产精品综合久久久久久久免费| 国产一区二区在线观看日韩 | 国产99白浆流出| 两个人的视频大全免费| 村上凉子中文字幕在线| 麻豆国产97在线/欧美| 黄片大片在线免费观看| 美女高潮的动态| 天堂av国产一区二区熟女人妻| 一个人看的www免费观看视频| 国产午夜福利久久久久久| 成人国产一区最新在线观看| 97人妻精品一区二区三区麻豆| 1024香蕉在线观看| 日韩三级视频一区二区三区| 国产精品久久久久久人妻精品电影| 国产精华一区二区三区| 精品久久蜜臀av无| 99热这里只有是精品50| 国产成+人综合+亚洲专区| 国产亚洲av高清不卡| 欧美日本亚洲视频在线播放| 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 欧美色视频一区免费| 国产精品永久免费网站| 好看av亚洲va欧美ⅴa在| 亚洲 国产 在线| 亚洲av成人av| 国产成人aa在线观看| 99精品在免费线老司机午夜| 亚洲av成人不卡在线观看播放网| 俄罗斯特黄特色一大片| 国产v大片淫在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品,欧美在线| 亚洲国产精品成人综合色| 99riav亚洲国产免费| 免费看日本二区| 久久久久久久久免费视频了| 亚洲无线观看免费| 国产亚洲精品久久久久久毛片| 国产乱人视频| 午夜福利欧美成人| 亚洲精品美女久久av网站| 搡老妇女老女人老熟妇| 国产精华一区二区三区| 中文字幕精品亚洲无线码一区| 成人亚洲精品av一区二区| 日本一本二区三区精品| 国产亚洲精品一区二区www| 日本成人三级电影网站| www.自偷自拍.com| 精华霜和精华液先用哪个| 在线免费观看的www视频| 十八禁网站免费在线| 男女之事视频高清在线观看| 亚洲成av人片免费观看| 久久久国产成人精品二区| 淫妇啪啪啪对白视频| 亚洲午夜理论影院| 曰老女人黄片| 亚洲一区二区三区色噜噜| 亚洲自拍偷在线| 国产黄a三级三级三级人| 国产成人精品无人区| 好男人在线观看高清免费视频| 成人三级做爰电影| 波多野结衣巨乳人妻| 白带黄色成豆腐渣| 日本五十路高清| 免费av毛片视频| 99国产综合亚洲精品| 国产精品久久视频播放| 一进一出抽搐动态| 欧美大码av| 99国产精品一区二区三区| 99久久久亚洲精品蜜臀av| 国产主播在线观看一区二区| 变态另类丝袜制服| 国产精品 国内视频| 99久久精品国产亚洲精品| 久久国产精品人妻蜜桃| 香蕉av资源在线| 亚洲性夜色夜夜综合| 久久九九热精品免费| 亚洲国产中文字幕在线视频| 99久久99久久久精品蜜桃| 国产一级毛片七仙女欲春2| 国产一区在线观看成人免费| 欧美性猛交黑人性爽| 亚洲精品久久国产高清桃花| 成人av在线播放网站| 欧美成人免费av一区二区三区| 国产欧美日韩精品一区二区| 欧美日韩瑟瑟在线播放| 此物有八面人人有两片| 宅男免费午夜| 成人午夜高清在线视频| 国产麻豆成人av免费视频| 亚洲在线自拍视频| 特级一级黄色大片| 婷婷精品国产亚洲av在线| 国产精品一及| 日韩av在线大香蕉| 国产精品亚洲av一区麻豆| 免费一级毛片在线播放高清视频| 国产亚洲av高清不卡| 午夜福利欧美成人| 极品教师在线免费播放| 18禁黄网站禁片午夜丰满| 最近最新中文字幕大全电影3| 国产精品99久久99久久久不卡| 亚洲成人中文字幕在线播放| 国产精品爽爽va在线观看网站| 午夜精品在线福利| 天堂影院成人在线观看| 免费无遮挡裸体视频| 亚洲成人精品中文字幕电影| 久久久久久国产a免费观看| tocl精华| 最新美女视频免费是黄的| 久久国产精品影院| 欧美在线黄色| 国内毛片毛片毛片毛片毛片| 国产美女午夜福利| 在线十欧美十亚洲十日本专区| 少妇丰满av| 久久久久九九精品影院| 国产av一区在线观看免费| 亚洲最大成人中文| 97碰自拍视频| 亚洲熟妇熟女久久| 久9热在线精品视频| 国产aⅴ精品一区二区三区波| 国产精品99久久99久久久不卡| 美女大奶头视频| 男女之事视频高清在线观看| aaaaa片日本免费| 一本综合久久免费| 99国产极品粉嫩在线观看| 成年女人永久免费观看视频| 精品熟女少妇八av免费久了| 国内精品久久久久精免费| 国产99白浆流出| 亚洲,欧美精品.| a在线观看视频网站| 国产人伦9x9x在线观看| 一级黄色大片毛片| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 中文字幕最新亚洲高清| 日本成人三级电影网站| 久久精品国产亚洲av香蕉五月| 男女视频在线观看网站免费| 欧美日韩亚洲国产一区二区在线观看| 日韩国内少妇激情av| 久久国产乱子伦精品免费另类| 久久精品综合一区二区三区| 亚洲人与动物交配视频| 99久久精品一区二区三区| 国产乱人视频| 最近视频中文字幕2019在线8| 久久这里只有精品19| 少妇的逼水好多| www.自偷自拍.com| 九九热线精品视视频播放| 99热6这里只有精品| 91九色精品人成在线观看| 国产av在哪里看| 国产精品女同一区二区软件 | 免费搜索国产男女视频| 精品久久久久久成人av| 亚洲国产精品合色在线| 美女免费视频网站| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| 天天一区二区日本电影三级| 一级a爱片免费观看的视频| 免费观看人在逋| 精品日产1卡2卡| 亚洲午夜精品一区,二区,三区| av欧美777| 在线观看午夜福利视频| 中文在线观看免费www的网站| 欧美日韩黄片免| 国产精品九九99| 桃色一区二区三区在线观看| 欧美日本视频| 午夜福利18| 69av精品久久久久久| 熟女少妇亚洲综合色aaa.| 国语自产精品视频在线第100页| 婷婷六月久久综合丁香| 国产亚洲av嫩草精品影院| 人妻丰满熟妇av一区二区三区| 嫩草影视91久久| 叶爱在线成人免费视频播放| 欧美成人免费av一区二区三区| 亚洲色图av天堂| 99热只有精品国产| 欧美成狂野欧美在线观看| 久久久久久久久免费视频了| 黄色 视频免费看| 亚洲av电影在线进入| 久久午夜综合久久蜜桃| 国产一区二区在线观看日韩 | 久久国产乱子伦精品免费另类| 在线观看美女被高潮喷水网站 | 国产真人三级小视频在线观看| 久久亚洲真实| 中文字幕最新亚洲高清| 999久久久精品免费观看国产| 久久久水蜜桃国产精品网| 亚洲欧美激情综合另类| 免费在线观看视频国产中文字幕亚洲| 国内揄拍国产精品人妻在线| 久久性视频一级片| 国产av麻豆久久久久久久| 精品人妻1区二区| 亚洲av第一区精品v没综合| 亚洲精品乱码久久久v下载方式 | 国内精品一区二区在线观看| 黄色女人牲交| 欧美色视频一区免费| 国产伦在线观看视频一区| 亚洲片人在线观看| 亚洲成av人片免费观看| 午夜日韩欧美国产| 美女免费视频网站| 18美女黄网站色大片免费观看| 久久热在线av| 国产精品久久久久久亚洲av鲁大| 亚洲欧美一区二区三区黑人| 99热这里只有精品一区 | 国产精品香港三级国产av潘金莲| 美女cb高潮喷水在线观看 | 久久中文看片网| 国产美女午夜福利| 国产精品乱码一区二三区的特点| 精品久久久久久久久久免费视频| 老熟妇仑乱视频hdxx| 中文字幕熟女人妻在线| 校园春色视频在线观看| 成人av一区二区三区在线看| 偷拍熟女少妇极品色| 99精品在免费线老司机午夜| 一区福利在线观看| 在线观看日韩欧美| 成人高潮视频无遮挡免费网站| 1024香蕉在线观看| 亚洲人与动物交配视频| 亚洲国产精品合色在线| 在线免费观看的www视频| 熟女人妻精品中文字幕| 亚洲最大成人中文| 国产高清视频在线观看网站| 欧美黄色淫秽网站| 日本黄色视频三级网站网址| av在线天堂中文字幕| 99久久成人亚洲精品观看| 日韩免费av在线播放| 岛国在线免费视频观看| 一级毛片女人18水好多| 18禁黄网站禁片免费观看直播| 一本久久中文字幕| 在线观看舔阴道视频| 国产主播在线观看一区二区| 99精品久久久久人妻精品| 757午夜福利合集在线观看| 十八禁网站免费在线| 亚洲一区高清亚洲精品| 亚洲国产高清在线一区二区三| 国产精品自产拍在线观看55亚洲| 99久久成人亚洲精品观看| 伊人久久大香线蕉亚洲五| 欧美zozozo另类| 99精品在免费线老司机午夜| 一a级毛片在线观看| 国产精品美女特级片免费视频播放器 | 久久精品亚洲精品国产色婷小说| 亚洲激情在线av| 国产日本99.免费观看| 国产亚洲精品久久久久久毛片| 国产精品自产拍在线观看55亚洲| 丰满人妻熟妇乱又伦精品不卡| 日本五十路高清| 高清毛片免费观看视频网站| 欧美激情久久久久久爽电影| 床上黄色一级片| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| 最近最新中文字幕大全免费视频| 男女那种视频在线观看| 无遮挡黄片免费观看| 国产亚洲av嫩草精品影院| 日韩国内少妇激情av| 国产高清有码在线观看视频| 特级一级黄色大片| 天天添夜夜摸| 午夜日韩欧美国产| 亚洲性夜色夜夜综合| 国产伦在线观看视频一区| 久久久国产精品麻豆| 日本在线视频免费播放| 女人高潮潮喷娇喘18禁视频| 这个男人来自地球电影免费观看| 又黄又粗又硬又大视频| 亚洲国产欧洲综合997久久,| 亚洲男人的天堂狠狠| 久久欧美精品欧美久久欧美| 亚洲五月婷婷丁香| 岛国在线观看网站| 90打野战视频偷拍视频| xxx96com| 搡老岳熟女国产| 在线免费观看不下载黄p国产 | 白带黄色成豆腐渣| 久久久久久国产a免费观看| 日韩欧美国产在线观看| 国产精品久久久av美女十八| 精品熟女少妇八av免费久了| 精品不卡国产一区二区三区| 51午夜福利影视在线观看| 看片在线看免费视频| 亚洲av美国av| 欧美中文综合在线视频| 久久久国产欧美日韩av| 无人区码免费观看不卡| 全区人妻精品视频| 啪啪无遮挡十八禁网站| netflix在线观看网站| 久久精品影院6| svipshipincom国产片| 村上凉子中文字幕在线| 欧美日韩一级在线毛片| 黄色 视频免费看| 麻豆成人午夜福利视频| 亚洲一区二区三区色噜噜| 18禁裸乳无遮挡免费网站照片| 老熟妇乱子伦视频在线观看| 日日干狠狠操夜夜爽| 搡老熟女国产l中国老女人| 亚洲av日韩精品久久久久久密| 男女那种视频在线观看| 久久精品aⅴ一区二区三区四区| 老司机福利观看| 国产不卡一卡二| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 黄色丝袜av网址大全| 国产主播在线观看一区二区| 欧美另类亚洲清纯唯美| 亚洲成人精品中文字幕电影| 午夜激情欧美在线| 不卡一级毛片| 又大又爽又粗| 婷婷精品国产亚洲av| 欧美zozozo另类| www日本在线高清视频| 日韩欧美免费精品| 亚洲av片天天在线观看| 国产三级黄色录像| 免费在线观看日本一区| 精品乱码久久久久久99久播| 亚洲天堂国产精品一区在线| 男人舔女人下体高潮全视频| 久久草成人影院| 欧美在线黄色| 亚洲av电影不卡..在线观看| av国产免费在线观看| 欧美成人免费av一区二区三区| 美女大奶头视频| 国产成人欧美在线观看| 成人三级黄色视频| 精品久久蜜臀av无| 无遮挡黄片免费观看| 桃色一区二区三区在线观看| 日韩精品青青久久久久久| 亚洲人成电影免费在线| www.自偷自拍.com| 精品熟女少妇八av免费久了| 国内精品美女久久久久久| 国产成人影院久久av| 日本成人三级电影网站| 久久精品91无色码中文字幕| 欧美不卡视频在线免费观看| 一个人免费在线观看电影 | 午夜免费观看网址| 伦理电影免费视频| 国产私拍福利视频在线观看| 国产人伦9x9x在线观看| 国产私拍福利视频在线观看| 九色国产91popny在线| 宅男免费午夜| 久久久精品大字幕| 亚洲熟妇中文字幕五十中出| 成人一区二区视频在线观看| 亚洲精品粉嫩美女一区| 国产伦精品一区二区三区四那| 亚洲av日韩精品久久久久久密| 国产1区2区3区精品| 国产精品爽爽va在线观看网站| 韩国av一区二区三区四区| 少妇裸体淫交视频免费看高清| 99国产精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 免费观看的影片在线观看| 亚洲中文字幕日韩| 男人和女人高潮做爰伦理| 色哟哟哟哟哟哟| 后天国语完整版免费观看| 大型黄色视频在线免费观看| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 九色国产91popny在线| 久久精品影院6| 色综合婷婷激情| 亚洲av第一区精品v没综合| h日本视频在线播放|