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

    Enhancement of electrochemical performance in lithium-ion battery via tantalum oxide coated nickel-rich cathode materials

    2022-05-16 07:12:08FenglingChen陳峰嶺JiannanLin林建楠YifanChen陳一帆BinbinDong董彬彬ChujunYin尹楚君SiyingTian田飔瑩DapengSun孫大鵬JingXie解婧ZhenyuZhang張振宇HongLi李泓andChaoboLi李超波
    Chinese Physics B 2022年5期
    關(guān)鍵詞:大鵬

    Fengling Chen(陳峰嶺) Jiannan Lin(林建楠) Yifan Chen(陳一帆) Binbin Dong(董彬彬)Chujun Yin(尹楚君) Siying Tian(田飔瑩) Dapeng Sun(孫大鵬) Jing Xie(解婧)Zhenyu Zhang(張振宇) Hong Li(李泓) and Chaobo Li(李超波)

    1Institute of Microelectronics,Chinese Academy of Sciences,Beiing 100029,China

    2University of Chinese Academy of Sciences,Beijing 100049,China

    3Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

    4Beijing Welion New Energy Technology Co.,Ltd,Beijing 102402,China

    Keywords: LiNi0.88Mn0.03Co0.09O2,tantalum oxide,surface coating,lithium-ion battery,cathode material

    1. Introduction

    As energy storage devices, lithium-ion batteries are thought to be able to change the traditional energy supply structure by harnessing renewable energy sources. Nowadays, commercialized lithium-ion batteries have been widely applied in the markets of 3C (consumer, computer, and communication), electric tools, and energy storage, which is benefited from their excellent performance, such as high energy density and long cycle life.[1]Among all the components of the lithium-ion battery, cathode material has the highest cost ratio,and it is also one of the important parts that will greatly limit battery performance.[2—8]Cathode material can be divided into three types from the perspective of structure,that is layered, spinel and polyanion originated from Goodenough’s group.[2]And attracted by higher reversible capacity,the layered nickel-rich cathode material (NRM) has been chosen to achieve a longer range for electric vehicles.

    The mentioned NRM is generally referred to a polynary derivative of LiNiO2, and in which, Ni ions are partially replaced by other ions (Co, Mn, Al,etc.). And the high reversible capacity of NRM is stemmed from the slight overlaps of Ni3+/4+egand O 2p band.[9]During the charging process, the NRM undergoes a series of phase transitions from the original hexagonal phase (H1), monoclinic (M), another hexagonal phase (H2), and finally to a third hexagonal phase(H3).[10]From the structural perspective, the great change of anisotropic lattice, a and c axis during the H2 and H3 phase transition will cause large volume change,resulting in cathode particle microcracks.[11]The generation of microcracks facilitates electrolyte penetration and cathode particle pulverization,which will further aggravate the degradation of long-term cycle performance.[12]Additionally,the similar ionic radius of Ni2+(0.69 °A)and Li+(0.76 °A),as well as the magnetic frustration of Jahn—Teller active Ni3+(t62ge1g) has been believed to be the main culprit of cation mixing.[11,13,14]These factors contribute to voltage decay and capacity fading at the cathode level, especially under high temperatures.[15]Nevertheless, the surface deterioration of NRM is normally attributed to the presence of a large amount of Ni ions with high valance(Ni3+and Ni4+), which has high oxidation activity and will further interact with moisture in the air during storage and electrolyte in delithiation state.[16]It is worth noting that all those factors interplay with each other and contribute to the failure of NRM.

    A number of improved methods acted on the cathode material bulk and surface has been developed to mitigate these degradation issues of NRM,such as bulk/surface doping,surface coating, single crystal, heterostructure, and concentration gradient design.[15,17,18]Among the aforementioned strategies,surface modification is one of the most straightforward approaches. The surface state of the cathode material is closely related to many undesired reactions, which further lead to surface reconstruction,formation of the resistive layer,and transition metal desolution.[19]And a reasonable cathode—electrolyte interface has a decisive impact on the overall efficiency of the batteries. Therefore,the surface design of cathode material is very necessary. Recently, tantalum has been found to be a very useful element for bulk/surface doping and surface coating of cathode material.[20—22]And according to the research results, tantalum possesses the strongest binding energy with oxygen (839 kJ/mol), which can inhibit oxygen release during lithium extraction.[23]The similar radius of Ta5+(0.64 °A) with Ni2+(0.69 °A) can promote the occupation of lithium site instead,[24]and Ta5+with higher valance in lithium sites can also compensate for the charge variation during lithium migration via screening oxygen ion repulsion.[24,25]Whereas there are few studies about surface coating by tantalum element,and further investigation is needed. In this work, the electrochemical performance of LiNi0.88Mn0.03Co0.09O2(NMC)was well improved via functional surface modification of tantalum oxide,and the reasons were also studied.

    2. Experimental details

    2.1. Materials

    The NMC powder was made by Beijing Welion New Energy Technology Company Limited. The detailed preparation process was described in supporting information.A wet chemical method was chosen to perform the coating process.An optimized amount of tantalum ethoxide (Ta(C2H5O)5, 99.99%,Aladdin) was firstly dissolved into ethanol and then NMC powder was added under stirring, and after the ethanol evaporated in oil bath followed by further heating at 120°C for 30 min, the coated LiNi0.88Mn0.03Co0.09O2(C-NMC) power was obtained. The C-NMC powder was further annealed at 550°C for 12 hours in an oxygen atmosphere to get the final sample, that is CA-NMC powder. The mass fraction of tantalum elements in CA-NMC is 0.4%. Additionally, the bare NMC (B-NMC) powder was prepared by the same process without adding tantalum ethoxide. All the above powder was stored in a vacuum. The coating amount is optimized in advance(0.1%,0.2%,0.3%,0.4%,and 0.5%)and choose 0.4%as research target due to its superior overall performance.

    2.2. Coin cell assembling

    To prepare the positive electrode,firstly,the cathode materials(B-NMC,C-NMC,and CA-NMC),carbon black(super P,Timcal)and polyvinylidenefluoride(PVDF,Sovlay Group)with a weight ratio of 90:5:5 were dispersed in N-methyl pyrrolidone(NMP,MYJ Chemical Co.,Ltd.)solvent,and then the resulting slurry was coated on aluminum foil followed by drying at 120°C for~12 hours in a vacuum oven. The test cells were assembled into CR2032 coin-type cells in an argonfilled glove box,with Li foil as the anode and Celgard 3501 as the separator.The electrode diameter was 12 mm and the cathode material load was around 12.0 mg/cm2. The amount of electrolyte(1-M LiPF6 dissolved in EC/DEC/EMC with volume ratio=1:1:1) used in each coin cell was 90 μL. All the cells were placed for 6 hours before the test to ensure sufficient infiltration of the electrolyte.

    2.3. Electrochemical test

    All the cells were activated at 0.1 C (1 C=220 mAh/g)for two cycles before the subsequent cycle(1 C/1 C)and rate(0.1 C/0.1 C,0.2 C/0.2 C,0.5 C/0.5 C,1 C/1 C,5 C/5 C,and 0.1 C/0.1 C)tests,and those electrochemical tests were carried out using galvanostatic charge/discharge method (LANHE,CT2001A) between 3.0 V—4.4 V (versusLi/Li+) at 25°C or 45°C. Cyclic voltammetry (CV) tests were performed on an RST 5200F electrochemical workstation at a scan rate of 0.1 mV and a voltage range of 3.0 V—4.4 V. Electrochemical impedance spectroscopy (EIS) was conducted on an RST 5200F electrochemical workstation with a voltage amplitude of 5 mV and frequency range of 0.1 Hz—106Hz.

    2.4. Characterization

    The cathode materials were measured by powder xray diffraction (XRD, Brucker D8 Advanced) to identify the sample phase. The sample morphology and microstructure were acquired on scanning electron microscopy(SEM,Hitach SU8200) and transmission electron spectroscopy (TEM, FEI Talos F200X). Besides, energy dispersive x-ray spectroscopy(EDS) mapping was implemented by using EDS equipment connected to SEM. X-ray photoelectron spectroscopy (XPS,Thermo escalab 250Xi)analysis was conducted to investigate the valance state of the main element composition on the surface of materials. The cross-section of cycled electrodes was prepared by an ion milling machine(Gatan 691). DSC(differential scanning calorimetry, Jupiter Netzsch) was performed from 25°C to 500°C at a heating rate of 5°C/min under a nitrogen atmosphere.

    3. Results and discussion

    3.1. Morphology and structure

    The SEM images of NMC materials before and after coating are shown in Figs.1(a)—1(c). The spherical cathode aggregates are composed of primary particles below 400 nm, and the aggregates remain intact after the wet chemical treatment and calcination process. And according to the inset images of higher magnification,there is no apparent difference within the surface of NMC samples before and after coating, indicating that the coating or annealing process induces a minimal change in the morphology. The EDS mapping images(Fig. 1(d)) present the existence and uniform distribution of Ta element in the CA-NMC sample, which means that the modification of tantalum oxide is believable and successful.High-resolution transmission electron microscopy (HRTEM)and fast Fourier transform(FFT)measurements were utilized to estimate the coating thickness and the surface structure change after coating or heat treatment,and the results are displayed in Figs.1(e)—1(g). The coating layer with a thickness of~8 nm on the surface of the C-NMC sample can be observed in Fig.1(e),and it is decreased to~6 nm after calcination treatment(Fig.1(f)). However,the surface of the B-NMC sample shows a smooth edge without any coatings(Fig.1(g)).Meanwhile, it can be seen that all these three samples possess a layered structure near the surface. The crystal faces of(006),(113),(104),(107),(101)are well-identified in the FFT patterns of the three samples, which implies that the coating or calcination process does not change the hexagonal layered crystalline structure of the NMC materials.[22,26]All the above results demonstrate that the surface structure is retained and consistent with the bulk after surface coating of tantalum oxide.

    As shown in Fig.2(a),the XRD patterns of C-NMC,CANMC,and B-NMC indicate that all the cathode samples have a hexagonalα-NaFeO2structure with good crystallinity (R-3mspace group, PDF#97-016-9367). Arising from the low amount of coating materials,the diffraction peaks of tantalum oxide is absent in C-NMC and CA-NMC sample (Figs. 2(a)and 2(b)). The typical (003) peak of C-NMC and CA-NMC slightly shift to a lower degree in comparison with B-NMC(Fig. 2(b)). The slight shift cannot exactly indicate the increase ofcaxis and surface doping of coating materials. And the XRD fitting results are summarized in Table 1,in which the Rwp represent overall error in fitting. Moreover, all of thea,caxis and the volume increase after coating,which may stem from the larger ionic radius of Ta5+(0.64 °A)compared to Ni3+(0.56 °A).[27]All three samples exhibit good layered structure and a low degree of cation mixing,confirmed by the value of peak intensity ratio (I003/I104>1.2) andc/a(>4.9).[28,29]The distinct splitting of (006)/(012) and (108)/(110) peaks point out that the materials have a good hexagonal structure(Figs.2(c)and 2(d)).[30]It can be seen that the introduction of tantalum oxide has no bad effect on the layered nature of the bulk phase.

    To study the influence of coated tantalum oxide on the elemental states of the materials,XPS measurements were conducted and the results are displayed in Figs. 2(e)—2(g). The major peaks of the Ni element can be found in the whole spectra(Fig.2(e)),and figure 2(f)illustrates that the Ni 2p3/2peak is consisted of Ni2+and Ni3+with different proportions,locating at 854.9 eV and 856.0 eV respectively.[31]The Ni 2p3/2peak of B-NMC at 855.98 eV is shifted to the higher binding energy of 856.04 eV for CA-NMC, which is due to the increase of valance state.[32]After the calculation of peak area,the area ratios of the Ni2+/Ni3+are 0.99,0.68,and 1.27 for C-NMC, CA-NMC, and B-NMC respectively. The lower value of the Ni2+/Ni3+ratio means the lower level of cation mixing,[22]and the increased amount of Ni3+on the surface of CA-NMC is mainly owing to the charge compensation mechanism of Ta5+,which is beneficial for the cation ordering and the electrochemical activity of the Ni-rich cathode.[22]Additionally,it is also believed that the coating layer could stabilize the surface of the CA-NMC sample,allowing more Ni3+to be retained. The XPS peaks in Fig. 2(g) prove the existence of Ta5+in CA-NMC.[33—35]

    Fig.2. XRD patterns of(a)C-NMC,CA-NMC,and B-NMC.(b)—(d)Magnified view of(003)peak(b),(006)/(012)peaks(c),and(018)/(110)peaks(d). XPS results for(e)survey scan,(f)Ni 2p3/2,and(g)Ta 4d CA-NMC.

    3.2. Room temperature test

    Electrochemical tests were carried out at room temperature to investigate the effects of surface modification. Figure 3(a) shows the charge—discharge profiles of C-NMC, BNMC,and CA-NMC samples at 0.1 C(1 C=220 mAh/g)with a voltage range of 3.0 V—4.4 V.No new platform is found after the coating of tantalum oxide,but the initial discharge capacity of C-NMC sample (217 mAh/g) is significantly lower than that of the other two samples(B-NMC:224 mAh/g;CANMC: 225 mAh/g). Meanwhile, the initial coulombic efficiencies of the C-NMC, B-NMC, and CA-NMC samples are~91%,~93%, and~93% respectively. After 150 cycles at 1 C,the CA-NMC sample shows the highest capacity retention of 83%,and the cycle performance of C-NMC is the poorest,as presented in Fig.3(b). And according to the results of rate capability, a similar conclusion can be drawn(Fig.3(c)).The reason for the superior performance of the CA-NMC sample might be that a beneficial and protective coating layer is formed via the coating process with further annealing, which can not only reduce the side reaction between the cathode material surface and electrolyte,but also will not hinder the surface transport of Li+. Although the surface modification is performed for the C-NMC sample,a coating layer without calcination may lead to slow Li+diffusion kinetics. Actually,the annealing treatment is a process that can repair the surface damage caused by the wet chemical coating method and also make the coating layer tighter,which will generate an interface with Li+transport activity. For the uncoated B-NMC sample,where there is no protective layer on the surface,poor performance is reasonable.

    The curves of differential capacity as a function of voltage are shown in Figs. 3(d)—3(f), containing several redox peaks.The nature of the redox species is concerned with the Ni2+/4+and Co3+/4+oxidation/reduction process accompanied with phase transition(hexagonal(H1),monoclinic(M),hexagonal(H2), hexagonal (H3)).[36]With the increased cycling number, the peak height at~3.7 V decreases, standing for a kinetic hindrance to Li+diffusion and loss of reversible capacity originated from structural damage.[37]For C-NMC,the H1—M peak heights decrease from 785 mAh/gV to 398 mAh/gV after 100 cycles corresponding to 49% height loss. This value is smaller for CA-NMC(43%)and B-NMC(48%). Moreover,a more pronounced peak shift is observed in C-NMC and BNMC samples, especially at the start of discharge (H3—H2),which indicates greater impedance growth during a prolonged cycle.[38]The C-NMC shifted 0.1 V to lower voltage,which is larger than CA-NMC(0.06 V)and B-NMC(0.06 V).The shift to higher polarization can be attributed to higher overpotentials for Li+intercalation/extraction into/from the host lattice.[39]A progressively increasing charge transfer resistance is a possible explanation for this observation and it will be further studied by subsequent EIS measurement.

    Fig. 3. Initial charge—discharge profiles at 0.1 C (a), cycling performance of 150 cycles at 1 C (b), rate capability (c) of C-NMC, B-NMC, and CA-NMC,dQ/dV profiles obtained by differentiating charge—discharge curves of every 25 cycles of C-NMC(d),CA-NMC(e),and B-NMC(f). All the measurements were tested under a voltage range of 3.0 V—4.4 V at room temperature.

    3.3. High temperature test

    High temperature tests at 45°C were also performed to study the electrochemical performance,as presented in Fig.4.Compared to the room temperature performance, the three samples exhibit higher initial discharge capacities (B-NMC:229 mAh/g; C-NMC: 229 mAh/g; CA-NMC: 233 mAh/g)and initial coulombic efficiencies (~95% of all) and lower potential polarization (Fig. 4(a)). It is worth noting that the CA-NMC sample can bring out the highest capacity, and as expected, the best cycle performance is found for the CANMC sample (Fig. 4(b)). The initial capacity is even parallel to those Ni90 or Ni95 cathode.[31,40]With the increase of cycling number, the voltage decay is noticeable for C-NMC and B-NMC samples,which is much slower for the CA-NMC sample (Figs. 4(c)—4(e)). As is well known, the voltage decay for cathode material is partly attributed to element dissolution and structure degradation, and surface coating is a very helpful strategy due to the formation of a protective layer,which is proved by the above results. The NRM suffers from the generation of microcracks, and the cross-sectional SEM images of B-NMC and CA-NMC samples after 150 cycles(Fig. S1 in supporting information) show that a stable coating layer is beneficial to reduce the risk of microcrack formation and the further side reaction and structure degradation,that means of the surface stability is very important for NRM.The allevation of microcrack formation is may be related to the buffering effect of coating layer and accommodation of internal stress.[41]The dQ/dVprofiles of the three samples obtained at 45°C exhibits a similar result with the room temperature test,as displayed in Figs.4(f)—4(h). Higher peak shift and polarization are observed for B-NMC and C-NMC samples,and the modified one is more stable. More specifically, the peak around 4.2 V during charging process decrease 76%and 70% for C-NMC and B-NMC after 100 cycles respectively,whereas, the CA-NMC sample presents moderate height loss of 59%. Meanwhile, the H3—H2 peak shifted 0.12 V for CNMC and B-NMC after 100 cycles,which is larger than CANMC(0.07 V).What is more,the plots of the average voltageversuscycle number illuminate a smaller voltage hysteresis of 0.18 V for the CA-NMC sample(Fig.S2 in supporting information), and the better thermal stability of the CA-NMC sample is confirmed by the DSC curves in Fig.S3 in supporting information. According to the results of room and high temperature tests, it can be concluded that surface coating of tantalum oxide is a useful way to improve the electrochemical performance of NRM.To gain a better understanding of the coating layer on electrochemical polarization during the charge/discharge process,the CV profiles of C-NMC,CA-NMC,and B-NMC samples for the initial three cycles at a scanning rate of 0.1 mV/s between 3.0 V—4.4 V are presented in Fig. 5. Several anodic peaks (oxidation) and cathodic peaks (reduction) resulting from multiple phase change with Li+extraction and insertion process consist in the CV curves. The potential gaps between the redox peaks of the CA-NMC sample marked in the CV profiles are narrower than that of both B-NMC and C-NMC samples, indicating that the polarization becomes smaller via functional coating of tantalum oxide,and the results are summarized in Table 2.

    Fig.4. Initial charge—discharge profiles at 0.1 C(a),cycling performance of 150 cycles at 1 C(b)of C-NMC,B-NMC,and CA-NMC;voltage-discharge capacity curves at different cycle[(c)—(e)]and dQ/dV profiles obtained by differentiating charge-discharge curves of every 25 cycle for C-NMC[(c)and(f)],CA-NMC[(d)and(g)]and B-NMC[(e)and(h)]. All the measurements were tested under a voltage range of 3.0 V—4.4 V at 45 °C.

    Table 2. The potential gaps between the first anodic and cathodic peaks marked in Fig.5.

    Electrochemical impedance spectroscopy was employed to offer insight into the resistance evolution for the three samples at a frequency ranging from 106Hz to 0.1 Hz after different charge/discharge numbers (Figs. 5(d)—5(f)).Z-view software is used to fit the impedance of the electrode and the simulated impedance parameters are listed in Table S1. The simulated circuit is presented in inset of Fig. 5(d). Two semicircles at the high and middle-frequency regions are consisted in each spectrum, and the semicircle at high frequency is related to the resistance of CEI(cathode electrolyte interphase)film on the surface,without a noticeable change in the whole cycling range. The C-NMC and B-NMC samples exhibit a larger radius for the semicircle at middle frequency, and the radius notably increases with the cycle number,implying that the charge transfer resistance is larger and the growth is more significant after the cycle. This may be the consequence of electrolyte consumption and accumulation of poorly conductive by-products, which hinder the fast transport of Li+.[42]However,the CA-NMC sample shows smaller charge transfer resistance and suppressed impedance growth, suggesting that the coating layer of tantalum oxide with further calcination is very stable and helpful in providing and maintaining fast charge transfer kinetics for NRM upon repeated cycle. Besides,the CA-NMC sample exhibits a larger initial Li+diffusion coefficient(D+Li),and the decrease of D+Liafter 150 cycles is insignificant compared with that of the other two samples(Table S1,supporting information).The higher D+Liof the CANMC sample is also responsible for its good electrochemical performance. The calculation process is presented in supporting information(Fig.S4,supporting information).

    Fig.5. Cyclic voltammograms of C-NMC(a),CA-NMC(b),and B-NMC(c)for initial three cycles between 3.0 V—4.4 V at a scanning rate of 0.1 mV/s.Nyquist plots of the first cycle(d),the 50th cycle(e),and 150th cycle(f)at charged state.

    In view of the great importance of CEI film on the electrochemical behavior of cathode materials, it is necessary to explore the chemical composition of CEI film generated on the surface of cycled samples with or without tantalum oxide coating. Therefore,the electrodes with different samples after 150 cycles at 45°C were disassembled and measured by XPS,and the results reveal that the CEI layer is mainly composed of C—O, C=O, C—F, and LiF as the decomposition products of electrolytes (Fig. 6). In the C 1s spectra (Fig. 6(a)), the peaks located at 284.8 eV correspond to the acetylene black,while the peak located at 290.8 eV is assigned to the CF2group from PVDF.[43]The peaks at 286.1 eV and 289.1 eV are associated with the C—O and C=O.[42,44]In Fig. 6(b),the O 1s spectra indicate the presence of C—O (531.7 eV)and C=O (533.1 eV).[42,43,46]And these species can be attributed to the decomposition of the carbonate electrolyte,leading to the formation of organic lithium alkyl carbonates(ROCO2Li) and lithium carbonates (Li2CO3). The tantalum oxide layer could partly reduce the electrolyte decomposition as seen from the proportion of surface chemical composition in Fig. 6(d). Compared to the B-NMC sample, CA-NMC possesses fewer organic components and a large proportion of inorganic Li—F. This may be originated from the formation of more beneficial electrolyte decomposition products LiF(685.03 eV) on the CA-NMC sample (Fig. 6(c)).[43,47]It has been reported that the LiF film is Li+conductive and can act as HF scavenger.[48,49]Overall,there are more by-products of electrolyte decomposition for B-NMC,and it implies that tantalum oxide modification could provide effective protection from side reactions.

    Fig.6. XPS spectra of(a)C 1s,(b)O 1s,and(c)F 1s. (d)The substance concentration for B-NMC and CA-NMC after 150 cycles at 45 °C between 3.0 V—4.4 V at 1 C.

    4. Conclusion and perspectives

    In the present study, ultrathin tantalum oxide coated LiNi0.88Mn0.03Co0.09O2cathode material was successfully prepared via wet chemical method followed by annealing treatment, and the structural characteristics and electrochemical performance were investigated. With the functional modification of tantalum oxide, the surface stability of LiNi0.88Mn0.03Co0.09O2is enhanced without hindering Li+transport and diffusion, and the side reaction with electrolyte is also suppressed. Furthermore, the capacity, rate and cycle performance are found to be improved along with low polarization and voltage decay. It is worth noting that the further calcination treatment is favorable for the formation of a tight coating layer and fast Li+transport interface. These findings of this study provide a beneficial and practical approach to protect the surface of nickel-rich cathode materials,which offers helpful insight for the optimization of cathode materials in the industry.

    Acknowledgements

    Zhenyu Zhang, Hong Li and Chaobo Li conceived and directed the project. Fengling Chen performed synthesis,performance testing and characterization. Jiannan Lin, Yifan Chen and Binbin Dong assisted the synthesis. Siying Tian,Dapeng Sun,Jing Xie and Chujun Yin collected and analyzed the XRD, SEM, and TEM data. Fengling Chen wrote the manuscript with assistance from all the authors.

    Project supported by the Key Laboratory Fund (Grant No.6142804200303)from Science and Technology on Microsystem Laboratory,the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences: Original Innovation Projects from 0 to 1 (Grant No. ZDBS-LY-JSC010), the Key Research and Development Project of the Department of Science and Technology of Jiangsu Province, China (Grant No. BE2020003), the Beijing Municipal Science and Technology Commission(Grant No.Z191100004719001),and the National Key Research and Development Program of China(Grant No.2017YFB0405400).

    猜你喜歡
    大鵬
    周鵬飛:大鵬展翅 跨界高飛
    Beating standard quantum limit via two-axis magnetic susceptibility measurement
    看圖紙
    三棱錐中的一個(gè)不等式
    Vorticity vector-potential method based on time-dependent curvilinear coordinates for two-dimensional rotating flows in closed configurations *
    Proton Beam Generated by Multi-Lasers Interaction with Rear-Holed Target
    李大鵬:打造縱向、橫向全域發(fā)展的蘇交科
    劉業(yè)偉、王大鵬設(shè)計(jì)作品
    非誠(chéng)勿擾
    AComparingandContrastingAnalysisofCooperationandPoliteness
    在线天堂中文资源库| 叶爱在线成人免费视频播放| 女人被躁到高潮嗷嗷叫费观| 亚洲成人手机| 午夜日韩欧美国产| 91精品国产国语对白视频| 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 亚洲欧美精品综合一区二区三区| 久久久久国产精品人妻一区二区| 91成人精品电影| 在线精品无人区一区二区三| 国产精品久久久久久人妻精品电影 | 国产精品偷伦视频观看了| 日本wwww免费看| e午夜精品久久久久久久| 国产成人免费观看mmmm| 亚洲综合色网址| 淫妇啪啪啪对白视频 | 精品国产超薄肉色丝袜足j| 成人av一区二区三区在线看 | 91国产中文字幕| 中文字幕人妻熟女乱码| 久久精品久久久久久噜噜老黄| 国产极品粉嫩免费观看在线| 国产成人一区二区三区免费视频网站| 中文字幕色久视频| 两个人看的免费小视频| 搡老熟女国产l中国老女人| 国产老妇伦熟女老妇高清| 老熟妇乱子伦视频在线观看 | 欧美97在线视频| 成年动漫av网址| 人人妻人人爽人人添夜夜欢视频| 成人亚洲精品一区在线观看| cao死你这个sao货| 男人操女人黄网站| 在线av久久热| 久久人人爽av亚洲精品天堂| 看免费av毛片| 亚洲中文日韩欧美视频| 色婷婷av一区二区三区视频| 一区二区日韩欧美中文字幕| 久久狼人影院| 十八禁网站免费在线| 日韩,欧美,国产一区二区三区| 日韩欧美一区二区三区在线观看 | 欧美精品人与动牲交sv欧美| 国产精品九九99| 欧美日本中文国产一区发布| 国产在线视频一区二区| 精品免费久久久久久久清纯 | 国产精品麻豆人妻色哟哟久久| 国产欧美日韩一区二区三 | 性少妇av在线| 亚洲综合色网址| 精品国产乱子伦一区二区三区 | 男女之事视频高清在线观看| 老鸭窝网址在线观看| 久久这里只有精品19| 久久久精品94久久精品| 日日爽夜夜爽网站| 国产精品熟女久久久久浪| 久久久精品区二区三区| 成人免费观看视频高清| 爱豆传媒免费全集在线观看| 国产精品免费视频内射| 日韩欧美一区视频在线观看| 一区二区三区精品91| 久久热在线av| 嫩草影视91久久| 色94色欧美一区二区| 热99re8久久精品国产| 国产成人精品久久二区二区免费| 天天影视国产精品| 亚洲七黄色美女视频| 国产亚洲欧美精品永久| 一级毛片女人18水好多| 亚洲国产精品成人久久小说| 99久久国产精品久久久| 亚洲av男天堂| 免费久久久久久久精品成人欧美视频| 夫妻午夜视频| 岛国在线观看网站| 亚洲国产毛片av蜜桃av| 人人妻人人爽人人添夜夜欢视频| a级毛片在线看网站| 老司机午夜福利在线观看视频 | 免费人妻精品一区二区三区视频| 18禁黄网站禁片午夜丰满| 老司机影院毛片| 纯流量卡能插随身wifi吗| 少妇粗大呻吟视频| 大型av网站在线播放| 久久国产精品人妻蜜桃| a 毛片基地| 国产成人系列免费观看| 欧美激情 高清一区二区三区| 男女免费视频国产| 国产成人免费观看mmmm| 日日摸夜夜添夜夜添小说| 午夜两性在线视频| 精品国产一区二区三区四区第35| 久久毛片免费看一区二区三区| 日韩,欧美,国产一区二区三区| 亚洲国产欧美一区二区综合| 亚洲天堂av无毛| 成人三级做爰电影| 黄片播放在线免费| 五月天丁香电影| 在线亚洲精品国产二区图片欧美| 欧美97在线视频| 91成年电影在线观看| 亚洲五月婷婷丁香| 大码成人一级视频| 99久久综合免费| 亚洲中文字幕日韩| 国产日韩欧美在线精品| 色94色欧美一区二区| 后天国语完整版免费观看| 免费不卡黄色视频| 国产老妇伦熟女老妇高清| 午夜福利一区二区在线看| 大香蕉久久成人网| 亚洲精品国产av成人精品| 男女无遮挡免费网站观看| 久久久久国内视频| 91精品伊人久久大香线蕉| 午夜福利视频在线观看免费| 亚洲一区二区三区欧美精品| 中文字幕高清在线视频| 考比视频在线观看| 一区二区av电影网| 水蜜桃什么品种好| 黑人欧美特级aaaaaa片| 国产在线一区二区三区精| 亚洲七黄色美女视频| 伦理电影免费视频| 国产野战对白在线观看| 十八禁人妻一区二区| 国产99久久九九免费精品| 亚洲精品一区蜜桃| 亚洲人成电影观看| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 亚洲av欧美aⅴ国产| 两个人看的免费小视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲五月婷婷丁香| 大片免费播放器 马上看| 欧美在线黄色| 欧美日韩福利视频一区二区| 日韩 亚洲 欧美在线| 亚洲av美国av| 日韩大片免费观看网站| 久久久久久免费高清国产稀缺| 精品久久久久久久毛片微露脸 | 午夜视频精品福利| 欧美精品一区二区大全| 丁香六月欧美| 欧美黑人精品巨大| 精品国产一区二区久久| 亚洲专区国产一区二区| 亚洲 国产 在线| 十八禁高潮呻吟视频| 国产精品久久久久久人妻精品电影 | 搡老乐熟女国产| 成人av一区二区三区在线看 | 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 亚洲国产欧美在线一区| 女人被躁到高潮嗷嗷叫费观| 电影成人av| 一区二区三区四区激情视频| 国产一区有黄有色的免费视频| 十八禁人妻一区二区| 午夜日韩欧美国产| 午夜福利,免费看| 免费不卡黄色视频| 老汉色∧v一级毛片| 国产av精品麻豆| 国产精品 欧美亚洲| 飞空精品影院首页| 一级,二级,三级黄色视频| 脱女人内裤的视频| 国产主播在线观看一区二区| 国产麻豆69| 日韩 亚洲 欧美在线| 欧美激情极品国产一区二区三区| svipshipincom国产片| 欧美日韩亚洲高清精品| 欧美精品人与动牲交sv欧美| 欧美乱码精品一区二区三区| 人妻一区二区av| 国产麻豆69| 国产精品免费视频内射| 欧美老熟妇乱子伦牲交| 99热国产这里只有精品6| 中文字幕精品免费在线观看视频| 亚洲精品国产精品久久久不卡| 两人在一起打扑克的视频| 啦啦啦视频在线资源免费观看| 乱人伦中国视频| 精品久久久久久久毛片微露脸 | 黄网站色视频无遮挡免费观看| 在线天堂中文资源库| 国产欧美日韩综合在线一区二区| www.熟女人妻精品国产| 不卡一级毛片| 久久精品人人爽人人爽视色| 一进一出抽搐动态| xxxhd国产人妻xxx| 亚洲精品国产区一区二| 蜜桃在线观看..| 亚洲美女黄色视频免费看| 两个人看的免费小视频| 欧美日韩黄片免| 另类亚洲欧美激情| 啪啪无遮挡十八禁网站| 欧美另类一区| 水蜜桃什么品种好| 午夜免费观看性视频| 亚洲av欧美aⅴ国产| 婷婷丁香在线五月| 国产老妇伦熟女老妇高清| 久久精品熟女亚洲av麻豆精品| 亚洲九九香蕉| 国产精品秋霞免费鲁丝片| 啦啦啦中文免费视频观看日本| 国产亚洲一区二区精品| 80岁老熟妇乱子伦牲交| 1024视频免费在线观看| 高清在线国产一区| 亚洲欧美激情在线| 宅男免费午夜| tocl精华| 精品少妇久久久久久888优播| 久久久久久人人人人人| 久久精品人人爽人人爽视色| 亚洲国产av新网站| 咕卡用的链子| 日韩视频在线欧美| 国产精品 国内视频| 伦理电影免费视频| 一区二区三区激情视频| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 黑人操中国人逼视频| 国产一区有黄有色的免费视频| 美女大奶头黄色视频| 夫妻午夜视频| tube8黄色片| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产日韩一区二区| 亚洲精品国产一区二区精华液| 亚洲 国产 在线| 这个男人来自地球电影免费观看| 五月天丁香电影| 视频区欧美日本亚洲| 欧美大码av| av在线老鸭窝| 一区二区av电影网| 欧美激情久久久久久爽电影 | 99热全是精品| 国产欧美日韩精品亚洲av| 亚洲天堂av无毛| 国产成人免费无遮挡视频| 天天躁夜夜躁狠狠躁躁| 亚洲欧美精品自产自拍| 亚洲中文字幕日韩| 岛国毛片在线播放| 爱豆传媒免费全集在线观看| 久久国产精品影院| 国产欧美亚洲国产| 久久久国产成人免费| 久久人人爽av亚洲精品天堂| 美女视频免费永久观看网站| 一个人免费在线观看的高清视频 | 老司机深夜福利视频在线观看 | 久久久久久亚洲精品国产蜜桃av| 手机成人av网站| 国产av一区二区精品久久| 视频在线观看一区二区三区| 国精品久久久久久国模美| 好男人电影高清在线观看| videosex国产| 欧美日韩亚洲国产一区二区在线观看 | 亚洲人成电影观看| 一二三四在线观看免费中文在| 欧美激情高清一区二区三区| 亚洲欧美一区二区三区久久| 亚洲av日韩精品久久久久久密| 国产日韩一区二区三区精品不卡| 欧美黑人精品巨大| 欧美另类亚洲清纯唯美| 日本欧美视频一区| 欧美黄色淫秽网站| 亚洲精品一区蜜桃| 精品国产一区二区久久| 成年av动漫网址| 久久性视频一级片| 午夜两性在线视频| 亚洲激情五月婷婷啪啪| 一边摸一边抽搐一进一出视频| 天天添夜夜摸| 王馨瑶露胸无遮挡在线观看| 亚洲成人国产一区在线观看| 久久久精品国产亚洲av高清涩受| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一区二区 视频在线| 黑人巨大精品欧美一区二区蜜桃| 黄色视频不卡| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 热99久久久久精品小说推荐| 国产成人av教育| 热99国产精品久久久久久7| 宅男免费午夜| 黄色片一级片一级黄色片| 一本大道久久a久久精品| 久久天堂一区二区三区四区| 午夜精品久久久久久毛片777| 亚洲精品国产av成人精品| 欧美变态另类bdsm刘玥| 美女主播在线视频| 午夜福利影视在线免费观看| 欧美变态另类bdsm刘玥| 美女主播在线视频| 99国产精品一区二区三区| 深夜精品福利| 桃花免费在线播放| 又紧又爽又黄一区二区| 亚洲国产日韩一区二区| 大码成人一级视频| 久久久久久免费高清国产稀缺| 国产精品av久久久久免费| √禁漫天堂资源中文www| 中文字幕制服av| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 亚洲国产精品999| av一本久久久久| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 精品国产一区二区三区四区第35| 欧美另类亚洲清纯唯美| 日韩精品免费视频一区二区三区| 亚洲国产中文字幕在线视频| 亚洲欧美色中文字幕在线| 亚洲免费av在线视频| 韩国精品一区二区三区| 性少妇av在线| 久久久国产欧美日韩av| 成年女人毛片免费观看观看9 | 国产免费一区二区三区四区乱码| 黑丝袜美女国产一区| 啦啦啦在线免费观看视频4| 男女床上黄色一级片免费看| 久久99热这里只频精品6学生| 亚洲,欧美精品.| 国产高清videossex| 久久狼人影院| 免费av中文字幕在线| 久久亚洲精品不卡| av一本久久久久| 少妇被粗大的猛进出69影院| 91大片在线观看| 精品亚洲成a人片在线观看| 亚洲激情五月婷婷啪啪| 一本久久精品| 91麻豆精品激情在线观看国产 | 亚洲色图综合在线观看| 久热这里只有精品99| 欧美日韩福利视频一区二区| 一级黄色大片毛片| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 欧美变态另类bdsm刘玥| 一本—道久久a久久精品蜜桃钙片| 秋霞在线观看毛片| 天堂中文最新版在线下载| 黄片播放在线免费| 一区福利在线观看| 99热网站在线观看| 亚洲国产欧美日韩在线播放| 最近最新免费中文字幕在线| 丰满少妇做爰视频| 成年人免费黄色播放视频| 老司机午夜福利在线观看视频 | 老司机影院成人| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 久久九九热精品免费| 久久久久国内视频| 午夜日韩欧美国产| 男女床上黄色一级片免费看| 亚洲精品国产色婷婷电影| av在线老鸭窝| 免费在线观看黄色视频的| 少妇 在线观看| 免费高清在线观看日韩| 国产成人欧美在线观看 | 亚洲五月婷婷丁香| 一本久久精品| 欧美黄色片欧美黄色片| 成人国产av品久久久| 热re99久久国产66热| 18禁国产床啪视频网站| 在线观看人妻少妇| av又黄又爽大尺度在线免费看| 韩国精品一区二区三区| 69精品国产乱码久久久| 秋霞在线观看毛片| 日韩视频一区二区在线观看| 久久久久国产一级毛片高清牌| 一级黄色大片毛片| 日韩人妻精品一区2区三区| 超碰成人久久| 一个人免费在线观看的高清视频 | 久久精品久久久久久噜噜老黄| 精品高清国产在线一区| 亚洲,欧美精品.| 91麻豆av在线| www.自偷自拍.com| 一区二区三区精品91| 国产精品一区二区精品视频观看| 国产亚洲精品久久久久5区| 亚洲欧美日韩高清在线视频 | 久久精品国产亚洲av香蕉五月 | 老鸭窝网址在线观看| 国产精品久久久久久精品电影小说| 国产伦理片在线播放av一区| 成人国产av品久久久| 亚洲va日本ⅴa欧美va伊人久久 | 久久精品亚洲av国产电影网| 欧美性长视频在线观看| 五月开心婷婷网| 国产一区二区 视频在线| 日韩精品免费视频一区二区三区| 国产成人免费观看mmmm| 18禁观看日本| 亚洲美女黄色视频免费看| 国产精品一区二区精品视频观看| 中文字幕另类日韩欧美亚洲嫩草| 国产成人免费无遮挡视频| 男女床上黄色一级片免费看| 日韩一区二区三区影片| 黑人操中国人逼视频| 精品免费久久久久久久清纯 | 国产精品熟女久久久久浪| av超薄肉色丝袜交足视频| 国产精品一二三区在线看| 久9热在线精品视频| 精品亚洲乱码少妇综合久久| 熟女少妇亚洲综合色aaa.| 亚洲成人国产一区在线观看| 精品高清国产在线一区| 99国产精品一区二区蜜桃av | 老司机在亚洲福利影院| 成年人黄色毛片网站| 亚洲精品国产区一区二| 桃红色精品国产亚洲av| 在线十欧美十亚洲十日本专区| 久久中文字幕一级| 脱女人内裤的视频| 欧美在线黄色| 中文精品一卡2卡3卡4更新| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 精品国产超薄肉色丝袜足j| 丝袜人妻中文字幕| 岛国在线观看网站| 女人久久www免费人成看片| 久久久久久久久免费视频了| 午夜影院在线不卡| 十分钟在线观看高清视频www| 亚洲中文字幕日韩| 国产精品 国内视频| 亚洲国产精品成人久久小说| 欧美xxⅹ黑人| 美国免费a级毛片| 亚洲av电影在线进入| 在线观看免费高清a一片| 蜜桃在线观看..| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 桃花免费在线播放| av天堂在线播放| 久久这里只有精品19| 免费在线观看影片大全网站| 如日韩欧美国产精品一区二区三区| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯 | 91国产中文字幕| 五月天丁香电影| 两个人免费观看高清视频| 亚洲av男天堂| 国产在线一区二区三区精| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| 欧美在线一区亚洲| 久久国产亚洲av麻豆专区| 欧美日韩亚洲综合一区二区三区_| 我要看黄色一级片免费的| 狠狠婷婷综合久久久久久88av| 精品少妇内射三级| 最近最新免费中文字幕在线| 亚洲专区中文字幕在线| 亚洲国产中文字幕在线视频| 国产在线免费精品| 丝袜脚勾引网站| 波多野结衣一区麻豆| 国产一区有黄有色的免费视频| 亚洲国产精品一区三区| 99国产综合亚洲精品| 亚洲专区国产一区二区| 久久久精品国产亚洲av高清涩受| 一区二区日韩欧美中文字幕| 国产免费视频播放在线视频| 69av精品久久久久久 | av有码第一页| 老司机在亚洲福利影院| 99热国产这里只有精品6| 久久狼人影院| 啦啦啦视频在线资源免费观看| 男女下面插进去视频免费观看| 美女视频免费永久观看网站| 中文精品一卡2卡3卡4更新| 亚洲久久久国产精品| 中文字幕精品免费在线观看视频| 中文欧美无线码| av在线老鸭窝| 人人妻人人添人人爽欧美一区卜| 国产伦人伦偷精品视频| 蜜桃国产av成人99| 久久久久视频综合| 亚洲男人天堂网一区| 我要看黄色一级片免费的| 久久久水蜜桃国产精品网| 超碰成人久久| 欧美人与性动交α欧美精品济南到| 天天影视国产精品| 高潮久久久久久久久久久不卡| 肉色欧美久久久久久久蜜桃| 国产无遮挡羞羞视频在线观看| 18禁裸乳无遮挡动漫免费视频| 美女高潮喷水抽搐中文字幕| 国产在视频线精品| 成人手机av| 无遮挡黄片免费观看| 黄色a级毛片大全视频| 99国产综合亚洲精品| 韩国高清视频一区二区三区| 中文字幕色久视频| 日本精品一区二区三区蜜桃| 欧美精品啪啪一区二区三区 | 久久久久国产精品人妻一区二区| 老司机午夜十八禁免费视频| 久久久欧美国产精品| 91九色精品人成在线观看| 亚洲人成电影观看| 欧美精品高潮呻吟av久久| 看免费av毛片| 亚洲黑人精品在线| 国产精品影院久久| 国产一区二区 视频在线| 精品免费久久久久久久清纯 | 亚洲欧美激情在线| 在线av久久热| 捣出白浆h1v1| 欧美日韩精品网址| 日本vs欧美在线观看视频| 黄片小视频在线播放| 高清欧美精品videossex| 国产av又大| 日韩 欧美 亚洲 中文字幕| 真人做人爱边吃奶动态| 18在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产av国产精品国产| 中亚洲国语对白在线视频| 精品亚洲成a人片在线观看| 我的亚洲天堂| 亚洲黑人精品在线| av天堂在线播放| 淫妇啪啪啪对白视频 | 亚洲国产欧美网| 桃红色精品国产亚洲av| 久久国产精品男人的天堂亚洲| 啦啦啦视频在线资源免费观看| 麻豆av在线久日| 丝袜喷水一区| av在线播放精品| 1024视频免费在线观看| 免费黄频网站在线观看国产| 啦啦啦在线免费观看视频4| 亚洲激情五月婷婷啪啪| 午夜福利视频在线观看免费| 久久99热这里只频精品6学生| 一级片免费观看大全| 涩涩av久久男人的天堂| 精品国内亚洲2022精品成人 | 久久中文看片网| 欧美精品啪啪一区二区三区 | 老汉色∧v一级毛片| 国产亚洲午夜精品一区二区久久| 亚洲性夜色夜夜综合| 精品国产超薄肉色丝袜足j| 亚洲avbb在线观看| 黄色视频,在线免费观看| 亚洲国产欧美日韩在线播放| 高清在线国产一区| 日本a在线网址| 99久久人妻综合| 国产成人av激情在线播放|