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

    In-situ polymerized carbonate induced by Li-Ga alloy as novel artificial interphase on Li metal anode

    2023-10-14 03:03:04ZipingWngShuyunXieXuejieGoXinyngChenLinCongJunLiuHimingXieChungYuYulongLiu
    Chinese Chemical Letters 2023年9期

    Ziping Wng, Shuyun Xie, Xuejie Go, Xinyng Chen, Lin Cong, Jun Liu,Himing Xie,?, Chung Yu, Yulong Liu,?

    a School of Chemistry, Northeast Normal University, Changchun 130024, China

    b Center for Lignocellulosic Chemistry and Biomaterials, College of Light Industry and Chemical Engineering, Dalian Polytechnic University, Dalian 116034,China

    c State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

    Keywords:Lithium metal Alloy Polymer film High energy density Self-healing

    ABSTRACT Li metal is considered an ideal anode material because of its high theoretical capacity and low electrode potential.However, the practical usage of Li metal as an anode is severely limited because of inevitable parasitic side reactions with electrolyte and dendrites formation.At present, single-component artificial solid electrolyte interphase cannot simultaneously meet the multiple functions of promoting ion conduction, guiding lithium ion deposition, inhibiting dendrite growth, and reducing interface side reactions.Therefore, multi-component design on Li metal surface is widely investigated to achieve long-term cycling.Herein, we report a Li2Ga-carbonate polymer interphase layer to solve volume changes, Li dendrites formation and side-reactions.As a result, the Li symmetric cell can be stabilized at 3.0 mA/cm2 in carbonate electrolyte with limited volume of 20 μL.Coupled with 13.6 mg/cm2 (loading of 2 mAh/cm2) LiFePO4 cathode, discharge capacity retains at 90% for over 150 cycles under limited electrolyte conditions.With such an alloy-polymer interphase layer, higher energy density Li metal batteries become prominent in the near future.

    Anode materials with higher specific capacity are the most challenge to the next-generation battery systems [1-4].Among them, Li metal is considered as the ideal material for anodes because of the higher energy density of 3860 mAh/g and the lower potential of -3.04 V (vs.H2/H+) [5-7].However, organic electrolytes would be severely reduced by Li metal, resulting in formation of dendrites.Dendrites are lost from Li anode and become dead lithium during lithium stripping, causing low Columbic effi-ciency (CE).Practical application of lithium metal battery (LMB) is still an impossible at current stage [8,9].Therefore, two strategies have been proposed including construct conductive skeletons for Li deposition and an artificial interface protective layer.A hierarchical structure host/framework composed of carbon-based and metallic materials has a large surface area and exhibits good conductivity,which cannot only reduce the local current density, but also relieve lithium dendrites and accommodate volume change of Li [10-18].However, due to its increased inert weight and volume, the use of additional materials significantly reduces the energy density.Additionally, large specific areas consume more electrolytes and active Li, often causing safety hazards [19].

    In order to prevent inhomogeneous Li plating/stripping of metal anodes, a protection layer is crucial [20].Thereby, Li metal maintains at its specific capacity, but also decrease reactions between it and its electrolyte.The intrinsic solid electrolyte interphase (SEI) layer is usually mechanical instability during Li deposition/dissolution processes [21].The properties of SEI formation can be modulated by optimizing electrolyte compositions such as solvents, Li salts, and fluorinated compounds [22-25].Besides, artificial SEI layers can protect the Li metal [26], including Li+conductor [27,28], carbon nanofibers [29,30], metal oxides [31,32] and Li alloy [33-37].However, artificial SEI normally cannot heal itself,leading to gradual deterioration of its performance over time.Recently, liquid metals based SEI (self-healing [38], damage repairing[39]) shows great potential for Li metal anode modification.Among them, Ga-based liquid metals are more attractive due to their nontoxic, chemical stability, and strong bonding interactions [40-43].Self-healing ability of metal Ga in liquid electrolyte is limited, and the performance of self-healing electrodes remains unsatisfactory.Especially, other metallic elements are added to adjust surface tension and melting point.Wanget al.reveal that adopting soft substrates to alleviate the residual stress can suppress Li dendrites[44].In addition, commercial electrolyte without initiator forming a flexible membrane on Li electrode with simple process is important to achieve improved electrochemical properties.

    Fig.1.The mechanism of alloying and in-situ multicomponent SEI formation.(a) Proposed mechanism of SEI formation on Li anode.(b) Bare Li.(c) Ga dispersion on the surface of Li metal.(d) Li-Ga alloy on Li metal.(e) Li-Ga/polymer coated Li metal.(f) The proposed dimethyl carbonate polymerization reactions of EC via the ring-opening reaction.Lithium ethylene dicarbonate (LEDC), polycarbonate diols (PCDL).

    In this work, nanoparticle dispersion of gallium (Ga) was applied to Li foil by dropping it on the surface, and alloy layer (Li-Ga) was formed.Taking advantage of the Ga catalytic properties,carbonate polymer film was formed on the surface of Li-Ga alloy.The composite artificial interphase suppress the dendrites of Li and promotes uniform Li nucleation and deposition through the synergetic function of alloy and polymer.As a result, artificialinsituformed polycarbonate (AIFP) Li-Ga symmetric cells show excellent cycling stability for over 1300 h and hysteresis of only 3 mV at 1 mA/cm2.Besides, by matching AIFP Li-Ga with LiFePO4(LFP),LiNi0.8Co0.1Mn0.1O2(NCM811) and LiNiO.6Co0.2Mn0.2O2(NCM622)cathodes, the full cells achieve greater cycling stability, higher discharge capabilities and stable CE than bare Li anodes.A longterm at 0.5 C revealed that the AIFP Li-Ga||LFP cell delivers a high discharge capacity of 136.77 mAh/g after 150 cycles.This work presents a simple and effective method for the development of Li metal anodes with high performance.

    The proposed mechanism of alloy andin-situformed polycarbonate multicomponent SEI was shown in Fig.1a.Ga nanoparticles dispersed in tetrahydrofuran (THF) liquid (Fig.S1 in Supporting information), and was dropped onto the surface of Li foil (Fig.1b).Li alloy was in the shape of sphere of a diameter of 2-10 μm (Fig.1c).Alloy layers typically measure between 1 and 2 μm thick (Fig.S2 in Supporting information).Then, an alloy layer was formed, as seen from the darkness of Li metal surface (Fig.S3 in Supporting information).The alloy reaction was uniform, as evidenced by the energy dispersive spectroscopy (EDS) mapping of Ga element and backscattering electron image (Fig.1d).Dense and darker film was formed on the surface of Li-Ga alloy after the electrochemical polymerization process at 50 °C (Fig.1e).Under the same conditions,no film is found on the surface when bare lithium used (Fig.S4 in Supporting information).Polycarbonate film was formed through carbonate solvent ring-opening reaction, which is catalyzed by Li-Ga alloy (Fig.1f).

    To testify our assumption, X-ray diffraction (XRD) patterns of Li metal, Li-metal after alloy reaction and Li metal after polymerization reaction were recorded (Fig.2a).Li2Ga diffraction peaks replaced the pristine Li peaks after the first reaction, which proved the formation of Li-Ga alloy on the surface.A small amount of unreacted Li peak was observed in the sample after polymerization process, which was due to the consumption of partial Ga.The multicomponent SEI structure was further analyzed using transmission electron microscopy (TEM).From the TEM image, a Li-Ga nanoparticle of 100 nm length and 120 nm width was covered by a dense film as shown in Fig.2b.The nanoparticle was enlarged with higher magnification, and the lattice fringes in Fig.2c was assigned to Li2Ga with a spacing of 0.32 nm, corresponding to the Li2Ga (060) facet.To understand the affinity between electrolyte and anode, the contact angle tests were conducted (Fig.2d).For Li-Ga anode, the 0° at 0 s indicates that the electrolyte was fully wetted.The bare Li shows bad electrolyte wetting phenomenon.Even after 30 s, the contact angle was still as high as 38.5°.Li-Ga had better electrolyte wetting, which will facilitate the mass transport during reaction.

    To understand the polymer film composition, the13C and1H NMR spectra was presented in Figs.2e and f.Apart from the pristine peaks assigned to dimethyl carbonate (DMC, 3.72 ppm) and ethylene carbonate (EC, 4.50 ppm), several additional peaks were observed in the1H spectroscopy between 1 ppm and 4.5 ppm (Fig.2e) that were not present in the raw electrolyte (Fig.S5 in Supporting information).The results confirmed that the main products were Lithium ethylene di-carbonate (LEDC), 2-propoxyethyl hydrogen carbonate (PCDL), ethylene glycol and methanol.In the cycled electrolyte samples, these included an intense singlet at 4.36 (h proton), 3.69 (f proton), 3.66 (a proton), 3.52 (b proton), 3.45 (d proton) ppm, 3.18 (g proton) ppm that appears after the cycle and several multiple peaks at 4.12-4.05 (c proton), 1.38-1.37 (e proton) ppm.The h proton at 4.36 ppm was assigned to LEDC, the labeled f proton and g proton were assigned to ethylene glycol, the singlet a (ppm) and e (ppm) protons were assigned to methanol.The single d (ppm) proton, b (ppm) proton and c (ppm) proton analysis were the characteristic peak of PCDL.Moreover, as further support for the structure conjecture,13C NMR spectra (Fig.2f and Fig.S6 in Supporting information) was also obtained, and these results were consistent with NMR spectrum results obtained for1H NMR.The molecular weight content of the carbonate polymer was characterized by Mass Spectrometry, showing different molecular weight (Mw=344, 681 and 1019 g/mol), suggesting the carbonate polymer was an oligomer (Fig.S7 in Supporting information).Considering the mechanism of anionic ring-opening polymerization of caprolactone initiated by active alkali metals (such as K and Na)[45], this ring-opening polymerization is also expected to follow the anionic polymerization mechanism initiated by Li-Ga.Due to LEDC’s poor stability, the carbonate polymer with excellent flexibility properties was formed after blending with PCDL.The AIFP skeleton with a linear chain structure rich in ether bonds shows fast lithium ions (Li+) transportation, which is proved in polymers with similar structures, such as polyethylene oxide [46].Different polymerization pathways were proposed to synthesize EC and DMC(Fig.S8 in Supporting information), calculation result indicates that the polymerization process of electrolyte is spontaneously.

    To identify the detailed chemical composition of Li metal surface, XPS analysis with different etching time were recorded (Figs.2g-o).The XPS results proved that the Li2Ga-carbonate polymer interphase layer was divided into outer carbonate polymer at the top surface and inner Li2Ga alloy close to Li metal.Li2Ga alloy layer was composed of Li2Ga alloy (Figs.2g, j and m), and there was lithium fluoride at the bottom of the metal layer, which was the early product of the decomposition of lithium hexafluorophosphate in the electrolyte.Outer carbonate polymer was evidenced by carbonate polymer groups the C-C (284.7 eV), C-O(286.8 eV), and polycarbonate (290.0 eV) in C 1s high resolution XPS [47].After contacted with Li2Ga alloy, the polycarbonate content in the inner layer increases significantly according to C 1s XPS due to the decomposition of the electrolyte.(Figs.2i, l and o).Besides, no elemental gallium was detected in the surface layer, and Ga-Li (19.6 eV) was detected with increasing etching depth, indicating that the polycarbonate was covered with the surface of the metal layer from the Ga 3d spectrum (Figs.2h, k and n).Therefore,on top of the Li metal surface, we constructed a carbonate polymer layer based on Li2Ga alloy, which may maintains high stability over extended cycling.

    To understand the lithium plating behavior of protected Li metal, a symmetrical battery was assembled under 1 mA/cm2and 1 mAh/cm2As presented Fig.3, AIFP Li-Ga modified Li symmetrical cell showed smooth voltage plateaus and maintained a stable overpotential of 7 mV until 10 cycles (Fig.3g), implying that Li dissolution was uniform.At first, the abrupt voltage increased to 100 mV, corresponding to the big resistance of the initial stripping from the uncycled Li.At the end of the initial stripping, pits can be observed, implying that Li dissolution was not uniform.The bare lithium exhibits large voltage fluctuations, mainly arcing voltage profiles, with an overpotential of more than 60 mV in 10 cycles(Fig.3a).For negative electrode protected by Li-Ga alloy, the initial stripping process was optimized, and lithium was easier to be plated.But the voltage fluctuation increased gradually, and the arc voltage curve was still dominant.This was ascribing to the increase in the migration resistance of Li ions at an electrode surface.The polarization voltage reached 50 mV after 10 cycles (Fig.3d).

    The lithium surface morphology was checked after 10 cycles.Many holes with a diameter of 2 μm were observed on bare Li after the plating process (Fig.3b).For the Li-Ga protected anode, Ga may flow into the void where Li was extracted, so the surface became flat and only a few voids can be detected (Fig.3e).As for a hybrid layer coated Li, it exhibits a very smooth surface without any breakage on the surface (Fig.3h).The polycarbonate layer on the surface was retained after cycling.Moreover, under layer Li-Ga alloy was smooth with breakage as well(Figs.S9 and S10 in Supporting information).This proves that uniform Li plating was ascribed to the synergistic work of Li2Ga and polycarbonate.

    Fig.3.Galvanostatic cycling of the symmetric Li||Li, Li-Ga||Li-Ga and Li-Ga AIFP||AIFP Li-Ga symmetric cells at a controlled current density and fixed capacity from 1st to 10th cycles and the morphology after the cycle: (a) Li||Li, (b) SEM image of Li, (c) AFM image of Li, (d) Li-Ga||Li-Ga, (e) SEM image of Li-Ga, (f) AFM image of Li-Ga, (g) Li-Ga AIFP||AIFP Li-Ga, (h) SEM image of AIFP Li-Ga, (i) AFM image of AIFP Li-Ga.

    The atomic force microscopy (AFM) information could also reveal information about the surface morphology after 10 cycles(Figs.3c, f and i), AIFP Li-Ga exhibited smooth surface morphology with surface roughness approximately 35 nm (Fig.3i).However, due to a large number of holes, there was a rough morphology on the surface of bare Li, and over 350 nm of roughness was obtained (Fig.3c).Li-Ga protected Li surface was full of bumps,with over 200 nm of roughness (Fig.3f).Therefore, cycling consumes a great deal of active lithium.

    To assess the stability of Li during cycling, we performed a plating/stripping test for Li||Cu batteries.In Li||Cu batteries, the peak current was related to the lithium deposition capacity and electrochemical reaction activity.As shown in Figs.4a-c, due to the continuous cracking and repairing of the SEI layer, the peak current of Li||Cu cells with bare Li gradually decreases.It implied the gradual decline of Columbic efficiency, in line with Li symmetrical cell tested results.The peak current of Li||Cu cell with Ga-alloy gradually increased and remained stable owing to Li-Ga protection.Similarly, peak current of the AIFP Li-Ga cell protected by composite layer also increased gradually, and peaks in different cycles were well-overlapped, confirming the excellent plating/string of Li on Cu foil.In Li-Ga cell, the lithium insertion for Ga took place between 0.55 V and 0.73 V, while the anodic lithium extraction took place among 0.4, 0.75 and 0.92 V [48].Furthermore, oxidation peak appeared in Li-Ga coated cell at 1.5-2 V was assigned to electrolyte decomposition.However, AIFP could inhibit the decomposition reaction through the protection of polycarbonate film.

    The Li symmetrical cells (electrolyte dosage 20 μL) were assembled to evaluate the lithium cycling behavior under 1 mAh/cm2(Fig.4d).The overpotential for bare lithium reaches 100 mV after 225 h, and the voltage suddenly dropped after 380 h due to the formation of Li dendrite.Li-Ga-protected Li anode presented a small overpotential of 68 mV.However, a soft short circuit occurred after 400 h.AIFP Li-Ga cell delivers the lowest overpotential (about 3 mV) and a lifespan of more than 1300 h.When the current was 2 mA/cm2, AIFP Li-Ga symmetric cell also showed low overpotential fluctuation of 10.2 mV after 1640 h.To the contrast, Li-Ga and Li symmetric cells could only cycle for 213 h with high current density (Fig.S12 in Supporting information).

    To analyze the interfacial stability, EIS was carried out using symmetric Li cells (Fig.4e and Fig.S11 in Supporting information).The interfacial impedance of AIFP Li-Ga (26.6Ω) was comparatively smaller than the bare Li (51.1Ω) and the Li-Ga protected Li (57.4Ω) before the cycle due to double SEI layer.With the advance of the lithiation process, the impedance of Li-Ga cells decrease from 13.8Ωat 60 cycles to 9.1Ωat 100 cycles, respectively.The Li-Ga dropped to 29.9Ωfor the first 40 cycles and then increasing to 231.2Ωin the 100thcycle.In the case of bare Li anodes, impedance increased continuously as cycling progresses of 9.6Ωat 20th cycle to and 269.7Ωat 100thcycle.Especially, the resistance was 20 times higher than AIFP Li-Ga, which was in consistent with Fig.4a.Li-Ga protected anode showed a relative stable resistance in the range of 20thto 60thcycle, the resistance began to rise afterward.This indicated that Li-Ga alloy alone cannot inhibit the side-reactions between electrolyte and Li.

    Various current densities were used to measure the rate capability of AIFP Li-Ga, Li-Ga, and bare Li symmetric cells (Fig.4f).Even increasing the current density to 3 mA/cm2with constant area capacity (1 mAh/cm2), the overpotential of AIFP Li-Ga was stable at all current densities (9 mV at 1 mAh/cm2, 12 mV at 3 mAh/cm2).For bare Li, the overpotential increased from 30 mV at 1 mAh/cm2to 130 mV at 3 mAh/cm2.Li-Ga protected Li showed similar trend with different current density (from 25 mV at 1 mAh/cm2to 92 mV at 3 mAh/cm2).When the current was returned to 1 mA/cm2, the average overpotential for AIFP Li-Ga was retained(16 mV), in contrast to that of bare Li (381 mV) and Li-Ga (35 mV).

    Fig.4.(a-c) CV scanning of the carbonate electrolytes of different Li metal.(d) Galvanostatic cycling of the symmetric Li||Li, Li-Ga||Li-Ga and Li-Ga AIFP||AIFP Li-Ga symmetric cells at a controlled current density and fixed capacity: 1 mA/cm2, 1 mAh/cm2.(e) EIS evolution of Li symmetrical cell during cycling.(f) Rate performances of electrode symmetric cells at different current densities.Electrochemical performances of the battery with bare Li, Li-Ga and AIFP Li-Ga: (g) Li||LFP; (h) Cycling rate performance; (i)Li||NCM811.

    The stability of Li||LFP full cells (electrolyte usage=20 μL and cathode loading of 2 mAh/cm2) with bare Li, Li-Ga and AIFP Li-Ga was shown in Fig.4g.In comparison with the other two anodes,bare Li had the worst cycling stability.After 110 cycles, its capacity was only 26.88 mAh/g, with a capacity decay of 82.08%.The cycling stability of the LiFePO4full cells was improved on Li-Ga protected Li metal.Nevertheless, after 100 cycles, the capacity began to decline rapidly.Eventually, it reached 80.35 mAh/g after 150 cycles, and there was a 46.43% decay in capacity.However, the cell containing the AIFP Li-Ga was extremely stable, with the capacity of 136.77 mAh/g after 150 cycles, resulting in a capacity decay of only 8.82% (0.09 mAh/g decay per cycle).With cathode loading of 1 mAh/cm2and electrolyte usage=20 μL, AIFP Li-Ga still showed capability with the capacity of 126.58 mAh/g after 220 cycles.In comparison, batteries with Li and Li-Ga protected Li were dead after 160 cycles and 180 cycles as a result of the accumulation of dead lithium (Fig.S13 in Supporting information).

    The AIFP Li-Ga||LFP full cell offers superior rate performance(Fig.4h), displaying discharge capacities between 153.3, 144.4,120.1, 97.1 and 153.4 mAh/g at ampere densities of 0.5, 1, 2, 3 and 0.5 C, respectively.Discharge capacities of the Li-Ga||LFP full cells are 148.2, 136.5, 93.3, 0.8 and 149.4 mAh/g, and the Li||LFP full cell is only 148.4, 139.9, 112.4, 0.9 and 148.5 mAh/g, respectively.Especially, there is a significant difference in capacity at the three types of electrodes (AIFP Li-Ga, Li-Ga and bare Li) at high current densities.Moreover, this AIFP Li-Ga alloy structure exhibits stability even during rapid charge/discharge at high current density, indicating AIFP Li-Ga alloy structure is also beneficial to fast transfer of Li+.

    The full cells with LiNi0.6Co0.2Mn0.2O2(NCM622, loading of 0.8 mAh/cm2, electrolyte usage=20 μL) and LiNi0.8Co0.1Mn0.1O2(NCM811, loading of 0.6 mAh/cm2and electrolyte usage=20 μL)were prepared.For AIFP Li-Ga||NCM622cell, an initial capacity of 143.5 mAh/g was obtained at 0.5 C.At the end of 200 cycles, 129 mAh/g and 89.9% capacity retention were obtained, (0.07 mAh/g decay per cycle).Using Li-Ga and bare Li as anodes, the capacity decreased obviously after 120 cycles, and the capacity decreased to 43 mAh/g and 36 mAh/g after 200 cycles (Fig.S14 in Supporting information).Due to the side-reaction between Li and carbonate electrolyte, the liquid electrolyte dries out and dead lithium accumulated.With NCM811as cathode and AIFP Li-Ga as anode, the initial capacity was 152.6 mAh/g.The capacity was 180 mAh/g after 14 cycles and dropped to 123.7 mAh/g after 200 cycles (0.3 mAh/g dropped per cycle), corresponding to a 68.7% capacity retention(Fig.4i).Without exception, cell with bare Li decayed quickly after 80 cycles and dead after 100 cycles.With Li-Ga as the anode, the capacity decay was still higher than that of AIFP Li-Ga within 200 cycles, and dropped to 100 mAh/g after 200 cycles (0.43 mAh/g drop per cycle).

    The SEM images after cycling were analyzed, there were lithium dendrites structure on the bare Li metal (Fig.5a).The surface of Li-Ga protected Li became rough after cycling (Fig.5b).As expected,the surface of AIFP Li-Ga is uniform without any breakage.The film keeps its original morphology after repeated plating and stripping(Fig.5c).Even peeling off the polymer film by solvent, the underlying Li-Ga layer keeps smooth (Fig.S15 in Supporting informaation).Furthermore, the composition of Li-Ga protected Li and AIFP Li-Ga were checked by EDS (Figs.5d and e).For Li-Ga protected Li, the mass ratio of elements F to O was high, suggesting an inorganic compound was formed into large quantities (Fig.5d).In AIFP Li-Ga, large amount of C and O was detected due to the appearances of organic SEI layer (Fig.5e).Because of the formation of LiF,high content of F can also be detected.Obviously, composite artificial layer can prevent the Li preferential deposition and result in smooth morphology.

    In summary, we report a Li2Ga-carbonate polymer interphase layer to solve the volume changes, side reactions and Li dendrite formation of Li metal anode.Carbonate polymer is closely covered on the Li-Ga alloy layer, Li dendrites are effectively suppressed due to the soft carbonate polymer.Li-Ga can fill the small gaps during plating process timely, thus preventing microcracks formation.As a result, Li symmetrical cells with AIFP Li-Ga keep stable at 1 mAh/cm2for 1300 h (electrolyte usage=20 μL).When AIFP Li-Ga was coupled with LFP cathode (cathode loading 2 mAh/cm2),reversible capacity delivers 137 mAh/g after 150 cycles.Using NCM622and NCM811, discharge capacity and capacity of AIFP Li-Ga shows better performance than Li-Ga and bare Li over 200 cycles at 0.5 C rate.

    Fig.5.Bird-view SEM images of Li-stripping electrode.(a) Li.(b) Li-Ga.(c) AIFP Li-Ga.Characterizations of (d) Li-Ga and (e) AIFP Li-Ga with TESM.

    Declaration of competing interest

    The authors declare that they have no conflict of interest.

    Acknowledgments

    This work was supported by Jilin Province Science and Technology Department Major Science and Technology Project (Nos.20220301004GX, 20220301005GX), Key Subject Construction of Physical Chemistry of Northeast Normal University, National Natural Science Foundation of China (Nos.21905110, 22102020); National Natural Science Foundation of China (No.21905041), Special foundation of Jilin Province Industrial technology Research and Development (No.2019C042), the Fundamental Research Funds for the Central Universities (No.2412020FZ008).

    Supplementary materials

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

    中文欧美无线码| 日本wwww免费看| 亚洲国产精品一区三区| 亚洲成人手机| 日韩中文字幕欧美一区二区 | 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| 女警被强在线播放| 日韩大片免费观看网站| 欧美另类一区| 亚洲精品成人av观看孕妇| 日韩视频在线欧美| 黄频高清免费视频| 国产黄色免费在线视频| 午夜激情久久久久久久| 亚洲欧美一区二区三区久久| 久久国产精品影院| 国产成人欧美| 99久久99久久久精品蜜桃| 日日摸夜夜添夜夜爱| 亚洲av片天天在线观看| 在线观看免费高清a一片| 高清不卡的av网站| 99国产精品99久久久久| 丰满少妇做爰视频| 国语对白做爰xxxⅹ性视频网站| 欧美日韩国产mv在线观看视频| 欧美 亚洲 国产 日韩一| 久久久久网色| 手机成人av网站| 老汉色∧v一级毛片| 亚洲 国产 在线| 日韩制服丝袜自拍偷拍| 男女午夜视频在线观看| 亚洲国产精品国产精品| 97精品久久久久久久久久精品| 国产成人欧美| 涩涩av久久男人的天堂| 男女无遮挡免费网站观看| 国产精品国产av在线观看| www.999成人在线观看| 精品人妻熟女毛片av久久网站| 亚洲av日韩精品久久久久久密 | 日本wwww免费看| 啦啦啦视频在线资源免费观看| 飞空精品影院首页| 夫妻性生交免费视频一级片| 日韩中文字幕视频在线看片| 免费在线观看视频国产中文字幕亚洲 | 精品国产国语对白av| 亚洲专区中文字幕在线| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 男人爽女人下面视频在线观看| 亚洲第一av免费看| 日韩人妻精品一区2区三区| 日韩中文字幕视频在线看片| 少妇 在线观看| 赤兔流量卡办理| 色精品久久人妻99蜜桃| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩av久久| 男人添女人高潮全过程视频| 亚洲精品自拍成人| 国产片内射在线| 久久久久国产精品人妻一区二区| 日韩欧美一区视频在线观看| 免费人妻精品一区二区三区视频| 高清不卡的av网站| 国产在线观看jvid| 久久精品国产亚洲av高清一级| 首页视频小说图片口味搜索 | 汤姆久久久久久久影院中文字幕| 久久久久久久精品精品| 人人妻,人人澡人人爽秒播 | 久久精品国产综合久久久| 韩国精品一区二区三区| 欧美xxⅹ黑人| 亚洲精品国产色婷婷电影| 亚洲五月色婷婷综合| 国产一卡二卡三卡精品| 十八禁人妻一区二区| 国产高清国产精品国产三级| 91字幕亚洲| 久久国产精品大桥未久av| 日韩大片免费观看网站| 丝袜脚勾引网站| 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 国产人伦9x9x在线观看| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 国产亚洲欧美精品永久| 久久精品久久久久久久性| 国产一区二区三区综合在线观看| 性少妇av在线| 下体分泌物呈黄色| 精品少妇黑人巨大在线播放| 十分钟在线观看高清视频www| 一本综合久久免费| 国产精品久久久人人做人人爽| 欧美亚洲 丝袜 人妻 在线| 9色porny在线观看| 欧美日韩国产mv在线观看视频| 新久久久久国产一级毛片| 国产精品av久久久久免费| 少妇的丰满在线观看| 在现免费观看毛片| 亚洲,欧美,日韩| 精品熟女少妇八av免费久了| 老司机深夜福利视频在线观看 | 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 一区福利在线观看| 中文字幕人妻熟女乱码| 久久精品亚洲熟妇少妇任你| 成年动漫av网址| 中国美女看黄片| 久久综合国产亚洲精品| 男女午夜视频在线观看| 久久99精品国语久久久| 国产成人a∨麻豆精品| 新久久久久国产一级毛片| 满18在线观看网站| 啦啦啦中文免费视频观看日本| 久久久久久久国产电影| 欧美精品一区二区大全| 高清不卡的av网站| 国产亚洲av片在线观看秒播厂| 大型av网站在线播放| 中文字幕最新亚洲高清| 男女午夜视频在线观看| 一区二区三区精品91| 日日爽夜夜爽网站| 一区二区日韩欧美中文字幕| 亚洲熟女精品中文字幕| 精品国产国语对白av| 青春草视频在线免费观看| 又大又爽又粗| 99九九在线精品视频| 在线观看一区二区三区激情| 免费观看a级毛片全部| 亚洲中文日韩欧美视频| 波多野结衣av一区二区av| 亚洲av日韩精品久久久久久密 | 伦理电影免费视频| 三上悠亚av全集在线观看| 叶爱在线成人免费视频播放| 久久国产精品人妻蜜桃| 黑人巨大精品欧美一区二区蜜桃| av网站免费在线观看视频| 日本猛色少妇xxxxx猛交久久| 99久久综合免费| 国产高清视频在线播放一区 | 天天躁夜夜躁狠狠躁躁| 亚洲第一青青草原| 日韩制服骚丝袜av| 一级黄片播放器| 在线天堂中文资源库| 一级毛片 在线播放| 丝袜脚勾引网站| 亚洲久久久国产精品| 日韩一卡2卡3卡4卡2021年| 中文乱码字字幕精品一区二区三区| 免费看不卡的av| 日本一区二区免费在线视频| 亚洲欧美一区二区三区国产| 精品久久久久久久毛片微露脸 | 精品一区二区三区av网在线观看 | 一本久久精品| 午夜免费观看性视频| 波多野结衣一区麻豆| 天天躁日日躁夜夜躁夜夜| 在线观看免费视频网站a站| 后天国语完整版免费观看| 91九色精品人成在线观看| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 国产免费一区二区三区四区乱码| 亚洲视频免费观看视频| 少妇人妻久久综合中文| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 国产精品 欧美亚洲| 美女视频免费永久观看网站| 国产亚洲精品久久久久5区| www.熟女人妻精品国产| 国产野战对白在线观看| 国产一区二区三区av在线| 午夜91福利影院| 一区二区日韩欧美中文字幕| 亚洲免费av在线视频| 亚洲欧美精品综合一区二区三区| 国产成人精品在线电影| 丰满迷人的少妇在线观看| 精品亚洲乱码少妇综合久久| 电影成人av| 久久国产亚洲av麻豆专区| 男人添女人高潮全过程视频| 国产精品秋霞免费鲁丝片| 欧美日韩精品网址| 人体艺术视频欧美日本| 欧美精品av麻豆av| 午夜福利视频在线观看免费| 99精品久久久久人妻精品| 久久久久网色| 人人妻人人澡人人爽人人夜夜| 欧美少妇被猛烈插入视频| 成人国产一区最新在线观看 | 国产成人欧美| 日韩av不卡免费在线播放| 在线观看免费高清a一片| 日韩精品免费视频一区二区三区| 日韩伦理黄色片| 欧美精品亚洲一区二区| 成人亚洲精品一区在线观看| 久久热在线av| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线| 99热全是精品| 交换朋友夫妻互换小说| 最近手机中文字幕大全| 视频区欧美日本亚洲| 男女高潮啪啪啪动态图| 人人妻,人人澡人人爽秒播 | 一边摸一边做爽爽视频免费| 51午夜福利影视在线观看| 精品一区二区三区av网在线观看 | a级毛片黄视频| 精品人妻熟女毛片av久久网站| 午夜激情久久久久久久| 丁香六月天网| 一区二区av电影网| 悠悠久久av| 欧美精品一区二区大全| 国产精品久久久久成人av| 老汉色av国产亚洲站长工具| 男女高潮啪啪啪动态图| 美女中出高潮动态图| 女警被强在线播放| 黄色 视频免费看| 男女下面插进去视频免费观看| 日日摸夜夜添夜夜爱| 午夜福利一区二区在线看| 别揉我奶头~嗯~啊~动态视频 | 久久精品aⅴ一区二区三区四区| 亚洲欧洲国产日韩| 建设人人有责人人尽责人人享有的| 波野结衣二区三区在线| 亚洲,一卡二卡三卡| 久久久精品免费免费高清| 99久久精品国产亚洲精品| 一区二区三区激情视频| 老鸭窝网址在线观看| 大片免费播放器 马上看| 国产又色又爽无遮挡免| 国产男女内射视频| 国产成人系列免费观看| 极品人妻少妇av视频| 妹子高潮喷水视频| 久久精品成人免费网站| 亚洲精品日本国产第一区| 亚洲av日韩精品久久久久久密 | h视频一区二区三区| 国产麻豆69| 热re99久久精品国产66热6| 丝袜美腿诱惑在线| 热re99久久国产66热| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站| 三上悠亚av全集在线观看| 久久精品成人免费网站| 日韩 亚洲 欧美在线| 你懂的网址亚洲精品在线观看| 啦啦啦在线免费观看视频4| 免费看不卡的av| 婷婷丁香在线五月| 国产精品麻豆人妻色哟哟久久| 香蕉国产在线看| 亚洲欧美一区二区三区国产| 精品福利永久在线观看| 五月天丁香电影| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 久久久精品免费免费高清| 国产成人欧美在线观看 | 大码成人一级视频| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 老司机靠b影院| 亚洲精品自拍成人| 一二三四在线观看免费中文在| 久久人妻熟女aⅴ| 中文字幕高清在线视频| 大香蕉久久网| 久久精品久久精品一区二区三区| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 美女午夜性视频免费| 免费不卡黄色视频| 男人舔女人的私密视频| 久久精品久久久久久噜噜老黄| 啦啦啦啦在线视频资源| 十八禁人妻一区二区| 老司机深夜福利视频在线观看 | 亚洲,欧美,日韩| 国产成人啪精品午夜网站| 91成人精品电影| 亚洲欧美成人综合另类久久久| 热re99久久精品国产66热6| 色精品久久人妻99蜜桃| 一二三四在线观看免费中文在| 成人亚洲欧美一区二区av| 国产精品三级大全| 亚洲专区国产一区二区| 国产欧美日韩一区二区三 | 免费女性裸体啪啪无遮挡网站| 国产一区二区 视频在线| 国产成人系列免费观看| 亚洲黑人精品在线| 黄色a级毛片大全视频| 国产欧美日韩一区二区三 | 免费在线观看日本一区| 欧美人与善性xxx| www.熟女人妻精品国产| 免费看不卡的av| √禁漫天堂资源中文www| 亚洲欧美色中文字幕在线| 亚洲伊人久久精品综合| 女人高潮潮喷娇喘18禁视频| 99热网站在线观看| 丝袜喷水一区| 18禁裸乳无遮挡动漫免费视频| 黑人巨大精品欧美一区二区蜜桃| 精品国产一区二区久久| 国产在线一区二区三区精| 制服诱惑二区| 老司机靠b影院| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 十八禁人妻一区二区| 亚洲专区国产一区二区| 国产精品国产三级专区第一集| 老司机亚洲免费影院| 久久中文字幕一级| 我要看黄色一级片免费的| 黑丝袜美女国产一区| 亚洲av片天天在线观看| 99精国产麻豆久久婷婷| 欧美日韩视频高清一区二区三区二| 亚洲精品一区蜜桃| 久久久久久亚洲精品国产蜜桃av| 飞空精品影院首页| 色精品久久人妻99蜜桃| 下体分泌物呈黄色| 亚洲一区二区三区欧美精品| 国产亚洲av高清不卡| 日韩人妻精品一区2区三区| 老司机在亚洲福利影院| 久久免费观看电影| 国产精品国产三级国产专区5o| 中国美女看黄片| 一级黄色大片毛片| 国产色视频综合| 妹子高潮喷水视频| 欧美日韩福利视频一区二区| 欧美+亚洲+日韩+国产| bbb黄色大片| 日韩制服丝袜自拍偷拍| 女人精品久久久久毛片| 免费少妇av软件| 少妇的丰满在线观看| 国产成人系列免费观看| 国产激情久久老熟女| 久久精品aⅴ一区二区三区四区| www.精华液| 在线观看免费视频网站a站| 欧美中文综合在线视频| 国产午夜精品一二区理论片| 黑人欧美特级aaaaaa片| 欧美日本中文国产一区发布| 男男h啪啪无遮挡| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 黑丝袜美女国产一区| 欧美乱码精品一区二区三区| 男女免费视频国产| 午夜激情av网站| 亚洲精品自拍成人| 亚洲精品中文字幕在线视频| 精品熟女少妇八av免费久了| 十八禁高潮呻吟视频| 一级毛片 在线播放| 欧美精品高潮呻吟av久久| 91麻豆av在线| av视频免费观看在线观看| 黑人猛操日本美女一级片| 国产精品国产av在线观看| 制服人妻中文乱码| 看免费成人av毛片| 亚洲国产最新在线播放| 高清黄色对白视频在线免费看| 欧美日韩av久久| 中国美女看黄片| 一区二区三区激情视频| 欧美日韩福利视频一区二区| 涩涩av久久男人的天堂| 一级片'在线观看视频| 国产主播在线观看一区二区 | av有码第一页| 亚洲欧美一区二区三区久久| 男女下面插进去视频免费观看| 日韩 亚洲 欧美在线| 最新的欧美精品一区二区| 亚洲成人免费电影在线观看 | av在线老鸭窝| 黑人猛操日本美女一级片| 熟女av电影| 一级毛片女人18水好多 | 美女国产高潮福利片在线看| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 精品国产乱码久久久久久小说| 亚洲伊人色综图| 狂野欧美激情性xxxx| 日韩视频在线欧美| 亚洲av电影在线观看一区二区三区| 在线精品无人区一区二区三| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看 | avwww免费| 欧美精品av麻豆av| 九草在线视频观看| 夜夜骑夜夜射夜夜干| 亚洲黑人精品在线| 妹子高潮喷水视频| 少妇人妻久久综合中文| 日韩精品免费视频一区二区三区| 国产精品一区二区免费欧美 | 国产精品秋霞免费鲁丝片| 一个人免费看片子| 国产片内射在线| 91精品三级在线观看| 午夜日韩欧美国产| cao死你这个sao货| 精品国产超薄肉色丝袜足j| 99热国产这里只有精品6| videos熟女内射| 久久久国产一区二区| 宅男免费午夜| 亚洲中文日韩欧美视频| 国产成人精品久久二区二区91| 久久人人爽人人片av| 69精品国产乱码久久久| 精品人妻一区二区三区麻豆| 欧美变态另类bdsm刘玥| 精品一区二区三区av网在线观看 | 黑人巨大精品欧美一区二区蜜桃| 久久久久国产一级毛片高清牌| 久久毛片免费看一区二区三区| 校园人妻丝袜中文字幕| 久久中文字幕一级| 又紧又爽又黄一区二区| 国产欧美日韩综合在线一区二区| 日本91视频免费播放| 男女边吃奶边做爰视频| 成年人免费黄色播放视频| 一区二区三区乱码不卡18| av网站免费在线观看视频| 热re99久久精品国产66热6| 国产精品秋霞免费鲁丝片| 久久国产精品男人的天堂亚洲| 中文字幕最新亚洲高清| 久久精品人人爽人人爽视色| 亚洲国产精品成人久久小说| 后天国语完整版免费观看| 日本五十路高清| 成人三级做爰电影| 十八禁网站网址无遮挡| 婷婷色综合www| 日韩av不卡免费在线播放| 亚洲国产欧美日韩在线播放| 亚洲精品乱久久久久久| a级片在线免费高清观看视频| 午夜福利视频精品| 男的添女的下面高潮视频| 蜜桃国产av成人99| 18禁裸乳无遮挡动漫免费视频| 午夜免费男女啪啪视频观看| 中文字幕高清在线视频| 伊人亚洲综合成人网| 妹子高潮喷水视频| 成人国语在线视频| 亚洲国产精品999| 99久久99久久久精品蜜桃| av网站免费在线观看视频| 美女视频免费永久观看网站| 狂野欧美激情性xxxx| 亚洲第一青青草原| 天天躁日日躁夜夜躁夜夜| 男女边摸边吃奶| 午夜精品国产一区二区电影| xxx大片免费视频| 一本一本久久a久久精品综合妖精| 亚洲色图综合在线观看| 亚洲一码二码三码区别大吗| 国产一级毛片在线| 成年女人毛片免费观看观看9 | 一边亲一边摸免费视频| 亚洲国产日韩一区二区| 欧美国产精品va在线观看不卡| 久久99精品国语久久久| 精品亚洲成a人片在线观看| 久久久精品94久久精品| 视频区欧美日本亚洲| a 毛片基地| 日韩伦理黄色片| 精品国产一区二区三区四区第35| 人人澡人人妻人| 欧美激情高清一区二区三区| 久久精品久久久久久久性| 一级毛片 在线播放| 国产91精品成人一区二区三区 | 少妇人妻 视频| 悠悠久久av| 国产精品成人在线| 99热全是精品| 两性夫妻黄色片| 好男人电影高清在线观看| 国产精品亚洲av一区麻豆| 亚洲精品久久成人aⅴ小说| 国产人伦9x9x在线观看| h视频一区二区三区| www.熟女人妻精品国产| 亚洲av美国av| 国产男女内射视频| 国产在线观看jvid| 爱豆传媒免费全集在线观看| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 欧美激情极品国产一区二区三区| 热re99久久国产66热| 久久人妻福利社区极品人妻图片 | 久久久国产精品麻豆| 99九九在线精品视频| 国产熟女欧美一区二区| 亚洲五月色婷婷综合| 性少妇av在线| 亚洲av美国av| 99国产精品一区二区三区| 国产熟女午夜一区二区三区| 久久国产精品男人的天堂亚洲| 免费在线观看日本一区| 男的添女的下面高潮视频| 日韩制服丝袜自拍偷拍| 欧美精品av麻豆av| 精品国产乱码久久久久久男人| 另类精品久久| 一级片'在线观看视频| 国产三级黄色录像| 亚洲精品第二区| 久久久久久亚洲精品国产蜜桃av| 欧美黑人欧美精品刺激| 久久久久网色| 亚洲av成人精品一二三区| kizo精华| 一本久久精品| 欧美人与性动交α欧美精品济南到| 男女无遮挡免费网站观看| 99国产精品99久久久久| 日韩大码丰满熟妇| 亚洲欧美色中文字幕在线| 一区二区三区四区激情视频| 天天躁夜夜躁狠狠久久av| av天堂在线播放| 久久久久久久久免费视频了| 大型av网站在线播放| 嫩草影视91久久| 亚洲av电影在线观看一区二区三区| 性色av乱码一区二区三区2| 久久久久久久大尺度免费视频| 多毛熟女@视频| 人体艺术视频欧美日本| 美女视频免费永久观看网站| 国产成人啪精品午夜网站| 国产高清不卡午夜福利| 欧美日韩黄片免| 又粗又硬又长又爽又黄的视频| videosex国产| 国产亚洲一区二区精品| 久久久国产一区二区| 亚洲成色77777| 久久久精品区二区三区| 大片电影免费在线观看免费| 国产精品.久久久| 黄网站色视频无遮挡免费观看| xxx大片免费视频| 亚洲av日韩精品久久久久久密 | 大型av网站在线播放| 亚洲图色成人| 欧美日韩视频高清一区二区三区二| 欧美精品一区二区大全| 精品人妻一区二区三区麻豆| 一级毛片女人18水好多 | 91字幕亚洲| 久久国产精品人妻蜜桃| 国产福利在线免费观看视频| 亚洲精品一区蜜桃| 校园人妻丝袜中文字幕| 国产精品 国内视频| 午夜福利视频在线观看免费| 中文字幕av电影在线播放| 亚洲人成电影观看| 日韩 欧美 亚洲 中文字幕| 99九九在线精品视频|