• <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.

    亚洲国产精品999| 欧美激情极品国产一区二区三区| 国产精品一区二区在线观看99| 国产精品国产三级专区第一集| 国产成人91sexporn| 黄片播放在线免费| 丝袜在线中文字幕| 国产在线免费精品| 亚洲欧美成人综合另类久久久| 国产一区二区 视频在线| 国产片特级美女逼逼视频| 97在线人人人人妻| 欧美日韩亚洲综合一区二区三区_| 日本色播在线视频| 国产1区2区3区精品| 一级黄片播放器| 一区福利在线观看| 久久久亚洲精品成人影院| av又黄又爽大尺度在线免费看| 亚洲国产欧美网| 欧美日韩亚洲高清精品| 国产免费现黄频在线看| 老司机午夜十八禁免费视频| 最近手机中文字幕大全| 天天躁狠狠躁夜夜躁狠狠躁| 日本a在线网址| 免费不卡黄色视频| 五月天丁香电影| 少妇精品久久久久久久| 人成视频在线观看免费观看| 欧美日韩亚洲国产一区二区在线观看 | 国精品久久久久久国模美| 亚洲熟女毛片儿| 另类亚洲欧美激情| 日韩大码丰满熟妇| 青青草视频在线视频观看| 亚洲中文字幕日韩| 久久鲁丝午夜福利片| 欧美日韩视频精品一区| 在线观看免费高清a一片| 丝袜在线中文字幕| 国产高清不卡午夜福利| 性少妇av在线| 成年人午夜在线观看视频| av在线老鸭窝| 精品福利永久在线观看| 多毛熟女@视频| 99国产精品一区二区三区| 97在线人人人人妻| 男女下面插进去视频免费观看| 18禁裸乳无遮挡动漫免费视频| 午夜91福利影院| 韩国精品一区二区三区| 五月开心婷婷网| 香蕉国产在线看| 黄片播放在线免费| 少妇粗大呻吟视频| 岛国毛片在线播放| 国产免费福利视频在线观看| 亚洲成色77777| 97精品久久久久久久久久精品| 国产精品99久久99久久久不卡| 亚洲精品国产av成人精品| 在线 av 中文字幕| 欧美成人午夜精品| av天堂在线播放| 在线观看免费视频网站a站| 国产伦理片在线播放av一区| 国产男女超爽视频在线观看| 欧美精品一区二区大全| 一区在线观看完整版| 十分钟在线观看高清视频www| 男女高潮啪啪啪动态图| 人妻人人澡人人爽人人| 亚洲精品在线美女| 成人国产一区最新在线观看 | av福利片在线| 校园人妻丝袜中文字幕| 乱人伦中国视频| 国产成人精品在线电影| 色婷婷av一区二区三区视频| 老司机在亚洲福利影院| 欧美激情 高清一区二区三区| 只有这里有精品99| 高清不卡的av网站| 丝袜喷水一区| 久久99一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产片特级美女逼逼视频| 日韩欧美一区视频在线观看| 亚洲色图 男人天堂 中文字幕| 日本av手机在线免费观看| 亚洲av综合色区一区| 欧美激情极品国产一区二区三区| 久久久久久亚洲精品国产蜜桃av| 王馨瑶露胸无遮挡在线观看| 男人添女人高潮全过程视频| 久久99精品国语久久久| 亚洲成人手机| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 午夜久久久在线观看| 久久久久精品国产欧美久久久 | 亚洲自偷自拍图片 自拍| 亚洲中文字幕日韩| 国产真人三级小视频在线观看| 人人妻人人爽人人添夜夜欢视频| 免费女性裸体啪啪无遮挡网站| 美女大奶头黄色视频| 两性夫妻黄色片| 国产亚洲av片在线观看秒播厂| 1024香蕉在线观看| 免费女性裸体啪啪无遮挡网站| 日韩一本色道免费dvd| 久久久精品国产亚洲av高清涩受| 伦理电影免费视频| 欧美人与性动交α欧美软件| 黄色 视频免费看| 人妻 亚洲 视频| 国产精品一区二区在线观看99| 亚洲国产欧美一区二区综合| 性少妇av在线| 777久久人妻少妇嫩草av网站| 国产又色又爽无遮挡免| 18禁黄网站禁片午夜丰满| 99香蕉大伊视频| 搡老乐熟女国产| 啦啦啦中文免费视频观看日本| 亚洲精品日本国产第一区| 亚洲图色成人| 国产精品九九99| 成人手机av| 日韩大片免费观看网站| 高清不卡的av网站| a 毛片基地| 国产淫语在线视频| 在线天堂中文资源库| 久久天躁狠狠躁夜夜2o2o | 天天躁日日躁夜夜躁夜夜| 美女视频免费永久观看网站| 国产av一区二区精品久久| 日韩av不卡免费在线播放| 国产视频首页在线观看| 亚洲三区欧美一区| 91精品三级在线观看| 国产欧美日韩综合在线一区二区| 日日夜夜操网爽| 亚洲国产成人一精品久久久| 亚洲成色77777| 免费在线观看黄色视频的| 天堂8中文在线网| 啦啦啦啦在线视频资源| 久久亚洲国产成人精品v| 91国产中文字幕| 精品第一国产精品| 国产精品免费视频内射| 欧美黑人精品巨大| 女人久久www免费人成看片| 91国产中文字幕| 一级毛片女人18水好多 | cao死你这个sao货| xxx大片免费视频| av不卡在线播放| 久久狼人影院| 满18在线观看网站| 亚洲精品乱久久久久久| 欧美精品啪啪一区二区三区 | 欧美日韩亚洲高清精品| 欧美黄色片欧美黄色片| 免费观看人在逋| 少妇人妻久久综合中文| 一本久久精品| 精品一区二区三卡| 母亲3免费完整高清在线观看| 两人在一起打扑克的视频| 亚洲精品av麻豆狂野| 欧美日韩黄片免| 岛国毛片在线播放| 黄色片一级片一级黄色片| 午夜91福利影院| 久久狼人影院| 美女大奶头黄色视频| 欧美黄色淫秽网站| 91精品伊人久久大香线蕉| 97在线人人人人妻| 人人妻人人澡人人爽人人夜夜| 亚洲色图 男人天堂 中文字幕| 人人澡人人妻人| 天天躁夜夜躁狠狠躁躁| 久久女婷五月综合色啪小说| 国产精品.久久久| 汤姆久久久久久久影院中文字幕| bbb黄色大片| 9色porny在线观看| 黄色 视频免费看| 中文字幕人妻丝袜一区二区| 日韩人妻精品一区2区三区| 制服人妻中文乱码| 老司机靠b影院| 久久久欧美国产精品| 日本a在线网址| 一级毛片 在线播放| 亚洲欧美日韩高清在线视频 | 视频在线观看一区二区三区| 女人精品久久久久毛片| 高清不卡的av网站| 亚洲精品中文字幕在线视频| 宅男免费午夜| 国产精品 国内视频| 久久久国产一区二区| 2021少妇久久久久久久久久久| 国产在线一区二区三区精| 国产一区二区激情短视频 | 50天的宝宝边吃奶边哭怎么回事| 美女午夜性视频免费| 成年人黄色毛片网站| 国产在视频线精品| 一二三四社区在线视频社区8| 一本—道久久a久久精品蜜桃钙片| 韩国高清视频一区二区三区| 中文字幕人妻丝袜制服| 久久国产亚洲av麻豆专区| 国产主播在线观看一区二区 | 日本欧美视频一区| bbb黄色大片| 亚洲欧美精品综合一区二区三区| 国产精品av久久久久免费| 国产精品一区二区在线不卡| 国产日韩欧美在线精品| 在线观看免费午夜福利视频| 亚洲五月色婷婷综合| 黑丝袜美女国产一区| 亚洲欧洲精品一区二区精品久久久| av天堂在线播放| 久久99精品国语久久久| 777久久人妻少妇嫩草av网站| 19禁男女啪啪无遮挡网站| 欧美中文综合在线视频| 亚洲图色成人| 99精品久久久久人妻精品| 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| 青春草亚洲视频在线观看| 一个人免费看片子| 麻豆国产av国片精品| 下体分泌物呈黄色| 精品免费久久久久久久清纯 | 国产av精品麻豆| 久久热在线av| 成人午夜精彩视频在线观看| 中文字幕色久视频| 黑人巨大精品欧美一区二区蜜桃| 超色免费av| 久久精品亚洲熟妇少妇任你| 五月开心婷婷网| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 久久精品久久久久久噜噜老黄| 久久国产精品影院| 久久99一区二区三区| 国产免费视频播放在线视频| 国产午夜精品一二区理论片| av网站在线播放免费| 男人舔女人的私密视频| 超色免费av| www日本在线高清视频| 国产一区二区激情短视频 | 亚洲熟女毛片儿| 国产免费福利视频在线观看| av电影中文网址| 久久鲁丝午夜福利片| 激情五月婷婷亚洲| 亚洲精品久久成人aⅴ小说| 精品国产一区二区久久| 久久久久久人人人人人| 久久九九热精品免费| 手机成人av网站| 午夜福利在线免费观看网站| 亚洲精品av麻豆狂野| 91麻豆精品激情在线观看国产 | 一本久久精品| 精品久久久久久久毛片微露脸 | 精品国产国语对白av| 国产欧美日韩一区二区三 | 免费av中文字幕在线| 日本五十路高清| 香蕉国产在线看| 亚洲欧美一区二区三区国产| 亚洲国产精品一区二区三区在线| 亚洲欧美一区二区三区久久| 亚洲久久久国产精品| 免费久久久久久久精品成人欧美视频| 免费少妇av软件| 国产成人av教育| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 美女大奶头黄色视频| 久久精品国产亚洲av涩爱| 制服诱惑二区| 天天躁日日躁夜夜躁夜夜| 欧美97在线视频| 国产成人精品久久二区二区91| 国产片内射在线| 国产一级毛片在线| 精品国产一区二区三区久久久樱花| 大片电影免费在线观看免费| 久久人人爽人人片av| 成人亚洲欧美一区二区av| 丰满人妻熟妇乱又伦精品不卡| 性色av一级| 亚洲成色77777| 2018国产大陆天天弄谢| 亚洲中文字幕日韩| 国产精品九九99| 精品久久蜜臀av无| 丰满饥渴人妻一区二区三| 久久久国产精品麻豆| 色综合欧美亚洲国产小说| 亚洲精品国产区一区二| 亚洲av国产av综合av卡| 国产片特级美女逼逼视频| 少妇裸体淫交视频免费看高清 | 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 大码成人一级视频| 日本五十路高清| 国产亚洲精品久久久久5区| 99国产精品99久久久久| 看十八女毛片水多多多| 国精品久久久久久国模美| 欧美人与善性xxx| 在线观看国产h片| 考比视频在线观看| 国产精品一区二区在线不卡| 亚洲第一青青草原| 青春草亚洲视频在线观看| 午夜福利一区二区在线看| 国产精品免费大片| av视频免费观看在线观看| av网站在线播放免费| 欧美精品高潮呻吟av久久| 99热网站在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产av影院在线观看| 自线自在国产av| 99热国产这里只有精品6| 晚上一个人看的免费电影| 久久性视频一级片| 免费在线观看影片大全网站 | 欧美日韩av久久| 最近最新中文字幕大全免费视频 | 又大又爽又粗| 久久人人爽人人片av| 亚洲激情五月婷婷啪啪| 侵犯人妻中文字幕一二三四区| 国产精品 国内视频| 777久久人妻少妇嫩草av网站| 永久免费av网站大全| a级毛片黄视频| 人成视频在线观看免费观看| 狂野欧美激情性bbbbbb| 亚洲精品自拍成人| 亚洲熟女毛片儿| 国产精品国产三级专区第一集| 久久人人爽人人片av| 中国美女看黄片| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 9热在线视频观看99| 亚洲av日韩在线播放| 无限看片的www在线观看| 99精国产麻豆久久婷婷| 日日爽夜夜爽网站| 一区二区日韩欧美中文字幕| 亚洲九九香蕉| 叶爱在线成人免费视频播放| 国产欧美日韩综合在线一区二区| 国产免费视频播放在线视频| 久久精品国产a三级三级三级| 国产成人精品久久二区二区91| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久精品古装| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 国产成人系列免费观看| 后天国语完整版免费观看| 老司机亚洲免费影院| 欧美人与善性xxx| 国产亚洲欧美精品永久| 久久人人97超碰香蕉20202| 中国美女看黄片| 一本大道久久a久久精品| 91精品三级在线观看| 国产精品免费大片| 老熟女久久久| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 国产一区二区三区综合在线观看| 国产成人一区二区三区免费视频网站 | 19禁男女啪啪无遮挡网站| 欧美另类一区| 丁香六月天网| 精品国产一区二区久久| 女人被躁到高潮嗷嗷叫费观| 嫁个100分男人电影在线观看 | 在线看a的网站| 91精品伊人久久大香线蕉| 久久亚洲国产成人精品v| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 桃花免费在线播放| 亚洲av成人不卡在线观看播放网 | 中文精品一卡2卡3卡4更新| 国产亚洲精品第一综合不卡| 国产国语露脸激情在线看| 亚洲精品中文字幕在线视频| 久久精品成人免费网站| 手机成人av网站| 丝袜美腿诱惑在线| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜一区二区| 中国国产av一级| 在线观看一区二区三区激情| 久久九九热精品免费| 97人妻天天添夜夜摸| 成人亚洲欧美一区二区av| 国产成人免费无遮挡视频| 人妻一区二区av| www.熟女人妻精品国产| 精品视频人人做人人爽| 欧美黄色淫秽网站| 亚洲av综合色区一区| 天天躁夜夜躁狠狠久久av| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 欧美另类一区| 国产xxxxx性猛交| 精品国产一区二区久久| 亚洲精品中文字幕在线视频| 一级,二级,三级黄色视频| 超碰成人久久| 免费高清在线观看日韩| 国精品久久久久久国模美| 午夜福利,免费看| 亚洲精品国产一区二区精华液| 精品亚洲成a人片在线观看| 最黄视频免费看| 免费看av在线观看网站| 久久久亚洲精品成人影院| 老司机亚洲免费影院| 国产熟女欧美一区二区| 亚洲成人手机| 女性生殖器流出的白浆| 国产免费视频播放在线视频| 国产免费福利视频在线观看| 亚洲中文字幕日韩| 国产精品秋霞免费鲁丝片| 久久国产精品男人的天堂亚洲| 男女国产视频网站| 国产在线视频一区二区| 中文字幕精品免费在线观看视频| 欧美黄色片欧美黄色片| 国产黄色视频一区二区在线观看| 男女免费视频国产| 精品少妇内射三级| 久久av网站| 日本猛色少妇xxxxx猛交久久| 日本五十路高清| 黄色片一级片一级黄色片| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 少妇 在线观看| 男女下面插进去视频免费观看| 又黄又粗又硬又大视频| 亚洲成色77777| 91精品伊人久久大香线蕉| 国产精品三级大全| 欧美精品高潮呻吟av久久| 亚洲av男天堂| av在线app专区| 极品人妻少妇av视频| 天天躁夜夜躁狠狠躁躁| 最黄视频免费看| 美女视频免费永久观看网站| 欧美精品一区二区免费开放| 国产精品一二三区在线看| 又大又爽又粗| 婷婷成人精品国产| 少妇人妻久久综合中文| 精品人妻在线不人妻| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 日韩一区二区三区影片| 日本vs欧美在线观看视频| 一区二区日韩欧美中文字幕| 亚洲精品av麻豆狂野| 51午夜福利影视在线观看| 韩国高清视频一区二区三区| av片东京热男人的天堂| 午夜免费男女啪啪视频观看| 热re99久久国产66热| 国产女主播在线喷水免费视频网站| 免费看av在线观看网站| 亚洲成国产人片在线观看| 宅男免费午夜| 女人爽到高潮嗷嗷叫在线视频| 日韩一本色道免费dvd| 纯流量卡能插随身wifi吗| 色综合欧美亚洲国产小说| 欧美人与善性xxx| 最近手机中文字幕大全| 国产黄频视频在线观看| 亚洲少妇的诱惑av| 男女高潮啪啪啪动态图| 亚洲成色77777| 九色亚洲精品在线播放| 亚洲欧洲国产日韩| 80岁老熟妇乱子伦牲交| 99久久综合免费| a级毛片在线看网站| 成年av动漫网址| www.精华液| 岛国毛片在线播放| 精品福利观看| 少妇人妻 视频| 久热这里只有精品99| 91精品国产国语对白视频| 欧美精品人与动牲交sv欧美| 自线自在国产av| 免费av中文字幕在线| 久久久国产精品麻豆| 亚洲精品一卡2卡三卡4卡5卡 | a级毛片在线看网站| 三上悠亚av全集在线观看| 精品久久蜜臀av无| 在线精品无人区一区二区三| 我要看黄色一级片免费的| 嫁个100分男人电影在线观看 | av一本久久久久| 午夜福利免费观看在线| 久久久久国产精品人妻一区二区| 免费在线观看黄色视频的| 成年动漫av网址| 久久99热这里只频精品6学生| 欧美xxⅹ黑人| 亚洲国产欧美日韩在线播放| av网站在线播放免费| 国产欧美亚洲国产| 国产精品久久久久久精品电影小说| 观看av在线不卡| 午夜影院在线不卡| 国产99久久九九免费精品| 亚洲欧美一区二区三区黑人| 丝袜在线中文字幕| 又粗又硬又长又爽又黄的视频| 亚洲精品美女久久久久99蜜臀 | 成年人黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 青青草视频在线视频观看| 黄频高清免费视频| 欧美成人午夜精品| 国产精品亚洲av一区麻豆| 在线精品无人区一区二区三| 男人爽女人下面视频在线观看| 日本一区二区免费在线视频| 国产精品国产av在线观看| 精品久久久精品久久久| 日本av免费视频播放| 欧美在线黄色| 精品高清国产在线一区| 汤姆久久久久久久影院中文字幕| 波多野结衣av一区二区av| 免费在线观看日本一区| 亚洲国产欧美一区二区综合| 97在线人人人人妻| 欧美日韩视频高清一区二区三区二| 满18在线观看网站| 韩国精品一区二区三区| 精品第一国产精品| 高清黄色对白视频在线免费看| 91字幕亚洲| 久久午夜综合久久蜜桃| 国产成人精品久久久久久| 狠狠精品人妻久久久久久综合| 一级毛片电影观看| 一级毛片黄色毛片免费观看视频| 悠悠久久av| 欧美国产精品一级二级三级| 欧美日韩亚洲综合一区二区三区_| 亚洲精品在线美女| 亚洲精品国产一区二区精华液| 亚洲人成77777在线视频| 啦啦啦啦在线视频资源| 亚洲精品国产一区二区精华液| 亚洲av电影在线进入| 久久久久久免费高清国产稀缺| 高清av免费在线| 亚洲av电影在线进入| 精品国产超薄肉色丝袜足j| 2018国产大陆天天弄谢| 国产男女超爽视频在线观看| 亚洲精品国产区一区二| svipshipincom国产片| 一级毛片黄色毛片免费观看视频| 久久久久久人人人人人| 国产高清videossex| 男女之事视频高清在线观看 | 久久性视频一级片| 久久精品熟女亚洲av麻豆精品| 热re99久久国产66热| 成年人黄色毛片网站| 久久精品aⅴ一区二区三区四区| av电影中文网址|