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

    Characteristics of secondary electron emission from few layer graphene on silicon(111)surface

    2022-10-26 09:53:02GuoBaoFeng封國(guó)寶YunLi李韻XiaoJunLi李小軍GuiBaiXie謝貴柏andLuLiu劉璐
    Chinese Physics B 2022年10期
    關(guān)鍵詞:劉璐國(guó)寶

    Guo-Bao Feng(封國(guó)寶) Yun Li(李韻) Xiao-Jun Li(李小軍) Gui-Bai Xie(謝貴柏) and Lu Liu(劉璐)

    1National Key Laboratory of Science and Technology on Space Microwave,China Academy of Space Technology,Xi’an 710100,China

    2School of Computer Science and Engineering,Xi’an University of Technology,Xi’an 710048,China

    Keywords: secondary electron emission,graphene on silicon,numerical simulation

    1. Introduction

    Secondary electron emission from the material surface plays an essential role in multipactor,[1–4]which will seriously deteriorate the performances of spacecraft microwave components,[5–7]such as the fluctuation of cavity tuning, parameter coupling, waveguide loss, and phase constant, producing harmonics, leading to out-of-band interference and passive intermodulation products, and corroding the surface of components.[8–11]In recent years, since deep space exploration requires more miniaturization and high-power microwave components, the issue,i.e.the effective reduction of secondary electron emission and multipactor suppression,has become a significant challenge.[12–14]Nowadays,more attention has been focused on the surface coating treatment with a thin film, like low-dimensional carbon-based materials. In other fields,the studying of the electronic transport characteristics of 2D system heterostructures is also a critical topic,particularly the future application of these 2D materials in manufacturing microelectronic devices. Actually,considering that the integrated circuits keep miniaturizing on, physical limitations to effective carrier transport are experienced, such as secondary electron emission, which can dramatically hamper the performances of electronic devices. Furthermore, the investigation of interface formation for the reliable control of heteroboundary quality is of great significance in producing device-quality graphene with a realistic chance of entering into the market.

    As a typical 2D material, graphene presents a remarkable inhibitory effect on secondary electron emission.[15–17]The graphene-coated microwave components can remarkably enhance the multipactor threshold when suffering high power and present a vast application prospect in several areas,ranging from scanning electron microscopy to spacecraft materials engineering to hadrontherapy. Nevertheless,owing to the lack of an integrated theoretical analysis method, the mechanism of suppressing the secondary electron emission of graphene is still under debate. By investigating the secondary electron emission of graphene on copper, Caoet al.thought that the suppression of graphene is mainly due to the greater considerable surface work function.[18]With statistically analyzing the correlations between the maximum values of secondary electron yield and surface work function of a variety of materials, Wanget al.found that a material with a lower work function always comes with a less secondary electron yield,which is mainly because most of conduction band electrons can be act as emitted secondary electrons.[19]Many researches mainly focused on surface geometric construction and adsorption condition.[20,21]Actually,the secondary electron emission should be related to the effects of both inner scattering and surface work function. In theoretical research, although the classical formula, such as the Vaughan and Furman model,can satisfy the measurement results well,[22–24]those empirical formulas can hardly help reveal the physical mechanism from a microcosmic perspective. Monte Carlo simulation is the appropriate method for describing the micro-process of electrons inside the materials by considering elastic and inelastic scattering.[25–27]Maoet al.[28]and Pennet al.[29]introduced the dielectric loss of optical frequency to describe energy loss during electron–electron inelastic scattering. However, the critical optical constants, the extinction and refraction coefficients, rely on experiment data, making the micromechanism of secondary electron generation ambiguous, especially for 2D graphene-coated materials.[30]Without relying on experiment testing data,the first-principles calculation method can offer a practical method of obtaining the dielectric response characteristics of materials and the interface energy states of multilayer materials.[31–33]It can also offer a viable approach to figuring out the micro-mechanism of secondary electron emission from graphene.

    In this work, to reveal the micro-mechanism for coating graphene to suppress secondary electron emission, we propose a numerical simulation method of the first-principles based on Monte Carlo to investigate the characteristics of secondary electron emission from the graphene on a silicon surface. Electrons inside material suffering elastic scattering are calculated with Mott scattering cross-section. The full Penn algorithm (FPA) is used to calculate the inelastic scattering,whose energy loss function is obtained by the first-principles simulation.We analyze the energy loss function from the electron energy level transition via the partial density of states(PDOS).Simulations by using this method are compared with experimental results from recent researches and our measurement system. One to four graphene layers are set on the silicon(111)surface to investigate suppression characteristics of secondary electron emission. Considering the Van der Waals’interaction between multilayer graphenes,the electron density differences between few-layer graphenes are analyzed. In addition, considering that the high conductivity graphene layer will remarkably change the dissipation form of internal deposited charge,we also analyze the internal charging states in different layers coating situations.

    2. Simulation and measurement method

    Incident primary electron will suffer a series of collisions with material atoms and the stimulation of outer shell electrons. Part of inner electrons can escape from the surface after traversing the surface potential to form secondary electrons(SE),while other parts will remain inside the dielectric,forming an internal charged states. Here, a series of inner electrons collisions is simulated with the Monte Carlo method,and the critical energy loss function is derived via the firstprinciples method. The secondary electron emission simulations are compared with recently reported results and experimental results in our measurement system.

    2.1. Inner electron scattering

    The scattering process can be divided into elastic scattering and inelastic scattering based on whether the collision of inner electrons comes with energy loss. The elastic scattering always happens between electrons and the atomic nucleus because the nuclear mass is much larger than the electron mass. The differential cross-section of elastic scattering dσe(θ)/dΩcan be calculated by using the Mott scattering model via incident partial wave function and scattering partial wave function.[34]The total scattering cross-sectionσTecan be obtained by integrating scattering angleθbetween limits 0 andπas follows:

    whereλeis the mean free path of elastic scattering in units of cm andNis the atom density in units of cm-3.

    For the inelastic scattering of energetic electron–electron interaction (apart from electron–phonon and electron–polaron),part of the energy will transform from energetic electron to extranuclear electron due to energy level transition or free electron excitation.

    This energy loss process can be described via the dielectric function model. Then,Yanet al. developed an algorithm to derive this function from optical data,[35]which takes into consideration the energy loss ˉhωand momentum conservation.

    In this work,we use the full Penn algorithm(FPA)to derive the energy loss function.[36–38]Compared with the single pole approximation(SPA),the FPA considers both the single electron excitation and the plasmon excitation,making the Monte Carlo simulation of secondary electrons in a low energy range more accurate.

    The energy loss per length-dE/dsand mean free pathλincan be calculated as follows:

    In FPA,the energy loss function can be described as two contributions,i.e.the single electron excitation part and plasmon excitation part:

    Here,kis the Boltzmann constant, andTis the absolute Kelvin temperature. The transferred energy from electron to phonon is set to beWph= ˉhω=0.1 eV. In this simulation,the process of phonon annihilation is ignored because its probability is much less than that of phonon generation. For materials with electric polarization,free electrons with low energy may also be trapped by the polaron. Here,considering that the simulated materials are elementary substances, we ignore the interaction between electron and polaron.

    2.2. Energy loss function

    The reaction between energetic electron–extranuclear electrons can be described as a dielectric related energy loss function,which is always hardly obtained. Dielectric functionε(q,ˉhω)describes the material response to inner point charge related to momentum transformationqand energy loss ˉhω.In the limited optical conditionq →0, the dielectric function can be described asε(0,ˉhω)=ε1+iε2, whereε1andε2are the real part and the imaginary part of the dielectric function,respectively. Then the energy loss function can be described as

    wherenandkare the refraction index and the extinction coefficient of the material,respectively. Now,the primary method to obtain the energy loss function is based on experiment measurements ofnandk. When the electron energy is larger than 72 eV,thenandkdata can be calculated by component atoms based on Henke.[42]Acutely, not allnandkdata of material are available from the published experimental data in lower energy range<72 eV.This paper uses the first-principles method to derive the dielectric functionε1andε2of graphene on silicon(111)in the lower energy range.

    The three-dimensional(3D)lattice structure of graphene on silicon is shown in Fig.1. The silicon unit cell is set to be in the 2×2 inu–vplane in crystal orientation of (111), and its lattice parameters are (u,v,θ): 7.68017 ?A, 7.68017 ?A,120°. The surface graphene unit cell is set to be 3×3 inu–vplane with the lattice parameters are (u,v,θ): 7.379997 ?A,7.379997 ?A, 120°. The average lattice parameters of build layers are set to be(u,v,θ): 7.530083 ?A,7. 530083 ?A,120°with the mismatching rate 1.99%. The thickness of the vacuum layer between the layers is set to be larger than 20 ?A for building a less interaction vacuum boundary.

    Fig.1. Vertical view and front view of 3D lattice structure of graphene on silicon.

    The energy loss related dielectric function is derived with density functional theory (DFT) based on the first-principles method. In this work, the DFT is calculated via the Cambridge Sequential Total Energy Package (CASTEP, version 2019) software. We first use Broyden–Fletcher–Goldfarb–Shanno (BFGS) for the unit cell structure optimization and select the TS parameters for the van der Waals dispersion correction. Then,we calculate the energy band structure,density of states,and optical properties.[43]The exchange–correlation function in the Kohn–Sham equation is described by LDA(Local Density Approximation).[44,45]The electron wave function is expanded in basis plane wave. The energy cutoff is set to be 400 eV,the convergence precision of the SCF iterative process is set to be 5×10-6eV/atom. Thek-points of monolayer graphene and multilayer graphene are set to be 30×30×1 and 20×20×6. The actual minimum spacing reaches 0.016 1/?A. The pseudopotential is described with OTFG ultra-soft for calculating electron structure.[46]The imaginary part of the dielectric constant can be described via electron transition probability as follows:

    Ultra-soft pseudopotential produces an additional contribution to optical matrix elements that are also included in results. Then,based on the Kramers–Kronig dispersion relation,the real part of the dielectric constantε1(ω)can also be calculated. The energy loss function with electron energy loss ˉhωcan be derived from Eq.(7).

    When the loss energy is larger than a hundred eV,the inelastic scattering is mainly dominated by the interaction between free electrons and inner shell electrons. The inner shell interaction has an insignificant effect on the exchange–correlation field and is basically independent of each other for atoms. Here,in this simulation,for the situation of inner shell inelastic scattering,we can use the energy loss data by Henkeet al.[42]And the integrated energy loss function can be related to low energy loss from the outer shell and to high energy loss from the inner shell.

    2.3. Monte Carlo simulation

    A series of interactions,including the elastic and inelastic scattering after the primary electron irradiation on the material surface, can be described via the Monte Carlo simulation method. By using random numbers for describing the collective behavior of a large number of particles,Monte Carlo simulation is widely used in particle physics fields and engineering discharge simulation.[47,48]

    In general, the electron transport process inside the material includes the primary electrons and the generated secondary electrons. As shown in Fig.2, the simulation process begins with the initial parameters setting, and then calculates each step length via the MFP and a random number. Besides,the classification of scattering types and the judgement of the generation and emission of SEs generation and electrons emission are both determined by random numbers. The generated internal SEs will be treated as the initial electrons of the next step and added to the tracking electrons. Finally,by counting all the emission electrons,we can obtain the information about the secondary electron emission.

    Fig.2. Flow diagram of secondary electron emission for MC simulation.

    The azimuth angleφsatisfies the uniform distribution in a range of 0–π

    In FPA,the internal secondary electron generation should also be treated based on the single electron excitation part and the plasmon excitation part.[36]For the plasmon excitation part,the secondary electron is excited from the Fermi sea by loss energy ˉhω,the excitation probability of initial energy of an excited electronE′is

    Monte Carlo simulation needs to track each electron’s movement, and the step between each collision can be described via a random numberR1, Δs=-λlnR1. Here,λis the total mean free path covering all the mentioned processes,

    Here,Ξis the impulse function.

    During inelastic scattering,when the transforming energy is large enough for extranuclear electrons to escape from the nuclear constraint, the loss energy can generate internal free secondary electrons. We still need to follow their trajectories with the same method. For the electron–phonon process,based on Llacer theory,[49]the scattering angleθcan be described with random numberR5as follows:

    Here,EandE1are electron energy before and after the electron–phonon process,respectively.

    When the electrons move across the interface between different materials(including the inner interface of graphene–silicon and the surface interface of solid–vacuum),the energy and moving direction will be affected by the interface potential barrier(work functionWF)χ=EF+WF.Based on the solution of the incident and reflected electron wave functions,[50,51]the penetrating coefficient of electron crossing the interfaceTinterfacecan be given below:When the electron energyEand the normal angleφsatisfy the inequality:Ecos2φ <χ, the penetrating coefficientTinterface=0. As the electron escapes from the surface,it may be treated as the secondary electron emission.

    It should be pointed out that although the upper film thickness is comparable to or even less than the electron mean free path (MFP), this MC method may still be used for calculating the SE emission. For monolayer graphene,its thickness is about 0.3 nm, while the electron mean free path ranges from several angstroms to tens of nm, depending on energy. Even so,the Monte Carlo simulation method is still available for describing electron transport within the thin film.In Monte Carlo simulation,each step length is generated by a random number based on the MFP.The actual step length is not a constant of MFP but it is of a probability distribution in a specific range.Hence, the actual step length can be less than film thickness with scattering processes occurring. Owing to the fact that the Monte Carlo simulation relies on statistics of plentiful particle events,the macroscopic result still obeys the physical process.

    3. Result and analysis

    3.1. Monolayers of GoSi

    In the inelastic scattering process,both the scattering angle and internal free secondary electron excitation are determined by the energy loss function. Figure 3 shows the curves of density of states(DOS),dielectric constant and energy loss function calculated with the first-principles method as mentioned before. The curve of DOS of graphene is continuous at fermi energy due to the contribution of p orbital electrons,while for a typical semiconductor,the DOS of silicon presents a forbidden band above the fermi energy as shown in Fig.3(a).

    For the 2D graphene, the different hybridized orbitals show different behaviors in different directions, resulting in the discrepancy between in-plane and out-of-plane energy loss. For the electrons across the graphene layer, the anisotropy of energy loss should also be taken into consideration. In this simulation,we calculate the energy loss function of graphene along thecaxis (0, 0, 1) and theaaxis (1,0, 0), while the silicon energy loss function is treated as being isotropic. The composited contribution of the energy loss function can be integrated over all theta directions (-π/2 toπ/2),

    Here,θsurfis the surface normal angle,Psurf-ang(θsurf)is the calculated angle distribution of inner electrons when they reach the surface thin film from sub-material. The energy loss of sub-material silicon is treated as being isotropic. Considering that the energy loss in the plane in each direction is equal,the composited contribution of electron energy loss can be divided intocaxis part andaaxis part.

    For graphene,the real partε1and the imaginary partε2of dielectric constant in thecaxis part and theaaxis part cross around at the pointP(with energy 5.3 eV)andQ(with energy 15.2 eV),respectively,where the energy loss function present two peaks in Fig.3(b). The peak around 5.3 eV incaxis part results from the plasmon resonance in theπ-bond of graphene,while the 15.2 eV peak inaaxis part is due to the plasmon resonance in theσ+π-bond of graphene. Those two peaks of the energy loss function also result in the two peaks of inelastic scattering angle distribution. Peaks of imaginary partε2of graphene in thecaxis part(at 11.2 eV and 14.5 eV)and in theaaxis part(at 0.79 eV,4.0 eV,and 13.8 eV)correspond to different electron energy level transitions in Fig.3(a).

    Comparing the dielectric constants in Fig. 3(b) with the PDOS in Fig.3(a), we can find that peaks at 0.79 eV,4.0 eV,and 13.8 eV of the imaginary partε2ofaaxis part in Fig.3(b)may primarily result from the electrons transiting from valence band p orbital to conduction band p orbital in Fig.3(a). And,peaks at 11.2 eV and 14.5 eV of the imaginary partε2ofcaxis part may primarily result from the electrons transiting from valence band p orbital to conduction band s orbital. For silicon,the real partε1and imaginary partε2of the dielectric constant cross around the loss energy 17.4 eV,corresponding to a single peak of the energy loss function. The prominent peak of the imaginary partε2of silicon is around 7.0 eV,also denoting the dominated electron energy level transitions. After obtaining energy loss functions inaaxis part andcaxis part,we can calculate the composited contribution of the energy loss function based on Eq.(19)in Fig.3(c). On the whole,the single peak of silicon energy loss function presents much higher intensity thanπandσ+πpeak of graphene.

    Figures 4(a)and 4(b)show the energy loss function of Si and monolayer graphene by using the mentioned full Penn algorithm. The energy loss function is affected by energy loss,ˉhω,and momentum transferˉhq. For Si in Fig.4(a),the energy loss function forms the sharp Bethe ridge along with both energy loss and momentum transfer directions, while the Bethe ridge appears multimodal and gentler. The entire ˉhω–ˉhqcoordinate region includes the single electron excitation part and plasmon excitation part. When the momentum transfer equals 0, the energy loss function degrades into the optical limiting energy-loss function.

    Fig.3. Curves of(a)graphene density of states and(b)graphene dielectric constants(ε1 and ε2)versus energy,and(c)curves of graphene loss function versus loss energy.

    Fig. 4. Energy loss function of (a) Si and (b) monolayer graphene via full Penn algorithm.

    In the process of elastic scattering, only the motion direction of the energetic electron changes under the force between the electron and atomic nucleus, while both direction and energy of electron change during inelastic scattering. Figures 5(a) and 5(b) are electron scattering angle distributions at different electron energy of silicon and graphene. When the electron energy is low,such as 20 eV,the scattering angle presents two peaks (at~23°and~96°for silicon, at~34°and~148°for graphene). The secondary peak recovers and vanishes gradually with electron energy increasing, while the prominent peak increases and becomes sharp towards a slight scattering angle. It shows that the influence of elastic scattering on motion direction is less for high energy electrons.

    To verify the secondary electron emission by using the above-mentioned method, we compare the calculated secondary electron yield with experiment results from recent reports and our measurements. Figures 6(a) and 6(b) show our secondary electron measurement system and testing schematic diagram, respectively. The system consists of the preprocessing chamber(on the left)and the analysis chamber(on the right). The sample placed on the holder is covered with a three-layer shield grid to collect the emitted electrons more effectively in the analysis chamber. The background vacuum degree of the analysis chamber is maintained via a 3-level vacuum pump within 10-8Pa. The pulsed E-beam gun works in a range from 30 eV to 3 keV, with the emission current ranging from 50 pA to 50 μA. To avoid strongly charging state on the dielectric surface, the system sets a neutralization low energy E-beam gunin situand works in a range of 2 eV–100 eV. A positive +38 V bias voltage is added to the collection grid, and the collected current is obtained via a picoammeters Keithley 6487E.The secondary electron yield can be measured via the E-beam gun currentJPEand collected currentJCas SEY=JC/JPE, where the SEY stands for the secondary electron yield.

    Fig.5. Elastic scattering probability and angle of(a)silicon and(b)graphene,(c)variations of average scattering angle with electron energy of elastic and inelastic scattering in graphene and silicon, (d)variations of elastic and inelastic scattering probability with electron energy of silicon and graphene.

    Fig.6. (a)Diagram of secondary electron generation,(b)measurement system,and(c)comparison among secondary electron yields.

    In measuring the secondary electrons emitted from a dielectric material, factors of surface microstructure morphology, adsorption states, and charge accumulation will significantly affect the measured results,which results in experiment data varying in a significant range.Here,we compare this simulation results of secondary electron yield with the counterparts for recent references and our experiment measurements as shown in Fig. 6(c).[52,53]The deviations between simulation results and different experiments are existent, especially for monolayer graphene on silicon. Apart from measurement uncertainty, interface contamination and defect may both increase the deviation.

    3.2. Few layers of GoSi

    To further investigate the secondary electron emission characteristics in graphene layers,we also investigate the electron interaction with one to four layers of graphene on silicon. Since the potential field in multilayer graphene may be changed due to the exchange–correlation effect across the layers.We calculate the surface work function and the energy loss function in one to four layers of graphene on silicon. Here,in order to avoid the fringe effect in periodic boundary as much as possible, the sub-based material is set to be of four layers,and the vacuum layer is set to be 20-?A thick.The surface work functions and the average potential profiles of one to four layer GoSi are shown in Fig.7. Under the van der Waals force between silicon and graphene,for the 1 layer of GoSi as shown in Fig.7(a),the maximum average potential of interface rises to 5.43 eV in comparison with the potential in periodic structural gap of silicon. With the number of graphene layers increasing, the van der Waals force tends to be balanced, and this average potential difference decreases to 3.77 eV,3.16 eV,and 3.06 eV for 2–4 layers of GoSi respectively.The values of surface work functionWfof 1–4 layers of GoSi are determined to be 4.22 eV,4.37 eV,4.48 eV,and 4.41 eV,respectively,by the difference between Fermi energy and vacuum energy, under the combined effect of graphene potential and interface potential. Compared with Si (111) surface, as shown in Fig. 7(b),the work functionWfof the monolayer of GoSi decreases from 4.43 eV to 4.22 eV.With layer number increasing,work functionWfof GoSi increases to that of the graphite. The dotted line in Fig. 7(b) shows the ultimate value of the work function for a giant layer number graphene about 4.43 eV,which is close to the value reported previously.[54]

    In the microscopic perspective,the geometrical parameter and the electron density distribution may vary for multilayer graphene on silicon as shown in Fig.8.The GoSi system is optimized and its total energy reaches a minimum value in geometry optimization.To obtain the minimum energy of system,in the geometry optimization process, we relax silicon atoms in surface four layers and all of coating graphene C atoms. Since the strength of graphene CC bond is much stronger than that of the Si–Si bond,the C atom in a layer graphene presents a slight displacement. Owing to a weak Van der Waals interface force for the monolayer graphene on silicon,the minimum interface distance reaches 2.89 ?A, which is larger than the counterpart in other multilayer graphene situations. However, under the imbalanced Van der Waals forces from two sides for a fewerlayer GoSi(1–2 layers),the electron density difference of top silicon atoms(as the denoted district in Fig.5)presents a noticeable positive expansion. For the situation of more layers GoSi(3–4 layers),Van der Waals forces from two sides come to be more balanced,and the electron density difference deviation gradually vanishes.It is also the reason that the difference of interlayer peak potential decreases with graphene layers increasing as shown in Fig.8(a). And the interlayer distance of graphene layers keeps a constant of 3.30 ?A,while the interval of the upper graphene layer enlarges slightly, from 3.31 ?A to 3.33 ?A of 2 layers to 4 layers. Finally,under the combined effects of both the deviation of interface electron density difference and the variation of upper layer interval, the work function value of GoSi changes from 4.22 eV of 1 layer to that of graphite situations as the number of graphene layers increases.The apparent positive deviation of electron density difference leads the work function of 1 layer of GoSi to be less than that of non-coated Si(111)and graphite situation.

    Fig. 7. (a) Average potential profiles against normolized z, and (b) work function versus layer number.

    Fig.8. Electron density differences of multilayer(1–4 layers)graphene on silicon.

    Fig.9. Energy loss functions of one to four layer graphene.

    Figure 9 shows the electron energy loss function inside one to four layers of graphene. We can find that with the increase of layer number,theπpeak and theσ+πpeak present the tendency of moving toward higher energy and enhancing.Compared with the scenario in Ref. [55], the tendencies of the right-movement ofπpeak and the enhancement ofσ+πpeak basically agree with the simulation results. The value difference may result from the charge migration and structure variation due to the sub silicon interaction. Theπpeak energy moves from 5.3 eV of 1 layer to 6.3 eV of 4 layers, and theσ+πpeak energy moves from 15.2 eV of 1 layer to 16.6 eV of 4 layers. It indicates that the coupling effect of multilayer increases the out-of-planeπpeak energy and the in-planeσ+πpeak energy,and also enhances their intensities of energy loss.Because the interaction from sub-layer material on the surface decreases with distance increasing, the effect of layer numbers on the surface work function decreases rapidly. Based on previous researaches,[56,57]the work function reaches its stable value after layer number becomes larger than 4. For the body properties, the range affected by layer number is larger.For example, for the energy loss function, it is found that the EELS of 5-layer graphene is obviously different from that of graphite until the layer number is larger than 10.[55]

    Fig. 10. Variations of electron permeability coefficient Tsurf with electron energy for different number layers GoSi.

    Figure 10 shows the penetration coefficientTsurfwhen electrons pass through the solid-vacuum’s surface potential barrier, obtained from Eq. (18a). We can find that only the free electrons whose energy is larger than the surface work functionWFmay pass through the surface to form an outgoing secondary electron. Even so,part of electrons with energy larger thanWFmay still reflect back. For instance,the probability of 5 eV electrons reflecting back from 1 layer GoSi surface still reaches 19%. A low surface work function, such as 1 layer GoSi,still comes with a greater electron emission than Si(111)surface work function. Here in Fig.10,the electrons are assumed to pass through the surface vertically. For other directions, the electron energy in thex-axis direction should still be multiplied by cos2φ.

    After passing through the surface interface,the secondary electron emission present various characteristics in different layers of GoSi as shown in Fig. 11. Comparing with the original silicon (111) surface, the secondary electron yields from 0–4 layers of graphene on silicon decrease from 1.65 to 1.11–1.12. It is mainly because, on the one hand, as denoted in Fig. 10, the surface coating graphenechanges the electron emission potential barrier due to the fact that a lot of surfacearound low energy electrons still keep on escaping. Owing to a lower work function of 1-layer GoSi, the surface potential barrier of fewer layers (1–2 layers) graphene coatings cannot contribute to the suppression of secondary electron emission.Conversely, the energetic freemoving electrons can easily go to the inner with an apparently less elastic scattering angle in graphene. Electrons inside the graphene layers are accompanied with a more extensive in-plan energy loss with a prominent peak of 15.2 eV,while the energy loss peak of silicon is around 17.4 eV. In addition, the smaller scattering cross section of graphene also results in a lower collision probability,which means that the electrons tend to dissipate more deeply inside the graphene than inside silicon. So,under the two inverse factors of surface potential barrier and internal scattering process,the total emission secondary electron yield of 1-layer GoSi (the maximum secondary electron yieldδmax= 1.38)is obviously lower than that in the silicon (111) situation(δmax=1.65) as shown in Fig. 11. But, with the increase of graphene layer number, the recede of interface electron density difference deviation promotes the recovery of surfaceWFof GoSi,and also enlarges the range of electron scattering process in graphene. Those two factors both inhibit the secondary electron emission as shown in Fig. 11(a). The value ofδmaxdecreases from 1.38 to 1.12,and the primary energy of maximum secondary electron yieldEmaxalso recedes from 229 eV to 165 eV for 1–4-layers GoSi. In Fig. 11(b), we find that theδmaxacts as an exponential function of the graphene layers numberNlayer,δmax(Nlayer)=a·exp((1-Nlayer)/b)+δlim.Here,ais the parameter related to coating material, which is fitted as 0.28, andbis the parameter related to single-layer thickness,which is fitted as 1.34. Theδmaxis indicated as the limit value of the maximum secondary electron yield increasing with layer number increasing,and it is 1.06 in this case.

    Fig.11. (a)Curves of secondary electron yield versus primary energy,and curves of maximum SEY and the Max SEY energy versus number of layers from 0-to 4-layer graphene on silicon(111)surface.

    Under the two effects,i.e.the surface work function barrier and the inner scattering intrinsic properties, the graphene layer can effectively reduce the emitted secondary electrons.The dominant effect causing such a reduction in the emission process remains controversial. Researchers always use different factors to explain phenomena,thus drawing conflicting conclusions sometimes. In order to more directly demonstrate which of those two factors(the surface work function and inner intrinsic properties)makes a dominant contribution to the secondary electron yield, we compare the situation with the case of four GoSi samples with different layers, and the results are shown in Fig.12. The contribution of work function to GoSi SEY suppression is defined as

    where Δδmaxis defined as the reduction value by comparing GoSi with Si(111),i.e.Δδmax=δmax|GoSi-δmax|Si(111).The subscript WFGoSi=WFSiin Δδmax|WFGoSi=WFSidenotes that the GoSi work function is supposed to be the same as Si (111). For the situations of 1–2-layers GoSi, the values of Δδmax|WFGoSi=WFSi(0.32 and 0.50)are larger than true values of Δδmax(0.28 and 0.493), thus, the values ofCWFare negative,i.e.-13.2%and-1.4%for SEY suppression. With the increase of graphene layer number (3–4 layers GoSi), although the Δδmax|WFGoSi=WFSiis larger than Δδmax,the values ofCWFare merely 1.0%and-0.3%for the SEY suppression.It suggests that the surface potential barrier of fewer graphene layer GoSi(1–2 layers)may inversely promote the secondary electron emission, which hardly contribute to SEY suppression even for more layers of GoSi(3–4 layers). The inhibition effect of graphene on SEY is mainly by inducing the deeper electrons in the scattering process. Besides, what should be noted is that for other base materials, various Van der Waals interface forces may result in the effects on work function to varying degrees or even leads to the converse effects on work function.

    Fig. 12. Contributions of work function to SEY suppression in different layer numbers of GoSi.

    Fig. 13. Angle distributions from different layer number GoSi (a) before and(b)after emission.

    Since moving electrons exhibit various behaviors when residing inside the materials with different intrinsic properties or passing through the interface, the distribution of deposited electrons also presents various characteristics inside different layer-number GoSi,which can also explain the corresponding various SEYs from the inside viewpoint. Figure 14(a) shows the true SE emission proportions and internal deposited SE distributions in depth of 2-layer GoSi under 100 eV, 300 eV,500 eV, and 800 eV electron-beam irradiation. We can find that most of the true SE emissions exit from 20-?A depth.Hence,the multiplication of internally deposited electron distribution and emission electron proportion in this region correlate with the SEY.As shown in Fig.14(a),when the primary energy is 300 eV, the emission electron proportion becomes flatter in depth. It is mainly because the situation with primary energy 300 eV comes with a larger inelastic scattering section around the surface than the situation with higher 800 eV.Meanwhile, it can generate more energy loss than the situation with the lower 100 eV as indicated by the dash lines in Fig.14(a). It results in the SEY for 300 eV being larger than for 100 eV and 800 eV.In addition,comparing with graphene layers and sub-silicon, the internal electron distribution appears to have a remarkable sharp rise,resulting from the difference in scattering process and the defect and potential barrier around the interface.

    For different layer-number GoSi as shown in Fig. 14(b)in primary energyEPE=300 eV, as the number of layers increases,the solid red line(4 layers)decreases with depth faster than the blue and black lines (1–2 layers), which means that more SEs are emitted from the surface region of the biggerlayer-number graphene. Meanwhile, for the graphene layer with a smaller elastic scattering cross section,deposited SEs in graphene region decrease with layer number increasing from 1 to 4,indicating the excited SE quality in graphene region decreases. Consequently, under the combining effects of those two distributions(the emission SE position and the excited SE position), internal SEs that can reach the sample surface become less for a thicker graphene-coated GoSi. Moreover, a larger surface work function further results in less SE emission.

    Fig.14. Internal electrons distributions of graphene on silicon for(a)under different incident energies and(b)different layer numbers of GoSi.

    In addition,when calculating interface with DFT,in order to reduce the cell size,a mismatch is alway created due to the uniformity of the material lattice constants on both sides of the interface, which may have an influence on the absolute value of material properties. Nevertheless,in this study,the intrinsic mechanisms and characteristics are still kept consistently since the cases considered are all under the same mismatch parameters. Here, we only investigate the situation of graphene on silicon (111) crystal face. As for other crystal faces, such as(110) or (100), the method proposed in this study is still applicable to the calculating of the SEY.It should be noted that the interface matching requirement may make the calculation much longer due to a larger size unit cell.

    Here, we choose a large atom-density crystal face at the interface Si (111) as a typical semiconductor substrate for coating graphene layers in this research. Although owing to the different surface atom densities and the interface matching state, the work function, and interface potential could be different from those of other directional crystal surface, this method of calculating SE emission is also appropriate to other situations. In practice,the substrate silicon wafer is easily oxidized to form a dangling bond with the O atom. Hence, for forming a pure silicon surface,the silicon substrate should first experience the process of outgassing at a temperature higher than 600°C in an UHV environment less than 10-10mbar(1 bar=105Pa) UHV for 24 h. Then, after repeatedly flashheating in 5×10-9-mbar UHV to over 1200°C, the silicon substrate could present the equilibrium structure of Si (111)without dangling bond. Actually, the dangling bond will change the surface binding energy and the interface potential.A thicker interface oxide layer can also be treated as an extra layer to influence the secondary electron emission bychanging the inner electron collision process. Those related characteristics will be investigated in further researches.

    4. Conclusions

    In this work,we investigated the electron interaction and SE emission of few layer graphene coated on silicon with a first-principles method based on Monte Carlo numerical simulation. The electron interaction inside the GoSi is calculated with considering elastic scattering via the Mott model,and inelastic scattering via the Full Penn Algorithm of the dielectric function model. The energy loss function is calculated via the first-principles method based on the energy level transition theory. From this study some conclusions can be obtained below. (i) With reverse effects of the surface potential barrier and internal scattering process on SE emission suppression,monolayer graphene can still reduce the maximum SEY to 1.38 from the original silicon 1.65.(ii)When the layer number of graphene increases from 1 to 4,because of the coupling effect of multilayers,both the out-of-planeπenergy loss peak and the in-planeσ+πenergy loss peak shift toward higher energy. The maximum SEY further decreases from 1.38 to 1.12 of 1 to 4 layers,behaves as a negative exponential function of graphene layer number.(iii)Owing to the positive deviation of electron density difference at the graphene–silicon interface,the surface work function of monolayer GoSi is pulled down to 4.22 eV from Si (111) 4.43 eV. The work function contribution of monolayer GoSi to the SEY suppression reversely reaches-13.2%, and even for 3-layer GoSi, the work function contribution is just positive 1.0%. The dominant factor of the SEY inhibition for GoSi is the mechanism by inducing electrons more deeply in the scattering process. This research can contribute to the understanding of the microscopic mechanism of electron emission from the coated graphene and has reference significance in suppressing the secondary electron multiplication effect in many engineering fields.

    Acknowledgements

    Project supported by the National Natural Science Foundation of China (Grant Nos. 61901360 and 12175176), the Natural Science Foundation of Shaanxi Province, China(Grant No. 2020JQ-644), and the Scientific Research Projects of the Shaanxi Education Department, China (Grant No.20JK0808).

    猜你喜歡
    劉璐國(guó)寶
    由年齡推斷屬相
    磨刀不誤砍柴工
    誰(shuí)算得對(duì)
    劉璐規(guī)劃設(shè)計(jì)作品
    劉璐規(guī)劃設(shè)計(jì)作品
    “國(guó)寶”競(jìng)拍會(huì)
    《珍惜點(diǎn)滴》
    漫話國(guó)寶
    我們是國(guó)寶
    國(guó)寶求學(xué)記
    22中文网久久字幕| 最好的美女福利视频网| 国内精品一区二区在线观看| 国产亚洲精品综合一区在线观看| 一级黄色大片毛片| 男女啪啪激烈高潮av片| 免费观看在线日韩| 天天一区二区日本电影三级| 亚洲精品乱码久久久v下载方式| 一本精品99久久精品77| 免费一级毛片在线播放高清视频| 久久久久国产网址| 国产免费一级a男人的天堂| 三级国产精品欧美在线观看| 日韩人妻高清精品专区| 97人妻精品一区二区三区麻豆| 国产一区二区在线观看日韩| 色噜噜av男人的天堂激情| 老熟妇仑乱视频hdxx| 色在线成人网| 日日啪夜夜撸| 精品一区二区三区人妻视频| 在线天堂最新版资源| 国产伦一二天堂av在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产在线精品亚洲第一网站| 亚洲人成网站在线播| 国产成人aa在线观看| 亚洲成人久久爱视频| 在线免费观看不下载黄p国产| 国产精品伦人一区二区| av黄色大香蕉| 看片在线看免费视频| 国产精品乱码一区二三区的特点| 精品不卡国产一区二区三区| 一进一出抽搐gif免费好疼| 亚洲成人精品中文字幕电影| 久久99热这里只有精品18| 亚洲精品粉嫩美女一区| 精品免费久久久久久久清纯| 插阴视频在线观看视频| 欧美在线一区亚洲| 淫秽高清视频在线观看| 九九热线精品视视频播放| 狠狠狠狠99中文字幕| 久久久久九九精品影院| 久久人人爽人人片av| 亚洲国产精品久久男人天堂| 又爽又黄a免费视频| 中文在线观看免费www的网站| 成年女人毛片免费观看观看9| 成人鲁丝片一二三区免费| 亚洲图色成人| 久久久欧美国产精品| h日本视频在线播放| 两个人的视频大全免费| 午夜福利18| 我的女老师完整版在线观看| 激情 狠狠 欧美| 美女内射精品一级片tv| 国产精品嫩草影院av在线观看| 老女人水多毛片| 高清午夜精品一区二区三区 | 日产精品乱码卡一卡2卡三| 一进一出抽搐gif免费好疼| 欧美人与善性xxx| 麻豆一二三区av精品| 中文字幕久久专区| 国产精品99久久久久久久久| 人人妻人人澡人人爽人人夜夜 | 精品熟女少妇av免费看| 午夜福利18| 亚洲精品色激情综合| 亚洲av一区综合| 亚洲最大成人手机在线| 中出人妻视频一区二区| 免费观看精品视频网站| 精品一区二区三区人妻视频| 久久精品影院6| 国产片特级美女逼逼视频| 国内精品久久久久精免费| 全区人妻精品视频| 日本精品一区二区三区蜜桃| 国产淫片久久久久久久久| 久久久精品大字幕| 国产成人精品久久久久久| 亚洲综合色惰| 亚洲国产精品成人综合色| 亚洲国产精品久久男人天堂| 国产高清不卡午夜福利| 久久婷婷人人爽人人干人人爱| 99久久无色码亚洲精品果冻| 色av中文字幕| 精品福利观看| 99久久精品热视频| 在线观看午夜福利视频| 别揉我奶头 嗯啊视频| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久伊人网av| 一本久久中文字幕| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站 | 久久久久久久午夜电影| 美女 人体艺术 gogo| 久久亚洲国产成人精品v| 无遮挡黄片免费观看| 最近中文字幕高清免费大全6| 国产精品伦人一区二区| 少妇高潮的动态图| 草草在线视频免费看| 国产成人影院久久av| 老女人水多毛片| 身体一侧抽搐| 九色成人免费人妻av| 国产精品美女特级片免费视频播放器| 小蜜桃在线观看免费完整版高清| 蜜臀久久99精品久久宅男| 超碰av人人做人人爽久久| 国产精品伦人一区二区| 少妇高潮的动态图| 亚洲综合色惰| 欧美zozozo另类| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 日韩精品有码人妻一区| 一区二区三区免费毛片| 免费黄网站久久成人精品| 精品不卡国产一区二区三区| 免费高清视频大片| 热99re8久久精品国产| 亚洲成av人片在线播放无| 亚洲精品一区av在线观看| 日韩成人av中文字幕在线观看 | 国产欧美日韩精品亚洲av| 乱系列少妇在线播放| 亚洲人成网站在线观看播放| 免费看a级黄色片| 久久人人爽人人爽人人片va| 99视频精品全部免费 在线| 日本黄大片高清| 亚洲自拍偷在线| 亚洲精品色激情综合| 亚洲aⅴ乱码一区二区在线播放| 最近的中文字幕免费完整| 男女那种视频在线观看| 久久亚洲国产成人精品v| 午夜视频国产福利| 免费看光身美女| 少妇人妻一区二区三区视频| 久久久久久久久大av| 久久中文看片网| 久久精品夜夜夜夜夜久久蜜豆| 亚洲一级一片aⅴ在线观看| 国产高清视频在线观看网站| 观看美女的网站| 精品人妻熟女av久视频| 99热这里只有是精品50| 免费av不卡在线播放| 精品人妻视频免费看| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 亚洲av不卡在线观看| 欧美日本视频| 日本黄大片高清| 91麻豆精品激情在线观看国产| 女生性感内裤真人,穿戴方法视频| 性插视频无遮挡在线免费观看| 午夜精品国产一区二区电影 | 国产欧美日韩精品一区二区| 日产精品乱码卡一卡2卡三| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 1024手机看黄色片| 色av中文字幕| 欧美成人a在线观看| 日韩亚洲欧美综合| 级片在线观看| 成人性生交大片免费视频hd| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 精品人妻视频免费看| 亚洲av.av天堂| 嫩草影院新地址| 国产欧美日韩一区二区精品| 国产精品久久久久久久久免| 亚洲精华国产精华液的使用体验 | 国产免费男女视频| 悠悠久久av| 日本精品一区二区三区蜜桃| videossex国产| 国产精品久久电影中文字幕| 内射极品少妇av片p| 人人妻人人澡欧美一区二区| 激情 狠狠 欧美| 成人性生交大片免费视频hd| 久久久精品大字幕| 老熟妇仑乱视频hdxx| 一夜夜www| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 一夜夜www| 人妻丰满熟妇av一区二区三区| 自拍偷自拍亚洲精品老妇| 久久午夜福利片| a级一级毛片免费在线观看| avwww免费| 亚洲欧美中文字幕日韩二区| 久久久久免费精品人妻一区二区| 日日啪夜夜撸| 欧美xxxx黑人xx丫x性爽| 午夜免费男女啪啪视频观看 | 欧美高清性xxxxhd video| 亚洲第一区二区三区不卡| 在线观看美女被高潮喷水网站| av黄色大香蕉| 国产av在哪里看| 久久6这里有精品| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡| 欧美一区二区亚洲| 亚洲在线观看片| 午夜影院日韩av| 亚洲欧美清纯卡通| 中国美白少妇内射xxxbb| 国产精品,欧美在线| 99热这里只有是精品在线观看| 精品不卡国产一区二区三区| 18禁在线播放成人免费| 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 亚洲高清免费不卡视频| 97超碰精品成人国产| 久久精品国产清高在天天线| 最近的中文字幕免费完整| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出| 欧美极品一区二区三区四区| 国产高清视频在线播放一区| 91久久精品国产一区二区成人| av专区在线播放| 99久久精品热视频| 热99re8久久精品国产| 免费av不卡在线播放| 一级毛片电影观看 | 免费观看人在逋| 国产 一区 欧美 日韩| 亚洲成人精品中文字幕电影| 别揉我奶头~嗯~啊~动态视频| 日韩精品中文字幕看吧| 亚洲熟妇熟女久久| 99热6这里只有精品| 婷婷精品国产亚洲av| 久久久精品94久久精品| 国产色爽女视频免费观看| 美女被艹到高潮喷水动态| 中文亚洲av片在线观看爽| 亚洲精品国产av成人精品 | 久久人人爽人人片av| 内地一区二区视频在线| 日日摸夜夜添夜夜添av毛片| 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 给我免费播放毛片高清在线观看| 亚洲在线自拍视频| 亚洲中文字幕日韩| 色吧在线观看| 97热精品久久久久久| 精品人妻一区二区三区麻豆 | 日本精品一区二区三区蜜桃| 久久久色成人| 免费观看精品视频网站| 国产免费一级a男人的天堂| 日韩强制内射视频| 国产aⅴ精品一区二区三区波| 精品一区二区三区人妻视频| 超碰av人人做人人爽久久| 大又大粗又爽又黄少妇毛片口| 看非洲黑人一级黄片| 亚洲国产色片| 国产黄片美女视频| 国产精品精品国产色婷婷| 亚洲国产精品国产精品| 亚洲真实伦在线观看| 一本久久中文字幕| 久久久久国产网址| 国产一区二区亚洲精品在线观看| 别揉我奶头 嗯啊视频| av.在线天堂| 少妇的逼好多水| 午夜福利视频1000在线观看| 欧美另类亚洲清纯唯美| 美女被艹到高潮喷水动态| 国产精品福利在线免费观看| 日韩欧美一区二区三区在线观看| 日本爱情动作片www.在线观看 | 亚洲电影在线观看av| 日本 av在线| 搡老妇女老女人老熟妇| 美女高潮的动态| 国产成人a区在线观看| av福利片在线观看| 婷婷六月久久综合丁香| 亚洲av一区综合| 久久精品国产亚洲网站| 午夜日韩欧美国产| 中文亚洲av片在线观看爽| 内射极品少妇av片p| 久久国产乱子免费精品| 久久人人爽人人爽人人片va| 国产69精品久久久久777片| 免费黄网站久久成人精品| 小蜜桃在线观看免费完整版高清| 国产精品一区二区三区四区免费观看 | 免费高清视频大片| 美女内射精品一级片tv| 国产精品久久久久久久久免| 人妻制服诱惑在线中文字幕| 最好的美女福利视频网| 日韩三级伦理在线观看| 欧美成人一区二区免费高清观看| 能在线免费观看的黄片| 尾随美女入室| 国产高清视频在线播放一区| 日本色播在线视频| 亚洲五月天丁香| 男人狂女人下面高潮的视频| 看片在线看免费视频| 欧美日韩综合久久久久久| 国产精品一及| 丝袜美腿在线中文| 久久精品国产亚洲网站| 伦精品一区二区三区| 午夜福利18| 波多野结衣高清作品| 可以在线观看的亚洲视频| 国产欧美日韩精品一区二区| 99久久九九国产精品国产免费| 91久久精品电影网| 在线观看av片永久免费下载| 免费高清视频大片| 欧美xxxx性猛交bbbb| 久久久久国产网址| 亚洲综合色惰| 最后的刺客免费高清国语| 国产在视频线在精品| 成人毛片a级毛片在线播放| 国产精品精品国产色婷婷| 国产精品永久免费网站| 人妻夜夜爽99麻豆av| 我要搜黄色片| 欧美国产日韩亚洲一区| 三级国产精品欧美在线观看| 久久久精品欧美日韩精品| 精品日产1卡2卡| 全区人妻精品视频| 综合色av麻豆| 免费看日本二区| 国产老妇女一区| 色综合站精品国产| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 亚洲最大成人手机在线| 在线观看免费视频日本深夜| 国产高清有码在线观看视频| 国产在线男女| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 婷婷六月久久综合丁香| 成人亚洲精品av一区二区| 免费看a级黄色片| 最近的中文字幕免费完整| 1000部很黄的大片| 国产精华一区二区三区| .国产精品久久| 日日干狠狠操夜夜爽| 日韩在线高清观看一区二区三区| 国产一区二区亚洲精品在线观看| 中文字幕久久专区| 亚洲av免费在线观看| av在线老鸭窝| av免费在线看不卡| 一本精品99久久精品77| av免费在线看不卡| 天堂影院成人在线观看| 久久久成人免费电影| 亚洲精品乱码久久久v下载方式| 婷婷色综合大香蕉| 国产单亲对白刺激| 嫩草影视91久久| 好男人在线观看高清免费视频| 女同久久另类99精品国产91| 天天躁日日操中文字幕| 成人高潮视频无遮挡免费网站| 在线观看美女被高潮喷水网站| 少妇的逼水好多| 久久精品影院6| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 欧美最黄视频在线播放免费| 亚洲自偷自拍三级| 亚洲欧美日韩高清在线视频| 美女黄网站色视频| 最近最新中文字幕大全电影3| 色综合亚洲欧美另类图片| 久久草成人影院| 深夜a级毛片| 中出人妻视频一区二区| 最新中文字幕久久久久| 老师上课跳d突然被开到最大视频| 久久中文看片网| 国产在线精品亚洲第一网站| 亚洲国产欧洲综合997久久,| 男女那种视频在线观看| 亚洲自偷自拍三级| 天堂网av新在线| 五月玫瑰六月丁香| 一级毛片aaaaaa免费看小| 国产国拍精品亚洲av在线观看| 欧美潮喷喷水| 亚洲精品国产成人久久av| 国产精品不卡视频一区二区| 18禁在线播放成人免费| 国产精品99久久久久久久久| 久久6这里有精品| 大香蕉久久网| 亚洲五月天丁香| 欧美区成人在线视频| 嫩草影院入口| 国产高清不卡午夜福利| 午夜福利视频1000在线观看| 亚洲成人久久性| 韩国av在线不卡| 亚洲自拍偷在线| 精品午夜福利视频在线观看一区| 国内揄拍国产精品人妻在线| 欧美国产日韩亚洲一区| 日日撸夜夜添| 精品久久久久久久久久免费视频| 少妇的逼好多水| 国产精品野战在线观看| 乱人视频在线观看| 国产免费男女视频| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 性插视频无遮挡在线免费观看| av福利片在线观看| 午夜福利视频1000在线观看| 色哟哟哟哟哟哟| 欧美xxxx性猛交bbbb| 亚洲精品在线观看二区| 国产视频内射| 成人国产麻豆网| 最近最新中文字幕大全电影3| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 免费观看人在逋| 亚洲色图av天堂| 久久九九热精品免费| 日韩人妻高清精品专区| 欧美色欧美亚洲另类二区| 在线观看66精品国产| 亚洲精品影视一区二区三区av| 一级黄色大片毛片| 一本久久中文字幕| 大型黄色视频在线免费观看| 国产精品人妻久久久影院| 赤兔流量卡办理| 观看美女的网站| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩一区二区精品| 日韩欧美精品v在线| 欧美高清性xxxxhd video| 国产一区亚洲一区在线观看| 国产亚洲精品综合一区在线观看| 欧美激情在线99| 成年女人看的毛片在线观看| 在线看三级毛片| 一a级毛片在线观看| 国产精品一区二区免费欧美| 亚洲最大成人av| 白带黄色成豆腐渣| 春色校园在线视频观看| 在线观看av片永久免费下载| 欧美另类亚洲清纯唯美| 国产精品国产三级国产av玫瑰| 看黄色毛片网站| 两个人视频免费观看高清| 欧美日韩一区二区视频在线观看视频在线 | 日本爱情动作片www.在线观看 | 午夜福利视频1000在线观看| 97碰自拍视频| 国产91av在线免费观看| 日韩三级伦理在线观看| 国产伦在线观看视频一区| 欧美性感艳星| 亚洲精华国产精华液的使用体验 | 午夜免费男女啪啪视频观看 | 免费看日本二区| 成人三级黄色视频| 亚洲精品456在线播放app| 嫩草影视91久久| 1024手机看黄色片| 全区人妻精品视频| 一级毛片我不卡| 国语自产精品视频在线第100页| 亚洲七黄色美女视频| 久久久久免费精品人妻一区二区| 淫秽高清视频在线观看| 中文字幕久久专区| 美女黄网站色视频| 中文字幕免费在线视频6| 亚洲成人中文字幕在线播放| 免费在线观看成人毛片| 国产乱人视频| 人妻少妇偷人精品九色| 国产高潮美女av| 亚洲av免费高清在线观看| 九九在线视频观看精品| 一级毛片我不卡| 在线观看免费视频日本深夜| 欧美xxxx黑人xx丫x性爽| 三级经典国产精品| 69人妻影院| 国产伦在线观看视频一区| 国产精品爽爽va在线观看网站| 美女被艹到高潮喷水动态| 亚洲熟妇熟女久久| 啦啦啦韩国在线观看视频| 日韩三级伦理在线观看| 99热这里只有是精品50| 99九九线精品视频在线观看视频| 三级男女做爰猛烈吃奶摸视频| 国产成人aa在线观看| 黄片wwwwww| 日韩大尺度精品在线看网址| 午夜亚洲福利在线播放| 夜夜看夜夜爽夜夜摸| 成人鲁丝片一二三区免费| 国产精品,欧美在线| 在线播放无遮挡| 欧美xxxx黑人xx丫x性爽| 亚洲激情五月婷婷啪啪| 日本在线视频免费播放| 国产精品日韩av在线免费观看| 国产黄片美女视频| 国产美女午夜福利| 热99在线观看视频| 成人精品一区二区免费| 91久久精品国产一区二区成人| 国产女主播在线喷水免费视频网站 | 精品不卡国产一区二区三区| 国产美女午夜福利| 一a级毛片在线观看| 欧美一级a爱片免费观看看| 美女xxoo啪啪120秒动态图| 一级黄色大片毛片| 热99re8久久精品国产| 久久精品久久久久久噜噜老黄 | 天堂影院成人在线观看| 欧美性猛交黑人性爽| 波野结衣二区三区在线| 一进一出抽搐动态| 一区福利在线观看| 免费av观看视频| 欧美在线一区亚洲| 国产乱人视频| 内地一区二区视频在线| 99在线人妻在线中文字幕| 免费看光身美女| 麻豆成人午夜福利视频| 草草在线视频免费看| 天堂av国产一区二区熟女人妻| 女人十人毛片免费观看3o分钟| 色哟哟哟哟哟哟| 人妻久久中文字幕网| 亚洲自偷自拍三级| 国产免费男女视频| 国产片特级美女逼逼视频| 国产高清激情床上av| 久久久欧美国产精品| 最后的刺客免费高清国语| 午夜福利在线在线| 国产色爽女视频免费观看| 日本a在线网址| 日韩欧美 国产精品| 老司机福利观看| 成年女人永久免费观看视频| 三级经典国产精品| 色视频www国产| 国产精品一区二区三区四区久久| 精品少妇黑人巨大在线播放 | 男女视频在线观看网站免费| 国产三级中文精品| 亚洲第一电影网av| 一级毛片电影观看 | 一本精品99久久精品77| 在线观看美女被高潮喷水网站| 身体一侧抽搐| 女的被弄到高潮叫床怎么办| 国产精品日韩av在线免费观看| 三级男女做爰猛烈吃奶摸视频| 又爽又黄无遮挡网站| 97超视频在线观看视频| 国产免费男女视频| 欧美色视频一区免费| 亚洲美女视频黄频| 真人做人爱边吃奶动态| 六月丁香七月| 日韩精品青青久久久久久| 国产aⅴ精品一区二区三区波| 露出奶头的视频| 又黄又爽又免费观看的视频| 精品久久久久久久末码| 亚洲在线观看片| 国产高清不卡午夜福利|