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

    Ce-Doped Smart Adsorbents with Photoresponsive Molecular Switches for Selective Adsorption and Efficient Desorption

    2020-09-12 03:11:18PengTanYaoJiangShiChaoQiXiaJunGaoXiaoQinLiuLinBingSun
    Engineering 2020年5期

    Peng Tan, Yao Jiang, Shi-Chao Qi, Xia-Jun Gao, Xiao-Qin Liu*, Lin-Bing Sun*

    State Key Laboratory of Materials-Oriented Chemical Engineering, Jiangsu National Synergetic Innovation Center for Advanced Materials, College of Chemical Engineering, Nanjing Tech University, Nanjing 211816, China

    Keywords:Azobenzene Ce-doped mesoporous silica Molecular switches Selective adsorption Efficient desorption

    A B S T R A C T Achieving efficient adsorption and desorption processes by controllably tuning the properties of adsorbents at different technical stages is extremely attractive. However, it is difficult for traditional adsorbents to reach the target because of their fixed active sites. Herein, we report on the fabrication of a smart adsorbent, which was achieved by introducing photoresponsive azobenzene derivatives with cis/trans isomers to Ce-doped mesoporous silica. These photoresponsive groups serve as ‘‘molecular switches”by sheltering and exposing active sites,leading to efficient adsorption and desorption.Ce is also doped to provide additional active sites in order to enhance the adsorption performance. The results show that the cis isomers effectively shelter the active sites, leading to the selective adsorption of methylene blue(MB)over brilliant blue(BB),while the trans isomers completely expose the active sites,resulting in the convenient release of the adsorbates. Both selective adsorption and efficient desorption can be realized controllably by these smart adsorbents through photostimulation. Moreover, the performance of the obtained materials is well maintained after five cycles.

    1. Introduction

    Adsorptive separation utilizing solid adsorbents has attracted increasing attention in the past decades due to the low equipment investment required, mild operating conditions, and efficient utilization of energy[1-4].During adsorption processes,it is essential to be able to selectively remove the target substances from mixtures and to efficiently release the adsorbates for regeneration,since doing so is strongly associated with product quality and with the input time and energy [5]. As adsorbents play a critical role in adsorption processes, achieving efficient adsorption and desorption processes by controllably tuning the properties of adsorbents at different technical stages is extremely attractive[6-9].However,it is difficult for traditional adsorbents to reach the target due to their fixed active sites.To be specific,the active sites of traditional adsorbents can be easily accessed by all adsorbates when the pore sizes are sufficiently large.This results in low selectivity,such that the product purity cannot meet the increasing industrial requirements. Modulating the pore sizes or constructing target-specific sites in traditional adsorbents can promote selectivity, but results in a problem: The release of adsorbates becomes difficult. This problem is caused by increased diffusion resistance or enhanced adsorbate-adsorbent interactions,and inevitably gives rise to difficulty in desorption and intensive energy consumption.The tradeoff between selective adsorption and efficient desorption poses enormous challenges when devising and synthesizing desirable adsorbents [10]. As adsorption and desorption are sequential stages in adsorption processes, the technical efficiency can be remarkably improved if adsorbents can be made to be responsive to external stimuli and to shelter or expose their active sites on demand. This concept is very different from the traditional concept, in which variation of the adsorption capacity depends on the cyclic variation of pressure and/or temperature, which requires massive energy consumption.

    Inspired by living systems that can adjust themselves to changeable environments for survival, functional materials that are responsive to external stimuli with structurally or chemically switchable properties have been engineered to perform various roles for practical application in recent years[11-15].The external stimuli include light [16-19], temperature [20-23], pH [24-28],and redox reactions[29-32].Among these diverse forms of stimulation, light has been given particular attention because it permits fast and precise control and generates few undesired byproducts.As typical photoresponsive molecules, azobenzenes can reversibly transform between cis and trans isomerization through stimulation with ultraviolet (UV) light and visible (Vis) light [33-35]. This isomerization causes azobenzene and its derivatives to change conformation between the rod-like and curved states. Based on this conformational change,azobenzenes are promising candidates for sheltering and exposing active sites to adsorbates in a manner that is well tuned by UV/Vis irradiation on demand [36-39]. It is worth noting that this target requires suitable pore sizes in the supports.Macroporous/mesoporous materials have relatively large pore sizes; thus, both large and small molecules can pass through the pore channels and come into contact with active sites,resulting in the failure of selective adsorption [40-42]. Although microporous materials with relatively small pore sizes can hinder the entrance of large molecules, it is difficult for the adsorbed substances to be released after adsorption [43,44]. Consequently,mesoporous materials with pore sizes of around 2-5 nm are good candidates for capturing both small and large molecules. The rational introduction of photoresponsive molecules can endow mesoporous materials with tunable selectivity and decrease the energy consumption required for desorption.

    In this work, molecular switches—namely, 4-(3-triethoxysilyl propyl-ureido)azobenzene (AB-TPI)—were introduced into mesopores during the synthesis of Ce-doped mobile composition of matters (MCM)-41 (denoted herein as CeM). Cerium (Ce) was doped in order to generate additional active sites to enhance the adsorption performance [45-48]. The relatively large molecule methylene blue (MB) and the relatively small molecule brilliant blue (BB) were employed to test the adsorption performance. The results show that Ce doping leads to a remarkable increase of 13.5% and 65.5%, respectively, in the capacity of MCM-41 to separately capture MB and BB. AB-TPI molecules act as molecular switches that alternately shelter and expose active sites through photoisomerization. To be specific, the cis isomers with curved states are close to the active sites in space, and thus have steric effects on relatively large molecules while negligibly affecting the approach of relatively small molecules (Fig. 1(a)). When AB-TPI transforms into the trans configuration, it fully exposes the active sites, allowing the attached adsorbates to be easily released, and achieving efficient desorption (Fig. 1(b)). As a result, these molecular switches can reversibly regulate adsorption and desorption. Triggered by UV/Vis light, such smart adsorbents display excellent performance in selective adsorption and efficient desorption, and thus provide clues for the development of controllable adsorbents for energy efficiency.

    2. Experimental section

    2.1. Chemicals

    3-(Triethoxysilyl)propyl isocyanate(TESPIC),cetyltrimethylammonium bromide (CTAB), tetraethylorthosilicate (TEOS), and Ce(SO4)2·4H2O were purchased from Sigma-Aldrich Co., USA. 4-Phenylazoaniline(PAA)was produced by Tokyo Chemical Industry(TCI)Co.,Ltd.,Japan.Methanol,hexane,toluene,and tetrahydrofuran (THF) were purchased from Aladdin, China. The trace amount of water in the toluene and THF was removed by a 4A molecular sieve before using;the other chemicals were used without further purification.

    2.2. Materials synthesis

    The photoresponsive molecule AB-TPI was prepared using a procedure similar to that reported in the literature(Supplementary data Fig. S1) [49]. 1.58 g of PAA and 2.05 g of TESPIC were first mixed in 12 mL of anhydrous THF. The mixture was then transferred to a Pyrex tube under a nitrogen (N2) atmosphere. After heating at 65 °C for 12 h, 40 mL of hexane was added to the mixture. The parent materials were then incubated overnight at 20 °C. Orange crystals were obtained after washing the solid products three times with hexane.

    MCM-41 was prepared based on a previous method[50].3.5 mL of NaOH and 1.00 g of CTAB were mixed in 480 mL of water, and the solution was heated to 80 °C. Next, TEOS (5 mL) was added dropwise to the solution with vigorous stirring. After reacting for 2 h, the white silica colloidal suspension was filtered, and the obtained solids were washed with water several times. The obtained solids were then dispersed in ethanol (100 mL) and HCl(1 mol·L-1, 1 mL), followed by refluxing at 80 °C for 6 h. Finally,the extracted template was removed by washing the products with ethanol several times.The process of preparing CeM was similar to the process used to prepare MCM-41, but 0.18 g of Ce(SO4)2·4H2O was initially added to the mixture.

    Fig. 1. (a) Illustration of the selective adsorption of smart adsorbents with AB-TPI in the cis configuration; (b) efficient desorption with AB-TPI in the trans configuration;(c) cis/trans isomerization of AB-TPI by UV/Vis light irradiation.

    The smart azobenzene-functionalized adsorbents were prepared by in situ synthesis. 3.5 mL of NaOH with a concentration of 2 mol·L-1and 1.00 g of CTAB were mixed in 480 mL of water.Next, TEOS (5 mL) and a certain amount of AB-TPI (see below)was added dropwise to the solution with vigorous stirring. The subsequent steps were the same as those used with MCM-41.The obtained samples were denoted as Azo-CeM-1, Azo-CeM-2,Azo-CeM-3, and Azo-CeM-4, corresponding to the added AP-TPI amounts of 34.7, 69.4, 138.7, and 218.0 mg, respectively.

    2.3. Materials characterization

    1H and13C nuclear magnetic resonance(NMR)was performed on a Bruker AVANCE NEO spectrometer(400 MHz,Bruker Corporation,USA) in dimethyl sulfoxide. X-ray diffraction (XRD) patterns were collected by a Rigaku D/MAX-γA with CuKα radiation(Rigaku Corporation, Japan). Scanning electron microscopy (SEM) was performed on a Hitachi New Generation SU8010 electron microscope(Hitachi, Ltd., Japan). The metal content of different samples was determined using energy-dispersive X-ray (EDX) analyses with Kevex equipment attached to the scanning microscope. Transmission electron microscope (TEM) images were collected using an FEI Tecnai G2 F20 electron microscope(FEI Company,Japan)operated at an acceleration voltage of 200 kV.N2adsorption-desorption tests were performed on an ASAP 2020 analyzer(GAT Scientific Sdn Bhd,Malaysia)at-196°C,and samples were pretreated at 100°C in a high vacuum for 4 h.The Brunauer-Emmett-Teller(BET)surface area was obtained according to the adsorption data, with relative pressure(p/p0)ranging from 0.04 to 0.20.The Barrett-Joyner-Hale nda(BJH)method was employed to calculate the pore size distributions, based on the adsorption branches of the isotherms. UV-Vis spectra were recorded on a PerkinElmer Lambda 35 spectrometer(PerkinElmer Inc., USA) with wavelengths from 300 to 600 nm.Fourier-transform infrared (IR) spectra were obtained with wavelengths from 550 to 4000 nm(Nicolet Nexus 470,Thermo Fisher Scientific Inc.,USA).The content of the elements carbon(C),hydrogen(H), and nitrogen(N) was obtained by an elemental analyzer from Elementar Analysensysteme GmbH, Germany. Thermogravimetric(TG)analysis was carried out on an STA-499C(NETZSCH,Germany)with a heating rate of 20°C·min-1.

    2.4. Adsorption experiments

    For the selective adsorption study,the azobenzene groups of the adsorbents remained as cis isomers; thus, the active sites were sheltered from the relatively large molecules. Two dye molecules,MB and BB,were chosen as the adsorbates.In a typical experiment,one of the adsorbates (3 mL, 25 mg·L-1) was added in a cuvette containing 2 mg of adsorbent. The system was treated with UV light(365 nm)for 3 h.The amount adsorbed(Qe)was determined according to Eq. (1).

    where ciis the initial concentration,ceis the equilibrium concentration, V is the volume of water, and m is the mass of adsorbent.

    For the desorption study, the azobenzene groups of the adsorbent were transformed from cis to trans isomers; thus, the active sites were exposed. In a typical experiment, the saturated adsorbent was irradiated with Vis light with a wavelength of 450 nm for 180 min. An ethanol/water mixture was used to wash the adsorbent.The concentration of the desorbed adsorbate was determined at regular intervals. The amount desorbed (Qd) was calculated based on Eq. (2):

    where cdrefers to the equilibrium concentration during desorption.

    3. Results and discussion

    3.1. Textural and structural properties

    AB-TPI was first analyzed by1H and13C NMR (Figs. S2(a) and(b)). The characteristic resonances clearly showed the successful synthesis of AB-TPI [19]. No other undefined resonances were observed, suggesting few impurities in the products. After removing the template, the photoresponsive groups were modulated to shelter and expose the active sites. The grafted amounts of the photoresponsive groups on Azo-CeM-1, Azo-CeM-2, Azo-CeM-3,and Azo-CeM-4 were calculated as 5.6 wt%, 7.8 wt%, 12.6 wt%,and 17.8 wt%, respectively (Table 1).

    The structural information of the materials was analyzed by means of XRD patterns.Fig.2(a)shows the low-angle XRD pattern of CeM.The prominent peak at the 2θ angle of 2.4 can be assigned to the reflection of(100)interplanar spacing(d100).This represents the two-dimensional (2D) hexagonal pore symmetry of MCM-41(Fig. S3(a)) [51,52]. The diffraction peaks of Azo-CeM-1, Azo-CeM-2, Azo-CeM-3, and Azo-CeM-4 are the same as those of CeM, indicating that the azobenzene groups have little effect on the structure of CeM.The wide-angle XRD patterns of the different samples were also recorded. MCM-41 displays only one broad diffraction peak centered around 23° (Fig. S3(b)) [53-55]. After Ce is introduced, three diffraction peaks at the 2θ angle of 28°,47°, and 56° can be observed (Fig. 2(b)). These peaks can be ascribed to the d111, d220, and d311reflections of a typical CeO2phase (Joint Committee on Powder Diffraction Standards (JCPDS)no. 034-0394) [56]. It should be stated that these new peaks are weak and broad, implying that the CeO2is well dispersed.

    Table 1 Physicochemical parameters and elemental composition of different samples.

    Fig. 2. (a) Low-angle and (b) wide-angle XRD patterns of different samples.

    N2adsorption-desorption tests were conducted to study the pore structures.Figs. 3(a) and S4(a) show that all the samples displayed typical type-IV isotherms. A remarkable increase can be observed at relative pressure from 0.3 to 0.4, indicating the abundant mesopores of different samples. MCM-41 has a BET surface area of 1224 m2·g-1, suggesting that it has highly developed open pores (Table S1). After the introduction of Ce, this value decreases to 971 m2·g-1(Table 1). The introduction of photoresponsive groups further decreases the surface area,although it is still maintained at a high level.For example,the surface area of Azo-CeM-4,which has the greatest amount of grafting,is as high as 825 m2·g-1.Figs.3(b)and S4(b)show that the pore size distributions of different samples are narrow. The apertures of most mesopores are around 2.9 nm, which is suitable for the stretching and bending of the azobenzene groups. In addition, the mesoporous structure can be well maintained after Ce doping and the AB-TPI grafting reaction. However, compared with MCM-41, CeM exhibits a decrease of pore volume, and the introduction of AB-TPI also results in a declining trend in pore volume with increasing amounts of AB-TPI. This finding confirms that the active sites and AB-TPI were introduced into the pore channels of MCM-41. Azo-CeM-3 presented identical N2adsorption behavior in the trans and cis configurations (Fig. S5). This is because the size of an N2molecule is relatively small.

    TEM and SEM images were employed to observe the morphology of the adsorbents. As shown in the SEM images (Fig. S6), the MCM-41 particles were regular spheres. CeM had a suborbicular morphology with relatively large particles, and some pieces were packed on the surface. This is probably due to the slight aggregation of particles caused by Ce doping. No further aggregation was observed for Azo-CeM-3, indicating that the AB-TPI grafting had little impact on the morphology. The TEM images clearly show the pore channels of Azo-CeM-3 (Fig. 4(a)). The uniform and straight pore channels of Azo-CeM-3 are important for the capture of adsorbates. This observation supports the XRD and N2adsorption results,which indicate that the pore structures were well preserved during synthesis, despite some aggregation of adsorbent particles.

    The identification of AB-TPI and the thermal stability of the adsorbents were demonstrated by the following techniques. As shown by the EDX analysis (Fig. S7), Azo-CeM-3 was composed of the elements C,N,oxygen(O),silicon(Si),and Ce;the EDX mapping (Fig. 4(b)) clearly shows that these elements were uniformly dispersed in the frameworks. The main organic groups and covalent bonds of the adsorbents were measured by means of IR spectra(Figs. 5(a) and S8(a)). For MCM-41, the band at 957 cm-1can be ascribed to the stretching vibrations of the silanol groups. The bands at 1086 and 797 cm-1can be respectively assigned to the asymmetric and symmetric stretching vibrations of the Si-O-Si frameworks [57,58]. For the other samples, the characteristic bands of MCM-41 can be clearly observed. After introducing ABTPI, three new bands appear. The band at 1548 cm-1can be ascribed to the stretching vibration of -NH-CO-NH-, the band at 1502 cm-1corresponds to the stretching vibration of -NH-, and the band at 688 cm-1can be assigned to the vibration of the aromatic amine[49].These results confirm the successful introduction of AB-TPI; furthermore, the intensity of these characteristic bands increases as more AB-TPI groups are introduced. The weak vibrational band at around 2950 cm-1is attributed to the vibration of the carbon chains of CTAB, and derives from its incomplete removal by extraction.TG analysis was performed to test the thermal stability of different samples(Figs.5(b)and S8(b)).The results show that AB-TPI decomposed in the temperature interval from 200 to 700 °C, and the residual template CTAB decomposed in the temperature interval from 500 to 600 °C.

    Fig. 3. (a) N2 adsorption-desorption isotherms and (b) pore size distributions of different samples. The curves are offset on the y axis for clarity.

    Fig. 4. (a) TEM image of Azo-CeM-3 and (b) its EDX mapping images of the C, N, O, Si, and Ce elements.

    Fig. 5. (a) IR spectra and (b) TG curves of different samples.

    3.2. Photoresponsive properties

    Fig. 6. Alteration in the UV-Vis spectra of AB-TPI upon UV/Vis light irradiation.

    The photoisomerization of AB-TPI (dissolved in ethanol) was recorded by means of UV/Vis spectroscopy. As shown in Fig. 6,two bands around 350 and 450 nm were observed, which correspond to the π-to-π* and n-to-π* transitions of the π-conjugated system for AB-TPI. UV irradiation resulted in an intensity decline of the broad band at 360 nm,while leading to an intensity increase of the broad band at 450 nm.This result corresponds to the isomerization of AB-TPI from the trans to cis isomer.When the irradiation time reached 60 s,AB-TPI reaches a photostationary state.UV irradiation recovered the adsorption intensity of the two bands,which corresponds to the conformational change of AB-TPI from the cis to the trans configuration. The intensity of the 350 nm band became stronger than that of the pure AB-TPI when the irradiation time was longer than 90 s. This result indicates that some azobenzene groups initially exist as cis isomers. Through cyclic variations of UV and Vis, the AB-TPI is able to switch between the trans and cis configurations (Figs. S9(a) and (b)), reversibly acting as molecule switches modulated by photostimulation. After Ce and ABTPI are introduced into the pore channels of MCM-41, photoisomerization of the azobenzene groups can be achieved. Taking Azo-CeM-3 as an example, its UV-Vis spectrum (Figs. 7(a) and(b))exhibited the same change as pure AB-TPI upon UV or Vis irradiation. In addition, the photoisomerization was reversible after three cycles. Thus, it is confirmed that the photoresponsiveness of the azobenzene groups was well maintained after immobilization on MCM-41. On this basis, switching of the active sites by means of AB-TPI is likely, and may have significant effects on adsorption and desorption.

    3.3. Adsorption performance

    Two probing molecules were employed to demonstrate the promoting effects of doping with Ce and the influence of the photoresponsive groups on adsorption. MB was used as a relatively small molecule,with a size of 1.26 nm,while BB was used as a relatively large molecule, with a size of 1.98 nm (Fig. S10). The effects of Ce doping on adsorption were investigated first. Compared with MCM-41 (Fig. S11(a)), CeM exhibited an obviously larger adsorption capacity for both MB and BB(Figs.S12(a)and S13(a)),indicating that Ce doping resulted in an increased number of active sites.As MB is a weak Lewis base [59], it can be captured through acidbase interactions between the dye molecules and the active sites on the adsorbents’ surfaces. Ce doping increases the electron-pair donation to the inner-side pore wall surface, forming Lewis acid centers that promote the Lewis acidity of the adsorbents [60].For the capture of BB, which is an anionic dye, Ce doping can promote cationic-anionic interactions between BB molecules and the modified surfaces.The positively charged surface of the adsorbents contributes to the electrostatic attraction between the adsorbates and adsorbents [61-63]. Consequently, Ce doping provided additional active sites for capturing both MB and BB.

    Fig.7. (a)Alteration in the UV-Vis spectra of Azo-CeM-3 upon UV/Vis light irradiation;(b)reversible changes in absorbance at 353 nm(A353)as a function of cycles for Azo-CeM-3.

    Next, the influence of the photoresponsive groups was investigated. Upon UV light irradiation, CeM displayed an inconspicuous change in the amount adsorbed of both MB and BB,which indicates that external photostimulation cannot interfere with the adsorption performance of CeM. After sheltering the active sites, the capacities of the photoresponsive samples exhibited insignificant differences for MB (Figs. S12(b), (c), and (d)), but showed remarkable declines for BB(Figs.S13(b),(c),and(d)).Taking Azo-CeM-3 as an example, its capacity for MB slightly decreased from 22.4 to 20.5 mg·g-1(Fig.8(a)).Interestingly,an obvious decrease in capacity of around 43%(from 13.8 to 7.8 mg·g-1,Fig.8(b))was observed for BB after UV light irradiation.This result suggests that Azo-CeM-3 has excellent selectivity for MB over BB after photostimulation.We propose that steric hindrance derived from the photoisomerization of AB-TPI in MCM-41 explains this phenomenon. After UV light irradiation,the active sites in the pore channels of the adsorbents are switched off by the AB-TPI in the cis configuration.In this case, the active sites are sheltered by the curved structure of the photoresponsive groups, making it difficult for the relatively large BB molecules to access them. In contrast, the relatively small MB molecules can pass freely through the confined spaces and be adsorbed by the active sites.

    The loading amount of AB-TPI was further modulated, and the corresponding effects on adsorption were investigated. As shown in Fig. 9, although more AB-TPI groups were introduced to CeM,a change in isotherms for MB was not obvious before and after UV light irradiation.This result confirms that MB molecules cannot be excluded by photoresponsive groups. It was notable that BB capture became difficult after UV light irradiation,and the adsorption capacity decreased as the loading amount of AB-TPI increased.It should be stated that when excessive AB-TPI groups are loaded,the adsorption capacity is adversely affected because the ordering of MCM-41 decreases and the pore channels can be blocked. The optimal sample was Azo-CeM-3, with a variation of 43% in capacity. This finding indicates that the variation in the spatial configuration of AB-TPI greatly interferes—through steric hindrance—with the approach of BB molecules to the active sites,thus providing an efficient and energy-saving method to vary the adsorption capacity.

    The desorption performance was explored on Azo-CeM-3 saturated by MB.As shown in Fig.10,the obtained desorption amount reached as high as 97.4% with AB-TPI in the trans configuration.However, only 68.3% of the adsorbates were released with AB-TPI in the cis configuration. This finding indicates that the exposure of active sites contributes to desorption when AB-TPI is transformed to the trans configuration. Finally, the reusability of Azo-CeM-3 was studied. As shown in Fig. 11, there was insignificant loss in capacity after five cycles,suggesting excellent recyclability.

    Fig. 8. Dynamic adsorption curves of (a) MB and (b) BB on Azo-CeM-3 before and after UV light irradiation.

    Fig. 9. Change in the adsorption capacity on MB and BB over different adsorbents before and after UV light irradiation.

    Fig. 10. Desorption curves of MB on Azo-CeM-3 before and after UV light irradiation.

    Fig. 11. Adsorption/desorption cycles of MB on Azo-CeM-3 under UV/Vis light irradiation. Ads: adsorption; Des: desorption.

    4. Conclusions

    In summary, Ce-doped smart adsorbents with photoresponsive molecular switches were successfully fabricated. Ce was doped to provide additional active sites for enhancing adsorption performance.Through the precise sheltering and exposing of active sites by photoresponsive groups, the smart adsorbents exhibited excellent performance in selective adsorption and efficient desorption.Moreover,the adsorption and desorption processes could be regulated reversibly by these molecular switches. These smart adsorbents displayed a remarkable performance in selective adsorption and efficient desorption, giving clues for the development of controllable adsorbents for energy efficiency.Furthermore,this new conception with the introduction of stimuli-responsive molecules should be transferrable to other adsorbents, thus extending the application of light energy in adsorption to protect the environment and increase energy efficiency.

    Acknowledgements

    This work was supported by the National Science Fund for Excellent Young Scholars(21722606),the National Natural Science Foundation of China (21676138, 21878149, 21808110, and 21576137), the China Postdoctoral Science Foundation(2018M632295), and the Six Talent Plan (2016XCL031).

    Compliance with ethics guidelines

    Peng Tan,Yao Jiang,Shi-Chao Qi,Xia-Jun Gao,Xiao-Qin Liu,and Lin-Bing Sun declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2020.03.005.

    av福利片在线观看| 99久久精品一区二区三区| 最后的刺客免费高清国语| 午夜亚洲福利在线播放| 免费观看人在逋| 成人漫画全彩无遮挡| 亚洲国产欧美人成| 久久精品综合一区二区三区| 亚洲精品粉嫩美女一区| 波多野结衣高清作品| 丰满乱子伦码专区| 在线观看美女被高潮喷水网站| 国产淫片久久久久久久久| 在线播放国产精品三级| 免费av不卡在线播放| 国产精品无大码| 床上黄色一级片| 一边摸一边抽搐一进一小说| 色尼玛亚洲综合影院| 色5月婷婷丁香| 12—13女人毛片做爰片一| 久久中文看片网| 91久久精品国产一区二区成人| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 欧美日本视频| 男女啪啪激烈高潮av片| 精品久久国产蜜桃| 最近的中文字幕免费完整| 高清在线视频一区二区三区 | 麻豆国产97在线/欧美| 久久久精品欧美日韩精品| 午夜爱爱视频在线播放| 大又大粗又爽又黄少妇毛片口| 一级黄色大片毛片| 久久婷婷人人爽人人干人人爱| 搡老妇女老女人老熟妇| 久久草成人影院| 午夜免费男女啪啪视频观看| 久久欧美精品欧美久久欧美| 激情 狠狠 欧美| 成人av在线播放网站| 老女人水多毛片| 黄色一级大片看看| 久久99蜜桃精品久久| 欧美一区二区国产精品久久精品| 精品无人区乱码1区二区| 国产精品野战在线观看| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 九九爱精品视频在线观看| 天天一区二区日本电影三级| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 青春草亚洲视频在线观看| 午夜激情福利司机影院| 国产伦在线观看视频一区| 久久99精品国语久久久| а√天堂www在线а√下载| 性欧美人与动物交配| 欧美日韩国产亚洲二区| 又爽又黄a免费视频| 国产精品国产高清国产av| 26uuu在线亚洲综合色| 亚洲av熟女| 晚上一个人看的免费电影| 久久久久国产网址| 色噜噜av男人的天堂激情| av视频在线观看入口| 桃色一区二区三区在线观看| 精品少妇黑人巨大在线播放 | 美女xxoo啪啪120秒动态图| 久久99热6这里只有精品| 在线播放无遮挡| 国产av在哪里看| 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 亚洲性久久影院| 国产精品无大码| 亚洲欧美日韩卡通动漫| 直男gayav资源| 三级男女做爰猛烈吃奶摸视频| 色综合亚洲欧美另类图片| 国产精品久久视频播放| 亚洲丝袜综合中文字幕| 国产亚洲5aaaaa淫片| 久久综合国产亚洲精品| 国产久久久一区二区三区| 看免费成人av毛片| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 国产成人精品一,二区 | 亚洲第一电影网av| 一进一出抽搐gif免费好疼| 色哟哟·www| 久久人人爽人人片av| 欧美日韩国产亚洲二区| 久久精品国产亚洲av香蕉五月| 中国美女看黄片| 日本一本二区三区精品| 亚洲精品久久久久久婷婷小说 | 亚洲精品久久国产高清桃花| 国产欧美日韩精品一区二区| 少妇人妻一区二区三区视频| 日日摸夜夜添夜夜爱| 日本在线视频免费播放| 男女啪啪激烈高潮av片| 免费搜索国产男女视频| 国产日韩欧美在线精品| 又粗又爽又猛毛片免费看| 18禁在线播放成人免费| 国产91av在线免费观看| 国产极品精品免费视频能看的| 国产国拍精品亚洲av在线观看| 成人性生交大片免费视频hd| 看免费成人av毛片| АⅤ资源中文在线天堂| 桃色一区二区三区在线观看| 最近2019中文字幕mv第一页| 中国美白少妇内射xxxbb| 深夜a级毛片| 搡老妇女老女人老熟妇| 级片在线观看| 神马国产精品三级电影在线观看| 99久久无色码亚洲精品果冻| 精品久久久噜噜| 大型黄色视频在线免费观看| 91久久精品国产一区二区成人| 精品人妻一区二区三区麻豆| 国内精品美女久久久久久| 国国产精品蜜臀av免费| 久久精品国产99精品国产亚洲性色| 大香蕉久久网| 久久人人爽人人爽人人片va| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 欧美性猛交黑人性爽| 久久精品久久久久久久性| 岛国毛片在线播放| 国产精品人妻久久久影院| 日韩 亚洲 欧美在线| 国产伦在线观看视频一区| 最新中文字幕久久久久| 一本久久精品| 久久精品国产亚洲av涩爱 | 精品国内亚洲2022精品成人| 亚洲精品成人久久久久久| 精品久久久久久久末码| 一级毛片aaaaaa免费看小| 国产av在哪里看| 色播亚洲综合网| 久久久a久久爽久久v久久| 91午夜精品亚洲一区二区三区| 黄色视频,在线免费观看| 女的被弄到高潮叫床怎么办| 两个人的视频大全免费| 搡老妇女老女人老熟妇| 最近最新中文字幕大全电影3| 日韩一本色道免费dvd| 午夜亚洲福利在线播放| 99精品在免费线老司机午夜| 国内揄拍国产精品人妻在线| 一进一出抽搐动态| 免费av不卡在线播放| 色哟哟·www| 中文字幕制服av| 日本-黄色视频高清免费观看| 少妇的逼好多水| 性欧美人与动物交配| 亚洲成人av在线免费| 国产av不卡久久| 麻豆av噜噜一区二区三区| 国产黄片美女视频| 一级二级三级毛片免费看| 国产91av在线免费观看| 国产极品天堂在线| 一级毛片我不卡| 免费一级毛片在线播放高清视频| 成人国产麻豆网| 精品欧美国产一区二区三| 日日摸夜夜添夜夜添av毛片| 亚州av有码| 久久精品国产亚洲网站| 一级黄色大片毛片| 校园人妻丝袜中文字幕| av在线老鸭窝| 精品国产三级普通话版| 国产中年淑女户外野战色| 九草在线视频观看| or卡值多少钱| 成人综合一区亚洲| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类| av在线蜜桃| 久久综合国产亚洲精品| 又黄又爽又刺激的免费视频.| 此物有八面人人有两片| 久久久久九九精品影院| 国产毛片a区久久久久| 国产av不卡久久| 如何舔出高潮| 国产乱人视频| 嫩草影院入口| 亚洲美女搞黄在线观看| 国产三级在线视频| av女优亚洲男人天堂| 内射极品少妇av片p| 国产精品人妻久久久影院| 一级黄片播放器| 国产精品人妻久久久久久| 日本免费一区二区三区高清不卡| 国产精品久久久久久精品电影| 国产精品三级大全| 亚洲欧美中文字幕日韩二区| 国产高清三级在线| 在线天堂最新版资源| 久久久久久伊人网av| 性插视频无遮挡在线免费观看| 九九爱精品视频在线观看| 亚洲成人av在线免费| 亚洲av男天堂| 一本久久精品| 久久久精品大字幕| 一区二区三区高清视频在线| 又粗又硬又长又爽又黄的视频 | 赤兔流量卡办理| 亚洲精品乱码久久久久久按摩| 一级毛片久久久久久久久女| 尤物成人国产欧美一区二区三区| 日韩中字成人| 干丝袜人妻中文字幕| 九九在线视频观看精品| 久久韩国三级中文字幕| 欧美人与善性xxx| 亚洲18禁久久av| 91精品国产九色| 一进一出抽搐gif免费好疼| 卡戴珊不雅视频在线播放| 看黄色毛片网站| 99久久精品热视频| 人妻系列 视频| 99视频精品全部免费 在线| 99九九线精品视频在线观看视频| 国产精品电影一区二区三区| 久久中文看片网| 免费看日本二区| 国产免费一级a男人的天堂| 岛国在线免费视频观看| 日韩成人伦理影院| a级毛片免费高清观看在线播放| 国内久久婷婷六月综合欲色啪| 亚洲精品日韩在线中文字幕 | 男人舔奶头视频| 一边摸一边抽搐一进一小说| 国产黄片美女视频| 久久6这里有精品| 久久久欧美国产精品| 99热网站在线观看| 成熟少妇高潮喷水视频| 久久99蜜桃精品久久| 久久久国产成人免费| 国产精品1区2区在线观看.| 乱人视频在线观看| 可以在线观看毛片的网站| www.色视频.com| 欧美3d第一页| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 91久久精品电影网| 亚洲av免费高清在线观看| 国产毛片a区久久久久| 国产视频首页在线观看| 成人一区二区视频在线观看| 亚洲av免费在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日韩人妻高清精品专区| 成人av在线播放网站| 熟女电影av网| 18禁裸乳无遮挡免费网站照片| 国产免费男女视频| 青青草视频在线视频观看| 久久久精品大字幕| 国产成人91sexporn| 国产单亲对白刺激| 2021天堂中文幕一二区在线观| 精品欧美国产一区二区三| av天堂中文字幕网| 免费搜索国产男女视频| 国产精品.久久久| 99九九线精品视频在线观看视频| kizo精华| 中国国产av一级| av在线老鸭窝| 日韩制服骚丝袜av| 内地一区二区视频在线| 男女下面进入的视频免费午夜| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 久久99热6这里只有精品| 日本免费一区二区三区高清不卡| 久99久视频精品免费| 高清午夜精品一区二区三区 | 国产女主播在线喷水免费视频网站 | 在线观看66精品国产| а√天堂www在线а√下载| 久久亚洲国产成人精品v| 少妇的逼好多水| 国产高清三级在线| 亚洲精品色激情综合| 在线播放国产精品三级| 久久久久久伊人网av| 精品久久久久久成人av| 国内揄拍国产精品人妻在线| 欧美高清性xxxxhd video| 美女大奶头视频| 99热全是精品| 日日啪夜夜撸| 久久精品久久久久久久性| 亚洲av中文av极速乱| 国产精品一区www在线观看| 欧美日韩一区二区视频在线观看视频在线 | 又爽又黄无遮挡网站| 亚洲精品自拍成人| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 亚洲欧美日韩高清专用| 国内精品美女久久久久久| 亚洲精品国产av成人精品| 三级经典国产精品| 午夜福利在线观看吧| 99久久九九国产精品国产免费| 青青草视频在线视频观看| 国产精品不卡视频一区二区| 嫩草影院精品99| 国产激情偷乱视频一区二区| 亚洲欧洲日产国产| 人人妻人人看人人澡| 床上黄色一级片| 青春草视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最好的美女福利视频网| 99热这里只有是精品在线观看| 大型黄色视频在线免费观看| 少妇熟女欧美另类| 国产av在哪里看| 久久久久性生活片| 久久久午夜欧美精品| 美女 人体艺术 gogo| 91午夜精品亚洲一区二区三区| 久99久视频精品免费| 久久热精品热| 精品无人区乱码1区二区| 久久精品人妻少妇| a级毛色黄片| 亚洲精品日韩在线中文字幕 | 亚洲美女视频黄频| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 精品人妻一区二区三区麻豆| 日韩三级伦理在线观看| 日韩视频在线欧美| 色综合站精品国产| 日韩在线高清观看一区二区三区| 偷拍熟女少妇极品色| 国产成人一区二区在线| av视频在线观看入口| 九草在线视频观看| 欧美xxxx黑人xx丫x性爽| 网址你懂的国产日韩在线| 丝袜美腿在线中文| 国产一区亚洲一区在线观看| 中文字幕久久专区| 丰满乱子伦码专区| 久久这里有精品视频免费| 男女边吃奶边做爰视频| 成人欧美大片| 九九热线精品视视频播放| 亚洲av中文字字幕乱码综合| 一级二级三级毛片免费看| 亚洲性久久影院| 午夜亚洲福利在线播放| 国产白丝娇喘喷水9色精品| 色播亚洲综合网| 人人妻人人澡人人爽人人夜夜 | 欧美bdsm另类| 久久精品国产亚洲av天美| 日韩欧美一区二区三区在线观看| 99热这里只有是精品50| 97在线视频观看| 国产精品人妻久久久久久| 99九九线精品视频在线观看视频| 亚洲最大成人中文| 欧美成人一区二区免费高清观看| 高清午夜精品一区二区三区 | 国产精品无大码| 特大巨黑吊av在线直播| 国产高清激情床上av| 级片在线观看| АⅤ资源中文在线天堂| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 综合色丁香网| 三级经典国产精品| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 成人av在线播放网站| 亚洲精品日韩在线中文字幕 | 99久久中文字幕三级久久日本| 日本一本二区三区精品| 国产精品日韩av在线免费观看| 亚洲人成网站在线观看播放| 日本免费a在线| 九草在线视频观看| 高清毛片免费观看视频网站| 1000部很黄的大片| 国产精品一二三区在线看| 亚洲国产精品sss在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日韩制服骚丝袜av| 亚洲三级黄色毛片| 亚洲最大成人av| 久久久欧美国产精品| 免费大片18禁| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频| 我的老师免费观看完整版| 非洲黑人性xxxx精品又粗又长| 国产亚洲5aaaaa淫片| 如何舔出高潮| 久久中文看片网| 成人欧美大片| 一夜夜www| .国产精品久久| 午夜福利视频1000在线观看| 国内少妇人妻偷人精品xxx网站| 在线免费观看不下载黄p国产| 日产精品乱码卡一卡2卡三| 久久久国产成人精品二区| 男人舔女人下体高潮全视频| 国产又黄又爽又无遮挡在线| 色综合站精品国产| 亚洲国产精品合色在线| 日韩高清综合在线| 91久久精品电影网| 男人狂女人下面高潮的视频| 久久久久久久久大av| 老女人水多毛片| 亚洲欧洲日产国产| 不卡视频在线观看欧美| 最近手机中文字幕大全| 亚洲精品乱码久久久久久按摩| 国产单亲对白刺激| 国产精品一及| 中文字幕av在线有码专区| 亚洲乱码一区二区免费版| 男人狂女人下面高潮的视频| 久久精品国产清高在天天线| 麻豆一二三区av精品| 天美传媒精品一区二区| 亚洲av二区三区四区| 久久国内精品自在自线图片| 亚洲激情五月婷婷啪啪| 老司机福利观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲综合色惰| 天堂av国产一区二区熟女人妻| 欧美性感艳星| 精品不卡国产一区二区三区| 久99久视频精品免费| 国产真实伦视频高清在线观看| 亚洲第一电影网av| 国国产精品蜜臀av免费| www.av在线官网国产| 午夜免费激情av| av在线蜜桃| 美女高潮的动态| 精品一区二区免费观看| 精品久久久久久久人妻蜜臀av| 最近中文字幕高清免费大全6| 九九在线视频观看精品| 国产中年淑女户外野战色| 在线观看一区二区三区| 成人美女网站在线观看视频| 18禁在线播放成人免费| 国产老妇女一区| 久99久视频精品免费| 免费观看a级毛片全部| 最近最新中文字幕大全电影3| 九九热线精品视视频播放| 国产精品一区二区三区四区免费观看| 国产亚洲av嫩草精品影院| av在线观看视频网站免费| 别揉我奶头 嗯啊视频| 国产成人精品婷婷| 麻豆成人av视频| 亚洲婷婷狠狠爱综合网| 男的添女的下面高潮视频| 国产精品久久久久久久久免| 日韩欧美三级三区| 91狼人影院| 精品久久久久久久末码| 亚洲精品乱码久久久久久按摩| 午夜精品国产一区二区电影 | 国产大屁股一区二区在线视频| 欧美性猛交黑人性爽| 国产精品.久久久| 麻豆成人午夜福利视频| 婷婷色综合大香蕉| 国产日韩欧美在线精品| 一夜夜www| 国产午夜福利久久久久久| 国产69精品久久久久777片| 久久久成人免费电影| 成人欧美大片| 男女视频在线观看网站免费| 日韩一本色道免费dvd| 特大巨黑吊av在线直播| 哪里可以看免费的av片| 国产精品一区二区三区四区免费观看| 亚洲国产精品成人久久小说 | 国内精品久久久久精免费| 18禁在线播放成人免费| 你懂的网址亚洲精品在线观看 | 欧美日韩国产亚洲二区| 美女cb高潮喷水在线观看| av天堂在线播放| 久久精品国产鲁丝片午夜精品| 亚洲无线在线观看| 日本成人三级电影网站| 麻豆精品久久久久久蜜桃| 高清日韩中文字幕在线| 日日摸夜夜添夜夜添av毛片| 男插女下体视频免费在线播放| 国产极品天堂在线| 国产又黄又爽又无遮挡在线| 一级毛片久久久久久久久女| 永久网站在线| 极品教师在线视频| 亚洲国产欧洲综合997久久,| 一个人看视频在线观看www免费| 国产精品嫩草影院av在线观看| 成人三级黄色视频| 天堂网av新在线| 一区二区三区四区激情视频 | 亚洲精品456在线播放app| 亚洲18禁久久av| 人妻制服诱惑在线中文字幕| 欧美色视频一区免费| 日本色播在线视频| 国内精品久久久久精免费| av视频在线观看入口| 日本免费一区二区三区高清不卡| 直男gayav资源| 中国国产av一级| 成人特级黄色片久久久久久久| 成人二区视频| 国产黄a三级三级三级人| 黄色一级大片看看| 亚洲人与动物交配视频| 99热网站在线观看| 久久久久久久久久成人| 99在线视频只有这里精品首页| avwww免费| 18+在线观看网站| 美女内射精品一级片tv| 免费观看在线日韩| 亚洲欧美中文字幕日韩二区| 亚洲久久久久久中文字幕| 欧美日韩综合久久久久久| 18禁裸乳无遮挡免费网站照片| 天堂中文最新版在线下载 | 免费人成视频x8x8入口观看| 伊人久久精品亚洲午夜| 欧美激情在线99| 婷婷精品国产亚洲av| 韩国av在线不卡| 日韩欧美精品v在线| 一级毛片我不卡| av专区在线播放| 国产精品久久久久久精品电影| 成人美女网站在线观看视频| 九九热线精品视视频播放| 免费av观看视频| 岛国在线免费视频观看| av.在线天堂| 亚洲国产精品成人久久小说 | 亚洲婷婷狠狠爱综合网| 99视频精品全部免费 在线| 国产精品久久久久久精品电影小说 | 麻豆精品久久久久久蜜桃| 少妇被粗大猛烈的视频| 国产色爽女视频免费观看| 一级毛片电影观看 | 成人美女网站在线观看视频| 国产视频首页在线观看| 国产毛片a区久久久久| 国内精品宾馆在线| 亚洲精品久久久久久婷婷小说 | 狂野欧美白嫩少妇大欣赏| 91aial.com中文字幕在线观看| 国产91av在线免费观看| 国产一区二区激情短视频| 九九在线视频观看精品| 国产精品蜜桃在线观看 | 99热这里只有是精品50| 卡戴珊不雅视频在线播放| 久久99蜜桃精品久久| 美女被艹到高潮喷水动态| 午夜免费激情av|