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

    Collision off-axis position dependence of relativistic nonlinear Thomson inverse scattering of an excited electron in a tightly focused circular polarized laser pulse

    2024-01-25 07:27:56YuboWang王禹博QingyuYang楊青嶼YifanChang常一凡ZongyiLin林宗熠andYouweiTian田友偉
    Chinese Physics B 2024年1期

    Yubo Wang(王禹博), Qingyu Yang(楊青嶼), Yifan Chang(常一凡), Zongyi Lin(林宗熠), and Youwei Tian(田友偉),?

    1Bell Honors School,Nanjing University of Posts and Telecommunications,Nanjing 210023,China

    2College of Science,Nanjing University of Posts and Telecommunications,Nanjing 210023,China

    Keywords: relativistic nonlinear Thomson inverse scattering, off-axis collision, radiation angle distribution,tightly focused laser pulse

    1.Introduction

    Recent years have witnessed significant advancement in the study of high-power laser–matter interactions,[1–6]particularly within the field of ultrafast laser physics.[7–10]A notable area of focus within this realm is the investigation of relativistic nonlinear Thomson inverse scattering (RNTIS),an important mechanism behind the modulated generation of high-energy x-rays[11,12]and gamma rays.[13–15]The study of RNTIS not only contributes to the broader understanding of laser–matter interactions but also has far-reaching applications in various disciplines,such as biomedicine,[16,17]atomic physics,[18,19]and laser-constrained nuclear fusion.[20]

    RNTIS involves the elastic scattering of photons and electrons, as depicted in Fig.1.Its nonlinear nature is due to the inclusion of the Doppler shift.Earlier studies have examined various aspects of radiation produced by laser pulses colliding with electrons.Zhuanget al.[11]focused on the quasi-monochromatic x-ray spectra generated by the interaction of low-energy electrons and lasers.The work of Tairaet al.[15,21,22]probed the angular distribution of radiative harmonic energy in RNITS.Significant contributions to understanding the dynamics of tightly focused laser-driven electrons were made by Liuet al.[23]Changet al.[24]explored the properties of radiation from tightly focused center colliding RNITS, delving into methods of modulating high-energy gamma rays.

    Most previous research primarily investigates the direct collision of electrons with tightly focused laser pulses.However, in practical settings, it is challenging to constrain the electron to the axis along which the laser pulse is located.Consequently,the characteristics of the electron’s interaction with the laser under off-axis conditions urgently need exploration and understanding.Yuet al.[25]studied the effect of the initial off-axis position of the electron acceleration driven by ponderomotive force in a plane wave laser field.Nonetheless, in the tightly focused Gaussian polarized laser field,the dynamic and nonlinear properties of the electron are significantly influenced by the tight focus characteristics of the strong laser field.Furthermore,the impact of electron off-axis position on space radiation has not been previously investigated.

    This paper presents the first-ever study of the dynamic properties,spatial radiated power,and spectral characteristics of electrons under off-axis collision conditions within a circularly polarized, tightly focused laser field.We investigate the effects of the peak amplitude of the laser pulse and the initial off-axis position of the electron on its dynamic properties,the angular spatial distribution of radiated power,the peak power of spatial radiation,and radiation collimation.

    The research findings suggest that the trajectory of an excited electron’s motion is influenced by both laser intensity and its initial off-axis position.The electron eventually deviates from the laser pulse’s central axis at a constant speed,with the deviation angle’s magnitude affected by the laser intensity and the initial off-axis position of the electron.With relatively low-intensity laser pulses,the radiation radial direction shows significant symmetry aroundφ=0°.As the electron’s initial off-axis distance increases, the spatial radiation peak power decreases, and collimation improves.Additionally, the multimodality of the electron radiation spectrum lessens, supercontinuum phenomena emerge,and the harmonic peak of the radiation spectrum exhibits red-shifts, resulting in a decrease in radiation light intensity.Conversely,under the influence of a relatively high-intensity laser pulse,the spatial radiation power angle distribution of the electron initially increases then decreases as the initial off-axis distance of the electron increases.Similarly, collimation initially weakens and then strengthens.The initial off-axis distance that corresponds to the highest radiation peak power point is identical to that corresponding to the lowest radiation peak angle, under the same laser pulse intensity.The multimodality of the electron radiation spectrum varies in different directions with the change in the electron’s initial off-axis distance.The radiation spectrum harmonic peak shows a distinct blue-shift phenomenon in specific directions.However, with a further increase in the electron’s initial off-axis distance,the radiation spectrum harmonic peak undergoes a red-shift phenomenon in all directions.

    Fig.1.Schematic diagram of relativistic nonlinear Thomson inverse scattering process.

    The remainder of this paper is structured as follows.In Section 2, we derive analytical expressions for the vector potential of the laser pulse,the laws governing electron motion,and the radiation spectrum, using principles of classical electrodynamics.In Section 3,we explore the impact of the peak amplitude of the laser pulse and the initial off-axis position of the electron on the spatial radiation power and radiation spectrum of the electron.Section 4 presents a summary of the effects of the electron’s initial off-axis position and the peak amplitude of the laser pulse on electron dynamics, spatial radiation power, and the radiation spectrum.Additionally, we discuss how to modulate radiation rays to achieve desirable monochromaticity, intensity, and high frequency or high collimation by controlling laser intensity and the initial off-axis position of the electron.

    2.Theory and formula

    To begin,we must note that all subsequent formula definitions have spatial and temporal coordinates normalized by the wave number of the laser,denoted as10 =λ0/(2π),and the frequency of the laser, denoted as=λ0/(2πc).Hereλ0represents the wavelength of the laser,corresponding to 1μm,cis the speed of light,εis the dielectric constant,manderepresent the mass and charge of the electron,respectively.

    We consider a Laguerre–Gaussian(LG)laser pulse propagating along thez-axis at an angle of incidenceσin=0,and with the dielectric media being isotropic,homogeneous,nonmagnetic, and nonconducting.In a tightly focused Gaussian laser field, the electric fieldEand magnetic fieldB, satisfying Maxwell’s equations,can be expressed as follows:[26,27]

    whereAis the solution to the Helmholtz equation

    Within the laser field, circularly polarized laser light can be decomposed into a pair ofx-axis andy-axis linearly polarized lasers with a phase difference ofπ/2.Thus,the electromagnetic field of the circularly polarized laser can also be decomposed similarly,that is,E=Exp+Eyp,B=Bxp+Byp.Yousefet al.[28]derivedA,B, andEin thex-axis linearly polarized laser field from Eqs.(1)–(3).Zhang[29]and Barton[26]proposed the electromagnetic field expression of they-axis linearly polarized laser field using symmetry expression.Combining the above, the electric field componentE={Ex,Ey,Ez}can be derived as

    and the magnetic field components are described as follows:

    whereξ=x/w0,β=y/w0,r=ρ/w0,w0represents the minimum waist radius.The electromagnetic field is fifth-order expanded accurate to the diffraction angleε=w0/zR,andELcan be expressed as follows:

    The electromagnetic field is fifth-order expanded accurate to the diffraction angleε=w0/zR,and among them,SnandCnare shown as follows:

    where?=?0+η ??R+?G,the initial phase?0is the phase of plane wave.The phase related to the curvature of wave fronts?Ris defined as?R=ρ2/(2R(z)),andR(z)=z+z/zindicates the radius of curvature of a wave front intersecting the beam axis at the coordinatez.The Gouy phase is?G=tan?1(z/zr).

    In the calculation process, we set the focus of the laser pulse as the origin,and electrons move toward the laser from{x0,0,z0},wherez0denotes a sufficiently far position,andx0denotes the initial off-axis position of the electron.

    The momentum of the electron in an intense laser pulse can be determined by Lorentz and energy equations given below:

    Here,Γ=γmc2represents the electron energy defined by the Lorentz factorγ=[1?(υ/c)2]?1/2.p=Γ u/c2,υis the electron velocity, andu=υ/c.By solving the ordinary differential Eqs.(4)–(9) with Runge–Kutta method through MATLAB,the numerical solution of the electron’s position,velocity,and acceleration are recorded for each step.

    The electron dynamics can be explained using the ponderomotive potential model.[30]Linet al.[31]suggested that for free electrons moving in a laser pulse, the mass force is the time average force experienced by the electron,that is,the low-frequency force generated by the high-frequency field:[31]

    whereVPondis the ponderomotive potential and can be expressed as follows:

    Here,a=eE/(mc)is the normalized amplitude of the electromagnetic field of the laser.As this paper is mainly a qualitative analysis,in order to simplify the calculation,the first-order approximation of the electric field is substituted intoato obtain the first-order approximate direct proportion of the transverse and longitudinal ponderomotive forces:

    From knowledge of electrodynamics we know that the electron in relativistic acceleration emits electromagnetic radiation, and the radiated power per unit steradian angle can be expressed as follows:[32,33]

    Here, the power radiated per unit solid angleP?is normalized and the direction of radiationn={sin(θ)cos(φ),sin(θ)sin(φ),cos(θ)}, whereθrepresents the polar angle to the laser movement direction, andφrepresents the azimuth angle on the plane perpendicular to origin point,t=t'+R,R~R0?n·r,wheret'is the time at which the electron interacts with the laser pulse, andR0denotes the distance from the origin to the observer,ris the position vector of the electron.The formula for the radiation energy per unit solid angle per unit frequency interval during the interaction of the electron with the laser pulse can be expressed using the following equation:[34,35]

    where d2I/dωd?is normalized usinge2/(4π2c),s=ωs/ω0,andωsis the frequency of the harmonic radiation.By solving Eqs.(20) and (21), the full time, full space and full spectral properties of the electron harmonic radiation can be obtained.

    3.Numerical results

    In our observations,we focused on radiation emitted from a sphere with a radius of 1 m,centered at the origin of the coordinates.The beam waist radiusw0of the laser pulse is set to 2λ0, which corresponds to 2 μm, and the pulse durationLis 6.6 fs.The initial phase?0=0.The electron initially carries energy denoted byγ=5(equivalent to 2.56 MeV),moving in the direction opposite to thez-axis.This energy level can be achieved through acceleration via a linear accelerator.Given the characteristics of the circularly polarized laser pulse, the effects of the electron’s off-axis collision with the laser pulse on the RNITS are modulated by the initial phase?0of the laser pulse.Specifically, the impact on the electron from a pulse with an initial phase?0at the radial off-axis angleα0is equivalent to the effect experienced by an electron subjected to a pulse with an initial phase?0+α1?α0at the radial offaxis angleα1.As such, the point of collision between the electron and the laser pulse, that is, the maximum value of the electron motion amplitude is positioned on thex-axis of the Cartesian coordinate system.For the purposes of our discussion, the radiated peak power per unit solid angle will be denoted as dp/d?.It is important to note that all numerical results presented here are derived based on the physical laws and equations outlined in Section 2.

    3.1.Electronic motion properties

    As illustrated in Figs.2(a)–2(d), with increasing laser pulse intensity, the final axial deviation extent of the off-axis electron’s motion trajectory is more pronounced.Additionally,an electron with an initial off-axis position ofx0/=0 deviates from thez-axis direction at a constant speed after the primary interaction stage with the laser pulse.The final deviation speed is closely tied to the intensity of the laser pulse,as depicted in Figs.2(e)–2(h).

    This phenomenon can be attributed to the fact that,when the electron centrally collides with the laser pulse, the transverse ponderomotive force experienced by the electron at each directional angle is equal over a macroscopic time scale.As a result,the electron follows a circular-like spiral motion,observable in Figs.2(a)–2(d).

    When an electron moves off-axis,it is subjected to a slight decay of the transverse ponderomotive force in the off-axis direction, along with a relatively larger ponderomotive force in other directions.This can be seen in Figs.2(i)–2(l).In extreme cases(e.g.,a0=10 as shown in Fig.2(l)),the ponderomotive force experienced by the electrons becomes highly irregular.This irregularity in force induces a change in the trajectory and transverse velocity of the electron’s motion.As the laser intensity increases,so too does the transverse ponderomotive force on the electron, amplifying the aforementioned effects.Consequently, the degree of axial deviation of the electron with respect to its initial off-axis position naturally escalates.

    Moreover, it is noteworthy from Figs.2(a)–2(d) that the final axial deviation angle of the electron’s trajectory does not have a direct correlation with the electron’s initial off-axis positionx0.Specifically, when thea0<10, the final axial deviation angle of the electron’s trajectory is at its maximum when the initial off-axis positionx0=1λ0.However, whena0=10, the maximum final axial deviation angle of the trajectory occurs when the initial off-axis position of the electronx0=0.5λ0.

    This phenomenon can be attributed to the limitation of the laser pulse’s beam waist on the effective transverse ponderomotive action of the laser pulse on the off-axis electron.In this study,the beam waist radius of the laser pulseb0=2λ0.The ponderomotive force’s impact on the electron in the longitudinal direction is significantly smaller than that in the transverse direction,as depicted in Figs.2(i)–2(p).Hence,the electron’s longitudinal motion speed is primarily influenced by its initial energy.

    Fig.2.Effect of the electron’s initial off-axis position x0 on the electron trajectories under the action of laser pulses with different laser peak amplitudes a0,where(a)–(d)illustrate the longitudinal changes of the electron motion compared to the initial position,along with the projection deviation angle (the connection between the origin and the point z ≤?100λ0), and (e)–(h) represent the distances of the electrons from the center axis of the laser pulse at different moments.(i)–(l)represent the transversal ponderomotive factors(the term after the proportional sign in Eqs.(17)and(18)),(m)–(p)represent the longitudinal ponderomotive factors(the term after the proportional sign in Eq.(19)).

    Indeed, this study primarily conducts a qualitative analysis of the ponderomotive force, and the results shown in Figs.2(i)–2(p)represent the first-order approximate direct proportional quantitiesαxPond,αyPondandαzPondin Eqs.(14)–(16).These can be interpreted as normalized first-order approximations of the ponderomotive force.This approach allows us to compare the magnitudes of the ponderomotive forces experienced by electrons under different conditions.

    The degree of transverse deviation of the electron can be characterized by the following two measures:

    (i) The ratio of the transverse ponderomotive force exerted on the electron in the off-axis direction to that in the on-axis direction.

    (ii)The magnitude of the total transverse ponderomotive force applied to the electron.

    As can be observed from Figs.2(i)–2(l), the first aspect demonstrates a positive correlation with the electron’s initial off-axis positionx0,while the latter shows a negative correlation with the same.This contrasting behavior creates a situation where the electron’s initial off-axis positionx0does not exhibit a straightforward positive or negative correlation with the final axial deviation of the electron’s trajectory.Instead,this relationship is significantly influenced by the intensity of the laser pulse.

    3.2.Electron space relativistic nonlinear Thomson inverse scattering properties

    3.2.1.Peak power properties of electron space radiation

    It is readily apparent from Figs.3(a) and 3(b) that when the intensity of the laser pulse is relatively low(a0≤3),the angular distribution of the electron radiation power displays an annular shape with a hollow center.This finding aligns with Tairaet al.’s[19]studies of higher harmonics.

    In this study, the beam waist radius of the laser pulse,b0=2λ0, placing it in the category of highly tightly focused laser pulses.The action of the transverse ponderomotive force of the laser field is tightly confined to a relatively close distance from the laser pulse axis.Additionally,the electron under examination does not fall into the category of high-energy particles, as it has a relatively low longitudinal velocity compared to high-energy electrons.This results in the laser field’s period of action being relatively long, which means that the overall action of the transverse ponderomotive force on the electron is considerably higher than that applied to ultra-highenergy electrons in a plane wave laser field.Consequently,the divergence of its radiation, or the angular distribution of the first harmonic peak of the radiation,is quite distinct from the distribution of higher harmonics in the plane wave laser field,as studied by Tairaet al.[19]

    As shown in Figs.2(i)–2(l),when the laser pulse peak amplitudea0≤5,the longer the electron’s initial off-axis positionx0,the lower the transverse ponderomotive force experienced by the electron,and the better the collimation of the electron’s spatial radiation.This is evident even whenx0= 2.8λ0as shown in Fig.3(a), where the annular shape distribution of the electron radiation power disappears and reverts to a single peaked distribution.However, as the laser pulse intensity increases,the collimation of the electron’s spatial radiation deteriorates and the asymmetry of the radiation becomes increasingly pronounced for the samex0.

    As shown in Figs.4(a)and 4(b),when the peak amplitude of the light pulsea0≥5,and the initial off-axis position of the electronx0<1λ0,the spatial radiation collimation of off-axis electrons is weaker than when the initial position of the electron is on the central axis of the laser pulse.This indicates that,under the influence of a strong laser pulse,the spatial radiation collimation of the electron is not positively correlated with the initial off-axis position of the electron,x0.The explanation for this phenomenon will be discussed subsequently.

    Fig.3.Effect of the distance x0 of the electron initial position collision from the central axis of the laser pulse on the maximum power distribution per unit angle of spatial radiation under the action of relatively low laser pulses(a0=1(a)and a0=3(b)).

    Fig.4.Effect of the distance x0 of the electron initial position collision from the central axis of the laser pulse on the maximum power distribution per unit angle of spatial radiation under the action of relatively high laser pulses(a0=5(a)and a0=10(b)).

    When the laser pulse intensity surpasses the relativistic intensity,the spatial radiation symmetry of the electron breaks as the initial off-axis positionx0increases,exhibiting a vortex state for the first time in the collision scenario.This occurs due to the instability of the transverse ponderomotive force acting on the electron during off-axis motion, which leads to an unstable transverse drift motion.The radius of curvature frequently changes in the strong laser field, which subsequently causes significant alterations to the electron’s spatial radiation and disrupts the symmetry of the power distribution.

    When the initial off-axis positionx0of the electron exceeds the beam waist radius of the laser pulse, the transverse ponderomotive force acting on the electron becomes too small.Consequently,the unstable transverse drift motion is no longer as pronounced compared to the longitudinal motion, and the instability of the electron’s transverse angular radiation weakens.This results in the disappearance of its vortex state and an improvement in collimation.

    When the laser pulse peak amplitudea0=10 and the initial off-axis positionx0just begins to deviate from the central axis of the laser pulse, the primary radiation direction of the electron surprisingly shifts to the forward direction.This implies that the lateral drift motion of the electron is highly unstable in an ultra-intense laser field(a0≥10),and the asymmetry exacerbates to the extent that it alters the basic characteristic of backward radiation.Therefore,in practical experiments,when the peak amplitude of the laser pulsea0≥10,researchers must consider observing the spatial radiation of off-axis electrons from the forward direction.Failing to do so may lead to incorrect conclusions regarding the angular power of the radiation unit.

    Next, we delve deeper into the instantaneous maximum power represented by max(dp/d?)and the off-axis angleθ0at the extreme value of the electron’s spatial radiation.These parameters respectively reflect the instantaneous maximum power of the electron’s spatial radiation and the collimation of the radiation, both of which are crucial for experimental purposes.Figures 5(a)and 5(b)depict the relationships between these parameters and the laser intensity as well as the initial off-axis position of the electron.

    Fig.5.Effect of the distance x0 of the initial position of the electron collision off the central axis of the laser pulse on the logarithmic value of the peak unit angular power of radiation log10[max(d p/d?)] (a) and the peak angle of radiation θ0 (b) under the action of laser pulses with different peak amplitudes a0,where the red lines mark the initial off-axis position x0 and value corresponding to the maximum peak radiation(a)and the minimum peak angle of radiation(b),respectively,under the same a0 condition.

    One key point to note is that as the laser pulse intensity intensifies, the relationship between the peak power magnitude of the electron’s spatial radiation and the electron’s initial off-axis position,x0, shifts from being negatively correlated to first increasing and then decreasing.This suggests that the maximum power radiation of the electron under the influence of a high-energy laser pulse is not associated with an on-axis collision,but instead with an off-axis collision.

    The explanation for this phenomenon can be drawn from the discussion about the ponderomotive force on off-axis electrons in Subsection 3.1.When the laser pulse is relatively weak (a0<5 in this case, but more specificallya0<4.5),the transverse ponderomotive force of the laser pulse is rather small,as shown in Figs.2(i)and 2(j).Given the electron’s high longitudinal velocity, the impact of the asymmetric ponderomotive force on its motion is relatively weak.As the electron strays further and further from the laser pulse’s central axis,its transverse ponderomotive force decreases,leading to a corresponding reduction in the power of the electron’s excited radiation.

    When the laser pulse is relatively intense (a0>4.5), as studied in this paper), the asymmetry of the ponderomotive force on the electron moving off-axis cannot be ignored.As illustrated in Figs.2(k) and 2(l), the magnitude of the ponderomotive force in the off-axis direction is greater than that in the direction opposite to the off-axis.Consequently, over the macroscopic time scale,the electron experiences an additional net ponderomotive force in the off-axis direction.The ratio of this force to the combined ponderomotive force rises with an increase inx0.

    However, due to the constraints of the laser beam waist,the combined ponderomotive force diminishes asx0increases.Under the combined effect of both the ratio and combined force, when the initial off-axis positionx0is relatively small butx0/=0,the net effect of the combined transverse ponderomotive force over the macroscopic time scale still exhibits an increasing trend asx0gradually increases.This is despite the rising transverse asymmetric ponderomotive force ratio and the declining combined ponderomotive force.

    Asx0increases further, the net effect of the combined transverse ponderomotive force over the macroscopic time scale diminishes due to the even faster decline of the transverse ponderomotive force from all directions.Although the ratio of the transverse asymmetric ponderomotive force on the electron continues to increase,it cannot compensate for the drastic decrease in the overall ponderomotive force.This results in a decline in the net combined transverse ponderomotive force,which in turn leads to a decrease in the peak power radiated from the electron in space.

    Additionally, a detailed examination of Fig.5(b) reveals that the numerical distribution in Figs.5(b) and 5(a) exhibits clear symmetry.Even the distribution of the extreme point within the unit laser peak amplitude is the same,which means that the patterns in these two aspects correspond to each other.

    As the intensity of the laser pulse increases,the relationship between the peak power angle of the electron’s spatial radiation and the electron’s initial off-axis positionx0gradually transitions from a positive correlation to a pattern that first decreases and then increases.This implies that the collimation of the electron’s spatial radiation is negatively correlated with the magnitude of the maximum power under identical initial electron energy,laser pulse width,and beam waist radius conditions.This means we cannot achieve electron radiation with both high collimation and high power simultaneously by simply modifying the laser pulse intensity and the initial off-axis position of the electron.The explanation for this can be found by considering electron dynamics.The electron, subjected to the action of the transverse ponderomotive force, undergoes energy changes and subsequently generates radiation.The power magnitude of the radiation is positively correlated with the degree of action of the transverse ponderomotive force.However, the action of the transverse ponderomotive force also leads to instability in the electron’s transverse drift.Therefore,under the same basic conditions of initial electron energy,laser pulse width,and beam waist radius,the collimation of the electron’s spatial radiation and its power magnitude are mutually exclusive.

    3.2.2.Spectral properties of electron radiation

    The deviation of the electron radiation spectrum from the symmetry axis of the laser pulse,along with its light intensity and harmonic frequency characteristics, were assessed.The results are presented in Figs.6 and 7.

    Figures 6(a)and 6(b)clearly shows that at relatively low laser intensities(a0<5),the spatial radiation spectrum of the electron gravitates toward the central axis of the laser pulse as the initial off-axis positionx0increases.This suggests a strengthening of the collimation of the electron’s spatial radiation spectrum.Simultaneously,as the initial off-axis positionx0of the electron increases, the harmonic peak of the electron radiation spectrum starts to red-shift.The wave number is gradually increasing,and the harmonic light intensity of the spectrum also diminishes.

    This behavior arises from the relatively low joint transverse ponderomotive force acting on the off-axis electron in low-intensity laser fields.As the initial off-axis distance of the electron increases,the transverse ponderomotive force decreases, resulting in a reduction of the transverse drift motion.Consequently, the collimation of the electron spectrum improves,and the harmonic light intensity decreases.

    Moving on to Figs.7(a) and 7(b), we observe that in the case of relatively strong laser pulses, similar to the spatial radiation properties discussed earlier, the radiation spectral properties of the electron change with an increase in the initial off-axis positionx0.Initially, the spectrum shifts toward higher frequencies and energies, becomes multimodal,and subsequently transitions to a state of under-collimation,lower frequencies, lower energies, and super-continuum generation with higher collimation.

    Fig.6.Spatial radiation spectrum θ angle distribution when the distance x0 of the electron’s initial position of the collision deviates from the center axis of the laser pulse varies under the action of relatively low laser pulses(a0=1(a)and a0=3(b)).

    Fig.7.Spatial radiation spectrum θ angle distribution when the distance x0 of the electron’s initial position of the collision deviates from the center axis of the laser pulse varies under the action of relatively high laser pulses(a0=5(a)and a0=10(b)).

    In particular,when the laser peak amplitudea0=10,the first harmonic peak of the electron at an initial off-axis positionx0=2.4λ0(with an observedθangle of approximately 80°)can reach a frequency of 1.8×1020Hz.This frequency range corresponds to high-frequency gamma rays and holds potential for applications in cancer treatment through radiation therapy.Subsequently, asx0increases further, the first harmonic peak of the electron radiation starts to red-shift.Forx0≥2.4λ0,the first harmonic peak no longer falls within the range of gamma rays,but rather within the range of x-rays.

    At the electron’s initial off-axis positionx0=2.8λ0, the collimation of the radiation spectrum is notably enhanced compared to when the electron is closer to the axis.The peak harmonic distribution of the radiation is concentrated within the range of 180°±15°.Additionally, the attenuation of the intensity in the above electron radiation spectrum is considerably limited within this harmonic distribution range due to the superposition of interference between the harmonics.

    Based on the observations from Figs.6 and 7, we can draw the following conclusion: when the laser pulse peak amplitudea0is relatively low,the electron radiation spectrum exhibits a noticeable axial symmetry with respect to the laser pulse axis.However,as the laser pulse peak amplitudea0increases, this axial symmetry is broken, resulting in an asymmetric frequency distribution of the first harmonic peak even with a small increase in the initial off-axis positionx0.Specifically,the harmonic peak shows a blue-shifted trend in theX+direction and a red-shifted trend in theX?direction.Asx0continues to increase,the collimation of the electron radiation improves,leading to red-shifted harmonic peaks and a restoration of symmetry in the electron radiation spectrum about the laser pulse axis.

    Furthermore, the convergence of electron radiation harmonics toward the central axis of the laser pulse is evident in Figs.6 and 7.As this convergence occurs, the multimodal nature of the radiation spectrum diminishes, giving way to significant interference superposition between the harmonics.This results in the electron radiation harmonics displaying distinct super-continuum spectral properties.This effect is particularly prominent when the initial off-axis position of the electrons is beyond the waist of the laser pulse beam,that is,whenx0>2λ0.

    Our analysis included measurements of the radial distribution, light intensity, and harmonic frequency properties of the electron radiation spectrum,as depicted in Figs.8 and 9.

    From Figs.8(a) and 8(b) it is evident that the radial distribution of the electron radiation exhibits significant symmetry aroundφ=0°in the case of a low-intensity laser field.This symmetry arises due to the relatively small transverse ponderomotive force acting on the electrons and the resulting weak influence of transverse drift motion, which leads to pronounced transverse symmetry.

    Fig.8.Angular distribution of the spatial radiation spectrum φ for the variation of the distance x0 from the central axis of the laser pulse at the initial position of the electron collision under the action of the laser pulse with relatively low peak amplitude(a0=1(a),a0=3(b)).

    Furthermore, when the laser pulse intensity is relatively low, the harmonic peaks nearφ=0°andφ=180°exhibit a red-shift as the initial off-axis positionx0of the electron increases.However, there is no significant frequency change observed in the harmonic peaks nearφ=90°andφ=270°.Additionally,the intensity of the radiation spectrum decreases as the initial off-axis positionx0of the electron increases.The explanation for this behavior has been previously discussed and is not reiterated here.

    However,when the off-axis electron is subjected to an intense laser pulse, as shown in Figs.9(a) and 9(b), the radial distribution of the radiation spectrum loses its symmetry.The distribution of harmonic peaks becomes concentrated in the rangeφ ∈(?90°,90°),with a significant decrease in harmonic light intensity nearφ=180°.The explanation for this phenomenon can be linked to the earlier discussion in Figs.7(a)and 7(b),where it was observed that the harmonic peaks experience a blue-shift followed by a red-shift as the initial off-axis positionx0of the electron increases.Additionally, due to the asymmetry of the ponderomotive force exerted by the intense laser pulse, the electron undergoes unstable transverse drift motion.This instability contributes to the loss of symmetry in the radiation spectrum distribution, with dominance from the direction of transverse drift motion.However,as the initial position of the electron moves further away from the central axis of the laser pulse, the unstable transverse drift motion weakens, resulting in a restoration of symmetry in the distribution of the electron’s spatial radiation spectrum.

    In addition,figures 9(a)and 9(b)reveal that in the strong laser field, the off-axis electron radiation spectrum exhibits two main radiation directions.These directions are characterized by the main harmonic direction angles:φmain∈(?90°,0°), which are called referred to asY?direction harmonics,andφmain∈(0°,90°),referred to asY+direction harmonics in this paper.When the laser pulse peak amplitudea0<10,the harmonic peaks of the electron radiation are distributed aroundφ=0°due to the interference superposition of theY?andY+direction harmonics.However,when the laser pulse peak amplitudea0=10, the multimodality of the electron radiation spectrum becomes significantly enhanced.As a result, theY?andY+ direction harmonics no longer appear in interference superposition, leading to a different distribution of the harmonic peaks in the electron radiation spectrum away fromφ=0°.

    In addition, in Figs.9(a)and 9(b), theY?direction harmonics of the off-axis electron generally exhibit higher frequencies compared to theY+ direction harmonics.Furthermore, as the initial off-axis distancex0increases, the frequency difference between theY?andY+direction harmonics gradually decreases.This indicates that the red-shift and blue-shift effects on theY?direction harmonics with changes in the initial off-axis distancex0are more pronounced compared to the frequency changes observed in theY+ direction harmonics.

    Fig.9.Angular distribution of the spatial radiation spectrum φ for the variation of the distance x0 from the central axis of the laser pulse at the initial position of the electron collision under the action of the laser pulse with relatively high laser pulses(a0=5(a)and a0=10(b)).

    4.Conclusion and perspectives

    In this study,we have investigated the variations in electron dynamics,radiated power,and spectral properties in RNTIS.Specifically,we focused on the collisions between tightly focused circular polarized lasers and electrons considering the influence of laser intensity and the initial off-axis position of the electrons.From the perspective of electron dynamics,we observed that off-axis electrons experience deviations from axial motion due to the combined effect of the transverse ponderomotive force.The magnitude of the deviation angle and speed is influenced by the initial off-axis position of the electrons and the peak amplitude of the laser pulse.

    At relatively low laser pulse intensities(a0<4.5),we observed a monotonic variation in electron radiation properties with the initial off-axis position of the electron.The angular distribution of the spatially radiated power decreases as the initial off-axis distance of the electronx0increases, indicating an enhancement in collimation.The multimodality of the electron radiation spectrum weakens, and a phenomenon of super-continuous spectrum appears.Additionally, the harmonic peaks of the radiation spectrum are red-shifted, resulting in a decrease in radiation intensity.Notably,the radiation radial direction exhibits significant symmetry aboutφ=0°.

    At relatively high laser pulse intensities (a0>4.5), the electron spatial radiation properties exhibit a distinct variation with the initial off-axis position of the electron.This variation is attributed to the transverse asymmetric ponderomotive force acting on the electron.The spatial radiation power angle distribution of the electron initially increases and then decreases asx0increases,while the collimation first decreases and then increases.Interestingly, the value ofx0corresponding to the highest peak radiation power is the same as the value ofx0corresponding to the lowest peak radiation angle for a given laser pulse intensity.The multimodal nature of the electron radiation spectrum behaves differently in different directions asx0varies, and the radiation spectrum exhibits a blue-shift phenomenon in specific directions.Furthermore,asx0further increases,the radiation spectrum harmonic peaks experience a red-shift phenomenon in all directions.

    In summary, to obtain high-energy, highly monochromatic,and high-frequency gamma rays within certain parameters such as the laser pulse beam radius,pulse width,and electron initial energy, the laser pulse intensity can be increased.Whena0>4.5, it is possible to deliberately deviate the electron from the central axis of the laser during the initial phase.However, this results in relatively weak collimation of the electron radiation.On the contrary, if high collimation of the electron radiation is desired,the initial off-axis position of the electron can be significantly increased, while simultaneously reducing the laser pulse intensity.However, this approach leads to lower power, weaker monochromaticity, and lower frequency of the electron radiation.These findings provide a valuable method for utilizing RNITS as a radiation source for scientific experiments and applications, including cancer treatment.By optimizing the laser pulse intensity and the initial off-axis position of the electron, researchers can strike a balance between power, collimation, monochromaticity, and frequency to meet their specific requirements.

    Acknowledgments

    Project supported by the National Natural Science Foundation of China (Grant Nos.10947170/A05 and 11104291),the Natural Science Fund for Colleges and Universities in Jiangsu Province(Grant No.10KJB140006),the Natural Sciences Foundation of Shanghai(Grant No.11ZR1441300),and the Natural Science Foundation of Nanjing University of Posts and Telecommunications (Grant No.NY221098).We also thank the Jiangsu Qing Lan Project for their sponsorship.

    久久午夜亚洲精品久久| www.999成人在线观看| 新久久久久国产一级毛片| 十八禁高潮呻吟视频| 伊人久久大香线蕉亚洲五| 免费日韩欧美在线观看| 91老司机精品| 最新美女视频免费是黄的| 免费在线观看日本一区| 欧美乱色亚洲激情| 国产乱人伦免费视频| 日韩中文字幕欧美一区二区| 日韩成人在线观看一区二区三区| 久久精品aⅴ一区二区三区四区| 国产欧美日韩一区二区三区在线| 免费黄频网站在线观看国产| 男女床上黄色一级片免费看| 大型av网站在线播放| 在线观看免费日韩欧美大片| 欧洲精品卡2卡3卡4卡5卡区| 色综合婷婷激情| 曰老女人黄片| 一进一出抽搐gif免费好疼 | 窝窝影院91人妻| 欧美日韩国产mv在线观看视频| 夜夜躁狠狠躁天天躁| 在线av久久热| 一进一出抽搐gif免费好疼 | 视频在线观看一区二区三区| 国产午夜精品久久久久久| 午夜福利一区二区在线看| 在线观看一区二区三区激情| 精品国产亚洲在线| 在线观看免费高清a一片| 在线播放国产精品三级| 国产97色在线日韩免费| 国产精品秋霞免费鲁丝片| 一级,二级,三级黄色视频| 香蕉丝袜av| 成人免费观看视频高清| 久久亚洲精品不卡| 久久狼人影院| 亚洲国产毛片av蜜桃av| 欧美 日韩 精品 国产| 91成人精品电影| 日韩欧美一区视频在线观看| 日韩三级视频一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品在线美女| 亚洲精品久久成人aⅴ小说| 亚洲精品美女久久av网站| 香蕉丝袜av| 亚洲一码二码三码区别大吗| 国产成人精品久久二区二区91| 如日韩欧美国产精品一区二区三区| 在线av久久热| 伦理电影免费视频| 91在线观看av| 日本黄色视频三级网站网址 | 亚洲熟女精品中文字幕| 色婷婷av一区二区三区视频| 亚洲精品久久午夜乱码| 女同久久另类99精品国产91| 国产真人三级小视频在线观看| 亚洲精品美女久久av网站| 啪啪无遮挡十八禁网站| 久99久视频精品免费| 欧美中文综合在线视频| 无限看片的www在线观看| 激情视频va一区二区三区| 亚洲少妇的诱惑av| www.999成人在线观看| 大型av网站在线播放| 人人妻人人澡人人看| 欧美国产精品一级二级三级| 国产精品久久久久成人av| 少妇粗大呻吟视频| 精品乱码久久久久久99久播| 亚洲av日韩在线播放| 久久精品成人免费网站| 欧美日韩亚洲国产一区二区在线观看 | 欧美成狂野欧美在线观看| 又黄又粗又硬又大视频| 亚洲av第一区精品v没综合| 成人av一区二区三区在线看| 国产乱人伦免费视频| 黑人操中国人逼视频| 久久 成人 亚洲| 我的亚洲天堂| 中文字幕色久视频| 成人三级做爰电影| 亚洲五月婷婷丁香| 91国产中文字幕| 如日韩欧美国产精品一区二区三区| 欧美日韩亚洲高清精品| 亚洲午夜理论影院| 身体一侧抽搐| 18禁国产床啪视频网站| 麻豆av在线久日| 最新在线观看一区二区三区| 色播在线永久视频| 日日夜夜操网爽| 久久精品人人爽人人爽视色| 亚洲少妇的诱惑av| 久久久久久免费高清国产稀缺| 97人妻天天添夜夜摸| av一本久久久久| 欧美黑人精品巨大| 国产乱人伦免费视频| 亚洲一区中文字幕在线| 咕卡用的链子| 国产成+人综合+亚洲专区| 9191精品国产免费久久| 国产免费男女视频| 99精品在免费线老司机午夜| 性少妇av在线| 18禁裸乳无遮挡动漫免费视频| 亚洲人成电影免费在线| 中国美女看黄片| 18禁美女被吸乳视频| 免费人成视频x8x8入口观看| 中文字幕人妻丝袜制服| 99国产综合亚洲精品| 久久精品人人爽人人爽视色| videosex国产| 黄色视频不卡| 国产精品二区激情视频| 天堂动漫精品| 成人免费观看视频高清| 性色av乱码一区二区三区2| 岛国毛片在线播放| 午夜亚洲福利在线播放| 日韩一卡2卡3卡4卡2021年| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美精品综合一区二区三区| 12—13女人毛片做爰片一| 久久国产乱子伦精品免费另类| www.999成人在线观看| 1024视频免费在线观看| 国产精品 欧美亚洲| 黑人猛操日本美女一级片| 在线国产一区二区在线| 好男人电影高清在线观看| 久久精品国产清高在天天线| 久久亚洲真实| 一级片'在线观看视频| 亚洲成人手机| 国产99白浆流出| av欧美777| 中文字幕人妻丝袜制服| av片东京热男人的天堂| 亚洲成人手机| 国产亚洲精品第一综合不卡| 一进一出好大好爽视频| 欧美精品一区二区免费开放| 看免费av毛片| 午夜视频精品福利| 久久久水蜜桃国产精品网| 国产黄色免费在线视频| 国产高清国产精品国产三级| 一级黄色大片毛片| 日本精品一区二区三区蜜桃| 国产色视频综合| 久久ye,这里只有精品| 欧美性长视频在线观看| 国产男靠女视频免费网站| 久久中文看片网| 母亲3免费完整高清在线观看| 久久久久久亚洲精品国产蜜桃av| 搡老熟女国产l中国老女人| 每晚都被弄得嗷嗷叫到高潮| 一个人免费在线观看的高清视频| 精品国产一区二区久久| 在线国产一区二区在线| 日韩制服丝袜自拍偷拍| 露出奶头的视频| 亚洲精品中文字幕一二三四区| 日韩三级视频一区二区三区| 久久国产精品人妻蜜桃| 91大片在线观看| 国产在线精品亚洲第一网站| 精品少妇久久久久久888优播| 我的亚洲天堂| e午夜精品久久久久久久| 亚洲精品国产一区二区精华液| 免费观看人在逋| 色精品久久人妻99蜜桃| 涩涩av久久男人的天堂| 欧美中文综合在线视频| 国产单亲对白刺激| bbb黄色大片| 一区在线观看完整版| 国产97色在线日韩免费| 午夜久久久在线观看| 淫妇啪啪啪对白视频| 午夜91福利影院| 色老头精品视频在线观看| 久久国产精品影院| 999久久久精品免费观看国产| 天天影视国产精品| 好男人电影高清在线观看| 嫩草影视91久久| 韩国精品一区二区三区| 久久久水蜜桃国产精品网| 亚洲一卡2卡3卡4卡5卡精品中文| av一本久久久久| 18禁黄网站禁片午夜丰满| 丝瓜视频免费看黄片| 一级毛片精品| 窝窝影院91人妻| 亚洲一区二区三区不卡视频| 成年版毛片免费区| 久久香蕉精品热| 精品电影一区二区在线| 无遮挡黄片免费观看| 国产91精品成人一区二区三区| 国产一区二区三区综合在线观看| 亚洲性夜色夜夜综合| 久久亚洲真实| 69精品国产乱码久久久| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 大型黄色视频在线免费观看| 亚洲人成77777在线视频| 一级作爱视频免费观看| 国产日韩一区二区三区精品不卡| 国产精品综合久久久久久久免费 | 露出奶头的视频| 亚洲人成伊人成综合网2020| 热re99久久国产66热| 亚洲av成人不卡在线观看播放网| 国产精品自产拍在线观看55亚洲 | 麻豆成人av在线观看| 一a级毛片在线观看| 国产有黄有色有爽视频| 日本vs欧美在线观看视频| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 深夜精品福利| 91av网站免费观看| 啦啦啦 在线观看视频| 在线天堂中文资源库| 叶爱在线成人免费视频播放| 人人妻人人澡人人爽人人夜夜| 国产精品香港三级国产av潘金莲| 亚洲美女黄片视频| 国产淫语在线视频| 国产精品九九99| 国产成人精品久久二区二区91| 一区在线观看完整版| 大片电影免费在线观看免费| 成人18禁高潮啪啪吃奶动态图| 好看av亚洲va欧美ⅴa在| 又黄又爽又免费观看的视频| 久久久精品区二区三区| 身体一侧抽搐| 美女午夜性视频免费| 国产精品久久久人人做人人爽| 亚洲精品美女久久久久99蜜臀| 国产精品自产拍在线观看55亚洲 | av网站免费在线观看视频| 亚洲五月天丁香| 久久精品91无色码中文字幕| 国产欧美日韩一区二区精品| 精品久久久久久电影网| 欧美日韩视频精品一区| 首页视频小说图片口味搜索| 免费在线观看日本一区| 国产精品99久久99久久久不卡| 动漫黄色视频在线观看| 国产一区二区三区视频了| 亚洲成人免费av在线播放| 国产不卡av网站在线观看| 国产一区二区三区在线臀色熟女 | 欧美大码av| av电影中文网址| 黄色成人免费大全| 亚洲欧美一区二区三区久久| 欧美成人午夜精品| 久久久精品区二区三区| 亚洲一区二区三区不卡视频| 久久久久国产精品人妻aⅴ院 | 久久热在线av| 国产视频一区二区在线看| 精品视频人人做人人爽| 波多野结衣av一区二区av| 少妇裸体淫交视频免费看高清 | 久久久国产成人免费| 精品无人区乱码1区二区| 伦理电影免费视频| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 嫁个100分男人电影在线观看| 久久国产亚洲av麻豆专区| 国产一区有黄有色的免费视频| 免费看十八禁软件| 热re99久久国产66热| 亚洲av美国av| 欧美乱色亚洲激情| 国产主播在线观看一区二区| 亚洲五月婷婷丁香| 国产精品电影一区二区三区 | 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 国产三级黄色录像| 国产97色在线日韩免费| 国产99白浆流出| 亚洲九九香蕉| 视频在线观看一区二区三区| 午夜福利影视在线免费观看| 亚洲av美国av| 欧美日韩亚洲国产一区二区在线观看 | 午夜视频精品福利| www.自偷自拍.com| 天天躁夜夜躁狠狠躁躁| 国产精品影院久久| 国产亚洲欧美在线一区二区| 久99久视频精品免费| 丰满迷人的少妇在线观看| 岛国在线观看网站| 亚洲成人免费av在线播放| 99精国产麻豆久久婷婷| 女警被强在线播放| 两个人看的免费小视频| 黄色成人免费大全| 午夜福利乱码中文字幕| 午夜福利一区二区在线看| 电影成人av| 99久久综合精品五月天人人| av电影中文网址| 成人特级黄色片久久久久久久| 国产免费av片在线观看野外av| 日韩欧美免费精品| 欧美日韩视频精品一区| √禁漫天堂资源中文www| 老司机深夜福利视频在线观看| 欧美精品亚洲一区二区| 欧美另类亚洲清纯唯美| 亚洲午夜理论影院| 欧美日韩亚洲高清精品| 亚洲精品一二三| 丝袜人妻中文字幕| 国产精品免费大片| av不卡在线播放| tube8黄色片| 中文字幕另类日韩欧美亚洲嫩草| 中文欧美无线码| 美女扒开内裤让男人捅视频| netflix在线观看网站| 18禁观看日本| 欧美激情极品国产一区二区三区| 99精品久久久久人妻精品| 黄色a级毛片大全视频| 又黄又爽又免费观看的视频| 丁香欧美五月| 免费观看人在逋| 国产成人av教育| 在线国产一区二区在线| 亚洲国产毛片av蜜桃av| 精品免费久久久久久久清纯 | 91精品国产国语对白视频| 欧美在线黄色| 亚洲av欧美aⅴ国产| 淫妇啪啪啪对白视频| 国产在线观看jvid| 亚洲精品中文字幕一二三四区| xxx96com| 免费黄频网站在线观看国产| 国产色视频综合| 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| 欧美中文综合在线视频| 精品人妻1区二区| 大香蕉久久成人网| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 精品久久蜜臀av无| 动漫黄色视频在线观看| 久久中文字幕一级| 黄色怎么调成土黄色| 天堂√8在线中文| 美女扒开内裤让男人捅视频| av线在线观看网站| 欧美亚洲 丝袜 人妻 在线| 黑人巨大精品欧美一区二区mp4| 精品久久久久久久毛片微露脸| 国产精品欧美亚洲77777| 国产欧美日韩一区二区三区在线| 1024视频免费在线观看| 国产精品乱码一区二三区的特点 | 色播在线永久视频| 午夜福利,免费看| 久久久久视频综合| 亚洲精品中文字幕在线视频| 天堂中文最新版在线下载| 啦啦啦在线免费观看视频4| avwww免费| 日本vs欧美在线观看视频| 999久久久精品免费观看国产| 亚洲黑人精品在线| 亚洲中文日韩欧美视频| 欧美黄色片欧美黄色片| 亚洲avbb在线观看| 国产成人av激情在线播放| 18在线观看网站| 精品第一国产精品| 欧美性长视频在线观看| 亚洲国产精品sss在线观看 | 涩涩av久久男人的天堂| 亚洲av第一区精品v没综合| 国产激情久久老熟女| 久久青草综合色| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦 在线观看视频| 久久久精品国产亚洲av高清涩受| 黑人操中国人逼视频| 国产精品av久久久久免费| www.熟女人妻精品国产| 无遮挡黄片免费观看| av福利片在线| 一边摸一边做爽爽视频免费| 精品国产一区二区三区久久久樱花| 人人妻人人澡人人爽人人夜夜| 亚洲 欧美一区二区三区| 一级a爱片免费观看的视频| 国产成人精品久久二区二区91| 大型av网站在线播放| 亚洲av熟女| 一进一出抽搐动态| 国产高清videossex| 中文字幕精品免费在线观看视频| 大型黄色视频在线免费观看| 黄片大片在线免费观看| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 久久人妻av系列| 日韩欧美一区视频在线观看| 777米奇影视久久| 国产精品影院久久| 中文字幕人妻熟女乱码| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费视频网站a站| aaaaa片日本免费| 美女视频免费永久观看网站| 最新的欧美精品一区二区| a级毛片黄视频| 久久国产精品大桥未久av| 国产精品电影一区二区三区 | 久久影院123| 丰满的人妻完整版| 熟女少妇亚洲综合色aaa.| 91成人精品电影| 亚洲人成伊人成综合网2020| 午夜视频精品福利| 国产av又大| 国产视频一区二区在线看| 日韩人妻精品一区2区三区| 咕卡用的链子| 宅男免费午夜| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 伦理电影免费视频| 国产真人三级小视频在线观看| 99热网站在线观看| 999久久久国产精品视频| 男人操女人黄网站| 一边摸一边做爽爽视频免费| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| 天天添夜夜摸| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 精品一区二区三区av网在线观看| 91麻豆精品激情在线观看国产 | 精品国内亚洲2022精品成人 | 日本一区二区免费在线视频| 超碰97精品在线观看| 黄片大片在线免费观看| 亚洲熟女毛片儿| 最近最新免费中文字幕在线| 欧美日韩精品网址| 久久99一区二区三区| 岛国在线观看网站| 国产精品久久视频播放| 999久久久国产精品视频| 久久久久国产一级毛片高清牌| 亚洲成人国产一区在线观看| 99久久99久久久精品蜜桃| 久久亚洲真实| 精品电影一区二区在线| 亚洲avbb在线观看| 国产欧美亚洲国产| 99国产精品一区二区三区| 嫁个100分男人电影在线观看| 黑人巨大精品欧美一区二区蜜桃| 法律面前人人平等表现在哪些方面| 亚洲av成人av| 欧美乱色亚洲激情| 青草久久国产| 人人妻,人人澡人人爽秒播| 日韩欧美在线二视频 | 超碰成人久久| 香蕉国产在线看| 久久国产乱子伦精品免费另类| 久久久久精品国产欧美久久久| 欧美日本中文国产一区发布| 老司机亚洲免费影院| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看 | 国产在视频线精品| 黄片大片在线免费观看| 村上凉子中文字幕在线| 9色porny在线观看| 国产亚洲精品久久久久5区| 亚洲熟女毛片儿| 操美女的视频在线观看| 18禁国产床啪视频网站| 搡老熟女国产l中国老女人| 这个男人来自地球电影免费观看| 好男人电影高清在线观看| 国产精品久久电影中文字幕 | 9色porny在线观看| 精品国产乱码久久久久久男人| 亚洲精品一卡2卡三卡4卡5卡| 国产极品粉嫩免费观看在线| av线在线观看网站| av天堂久久9| 午夜福利免费观看在线| 侵犯人妻中文字幕一二三四区| 国产精品秋霞免费鲁丝片| 午夜福利一区二区在线看| 精品高清国产在线一区| 麻豆av在线久日| 久久精品国产a三级三级三级| 欧洲精品卡2卡3卡4卡5卡区| 国产精品.久久久| 国产精品99久久99久久久不卡| 18禁黄网站禁片午夜丰满| 久久亚洲精品不卡| 热re99久久国产66热| 一级,二级,三级黄色视频| av线在线观看网站| 搡老乐熟女国产| 国内毛片毛片毛片毛片毛片| 在线观看www视频免费| 国产单亲对白刺激| 久久久久精品人妻al黑| 老熟妇乱子伦视频在线观看| bbb黄色大片| 另类亚洲欧美激情| 好看av亚洲va欧美ⅴa在| 国产精品成人在线| 久久国产精品男人的天堂亚洲| 脱女人内裤的视频| 99国产综合亚洲精品| 精品久久久精品久久久| 日本a在线网址| 国产极品粉嫩免费观看在线| 欧美成人午夜精品| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 久久精品91无色码中文字幕| 亚洲三区欧美一区| 国产精品偷伦视频观看了| 纯流量卡能插随身wifi吗| av片东京热男人的天堂| 国产一区有黄有色的免费视频| 国产精品久久电影中文字幕 | 午夜福利在线免费观看网站| 国产成+人综合+亚洲专区| 69av精品久久久久久| 国产亚洲精品久久久久久毛片 | av电影中文网址| 国产精品久久久久久人妻精品电影| 一级黄色大片毛片| 国产熟女午夜一区二区三区| 日本五十路高清| 亚洲精品乱久久久久久| 成年动漫av网址| 国产精品一区二区免费欧美| 国产精品.久久久| 在线观看免费高清a一片| 少妇裸体淫交视频免费看高清 | 高潮久久久久久久久久久不卡| 80岁老熟妇乱子伦牲交| 精品第一国产精品| 亚洲熟妇熟女久久| 精品久久久精品久久久| 欧美在线一区亚洲| 成人特级黄色片久久久久久久| 在线观看日韩欧美| 亚洲伊人色综图| 人妻久久中文字幕网| 国产一卡二卡三卡精品| 中文字幕av电影在线播放| 少妇裸体淫交视频免费看高清 | 免费在线观看亚洲国产| 国产一区在线观看成人免费| 91大片在线观看| 99re6热这里在线精品视频| 美女福利国产在线| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 嫁个100分男人电影在线观看| 男女午夜视频在线观看| 99精品在免费线老司机午夜| 免费看十八禁软件| 免费少妇av软件| 在线国产一区二区在线| 午夜精品在线福利| 久久久精品免费免费高清| 亚洲欧美一区二区三区久久| 久久久国产成人免费|