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

    A 3D MIL-101@rGO composite as catalyst for efficient conversion of straw cellulose into valuable organic acid

    2022-06-20 08:00:10YuanSuMenghuLuRuidianSuWeizhiZhouXingXuQianLi
    Chinese Chemical Letters 2022年5期

    Yuan Su, Menghu Lu, Ruidian Su, Weizhi Zhou, Xing Xu,*, Qian Li,*

    a School of Mathematic and Quantitative Economics, Shandong University of Finance and Economics, Ji’nan 250014, China

    bShandong Key Laboratory of Water Pollution Control and Resource Reuse, School of Environmental Science and Engineering, Shandong University, Qingdao 266200, China

    c School of Civil Engineering, Shandong University, Ji’nan 250100, China

    Keywords:Straw cellulose MIL-101(Cr)Three-dimensional graphene Catalystic conversion Organic acid

    ABSTRACT Efficient conversion of straw cellulose to chemicals or fuels is an attracting topic today for the utilization of biomass to substitute for fossil resources.The development of catalysts is of vital importance.In this work, a composite catalyst metal-organic frameworks (MOFs) immobilized on three-dimensional reduced graphene oxide (3D-rGO) were synthesized by in situ growth of the MIL-101(Cr) within the 3DrGO matrix.The supporting of 3D-rGO guaranteed the dispersion and acid site density of MIL-101(Cr).The MIL-101(Cr)@3D-rGO nanocomposite possesses excellent catalytic activity, stability, recyclability and is an idea catalyst for the efficient degradation of straw cellulose into formic acid (FA), acetic acid (AA)and oxalic acid (OA).A maximum FA conversion rates of 95.36% was obtained by using MIL-101(Cr)@3DrGO(1:1) as catalyst and hydrothermal reaction at mild conditions of 200 °C for 1h in alkaline aqueous medium.The MIL-101(Cr)@3D-rGO nanocomposite can be reused with high catalytic activity without any collapse of structure or leaching of chromium.

    With growing concerns about energy shortages and environment degradation, finding renewable and green energy sources as alternatives to petroleum-based productions have received worldwide attention.Conversion of biomass to chemicals or liquid fuels are attracting much interest in recent years, among which the degradation of cellulose into high value-added platform compounds is regarded as one of the most important direction.Among various of valuable biomass-derived products such as hydroxymethylfurfural (HMF), aldols and small molecule acids [1,2].Formic acid (FA) has received much attention because it is closely related to fuel-cell-based hydrogen economy and has the potential to serve as a hydrogen generation/storage material with its high energy density.On the other hand, crops residue such as straw produced in agriculture exceed 72 millionton annually in China.Much of it was used as fuel for direct combustion which could not only arouse pollution of the environment but also cause waste of resources.Among them, wheat straw (WS), containing 40%-60%of natural cellulose, hemicelluloses and lignin, has the potential to become the plentiful raw materials source for chemical conversion to produce high value-added fuel [3].

    Generally, the conversion of cellulose into platform molecules is achieved by the oxidation of carbonhydrates in the presence of catalysts.Nanoparticle catalyst, solid acid catalyst and alkali metal catalyst have shown remarkable effect [4-6].For example, with the solid heteropoly acid Cs2.5H0.5PW12O40catalyst, fructose could produces HMF with high efficiency in a two-phase system [5].The Ru-based catalyst could convert the cellobiose and cellulose to a hexahydric alcohol with high selectivity [6].The composite catalyst of Ir/SiO2and zeolite can directly convert microcrystalline cellulose inton-butane [7].However, these reactions are generally heterogeneously catalyzed.And there are few studies on the effi-cient catalytic degradation of cellulose to small molecule acids.The presence of repeated 1,4-βglycosidic bonds and complex hydrogen bonding networks in cellulose make them very difficult to dissolve,thus reducing the efficiency of heterogeneous catalysis.Recently,studies have shown that cellulose can be dissolved in ionic liquid or alkaline aqueous solutions, and then can be homogeneously degraded to HMF or small molecular acids, which depends on the catalyst system [8].But the high conversion efficiency and selectivity of product should be improved.

    Metal-organic frameworks (MOFs) are a kind of organicinorganic hybrid materials with metal ions or metal clusters as nodes and organic ligands as skeletons.They have great advantages in controllable regular nanochannel structure, large specific surface area and porosity, and easy functionalization compared with traditional inorganic and organic materials.As a result, there has been widespread interest in the fields of gas storage and separation, chemical sensing, drug delivery, and catalysis [9-11].In particular, owing to the rich Lewis acid sites based on unsaturated metal centers, MOFs are regarded as promising catalytic materials [12].However, MOFs also have some shortcomings including poor water resistance and difficulty in recovery, which limit their widespread application.Among various MOFs, MIL-101(M) possesses a rigid mesoporous and microporous cage structure with a very large specific surface area, in which MO4(OH)2(M = Cr, Fe, Al) octahedron and terephthalic acid (BDC) are bridged in space.The special matrix of MIL-101(M) ensures its excellent hydrothermal stability and chemical stability, and has a greater advantage in adsorption and catalysis than other MOFs [13].Moreover, the catalytic performance of MIL-101 is mainly derived from the unsaturated sites of surface metals.Compared with the Fe-based and Al-based MILs materials, MIL-101(Cr) has more unsaturated sites and stronger stability, therefore, it has great application potential especially in the field of heterogeneous catalysis.

    In addition, in order to solve the defects of the material and ensure the catalytic activity of the material, it is particularly important to fix the MOFs on a suitable solid support [14].Graphene is an excellent carrier material with high specific surface area, stable thermochemical properties, strong mechanical strength and structural adjustability [15].However, the existence of van der Waals forces andπ-πsuperposition between graphene sheets tend to cause its recombination and aggregation, resulting in a significant decrease in specific surface area [16].Assembling two-dimensional graphene to create a three-dimensional graphene structure could overcome these defects and provide the possibility for immobilizing MOFs nanoparticles within its matrix [17].Meanwhile, the adhesion of the internal solid particles, especially MOFs with caged skeleton, can also support the 3D structure and prevent the agglomeration of the graphene sheets.Therefore, the combination of the two materials not only solves the problems of difficult recovery of MOFs and easy agglomeration of graphene, but also contributes to the dispersion of active sites on MOFs, which ensure the sufficient contact of catalyst and reactants to improve catalytic effi-ciency [18,19].

    In this paper, MIL-101(Cr) wasin-situsynthesized within the 3D-rGO matrix to form MIL-101(Cr)@3D-rGO nanocomposite.Then,it was used as catalyst to efficiently convert straw cellulose to valuable organic acid in a homogeneous alkaline solution.Furthermore, the appropriate hydrothermal conditions of the conversion reactions were systematically evaluated, and a high selectivity and yield of formic acid was obtained.

    The synthesis process of MIL-101(Cr)@3D-rGO(n) (n represents the mass ratio of MIL-101(Cr) and 3D-rGO) is shown in Fig.1a.The SEM image of MIL-101(Cr) (Fig.1b) showed that the prepared MIL-101(Cr) possessed a regular crystal morphology with a diameter of about 300–400 nm.The internal defects and wrinkles of 3D-rGO were clearly observed (Fig.1c).SEM image of MIL-101(Cr)@3D-rGO(1:1) (Fig.1d) clearly showed that plenty of MIL-101(Cr) nanocrystalline were formed inside the network, almost completely occupied the defects and wrinkles of the 3D-rGO surface, which demonstrated that MIL-101(Cr) was uniformly attached within the 3D-rGO grid.Compared with the SEM image of MIL-101(Cr)@3D-rGO(1:2) (Fig.S1 in Supporting information), it can be seen that MIL-101(Cr) in MIL-101(Cr)@3D-rGO(1:1) was more fully attached inside 3D-rGO.The element composition and distribution(C, O, N, Cr) of MIL-101(Cr)@3D-rGO(1:1) were detected in SEM element mapping (Fig.1e) and elemental analysis (Fig.1f).The C were observed uniformly distributed within the composite, showing the skeleton of 3D-rGO.The O and Cr images could fit into each other and O atoms were distributed near the Cr atom positions, pointing to Cr-O.The Cr/O atom ratio was about 1:1, which was corresponding to the elemental composition of MIL-101(Cr).The scattered N atoms were attributed to the Cr(NO3)·9H2O during the preparation.And visible images of EDS mapping results confirmed that MIL-101(Cr) were successfully synthesized and dispersed uniformly within the 3D-rGO matrix.The TEM image of MIL-101(Cr)@3D-rGO(1:1) (Fig.S2 in Supporting information) also showed the regular crystal morphology of MIL-101(Cr).

    Fig.1.Preparation flow chart of MIL-101(Cr)@3D-rGO (a).SEM images of MIL-101(Cr), 3D-rGO and MIL-101(Cr)@3D-rGO(1:1) (b-d); SEM element mapping of MIL-101(Cr)@3D-rGO(1:1) (e).Elemental analysis of MIL-101(Cr)@3D-rGO(1:1) (f).

    The XRD patterns of 3D-rGO, MIL-101(Cr) and MIL-101(Cr)@3DrGO composite were exhibited in Fig.2a.For 3D-rGO pattern, a distinct diffraction peak appeared near 21.50°, which was consistent with the corresponding surface (002) of 3D-rGO.The main peaks of MIL-101(Cr) were recorded at 5.82°, 8.33°, 9.02°, 16.39°, corresponding to (531), (882), (911), (311) of standard MOFs, respectively [20].For the MIL-101(Cr)@3D-rGO composite, both characteristic peaks of MIL-101(Cr) and 3D-rGO were detected indicating that the addition of 3D-rGO would not affect the crystal structure of MIL-101(Cr).A stronger characteristic peak was observed as the increasing of the content of MIL-101(Cr) in the composite.

    The FT-IR spectrum of the catalysts were depicted in Fig.2b.The bands around 748 cm-1were related to the deformation vibration of C–H.A adsorption band appeared around 1015 cm-1was assigned to the C–O stretching.And the bands at 1348 cm-1could be attributed to –OCO group from the dicarboxylate linker.The peaks observed around 1507 cm-1and 1620 cm-1were assigned to aromatic C=C stretches and structural water molecules,respectively [21].It can be seen that the composite material has a good characteristic peak and the addition of 3D-rGO does not affect the formation and structure of MIL-101(Cr).

    Fig.2.XRD patterns (a) and FTIR spectra (b) of MIL-101(Cr), 3D-rGO MIL-101(Cr)@3D-rGO(1:1) and MIL-101(Cr)@3D-rGO(1:2); NH3-TPD curve (c) of MIL-101(Cr)@3DrGO(1:1) and MIL-101(Cr)@3D-rGO(1:2); TGA curve (d) of MIL-101(Cr)@3D-rGO(1:1) and MIL-101(Cr)@3D-rGO(1:2).

    The NH3-TPD curve is used to characterize the acidity of the composite catalyst surface.From Fig.2c, it can be observed that two different types of 3D-rGO dosing compound catalysts had detected acid sites with different strengths, of which MIL-101(Cr)@3D-rGO(1:1) possessed a stronger Bronsted acidic site than MIL-101(Cr)@3D-rGO(1:2).The difference in the strength of the Lewis acid sites between the two composites was not large,which was attributed to the higher proportion of MIL-101(Cr) in MIL-101(Cr)@3D-rGO(1:1) composites, of which had more activated Cr acid sites.NH3-TPD curve confirmed that the composite material was an acidic catalyst and MIL-101(Cr)@3D-rGO(1:1) had a stronger acidic site.

    The TGA curve of the composite material is shown in Fig.2d,reflecting the thermal stability of the material.The weight loss of the composite material was mainly divided into three stages.The first stage was 20–360 °C, the weight loss at this stage was due to removal of trace water on the surface of the material (8 wt% for MIL-101(Cr)@3D-rGO(1:1), 14 wt% for MIL-101(Cr)@3D-rGO(1:2)).The second stage occurred at 360–550 °C, the weight loss at this stage is attributed to the collapse skeleton of the composite catalyst caused by the removal of the material skeleton linker (terephthalic acid).The third stage took place when the temperature was higher than 550 °C, of which the remaining weight was assigned to chromium oxide after combustion.In addition, the weight loss of MIL-101(Cr)@3D-rGO(1:1) and MIL-101(Cr)@3D-rGO(1:2) are 67 wt% and 49 wt%, respectively, which proved that as the proportion of 3D-rGO increased, the thermal stability of the composite material was enhanced.At the same time, the composite material can keep the skeleton intact under 360 °C, which can be used for the catalytic reaction of cellulose.

    The XPS survey spectra of the as-prepared MIL-101(Cr)@3DrGO(1:1) were detected in Fig.S3 (Supporting information).The photoelectron peaks of main elements on the surface of MIL-101(Cr)@3D-rGO(1:1) appeared at the binding energies of 284 eV(C 1s), 534 eV (O 1s) and 578 eV (Cr 2p), respectively, confirming the presence of MIL-101(Cr) on the framework of 3D-rGO.The binding energy at 289.1 eV and 284.5 eV of C 1s can be assigned to C=O and C=C, respectively.Among them, the C=C peak is mainly attributed to 3D-rGO, and the C=O peak is mainly MIL-101(Cr).The presence of hydroxyl and epoxy C–O was not observed at 286.5 eV,indicating that the reduction effect of 3D-rGO was relatively obvious and the network composition of graphene was mainly C=C.The Cr 2p XPS spectrum showed two main peaks appeared at 586.97 eV and 577.10 eV corresponding to Cr 2p1/2and Cr 2p2/3,respectively, confirming the embedded Cr in MIL-101(Cr) [22].The binding energy of 531.5 eV was due to the coordination of Cr atom with the carboxyl oxygen to form Cr–O.The XPS results showed that MIL-101(Cr) was successfullyin situsynthesized within 3DrGO matrix.The Raman spectroscopy spectrum of the as-prepared MIL-101(Cr)@3D-rGO(1:1) was shown in Fig.S4 (Supporting information).The characteristic D band at 1348 cm-1was the disordered vibration peak of graphene, indicating the exist of defect on the surface of 3D-rGO.The main characteristic peak G of graphene at 1583 cm-1was caused by the in-plane vibration of sp2carbon atoms.The peaks of 2590 cm-1and 2700 cm-1are 2D bands of graphene, which is a combination of G’peak and D’peak.TheID/IGof MIL-101(Cr)@3D-rGO is 1.04, indicating that there are more defects in 3D-rGO for MIL-101(Cr) synthesisin situ.

    The pore sizes and specific surface areas of different asprepared catalysts samples were measured and shown in Fig.S5(Supporting information) and Table S1 (Supporting information).MIL-101(Cr) had a high N2adsorption capacity and exhibited a mesoporous structure with a pore diameter of 4.728 nm, while 3DrGO was macroporous and the average pore diameter was 2.420 μm.Meanwhile, the MIL-101(Cr)@3D-rGO exhibited a significant decrease of porous surface area with the increase of MIL-101(Cr)loading levels compared with MIL-101(Cr) (2575 m2/g) and 3D-rGO(1774 m2/g).The appreciable decreases in surface area indicated that the pore surface of 3D-rGO and MIL-101(Cr) were presumably occupied by each other and the cage structure MIL-101(Cr) was successfully interwoven with the 3D network of rGO.The porous structure of MIL-101(Cr)@3D-rGO facilitated full exposure of the acidic sites of Cr, which greatly enhanced their contact with the cellulose solution penetrating into the 3D-rGO matrix and thus improved catalytic property for cellulose conversion.

    The cellulose conversion performance catalyzed by four composite catalysts were shown in Fig.3a.The target product was small molecule acid such as FA and AA, which could be directly detected by HPLC.It could be observed that the all as-prepared catalysts had obvious catalytic property for cellulose degradation,of which MIL-101(Cr)@3D-rGO(1:1) was the best.The product concentration of FA obtained from the degradation of cellulose for 1 h catalyzed by different materials was 0.0044 g/mL for no catalyst, 0.0067 g/mL for pure MIL-101(Cr), 0.0062 g/mL for pure 3DrGO, 0.0572 g/mL for MIL-101(Cr)@3D-rGO(1:1) and 0.0102 g/mL for MIL-101(Cr)@3D-rGO(1:2), respectively.According to the original content of cellulose in the uniform cellulose solution (0.06 g/mL), the FA conversion rates catalyzed by different catalysts were calculated to be 7.3%, 11.1%, 10.3%, 95.4% and 16.9%, respectively.Similarly, the product concentration and conversion rates of AA were 0.0013, 0.0014, 0.0015, 0.0019, 0.0026 g/mL and 2.1%, 2.3%,2.4%, 3.1%, 4.3%, respectively.Without catalyst, the total conversion rates of obtained FA and AA was only 9.40%, and the HPLC peak(Fig.S6 in Supporting information) was relatively complicated indicating that there were many other by-products.When the degradation was catalyzed by pure 3D-rGO and MIL-101(Cr), the amount of FA was still low, while the conversion rates of AA was increased slightly and the total conversion rates of FA and AA was 13.60% and 13.45%, respectively.For MIL-101(Cr)@3D-rGO(1:1), the conversion rates of FA had an amazing increase to reach a maximum of 95.35%while the total conversion rates of all the products reached 98.48%.This indicated that the composite material with proper amount of 3D-rGO and MIL-101(Cr) could present a synergistic effect, with which the cellulose solution would be converted to FA efficiently.In the composite catalyst, MIL-101(Cr) promoted cellulose conversion to FA and the addition of an appropriate amount of 3D-rGO inhibited the production of by-products [17].Moreover, an appropriate amount of 3D-rGO doping in the composite catalyst could contribute to the dispersion of MIL-101(Cr), making the acidic site of Cr in MIL-101(Cr) sufficiently contact with the cellulose alkaline solution.However, when the 3D-rGO was excessively doped, the acidic sites of Cr in the composite catalyst reduced, which is insufficient for the efficient degradation of cellulose and resulted in the conversion rates of FA and the total conversion rates of cellulose catalyzed by MIL-101(Cr)@3D-rGO(1:2) drastically dropped to 16.99% and 21.26%.Therefore, the MIL-101(Cr)@3D-rGO(1:1) composite possessed the best degradation rate and selectivity for the conversion of homogeneous cellulose.The comparation of the cellulose conversion performance catalyzed by the as-prepared MIL-101(Cr)@3D-rGO(1:1) with the other reported catalysts were listed in Table S2 (Supporting information).It can be seen that MIL-101(Cr)@3D-rGO(1:1) can efficiently decompose cellulose into a small molecular acid in an alkaline system at a relatively mid temperature and has great advantages in cellulose conversion reaction.

    Fig.3.Catalytic degradation of cellulose by different catalysts (reaction conditions were 200°C, 1 h, 0.1 g of catalyst dosage) (a).Effect of reaction time and reaction temperature of MIL-101(Cr)@3D-rGO(1:1) on cellulose conversion products (b).Effect of the catalyst dosage on the conversion products of cellulose (reaction conditions were 200°C, 1 h) (c).Circulation experiment of MIL-101(Cr)@3D-rRGO(1:1) catalyst (d).XPS spectrum of the recovered MIL-101(Cr)@3D-rGO(1:1) composite catalyst after 5 cycles of experiment (e).Cr content in leachate of different catalysts (f).

    The effect of reaction time and temperature on cellulose degradation products were described in Fig.3b.As the time increased from 0.5 h to 2 h, the conversion rates of FA at 50 °C and 100°C was continuously increased from 1.883% and 31.83% to 3.850%and 53.33%, respectively.When the degradation temperature rises to 150 °C and 200 °C, the produced FA conversion rates reached maximum at 1 h.Continued to prolong reaction time will result in a decrease in FA production, which was attributed to the decomposition of FA into carbon dioxide and hydrogen at temperatures above 150 °C.It is obvious that the hydrolysis temperature is the dominant influencing factor to convert cellulose into FA using MIL-101(Cr)@3D-rGO(1:1) as catalyst, whereas reaction time has little effect [23].But for AA, the conversion rates did not change significantly with hydrolysis temperature or time, which remained at around 4% and reached a maximum of 6.36% at 50 °C, 2 h (Figs.S7-S9 in Supporting information).It is clear that in the presence of the as-prepared composite catalyst, the conversion rate of FA could be greatly improved under certain time and temperature conditions, which implies that we can achieve high conversion selectivity of FA by adjusting hydrolysis temperature and time using MIL-101(Cr)@3D-rGO(1:1) as catalyst.

    The effect of the amount of catalyst on the degradation of cellulose was shown in Fig.3c.As the dosage of the catalyst increasing, the conversion rates of FA enhanced slightly, while the conversion rates of AA declined, resulting in the overall change exceed no more than 1%, which indicated that the dosage of the catalyst had slight impact on cellulose conversion.

    Recyclability is one of the criteria for evaluating the pros and cons of the catalyst.The composite of MIL-101(Cr)@3D-rGO not only showed excellent catalytic property, but also exhibited more conducive to recycling.Thus, cycle experiment were conducted to evaluate the reusability MIL-101(Cr)@3D-rGO and the results were shown in Fig.3d.After repeated usages for 3 times, the conversion rates of FA decreased slightly from 95.05% to 85.16% and the total conversion rates decreased by 9.89%.Repeated catalytic processes may cause the MIL-101(Cr) acidic sites to be entrapped by carbon deposits, resulting in a decrease in catalytic performance.However, after repeated use of four times, the total conversion rates of small molecule acid is still 83.26%, which exhibited the excellent recycling property of the as-prepared composite catalyst.Fig.3e showed the XPS spectra of the composite catalyst after five cycles experiments, compared to the XPS pattern of the catalyst before the catalytic experiment (Fig.S3), the Cr 2p2/3peak at 577.10 eV and the C=O peak at 289.1 eV showed a slight decrease.This is attributed to the fact that after the catalyst is catalyzed, some carbon deposits adhere to the surface of the material and affect the catalytic strength of the acidic sites of Cr.The newly occurring hydroxyloxy peak at 536.23 eV may be derived from a cellulose alkaline solution [24,25].

    To evaluated the stability and environmental safety of the composite catalyst, the emission of chromium during the reaction process was determined (Fig.3f).When the pure MIL-101(Cr) was used as catalyst, a large amount of chromium ion leakage was detected.But for MIL-101(Cr)@3D-rGO, chromium ion leaching of the composite catalyst is significantly reduced with the increase of the proportion of 3D-rGO, and the chromium concentrations were lower than the national wastewater discharge standard (GB-8978–1996).So, the composite with 3D-rGO not only enhances the catalytic performance of the catalyst, but also enhances the stability of MIL-101(Cr).

    Since cellulose is formed by the linkage of glucose through a 1,4-β-glycosidic bond through hydrogen bonding, there are repeated 1,4-β-glycosidic bonds and complex hydrogen bonding networks inside the cellulose, which limits the dissolution and degradation of cellulose.In this study, after dissolving cellulose in an alkaline aqueous solution with sodium hydroxide solution, each cellulose chain is surrounded by sodium hydroxide, which makes the hydrogen bond network inside the cellulose being destroyed[26].Thus, as the catalyst is present, the cellulose macromolecule is decomposed into small molecule glucose under the action of Cr acidic site.In the catalytic process, Cr acid site of the composite catalyst accepts electrons, the glucose molecule loses electrons and the chain is progressively decomposed to glycolaldehyde, glyoxylic acid, oxalic acid and finally obtain FA.Or, glucose is converted to glycerol by hydrolysis, then AA was obtained by further hydrolysis through levulinic acid.The potential conversion pathways of cellulose was shown in Fig.4 and Fig.S10 (Supporting information).

    Fig.4.The potential degradation pathways of cellulose.

    In the present research, a new composite catalyst MIL-101(Cr)@3D-rGO(1:1) was synthesized byin-situsynthesis and used for the efficient degradation and conversion of wheat straw cellulose.The MIL-101(Cr) with about 300–400 nm were tightly wrapped inside the 3D-rGO matrix and the acidic site of Cr provided a catalytic possibility for cellulose degradation.When the cellueous was catalyzed by MIL-101(Cr)@3D-rGO(1:1) at 200 °C for 1 h, an ultra-high FA conversion rate (95.36%) could be obtained.Compared with the present catalysts, MIL-101(Cr)@3D-rGO prepared in this study could possess great potential in the energization of biomass.

    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 National Natural Science Foundation of China (Nos.U1906221 and 52070121).This work was also supported by grants from the Major Program of Shandong Province Technological Innovation Project (No.2020CXGC011403).

    Supplementary materials

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

    国精品久久久久久国模美| 亚洲av在线观看美女高潮| 日本猛色少妇xxxxx猛交久久| 夫妻性生交免费视频一级片| 欧美3d第一页| 成人影院久久| 一区在线观看完整版| 国产xxxxx性猛交| 热re99久久国产66热| 日产精品乱码卡一卡2卡三| videos熟女内射| 丰满少妇做爰视频| 自拍欧美九色日韩亚洲蝌蚪91| 肉色欧美久久久久久久蜜桃| 亚洲美女搞黄在线观看| av不卡在线播放| 亚洲综合色惰| kizo精华| 国产精品一区www在线观看| 最近最新中文字幕免费大全7| 亚洲精品美女久久久久99蜜臀 | 国产亚洲精品第一综合不卡 | 久久热在线av| 90打野战视频偷拍视频| 女人被躁到高潮嗷嗷叫费观| 搡女人真爽免费视频火全软件| 亚洲成av片中文字幕在线观看 | 久久精品国产a三级三级三级| 一级a做视频免费观看| 欧美最新免费一区二区三区| 国产精品久久久久成人av| 22中文网久久字幕| 国产一区有黄有色的免费视频| 亚洲综合精品二区| 国产精品三级大全| 波多野结衣一区麻豆| 国产成人精品在线电影| 国产成人精品一,二区| 国产精品久久久久久精品电影小说| 高清av免费在线| 国产极品天堂在线| 欧美日本中文国产一区发布| 高清在线视频一区二区三区| 老司机亚洲免费影院| 成人漫画全彩无遮挡| 99香蕉大伊视频| 免费av不卡在线播放| 国产高清三级在线| 亚洲精品成人av观看孕妇| 久久精品久久精品一区二区三区| 欧美人与性动交α欧美精品济南到 | 午夜视频国产福利| 亚洲图色成人| 亚洲高清免费不卡视频| 欧美精品av麻豆av| 久久精品国产亚洲av涩爱| 热99国产精品久久久久久7| 99热网站在线观看| 少妇高潮的动态图| 丰满饥渴人妻一区二区三| 国产精品免费大片| 亚洲激情五月婷婷啪啪| 精品人妻在线不人妻| av女优亚洲男人天堂| 男女边吃奶边做爰视频| 汤姆久久久久久久影院中文字幕| 国产探花极品一区二区| 高清在线视频一区二区三区| 国产亚洲精品久久久com| 亚洲内射少妇av| av播播在线观看一区| 国产白丝娇喘喷水9色精品| 免费少妇av软件| 久久狼人影院| 久久久精品区二区三区| 国产成人免费无遮挡视频| 中文精品一卡2卡3卡4更新| 精品一区二区三区四区五区乱码 | 国产精品一二三区在线看| 好男人视频免费观看在线| 自拍欧美九色日韩亚洲蝌蚪91| av播播在线观看一区| 国产精品久久久久久久电影| 成年动漫av网址| 中文字幕制服av| 欧美人与性动交α欧美软件 | 国产在线一区二区三区精| 丝袜脚勾引网站| 两个人看的免费小视频| 日韩伦理黄色片| 99精国产麻豆久久婷婷| kizo精华| 久久精品久久久久久久性| 国产成人a∨麻豆精品| 最近最新中文字幕免费大全7| 日韩av不卡免费在线播放| 90打野战视频偷拍视频| 精品一品国产午夜福利视频| 99久久人妻综合| 韩国av在线不卡| 嫩草影院入口| 多毛熟女@视频| 美女主播在线视频| 永久免费av网站大全| 国产欧美亚洲国产| av视频免费观看在线观看| 18禁观看日本| 热99国产精品久久久久久7| 亚洲色图综合在线观看| 激情五月婷婷亚洲| 国产国语露脸激情在线看| 亚洲国产精品一区二区三区在线| 国产精品国产av在线观看| 国产一区有黄有色的免费视频| 丝袜在线中文字幕| 99久国产av精品国产电影| 少妇猛男粗大的猛烈进出视频| 精品久久蜜臀av无| 超色免费av| 天堂中文最新版在线下载| 久久鲁丝午夜福利片| 考比视频在线观看| 十八禁高潮呻吟视频| 欧美3d第一页| 51国产日韩欧美| 日日爽夜夜爽网站| 制服人妻中文乱码| 国产成人免费观看mmmm| 日产精品乱码卡一卡2卡三| 国产有黄有色有爽视频| 国产在线视频一区二区| 天天影视国产精品| 中文字幕亚洲精品专区| 国产黄色视频一区二区在线观看| 丝袜美足系列| 在线观看美女被高潮喷水网站| 九色成人免费人妻av| 欧美日韩视频精品一区| 国产精品一区www在线观看| 亚洲欧洲日产国产| 精品久久久久久电影网| 91国产中文字幕| 成人毛片60女人毛片免费| 亚洲成人av在线免费| 国产不卡av网站在线观看| 欧美精品高潮呻吟av久久| av线在线观看网站| 另类精品久久| 51国产日韩欧美| 激情视频va一区二区三区| 黄色视频在线播放观看不卡| 又黄又粗又硬又大视频| 国产免费一区二区三区四区乱码| 人人澡人人妻人| 91久久精品国产一区二区三区| 91久久精品国产一区二区三区| 中文字幕最新亚洲高清| 日本猛色少妇xxxxx猛交久久| 在线观看三级黄色| 女人精品久久久久毛片| 国产精品免费大片| 国产成人免费无遮挡视频| 国产精品一区二区在线不卡| 国产成人精品婷婷| 亚洲国产精品专区欧美| 亚洲国产av新网站| 国产极品粉嫩免费观看在线| 美女主播在线视频| 看非洲黑人一级黄片| 街头女战士在线观看网站| 亚洲激情五月婷婷啪啪| 亚洲一区二区三区欧美精品| 男的添女的下面高潮视频| 色94色欧美一区二区| 国产免费现黄频在线看| 性高湖久久久久久久久免费观看| 亚洲欧美成人精品一区二区| 久久97久久精品| av女优亚洲男人天堂| 宅男免费午夜| 美国免费a级毛片| videos熟女内射| 美女内射精品一级片tv| 激情五月婷婷亚洲| 成人国语在线视频| 欧美bdsm另类| 国产永久视频网站| 精品国产一区二区三区久久久樱花| 国产精品国产av在线观看| 成年动漫av网址| 精品亚洲成a人片在线观看| 免费看av在线观看网站| 在线免费观看不下载黄p国产| 欧美日韩av久久| videossex国产| 男男h啪啪无遮挡| 久久综合国产亚洲精品| 国国产精品蜜臀av免费| 我的女老师完整版在线观看| 国产精品一国产av| 日韩av在线免费看完整版不卡| 日本-黄色视频高清免费观看| 成人免费观看视频高清| 亚洲国产日韩一区二区| 国产午夜精品一二区理论片| 熟女av电影| av网站免费在线观看视频| 免费不卡的大黄色大毛片视频在线观看| 中文字幕av电影在线播放| 久久久精品94久久精品| 黑人猛操日本美女一级片| 爱豆传媒免费全集在线观看| 久久精品国产自在天天线| 国产黄色免费在线视频| 亚洲av中文av极速乱| 国产成人精品无人区| 中文字幕人妻丝袜制服| 99国产综合亚洲精品| 男女免费视频国产| 亚洲av.av天堂| 亚洲伊人久久精品综合| 日本vs欧美在线观看视频| 我的女老师完整版在线观看| 国产成人免费观看mmmm| 久久久久久伊人网av| 午夜激情久久久久久久| 午夜激情久久久久久久| 2018国产大陆天天弄谢| 欧美另类一区| 国产国语露脸激情在线看| 成人手机av| 国产精品秋霞免费鲁丝片| 亚洲综合色惰| 亚洲av男天堂| 日韩三级伦理在线观看| 黑人猛操日本美女一级片| 丰满迷人的少妇在线观看| 亚洲欧美日韩卡通动漫| 日韩欧美精品免费久久| 国产精品嫩草影院av在线观看| 热re99久久精品国产66热6| 激情视频va一区二区三区| 欧美成人精品欧美一级黄| 亚洲精品一区蜜桃| 一级爰片在线观看| 日韩电影二区| 欧美激情 高清一区二区三区| 伦精品一区二区三区| 97超碰精品成人国产| 人成视频在线观看免费观看| 国产麻豆69| 亚洲精品乱久久久久久| 国产精品三级大全| 9色porny在线观看| 久久97久久精品| 亚洲高清免费不卡视频| 一级毛片我不卡| 麻豆乱淫一区二区| 久久久久久人人人人人| 久久精品国产鲁丝片午夜精品| 久久精品国产鲁丝片午夜精品| 久久这里有精品视频免费| 18禁国产床啪视频网站| 看免费成人av毛片| 亚洲成人av在线免费| 最近最新中文字幕大全免费视频 | 国产深夜福利视频在线观看| 新久久久久国产一级毛片| 两个人免费观看高清视频| 亚洲成av片中文字幕在线观看 | 三级国产精品片| 久久久久久久大尺度免费视频| 亚洲精品一区蜜桃| 亚洲人成网站在线观看播放| 亚洲性久久影院| 全区人妻精品视频| 国产精品国产三级专区第一集| 久久久久久伊人网av| 国产日韩一区二区三区精品不卡| 咕卡用的链子| 韩国高清视频一区二区三区| 久久精品国产a三级三级三级| 赤兔流量卡办理| 国产av精品麻豆| 18在线观看网站| 女人精品久久久久毛片| 免费观看性生交大片5| 久久久久精品性色| 18禁观看日本| h视频一区二区三区| 伦理电影大哥的女人| 永久免费av网站大全| 又大又黄又爽视频免费| 多毛熟女@视频| 丝袜美足系列| 激情视频va一区二区三区| 又黄又爽又刺激的免费视频.| 99国产精品免费福利视频| 男女免费视频国产| 少妇人妻久久综合中文| 视频中文字幕在线观看| 色婷婷久久久亚洲欧美| 欧美 日韩 精品 国产| 日韩视频在线欧美| 国产一区二区三区综合在线观看 | 观看美女的网站| 成年人午夜在线观看视频| 成年av动漫网址| 日韩在线高清观看一区二区三区| 精品国产露脸久久av麻豆| 在线观看免费日韩欧美大片| 女人精品久久久久毛片| 久久精品久久精品一区二区三区| 高清欧美精品videossex| 丝袜在线中文字幕| 国产淫语在线视频| 成人亚洲欧美一区二区av| 亚洲av国产av综合av卡| 少妇的丰满在线观看| 一二三四中文在线观看免费高清| 美女国产高潮福利片在线看| 啦啦啦啦在线视频资源| 国产色婷婷99| 欧美激情极品国产一区二区三区 | av国产久精品久网站免费入址| 国产日韩欧美视频二区| 亚洲欧美色中文字幕在线| 少妇高潮的动态图| 婷婷色麻豆天堂久久| 国产成人免费观看mmmm| 在线观看免费高清a一片| 亚洲欧美精品自产自拍| 日本午夜av视频| 亚洲国产欧美在线一区| 免费观看无遮挡的男女| 侵犯人妻中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 国产熟女欧美一区二区| 少妇的丰满在线观看| 亚洲四区av| 亚洲国产精品999| 男的添女的下面高潮视频| 久久久久精品久久久久真实原创| 天堂8中文在线网| 熟女人妻精品中文字幕| 久久久久久久亚洲中文字幕| 国产白丝娇喘喷水9色精品| 久久99精品国语久久久| 久久久国产一区二区| 亚洲高清免费不卡视频| 亚洲精品国产av蜜桃| 人人妻人人爽人人添夜夜欢视频| av网站免费在线观看视频| 亚洲国产欧美在线一区| 一本久久精品| 久久精品国产a三级三级三级| 大香蕉久久网| 亚洲熟女精品中文字幕| 9色porny在线观看| 99热6这里只有精品| h视频一区二区三区| 日韩人妻精品一区2区三区| 在线观看国产h片| 国产69精品久久久久777片| 欧美激情极品国产一区二区三区 | 少妇人妻精品综合一区二区| 免费黄频网站在线观看国产| 成人亚洲精品一区在线观看| 男人操女人黄网站| 亚洲欧洲精品一区二区精品久久久 | 美女脱内裤让男人舔精品视频| 一级片免费观看大全| 欧美精品国产亚洲| 好男人视频免费观看在线| 久久久久久人妻| 亚洲精品美女久久久久99蜜臀 | 国产不卡av网站在线观看| 免费大片黄手机在线观看| 日韩av不卡免费在线播放| 久久久久久伊人网av| av天堂久久9| 少妇猛男粗大的猛烈进出视频| 9热在线视频观看99| 天天躁夜夜躁狠狠躁躁| 在线精品无人区一区二区三| 亚洲激情五月婷婷啪啪| 日本av免费视频播放| 久久久久视频综合| 亚洲成人手机| 国产成人欧美| 欧美变态另类bdsm刘玥| 永久免费av网站大全| 亚洲伊人久久精品综合| 午夜免费鲁丝| 亚洲精品日本国产第一区| 丝袜喷水一区| www.av在线官网国产| 黄片无遮挡物在线观看| 一级a做视频免费观看| 亚洲综合精品二区| 国产精品一区www在线观看| 精品一区二区三卡| 在线天堂最新版资源| 亚洲av.av天堂| 日韩欧美精品免费久久| 99久久人妻综合| 久久久国产精品麻豆| 国产一区二区三区av在线| 亚洲 欧美一区二区三区| 如何舔出高潮| 女人久久www免费人成看片| 日本av免费视频播放| xxxhd国产人妻xxx| 中文字幕亚洲精品专区| 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| av在线播放精品| 如何舔出高潮| 高清视频免费观看一区二区| 最新中文字幕久久久久| 亚洲国产最新在线播放| tube8黄色片| 最近中文字幕2019免费版| 美女国产视频在线观看| 夜夜骑夜夜射夜夜干| 国产精品麻豆人妻色哟哟久久| 国产精品久久久久久久电影| 一本—道久久a久久精品蜜桃钙片| 在线观看免费高清a一片| 中文字幕免费在线视频6| 久久精品人人爽人人爽视色| 久久这里只有精品19| 午夜福利视频精品| 黄色视频在线播放观看不卡| 日韩一区二区三区影片| 曰老女人黄片| 最后的刺客免费高清国语| 久久人人爽av亚洲精品天堂| 免费久久久久久久精品成人欧美视频 | 国产日韩欧美视频二区| 日本91视频免费播放| 最近最新中文字幕大全免费视频 | 99热国产这里只有精品6| 国产1区2区3区精品| av线在线观看网站| 香蕉丝袜av| 午夜久久久在线观看| 亚洲精品第二区| 91精品三级在线观看| 菩萨蛮人人尽说江南好唐韦庄| av免费观看日本| 大香蕉久久网| 日本色播在线视频| 99久久精品国产国产毛片| 久热这里只有精品99| 黄片无遮挡物在线观看| 热re99久久国产66热| 高清毛片免费看| 一级毛片黄色毛片免费观看视频| 精品午夜福利在线看| 免费日韩欧美在线观看| 男女无遮挡免费网站观看| 2022亚洲国产成人精品| 久久久久久久精品精品| 亚洲av在线观看美女高潮| 寂寞人妻少妇视频99o| 91久久精品国产一区二区三区| 香蕉精品网在线| 成人亚洲欧美一区二区av| 国产精品三级大全| 一边摸一边做爽爽视频免费| 捣出白浆h1v1| 美女国产视频在线观看| av电影中文网址| 丝袜喷水一区| 丝袜人妻中文字幕| 视频中文字幕在线观看| 亚洲图色成人| 国产成人精品无人区| 日本猛色少妇xxxxx猛交久久| 最近2019中文字幕mv第一页| 捣出白浆h1v1| 日韩制服骚丝袜av| 亚洲国产色片| 国产伦理片在线播放av一区| 国产精品无大码| 精品久久久精品久久久| 精品人妻在线不人妻| 国产 一区精品| 久久午夜福利片| 成人午夜精彩视频在线观看| 69精品国产乱码久久久| 国产在视频线精品| 久久毛片免费看一区二区三区| 五月伊人婷婷丁香| 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 草草在线视频免费看| www日本在线高清视频| 国产精品不卡视频一区二区| 波多野结衣一区麻豆| 国产成人a∨麻豆精品| 午夜福利乱码中文字幕| 日韩伦理黄色片| 国产精品一区www在线观看| 国产一区二区激情短视频 | 国产 精品1| 国产免费现黄频在线看| 国产高清国产精品国产三级| av线在线观看网站| 亚洲av中文av极速乱| 亚洲婷婷狠狠爱综合网| 国产又色又爽无遮挡免| 草草在线视频免费看| 熟女人妻精品中文字幕| 亚洲成色77777| 亚洲精品国产色婷婷电影| 自线自在国产av| 这个男人来自地球电影免费观看 | 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| av国产久精品久网站免费入址| 宅男免费午夜| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩另类电影网站| 人体艺术视频欧美日本| 综合色丁香网| 日韩精品免费视频一区二区三区 | 国产亚洲最大av| av黄色大香蕉| 一区在线观看完整版| 在线免费观看不下载黄p国产| 伦理电影大哥的女人| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 亚洲精品久久久久久婷婷小说| 一本大道久久a久久精品| 久久久久国产网址| 国产视频首页在线观看| 十八禁高潮呻吟视频| av播播在线观看一区| av女优亚洲男人天堂| 亚洲成色77777| 久久人人爽av亚洲精品天堂| 乱码一卡2卡4卡精品| 婷婷色综合大香蕉| 亚洲国产色片| 最近中文字幕2019免费版| 天美传媒精品一区二区| 国产午夜精品一二区理论片| 亚洲精品色激情综合| 自线自在国产av| 久久99一区二区三区| 性色avwww在线观看| 伦理电影大哥的女人| 成人18禁高潮啪啪吃奶动态图| 国产xxxxx性猛交| 亚洲伊人久久精品综合| 成年av动漫网址| 波多野结衣一区麻豆| 性色av一级| 国产免费现黄频在线看| 国产综合精华液| 欧美老熟妇乱子伦牲交| 欧美国产精品一级二级三级| 久久久久久久亚洲中文字幕| 成人国产麻豆网| 亚洲国产欧美日韩在线播放| av网站免费在线观看视频| 欧美日韩国产mv在线观看视频| 亚洲在久久综合| 欧美精品av麻豆av| 成年女人在线观看亚洲视频| 99久久中文字幕三级久久日本| 女人被躁到高潮嗷嗷叫费观| 在线观看免费日韩欧美大片| 九色亚洲精品在线播放| 国产熟女欧美一区二区| 亚洲,一卡二卡三卡| 国产成人一区二区在线| 少妇的丰满在线观看| 中文字幕免费在线视频6| 亚洲精品,欧美精品| 人妻少妇偷人精品九色| 一级毛片黄色毛片免费观看视频| 91精品国产国语对白视频| 欧美激情 高清一区二区三区| 日韩av免费高清视频| 国产日韩欧美亚洲二区| 亚洲国产av影院在线观看| 香蕉精品网在线| 欧美亚洲 丝袜 人妻 在线| 午夜激情av网站| 亚洲人成77777在线视频| av在线播放精品| 欧美+日韩+精品| 精品一区在线观看国产| 日韩视频在线欧美| 国产成人91sexporn| 久久人人97超碰香蕉20202| 久久久久精品久久久久真实原创| 丝袜在线中文字幕| av在线观看视频网站免费| 丰满饥渴人妻一区二区三| 国产精品一国产av| 国产精品久久久久久精品古装| 男人操女人黄网站| 飞空精品影院首页| av免费观看日本| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 欧美丝袜亚洲另类| 两性夫妻黄色片 | 久久久国产欧美日韩av| 18禁裸乳无遮挡动漫免费视频|