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

    How to realize an ultrafast electron diffraction experiment with a terahertz pump: A theoretical study

    2022-05-16 07:10:48DanWang王丹XuanWang王瑄GuoqianLiao廖國前ZheZhang張喆andYutongLi李玉同
    Chinese Physics B 2022年5期
    關鍵詞:王丹李玉

    Dan Wang(王丹) Xuan Wang(王瑄) Guoqian Liao(廖國前)Zhe Zhang(張喆) and Yutong Li(李玉同)

    1Beijing National Laboratory for Condensed Matter Physics,Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

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

    3Songshan Lake Materials Laboratory,Dongguan 523808,China

    4Collaborative Innovation Center of IFSA,Shanghai Jiao Tong University,Shanghai 200240,China

    Keywords: terahertz pump,ultrafast electron diffraction,terahertz-electron interaction

    1. Introduction

    The rapid development of strong terahertz sources[1—10]opens up new opportunities for controlling and manipulating materials. Particularly, since the spectra of terahertz covers various low-energy excitations such as phonons in lattice,[11]rotations in molecular or hydration bonds in biological systems,[12]there is a great hope that with its electric field now exceeding 1 MV/cm,strong terahertz sources can directly perturb the nuclei degrees of freedom and even induce a structural phase transition,[13,14]free from the mediation of electronic excitations induced by infrared or visible light. Then,the subsequent structural dynamics can be monitored by a typical pump—probe experiment.Among various available probes that can extract the transient response of nuclei, diffraction methods such as ultrafast electron and x-ray diffraction techniques stand out, because they can follow the subtle changes of bond length or lattice constant together with the changes of structural symmetry through samples’ diffraction fingerprints with femtosecond temporal resolution.[15]So far, several terahertz-pump-x-ray-probe studies have been carried out which gain new insights of lattice dynamics.[14,16,17]However,the terahertz-pump-electron-probe studies on this topic are largely missing,[18]although terahertz sources have already been utilized to accelerate[19—28]or manipulate[19,22,23,29—38]electrons with many potential applications.

    The main obstacle for combining UED with terahertz sources lies in the fact that the electromagnetic field of the terahertz pump interacts strongly on the electron probe. Thus,if the terahertz pump and electron probe encounter each other before reaching the sample, the electron probe will be deflected preventing it from probing the structural dynamics afterwards. Figure 1 shows such deflection issue schematically.Suppose the incident angle of the electron probe in UED is initially set-up according to Bragg’s law

    to form diffraction patterns without terahertz-pumping. Here,dis the lattice constant,θis the diffraction angle,andλ0is the de Broglie wavelength of electrons. Then any pump-induced change indwill lead to a change inθcausing a relative shift in the position of the diffraction peak as Δd/d=Δθ/θaccording to Eq.(1). Meanwhile,the intensity of the diffraction peak may also alter either due to Δθ,[39]or due to other lattice dynamics such as lattice heating.[40]These fine structural changes can be resolved from the subtle differences of diffraction peaks between with and without pumping recorded by UED.However,if the electron probe encounters the terahertz pump before probing the sample as shown in Fig. 1(a), the pre-set incident angle (without pumping) is altered. Accordingly,the diffraction patterns are also altered or even disappear containing no information about lattice dynamics yet. As we will show in details later,such artifact varies dynamically with the probe’s time-delay according to the terahertz-pump,almost impossible to be corrected actively in a pump-probe experiment.And therefore,it must be carefully considered and eliminated when designing a terahertz-pump-electron-probe experiment. On the other hand,if such deflection happens after the sample is probed by the electrons as shown in Fig. 1(b), the diffraction patterns will still shift as a whole,yet the structural information has been recorded. Then, it could be recovered,in principle,even with improved temporal resolution since the terahertz can work as a streak camera.[22,30—33,35,41]

    Fig. 1. The electron probe encounters the terahertz pump (a) before the sample, which affects the probe’s pre-set incident angle, and (b) after the sample,which does not.

    In this article, we studied the terahertz-induced deflection effect on the electron probe of UED as shown in Fig.1(a)in details. We considered two colinear UED set-ups, namely the observe set-up and reversal set-up, in which the electron probe travels in the same direction as and in the opposite direction to the terahertz-pump,respectively. We found that for 100-keV electrons,no matter the sample is highly reflected or transmitted,significant deflection effect lasts so long that it is very difficult to carry out a pump—probe experiment with the observe set-up. To tackle this issue, we proposed two practical schemes by reshaping the terahertz waveform and by designing a noncolinear set-up. On the other hand, for highly reflected sample, a pump—probe experiment with the reversal set-up is largely not affected by such deflection effect, although the diffraction pattern will be shifted dynamically. It inspires us to use a metal-mesh as the sample holder to screen the terahertz pump for the reversal set-up. We also repeat all the above simulations for 3-MeV electrons. In general,the deflection effect for 3-MeV electrons is much smaller than for 100-keV electrons. With some minor adjustment, a pump—probe experiment can be realized in this case. We hope our study could arise people’s attention to the terahertz—electron interaction,and boost the multidisciplinary researches in both terahertz and UED societies.

    2. Methods

    Taking into account both the electric and magnetic effects of the terahertz field, the resultant force that acts on a single electron can be written as

    in whiche=1.60×10-19C is elementary charge,vrepresents velocity vector of the electron,EandBrepresent electric and magnetic terahertz field vectors,respectively.The acceleration of electron can be calculated as

    We considered two possible UED set-ups for terahertzpump—electron—probe experiment. In the obverse case as shown in Fig.2(a),the terahertz pump is focused on the sample by an off-axis parabolic mirror (OAP) with a horizontal hole through which the electron beam travels in the same direction as the terahertz pump.On the other hand,in the reverse case as shown in Fig.2(b),the terahertz pump and the electron probe travel collinearly in the opposite direction. The sample is assumed to be perpendicular to the electron beam.

    Fig.2. Experimental layout. (a)The obverse set-up, in which the electron moves in the same direction with focused terahertz pulse. (b)The reversal set-up,in which the electron moves in the opposite direction to the focused terahertz pulse.

    where “±” takes “+” in the obverse case and “-” in the reverse case. In Eq.(4),Erepresents the terahertz electric field

    Fig.3. The terahertz spectrum and waveform for a pump centered at fTHz =0.2 THz with fFWHM =π2/(16ln2)fTHz in panel(a)and in panel(b),respectively,and with fFWHM= fTHz/2 in panels(c)and(d),respectively.

    If we substitutet′=(±z/c)-t,then the terahertz waveform used in our simulation can be written as

    Since the terahertz beam size(several millimeters even at the focus)is much larger than the electron beam(several hundred microns), we use a single electron initially at (0, 0,-∞) att=-∞to characterize an electron pulse with finite transverse size. SinceEis set alongxaxis in our simulation(Ey=Ez=Bx=Bz=0),the electron accelerationa(x,y,z,t)can be derived according to Eqs.(2)and(3)as

    In our simulation the integration is achieved through step by step computing

    Table 1. Parameters used in the simulation.

    3. Results

    3.1. The obverse case

    We first consider the case when the terahertz pump completely transmits through the sample(T=100%). The corresponding evolution of deflection angleαand deflection distance Δxas a function of delay-timet0of a 100-keV electron probe is shown in Fig. 4(a). Depending on where the pump and probe encounter each other, their interaction can be divided into three regions as shown in Fig. 4(b). It can be seen that when the two beams overlap close to the sample (region I in Figs. 4(a) and 4(b)), electrons are strongly streaked. The overall time scale of this region is roughly the period of the terahertz pump. Also, in this region, only part of the terahertz waveform contributes onαand Δx, while the other part encounters electrons after the sample, which has no effect onαand Δx. Since the waveform of the terahertz beam starts with a quarter period of negative electric field followed by a half period of positive and another quarter period of negative field,αfirst increases, then decreases and increases again, accordingly. The streaking speed in the quasi-linear range (between the maximum and minimum ofα)is about-3.69×10-5rad/fs. Meanwhile,Δxalso changes several tenth micrometers. In region II, the two beams meet between the OAP and sample,so that the electron receives the interaction from the whole period of terahertz pulse. Since the terahertz beam is not focused transversely yet,the change ofαstarts to decrease quickly. And between 100 ps and 130 ps,the absolute change ofαis about 2.7157×10-4rad with an averaged streaking speed about 9.05×10-9rad/fs. On the other hand, Δxdoes not reduce likeαin this region, because the less focused terahertz beam the probe encounters, the longer distance it must travel to the sample. So, the reduced change inαis balanced by the increased travel distance. Finally, in region III,the two beams meet each other close to OAP,where only a portion of unfocused terahertz waveform encounters the probe. And therefore, a “wrinkle” appears in bothαand Δxjust like in region I,but with opposite variation trend,because the OAP cuts the wavefront with negative electric field first.

    Next, we consider the case when the terahertz pump is completely reflected by the sample(T=0%). And the corresponding evolution ofαand Δxis shown in Fig. 4(c). Since the electric field of the reflection beam points to the opposite direction of the original beam (see Fig. 4(d)), the deflection effect is partially cancelled. Also,because the reflection beam focuses tighter than the original beam when encounters the electrons,the deflection effect from the former is larger,leading to a positive value ofαand negative value of Δxin both regions I and II.In region III,electrons mainly encounter the original terahertz beam cut by OAP.So,a wrinkle appears inαand Δxcurves similar as the case of region III in Fig.4(a). In regions IV and V when the delay-time is long,electrons only meet the reflection beam, so the general variation trends of bothαand Δxcurves are just opposite to the ones in regions II and III in Fig.4(a), respectively. Comparing withT=100%case,the reflection of terahertz beam makes the amplitude ofαand Δxless but nonnegligible deflection effect lasts for longer time-scale.

    Fig.4. The deflection angle α and deflection distance Δx as a function of the delay-time t of 100-keV electrons in the obverse case when the sample transmittance T is(a)100%and(c)0%. The red solid line represents α,while the gray solid line represents Δx. (b)and(d): schematic diagrams for regions marked in panels(a)and(c).

    The continuous change inαand Δxmakes the above observe set-up unfavorable for a pump—probe experiment to study the transient structural dynamics in time domain. According to the two-beam dynamical diffraction theory under the assumption of no absorption, the intensity of Bragg peak can be described as[39,42]

    whereIrepresents the Bragg peak intensity,dsis the thickness of the sample,ξgis the extinction distance, andsis the excitation error that characterizes the deviation from the exact Bragg condition.[39]Since the electron probe always has certain transverse deviation,so the overall intensity of the Bragg peak should be considered as an integration through a range ofscovered by the deviation angle of the electron probe.Considering the initial energy spread of the photoelectrons is about 0.5 eV, the transverse deviation angle of the final 100-keV electron beam is about 0.002 rad, and this number decreases to about 0.001 rad for MeV electron beam because its longitudinal speed is doubled.

    To quantitatively determine how the intensity of Bragg peaks changes according toαis beyond the scope of this paper(see the appendix for more details),so here we use a simple criterion:if the change ofαis comparable to the deviation angle of the electron probe(in this case,0.002 rad),the intensity of the Bragg peaks are assumed to be attenuated to a degree that it couldn’t reflect the real lattice structural dynamics any longer. Judged by this criterion,region I in both Figs.4(a)and 4(c)cannot be used for a pump—probe experiment becauseαchanges drastically. However, region I is usually not counted in a pump—probe experiment, because the terahertz pump on the sample is not completed yet. Accordingly, the minimum delay-time of the electron probe should be at least one pulse width of the terahertz pump, roughly the same time scale of region I.And such streaking effect could be used for checking the alignment and synchronization between the electron and terahertz beam. Moving into region II, it can be seen that its early stage still could not be used for such experiments until about 50 ps, when the changing rate ofαbecomes much smaller. Even at this stage, the overall change ofαbetween 50 ps and 130 ps is still about 0.002 rad,comparable to the deviation angle of the electron beam. So,unless one can further divide this time range into several small sections and re-align the electron beam one section by one section, it is still difficult for a pump—probe experiment. All the other regions in Figs.4(a)and 4(c)are also hard to be used for a pump—probe experiment due to the existence of the wrinkles inαand Δxcurves. So, we may conclude that for the observe set-up, a pump—probe experiment could only be carried out whent0is longer than 150 ps forT=10%case,more than 10 periods of terahertz oscillation,and even longer forT=0%case when it reaches about 500 ps.

    To reduce the deflection from the terahertz pump in the observe set-up,here we propose two experimental schemes.

    The first one is to adjust the waveform of the terahertz field to make the contribution from the positive and negative fields cancel each other. For example,if we reduce the bandwidth of the terahertz beam and make the wave package contain several circles of oscillations, then the deflection from it will be significantly reduced. Figures 5(a)and 5(b) show the evolution ofαand Δxcaused by a terahertz beam that has roughly two circles of oscillations as shown in Fig. 3(d). It can be seen that for bothT=100%andT=0%cases,region II can be used for pump—probe experiment in which the amplitude ofαis much smaller than the deviation angle of the 100-keV electron beam. Another approach to obtain multi-cycle terahertz waveforms is to use higher frequency terahertz pulse.As shown in Fig.5(c),the higher frequency the terahertz pump has,the less it deflects the probe. For example,with the same electric field strength,the magnitude ofαis reduced by a factor of 592 forfTHz=10 THz compared tofTHz=0.2 THz. Another possible method to achieve the same goal is to change the carrier-envelope offset phase(CEP)of the terahertz pump to balance the positive and negative field strengths. For example,if the carrier-envelope offset phase of the terahertz pump shown in Fig.3(b)is shifted aboutπ/4,then the magnitude of positive and negative fields will be close to each other leading to a cancelation inαcurve. Now terahertz sources with different carrier-envelope offset phases and frequencies are available through different generation methods.[43—45]

    Fig.5. The deflection angle α and deflection distance Δx as a function of delay-time t0 for 100-keV electrons in obverse case when the waveform of terahertz field shown in Fig.3(d)is used with(a)T =100%and(b)T =0%. (c)The maximum of α and Δx as a function of fTHz with the waveform of terahertz field shown in Fig.3(b)and T =100%.

    Fig.6. Noncolinear scheme in the obverse case: (a)the experiment layout with(b)a highlight of the sample location. The corresponding deflection angle α and deflection distance Δx as a function of delay-time t0 for 100-keV electrons in obverse case with(c)T =100%and(d)T =0%. (e)Schematic diagram for a tilted sample. (f)The corresponding α and Δx as a function of t0 for 100-keV electrons in the tilting sample scheme with T =0%and γ =π/4.

    To further increase the time window for pump—probe experiment,we come up with a noncolinear scheme(Fig.6(a)),in which the OAP is tilted towards the sample so that the pump and probe travels noncolinearly. In this respect,the spacial overlap of the two beams is only achieved very close to the sample location (mm scale), and therefore, the deflection from the pump attenuated very quickly after about 10 ps. Figures 6(c) and 6(d) show the evolution ofαand Δxin such noncolinear scheme forT=100%andT=0%,respectively.In the simulation,the tilting angle of terahertz beam according to electron beam (γ) is set to beπ/4, while other parameters are the same as used in Fig. 4. It can be seen thatαreturns to zero about 10 ps and 20 ps forT=100%andT=0%,respectively. Also,the maximal amplitude of Δxis only several micrometers. Comparing Fig.6 with Fig.4,the advantage of noncollinear scheme on collinear one is clearly seen.

    One potential issue for such scheme is that the projection of the electric field on thezdirection will accelerate or decelerate electrons dynamically so that both the central energy and the energy spread of the probe will change as a function oft0.However, if we only use the time domain whenαis close to zero, then, such effect is negligible. Another issue for such scheme is that it will introduce the velocity and the geometrical mismatch[46]that worsen the overall temporal resolution as shown in Fig.6(b),because now the sample is not pumped and probed simultaneously alongxaxis. The overall temporal resolution could be estimated by[46]

    whereτpumpandτproberepresent the longitudinal durations of the terahertz pulse and electron beam,respectively.τjitterrepresents the timing jitter between the pump and the probe. Andτmisis the temporal smear due to the velocity and the geometric mismatch[46]as we just mentioned,which can be estimated by

    whereleis electron beam size,respectively.Assumingτpump=1 ps,τprobe=200 fs,τjitter=0 ps,ds=100 nm,le=100 μm,and even usingγ=π/2,we getτtotal=3.1 ps,not significantly lengthened if comparing with the time scale of region I. But asγincreases, the pumping area also increases as 1/cos(γ)which decreases the pumping fluence,another factor one must weight when designing the experiment.

    We also considered the possibility to tilt the sample(see Fig. 6(e)) rather than the pump. The difference between this scheme and the noncolinear scheme is that the former only reduces the spacial overlapping between the reflected pump beam and the probe. Since the original pump beam also deflects the electron probe significantly,this scheme does not extend the usable time window too much. Figure 6(f)shows the simulation result for this case withγ=π/4 andT=100%. It can be seen that in the first 20 ps(roughly region I in Fig.4),the deflection effect does reduce significantly if compared with the noncolinear scheme. However, it lasts for a longer time scale similar as the region II in Fig. 4(a) due to the contribution from the original pump beam. Also,the overall temporal resolution is reduced by the velocity and the geometric mismatch that can be estimated by

    As a result, it is not an effective approach for 100-keV electrons, But, as we will show later, it is very helpful for UED with 3-MeV electrons.

    3.2. The reversal case

    In the reversal case (Fig. 2(b)), the electron and the terahertz pulses travel in the opposite direction (γ=π) so that they only meet each other once. Figure 7(a)shows the evolution curves ofαand ΔxforT=100%,and figure 7(c)shows the two curves forT=0% and 1%, respectively. Here, we ignored the tinyτmis(less than 1 fs) in the simulation. ForT=100%, the interaction between pump and probe is similar to the situation described in Fig. 4(a) for the initial two regions,and both two curves also look very similar in the two cases. On the other hand,forT=0%,bothαand Δxare not affected by the terahertz beam because the pump and probe encounter each other after the sample and it does not affect the pre-set diffraction angle and probed position. ForT=1%,the deflection is just 1% ofT=100% case confirming that only the transmitted beam affectsα.

    Fig.7. The deflection angle α and deflection distance Δx as a function of delay-time t0 for 100-keV electrons in reversal case with(a)T =100%and(c)T =0%(horizontal lines)and 1%(solid curves). (b)and(d)Schematic diagrams for regions marked in panels(a)and(c). (e)Schematic diagram that shows how to screen the terahertz pump by a mesh in the reversal case.

    Compared with the obverse set-up, the reversal set-up withT=0% seems to be more situatable for a pump—probe experiment. To enhance the reflection of terahertz beam in this case, we proposed to use metal grid to hold the sample as shown in Fig. 7(e). So long as the wavelength of the terahertz beam is much larger than the pitch size of the mesh,terahertz wave will be screened out to achieveT=0%. In fact,these metal grids are widely used for transmission electron microscopy.The drawback for the reversal set-up is that OAP can block the diffracted electron beam. To solve this problem,one may apply the noncolinear scheme (see Fig. 6(a)) but with a negative titling angle.

    3.3. Using higher energy electrons

    According to Eqs. (6) and (7), bothvx(t) andqx(t) decreases asvincreases due to the relativistic effect, thus bothαand Δxdecreases around one order in magnitude for 3-MeV electrons if compared with 100-keV ones as shown in Fig. 8(also see Figs.4 and 7 for comparison). ForT=100%of the obverse set-up,the evolution ofαcurve mimics the waveform of terahertz beam(see Fig.3(b)),which is different from all the curves for 100-keV electrons. This is becausevnow is very close to the speed of light so that the electron interacts with terahertz mostly in pace. On the contrary, for lower energy(100 keV for example) electrons, terahertz beam surpasses electrons from behind so that a large portion of its waveform deflects electrons.At a first glance,this seems like a disadvantage because it prevents the contributions from the negative and positive field in one terahertz circle to cancel each other.However, such disadvantage is compensated by the cancelation between the forces from the electric and magnetic field,as the latter one increases withvto a comparable magnitude to the former (see Eq. (6)). So, comparing theαcurve in Fig. 8(a) with Fig. 4(a), its magnitude is largely reduced by a factor of 16. Also, as shown in Fig. 8(d), the maximum ofαdecreases when the energy of electron increases. Since the velocity of MeV electrons are close to the speed of light,αvery quickly decreases to zero around one period of terahertz,because longer than this time-delay,the terahertz beam is always ahead of electrons and they do not encounter each other any longer, another advantage for using MeV electrons in a terahertz—pump—electron—probe experiment. ForT=0%in obverse set-up, considerable amount ofα(in the order of 0.01 rad) persists to several tenth picoseconds due to the reflection of the terahertz pump(in this case,the forces from the electric and magnetic field add up) which impedes a pump—probe experiment. One possible method to solve this issue is to tilt the sample instead of terahertz pump to reduce the deflection from the reflected terahertz pump only.But it will sacrifice certain temporal resolution due to the geometrical mismatch(see Eq.(12)).

    In the reversal case,bothαand Δxcurves resemble those for 100-keV electrons(see Figs.7(a)and 7(c))but with almost one order smaller in magnitude,which make 3-MeV electrons better than 100-keV electrons for such a pump—probe experiment. But it still faces the blocking issue of diffracted electrons by the OAP,which could be improved by using the noncolinear pumping scheme(see Fig.7(a)).

    So far, we only consider the deflection effect of the terahertz pump on the intensity of Bragg peaks. In many UED experiments, the transient change in lattice constants is also monitored from the position shift of the Bragg peaks as a function of time,characterized by the relative peak position change Δθ/θ. Then,it is more proper to useα/θto characterize the deflection effect on Bragg peak positions. Figure 8(d) shows how the maximum ofα/θchanges as a function of electron energy in the obverse case. And we can see that the deflection effect on Bragg peak positions starts to decrease rapidly when the energy of electrons enters MeV range,making MeV UED a better choice for such study,although the much smaller Bragg angles in MeV UEDs may bring other challenges in resolving Bragg peaks.

    Fig.8. The deflection angle α and deflection distance Δx as a function of delay-time t0 for 3-MeV electrons in obverse case with(a)T =100%and(b)T =0%and in reversal case with(c)T =100%and 1%. (d)The maximum of α and deflection angle divided by diffraction angle α/θ for electrons with different energies for sample with d=5 °A and T =100%in the obversed case.

    3.4. The deflection from terahertz-induced non-propagating field

    Up to now, we only considered a very ideal case, where the sample is assumed to be perfectly flat and is perpendicular to the incident electron beam, so that the influence of the sample on terahertz beam is simplified by the reflection and negligible absorption. But in the real situation, the dielectric constant, the geometry, and the local structure of the sample and the sample mount can strongly attenuate the localEandBfields. If such near-field is evanescent (immediately disappears after the terahertz pump and roughly the same timescale as region I), we will show that its effect on the UED experiment is not so important. First,we argue that the effective interaction length alongzaxis between such evanescent near-field and electrons is much less than the wavelength of the pump beam. This argument comes from several superresolution techniques such as near-field scanning optical microscopy, where a change in the local field in the order of one tenth of the pump’s wavelength can be resolved. On the other hand, the effective interaction length between electrons and terahertz beam is around the Rayleigh length (in the order of the terahertz wavelength under the condition that the F-number of the OAP is around 1),within which the terahertz beam is tightly focused. As a result, unless the near-field is strongly enhanced such as in a terahertz resonator (then, it is not evanescent), its deflection effect is negligible if compared with that from the original and reflected beam. We think this is probably the case for a typical TEM sample,where the length scale of its roughness or the wrinkles on it are both much smaller than the terahertz wavelength, causing random scattering of light only. It is also true for our proposal where we use a metal grid as the sample holder as shown in Fig.7(e),where the terahertz wave is completely screened(reflected).

    However, if such terahertz-induced field is not evanescent, such as the dipole field due to the plasma,[47]the resonant field in a split-ring[48]and so on,the amplitude of such local field can be very large and its effective interaction length on electrons can be also significantly extended. As a result,their deflection effect dominates. In fact, ultrafast electron has already been used as a probe to map such transient field in several cases.[48—51]But such sample conditions should be avoided when designing a UED experiment.

    Finally,we want to give only a qualitative analysis on how such enhanced local field affects 3-MeV UED,for the cancellation betweenEandBis only important for such high energy electrons. We assume that the field does not propagate in pace with the electrons. Otherwise, the deflection effect will enhance a factor of one hundred without the cancellation betweenEandB. Now, the interaction between terahertz pulse and electron is very similar as for the low energy electron cases but the deflection force reduces a factor of 5.74 due to the relativistic mass of 3-MeV electrons. Given the field strength same as the terahertz pump and an effective interaction length of the Rayleigh length,from Figs.4(a)and 4(c),the maximal change ofαcan be estimated to be about 0.001 rad. This is also consistent with Fig.8(b),where the maximal change ofαis doubled to 0.002 rad,becauseEandBare added in the reflected terahertz beam. Based on these estimations,if the field can be localized in the length scale of at least one order less than the Rayleigh length,it will not affect the diffraction patterns significantly. However,if the strength of the local field is enhanced,then it must be confined further to reduce the effective interaction length on electrons.

    4. Conclusions

    In this article, we theoretically studied the terahertzinduced deflection effect on the electron probe of UED and how it affects the incident angle and the probe location on the sample. In general, when the electron energy reaches MeV range,the terahertz-induced deflection becomes small enough to realize a pump—probe experiment. However, for electron probe with its energy around 100 keV or less,such experiment is much harder to design,although this energy range is widely used in electron microscopy, and is more suitable for some fragile samples like 2D materials. Particularly,in the observe set-up, both original pump beam and its reflection from the sample deflects the electron probe significantly, while in the reverse set-up, only the original beam contributes. We introduced several schemes to tackle this issue. This first scheme is to reshape the waveform of the terahertz pump, so that the deflection from its positive and negative portion may cancel each other. Since the available tools to manipulate the optical parameters of a terahertz pulse such as its frequency,bandwidth and CEP are still sparce,the applicability of this scheme mostly depends on how the terahertz pump is generated. The other scheme is to use a noncolinear pump—probe set-up to reduce the spacial overlap between the pump and probe. As a result,the time-window suitable for a pump—probe experiment extends significantly. However,due to the velocity and the geometric mismatch,the overall temporal resolution and the energy stability of the electron probe is worsened, and the field amplitude of the terahertz pump is reduced. For the reverse set-up, we proposed to use a metal grid as the sample holder to screen the terahertz pump in order to eliminate the transmitted pump. We hope our study could arise people’s attention to the terahertz—electron interaction, and blend the merits of terahertz and UED techniques.

    Acknowledgements

    Project supported by the National Natural Science Foundation of China (Grant Nos. 11774409, 11827807, and 92050106). This work is supported by the Synergic Extreme Condition User Facility.

    Appendix A

    In this appendix,we would like to show howαaffects the intensity of Bragg peaks.The reason we use absolute angle instead of relative angle is that we mainly consider the intensity change of the Bragg peak due to the change of,which depends on the product ofξgands. Let us assume that the UED is initially perfectly aligned to obtain maximal diffraction intensity withξgs=0. Then, according to Eq.(9), the diffraction intensity will drop significantly whenξgs~1 due toα.

    Here,ξgcan be expressed as[42]and for the sameg(the same plane index),F(θ) ∝f(θ),where

    which is plotted as a function of electron energy as shown in Fig. A1. It can be seen thatξgincreases about a factor of 1.8 when the electron energy increases from 100 keV to MeV.This means that the diffraction intensity for MeV UED drops larger than 100-keV UED,given the sames. However,it does not depend onθ.Then, from Eqs. (A6) and (A8), we may conclude thatξgsonly depends onα.

    Fig.A1. The ξg as a function of electron energy Ee.

    Fig. A2. Schematic diagrams for the relation between s and α. kI is the wave vector of the incidence electron, and kD is the wave vector of the diffracted electron. k′I and k′D represent the wave vectors deflected by α due to the terahertz pump.

    In this article,we used a simple criterion:if the change ofαis comparable to the deviation angle of the electron probe,the intensity of the Bragg peaks is assumed to be attenuated to a degree that it could not reflect the real lattice structural dynamics any longer. The deviation angle for different ultrafast electron sources depends on many factors, which is hard to quantify. To make a fair comparison, here, we just simply assume that the initial energy spread of the photoelectrons(0.5 eV)is fixed. Since the velocity of MeV electrons is twice of 100-keV ones, the resultant deviation angle for MeV electrons is only half of 100 keV which balanced the factor of 1.8 inξg. So, our compassion between UEDs with different energy electrons based on this simple criterion is justified.

    猜你喜歡
    王丹李玉
    王丹美術作品欣賞
    大眾文藝(2023年18期)2023-10-13 02:07:16
    從《綴白裘》看李玉《一笠庵四種曲》的場上變遷
    戲曲研究(2023年1期)2023-06-27 06:55:48
    兩個“倍半”巧解題
    虛構的野心:“私生飯”追星自陷迷奸危局
    王丹作品
    藝術家(2019年8期)2019-12-16 22:35:07
    On the Category and Conceptual Process of Metaphor
    OntheEssentialIngredientsofModernCulturalSyllabus
    山東青年(2015年2期)2015-12-24 13:01:26
    Based on a literature review of rural finance development
    論李玉電影作品的“跳切”剪輯對其藝術風格的影響
    電影新作(2014年4期)2014-02-27 09:13:17
    捷徑是條不歸路
    做人與處世(2013年7期)2013-06-27 08:17:18
    一级作爱视频免费观看| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器 | 婷婷精品国产亚洲av在线| 神马国产精品三级电影在线观看 | 国产精品乱码一区二三区的特点| 美女国产高潮福利片在线看| 国产成人av激情在线播放| 50天的宝宝边吃奶边哭怎么回事| 精品卡一卡二卡四卡免费| 精品一区二区三区视频在线观看免费| 精品久久久久久,| 欧美+亚洲+日韩+国产| 国产三级在线视频| 亚洲,欧美精品.| 国产黄a三级三级三级人| 欧美性猛交黑人性爽| 夜夜爽天天搞| 人人妻人人澡人人看| 中文字幕精品亚洲无线码一区 | 91大片在线观看| 色尼玛亚洲综合影院| 高潮久久久久久久久久久不卡| 成人欧美大片| 激情在线观看视频在线高清| 国产爱豆传媒在线观看 | 在线观看免费视频日本深夜| 此物有八面人人有两片| 欧美午夜高清在线| 视频区欧美日本亚洲| 性色av乱码一区二区三区2| 午夜福利在线观看吧| 后天国语完整版免费观看| 在线观看日韩欧美| 一本一本综合久久| 精品午夜福利视频在线观看一区| 亚洲在线自拍视频| 亚洲自偷自拍图片 自拍| 亚洲人成网站在线播放欧美日韩| 99国产极品粉嫩在线观看| 亚洲国产精品sss在线观看| 国产精品久久久人人做人人爽| 日韩欧美三级三区| 91av网站免费观看| 欧美人与性动交α欧美精品济南到| 久99久视频精品免费| 成年女人毛片免费观看观看9| 色av中文字幕| 我的亚洲天堂| 成人国产综合亚洲| 最好的美女福利视频网| 欧美黑人欧美精品刺激| 人成视频在线观看免费观看| 中文字幕人妻熟女乱码| 亚洲成国产人片在线观看| 国产精品永久免费网站| 两性午夜刺激爽爽歪歪视频在线观看 | 两性夫妻黄色片| 99国产精品一区二区蜜桃av| 亚洲色图 男人天堂 中文字幕| 97超级碰碰碰精品色视频在线观看| 国产久久久一区二区三区| 怎么达到女性高潮| 白带黄色成豆腐渣| 亚洲欧美精品综合一区二区三区| 欧美又色又爽又黄视频| 在线免费观看的www视频| 欧美性猛交╳xxx乱大交人| 精品久久久久久久末码| 首页视频小说图片口味搜索| 国产在线精品亚洲第一网站| 婷婷精品国产亚洲av| 国内少妇人妻偷人精品xxx网站 | 窝窝影院91人妻| 欧美性长视频在线观看| 午夜福利高清视频| 午夜a级毛片| 91麻豆精品激情在线观看国产| 黄色视频,在线免费观看| www.www免费av| 99久久精品国产亚洲精品| 国产片内射在线| 亚洲av成人一区二区三| 午夜两性在线视频| 久久狼人影院| 免费女性裸体啪啪无遮挡网站| 好男人在线观看高清免费视频 | 老汉色∧v一级毛片| 2021天堂中文幕一二区在线观 | 日本一区二区免费在线视频| 色尼玛亚洲综合影院| 一级毛片精品| 侵犯人妻中文字幕一二三四区| 99久久久亚洲精品蜜臀av| 亚洲自拍偷在线| 黄片小视频在线播放| 久久久久久国产a免费观看| 黄色视频,在线免费观看| 日本精品一区二区三区蜜桃| 欧美国产精品va在线观看不卡| 女人被狂操c到高潮| 在线播放国产精品三级| 成人国产综合亚洲| 精品少妇一区二区三区视频日本电影| 亚洲精品在线观看二区| 亚洲精品久久成人aⅴ小说| 精品久久久久久久久久免费视频| 欧美不卡视频在线免费观看 | 亚洲天堂国产精品一区在线| 超碰成人久久| 久久欧美精品欧美久久欧美| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久av网站| 亚洲七黄色美女视频| 琪琪午夜伦伦电影理论片6080| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久久久电影| 丁香六月欧美| 妹子高潮喷水视频| 亚洲激情在线av| 19禁男女啪啪无遮挡网站| 一个人观看的视频www高清免费观看 | 婷婷精品国产亚洲av| 自线自在国产av| 午夜免费观看网址| 午夜免费激情av| 制服丝袜大香蕉在线| 一边摸一边抽搐一进一小说| 婷婷丁香在线五月| 欧美乱色亚洲激情| 精品电影一区二区在线| 久久国产精品男人的天堂亚洲| 国产精品日韩av在线免费观看| 亚洲精华国产精华精| 99热6这里只有精品| 亚洲国产中文字幕在线视频| 国产伦一二天堂av在线观看| 国产精品综合久久久久久久免费| 久久久久久久久久黄片| 18禁裸乳无遮挡免费网站照片 | 一级片免费观看大全| aaaaa片日本免费| 夜夜躁狠狠躁天天躁| 最近最新中文字幕大全电影3 | 国产精品电影一区二区三区| 99国产精品99久久久久| 国产午夜福利久久久久久| 日韩欧美一区视频在线观看| 69av精品久久久久久| 99久久精品国产亚洲精品| 久久国产精品人妻蜜桃| 亚洲国产日韩欧美精品在线观看 | 午夜视频精品福利| 好男人电影高清在线观看| 国产精品久久久人人做人人爽| 亚洲av第一区精品v没综合| aaaaa片日本免费| 中文字幕久久专区| 好看av亚洲va欧美ⅴa在| 99久久久亚洲精品蜜臀av| 午夜免费观看网址| 久久天躁狠狠躁夜夜2o2o| www国产在线视频色| 国产激情久久老熟女| 久久久久久久久中文| 久久精品亚洲精品国产色婷小说| 此物有八面人人有两片| 十八禁网站免费在线| 日韩欧美免费精品| 国产极品粉嫩免费观看在线| 精品卡一卡二卡四卡免费| 国产成人欧美| 免费在线观看完整版高清| 色播亚洲综合网| 久久久久国产一级毛片高清牌| 制服人妻中文乱码| 日韩大尺度精品在线看网址| 成年版毛片免费区| 久久久国产欧美日韩av| 久久精品影院6| 中亚洲国语对白在线视频| 老司机午夜福利在线观看视频| 观看免费一级毛片| 国产精品99久久99久久久不卡| 精品久久久久久,| www日本黄色视频网| 天堂动漫精品| 日本 欧美在线| www.www免费av| 俺也久久电影网| 特大巨黑吊av在线直播 | 国产在线精品亚洲第一网站| 动漫黄色视频在线观看| 久久天堂一区二区三区四区| 两性夫妻黄色片| 久久香蕉国产精品| 在线观看日韩欧美| 亚洲国产日韩欧美精品在线观看 | 69av精品久久久久久| 性色av乱码一区二区三区2| 99热6这里只有精品| 一级毛片精品| 亚洲精品一区av在线观看| 精品国产一区二区三区四区第35| av在线天堂中文字幕| 叶爱在线成人免费视频播放| 国产真实乱freesex| 高清在线国产一区| 99热这里只有精品一区 | 黄网站色视频无遮挡免费观看| 午夜福利欧美成人| cao死你这个sao货| 午夜久久久久精精品| 国产精品美女特级片免费视频播放器 | 久久久久久九九精品二区国产 | 99久久无色码亚洲精品果冻| 国产av一区在线观看免费| 人妻丰满熟妇av一区二区三区| 色综合欧美亚洲国产小说| 男女下面进入的视频免费午夜 | 老司机午夜十八禁免费视频| 一二三四在线观看免费中文在| 亚洲精品粉嫩美女一区| 久久这里只有精品19| 窝窝影院91人妻| 国产精品久久久久久亚洲av鲁大| 亚洲精品中文字幕一二三四区| 狂野欧美激情性xxxx| 日韩国内少妇激情av| 日韩成人在线观看一区二区三区| 色播亚洲综合网| 欧美国产日韩亚洲一区| 一本综合久久免费| 亚洲激情在线av| 国产色视频综合| 91老司机精品| 麻豆久久精品国产亚洲av| АⅤ资源中文在线天堂| bbb黄色大片| 成人精品一区二区免费| 国产又爽黄色视频| 久久性视频一级片| 国产成人av激情在线播放| 国产高清视频在线播放一区| 91av网站免费观看| 国产精品久久电影中文字幕| 亚洲电影在线观看av| 亚洲第一青青草原| 亚洲av熟女| 中国美女看黄片| 人人妻人人澡人人看| 亚洲 欧美一区二区三区| 手机成人av网站| 999久久久国产精品视频| 久久香蕉国产精品| 99精品久久久久人妻精品| 国产真实乱freesex| 老司机午夜福利在线观看视频| 一本精品99久久精品77| 99精品在免费线老司机午夜| 欧美久久黑人一区二区| 50天的宝宝边吃奶边哭怎么回事| 久久久久免费精品人妻一区二区 | 麻豆一二三区av精品| 久久精品国产综合久久久| 欧美成人午夜精品| 国产在线观看jvid| 少妇被粗大的猛进出69影院| 久久久久久免费高清国产稀缺| 欧美亚洲日本最大视频资源| 色播在线永久视频| 波多野结衣巨乳人妻| 欧美成人午夜精品| 久久婷婷成人综合色麻豆| 97碰自拍视频| 久久婷婷人人爽人人干人人爱| 午夜免费观看网址| 亚洲三区欧美一区| 日韩欧美国产一区二区入口| 中文字幕人妻丝袜一区二区| 亚洲av日韩精品久久久久久密| 亚洲黑人精品在线| 国产成人精品无人区| 国产欧美日韩一区二区精品| 日本a在线网址| 一本久久中文字幕| 日韩视频一区二区在线观看| 两个人视频免费观看高清| www.999成人在线观看| АⅤ资源中文在线天堂| 免费高清视频大片| 一本久久中文字幕| 正在播放国产对白刺激| 国产成人精品久久二区二区免费| 97人妻精品一区二区三区麻豆 | 人人妻人人澡人人看| 最新在线观看一区二区三区| 国产v大片淫在线免费观看| 中文字幕av电影在线播放| 国产伦一二天堂av在线观看| 精品电影一区二区在线| www日本黄色视频网| 男女做爰动态图高潮gif福利片| 啦啦啦 在线观看视频| 黄色视频,在线免费观看| 99国产精品一区二区蜜桃av| 女警被强在线播放| 婷婷丁香在线五月| 我的亚洲天堂| 欧美日韩瑟瑟在线播放| 免费无遮挡裸体视频| 亚洲国产欧美日韩在线播放| 精品卡一卡二卡四卡免费| www.www免费av| 99在线视频只有这里精品首页| 给我免费播放毛片高清在线观看| 国产极品粉嫩免费观看在线| 最新美女视频免费是黄的| 精品高清国产在线一区| 国产熟女xx| 女同久久另类99精品国产91| 久久国产乱子伦精品免费另类| 最好的美女福利视频网| 国产野战对白在线观看| 变态另类丝袜制服| 日韩 欧美 亚洲 中文字幕| 最近最新免费中文字幕在线| 久久 成人 亚洲| 欧美日韩瑟瑟在线播放| 老汉色av国产亚洲站长工具| 久久久久国内视频| 黑人操中国人逼视频| 国产在线观看jvid| 午夜精品在线福利| 欧美日韩亚洲国产一区二区在线观看| 好看av亚洲va欧美ⅴa在| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 成人亚洲精品av一区二区| 高清在线国产一区| 91字幕亚洲| 久久精品影院6| 少妇的丰满在线观看| 脱女人内裤的视频| 午夜影院日韩av| 18禁美女被吸乳视频| 久久久国产成人免费| 欧美乱色亚洲激情| 国产精品一区二区三区四区久久 | 国产亚洲精品久久久久久毛片| 国产亚洲欧美在线一区二区| x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 日韩欧美一区视频在线观看| 无人区码免费观看不卡| 成人永久免费在线观看视频| 一区福利在线观看| 中国美女看黄片| 精品熟女少妇八av免费久了| 亚洲中文日韩欧美视频| 国产免费av片在线观看野外av| 免费看十八禁软件| 黑人操中国人逼视频| 黄色成人免费大全| 啦啦啦 在线观看视频| 性欧美人与动物交配| 亚洲欧美精品综合一区二区三区| 精品不卡国产一区二区三区| 久久中文字幕人妻熟女| 特大巨黑吊av在线直播 | 久久久久久九九精品二区国产 | 国产真人三级小视频在线观看| 一进一出好大好爽视频| 色综合站精品国产| 丰满人妻熟妇乱又伦精品不卡| 午夜久久久久精精品| 日本在线视频免费播放| 好男人电影高清在线观看| 国产精品久久久久久人妻精品电影| 国产激情偷乱视频一区二区| 国产蜜桃级精品一区二区三区| 欧美色欧美亚洲另类二区| 欧美成人午夜精品| 国产高清有码在线观看视频 | 国产午夜精品久久久久久| 精品免费久久久久久久清纯| 757午夜福利合集在线观看| 99在线人妻在线中文字幕| 精品国产乱子伦一区二区三区| 99re在线观看精品视频| √禁漫天堂资源中文www| 韩国精品一区二区三区| 两人在一起打扑克的视频| 一区二区三区激情视频| 在线观看免费日韩欧美大片| 欧美精品亚洲一区二区| 免费看十八禁软件| 少妇熟女aⅴ在线视频| 变态另类丝袜制服| 狂野欧美激情性xxxx| 国内少妇人妻偷人精品xxx网站 | 午夜福利在线在线| 18美女黄网站色大片免费观看| 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 亚洲在线自拍视频| 免费电影在线观看免费观看| 精品午夜福利视频在线观看一区| 日韩欧美一区二区三区在线观看| 中出人妻视频一区二区| 日韩免费av在线播放| av有码第一页| 九色国产91popny在线| 美女免费视频网站| 99久久国产精品久久久| 色老头精品视频在线观看| 国产伦人伦偷精品视频| 成人欧美大片| 宅男免费午夜| 日日夜夜操网爽| 亚洲片人在线观看| e午夜精品久久久久久久| 搡老岳熟女国产| 久久香蕉激情| 国产亚洲精品久久久久久毛片| 国产三级黄色录像| 9191精品国产免费久久| 99热只有精品国产| 精品欧美一区二区三区在线| 在线观看www视频免费| 亚洲精品国产精品久久久不卡| www.熟女人妻精品国产| 97人妻精品一区二区三区麻豆 | 午夜久久久在线观看| 夜夜看夜夜爽夜夜摸| 精品欧美国产一区二区三| 欧美激情久久久久久爽电影| 久久久久久大精品| 怎么达到女性高潮| 俄罗斯特黄特色一大片| 可以在线观看毛片的网站| 老司机靠b影院| 在线观看日韩欧美| 国产成人欧美| 十八禁人妻一区二区| 国产高清视频在线播放一区| 在线观看免费日韩欧美大片| 777久久人妻少妇嫩草av网站| 又大又爽又粗| 亚洲色图av天堂| 亚洲av第一区精品v没综合| 欧美黄色淫秽网站| 丁香欧美五月| 国产区一区二久久| 此物有八面人人有两片| 欧美黑人巨大hd| 搡老熟女国产l中国老女人| 午夜免费观看网址| 欧美中文综合在线视频| 欧美黑人欧美精品刺激| 中文在线观看免费www的网站 | 久久精品91蜜桃| 人成视频在线观看免费观看| 色哟哟哟哟哟哟| 一区二区三区精品91| 欧美精品亚洲一区二区| 欧美成人性av电影在线观看| 亚洲熟女毛片儿| 国产成人精品无人区| 搞女人的毛片| 极品教师在线免费播放| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频| 欧美绝顶高潮抽搐喷水| 亚洲精品久久国产高清桃花| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 香蕉久久夜色| av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av| 少妇被粗大的猛进出69影院| 日日爽夜夜爽网站| 日本精品一区二区三区蜜桃| 黄片播放在线免费| 嫩草影院精品99| 国产1区2区3区精品| 亚洲av成人av| 亚洲av电影不卡..在线观看| 亚洲av成人av| 高潮久久久久久久久久久不卡| 成人国产一区最新在线观看| 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 久9热在线精品视频| 无遮挡黄片免费观看| x7x7x7水蜜桃| 午夜福利在线在线| 最近在线观看免费完整版| 身体一侧抽搐| 欧美成狂野欧美在线观看| 欧美三级亚洲精品| 色综合站精品国产| 婷婷亚洲欧美| 一级片免费观看大全| 亚洲一区中文字幕在线| 久久久久国内视频| 精品不卡国产一区二区三区| 久久久久国内视频| 韩国精品一区二区三区| 国产精品美女特级片免费视频播放器 | 欧美精品亚洲一区二区| 一个人观看的视频www高清免费观看 | 亚洲一区高清亚洲精品| 老汉色∧v一级毛片| 可以免费在线观看a视频的电影网站| 欧美又色又爽又黄视频| 97超级碰碰碰精品色视频在线观看| 黄色成人免费大全| 黄色片一级片一级黄色片| 身体一侧抽搐| av福利片在线| 美女大奶头视频| 特大巨黑吊av在线直播 | 热99re8久久精品国产| 精品人妻1区二区| 成年免费大片在线观看| 久久天堂一区二区三区四区| 久久久久精品国产欧美久久久| 午夜福利免费观看在线| 欧美乱码精品一区二区三区| 一区二区三区精品91| 亚洲熟女毛片儿| 国产区一区二久久| 亚洲国产精品999在线| 亚洲免费av在线视频| 国产又爽黄色视频| 午夜久久久在线观看| 亚洲中文字幕日韩| tocl精华| 在线观看www视频免费| 亚洲成av人片免费观看| 老司机在亚洲福利影院| 99久久国产精品久久久| 国内少妇人妻偷人精品xxx网站 | 欧美色欧美亚洲另类二区| av天堂在线播放| 成熟少妇高潮喷水视频| 香蕉av资源在线| 精品日产1卡2卡| 哪里可以看免费的av片| 中文字幕av电影在线播放| 亚洲国产欧美一区二区综合| 白带黄色成豆腐渣| 欧美日韩中文字幕国产精品一区二区三区| 成人国产综合亚洲| 人成视频在线观看免费观看| 亚洲色图av天堂| 中文字幕精品亚洲无线码一区 | 久久久久久大精品| 黑人操中国人逼视频| АⅤ资源中文在线天堂| 黄色毛片三级朝国网站| 成在线人永久免费视频| 黄网站色视频无遮挡免费观看| 俺也久久电影网| 日本成人三级电影网站| 免费看十八禁软件| 黄网站色视频无遮挡免费观看| 久久久精品国产亚洲av高清涩受| 男女午夜视频在线观看| 老司机午夜福利在线观看视频| av中文乱码字幕在线| 制服人妻中文乱码| 亚洲久久久国产精品| 一级毛片女人18水好多| 日本a在线网址| 精品国内亚洲2022精品成人| av在线播放免费不卡| 国产亚洲欧美精品永久| 午夜免费观看网址| 精品久久久久久久久久免费视频| 国产一区二区激情短视频| 婷婷精品国产亚洲av| 久久性视频一级片| 久久久久久免费高清国产稀缺| 午夜福利成人在线免费观看| 国产蜜桃级精品一区二区三区| 久久午夜综合久久蜜桃| 亚洲精品中文字幕一二三四区| 国产精品98久久久久久宅男小说| 久久久久国内视频| 亚洲成人久久性| 最近在线观看免费完整版| 99国产精品一区二区三区| 国产精品美女特级片免费视频播放器 | 丝袜人妻中文字幕| 久久国产精品影院| 99riav亚洲国产免费| 日本成人三级电影网站| 成年女人毛片免费观看观看9| 午夜a级毛片| 成人国语在线视频| 国产一区二区激情短视频| 久久人妻av系列| 欧美黄色片欧美黄色片| 好看av亚洲va欧美ⅴa在| 18禁美女被吸乳视频| 亚洲av成人一区二区三| 俺也久久电影网| 亚洲真实伦在线观看| 成人国语在线视频| 女性生殖器流出的白浆| 怎么达到女性高潮| 宅男免费午夜| 久久久久国产一级毛片高清牌|