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

    A compact electron storage ring for lithographical applications

    2021-10-18 01:48:30SiQiShenDaZhangHuangZhenTangZhaoQingLeiZhang
    Nuclear Science and Techniques 2021年9期

    Si-Qi Shen ? Da-Zhang Huang ? Zhen-Tang Zhao, ? Qing-Lei Zhang

    Abstract The physical design for a novel low-energy compact-storage-ring-based extreme ultraviolet (EUV)light source was systemically studied.The design process considers the linear and nonlinear beam optics,including transverse matching and the optimization of the dynamic aperture,momentum aperture,and beam lifetime.With a total circumference of 36.7 m and a beam energy of 400 MeV,the storage ring can operate with an average beam current of up to 1 A.With the undulator as the radiator,this facility has the potential to emit EUV radiation at 13.5 nm with an average power exceeding 10 W within the bandwidth.In addition,the collective instabilities of the lattice at high beam current were analyzed;it was found that the typical instabilities which may occur in an electron storage ring can be reasonably controlled in our design.With the advantages of variable beam energy and current,this design exhibits great promise as a new candidate for various EUV lithographical applications requiring tunable radiation power.

    Keywords Storage ring · Extreme ultraviolet (EUV) ·EUV lithography (EUVL)

    1 Introduction

    With the recent rapid development of the microelectronics industry,the need for extreme ultraviolet lithography (EUVL) metrology has increased rapidly.Among many challenging tasks in EUVL metrology such as measuring the reflectivity of EUV multilayer coatings,characterizing the figure and finish of EUV mirror surfaces with high accuracy,and inspecting EUV mask blank substrates for pits and bumps and determining their flatness [1],the technical threshold of conventional EUV light sources such as laser-produced plasma sources is remarkably high [2],and these sources have several disadvantages as well.For example,to achieve higher levels of in-band EUV brightness,time-multiplexing of several sources is required,which increases both the cost and the complexity [3].However,the recent development of electron-storage-ringbased light sources suggests a new approach to EUVL.New mechanisms such as stable-state microbunching and angular dispersion modulation have been studied [4–6],and machines have been designed and built,for example,the compact storage ring for actinic mask inspection(COSAMI) [7] and the Metrology Light Source [8].

    However,because the EUV light sources for lithography are ultimately intended for use in industrial mass production,the performance-to-cost ratio is an essential consideration.For example,because the electron beam emittance is not highly important for photolithography,it can be less important in our design.Consequently,the lattice complexity and the cost can be dramatically decreased compared with those of the schemes mentioned above.Moreover,a higher emittance enables a higher bunch charge,which supports a higher radiation power.Therefore,a design study of a compact storage ring lattice with reduced complexity and relatively low cost for EUVL applications has been initiated.

    In this study,our main goal was to systemically study the feasibility of a storage ring that can operate stably with an average beam current of up to 1 A and deliver EUV radiation with an average power suitable for lithographical applications such as wafer mask inspection.In addition,as mentioned above,to reduce the complexity of the machine and control the overall cost,little effort was made to decrease the emittance as long as the emittance is acceptable for the ring.The storage ring was optimized to obtain sufficient dynamic aperture for injection and ample momentum aperture to ensure the required beam lifetime.The working points were also analyzed to avoid the resonance lines as far as possible,and the bunch length was increased using the third harmonic cavity to increase the beam lifetime.

    The remainder of this paper is organized as follows.In Sect.2,we describe the lattice design,including the design and optimization of the linear and nonlinear beam optics.The lattice layout is presented,and the optimization of the dynamic and momentum apertures is discussed in detail.Frequency map analysis (FMA) and the evolution of the momentum aperture along the ring are also discussed.We estimate the beam lifetime when the third harmonic cavity is included.The injection scheme that will probably be used is also introduced.At the end of this section,the threshold of the typical collective instabilities is analyzed to estimate the stable operation of the machine,and more work is suggested to ensure good control of various collective effects.In Sect.3,the undulator parameters are calculated to obtain the required EUV light,and the key indexes of the EUV light are also presented.The conditions for achieving it are analyzed,and concluding remarks are given in the last section.

    2 Lattice design and optimization of the storage ring

    2.1 Linear optics and dynamic aperture

    The layout of the storage ring is illustrated schematically in Fig.1.The racetrack-like lattice consists of four double-bend achromat cells,which minimize the lattice complexity and meet the basic requirements of the storage ring for acceptable beta functions and dispersion,and two long straight sections each 8 m in length.These sections are reserved for insertion devices such as the undulator,also called a radiator,which will generate radiation.The total circumference of the ring is 36.7 m,and it will operate at the 61st harmonic of the 500 MHz fundamental radio frequency(RF).A storage ring energy of 400 MeV is chosen to reduce the required RF power and to obtain the desired EUV wavelength.The natural emittance is approximately 44 nm·rad in the horizontal direction,and the average beam current is 1 A so that photons with the necessary flux and power can be produced.The basic parameters of the ring are shown in Table 1.

    The integer of the fundamental RF wavelength was chosen as the total length of the ring.There are two bending magnets,four quadrupoles,and four sextupole magnets in each cell.The bending magnets are all 0.816 m in length and have a bending angle of 45°.The phase advance of the reference particle is adjusted by changing the strengths of the quadrupoles to optimize the working point,i.e.,by horizontal and vertical tuning,to avoid various types of resonance such as external resonance,parametric resonance,and sum,and difference resonance [9].The goal is to ensure the stability of the closed orbit of the particle and to obtain horizontal and vertical orbits that are closed to each other in the straight sections.After the linear optics were optimized,the main parameters of the storage ring were determined.The optimization was performed using the NSGA-II program.The typical optical functions of the lattice are shown in Fig.2.

    Fig.2 (Color online) Typical optical functions of the lattice.The blue curve is the horizontal beta function,the red curve is the vertical beta function,and the green curve is the horizontal dispersion function

    The sextupole magnets are used to correct the chromaticity and to optimize the dynamic and momentum apertures.There are four sextupoles in each cell.Two are used to obtain positive chromaticity;they are called the chromatic sextupoles and are located between the two bending magnets,where the dispersion is large.The other two sextupole magnets,which are called harmonic sextupoles,are employed to optimize the dynamic aperture,which should be large enough for injection,and to optimize the momentum aperture for a longer beam lifetime.Thesetwo magnets are placed at the ends of the bending magnets near the long straight section,where the dispersion is zero.

    The chromaticity is defined as the derivative of the transverse tune over the relative momentum spread of a particle.It arises because the quadrupoles exert different focusing forces on particles of different energies.Therefore,there is always a natural chromaticity (Eq.1) provided by the quadrupole magnets in the storage ring.The natural chromaticity is usually negative because the focusing force on higher-energy particles is weaker [9].

    where β is the beta function,Kis the strength of the quadrupole,andxandyrepresent the horizontal and vertical directions,respectively.Negative chromaticity will cause the difference Δν between particles with different energies to increase,and thus νx,ywill be close to the resonance value.To avoid negative chromaticity,the chromatic sextupoles are placed at locations where the dispersion is large,and the chromaticity can be reduced to zero by using the higher-order focusing term of the sextupole magnets (Eq.2) [9].

    whereSrepresents the sextupole strength,andDis the dispersion.During optimization,it is difficult to obtain zero chromaticity in the horizontal and vertical directions simultaneously.Therefore,the chromaticity is usually optimized to positive values slightly larger than zero to ensure that particles with higher energy receive a stronger focusing force so that Δν between particles is controllable.

    The dynamic aperture of the beam in a storage ring is usually defined as the phase space area in which particles cannot be lost after multiple revolutions.The area of the transverse cross-section of the beam is usually employed to define the dynamic aperture.Because of the tune spread of the beam and the nonlinear focusing force,the particles falling on the resonance lines will be lost.This determines the maximum transverse area in which beam particles can survive,which is the dynamic aperture.Because the focusing force is nonlinear,sextupole families are used for optimization.

    The chromaticity should remain constant during the maximization of the dynamic aperture.Thus,the harmonic sextupole magnets,used for optimization,are usually placed at locations where the dispersion is close to zero,and the dynamic aperture is optimized by adjusting the strengths of these magnets.For further optimization,octupole magnets could also be added in the future to eliminate higher-order terms.

    In the optimization process,we used the beam tracking code Elegant [10],and we used its built-in simplex algorithm for multi-objective optimization.The goal was to optimize the driving terms related to the chromaticity,tune,and particle momentum of the bunch (first and second order)[11].The optimized variables are the strengths of the harmonic sextupoles [12].The convergence of the results indicates that the optimal condition was reached.Moreover,radiation damping and quantum excitation were also taken into account to obtain more realistic results.

    Figure 3 shows the FMA diagram of the optimized storage ring.Colors represent the logarithm of the tune difference (Δν)between an arbitrary beam particle and the reference particle.The region,where Δν is relatively uniform,is considered the dynamic aperture of the beam.The figure shows that the dynamic aperture of the beam reached several centimeters in diameter in the transverse direction after optimization,which meets the beam injection requirements of the storage ring.Figure 4 shows the working points,which were obtained by taking the beam energy spread into account.The data points in the figure represent the fractional tunes in both the horizontal and vertical directions.The resonance lines up to the fourthorder resonance are drawn.Most of the working points within the dynamic aperture (blue) are not on resonance lines.

    Fig.3 (Color online) FMA diagram of the ring

    Fig.4 (Color online) Working points of the ring considering the beam energy spread

    Figure 5 shows the dynamic aperture,with error bars,before and after orbit correction.Each cell was divided into two parts on separate girders.The installation errors were set to 50 μm in both the horizontal and vertical directions within the girder and 80 μm between the girders,with 3 σ truncation.For orbit correction,four correctors and three beam position monitors (BPMs) were added to each cell,and all the correctors were combined with the sextupoles.Accelerator toolbox (AT) [13] was employed to track the beam for 5000 turns.The result shows that the average dynamic aperture of the beam can reach ± 10 mm in the horizontal direction and ± 7 mm in the vertical direction.In addition,the minimal dynamic aperture after orbit correction,considering the error bars,is approximately ± 7 and ± 5 mm in the horizontal and vertical directions,respectively.Therefore,the dynamic aperture after orbit correction meets the injection requirements.

    Fig.5 (Color online)Dynamic aperture,with error bars,before(red)and after (blue) orbit correction

    In the last part of this section,we briefly introduce the injection scheme that will probably be used for the main ring.First,a booster ring is needed to accumulate and inject the beam in top-up mode to achieve stable injection with a current variation of less than 0.1%.The booster will be placed immediately above or below the main storage ring to save space.In the scheme,the beam is extracted by a kicker magnet with a rapid rise time and immediately deflected toward the storage ring by a septum magnet[14,15].A slow orbit bump is employed to generate a displacement of <10 mm horizontally at the septum position.A linac injector with a photocathode electron gun will also be installed near the booster ring to supply the electron beam.

    2.2 Momentum aperture,higher harmonic cavity,and beam lifetime

    In a low-emittance storage ring,the Touschek lifetime dominates the beam lifetime [16].Touschek scattering causes the loss of electrons and determines the lifetime.Transverse scattering results in large energy kicks to particles,which undergo betatron oscillation and are lost rapidly on the physical aperture.The loss rate is proportional to the particle density in the bunch and inversely proportional to the momentum acceptance,which is also called the momentum aperture.The rate of beam loss due to the Touschek effect takes the form [17]

    whereTstands for the Touschek lifetime,and 〈...〉 represents the average throughout the entire ring.

    Similar to the dynamic aperture,the momentum aperture represents the momentum/energy tolerance of the beam.As above-mentioned,when the energy spread is too large,particle loss can occur.On the other hand,a larger momentum aperture corresponds to a longer Touschek lifetime,which means that the upper limit of the momentum aperture determines the Touschek lifetime.Touschek scattering is the large-angle scattering from the transverse to the longitudinal direction and excites the large-amplitude betatron oscillation of scattered particles,which will cause the beam to be lost if the oscillation exceeds the physical aperture of the machine,limiting the beam lifetime.The momentum aperture of the beam is the key factor that determines the Touschek lifetime.It can be estimated by numerical beam tracking.A positive momentum kick[δm(s)] is given to a particle in the bunch at a given position(s),causing it to undergo betatron oscillation.The maximum kick that the particle can withstand without loss after numerous turns is the positive momentum aperture.The negative momentum aperture can be obtained by the same method using a negative kick.

    We used Elegant to calculate and optimize the momentum aperture of the beam [18].The optimization method and algorithm are the same as those used to compute the dynamic aperture.In the optimization process,radiation damping by both the bending magnets and undulators was included,along with the RF cavities.In addition,the physical aperture was taken as ± 2.5 cm horizontally and ± 4 mm vertically in the straight section considering the undulator gap.To ensure that the damping time is sufficiently long for the beam to reach the stable state,we calculated the momentum aperture after numerous turns.We found that when the beam was tracked for more than 10,000 turns,the result remained essentially unchanged,as shown in Fig.6.We expect that additional work will further maximize the momentum aperture.

    Fig.6 (Color online) Variation of the momentum aperture throughout the lattice after 10,000 turns

    According to Eq.(3),the Touschek lifetime is determined by the number of particles in the bunch,beam energy,beta function,transverse emittance,and bunch length.Among these parameters,the transverse emittance and bunch length are most often tweaked to optimize the lifetime.The former can be adjusted by changing the horizontal/vertical coupling ratio,and the latter can be increased by inserting a higher-order harmonic cavity into the lattice.

    A cavity with a resonant frequency that matches the harmonics of the fundamental RF in a storage ring is called a higher-order harmonic cavity.Its main function is to adjust the beam energy to control the beam dynamics by adjusting either the energy of electrons at different locations in the beam or that of different bunches in the bunch train.The higher-order harmonic cavity interacts with the beam in different ways to stretch or compress the bunch,and the harmonics of the cavity determine the cavity size[19].In the storage ring,we designed the fundamental RF in the ring is approximately 500 MHz,and we chose a frequency of 1.5 GHz for the high-order harmonic cavity;thus,the harmonic number is 3.

    The bunch length in the time domain can be calculated as

    wherenis the harmonic number,αcis the momentum compaction factor,ωRFis the frequency of the fundamental RF cavity,σEis the energy spread,ωsis the angular frequency of synchrotron oscillation without the harmonic cavity,E0is the energy of the storage ring,cis the speed of light,and=3.6256.

    To illustrate the problem,we substitute the cavity voltage into the longitudinal motion equation of the beam to obtain the potential well as a function of RF phase when the harmonic cavity is present:

    For comparison,the potential well without the harmonic cavity reads

    where φ is the RF phase,and φsis the synchrotron phase.By comparing Eqs.(5) and (6),we can see that the potential well is lengthened,which provides more space to hold a longer bunch,as shown in Fig.7.

    Fig.7 (Color online)Potential well with(red)and without(blue)the third harmonic cavity.Note that the horizontal axis has units of the ratio of the RF phase over π

    Therefore,to obtain a sufficient Touschek lifetime for a high beam current,we introduce a third harmonic cavity into the lattice to increase the bunch length from 3.5 to 14.2 mm rms.Considering the natural coupling in the ring introduced by misalignments,magnet tilts,dispersion leaks,etc.,we take the coupling value to be 10%.From Fig.4,we can see that the minimum momentum aperture(-0.5%-0.3%) appears at the locations of the dipoles,and the maximum aperture(~±1.5%)almost fills the rest of the lattice,which is more than 80% of the total length.As a result,using Elegant calculations,we obtained average Touschek lifetimes of approximately 1.3 h at an average beam current of 1 A weighted over the length of regions with different momentum apertures with a 10%vertical/horizontal coupling ratio and of approximately 40 min at 2 A.These values are acceptable for this high beam current in top-up mode.Note that all the RF buckets will be filled with electrons;thus,the charge of a single bunch at 1 A is 2 nC at a fundamental RF of 500 MHz.The possibility of ion instability is discussed in Sect.2.3.

    2.3 Estimation of collective effects

    In this subsection,we roughly estimate the possible collective effects in our ring to ensure that the machine can operate stably.As noted above,the design goal is an average beam current of up to 1 A,which is quite high.Thus,the collective effects could play an important role in our machine.

    The sources of collective effects may include electron clouds,ions,and impedance elements in the ring.The impedance elements in the ring cause various instabilities,which often limit high-current operation,and interactions between an electron cloud and ions usually introduce beam-ion instabilities,including ion traps and fast ion instability (FII) at high beam currents.In this section,we discuss the typical instabilities that may appear in our ring.

    When a charged particle moves in a beam pipe,it will excite an electromagnetic field if the shape of the pipe changes or resistance is present in the pipe wall.Because this field is behind the source particle,it is called the wakefield.In many cases,the wakefield takes some time to decay;therefore,it affects the trailing particles and may introduce instabilities [20].In a real machine,we would observe the spectra of the beam.In addition,because the beam dynamics is usually analyzed in the frequency domain,we performed a Fourier transform of the wakefield to obtain the impedance,for example,s.

    wheresis the longitudinal coordinates of the witness particle.

    The narrowband impedance usually results in coupledbunch instability.The introduction of broadband impedance generally excites microwave instability in both the transverse and longitudinal directions.In the following,we outline both the narrowband and broadband impedance and their effects on the collective instabilities in our ring.

    The resistive-wall (RW) impedance can be both narrowband and broadband [16].In addition,because the resistive beam pipe occupies most of the length of the ring,the RW impedance plays a major role in our design.As an example,we estimated the longitudinal coupled-bunch instability by comparing the impedance threshold,as shown in Eq.(8).

    whereE0is the beam energy,τSis the longitudinal radiation damping time,I0is the average beam current,Qsis the synchrotron tune,ω is the frequency of the wakefield,η is the momentum compaction factor,and the longitudinal RW impedance is given by Eq.(9) [21].

    where/nis the effective longitudinal impedance,Z0≈377Ω is the vacuum impedance,ρcis the material resistance,bis the pipe radius,μ is the material permeability,and μ0is the free space permeability.Figure 8 compares the effective longitudinal RW impedance and the threshold of the longitudinal coupled-bunch instability as a function of frequency,which was calculated using the parameters in Table 1 with a 1 A average beam current.Copper was selected as the material of the vacuum pipe,which has an inner radius of 25 mm and a wall thickness of 1.5 mm.

    In Fig.8,we can see that the effective longitudinal RW impedance is below the threshold up to 20 GHz,which covers almost all the possible bandwidth of the instability.However,the longitudinal narrowband geometric impedance introduced by the RF cavities,kickers,BPMs,etc.,produces additional small peaks on the RW impedance curve in Fig.7.Among them,the higher-order modes of the RF cavity can be the most harmful sources of coupledbunch instability.Therefore,the cavities in our rings will be carefully designed.Nonetheless,although more work is necessary to investigate the details,on the basis of past experience and a comparison with the RW impedance given above,we are confident we can design and build devices that satisfy the threshold.

    Fig.8 (Color online) Comparison of impedance threshold of longitudinal coupled-bunch instability(red)and effective longitudinal RW impedance (blue)

    Next,as is well known,the RW instability usually refers to the transverse coupled-bunch instability arising from the long-range RW wakefields.Therefore,we also estimated the transverse coupled-bunch instability by calculating the current threshold at zero chromaticity [16],which is the worst case of the instability,and it reads

    where τ⊥is the transverse radiation damping time,β⊥is the average of the beta function along the ring,and ω0is the angular revolution frequency of beam particles.ΔQβis the fractional part of the betatron tune,and the other variables are defined above.The result obtained using the beam parameters in Table 1 shows that the average current threshold is approximately 80 mA.Although the nominal design current of our machine is 1 A,which gives a ratio of 12 between the two currents,the instability is still within the range that a standard transverse feedback system can currently handle [13,20].

    The longitudinal broadband impedance can introduce longitudinal microwave instability.When this instability occurs,the longitudinal current distribution of the beam is impaired,and the energy spread increases.As a result,the beam quality deteriorates,and beam loss may occur.We estimated the impedance threshold of a single bunch for the longitudinal microwave instability in a storage ring on the basis of the Keil–Schnell criterion for a coasting beam[9,22],which was extended to a bunched beam by Boussard [16].

    whereRis the average radius of the storage ring,E0is the energy of the storage ring,σE0is the energy spread,σl0is the bunch length,η is the phase slip factor,Ibis the bunch current,and|z/n|effis the effective longitudinal impedance.

    Our machine is designed for operation at an average beam current of 1 A with a harmonic number of 61.Thus,if we choose a single-bunch current of 16.4 mA,the total effective broadband impedance will be less than 0.39 Ω,according to Eq.(11).For comparison,the total effective broadband impedance of the SSRF is approximately 0.17 Ω,and those of Soleil and BEPC are approximately 0.21 and 0.23 Ω,respectively.Therefore,we believe that microwave instability will not be a problem in our machine.

    In addition to the longitudinal microwave instability,the transverse broadband impedance may also introduce another single-bunch instability called the transverse mode coupling instability(TMCI)or the fast head-tail instability,which is also detrimental to the machine.TMCI occurs when the frequencies of two neighboring head-tail modes approach each other because of detuning with increasing current during accumulation [21].To estimate the effect,we first compute the transverse broadband impedance from the longitudinal broadband impedance and the Panofsky–Wenzel theorem.

    Like that in the longitudinal direction,the single-bunch current threshold of the TMCI can be derived as [23]

    where β⊥represents the average beta function over the machine.In our design,the radius of the beam pipe is 25 mm;thus,the single-bunch current threshold of the transverse microwave instability is 82.6 mA,which is much larger than the longitudinal single-bunch current threshold of our machine.Therefore,the transverse microwave instability and TMCI will not be a problem either.Finally,it is necessary to note that the chromaticity is conventionally adjusted to a positive value to avoid the TMCI,whereas positive chromaticity usually results in certain higher-order modes of another transverse instability called the head-tail instability,and the threshold of the first unstable mode of the instability is even lower than that of the TMCI[16].The patterns of the instability are different at low and high positive chromaticity and have been studied at ESRF [24].Thus,caution should be used in designing the chromaticity in our machine.

    To achieve an average beam current of up to 1 A,all the RF buckets throughout the ring will be filled with electrons.Consequently,another collective effect called the ion instability,which includes ion trapping and the FII,had to be estimated because of the limited spacing between beam bunches.The instability arises from interactions between the electron beam and the ionized residual gas,and it induces transverse center-of-mass oscillation of the beam.The growth time of the ion trap instability is generally long,and it can be alleviated by introducing bucket gaps,which are usually approximately a few percent of the ring circumference,behind the bucket train during filling.However,in modern storage rings,the single-pass FII is more significant owing to the small emittance and multibunch operation.The asymptotic growth rate of the FII derived by Raubenheimer and Zimmermann for a flat beam is given by [25]

    whereNedenotes the number of electrons in a bunch,nbis the number of bunches in the train,γ is the relativistic factor,pgasis the residual gas pressure in Torr,andreandrpare the classical electron and proton radii,respectively.Ais the atomic mass number,andLsepis the bunch spacing.Because the average beam current is 1 A in our design,the charge of a single bunch is thus 2 nC for a fundamental RF of 500 MHz,and all the RF buckets are filled at a harmonic number of 61.CO is chosen as the residual gas in our calculation because its ionization cross section is 6 times that of H2,and the gas pressure is set to 10-9Torr,which is not hard to reach in storage rings[26,27].In addition to the parameters in Table 1,we have a rising time of 0.63 ms for the instability,which is within the bandwidth of the feedback system [27].Therefore,we conclude that the FII in our machine can be well-controlled with a relatively high but achievable vacuum level.

    Furthermore,intrabeam scattering and Touschek scattering are also important collective effects in high-current storage rings.As mentioned in the preceding section,the growth rates of both scattering effects are inversely proportional to the transverse emittance of the beam.Because the nominal transverse emittance in our design is one order larger than that of third-generation light sources (or even larger),the above two effects can also be controlled.

    As in other storage rings,the collective instabilities in our machine can be addressed by both chromaticity shifting and transverse feedback.However,as above-mentioned,because shifting the chromaticity to a positive value may also excite some higher-order modes of the head-tail instability,it alone is not sufficient for good control of the instabilities.Therefore,as in SOLEIL,transverse bunchby-bunch feedback will be combined with chromaticity[28] to relieve the collective instabilities in our machine.

    Finally,we must point out again that we use copper as the beam pipe material,and the inner radius of the vacuum chamber is 25 mm.According to the first-order geometric strengthK1of the quadrupoles and the second-order strengthK2of the sextupoles,the maximum integrated field for the strongest quadrupole is 1.36 T/m,and the maximum integrated field for the strongest sextupole is 163.3 T/m2,which are not hard to achieve for the 400 MeV electron beam in the vacuum pipe as described above.Moreover,more detailed work is needed to provide more precise estimations of the collective effects in our design,such as the coupled-bunch instability and beam loading.These effects are topics for further investigation.

    In conclusion,we believe that most of the typical collective instabilities that may arise in our machine can be controlled using feasible methods.However,as mentioned in Sect.1,the emittance of our ring is not of great importance;therefore,the tolerance of the collective instabilities of our machine could be slightly larger than those of other storage ring light sources,which can be considered to give our design an advantage over others.

    3 Considerations regarding EUV radiation

    The typical wavelength currently used for EUVL is 13.5 nm.As in other machines,an undulator is employed to provide the required radiation.To evaluate the performance of our design,the SPECTRA simulation code[29–31]was used.In the simulation,a photon beam with an emittance of 40 nm in the horizontal direction and 4 nm in the vertical direction is obtained by the convolution of the electron beam and diffraction phase spaces.Strong magnets are required to obtain a maximum field exceeding 0.3 T with a period length of 15 mm.The simulation result shows that the undulator parameter is approximately 0.47.To obtain the highest photon flux,the entire central cone of the photon emission must be focused on the sample,for example,a microchip.

    The parameters of the undulator as calculated by SPECTRA are shown in Table 2.According to the design parameters of the storage ring,the required 13.5 nm EUV radiation can be obtained using the undulator.According to the data in Table 2,we can see that for an average current of 1 A and a 7 m undulator,the total radiation power within 5% bandwidth (FWHM) at 13.5 nm is 4.1 W and can also exceed 10 W with a >2 A beam current.Note,however,that in these cases,the undulator gap is approximately 10 mm,and the wakefield induced by the undulator will be somewhat strong.A superconducting undulator[32]may be required to avoid this possibility,and the overall cost may rise.Alternatively,an economical method is to lower the average beam current to several hundred milliamps to reduce the wakefield dramatically and afford much easier control.Although there would be a corresponding reduction in the radiation power to a few watts,we can still use the machine for mask inspection,as at the COSAMI.On the other hand,because the radiation power is proportional to the undulator length,the required radiation power can also be obtained by decreasing the beam current and increasing the undulator length,which is easier to realize than a higher beam current in the real machine without excessive cost.Moreover,a long variable-gap undulator can be made by combining several short sections with a dedicated control and a diagnostic system,whereas there is no such problem for the fixed gap undulator.

    Table 2 Basic design parameters of undulator and key indexes of EUV light

    Furthermore,in our design,the bunch charge,duty factor of the RF bucket,and undulator gap can all be varied.Therefore,a wide radiation power range can be obtained to expand the use of our design under different conditions and at different costs.This characterizes our design as a versatile tool for various EUVL applications.

    4 Conclusion

    A systematic physical design study of a compact 400 MeV storage-ring-based EUV light source for industrial lithography was conducted.Owing to its compactness,reduced complexity,and relatively low cost,the storage ring can provide an average beam current of up to 1 A with a dynamic aperture of approximately ± 10 mm.If the third harmonic cavity is used,an average Touschek lifetime of 1.3 h can also be achieved.Moreover,the collective effects were also estimated at high average current.It was found that the typical collective instabilities can be controlled and relieved by implementing appropriate materials and methods,such as a copper beam pipe and modern feedback technology.With the designed undulator,the machine can radiate at 13.5 nm in the EUV regime with an average radiation power within the bandwidth of several watts to more than 10 W.Although the radiation power in the bandwidth appears to be inadequate for directphotolithography,our design can still be used for multiple lithographic applications,such as wafer mask inspection,EUV mirror characterization,and reflectivity measurements of EUV multilayer coatings.In conclusion,our study demonstrates the feasibility of a less complex,less costly compact-storage-ring-based EUV light source with tunable current and radiation power.As an alternative to the conventional EUVL technology,it provides a new option for multiple applications in EUVL and a valuable reference for the design of similar facilities in the future.

    AcknowledgementsThe authors would like to thank Dr.Bo-Cheng Jiang for useful discussions on storage ring physics,Dr.Zhen Wang for assistance with SPECTRA,Dr.Chao Feng for advice on the undulator,Dr.Xiao-Xia Huang for assistance with the impedance calculation,Dr.Yi-Yong Liu for the vacuum pipe calculation,and Dr.Wei Zhang for advice on the magnets.

    Author contributionsAll authors contributed to the study conception and design.Material preparation,data collection,and analysis were performed by Si-Qi Shen,Da-Zhang Huang,and Zhen-Tang Zhao.The first draft of the manuscript was written by Si-Qi Shen,and all authors commented on subsequent versions of the manuscript.All authors read and approved the final manuscript.

    国产成人a区在线观看| 午夜亚洲福利在线播放| 黄色一级大片看看| 日日摸夜夜添夜夜添av毛片| 精品久久久久久久久久久久久| 日韩 亚洲 欧美在线| 99热6这里只有精品| 免费在线观看成人毛片| 午夜老司机福利剧场| 日产精品乱码卡一卡2卡三| 中文字幕久久专区| 日韩视频在线欧美| 久久久久久久亚洲中文字幕| 最近最新中文字幕大全电影3| 亚洲在线观看片| 国产日本99.免费观看| 乱人视频在线观看| 久久久成人免费电影| 日日啪夜夜撸| 日本成人三级电影网站| 日韩在线高清观看一区二区三区| 日本三级黄在线观看| 九九热线精品视视频播放| 国产精品av视频在线免费观看| 亚洲欧美精品专区久久| 老女人水多毛片| 在线免费观看不下载黄p国产| 一个人观看的视频www高清免费观看| 免费不卡的大黄色大毛片视频在线观看 | 插阴视频在线观看视频| 国产精品日韩av在线免费观看| 亚洲自拍偷在线| 精品不卡国产一区二区三区| 免费一级毛片在线播放高清视频| 哪里可以看免费的av片| 日日摸夜夜添夜夜添av毛片| 午夜久久久久精精品| 免费在线观看成人毛片| 天天躁日日操中文字幕| 国产成人福利小说| www日本黄色视频网| 久久久久久九九精品二区国产| 赤兔流量卡办理| 国产探花在线观看一区二区| 国产精品野战在线观看| 日韩成人伦理影院| 97超碰精品成人国产| 成人二区视频| 好男人在线观看高清免费视频| 18禁在线播放成人免费| 久久婷婷人人爽人人干人人爱| 男女做爰动态图高潮gif福利片| 中文精品一卡2卡3卡4更新| 国产av不卡久久| 我要搜黄色片| 一个人免费在线观看电影| 久久久精品94久久精品| 日韩av在线大香蕉| 给我免费播放毛片高清在线观看| 国产精品一二三区在线看| 精品99又大又爽又粗少妇毛片| 国产又黄又爽又无遮挡在线| 午夜视频国产福利| 99热这里只有是精品50| 麻豆国产av国片精品| 天美传媒精品一区二区| 美女xxoo啪啪120秒动态图| 两个人视频免费观看高清| 欧美成人一区二区免费高清观看| 国产老妇伦熟女老妇高清| 亚洲国产色片| 成人午夜高清在线视频| 亚洲国产欧美人成| 免费电影在线观看免费观看| 夜夜看夜夜爽夜夜摸| 亚洲电影在线观看av| 精品人妻偷拍中文字幕| 久久99精品国语久久久| 国内揄拍国产精品人妻在线| 99国产极品粉嫩在线观看| 久久欧美精品欧美久久欧美| 自拍偷自拍亚洲精品老妇| 少妇人妻一区二区三区视频| 亚洲欧美日韩东京热| 久久亚洲国产成人精品v| 特大巨黑吊av在线直播| 亚洲人成网站高清观看| 日韩成人伦理影院| 中国美白少妇内射xxxbb| 99久国产av精品| 亚洲av成人精品一区久久| 亚洲人成网站在线观看播放| 国产精品电影一区二区三区| 日本黄色片子视频| 久99久视频精品免费| 日韩高清综合在线| 精品免费久久久久久久清纯| 97在线视频观看| 丰满人妻一区二区三区视频av| 最新中文字幕久久久久| 蜜臀久久99精品久久宅男| 日本黄大片高清| 高清午夜精品一区二区三区 | 精华霜和精华液先用哪个| 日韩人妻高清精品专区| 亚洲欧美精品自产自拍| 网址你懂的国产日韩在线| 亚洲,欧美,日韩| 在线观看一区二区三区| 特大巨黑吊av在线直播| 天堂影院成人在线观看| 大香蕉久久网| 夫妻性生交免费视频一级片| 国产精品野战在线观看| 日本一本二区三区精品| 国产日韩欧美在线精品| 国产在线精品亚洲第一网站| 在线国产一区二区在线| 你懂的网址亚洲精品在线观看 | 亚洲第一电影网av| 成人一区二区视频在线观看| 看十八女毛片水多多多| 国产亚洲精品久久久久久毛片| 欧美日韩在线观看h| 国产精品免费一区二区三区在线| 成人国产麻豆网| 色噜噜av男人的天堂激情| 精品一区二区三区人妻视频| 在线观看免费视频日本深夜| 亚洲人成网站在线播放欧美日韩| 2022亚洲国产成人精品| 国产精品,欧美在线| 青青草视频在线视频观看| 在线免费十八禁| 欧美成人免费av一区二区三区| 精品久久久久久成人av| 国产v大片淫在线免费观看| 午夜福利成人在线免费观看| 中文字幕免费在线视频6| 校园春色视频在线观看| 国产精品电影一区二区三区| 国产黄片美女视频| 国内揄拍国产精品人妻在线| 国产一区二区在线av高清观看| 欧美zozozo另类| 少妇人妻精品综合一区二区 | 能在线免费观看的黄片| 国产久久久一区二区三区| 久久久欧美国产精品| 国产伦理片在线播放av一区 | 日日啪夜夜撸| 亚洲av免费在线观看| 在线观看美女被高潮喷水网站| 国产真实乱freesex| 亚洲av不卡在线观看| 99热这里只有精品一区| 亚洲成a人片在线一区二区| 麻豆精品久久久久久蜜桃| 国产精品麻豆人妻色哟哟久久 | а√天堂www在线а√下载| 狠狠狠狠99中文字幕| 狠狠狠狠99中文字幕| 国产精品久久久久久av不卡| 能在线免费看毛片的网站| 一本一本综合久久| 少妇人妻精品综合一区二区 | 国产午夜精品一二区理论片| 国产成人一区二区在线| 波野结衣二区三区在线| 日本五十路高清| 国产亚洲91精品色在线| 国产精品免费一区二区三区在线| 美女黄网站色视频| 精品午夜福利在线看| 波多野结衣高清作品| 麻豆国产97在线/欧美| 成人特级av手机在线观看| 亚洲最大成人av| 免费看美女性在线毛片视频| 亚洲国产精品成人综合色| 久久久国产成人免费| 久久精品国产亚洲av天美| 中文资源天堂在线| 中文字幕av成人在线电影| 亚洲成人中文字幕在线播放| 自拍偷自拍亚洲精品老妇| 国产三级在线视频| 午夜a级毛片| 97在线视频观看| 免费观看的影片在线观看| 久久精品国产鲁丝片午夜精品| 国产探花极品一区二区| 身体一侧抽搐| 国产av在哪里看| 亚洲18禁久久av| 亚洲高清免费不卡视频| 麻豆乱淫一区二区| 精品久久久久久久人妻蜜臀av| 在线天堂最新版资源| 亚洲精品自拍成人| 精品不卡国产一区二区三区| 亚洲av二区三区四区| 亚洲激情五月婷婷啪啪| 亚洲av.av天堂| 午夜亚洲福利在线播放| 女人被狂操c到高潮| 国产激情偷乱视频一区二区| 久久精品91蜜桃| 免费观看在线日韩| av在线蜜桃| 亚洲高清免费不卡视频| 久久综合国产亚洲精品| 天天躁日日操中文字幕| 夜夜夜夜夜久久久久| 国产在线男女| 国产女主播在线喷水免费视频网站 | 久久中文看片网| 岛国在线免费视频观看| 非洲黑人性xxxx精品又粗又长| 成人欧美大片| 成年版毛片免费区| 成人毛片a级毛片在线播放| 亚洲精品色激情综合| 美女cb高潮喷水在线观看| 国产久久久一区二区三区| 久久久久国产网址| 国产极品天堂在线| 99久久精品热视频| 久久久久久久午夜电影| 能在线免费观看的黄片| 成人亚洲精品av一区二区| 免费观看精品视频网站| 国产伦理片在线播放av一区 | 日韩中字成人| 最近视频中文字幕2019在线8| 老司机影院成人| 欧美一级a爱片免费观看看| 黄色配什么色好看| 成人毛片a级毛片在线播放| av国产免费在线观看| 天堂√8在线中文| 一区二区三区四区激情视频 | 国产亚洲精品久久久久久毛片| 国产精品福利在线免费观看| 久久精品国产清高在天天线| 嫩草影院新地址| 亚洲精品456在线播放app| 在线国产一区二区在线| 国产黄片美女视频| 国产免费一级a男人的天堂| 亚洲成人精品中文字幕电影| 久久久成人免费电影| 亚洲无线在线观看| 卡戴珊不雅视频在线播放| 男人和女人高潮做爰伦理| 国产亚洲欧美98| 精品人妻一区二区三区麻豆| 成人永久免费在线观看视频| 久久精品91蜜桃| 国产 一区精品| 天美传媒精品一区二区| 99热这里只有精品一区| 欧美最新免费一区二区三区| 3wmmmm亚洲av在线观看| 色播亚洲综合网| 中出人妻视频一区二区| 精品久久久久久久久久免费视频| 欧美日韩国产亚洲二区| 午夜福利高清视频| 日韩欧美三级三区| 日韩高清综合在线| 一本精品99久久精品77| 两个人的视频大全免费| 久久久久久国产a免费观看| 可以在线观看的亚洲视频| 直男gayav资源| 熟妇人妻久久中文字幕3abv| 自拍偷自拍亚洲精品老妇| 在线观看美女被高潮喷水网站| 亚洲av电影不卡..在线观看| 三级经典国产精品| 91麻豆精品激情在线观看国产| 久久久国产成人免费| 99久久精品一区二区三区| av天堂中文字幕网| 日韩精品青青久久久久久| 黄片wwwwww| 成人特级av手机在线观看| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 国产极品精品免费视频能看的| 久久精品影院6| 亚洲成人久久爱视频| 午夜免费男女啪啪视频观看| 人体艺术视频欧美日本| 亚洲成人精品中文字幕电影| 丰满人妻一区二区三区视频av| 一区二区三区四区激情视频 | 成人一区二区视频在线观看| 91午夜精品亚洲一区二区三区| 久久久精品欧美日韩精品| 干丝袜人妻中文字幕| 国内精品一区二区在线观看| 成人国产麻豆网| 性色avwww在线观看| 精品久久久久久久久久久久久| 看十八女毛片水多多多| 男女做爰动态图高潮gif福利片| 好男人在线观看高清免费视频| 少妇的逼水好多| 男的添女的下面高潮视频| 亚洲av男天堂| 亚洲欧美清纯卡通| 日韩一区二区视频免费看| 神马国产精品三级电影在线观看| 国产高清有码在线观看视频| 尤物成人国产欧美一区二区三区| 色哟哟哟哟哟哟| 少妇裸体淫交视频免费看高清| 男人狂女人下面高潮的视频| 国产黄片美女视频| 黄片wwwwww| av在线天堂中文字幕| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 亚洲成人精品中文字幕电影| av视频在线观看入口| 国产av一区在线观看免费| 两个人的视频大全免费| 成人漫画全彩无遮挡| 亚洲av免费高清在线观看| 中出人妻视频一区二区| 国产伦精品一区二区三区四那| 三级国产精品欧美在线观看| 国产一区亚洲一区在线观看| 春色校园在线视频观看| 又爽又黄a免费视频| 小说图片视频综合网站| 久久这里只有精品中国| 精品久久久久久久末码| 日韩视频在线欧美| 波多野结衣高清无吗| 丝袜美腿在线中文| 亚洲在线观看片| 在线观看66精品国产| 菩萨蛮人人尽说江南好唐韦庄 | 欧美激情久久久久久爽电影| 国产精品无大码| 欧美一区二区亚洲| 国产中年淑女户外野战色| 国产精品99久久久久久久久| 欧美日韩精品成人综合77777| 亚洲av熟女| 性欧美人与动物交配| 亚洲国产日韩欧美精品在线观看| 国内少妇人妻偷人精品xxx网站| 97超视频在线观看视频| 嫩草影院入口| 欧美高清性xxxxhd video| 黄色一级大片看看| 婷婷色综合大香蕉| av卡一久久| 久久人妻av系列| 日韩av不卡免费在线播放| 午夜亚洲福利在线播放| 老女人水多毛片| 亚洲成人久久性| 国内揄拍国产精品人妻在线| 久久久精品欧美日韩精品| 久久99热6这里只有精品| 精品少妇黑人巨大在线播放 | 美女脱内裤让男人舔精品视频 | 麻豆成人午夜福利视频| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩高清在线视频| 欧美+亚洲+日韩+国产| 久久中文看片网| 人妻少妇偷人精品九色| 亚洲人成网站在线观看播放| 99久久九九国产精品国产免费| 婷婷六月久久综合丁香| 亚洲七黄色美女视频| 免费看日本二区| 成年版毛片免费区| 亚洲一级一片aⅴ在线观看| 简卡轻食公司| 久久精品国产亚洲av涩爱 | 亚洲在线自拍视频| 国产午夜精品久久久久久一区二区三区| 一级毛片我不卡| av在线亚洲专区| 亚洲精品色激情综合| 99热全是精品| 少妇高潮的动态图| 亚洲av男天堂| 一本一本综合久久| 天堂av国产一区二区熟女人妻| 91午夜精品亚洲一区二区三区| 少妇的逼好多水| 国产真实乱freesex| 伦理电影大哥的女人| 又粗又硬又长又爽又黄的视频 | 国产av麻豆久久久久久久| 久久亚洲精品不卡| 51国产日韩欧美| 老熟妇乱子伦视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产欧美在线一区| 美女高潮的动态| 成年版毛片免费区| 久久精品国产亚洲av香蕉五月| 国产精品美女特级片免费视频播放器| 国产伦理片在线播放av一区 | 99久久精品国产国产毛片| 成人永久免费在线观看视频| 国产精品人妻久久久影院| 99久久人妻综合| 最好的美女福利视频网| 91久久精品国产一区二区成人| 亚洲av电影不卡..在线观看| 你懂的网址亚洲精品在线观看 | 美女国产视频在线观看| 免费在线观看成人毛片| 三级经典国产精品| 久久久久久大精品| 国产色爽女视频免费观看| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 亚洲美女搞黄在线观看| 精品一区二区三区人妻视频| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 免费观看在线日韩| 中国国产av一级| 直男gayav资源| 日韩av不卡免费在线播放| 亚洲av熟女| 午夜精品一区二区三区免费看| 自拍偷自拍亚洲精品老妇| 国模一区二区三区四区视频| 欧美日本视频| 久久久久久久久久黄片| 内地一区二区视频在线| 亚洲精品乱码久久久久久按摩| 国产高清有码在线观看视频| 国产大屁股一区二区在线视频| 亚洲五月天丁香| 国产精品一区二区三区四区免费观看| 国产精品av视频在线免费观看| 国产精品嫩草影院av在线观看| 在线国产一区二区在线| 久久精品国产自在天天线| 亚洲最大成人手机在线| 欧美日本视频| 国产精品久久久久久精品电影| 久99久视频精品免费| 亚洲最大成人中文| 国产成年人精品一区二区| 欧美丝袜亚洲另类| 亚洲经典国产精华液单| 成年女人永久免费观看视频| 夜夜夜夜夜久久久久| 国产高清不卡午夜福利| 一个人免费在线观看电影| a级一级毛片免费在线观看| 久久久久久久久大av| 老司机福利观看| 中文欧美无线码| 乱系列少妇在线播放| 看十八女毛片水多多多| av又黄又爽大尺度在线免费看 | 亚洲精品粉嫩美女一区| 国产日本99.免费观看| 人体艺术视频欧美日本| 我要看日韩黄色一级片| 美女被艹到高潮喷水动态| 美女黄网站色视频| 国产91av在线免费观看| 欧美日韩在线观看h| 韩国av在线不卡| 一区二区三区高清视频在线| 日韩欧美精品免费久久| 欧美激情久久久久久爽电影| АⅤ资源中文在线天堂| 99久国产av精品| 国产极品精品免费视频能看的| 精品久久久久久久末码| 国产精品一区二区三区四区免费观看| 国产精品女同一区二区软件| 99久久精品一区二区三区| 亚洲自偷自拍三级| 国产大屁股一区二区在线视频| 搡女人真爽免费视频火全软件| 国产日本99.免费观看| 国产成人精品久久久久久| 欧美日韩乱码在线| 黄色欧美视频在线观看| 少妇丰满av| 久久精品国产亚洲网站| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| 亚洲最大成人中文| 综合色丁香网| 特级一级黄色大片| 青春草亚洲视频在线观看| 亚洲在线自拍视频| 国产精品一及| 久久久久久久久久久丰满| 在线观看一区二区三区| 日本黄色片子视频| 大型黄色视频在线免费观看| 亚洲欧美清纯卡通| 国产黄片美女视频| 搡老妇女老女人老熟妇| 人人妻人人看人人澡| 长腿黑丝高跟| 成人美女网站在线观看视频| 久久精品国产亚洲av香蕉五月| 久久草成人影院| 欧美日韩在线观看h| 美女高潮的动态| 欧美日韩在线观看h| 久久精品久久久久久久性| 波多野结衣高清无吗| 美女 人体艺术 gogo| 亚洲精品日韩在线中文字幕 | 国产精品一区二区在线观看99 | 中出人妻视频一区二区| 久久人人精品亚洲av| 精品无人区乱码1区二区| 国内精品美女久久久久久| 69av精品久久久久久| 99久久精品一区二区三区| 中文精品一卡2卡3卡4更新| 国产在线精品亚洲第一网站| 国产精品久久久久久精品电影| 国产精品永久免费网站| 日韩制服骚丝袜av| 午夜久久久久精精品| 日韩欧美三级三区| 亚洲在线观看片| 亚洲av成人av| 亚洲三级黄色毛片| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 亚洲精品久久久久久婷婷小说 | 99热网站在线观看| 国产亚洲av片在线观看秒播厂 | 如何舔出高潮| 国产伦精品一区二区三区视频9| 色尼玛亚洲综合影院| 亚洲精品成人久久久久久| 日本撒尿小便嘘嘘汇集6| 一进一出抽搐gif免费好疼| 亚洲欧美日韩高清专用| 欧美三级亚洲精品| 一本久久精品| 国产老妇女一区| 欧美xxxx黑人xx丫x性爽| 大香蕉久久网| 99国产精品一区二区蜜桃av| 日韩,欧美,国产一区二区三区 | 久久精品综合一区二区三区| 大又大粗又爽又黄少妇毛片口| 亚州av有码| 老女人水多毛片| 亚洲一级一片aⅴ在线观看| 最近最新中文字幕大全电影3| 中文精品一卡2卡3卡4更新| 亚洲三级黄色毛片| 九色成人免费人妻av| 久久亚洲国产成人精品v| 黄色视频,在线免费观看| 中文字幕av在线有码专区| 国产一区二区亚洲精品在线观看| 中文字幕久久专区| 国产男人的电影天堂91| av在线老鸭窝| 久久久久久久久久黄片| 国产高清三级在线| 国产精品国产三级国产av玫瑰| av黄色大香蕉| 久久久久久伊人网av| 亚洲国产精品久久男人天堂| 成人国产麻豆网| 亚洲图色成人| 丰满的人妻完整版| 日本三级黄在线观看| 99久久精品热视频| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 蜜桃久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 国产91av在线免费观看| 午夜a级毛片| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产自在天天线| 亚洲久久久久久中文字幕| 黄片wwwwww| 欧洲精品卡2卡3卡4卡5卡区| 久久亚洲精品不卡| 色尼玛亚洲综合影院| 丰满的人妻完整版| 国产亚洲av片在线观看秒播厂 | 变态另类丝袜制服| 99热这里只有是精品50| 免费电影在线观看免费观看| 久久精品国产鲁丝片午夜精品| 色噜噜av男人的天堂激情| 欧美成人精品欧美一级黄| 亚洲乱码一区二区免费版| 成人二区视频| 亚洲国产精品sss在线观看|