• <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网站的网址| 黄色欧美视频在线观看| 九色亚洲精品在线播放| 一级a做视频免费观看| 欧美日韩在线观看h| 精品少妇黑人巨大在线播放| 18禁观看日本| 插阴视频在线观看视频| 亚洲综合色网址| 这个男人来自地球电影免费观看 | 亚洲精品久久成人aⅴ小说 | 少妇精品久久久久久久| 下体分泌物呈黄色| 久久影院123| 日韩伦理黄色片| 国产男女超爽视频在线观看| 一个人免费看片子| 99热6这里只有精品| 国产av国产精品国产| 亚洲av男天堂| 91久久精品电影网| 人妻一区二区av| 99热这里只有精品一区| 最后的刺客免费高清国语| 又粗又硬又长又爽又黄的视频| 亚洲精品自拍成人| 免费大片黄手机在线观看| av黄色大香蕉| 人妻少妇偷人精品九色| 久久人人爽av亚洲精品天堂| 哪个播放器可以免费观看大片| 免费av不卡在线播放| 精品一品国产午夜福利视频| 亚洲国产欧美日韩在线播放| 最后的刺客免费高清国语| 黑人猛操日本美女一级片| 国产精品一国产av| 国产av精品麻豆| 精品人妻在线不人妻| 女的被弄到高潮叫床怎么办| 我要看黄色一级片免费的| 国国产精品蜜臀av免费| 婷婷色综合大香蕉| 一级毛片黄色毛片免费观看视频| 精品一区二区三卡| 啦啦啦啦在线视频资源| 18禁在线无遮挡免费观看视频| 夫妻午夜视频| 精品国产乱码久久久久久小说| 免费观看av网站的网址| 国产伦理片在线播放av一区| 久久精品国产亚洲av天美| 黄色视频在线播放观看不卡| 久久久久久久精品精品| 午夜激情久久久久久久| 欧美+日韩+精品| 国产爽快片一区二区三区| 一本一本综合久久| 视频中文字幕在线观看| 国产精品国产三级国产av玫瑰| 国产黄色免费在线视频| 久久精品人人爽人人爽视色| 婷婷色综合大香蕉| 你懂的网址亚洲精品在线观看| 国产精品国产av在线观看| 久久精品国产亚洲av天美| 亚洲欧美日韩另类电影网站| 亚洲综合精品二区| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品久久午夜乱码| 亚洲精品色激情综合| 欧美日韩在线观看h| 永久网站在线| 十八禁网站网址无遮挡| 夜夜爽夜夜爽视频| 国产精品蜜桃在线观看| 国产精品久久久久久精品电影小说| 91在线精品国自产拍蜜月| 亚洲综合色网址| 亚洲美女视频黄频| 最新的欧美精品一区二区| 免费av不卡在线播放| 熟女av电影| 精品一区二区三卡| 黑人巨大精品欧美一区二区蜜桃 | 成年人午夜在线观看视频| 国产黄片视频在线免费观看| 啦啦啦在线观看免费高清www| 精品国产乱码久久久久久小说| 天天躁夜夜躁狠狠久久av| 免费av不卡在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人午夜福利电影在线观看| 丝瓜视频免费看黄片| 建设人人有责人人尽责人人享有的| 欧美 亚洲 国产 日韩一| 国产免费视频播放在线视频| 亚洲欧洲精品一区二区精品久久久 | 99热这里只有是精品在线观看| 久久国内精品自在自线图片| 久久97久久精品| 国产亚洲午夜精品一区二区久久| 极品人妻少妇av视频| 国产一区二区三区av在线| 三级国产精品欧美在线观看| 夫妻性生交免费视频一级片| 久久av网站| 久久狼人影院| 精品亚洲成a人片在线观看| 蜜桃国产av成人99| 下体分泌物呈黄色| 我要看黄色一级片免费的| 一二三四中文在线观看免费高清| 97在线视频观看| 国产欧美另类精品又又久久亚洲欧美| 欧美人与性动交α欧美精品济南到 | 女人久久www免费人成看片| 中国美白少妇内射xxxbb| av专区在线播放| 一个人免费看片子| 亚洲av在线观看美女高潮| 香蕉精品网在线| 久久99热6这里只有精品| 两个人的视频大全免费| 国产精品一区二区在线观看99| 99国产精品免费福利视频| 免费人成在线观看视频色| 国产亚洲精品第一综合不卡 | 三级国产精品片| 日韩 亚洲 欧美在线| 一本一本综合久久| 亚洲综合精品二区| 精品国产国语对白av| 日日摸夜夜添夜夜爱| 免费大片18禁| 亚洲av电影在线观看一区二区三区| 一级片'在线观看视频| 女人精品久久久久毛片| 国产国语露脸激情在线看| 中文欧美无线码| av在线观看视频网站免费| 一级,二级,三级黄色视频| 寂寞人妻少妇视频99o| 99国产综合亚洲精品| 夫妻性生交免费视频一级片| 久久99蜜桃精品久久| 日本黄大片高清| 99视频精品全部免费 在线| 五月开心婷婷网| 欧美激情极品国产一区二区三区 | 精品视频人人做人人爽| 亚洲欧美色中文字幕在线| 国产亚洲av片在线观看秒播厂| 亚洲精品乱码久久久v下载方式| 飞空精品影院首页| 国产精品久久久久久精品古装| 99热6这里只有精品| 免费黄色在线免费观看| 午夜激情av网站| 狂野欧美激情性bbbbbb| 亚洲欧美中文字幕日韩二区| 国产老妇伦熟女老妇高清| 中文字幕精品免费在线观看视频 | 男男h啪啪无遮挡| 国产精品久久久久久av不卡| 毛片一级片免费看久久久久| 久久国产亚洲av麻豆专区| 久久国产亚洲av麻豆专区| 日日摸夜夜添夜夜爱| 卡戴珊不雅视频在线播放| 日日摸夜夜添夜夜爱| 国产精品久久久久久av不卡| 热re99久久国产66热| 美女福利国产在线| 热re99久久国产66热| 好男人视频免费观看在线| 热re99久久国产66热| 一区二区三区免费毛片| 一区二区三区免费毛片| 丰满少妇做爰视频| 高清欧美精品videossex| 日韩免费高清中文字幕av| 人妻人人澡人人爽人人| 夜夜骑夜夜射夜夜干| 亚洲av免费高清在线观看| 热99国产精品久久久久久7| 热99国产精品久久久久久7| 成年人午夜在线观看视频| 国产精品99久久99久久久不卡 | 美女国产高潮福利片在线看| 五月天丁香电影| 国产精品99久久99久久久不卡 | 久久久久国产精品人妻一区二区| 熟女人妻精品中文字幕| 欧美xxxx性猛交bbbb| 久久久精品区二区三区| 99热这里只有精品一区| 丰满饥渴人妻一区二区三| 日本黄色片子视频| 欧美日韩亚洲高清精品| 亚洲精品乱码久久久v下载方式| 亚洲精品一区蜜桃| 一级爰片在线观看| 男人添女人高潮全过程视频| 午夜视频国产福利| 少妇被粗大的猛进出69影院 | 成年人午夜在线观看视频| 最近的中文字幕免费完整| 欧美3d第一页| 亚洲不卡免费看| 日本91视频免费播放| 亚洲丝袜综合中文字幕| 久久亚洲国产成人精品v| 国产成人freesex在线| 国产一区二区三区综合在线观看 | 日日摸夜夜添夜夜添av毛片| 欧美日本中文国产一区发布| 久久女婷五月综合色啪小说| 久久人妻熟女aⅴ| 国产黄色免费在线视频| av免费观看日本| 精品久久久久久电影网| 自拍欧美九色日韩亚洲蝌蚪91| 国产视频首页在线观看| 亚洲色图综合在线观看| 好男人视频免费观看在线| 99久国产av精品国产电影| 97精品久久久久久久久久精品| 美女视频免费永久观看网站| 亚洲国产精品国产精品| 精品一品国产午夜福利视频| 亚洲无线观看免费| 色94色欧美一区二区| 国产成人精品无人区| 国产午夜精品久久久久久一区二区三区| 欧美bdsm另类| 亚洲精品色激情综合| 在线观看国产h片| 少妇猛男粗大的猛烈进出视频| 久久久午夜欧美精品| 如日韩欧美国产精品一区二区三区 | 欧美精品人与动牲交sv欧美| 欧美日韩在线观看h| 久久久久久伊人网av| 亚洲内射少妇av| 黄色怎么调成土黄色| 亚洲欧美清纯卡通| 狂野欧美白嫩少妇大欣赏| 国产一级毛片在线| 国产高清三级在线| kizo精华| 波野结衣二区三区在线| 国产永久视频网站| 日本午夜av视频| 成人国产av品久久久| 99热全是精品| 欧美成人午夜免费资源| 久久婷婷青草| 国产精品一区二区三区四区免费观看| 天堂中文最新版在线下载| 免费人妻精品一区二区三区视频| 亚洲伊人久久精品综合| 男女高潮啪啪啪动态图| 午夜福利视频精品| 亚洲国产精品成人久久小说| 久热久热在线精品观看| 免费少妇av软件| 国产一区二区三区综合在线观看 | 最近中文字幕2019免费版| 久久久午夜欧美精品| 日本vs欧美在线观看视频| 国产片内射在线| 精品少妇内射三级| 日韩 亚洲 欧美在线| 欧美成人午夜免费资源| 欧美变态另类bdsm刘玥| 国产白丝娇喘喷水9色精品| 熟女av电影| 精品人妻在线不人妻| 国产精品.久久久| 亚洲精品自拍成人| 国产成人91sexporn| 亚洲国产精品一区三区| 亚洲欧美清纯卡通| 国产精品久久久久久精品古装| 一级爰片在线观看| 亚洲精品美女久久av网站| 男的添女的下面高潮视频| 人人妻人人澡人人看| 国产视频首页在线观看| 春色校园在线视频观看| 妹子高潮喷水视频| 久久狼人影院| 老司机影院毛片| 日本与韩国留学比较| 大码成人一级视频| 国产免费视频播放在线视频| 少妇人妻精品综合一区二区| 在线观看人妻少妇| 国产一区二区三区综合在线观看 | 插逼视频在线观看| 一区二区三区精品91| 久久人人爽人人爽人人片va| 热re99久久国产66热| 亚洲精品第二区| 少妇精品久久久久久久| 九草在线视频观看| 日产精品乱码卡一卡2卡三| 最近最新中文字幕免费大全7| 天美传媒精品一区二区| 下体分泌物呈黄色| 精品亚洲成a人片在线观看| 国产成人freesex在线| 男女边摸边吃奶| 我的女老师完整版在线观看| 久久99热这里只频精品6学生| 欧美精品亚洲一区二区| 一个人免费看片子| videosex国产| 欧美日韩视频精品一区| 欧美另类一区| 国内精品宾馆在线| 亚洲图色成人| 国产一区二区三区综合在线观看 | 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 国产精品女同一区二区软件| 亚洲国产精品成人久久小说| videos熟女内射| 久久久久久久久久久免费av| 欧美精品一区二区免费开放| 欧美少妇被猛烈插入视频| 日本欧美国产在线视频| 97在线人人人人妻| 少妇高潮的动态图| 亚洲美女黄色视频免费看| 一级毛片电影观看| 男人爽女人下面视频在线观看| 亚洲国产精品国产精品| 国产欧美亚洲国产| 国产精品无大码| 亚洲精品国产色婷婷电影| 亚洲国产精品成人久久小说| tube8黄色片| 插逼视频在线观看| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 性色av一级| 91精品一卡2卡3卡4卡| 丝袜在线中文字幕| 91精品一卡2卡3卡4卡| 天天躁夜夜躁狠狠久久av| 亚洲av福利一区| 亚洲精品国产av成人精品| 久久久久久久久久久丰满| 黄色一级大片看看| 日本午夜av视频| 色视频在线一区二区三区| 三级国产精品片| 91精品一卡2卡3卡4卡| 亚洲一区二区三区欧美精品| 国产亚洲午夜精品一区二区久久| 久久久久网色| 男人添女人高潮全过程视频| 制服丝袜香蕉在线| 日韩三级伦理在线观看| 午夜福利影视在线免费观看| 97在线人人人人妻| 色网站视频免费| 久久久久久伊人网av| 18禁在线无遮挡免费观看视频| 亚洲人成网站在线播| 国产精品一区二区三区四区免费观看| 一区在线观看完整版| 久久青草综合色| 国产精品蜜桃在线观看| 精品熟女少妇av免费看| 草草在线视频免费看| a级毛片免费高清观看在线播放| 熟女人妻精品中文字幕| 午夜免费观看性视频| 久久久久久久久久久久大奶| 有码 亚洲区| 日韩 亚洲 欧美在线| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 国产精品一区二区在线不卡| 国产乱来视频区| 尾随美女入室| 欧美精品一区二区免费开放| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 岛国毛片在线播放| 草草在线视频免费看| 精品国产国语对白av| 91成人精品电影| 免费观看av网站的网址| 久久人人爽人人爽人人片va| a级毛片黄视频| av视频免费观看在线观看| 2022亚洲国产成人精品| 99九九线精品视频在线观看视频| 啦啦啦啦在线视频资源| 十八禁网站网址无遮挡| 国产极品粉嫩免费观看在线 | 成人二区视频| 九草在线视频观看| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 欧美97在线视频| a级毛片在线看网站| 亚洲国产精品一区三区| 如何舔出高潮| 国产男女超爽视频在线观看| 18禁裸乳无遮挡动漫免费视频| 中文乱码字字幕精品一区二区三区| 国产淫语在线视频| 国产精品 国内视频| 一级二级三级毛片免费看| 特大巨黑吊av在线直播| 久久国产精品大桥未久av| 一本—道久久a久久精品蜜桃钙片| 女性被躁到高潮视频| 国产精品国产av在线观看| 国产精品无大码| 大片免费播放器 马上看| 丁香六月天网| 久久久久久久久久久久大奶| 黑人欧美特级aaaaaa片| 三上悠亚av全集在线观看| 九九在线视频观看精品| 国产亚洲欧美精品永久| 国产乱来视频区| 女人久久www免费人成看片| 99久久中文字幕三级久久日本| 亚洲欧美中文字幕日韩二区| 亚洲国产精品999| 久久久久精品久久久久真实原创| 人妻人人澡人人爽人人| 日韩欧美一区视频在线观看| 黄片无遮挡物在线观看| 亚洲婷婷狠狠爱综合网| 在线 av 中文字幕| 一级毛片aaaaaa免费看小| 最新的欧美精品一区二区| 国精品久久久久久国模美| 91成人精品电影| 一级,二级,三级黄色视频| 免费观看性生交大片5| videos熟女内射| 亚洲精品亚洲一区二区| 国产高清有码在线观看视频| 极品少妇高潮喷水抽搐| 国产精品 国内视频| 91在线精品国自产拍蜜月| 在线观看人妻少妇| 国产黄色免费在线视频| 亚洲精品av麻豆狂野| 亚洲av成人精品一区久久| 2022亚洲国产成人精品| 亚洲成人手机| 久久久精品免费免费高清| 少妇被粗大猛烈的视频| 99久国产av精品国产电影| 国产伦理片在线播放av一区| 亚洲国产精品专区欧美| 简卡轻食公司| 九草在线视频观看| 啦啦啦在线观看免费高清www| 国产无遮挡羞羞视频在线观看| 国产av码专区亚洲av| 亚洲精品国产色婷婷电影| 老熟女久久久| 视频在线观看一区二区三区| 热re99久久国产66热| 两个人免费观看高清视频| 日韩精品有码人妻一区| 简卡轻食公司| 伦理电影大哥的女人| 国产精品三级大全| 久久韩国三级中文字幕| 精品视频人人做人人爽| 美女主播在线视频| 在线观看www视频免费| 国产黄色免费在线视频| 免费观看的影片在线观看| 午夜视频国产福利| 亚洲怡红院男人天堂| 狂野欧美激情性bbbbbb| 波野结衣二区三区在线| 男女国产视频网站| 精品人妻一区二区三区麻豆| 只有这里有精品99| 久久国内精品自在自线图片| 国产精品99久久久久久久久| 国产一区亚洲一区在线观看| 久热久热在线精品观看| 边亲边吃奶的免费视频| 午夜福利视频精品| 国产不卡av网站在线观看| av免费在线看不卡| 妹子高潮喷水视频| 亚洲国产精品一区二区三区在线| 久久青草综合色| 两个人免费观看高清视频| 亚洲国产精品专区欧美| 在线观看www视频免费| 国产免费福利视频在线观看| av在线播放精品| 成人午夜精彩视频在线观看| 99热网站在线观看| 欧美亚洲日本最大视频资源| 纵有疾风起免费观看全集完整版| 国产av国产精品国产| 赤兔流量卡办理| 亚州av有码| 精品99又大又爽又粗少妇毛片| 久久热精品热| 插阴视频在线观看视频| 亚洲国产日韩一区二区| 激情五月婷婷亚洲| tube8黄色片| 人人妻人人澡人人看| 国产精品一二三区在线看| 在线看a的网站| videossex国产| 免费看光身美女| 亚洲欧洲精品一区二区精品久久久 | 国产一区有黄有色的免费视频| av.在线天堂| 亚洲欧美日韩卡通动漫| 亚洲精品一区蜜桃| 精品熟女少妇av免费看| 久久韩国三级中文字幕| 色5月婷婷丁香| 多毛熟女@视频| 国产精品麻豆人妻色哟哟久久| 3wmmmm亚洲av在线观看| 免费黄色在线免费观看| av.在线天堂| 大片免费播放器 马上看| 麻豆成人av视频| 丰满乱子伦码专区| 日韩视频在线欧美| 七月丁香在线播放| 午夜精品国产一区二区电影| freevideosex欧美| 丝袜脚勾引网站| 大陆偷拍与自拍| av在线播放精品| 尾随美女入室| 精品久久久久久久久av| 在线观看人妻少妇| 一区二区av电影网| 欧美成人精品欧美一级黄| av在线观看视频网站免费| 国产精品一区www在线观看| av.在线天堂| 一区二区三区乱码不卡18| 日韩熟女老妇一区二区性免费视频| 天堂8中文在线网| 国产精品人妻久久久久久| 国产精品免费大片| 精品久久久噜噜| 成人国语在线视频| 十八禁高潮呻吟视频| 人人妻人人澡人人看| 晚上一个人看的免费电影| 街头女战士在线观看网站| 久久人人爽人人片av| a级毛片在线看网站| a级毛片黄视频| 三级国产精品欧美在线观看| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 精品卡一卡二卡四卡免费| 日韩熟女老妇一区二区性免费视频| av国产精品久久久久影院| 99re6热这里在线精品视频| 制服人妻中文乱码| 热99久久久久精品小说推荐| 黑人猛操日本美女一级片| 在线观看免费视频网站a站| 欧美另类一区| 考比视频在线观看| 多毛熟女@视频| 午夜视频国产福利| 97在线视频观看| 欧美变态另类bdsm刘玥| 多毛熟女@视频| 国产av一区二区精品久久| 夫妻午夜视频| 另类亚洲欧美激情| 日韩强制内射视频| 在线播放无遮挡| 一级二级三级毛片免费看| 色5月婷婷丁香| 青青草视频在线视频观看| 99re6热这里在线精品视频| 视频中文字幕在线观看| 91精品三级在线观看| 国产黄片视频在线免费观看| 免费大片黄手机在线观看| 久久 成人 亚洲| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 夜夜看夜夜爽夜夜摸| 国产精品熟女久久久久浪| 边亲边吃奶的免费视频| 最近2019中文字幕mv第一页|