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

    Strain-modulated ultrafast magneto-optic dynamics of graphene nanoflakes decorated with transition-metal atoms?

    2021-09-28 02:18:28YimingZhang張一鳴JingLiu劉景ChunLi李春WeiJin金蔚GeorgiosLefkidisandWolfgangubner
    Chinese Physics B 2021年9期
    關鍵詞:李春張一鳴

    Yiming Zhang(張一鳴),Jing Liu(劉景),Chun Li(李春),3,?,Wei Jin(金蔚),Georgios Lefkidis,,and Wolfgang H¨ubner

    1School of Mechanics,Civil Engineering and Architecture,Northwestern Polytechnical University,Xi’an 710072,China

    2Department of Physics and Research Center OPTIMAS,Technische Universit¨at Kaiserslautern,P.O.Box 3049,67653 Kaiserslautern,Germany

    3Department of Mechanical Engineering,University of Manitoba,Winnipeg MB R3T 5V6,Canada

    4School of Physics and Information Technology,Shaanxi Normal University,Xi’an 710119,China

    Keywords:graphene nanoflakes,straintronics,spin dynamics,strain effect,first principles

    1.Introduction

    Graphene nanoflakes(GNFs),arbitrarily shaped graphene fragments that are finite in both dimensions,have emerged as new materials for improving electronic,spintronic,optical,and sensing devices.[1–4]They have potential applications due to the fact that they can be cut into a large variety of different shapes,due to quantum confinement and variation of the edges.[5,6]Mono-layer and multi-layer GNFs can be generated with different techniques.[7–9]When a graphene sheet is cut into fragments,the created edges strongly affect the electronic structure of graphene depending on the geometry of edges(zigzag and armchair).The presence of a zigzag edge state,accompanied by a large local density of states(LDOS)and localized spins,plays an important role in the electronic,magnetic and chemical activities to graphene edges and nanosized graphene.[10]Furthermore,GNFs with all zigzag boundaries,such as triangulene,[11]rhombic GNFs[12,13]and Clargoblet nanographene,[14,15]have been synthesized with several scalable synthetic methods and possess excellent magnetic properties.Such graphene nanoflakes with zigzag edges attract interests for their applications in spintronic devices.[16,17]

    However,the gapless nature of graphene[18]limits its applications in fabrication of novel electronic devices.A substantial amount of theoretical and experimental works have predicted that adsorption of transition metal atoms on graphene can effectively induce and tune band gaps of GNFs.[19–23]It is worth noting that the transition metal atoms(TM-atoms)doped graphene structures have been synthesized experimentally with satisfactory stabilities.[24–27]Previous studies have,without exception,expounded that doped TM-atoms can also significantly manipulate the electronic and magnetic properties of graphene.[28,29]Actually,both the low binding energy of TM-doped graphene and the strong bonds between the adjacent carbon atoms and the TM atom greatly affect the magnetic properties.For instance,the strong bonds between the carbon atoms and the TM-atoms lead to the localization of theπorbits near the Fermi levels;[30,31]3d atoms such as Co,[32]Cr,[33]and Fe[28]give rise to magnetic moments in monolayer graphene.Combining GNFs and TMatoms,the GNFs can provide ultrafast electronic charge transport and the magnetic elements act as information carriers,which could introduce new two-dimensional(2D)spintronic devices and applications.

    Strain is another effective method to modulate physical properties of 2D materials.[34–37]As is well-known,strain in graphene has been proposed as a tool for modifying its electronic structure and physical properties.For instance,both experimental[39,40]and theoretical[41]studies demonstrate that the ultrafast magnetization dynamics in fcc Ni can be manipulated with pressure.The essence of strain is to change the geometry and thus affect the electronic structure of the system.[42]In TM-doped GNFs,the magnetism of the structure is mainly induced by the interaction between the TMatoms and the adjacent carbon atoms.Therefore,altering the bonds between the TM-atoms and the GNF can effectively adjust the electromagnetic characteristics of the system and more intuitively reflect the internal mechanism of the strain modulation.

    These investigations stimulate us to explore the local deformation effects on the laser-induced spin dynamics on TM&GNFs(TM=Fe,Co,Ni)structures.We find that the spin-flip processes on TM&GNFs take advantage of the synergy effect of the TM-atom species and the position where the TM-atoms are attached to the nanoflake(TM-adatoms on the edge or vertex of the rhombic graphene nanoflake).Moreover,by applying strain along the C–Ni bonds,we have calculated the strain effect on the spin-flip process in the Ni&GNF structure.It shows that the increasing strain along the C–Ni bonds can significantly promote the spin-flip processes in Ni&GNF structures owing to the better selection rules for increased bond lengths.The present results offer a promising way to actively control spin,which can lead to more flexibility in the functionalization of future straintronic devices.

    2.Theory and method

    Our calculations are based on an ab initio quantum many-body model termed“Λprocess”.[43–45]In our previous studies,we demonstrated thatΛprocesses can act as exact microscopic spin-transition mechanisms on ultrafast,laser-induced,magnetization-dynamics scenarios.[46–50]This theoretical model also applies to magnetic molecules and low-dimensional carbon structures.[51–54]Especially,a longdistance ultrafast spin transfer model over one-dimensional carbon chain structure has been proposed very recently.[55]

    In our calculations,the Hamiltonian of strongly correlated systems is solved in four distinct steps.[56,57](i)The first approximate geometry optimization of the system at the Hartree–Fock level.This serves as a starting geometry for the subsequent more detailed geometry adjustments,and produces the electronic reference states for the subsequent post-Hartree–Fock calculations.

    (ii)A series of highly correlational post-Hartree–Fock calculations,typically with the symmetry-adapted cluster configuration interaction(SAC-CI)method.[58]Such a method represents at the moment the best compromise between accurate calculations of magnetic properties and reasonable computation time,and at the same time provides a balanced way of including both dynamical and static correlations.

    (iii)Perturbative inclusion of relativistic effects(notably spin-orbit coupling).This step is at the very core of our methods,because all our suggested laser-induced scenarios necessitate spin-mixed intermediate electronic excited states,which are optically coupled to both the initial and final states(Λprocesses).At this step we concurrently include an infinitesimal static external magnetic field,which allows us to energetically differentiate between the spin-up and spin-down states(Zeeman splitting).Here we can already distinguish the initial and final states for local spin-flip process.The idea of spin-flip process is that,to generate a fast transition between two almost degenerate states(within Zeeman splitting)with opposite spin directions,we need a laser excitation to a highly excited state,which is a spin-mixed state due to spin–orbit coupling,and subsequently undergoes a de-excitation to the state of opposite spin.

    3.Adsorption of TM-adatoms

    3.1.Structural optimization

    Since the rhombic GNFs have been synthesized,[12,13]there have been three types of TM-doped monolayer GNFs reported.(i)The TM atoms are doped on the boundaries of GNFs,[61–63]in which the TM atoms are doped on the edge of monolayer GNFs,and the TM atoms prefer staying in the plane.(ii)The metal impurities are involved by substituting carbon atoms in graphene,[64]that is,several C atoms are replaced by magnetic atoms,and the magnetic atoms will keep a certain distance out of the plane.(iii)The magnetic atoms are placed perpendicular to the GNFs(π-conjugated graphene),where the magnetic elements interact with the C atoms on the hexagon carbon rings,and the TM atoms prefer positions above the hollow sites of the hexagon carbon rings.[65–67]Among these three types of TM-doped monolayer GNFs,we investigate the structures with TM atom decorated on the GNFs boundaries,and construct the TM&GNFs(TM=Fe,Co,Ni)structures to study the spin dynamics of the two types of TM&GNFs(TM=Fe,Co,Ni),owing to their high degree of spin-density localization and narrow band electronic structures,thus allowing for coherent spin manipulation.As shown in Fig.1,when a magnetic atom is absorbed on the graphene nanoflake,the obtained structures can be classified into two types with respect to the connecting position:(i)TM-adatoms connected to one of the C atoms at the zigzag edges;(ii)TM-adatom on the vertex of the rhombic graphene nanoflake.

    Fig.1.Optimized(a)V-type and(b)VT-type TM&GNFs structures.(c)Sketch of the strain model for the VT-type TM&GNF.

    Before we investigate the spin dynamics of the two types of TM&GNFs(TM=Fe,Co,Ni),all the structures are fully relaxed to obtain a balanced structure,the stability of each structure is confirmed by the calculated normal mode with the absence of imaginary frequencies.The binding energy of the TM&GNFs structures are listed in Table S1 in the supplementary materials.After optimization,two types of structures are further studied,noted as V-type[Fig.1(a)]and VT-type[Fig.1(b)],respectively.Furthermore,we probe the geometry deformation after doping the magnetic elements,which is discussed in Section 1 in the supplementary materials.It is noteworthy that since Li et al.found reduced resistance and more active electromagnetic properties by removing the H atoms on the boundaries of graphene nanoflakes,[12]we do not add H atoms on the boundaries in the present calculations.

    3.2.Spin density

    Before calculating the time-dependent spin-flip process,it is necessary to check the spin density distribution of each equilibrium structure.Based on theΛprocess theoretical model,we need to identify the states with adequate spin localization,which could be the initial and final states of the spinflip processes.To this end,we calculate the spin density on each many-body state using Mulliken population analysis.It is noteworthy that,although this analysis only gives an estimate of the spin distribution in our structures,it is sufficient to detect which atoms and states carry the spin.

    Table 1 shows the spin density in the TM&GNFs(TM=Fe,Co,Ni)configurations for the lowest five triplet states,while the initial and final states will be selected from one of the spin localized triplet states.The energy levels of the TM&GNFs structures are listed in Fig.S1 in the supplementary materials.In addition,the spin configuration of the lowest five triplet states of V-and VT-type structures are listed in Tables S2 and S3,respectively.An important advantage of triplet over doublet ground states is that a triplet is built from two electrons:using triplet states we can separate spin and charge dynamics.[42,48,51,53]Compared to the spin distribution,we find that the spin is localized mostly on the TM atoms for state 1.Although the spin distribution in the ground state of Fe&GNF is not as high as in Co&GNF or Ni&GNF,it is yet by far higher than on the C atoms.Moreover,compared the spin distribution with different TM-atoms,degree of spin localization in Ni&GNF is higher than for Co&GNF and Fe&GNF.For instance,in the V-type structures for the ground state,the spin density of the Ni atom in Ni&GNF is 1.830,which is higher than the values in Co&GNF(1.782)and Fe&GNF(1.793).As for the VT-type structures,Ni&GNF also attains the leading position of spin distribution,indicating higher spin polarization and more pronounced magnetism.In addition,for the purpose of experimental realization we want the initial and final states to lie energetically as low as possible.Therefore,we choose the spin-down(|↓〉)and spin-up(|↑〉)substates stemming from triplet state 1 as the initial and final states in the further spin dynamics study.

    Table 1.Spin density(without SOC)in the TM&GNFs(TM=Fe,Co,Ni)configurations for the lowest five triplet states.

    3.3.Spin flip on TM&GNFs(TM=Fe,Co,Ni)

    Based on the spin density analysis of the two types of TM-atom doped structures(V-and VT-type)in the previous section,the triplet ground state fulfills the requirement of aΛ process.The spin-down(|↓〉)and spin-up(|↑〉)states stemming from the triplet ground state are chosen as the initial and final states for spin-flip processes.Moreover,the spin angular momenta of states|↓〉and|↑〉have the same absolute values and the opposite directions,satisfying the selection requirements of the initial state and the final state.[68]

    Figures 2 and 3 show the spin-flip processes obtained in the two types of configurations with three different TMadatoms,all of which complete within subpicoseconds with high fidelity above 96%.For each scenario,the upper panel depicts the occupation probabilities of the involved states and the lower panel depicts time-resolved expectation values of the spin components.

    Fig.2.Spin-flip processes achieved on the V-type(a)Fe&GNF,(b)Co&GNF,and(c)Ni&GNF structures.For each scenario,the upper panel depicts the occupation probabilities of the involved states(initial:black dashed,final:red solid,intermediate:dash-dotted),and the lower panel shows the time-dependent expectation values of the spin components.The insets give the envelopes of the optimized laser pulses.

    Fig.3.Spin-flip processes achieved on the VT-type(a)Fe&GNF,(b)Co&GNF,and(c)Ni&GNF structures.For each scenario,the upper panel depicts the occupation probabilities of the involved states(initial:black dashed,final:red solid,intermediate:dash-dotted),and the lower panel shows the time-dependent expectation values of the spin components.The insets give the envelopes of the optimized laser pulses.

    Comparing the spin-flip processes achieved in V-and VTtype structures,first,the spin-flip processes completed in the V-type structures are obviously faster than those in the VTtype structures.For a more intuitive comparison,Table 2 lists the parameters for each spin-flip process,where the FWHM(full width at half maximum of the laser pulse)is close to half the time needed to complete the flip processes.Here we compare the FWHM to illustrate the time for spin-flip processes.The FWHMs for Fe-,Co-,and Ni-decorated V-type structures are 92.27,284.27,and 245.28 fs,respectively,while for the VT-type structures,the corresponding magnitudes climb to 225.91,510.38,and 473.28,respectively,almost twice the V-type ones.Moreover,it is noteworthy that in both the Vand VT-type structures,the time required to achieve a spinflip process for three different TM-atoms increases in the order of Fe,Ni,and Co.Santos et al.[64]have pointed out that the Fe atom in graphene holds different magnetic properties compared with other transition metal atoms due to the complex band structure.For Co and Ni doped graphene nanoflakes,the electronic levels of Co&GNF and Ni&GNF contain more contributions of carbon character,which lead to the small oscillatory spin moments.Owing to the different bonding states between the GNFs and Fe,Ni,Co atoms,the spin-flip processes in TM&GNFs structures can be effectively modulated by different TM adatoms.

    Furthermore,we compare the energy of the laser pulse(Table 2).With the same adatoms,the energy of of the laser pulse is influenced by the synergy effect of the number of oscillation cycles in the population propagation,the number of involved intermediate states and the occupation of intermediate states.For instance,when the TM adatom is Fe,the energy of laser pulse is 0.687 eV for the VT type and 0.685 eV for the V type.For the Ni&GNF systems,the energies of the laser pulse are 0.419 eV(V-type)and 0.593 eV(VT-type),which is mainly caused by the occupation of the intermediate state(37%for V-type and 75%for the VT-type systems).In the Codecorated systems,the energies of laser pulse increase from 0.451(V-type)to 1.696 eV(VT-type).

    Table 2.Optimized laser parameters for each spin-flip scenario.Hereθandφdenote the angles of incidence in spherical coordinates,andγis the angle between the polarization of the light and the optical plane.FWHM is the full width at half maximum of the laser pulse.

    In order to further reveal the underlying mechanism,we also examine the highest occupied molecular orbital(HOMO)energy levels of the system(Fig.4).It is obvious that the electron clouds are localized around the Ni atom in the VTtype structures,while the V-type structures produce a completely different picture with the electron cloud distributed in the GNF.However,from the analysis of the spin density in the equilibrium structures,the spin density for the ground state is always localized on the Ni atoms in both V-and VT-type structures.This indicates that in the V-type structures,the applied laser pulse mainly affects the spin properties of the system without electronic interference.For the VT-type structures,both the spin and the charge are excited during the local spinflip processes obtained in the Ni atom,which could significantly slow down spin-flip processes.The lower energy and shorter time in the V-type structures demonstrate low-energy consumption and efficient logic operation processes.

    Fig.4.The electron clouds of the HOMO of(a)V-type and(b)VT-type equilibrium Ni&GNFs.The isosurfaces of the wavefunctions are red and green for positive and negative values,respectively.The isovalue is 0.02.

    4.Local deformation effect on Ni&GNF

    4.1.Structural reoptimization

    We have revealed the temperature-dependent,phonongenerated geometric distortions of the symmetry,the optical effects,and the laser-inducedΛprocesses[45,69,70]in the past.In the present study,we introduce the local deformation to the magnetic molecule systems(Ni&GNF),to actively tune the spin-manipulation scenarios.As mentioned in the above section,among the three TM-atom doped GNFs,the total energy of the Ni&GNF is the lowest,and thus is the most stable equilibrium structure.Based on the spin density analysis,Ni&GNF holds the leading position of spin distribution in both V-type and VT-type structures,indicating higher spin polarization and magnetism,so Ni&GNF is studied in the following part for the local deformation effect.

    In order to apply the strain on the C–Ni bond of the Ni&GNF systems,we firstly vary the coordinates of the C and Ni atoms along the directions of the applied strain proportionally to their values[Fig.1(c)].Then,the structures are reoptimized to obtain the deformed geometries.To fix the local distortions on the C–Ni bonds,we keep the coordinates of the relevant atoms as constants during the optimization.Except for the fixed atoms,all the other atoms are allowed to geometrically relax to construct stable structures.

    It is known that the theoretical C–Ni bond length is 1.82?A.However,in chemical compounds,the bond length often varies greatly.Many theoretical and experimental studies show that the C–Ni bond length generally lies between 1.8 and 2.0?A.[71–74]In our present work,the maximum C–Ni bond length is 2.103?A(in 10%strain),slightly larger than the values in previous works.Nevertheless,as a theoretical prediction,and considering the possibility of the experiment to the greatest extent,it is still worth testing.For instance,Santos et al.increased the spin moment in TM-doped GNFs by stretching the C–Fe bond by 9%,[64]which is just similar to the present work.Based on the above data,we apply a strain from?10%to 10%with 2.5%increment on the C–Ni bonds.That is,the C–Ni bond lengths range between 1.655?A and 2.023?A in the V-type structures and from 1.721?A to 2.103?A in the VT-type systems.The bond lengths and total energies of the optimized structures are listed in Table S5 in the supplementary materials.Moreover,as is expected,the local deformation can slightly break the geometric stability of Ni&GNF,since the total energy increases with increasing tensile or compressive strains.

    4.2.Local deformation effect on spin density

    Table 3 lists the Mulliken spin density distribution of the lowest five triplet states in the equilibrium and distorted Ni&GNF structures.The detailed information of the initial and final states(including their energies and spin expectation values)for the spin-flip processes of the strained VT-and Vtype Ni&GNFs are listed in Table S4 in the supplementary materials.In the V-type structures,the spin density for the ground state is always localized on the magnetic elements.With strain increasing from?10% to 10%,the degree of spin-density localization gradually increases.However,for the VT-type structures,when the strain applied to the C–Ni bond is?10%,state 3 is the lowest localized triplet state,and the spin density on the ground state is almost homogeneously localized on the carbon atoms.When the applied strain ranges from?2.5%to 10%,the Ni atom carries most of the spin density,and its value is monotonously increasing.Note that with the same strain increment,in the ground state,the spin density increases more quickly in the VT-type structures,indicating a more obvious strain effect.According to the theoretical requirement of aΛprocess,the spin density distribution for initial and final states should be ultimately localized on the magnetic center.Therefore,in the subsequent spin-dynamics calculations,we still use the spin-down and spin-up states stemming from the triplet ground state 1 as initial and final states for each distorted structure.That is,the applied strains from?10%to?5%will not be further studied for the VT-type structures.

    Table 3.Mulliken spin density in the V-and VT-type equilibrium and distorted structures for the five low-lying states,the maximum value for each structure is in bold.

    4.3.Spin-flip processes

    Based on the analysis of the spin density distribution in the previous section,we focus on the modulation effect of local deformation on the spin-flip process of Ni&GNFs.Figures S3 and S4 in the supplementary materials show the spinflip processes of V-type and VT-type Ni&GNF structures under different magnitudes of strain along the C–Ni bond.The corresponding parameters of the laser pulses are listed in Table 2.It is found that the time is obviously reduced with strain,but the rate of reduction is not identical.For comparison,we show the relationship between the duration(time consumption for spin-flip processes)and strain in Fig.5.It is shown that increasing the distance between magnetic elements and GNF is beneficial to the spin-flip processes,while the strain along the bond in the vertex of GNF(VT-type)can manipulate the spin dynamics behavior even more efficiently.

    The strain-promoted spin-flip processes in Ni&GNFs structures convincingly show that the Ni-GNFs hybridization is the key parameter in explaining the magnetic state of Ni doped GNFs.Since the Ni doped GNFs present non-magnetic character,by increasing the C–Ni bond lengths,we firstly change the local electronic properties in Ni&GNFs structures,and then influence the energy levels and energy gap(including the HOMO-LUMO energy gap)of the systems.Thus,we reduce the energy penalty for promoting the electrons between the initial/final state to the intermediate states.Furthermore,increasing the C–Ni distance can manually reduce the interaction between the Ni and carbon atoms,so that the magnetic properties of Ni atom get significantly enhanced,and thus affect the spin-flip processes in Ni&GNFs structures.

    To further explore the physical mechanism of the strain effect on the V-and VT-type structures,we calculate the HOMO(highest occupied molecular orbital)and LUMO(lowest unoccupied molecular orbital)energies of the strained Ni&GNF structures(Table S5 in the supplementary materials).The lowest 31 energy levels are calculated and listed in Fig.S2 in the supplementary materials.It is interesting that the HOMO and LUMO energies change gradually with the increasing bond length,while the HOMO going up and LUMO going down.The HOMO energies change faster than the LUMO energies,for instance,in the VT-type structures,the energies of HOMO increase from?4.764 eV to?4.478 eV,with 0.286 eV increment,while the energies of LUMO change 0.033 eV(from 0.584 eV to 0.551 eV),indicating that the HOMO energy level determines most of the virtual excitations,constituting the correlations in the many-body wavefunctions.Moreover,since the HOMO energies go up and the LUMO energies go down with increasing C–Ni bond lengths,it is obvious that the HOMO-LUMO energy gap decreases with increasing C–Ni distances.The decreasing HOMO-LUMO energy gap is accomplished with the increased electronic activity of the system,which could facilitate the transitions between the initial/final state and the intermediate states,thus prominently influence the spin dynamic behavior.

    At the same time,it is noteworthy that the strain effect of HOMO-LUMO energy gap(ΔEHOMO?LUMO)is more obvious in the VT-type structures.TheΔEHOMO?LUMOdecreases 0.249 eV(from 5.278 eV to 5.029 eV)with the strain increases from?2.5% to 10%,while in the V-type structures,the value decreases 0.242 eV(from 5.004 eV with?10%strain to 4.762 eV with 10% strain).It is obvious that the HOMO-LUMO energy gap decreases with a higher slope in the VT-type structures,with the same trend as the interaction between the strain and the duration(Fig.5).The strain along the C–Ni bonds in the VT-type structures can efficiently manipulate the HOMO-LUMO energy,thus dominating the spin dynamics behavior more effectively.

    From the above results,we find that the spin-flip process in Ni&GNF can be modulated by applying uniaxial strain.Nevertheless,the symmetric tensile and compressive strains do not result in symmetric spin-flip scenarios.The tensile strain can speed up the spin-flip process,while the compressive strain slows down the process.

    Fig.5.Duration for spin-flip processes of VT-type(red inverted triangles)and V-type(black regular triangles)Ni&GNFs under different strains.

    5.Conclusions

    In summary,the laser-induced ultrafast spin flip on the TM&GNFs(TM=Fe,Co,Ni)is investigated using theΛprocess-based theoretical model and ab initio calculations,and the spin-dynamics behavior induced by the variation of the C–Ni bond length is investigated in detail.The following conclusions can be drawn in this work.

    (i)The spin-flip processes on TM&GNFs are strongly related to the connecting positions.When the magnetic center is a Ni atom,it takes 940 fs and 500 fs in the VT-type and V-type structures,respectively.Similar situations can also be found in Fe&GNF and Co&GNF structures.

    (ii)Different magnetic atoms(Fe,Co,and Ni)can significantly change the spin properties,which results from the different binding states of Fe,Co,and Ni.The spin density redistributes with different magnetic centers and the spin flip in Fe&GNF is the fastest for both V-and VT-type structures.

    (iii)In the Ni&GNF structures,the spin-flip process can be modulated by the local strain along the C–Ni bonds.The time scale for spin-flip processes decreases with the increasing C–Ni bond length,owing to the better selection rules for increased bond lengths.Especially,the effect of local strain on the spin flip of VT-type structures is more pronounced than that of the V-type structures.When the C–Ni bond of the VTtype structure is stretched to 110% of its equilibrium length,the fastest spin-flip process is achieved within 80 fs.

    The theoretical results obtained in the present work provide the precursory guidance for further study of ultrafast magnetization control in GNFs,which can lead to more potential applications in future integrated straintronic devices.

    Acknowledgment

    GL gratefully acknowledges the hospitality of the Northwestern Polytechnical University,Xi’an,China,during the summer semester 2019.

    猜你喜歡
    李春張一鳴
    Fast identification of mural pigments at Mogao Grottoes using a LIBS-based spectral matching algorithm
    嘿!來點陽光音樂
    張一鳴的9條經驗
    誰是我們的朋友
    張一鳴的星辰大海
    商界(2019年5期)2019-05-21 03:42:08
    撒謊的“證人”
    藍盾(2017年11期)2017-12-25 21:27:30
    鬧洞房
    故事會(2016年17期)2016-09-01 08:14:32
    化學階段性創(chuàng)新訓練
    今日頭條張一鳴
    海峽姐妹(2015年5期)2015-02-27 15:10:54
    今日頭條張一鳴:不只靠技術
    国产视频首页在线观看| 男女无遮挡免费网站观看| 最后的刺客免费高清国语| av免费在线看不卡| 久久久成人免费电影| 视频区图区小说| 少妇被粗大猛烈的视频| 草草在线视频免费看| 亚洲av男天堂| 亚洲怡红院男人天堂| 国产白丝娇喘喷水9色精品| 欧美xxⅹ黑人| 伦理电影大哥的女人| 色综合色国产| 国产毛片在线视频| 男女无遮挡免费网站观看| 极品少妇高潮喷水抽搐| 2021少妇久久久久久久久久久| 少妇裸体淫交视频免费看高清| 一区二区av电影网| 肉色欧美久久久久久久蜜桃| 一个人看的www免费观看视频| 亚洲人成网站在线播| 最近2019中文字幕mv第一页| 久久午夜福利片| 色视频在线一区二区三区| 各种免费的搞黄视频| 日本黄大片高清| 人妻一区二区av| 亚洲精品一区蜜桃| 丝袜脚勾引网站| 欧美日韩视频高清一区二区三区二| 日本wwww免费看| 91久久精品电影网| 久久精品国产亚洲av涩爱| 国产成人精品福利久久| 久久久欧美国产精品| 妹子高潮喷水视频| 亚洲va在线va天堂va国产| 男女国产视频网站| 亚洲欧洲国产日韩| 两个人的视频大全免费| 18+在线观看网站| 久久久久久人妻| 下体分泌物呈黄色| 视频中文字幕在线观看| 国产成人免费观看mmmm| 国产亚洲5aaaaa淫片| .国产精品久久| 亚洲av中文字字幕乱码综合| 亚洲国产成人一精品久久久| videossex国产| 丰满乱子伦码专区| 久久精品夜色国产| 伦理电影免费视频| 一区二区av电影网| 这个男人来自地球电影免费观看 | 简卡轻食公司| 丰满迷人的少妇在线观看| 久久久成人免费电影| 日韩人妻高清精品专区| 最近中文字幕2019免费版| 国产久久久一区二区三区| 美女福利国产在线 | 美女cb高潮喷水在线观看| 在线观看国产h片| 成人国产麻豆网| 中文天堂在线官网| 亚洲中文av在线| 人妻 亚洲 视频| 亚洲不卡免费看| 国产精品99久久99久久久不卡 | 精品久久国产蜜桃| 成人国产麻豆网| av在线蜜桃| 亚洲国产日韩一区二区| 大又大粗又爽又黄少妇毛片口| 国产精品一二三区在线看| 午夜福利视频精品| 成人影院久久| 国产成人免费观看mmmm| 亚洲伊人久久精品综合| 少妇人妻 视频| 免费观看av网站的网址| 中文字幕久久专区| 性色av一级| 国产免费一区二区三区四区乱码| 大片电影免费在线观看免费| 街头女战士在线观看网站| 欧美日韩一区二区视频在线观看视频在线| 男女无遮挡免费网站观看| 亚洲av不卡在线观看| 亚洲国产精品成人久久小说| 国精品久久久久久国模美| 各种免费的搞黄视频| 哪个播放器可以免费观看大片| 寂寞人妻少妇视频99o| av网站免费在线观看视频| 精品人妻偷拍中文字幕| 欧美精品人与动牲交sv欧美| 深夜a级毛片| 日日撸夜夜添| 黑人猛操日本美女一级片| 国产综合精华液| 18禁裸乳无遮挡动漫免费视频| 91aial.com中文字幕在线观看| 国产真实伦视频高清在线观看| 波野结衣二区三区在线| 久久久久久久亚洲中文字幕| 人妻少妇偷人精品九色| 久久久久久久精品精品| 国产成人免费观看mmmm| 男女啪啪激烈高潮av片| 亚洲国产精品999| 新久久久久国产一级毛片| 国产 精品1| 天堂俺去俺来也www色官网| 最近最新中文字幕大全电影3| 99热这里只有精品一区| 99热这里只有是精品在线观看| 777米奇影视久久| 男女边吃奶边做爰视频| 日本爱情动作片www.在线观看| 啦啦啦中文免费视频观看日本| 菩萨蛮人人尽说江南好唐韦庄| 日韩人妻高清精品专区| 看免费成人av毛片| tube8黄色片| 国产v大片淫在线免费观看| 热99国产精品久久久久久7| 一区二区三区四区激情视频| 欧美精品国产亚洲| 日产精品乱码卡一卡2卡三| 成人高潮视频无遮挡免费网站| 国产精品一二三区在线看| 91aial.com中文字幕在线观看| 国产在线男女| 久久久久久人妻| 简卡轻食公司| 啦啦啦视频在线资源免费观看| 久久人人爽人人爽人人片va| 搡女人真爽免费视频火全软件| 精华霜和精华液先用哪个| av播播在线观看一区| 亚洲精品久久久久久婷婷小说| 日韩人妻高清精品专区| 亚洲av欧美aⅴ国产| 日本黄大片高清| 亚洲欧美精品自产自拍| 欧美日韩一区二区视频在线观看视频在线| av福利片在线观看| 老司机影院毛片| 亚洲在久久综合| 国产精品99久久久久久久久| 人妻一区二区av| 国产伦在线观看视频一区| 80岁老熟妇乱子伦牲交| 男女无遮挡免费网站观看| 国产大屁股一区二区在线视频| 97在线人人人人妻| 高清在线视频一区二区三区| 直男gayav资源| 亚洲美女搞黄在线观看| 熟妇人妻不卡中文字幕| 九九爱精品视频在线观看| 久久精品国产a三级三级三级| 亚洲不卡免费看| 国产乱人偷精品视频| 中文在线观看免费www的网站| 尤物成人国产欧美一区二区三区| 新久久久久国产一级毛片| 精品亚洲成国产av| 国产成人a区在线观看| 男人添女人高潮全过程视频| 亚洲三级黄色毛片| 美女中出高潮动态图| 国产伦精品一区二区三区四那| 午夜免费观看性视频| 午夜视频国产福利| 97超碰精品成人国产| 中文资源天堂在线| 国产探花极品一区二区| 岛国毛片在线播放| 三级经典国产精品| 午夜福利在线观看免费完整高清在| 日日撸夜夜添| 丰满迷人的少妇在线观看| 热99国产精品久久久久久7| 国产久久久一区二区三区| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 自拍偷自拍亚洲精品老妇| 色吧在线观看| 少妇的逼好多水| 青春草亚洲视频在线观看| 啦啦啦中文免费视频观看日本| 一级毛片我不卡| 国产美女午夜福利| 男人狂女人下面高潮的视频| 国产91av在线免费观看| 精品久久久噜噜| 97超碰精品成人国产| 黄色一级大片看看| 插逼视频在线观看| 天天躁夜夜躁狠狠久久av| 精品人妻视频免费看| 色吧在线观看| 91久久精品国产一区二区三区| 美女xxoo啪啪120秒动态图| 欧美成人a在线观看| 在线 av 中文字幕| 亚洲四区av| 大香蕉97超碰在线| 久久久久性生活片| 亚洲国产毛片av蜜桃av| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 亚洲综合色惰| 中文字幕精品免费在线观看视频 | 国产成人免费无遮挡视频| 日韩一区二区三区影片| 国产久久久一区二区三区| 黄色欧美视频在线观看| 午夜视频国产福利| 男女免费视频国产| 日韩大片免费观看网站| 五月玫瑰六月丁香| 噜噜噜噜噜久久久久久91| 一区二区三区四区激情视频| 日本欧美国产在线视频| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 午夜福利网站1000一区二区三区| 久久久欧美国产精品| 三级国产精品欧美在线观看| 国产亚洲午夜精品一区二区久久| 三级国产精品片| 亚洲精品一二三| 男女啪啪激烈高潮av片| 亚洲最大成人中文| 天天躁夜夜躁狠狠久久av| 亚洲一级一片aⅴ在线观看| 搡老乐熟女国产| 欧美最新免费一区二区三区| 性高湖久久久久久久久免费观看| 亚洲精品,欧美精品| 成人国产麻豆网| 少妇高潮的动态图| 欧美区成人在线视频| 亚洲,欧美,日韩| 国产白丝娇喘喷水9色精品| av在线蜜桃| 欧美bdsm另类| 我要看日韩黄色一级片| 成人无遮挡网站| 久久久久精品性色| 91精品国产九色| 欧美丝袜亚洲另类| 国产在线一区二区三区精| 久久毛片免费看一区二区三区| 国产精品.久久久| 91久久精品电影网| 色婷婷av一区二区三区视频| xxx大片免费视频| 国产免费视频播放在线视频| 在线免费十八禁| 免费黄网站久久成人精品| 国产欧美日韩精品一区二区| 九草在线视频观看| 边亲边吃奶的免费视频| 日本wwww免费看| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久av不卡| 麻豆精品久久久久久蜜桃| 少妇精品久久久久久久| 日韩在线高清观看一区二区三区| 成人亚洲欧美一区二区av| 婷婷色麻豆天堂久久| 中国三级夫妇交换| 欧美性感艳星| 男男h啪啪无遮挡| 建设人人有责人人尽责人人享有的 | 亚洲国产高清在线一区二区三| 久久久亚洲精品成人影院| 久久精品久久久久久噜噜老黄| 国产精品人妻久久久影院| 女人十人毛片免费观看3o分钟| 国产精品国产三级国产专区5o| 成人一区二区视频在线观看| 免费av不卡在线播放| 亚洲成人中文字幕在线播放| 日韩强制内射视频| 成年人午夜在线观看视频| 99热这里只有精品一区| 人妻制服诱惑在线中文字幕| 久久精品人妻少妇| 精品午夜福利在线看| 91久久精品国产一区二区成人| 韩国高清视频一区二区三区| 小蜜桃在线观看免费完整版高清| 全区人妻精品视频| 91aial.com中文字幕在线观看| 成人特级av手机在线观看| 大香蕉97超碰在线| 精品亚洲乱码少妇综合久久| 人妻制服诱惑在线中文字幕| 国产亚洲5aaaaa淫片| 国产亚洲91精品色在线| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 亚洲国产精品999| 天天躁夜夜躁狠狠久久av| 一本色道久久久久久精品综合| 久久 成人 亚洲| 国产一级毛片在线| 亚洲精华国产精华液的使用体验| 亚洲精品自拍成人| 又粗又硬又长又爽又黄的视频| 最后的刺客免费高清国语| 婷婷色av中文字幕| 久久av网站| 日韩中字成人| 最近中文字幕高清免费大全6| 国产美女午夜福利| 成人一区二区视频在线观看| 午夜福利在线观看免费完整高清在| 日本wwww免费看| 亚洲欧美日韩卡通动漫| 天堂俺去俺来也www色官网| 国产伦精品一区二区三区四那| 蜜臀久久99精品久久宅男| 免费大片18禁| 人妻系列 视频| 十八禁网站网址无遮挡 | 午夜福利影视在线免费观看| 91精品一卡2卡3卡4卡| 久久久久久久久大av| 国产亚洲一区二区精品| 黑丝袜美女国产一区| 久久ye,这里只有精品| 高清毛片免费看| 日本黄大片高清| 免费观看性生交大片5| 亚洲av在线观看美女高潮| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 直男gayav资源| 日本黄色日本黄色录像| 亚洲精品色激情综合| 熟女人妻精品中文字幕| 国产精品不卡视频一区二区| 免费高清在线观看视频在线观看| .国产精品久久| av免费观看日本| 久久精品久久精品一区二区三区| 美女视频免费永久观看网站| 少妇的逼水好多| 午夜福利在线在线| 老司机影院成人| 成人美女网站在线观看视频| 人妻夜夜爽99麻豆av| 午夜激情久久久久久久| 老司机影院成人| 成人一区二区视频在线观看| 国产免费一级a男人的天堂| 中文精品一卡2卡3卡4更新| 免费观看的影片在线观看| 国产一区有黄有色的免费视频| 夜夜看夜夜爽夜夜摸| 日韩大片免费观看网站| 一本久久精品| 最近最新中文字幕大全电影3| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 国产精品久久久久成人av| 国产色婷婷99| 人妻夜夜爽99麻豆av| 国产精品免费大片| 久久久久精品性色| a级一级毛片免费在线观看| 亚洲精品亚洲一区二区| 欧美另类一区| 99热6这里只有精品| 尾随美女入室| 观看av在线不卡| 美女福利国产在线 | 欧美高清性xxxxhd video| 三级经典国产精品| 一本—道久久a久久精品蜜桃钙片| av不卡在线播放| 久久精品熟女亚洲av麻豆精品| 日韩av不卡免费在线播放| 国产精品一区二区在线不卡| 韩国高清视频一区二区三区| 亚洲最大成人中文| 99久久精品热视频| 丰满少妇做爰视频| 亚洲第一av免费看| 国产永久视频网站| 成人18禁高潮啪啪吃奶动态图 | 亚洲人成网站高清观看| 寂寞人妻少妇视频99o| 一边亲一边摸免费视频| 狂野欧美激情性bbbbbb| 91狼人影院| 久久久久精品久久久久真实原创| videossex国产| 国产成人精品一,二区| 国产乱人视频| 青青草视频在线视频观看| 亚洲国产成人一精品久久久| 人体艺术视频欧美日本| 久久人人爽人人爽人人片va| 狂野欧美激情性bbbbbb| 国产伦在线观看视频一区| 欧美bdsm另类| 亚洲精品国产av成人精品| 欧美97在线视频| 亚洲精品久久久久久婷婷小说| 免费观看性生交大片5| 丝瓜视频免费看黄片| 日本免费在线观看一区| 国产亚洲午夜精品一区二区久久| 啦啦啦啦在线视频资源| 亚洲国产日韩一区二区| 高清日韩中文字幕在线| 五月天丁香电影| 久久亚洲国产成人精品v| 黄片wwwwww| av不卡在线播放| 亚洲人成网站在线观看播放| 亚洲熟女精品中文字幕| 性色avwww在线观看| 午夜免费观看性视频| 中文欧美无线码| 国模一区二区三区四区视频| 成年女人在线观看亚洲视频| 亚洲美女黄色视频免费看| 王馨瑶露胸无遮挡在线观看| 日韩视频在线欧美| 黄片无遮挡物在线观看| 欧美xxxx黑人xx丫x性爽| 精品视频人人做人人爽| 免费人妻精品一区二区三区视频| 欧美成人a在线观看| 精品视频人人做人人爽| 99热国产这里只有精品6| 国产美女午夜福利| 精品国产一区二区三区久久久樱花 | 99热国产这里只有精品6| 中文乱码字字幕精品一区二区三区| 网址你懂的国产日韩在线| 国产亚洲午夜精品一区二区久久| 久久久久久久亚洲中文字幕| 久久久亚洲精品成人影院| 美女福利国产在线 | www.av在线官网国产| 亚洲av.av天堂| 亚洲精品成人av观看孕妇| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 热re99久久精品国产66热6| 国产成人免费无遮挡视频| 国产一级毛片在线| 日本-黄色视频高清免费观看| 亚洲国产精品一区三区| 又爽又黄a免费视频| 日韩欧美 国产精品| 国产黄片视频在线免费观看| 一级爰片在线观看| 国产精品无大码| 2021少妇久久久久久久久久久| www.av在线官网国产| 久久久国产一区二区| 久久久久久久大尺度免费视频| 成年美女黄网站色视频大全免费 | 久久久久久久亚洲中文字幕| 国产黄频视频在线观看| 舔av片在线| 秋霞在线观看毛片| 久久国产精品大桥未久av | 亚洲精品色激情综合| 国产精品伦人一区二区| 身体一侧抽搐| 91精品国产国语对白视频| a级一级毛片免费在线观看| 男女下面进入的视频免费午夜| 国产精品一区www在线观看| 精品酒店卫生间| 亚洲精品国产成人久久av| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片 | 国产精品.久久久| 男女边吃奶边做爰视频| av线在线观看网站| 只有这里有精品99| 一本一本综合久久| 日韩,欧美,国产一区二区三区| 一边亲一边摸免费视频| 欧美成人午夜免费资源| 日韩一本色道免费dvd| av国产久精品久网站免费入址| 另类亚洲欧美激情| 七月丁香在线播放| 免费看av在线观看网站| 身体一侧抽搐| 天天躁日日操中文字幕| 国产精品人妻久久久久久| 日韩大片免费观看网站| 如何舔出高潮| 久久精品人妻少妇| 一本久久精品| 久久鲁丝午夜福利片| 国产亚洲午夜精品一区二区久久| av线在线观看网站| 国产午夜精品久久久久久一区二区三区| 在线 av 中文字幕| 亚洲成色77777| 一级毛片久久久久久久久女| 亚洲丝袜综合中文字幕| 久久99热这里只有精品18| 国产视频内射| 成人毛片a级毛片在线播放| av国产久精品久网站免费入址| 亚洲在久久综合| 夫妻午夜视频| 国产精品麻豆人妻色哟哟久久| 精品国产乱码久久久久久小说| 99热全是精品| 日韩欧美 国产精品| 一区二区三区乱码不卡18| 人妻制服诱惑在线中文字幕| 日韩强制内射视频| 夜夜爽夜夜爽视频| 日本vs欧美在线观看视频 | a 毛片基地| 亚洲第一区二区三区不卡| 午夜福利视频精品| 久久久久久久久久成人| 伦理电影免费视频| 成人综合一区亚洲| 色婷婷av一区二区三区视频| 嘟嘟电影网在线观看| 这个男人来自地球电影免费观看 | 丰满少妇做爰视频| 免费看光身美女| 我要看日韩黄色一级片| 免费人成在线观看视频色| 久久久久久伊人网av| 熟妇人妻不卡中文字幕| 深夜a级毛片| 天堂俺去俺来也www色官网| 亚洲精品日韩在线中文字幕| 亚洲欧美中文字幕日韩二区| 亚洲国产高清在线一区二区三| 男男h啪啪无遮挡| 少妇高潮的动态图| 午夜免费鲁丝| 亚洲人成网站在线观看播放| 一级片'在线观看视频| 在线观看一区二区三区激情| 久久久成人免费电影| 国产91av在线免费观看| av免费观看日本| 亚洲精品中文字幕在线视频 | 免费大片黄手机在线观看| 国产高清国产精品国产三级 | 精品一区二区免费观看| 国产欧美日韩一区二区三区在线 | 在线免费十八禁| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 成人亚洲精品一区在线观看 | 人妻系列 视频| 亚洲av中文av极速乱| 少妇 在线观看| 91午夜精品亚洲一区二区三区| 久久久久精品性色| 永久免费av网站大全| 热99国产精品久久久久久7| 精品熟女少妇av免费看| 中文字幕精品免费在线观看视频 | 国产在线免费精品| 亚洲色图综合在线观看| 日本黄色日本黄色录像| 日韩成人av中文字幕在线观看| 人妻一区二区av| 在线观看美女被高潮喷水网站| 在线精品无人区一区二区三 | 另类亚洲欧美激情| 久久人人爽av亚洲精品天堂 | 日本爱情动作片www.在线观看| 伦理电影大哥的女人| 最后的刺客免费高清国语| 日本黄大片高清| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区欧美精品| 国产精品99久久99久久久不卡 | 肉色欧美久久久久久久蜜桃| 我要看黄色一级片免费的| 精品国产露脸久久av麻豆| 免费观看的影片在线观看| 午夜福利在线观看免费完整高清在| 国产在视频线精品| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看| 大话2 男鬼变身卡| 最近的中文字幕免费完整| 在线观看一区二区三区| 亚洲人成网站高清观看| 成人18禁高潮啪啪吃奶动态图 | xxx大片免费视频| 人妻 亚洲 视频| 青青草视频在线视频观看| 精品酒店卫生间|