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

    Noncollinear phase-matching geometries in ultra-broadband quasi-parametric amplification

    2022-05-16 07:10:08JiWang王佶YanqingZheng鄭燕青andYunlinChen陳云琳
    Chinese Physics B 2022年5期
    關(guān)鍵詞:陳云

    Ji Wang(王佶) Yanqing Zheng(鄭燕青) and Yunlin Chen(陳云琳)

    1Institute of Applied Mic-Nano Materials,School of Science,Beijing Jiaotong University,Beijing 100044,China

    2School of Materials Science and Chemical Engineering,Ningbo University,Ningbo 315211,China

    Keywords: quasi-parametric amplification,noncollinear phase-matching geometry,nonlinear optics,ultrafast optics

    1. Introduction

    Ultra-short high-intensity lasers based on chirped pulse amplification (CPA) technique have been widely studied. As a result, 10-petawatt-class lasers have been achieved, such as the Shanghai super-intense ultrafast laser facility (SULF,339 J, 21 fs),[1]the European Extreme Light Infrastructure(ELI, 250 J, 22.5 fs, 1/60 Hz) 10 PW facilities,[2]and the Apollon (150 J, 15 fs, 1/60 Hz) 10 PW facility.[3]Several PW-class lasers based on optical parametric chirped-pulse amplification (OPCPA) have been developed.[4–6]OPCPA has great potential in the next-generation hundred-petawatt class lasers because of its high gain, high contrast, and large gain bandwidth.[4]A picosecond-pumped OPCPA can provide a larger bandwidth than a nanosecond-pumped OPCPA because a thinner crystal can be applied with a larger pump intensity.[7]In a nanosecond-pumped OPCPA, the pump intensity is limited because of the damage threshold of nonlinear crystals.[8]Nanosecond-pumped OPCPA is widely used in high-power laser facilities,[4–6,9]whereas its bandwidth has been limited. Recently,a new scheme called quasi-parametric chirped pulse amplification (QPCPA) shows great potential in broadband amplification. It can amplify the signal to its theoretical conversion efficiency limit even with a large phase mismatch.[10,11]Here the theoretical conversion efficiency limit means that the conversion efficiency can reach the limit only constrained by the conservation of energy. The high tolerance towards the walk-off and thermal dephasing in QPCPA within a long propagation length is also discovered.[12]The difference between QPCPA and OPCPA is that the nonlinear crystal used in QPCPA has idler absorption. In OPCPA,the conversion efficiency is limited by the back conversion when there are spatial and temporal distributions.[10]However,this limitation does not exist in QPCPA because the idler absorption suppresses the back conversion process. Up to date,the signal spectrum together with the amplification process in QPCPA is not fully investigated yet.In this article,we describe the signal evolution process in QPCPA to show the origin of the stability towards phase mismatch and the origin of the spectral phase in QPCPA which is different from the spectral phase in OPCPA.We have also investigated the gain factor and the signal spectrum to show its potential in broadband amplification at a limited pump intensity. The numerical simulation results demonstrate that the pulse duration can achieve 7.92 fs at a pump intensity of 5 GW/cm2in a doped-YCa4O(BO3)3(YCOB)crystal by applying a novel triple-wavelength phasematching geometry.

    2. Description of the signal evolution in QPCPA

    Firstly, we analyze the amplification process by a set of simplified three-wave-mixing equations:

    Fig. 1. The description of signal amplification under extreme conditions to show the two different back conversion mechanisms. (a) The description to the amplification process when |Δk|?2G. (b) The corresponding intensity changes of signal, idler and pump in QPCPA calculated with an initial gain factor of 450 m-1, Δk=200 m-1, Ip/Is =3200 and α =200 m-1 with fourth order super-Gaussian profiles.(c)The description of signal amplification when|Δk|?2G. (d)The corresponding intensity changes in QPCPA,which is calculated with an initial gain factor of 450 m-1,Δk=2000 m-1,Ip/Is=3200 and α =200 m-1.

    As shown in Figs. 1(a) and 1(b), the two mechanisms dominate when|Δk|?2Gand|Δk|?2G,respectively.When|Δk|<2G, the signal can be similarly amplified in both OPCPA and QPCPA.For a specific condition Δk=0,the signal can be continually amplified until the pump is fully depleted. When the phase mismatch is much smaller than the gain factor and ignore the temporal and spatial distributions of the pump and the signal, the pump can be fully depleted or nearly depleted. Then a new pump field is generated by the signal and idler with an electric field direction opposite to the origin pump, which is equivalent to aπphase shift on the pump, as shown in Fig. 1(a). The newly generated pump gradually depletes the signal and idler and makes the signal intensity in OPCPA changing periodically. In QPCPA,the idler is partially absorbed during propagation. When the pump converts to signal and idler,it obeys the law of conservation of energy,which is also called the Manley–Rowe relations.Thus,if the temporal and spatial distributions are ignored,the intensity of the signal without phase mismatch can eventually reach the theoretical conversion efficiency limit no matter whether there is idler absorption or not. During the back conversion, the signal decreases until the idler fully converts to pump, which also obeys the Manley–Rowe relations. Hence,the remaining signal intensity in QPCPA is larger than in OPCPA,while the converted pump intensity in QPCPA is less than in OPCPA.After a long-distance propagation,both the idler and the pump are negligible while the signal is kept at the theoretical conversion efficiency limit,as shown in Fig.1(b).

    The relationship between the phase mismatch and the gain factor finally becomes|Δk|>2Gbecause the pump is gradually consumed with the amplification process. We shall consider the two kinds of processes that happen in the threewave mixing process after the signal and the pump produces a new idler component with an envelope dAiin a very short distance dz, and we temporally ignore the existed idler components. The newly generated idler component converts into the signal within a very short distance,or it propagates a long distancezbefore it converts to the signal. The envelope of the newly generated idler component carries a phaseφ(ApA*s)-Δkz0,z0is the position where the new idler component is generated. In the first case,the immediately generated new signal component carries a phaseφ(ApA*pAs)=φ(As), which is the same as the existing signal component. In the second case,the idler propagates a lengthzand then converts to signal at the positionz1=z0+z. In this article,we call these accumulated idler components‘the formerly generated idler’,and call the accumulated signal components (‘the formerly generated signal’). The envelope of the formerly generated idler component generated atz0carries a phaseφ(ApA*s)-Δkz0when it propagates toz1(notice thatAs,AiandApdo not carry the phase produced by the propagation, and Δkzterm carries this phase component is carried). The envelope of the signal component converted by the formerly generated idler component carries a phaseφ(ApA*pAs)-Δk(z1-z0)=φ(As)-Δkz. We want to discuss the condition|Δkz|>π, which leads to the back conversion.In OPCPA,the intensity of the amplified signal is nearly the same as the intensity of the formerly generated signal within the distancelefflimited by the parametric bandwidthleff<π/|Δk|. Then, the phase difference between the newly generated signal and the formerly generated signal gradually becomes larger thanπ/2. When the phase difference is larger thanπ/2, the idler and the signal back convert to the pump,and therefore the signal cannot be amplified. We name the phase of the signal,the idler,and the pump asφs,φi,andφp, respectively. The phase mismatch Δkzhas been considered inφs,i,p. Presume that the initial phase of the signal,the idler and the pump are the same,then when the phase mismatch is in the range of 0<|φs+φi-φp|<π/2,the signal and the idler is amplified together. When the phase mismatch is in the range ofπ/2<|φs+φi-φp|<π,the difference between the phase of formerly generated signal (idler) and the phase of newly generated signal(idler)is larger thanπ/2. Then the signal and the idler continually back convert to the pump until the later generated idler(signal)completely cancel out the formerly generated idler(signal). The difference between the phase of the newly generated signal and the phase of the existing signalφp-φi-φsvaries periodically aroundπ/2,and the signal intensity varies with the same period.

    In QPCPA, when the phase mismatch is in the range ofπ/2<|φs+φi-φp|<π, the phase mismatch of the existing signal and idler components is large,and the existing idler components are gradually absorbed. The phase mismatch of the newly generated signal and idler components is relatively smaller. The proportion of subsequently generated idler increases with the enhancement of the absorption of the idler,and the back conversion is inhibited. The amplification is suppressed in the next 0<|φs+φi-φp|<π/2 interval because the idler generated in the previousπ/2<|φs+φi-φp|<πinterval is not fully depleted. Finally,due to the idler absorption,the later generated idler intensity ratio to the existing idler intensity remains a constant greater than 1 in a|Δk|=2πperiod. The corresponding phase differenceφs+φi-φpremains at a value less thanπ/2, and the back conversion gradually disappears. The signal’s withdrawal in QPCPA is less than that in OPCPA, ensuring the signal increases. When the signal intensity is increasing, the phase differenceφs+φi-φpis further reduced and accelerates the amplification to exponential growth,as shown in Fig.1(c). The exponential growth speed,i.e.,the effective gain factor(the gain factor considering the phase mismatch and the absorption),depends on the ratio of the later generated idler and the existing idler. Therefore the effective gain factor increases with the absorption coefficient until the absorption displays more inhibition than promotion. This phenomenon has also been found and named the‘self-locked phase’.[11]This effect ensures that every component in the signal can be amplified to its theoretical conversion efficiency limit in a long enough nonlinear crystal, even if there is a large phase mismatch. An extreme situation is shown in Fig.1(d),which shows the back conversion gradually disappears in QPCPA even if there is a large phase mismatch Δk=2000 m-1and a short coherence lengthπ/Δk <2 mm.

    3. Broadband amplification and the phasematching strategy in QPCPA

    In a real system,the crystal length is limited. The crystal length required to ensure the conversion efficiency is longer than the crystal length required to maintain the parametric bandwidth in OPCPA.When the effective gain factor is close to zero due to the pump laser depletion,the signal components in the large phase mismatch regions are no longer amplified,as shown in Figs.2(a)and 2(b). Thus,the spectral bandwidth of the amplified signals may be narrower than the seed spectrum in a high efficient nanosecond-pumped OPCPA stage. However,the effective gain factor in the large phase mismatch region does not vanish in QPCPA. This property ensures ultrabroadband amplifications without spectrum narrowing.

    Fig. 2. The effective gain factor with the phase mismatch for (a)G=1000 m-1 and(b)G=500 m-1.

    Generally, the length of a nonlinear crystal is not unlimited. It is essential to investigate the gain factor in QPCPA and show its potential at limited crystal length. The exact expression of signal can be calculated with Eqs. (1)–(3)when the pump depletion can be ignored, and the coefficient of exponential term ofAsis (-j×Δk/2+(G-(Δk+j×α/2)2/4)1/2-α/4)z. Here the real part of (G-(Δk+j×α/2)2/4)1/2provides the gain and the imaginary part provides the phase, as described in our previous work.[12]When Δk ≈0, the QPCPA gain factor can be expressed as(G2-Δk2/4+α2/16)1/2-α/4,as described in Ref.[11],but this expression is lack of accuracy when Δkis large. When the phase mismatch|Δk|<2G, the OPCPA effective gain factorGOPCPA=(G2-Δk2/4)1/2is larger than the QPCPA effective gain factorGQPCPA=Re(G-(Δk+j×α/2)2/4)1/2-α/4.When|Δk|≈2G, the OPCPA effective gain factor is close to zero and shows a cutoff of the gain bandwidth. The QPCPA gain factor decreases to about(αG)1/2/2-α/4 and the maximum effective gain factor is 0.3Gatα=2G.When|Δk|>2G,the OPCPA cannot amplify the signal, and the signal intensity shows periodical oscillation. In contrast,the QPCPA can continually amplify the signal to the theoretical conversion efficiency limit. The results show that in QPCPA, the components with a large phase mismatch|Δk|≈2Gcan be amplified to the theoretical conversion efficiency limit using a crystal 3–4 times longer than the crystal used in OPCPA. In contrast,these components cannot be amplified in OPCPA. The phase mismatch increases rapidly with bandwidth extension so that a long crystal with a large absorption is required to amplify the components where|Δk|>2G. When the gain factorGdecreases,the QPCPA effective gain factorGQPCPAdecreases more than the factorGso that a longer crystal is required to amplify all possible components in QPCPA.If a QPCPA system is designed to have the same bandwidth as an OPCPA system but suppresses the back conversion,it just requires the absorption to be large enough to ensure that the residual idler is much smaller than the signal. After the pump has been totally converted to signal and idler, the idler continually decreases with a factor-α/4. We presume|Ei|/|Es|=1/e2≈e-α(z-ln(Ep/Es)/G0)/4, whereGis the initial gain factor without pump depletion, so that the required absorption can be expressed asα=8/(z-ln(Ep/Es)/G). For a 50-mm-long crystal with a gain factor of 450 m-1and|Ip|/|Is|~3200,the absorptionαshould be larger than 195 m-1. Suppose that a QPCPA is designed to provide an ultra-broadband amplification, then the maximum absorption should be at least 2G0=900 m-1to support the amplification at the point|Δk|≈2G,which puts forward a higher requirement to the crystal doping.

    A large gain factor can easily support an ultrabroadband spectrum in OPCPA or QPCPA towards few-cycle pulses. However, a high gain factor can hardly be achieved at nanosecond scale in the commonly used nonlinear crystals such as BBO (<1.2 GW/cm2@532 nm, 8 ns),[14]YCOB (7.24 GW/cm2@532 nm, 8.5 ns),[15]or Sm:YCOB(6.41 GW/cm2@532 nm, 8.5 ns)[16]because of the damage threshold. The pump intensity highly limits the bandwidth in nanosecond amplification stages,[6]and affects the pulse duration in petawatt laser facilities. In these crystals,the gain factorGcan hardly achieve 500 m-1, which limits the gain bandwidth in OPCPA. The development of petawatt lasers leads to higher intensities in amplifiers. Strong nonlinear effects in picosecond-pumped amplifiers such as selffocusing may become serious. Moreover, these picosecondpumped amplifiers may require extremely large optics in nextgeneration petawatt lasers. However, nanosecond-pumped amplifiers can avoid these nonlinear effects and reduce the size of optics. Therefore, it is reasonable to discuss the amplification at a limited pump intensity, e.g.,<5 GW/cm2in YCOB crystals. The noncollinear geometry (noncollinear OPCPA) can achieve a larger bandwidth in OPCPA than the collinear geometry. The noncollinear geometry in OPCPA aims to adjust the first-order derivation of wave-vector mismatch d(Δk)/dω|(ω=ωs)=0 while maintaining the wavevector mismatch at the central wavelength Δk|(ω=ωs)=0.By adjusting the signal wavelength in the system,the residual second-order wave-vector mismatch term becomes zero. The wavevector mismatch becomes the third-order mismatch term and a large bandwidth can be achieved. This scheme is called ultra-broadband OPCPA.[17]The phase-matching bandwidth in ultra-broadband OPCPA can be derived from Ref. [17]when ignoring the pump depletion as

    Fig.3.The Gaussian-shaped seed spectrum has a Fourier transform limit(F-T limit) of 10 fs. The triple-wavelength phase-matching QPCPA operates at a noncollinear angle of 2.68°,with a larger bandwidth than in the noncollinear OPCPA (a=2.81°). The noncollinear angles calculated with Eq. (14) use the refractive index of Sm:YCOB. When a=2.81°, GV=169 fs, GVD=-6×103 fs2 and k(3) =-3×105 fs3. When a=2.68°, GV=5×103 fs,GVD=3×103 fs2 and k(3) =-3×105 fs3. See the text for definitions of GV and GVD.

    When the maximum acceptable phase mismatch has been determined,the maximum GV of the triple-wavelength phasematching QPCPA is determined by|Δk(ω1)|<|Δkupper|,and the parametric bandwidth can be solved withωpm=±(6|GV|/|k(3)|)1/2. When the noncollinear angle is away from the broadband condition where the phase mismatch contributed by the GV is larger than the phase mismatch contributed by thek(3), the phase-matching bandwidth is limited by GV:

    whereωmaxis the largest frequency component in the signal spectrum, the bandwidth decreases rapidly away from the ultra-broadband noncollinear geometry (i.e.,|GVω|?|k(3)ω3/6| within the signal spectrum). Here we focus on the bandwidth in the ultra-broadband triple-wavelength phasematching geometry, and the corresponding phase-matching bandwidth can be estimated as follows:

    Equation (11) shows that a larger gain factor, longer crystal length, and higher signal intensity are helpful to get a larger bandwidth in QPCPA. When QPCPA is operated with the same geometry as an ultra-broadband noncollinear OPCPA,the bandwidth can be solved by

    The phase-matching bandwidth can be derived as follows:

    The phase-matching bandwidth in ultra-broadband noncollinear phase matching geometry is narrower than the bandwidth in triple-wavelength phase-matching geometry,ΔνQPCPA-NPM≈0.73ΔνQPCPA-TPM.

    Fig. 4. The parametric bandwidth changes with noncollinear angle when (a) G0 =450 m-1 and (b) G0 =320 m-1. Here, (a) and (b) are calculated with a crystal length of 50 mm for QPCPA and 10 mm for OPCPA, Ip/Is ≈3200 and α =2G. (c) The parametric bandwidth changing with the gain factor when z=50 mm for QPCPA and z=z0 for OPCPA,Ip/Is ≈3200 and α =2G. (d)The parametric bandwidth changing with the absorption when z=50 mm and Ip/Is ≈3200.

    Figure 4(c) shows the parametric bandwidth versus the gain factorG. The triple-wavelength phase-matching QPCPA can achieve a bandwidth larger than noncollinear QPCPA and noncollinear OPCPA. Only half of the pump intensity is required in triple-wavelength phase-matching QPCPA in order to get the same bandwidth with noncollinear QPCPA. The bandwidth increasing in QPCPA is limited when the gain factor is large enough because QPCPA cannot amplify the components where signal and idler are overlapped. In the triplewavelength phase-matching scheme,the best noncollinear angle is determined by GV. To simplify the question, we can ignore the changes of GVD andk(3)compared to the change of GV when changing the noncollinear angle (the angle isa) and phase-matching angleθaround Δk|(ω=ωs)=0 and GVD=0. The relationship between GV and noncollinear angleacan be expressed as

    whereb= arcsin(ks/ki×sin(a)).ksandkiare the wavevector of the signal and the idler, respectively. When the parametric bandwidth already determines the maximum acceptable phase mismatch atω1, GV can be calculated by Δkmax=±2(1/3)3/2(-6GV/k(3))1/2. And the noncollinear angle can be solved by Eq. (14). The parametric bandwidth changes with the noncollinear angle considering GV, GVD andk(3)are shown in Figs.4(a)and(b). When the gain factor is 450 m-1,the absorption is 900 m-1andIp/Is≈3200 with a 50 mm long crystal, the best noncollinear angle in QPCPA is between 2.6°and 2.7°. When the gain factor is reduced to 320 m-1and the absorption is 640 m-1,the best noncollinear angle is close to 2.7°. The noncollinear angle in OPCPA remains about 2.8°,which does not change with the gain factor.Figure 4(d) shows that a larger absorption coefficient is required to reach the parametric bandwidth maximum with the gain factor increasing,and the absorption required by the maximum parametric bandwidth is approximatelyα=2G.

    4. Results and discussion

    The triple-wavelength phase-matching QPCPA spectra have been simulated with three-wave mixing equations. In the following simulations, we use the refractive index of Sm:YCOB crystals considering up to the third-order material dispersion and assume ideal absorption spectra:

    whereτis the time normalized by the pulse duration of signalts,τ=t/ts;ks,ki,kpare the wave vectors of the signal, the idler and the pump,respectively.

    Fig.5.The parametric bandwidth and the simulated signal bandwidth changing with the propagation length in triple-wavelength phase-matching QPCPA when the gain factor G=450 m-1 [(a),blue]and G=320 m-1 [(b),orange]. The parametric bandwidth and the simulated signal bandwidth changing with the propagation length in(c)noncollinear phase-matching QPCPA(yellow)and in(d)noncollinear phase-matching OPCPA(purple). The initial signal is centered at 810 nm with an F-T limit of 10 fs. The pump is centered at 532 nm,Ip/Is ≈3200,the crystal length z=50 mm and the absorption in QPCPA is 900 m-1. The solid lines are the upper limits of the parametric bandwidths,while the dashed lines are the lower limits,the dotted lines are the simulated signal bandwidth,and the dash-dotted lines are the simulated whole spectrum range of the signal.

    The estimated parametric bandwidth and the simulated signal bandwidth in OPCPA and QPCPA are shown in Fig.5.In triple-wavelength phase-matching QPCPA,the signal bandwidth decreases whenz <z0G0/G′because of the gain distribution. The signal bandwidth rapidly increases to the parametric bandwidth maximum as shown in Fig.5(a)with a gain factorG0=450 m-1. Then the bandwidth increase is slowed down because the parametric bandwidth maximum requires a larger GV due to a longer propagation length. When the gain factor decreases to 320 m-1but the noncollinear angle remains 2.68°, the phase mismatch maximum located atω1is larger than the requirement of the parametric bandwidth so that the simulated signal bandwidth cannot reach the parametric bandwidth maximum,as shown in Fig.5(b). However,because the components at the phase-matching pointωpmcan be well amplified, the whole spectrum range can be larger than the parametric bandwidth maximum and reach the spectrum range which is supported by the noncollinear angle of 2.68°,as shown in Fig.5(a).When the QPCPA is operated ata=2.81°,the simulated signal bandwidth is close to the parametric bandwidth maximum as shown in Fig.5(c). In contrast,the signal bandwidth of noncollinear OPCPA can only reach the parametric bandwidth maximum in a small range aroundz=z0as shown in Fig.5(d).

    When the QPCPA works at a noncollinear anglea=2.81°as the experiment in Ref.[10],and assume a seed pulse duration of 10 fs and a pump intensity of about 5 GW/cm2, the amplified spectrum in QPCPA can cover 750–930 nm, which is slightly broader than in OPCPA,as shown in Fig.6(a). The OPCPA spectrum is coincident with the experiment and the theoretical results using YCOB crystals in Refs. [20,21] as shown in Fig. 6(b). The bandwidth in QPCPA is continually increasing with the propagation length because the components with larger phase mismatch can be amplified to their theoretical conversion efficiency limit. However, the bandwidth extension is small because the phase mismatch is increasing rapidly. In OPCPA,the signal intensity and the spectrum bandwidth cannot be kept to their maximum with an improper crystal length. In QPCPA,the signal intensity and the bandwidth always are kept at their maximum when the crystal is long enough. When the noncollinear angle is set toa=2.68°, the spectrum in QPCPA can be extended to 700–950 nm,as shown in Fig.6(c). The amplified spectrum has a Fourier transform limit of~8 fs,even shorter than the seed. It is found that more short-wavelength components can be amplified ata=2.68°compared toa=2.81°. Therefore, this geometry is suitable for amplifying seeds that come from Tisapphire lasers or the white-light super-continuum generated by self-phase modulation. The former (Ref. [6]) has a spectrum covering 700–900 nm,and the latter(Ref.[22])has more short-wavelength components than the long-wavelength components.When the OPCPA works ata=2.68°,different spectral components cannot be amplified to their maximum with the same crystal length as shown in Fig.6(d). The bandwidth of OPCPA cannot be extended to the same as QPCPA under the triple-wavelength phase-matching geometry with a moderate gain factor.

    Fig.6. The spectrum with propagation length in(a)QPCPA with a=2.81°,(b)OPCPA with a=2.81°,(c)QPCPA with a=2.68°,and(d)OPCPA with a=2.68°,respectively. The spectrum is calculated with the central wavelength of 810 nm,and the pump intensity of 5 GW/cm2.The seed is centered at 810 nm with an F-T limit of 10 fs. The pump is centered at 532 nm with a fourth order super-Gaussian profile and Ip/Is ≈3200. The idler absorption in QPCPA is fixed at 900 m-1.

    The conversion efficiency and the parametric bandwidth changes with the propagation length in QPCPA and OPCPA are shown in Fig. 7. The triple-wavelength scheme has the highest efficiency because the small components away from the central wavelength can be well amplified. It is impossible to get a large bandwidth by reducing the crystal length in OPCPA when the gain factor is limited to 450 m-1andIp/Is≈3200 in our simulation because the conversion efficiency is negligible. It is important that the bandwidth of QPCPA is increasing together with the conversion efficiency,which means that the signal bandwidth and the signal intensity can both reach their upper limits.

    Fig. 7. The parametric bandwidth (solid lines) and the conversion efficiency (dashed lines) changing with the propagation length of the triple-wavelength phase-matching QPCPA (blue), the noncollinear phase-matching QPCPA(orange)and the noncollinear phase-matching OPCPA (yellow). The pulse duration of the pump is 150 ps with a fourth-order super-Gaussian profile, and the signal is stretched from 10 fs to 50 ps. The initial gain factor G0 ≈450 m-1, Ip/Is ≈3200 and the absorption is fixed at 900 m-1.

    If the pump intensity changes in an established system,the spectrum and conversion efficiency may be changed dramatically. When the pump intensity changes from 5 GW/cm2to 2.5 GW/cm2, because of the unsuitable crystal length, the signal bandwidth in OPCPA is narrowed by 40 nm(as shown in Fig.8(a))and the conversion efficiency decreases to 42.5%compared to its maximum. In QPCPA, the narrowness of the signal bandwidth is much smaller than that in OPCPA,as shown in Figs. 8(b) and 8(c). When the noncollinear anglea=2.81°, the conversion efficiency is about 87% compared to the theoretical conversion efficiency limit. Whena=2.68°,the conversion efficiency remains about 80%,while the spectrum narrowing is negligible.

    It is essential to analyze the nonlinear spectral phase in triple-wavelength phase-matching QPCPA. The origin of the spectral phase accumulation induced by the idler absorption is the different phase accumulation during the pulse transmission between the idler and the signal. When there is a phase mismatch in the amplifying process,an additional phase term∫exp(-jΔkz)dzappears in both newly generated idler and signal components.In OPCPA,this additional phase is equally added to the newly generated signal and idler,which produces a spectral phase-Δkz/2.

    The phase of the amplified signal in OPCPA can be expressed as[23]

    When the absorption is large enough,the phase accumulation is close to the material dispersion from the beginning to the end. When there is a strong idler absorption, the formerly generated idler is quickly absorbed,the newly generated idler has a phase close toΦp-Φs-Δkz. When the newly generated idler converts to signal,the phase of the newly generated signal is nearly the same as that of the existing signal. The signal is only experiencing the material dispersion. If highorder dispersion in the crystals can be ignored,the production of the uncompensable phase from the amplification process is strongly suppressed. Therefore, the corresponding spectral phase can be compensated for, and only the high-order dispersion due to the non-uniformity of the absorption spectrum remains. The dispersion induced by absorption spectrum nonuniformity can be large if a strong absorption peak is placed at a large phase mismatch region and a weak absorption valley is placed at the small phase mismatch region. A uniform absorption spectrum is impossible to be achieved. Therefore strong absorption peaks should be ensured in the small phase mismatch region to decrease the phase difference between those components.

    When a uniform absorption spectrum is used, the spectral phase of signals after amplified by a 50-mm-long crystal can be well controlled by a high order dispersion control up to 6thorder, as shown in Fig. 8(c). However, if there is no idler absorption, the spectral phase cannot be well controlled even with high order dispersion control, as shown in Fig.8(d). In a triple-wavelength phase-matching QPCPA,the required absorption spectrum is possibly formed by multiple absorption peaks. Assuming an absorption spectrum as shown in Fig.8(e),the residual spectral phase is highly modulated by the nonuniformity of the absorption spectrum. In the large dispersion region,the nonlinear spectral phase in QPCPA with such an absorption spectrum is smaller than the nonlinear spectral phase in OPCPA.On the other hand, the remaining spectral phase using a non-uniform absorption spectrum is larger than the nonlinear spectral phase using a uniform absorption spectrum, as shown in Fig. 8(d). Therefore, it leads to minor side pulses,as shown in Fig.8(f).

    Fig. 8. The spectrum and spectral phase after dispersion control up to 6th order of (a) noncollinear OPCPA, (b) noncollinear QPCPA, (c)triple-wavelength phase-matching QPCPA, (d) triple-wavelength phase-matching QPCPA with multi-peaks absorption, respectively. (e) The absorption spectrum in(d). (f)The pulse duration after dispersion control. The seed is centered at 810 nm with an F-T limit of 10 fs. The pump is centered at 532 nm and Ip/Is ≈3200. The absorption in(b)and(c)is fixed at 900 m-1.

    The pulse duration after dispersion management is shown in Fig. 8(f). The pulse duration is 12.5 fs when amplified with the noncollinear OPCPA scheme at the pump intensity of 5 GW/cm2. When the pump intensity is changed to 2.5 GW/cm2, the pulse duration in OPCPA is expended to 14.5 fs. The pulse duration after amplifying by triplewavelength phase-matching QPCPA can achieve 7.92 fs,which is close to its Fourier transform limits. The shortest pulse duration, 7.92 fs, is obtained by a dispersion control up to 6thorder, which can be achieved by applying an acousto-optic programmable dispersive filter (AOPDF). The high-order dispersions in triple-wavelength phase-matching QPCPA are large because of the relatively larger phasemismatch and the absorption distributions in the broadband spectrum. The dispersions are group delay dispersions(GDD)=4418 fs2,third-order dispersions(TOD)=3767 fs3,fourth-order dispersions (FOD) = 5147 fs4, fifth-order dispersions (FiOD) = 104040 fs5and sixth-order dispersions(SOD)=615672 fs6when using a uniform absorption spectrum,respectively. Therefore,dispersion control up to the 3rdorder is required. Meanwhile, the dispersions are GDD =4430 fs2, TOD = 2662 fs3, FOD =-1168 fs4, FiOD =322680 fs5and SOD=578952 fs6when using a multi-peaks absorption spectrum,respectively. As a result,dispersion control up to the 5thorder is necessary. If the dispersions are controlled up to 3rdorder or 4thorder with a uniform absorption spectrum, the uncompensated phase within the spectrum is±0.6π,showing a pulse duration~7.92 fs with small side peaks (about 10% of the maximum intensity). The uncompensated phase within the spectrum is±0.2πwith a dispersion control up to 5thor 6thorder with no side peaks. If the dispersions are controlled up to 3rdorder or 4thorder with a multi-peaks absorption spectrum,the uncompensated phase within the spectrum is±1.8πwhile showing large side peaks(about 60% of the maximum intensity). The uncompensated phase within the spectrum is±0.6πwith a dispersion control up to 5thor 6thorder, side peaks almost disappear. The pump intensity changes almost do not change the pulse duration, which means triple-wavelength phase-matching QPCPA is very insensitive towards pump intensity fluctuations. These results show that QPCPA is tailored to broadband amplification at a limited pump intensity, which has great potential in full-OPCPA petawatt laser systems.

    5. Conclusion

    In conclusion,we have described the evolution of the signal in the quasi-parametric amplification process and proved that QPCPA is suitable for broadband amplification. We have modeled and developed a novel noncollinear phase-matching geometry in QPCPA, named as the triple-wavelength phasematching geometry, and compared this scheme with the noncollinear geometry commonly used in OPCPA.Adjusting the noncollinear angle to the triple-wavelength phase-matching point allows the signal bandwidth to about 260 nm when the pump intensity is limited to 5 GW/cm2in a doped YCOB crystal. This scheme provides a much broader bandwidth in QPCPA than the 150 nm bandwidth in noncollinear OPCPA.The results show that the pulse duration of the signal amplified by the triple-wavelength phase-matching QPCPA can achieve 7.92 fs, which is even shorter than the Fourier transform limit of the seed pulses. Furthermore,the pulse duration and the conversion efficiency are stable when the pump intensity changes from 5 GW/cm2to 2.5 GW/cm2, showing its great potential in few-cycle full-OPCPA petawatt lasers. The quasi-parametric amplifying process can overcome the high order dispersion induced by the phase mismatch. However,the nonuniformity of the absorption spectrum leads to another set of high order dispersion,which does not expend the pulse duration but brings several minor side peaks in the time domain.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China(Grant No.51832009)and the Fundamental Research Funds for the Central Universities, China(Grant No.2019YJS209).

    猜你喜歡
    陳云
    加快構(gòu)建旅游產(chǎn)業(yè)創(chuàng)新生態(tài)系統(tǒng)
    Biased random walk with restart for essential proteins prediction
    基于大數(shù)據(jù)分析與審計的關(guān)系研究
    My plan for new term
    向陳云學(xué)習(xí)錘煉“筆力”
    陳云:我黨干部的楷模
    A Study of ρ-ω Mixing in Resonance Chiral Theory?
    陳云貴:你是泥土你是光
    海峽姐妹(2016年2期)2016-02-27 15:15:59
    TransitivityandCharacterization:AnalysisonDickinTenderisTheNight
    略論陳云執(zhí)政黨黨風(fēng)建設(shè)的思想
    在线a可以看的网站| a级毛片a级免费在线| 久久久水蜜桃国产精品网| 久久天躁狠狠躁夜夜2o2o| x7x7x7水蜜桃| 亚洲av中文字字幕乱码综合| 国产三级中文精品| 国产三级黄色录像| 99在线视频只有这里精品首页| 久久伊人香网站| 中文资源天堂在线| av中文乱码字幕在线| 19禁男女啪啪无遮挡网站| 精品久久久久久久末码| 一卡2卡三卡四卡精品乱码亚洲| 国产精品亚洲美女久久久| 嫩草影视91久久| 亚洲精品美女久久av网站| 国产精华一区二区三区| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品女同一区二区软件 | 亚洲av片天天在线观看| 亚洲 国产 在线| 国产av不卡久久| 午夜精品在线福利| 久久久久久久久免费视频了| 午夜两性在线视频| 久久久久久久久久黄片| 免费观看的影片在线观看| 国内久久婷婷六月综合欲色啪| 国内精品一区二区在线观看| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 精品午夜福利视频在线观看一区| 亚洲精品一卡2卡三卡4卡5卡| 啦啦啦免费观看视频1| 我要搜黄色片| 成人亚洲精品av一区二区| 国产精品乱码一区二三区的特点| 日本 av在线| 国产成人精品无人区| 成年女人毛片免费观看观看9| 99久久99久久久精品蜜桃| 精品电影一区二区在线| 一个人免费在线观看的高清视频| 首页视频小说图片口味搜索| 色老头精品视频在线观看| 97人妻精品一区二区三区麻豆| 叶爱在线成人免费视频播放| 成人亚洲精品av一区二区| 成人国产综合亚洲| 中文字幕人成人乱码亚洲影| 免费看十八禁软件| 99热这里只有是精品50| 国产爱豆传媒在线观看| 噜噜噜噜噜久久久久久91| 18禁黄网站禁片午夜丰满| 一级毛片精品| 97碰自拍视频| 亚洲欧美一区二区三区黑人| 国产av一区在线观看免费| 老司机在亚洲福利影院| 精品国产美女av久久久久小说| 最近最新免费中文字幕在线| 精品一区二区三区av网在线观看| 91字幕亚洲| 一个人免费在线观看的高清视频| 婷婷丁香在线五月| 婷婷精品国产亚洲av| 久久亚洲精品不卡| 18禁裸乳无遮挡免费网站照片| 国产精品永久免费网站| 国产极品精品免费视频能看的| 亚洲av免费在线观看| 给我免费播放毛片高清在线观看| 成人精品一区二区免费| 亚洲美女黄片视频| 免费人成视频x8x8入口观看| 亚洲人成网站高清观看| 国内少妇人妻偷人精品xxx网站 | 国产精品亚洲一级av第二区| 久久婷婷人人爽人人干人人爱| 久久久久国内视频| 99精品在免费线老司机午夜| 日韩精品青青久久久久久| 噜噜噜噜噜久久久久久91| 成人鲁丝片一二三区免费| 欧美日韩瑟瑟在线播放| 男女下面进入的视频免费午夜| 欧美黑人巨大hd| 亚洲欧洲精品一区二区精品久久久| 女同久久另类99精品国产91| 床上黄色一级片| 热99re8久久精品国产| 午夜日韩欧美国产| h日本视频在线播放| 18禁黄网站禁片免费观看直播| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 亚洲国产高清在线一区二区三| 亚洲精品色激情综合| 久久久色成人| 黑人操中国人逼视频| 午夜福利18| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 在线观看免费视频日本深夜| 欧美黑人欧美精品刺激| 午夜激情福利司机影院| www日本在线高清视频| 国产一区在线观看成人免费| 男插女下体视频免费在线播放| 亚洲第一电影网av| 夜夜看夜夜爽夜夜摸| 波多野结衣高清作品| 欧美日韩一级在线毛片| 亚洲av成人一区二区三| 国产野战对白在线观看| a级毛片在线看网站| 亚洲欧洲精品一区二区精品久久久| 999久久久精品免费观看国产| 变态另类丝袜制服| 国内精品久久久久久久电影| 一级黄色大片毛片| 天堂动漫精品| 国产69精品久久久久777片 | 午夜激情福利司机影院| 欧美不卡视频在线免费观看| 国产av不卡久久| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲精品国产色婷小说| 亚洲国产看品久久| 黄片大片在线免费观看| 99精品在免费线老司机午夜| 国产精品香港三级国产av潘金莲| 嫩草影院入口| 日韩欧美在线乱码| 丰满的人妻完整版| 国产成人av教育| 又黄又粗又硬又大视频| 久久精品夜夜夜夜夜久久蜜豆| www.www免费av| 免费观看精品视频网站| 日本五十路高清| 狠狠狠狠99中文字幕| 国产精品 欧美亚洲| 日本黄大片高清| 90打野战视频偷拍视频| www.熟女人妻精品国产| 男人舔奶头视频| 制服丝袜大香蕉在线| 熟女少妇亚洲综合色aaa.| 国产亚洲精品久久久久久毛片| 99在线人妻在线中文字幕| 亚洲欧美激情综合另类| 又紧又爽又黄一区二区| 视频区欧美日本亚洲| 精品午夜福利视频在线观看一区| 成人av一区二区三区在线看| 精品免费久久久久久久清纯| 国产成人av激情在线播放| 男人舔女人的私密视频| cao死你这个sao货| 免费观看精品视频网站| 国产不卡一卡二| 禁无遮挡网站| 搞女人的毛片| 97超视频在线观看视频| e午夜精品久久久久久久| 欧美乱色亚洲激情| 美女高潮喷水抽搐中文字幕| aaaaa片日本免费| 黄色日韩在线| 欧美一区二区国产精品久久精品| 一本精品99久久精品77| 国产av在哪里看| 一个人观看的视频www高清免费观看 | 岛国在线观看网站| 麻豆国产97在线/欧美| 亚洲一区高清亚洲精品| 日本黄色视频三级网站网址| 搡老妇女老女人老熟妇| 亚洲在线观看片| 亚洲真实伦在线观看| 啪啪无遮挡十八禁网站| 丁香六月欧美| 九九久久精品国产亚洲av麻豆 | 国产视频一区二区在线看| 国产一区二区激情短视频| 老司机午夜福利在线观看视频| 99视频精品全部免费 在线 | 一区二区三区激情视频| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 此物有八面人人有两片| 99riav亚洲国产免费| 亚洲国产欧美网| 国内揄拍国产精品人妻在线| 好男人在线观看高清免费视频| 不卡av一区二区三区| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 国产av麻豆久久久久久久| 日韩 欧美 亚洲 中文字幕| 久久久国产成人免费| 少妇熟女aⅴ在线视频| 波多野结衣高清作品| 91麻豆av在线| 国产三级中文精品| a级毛片在线看网站| 国产成人精品久久二区二区免费| 国产精品久久久久久精品电影| 熟女人妻精品中文字幕| 国产亚洲av高清不卡| 亚洲av电影在线进入| 色尼玛亚洲综合影院| 午夜福利在线在线| 久久久久九九精品影院| 熟妇人妻久久中文字幕3abv| 亚洲精品美女久久久久99蜜臀| 色视频www国产| 国产精品99久久久久久久久| 少妇人妻一区二区三区视频| 国产精品爽爽va在线观看网站| 国产淫片久久久久久久久 | 天天一区二区日本电影三级| 国产精品亚洲一级av第二区| av国产免费在线观看| 国产精品av久久久久免费| 女警被强在线播放| 国产欧美日韩精品一区二区| 老司机在亚洲福利影院| 中文字幕人成人乱码亚洲影| 亚洲精品中文字幕一二三四区| 久久久久性生活片| 欧美最黄视频在线播放免费| 亚洲色图av天堂| 国产一区在线观看成人免费| 网址你懂的国产日韩在线| 日韩欧美三级三区| 午夜福利视频1000在线观看| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 老汉色av国产亚洲站长工具| av视频在线观看入口| 在线永久观看黄色视频| 99国产精品一区二区蜜桃av| 亚洲在线观看片| av中文乱码字幕在线| 色精品久久人妻99蜜桃| 成人精品一区二区免费| www日本黄色视频网| 国产成人影院久久av| 一级毛片高清免费大全| 久久国产乱子伦精品免费另类| 最近最新免费中文字幕在线| 精品久久久久久久人妻蜜臀av| 午夜福利18| 青草久久国产| 久久精品国产综合久久久| 国产精品av视频在线免费观看| 亚洲成av人片免费观看| 亚洲av中文字字幕乱码综合| 少妇丰满av| 久久午夜亚洲精品久久| 国产精品电影一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 他把我摸到了高潮在线观看| 久久午夜亚洲精品久久| 九色国产91popny在线| 精品福利观看| 国产成年人精品一区二区| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看| 变态另类成人亚洲欧美熟女| 香蕉av资源在线| 香蕉国产在线看| 亚洲 欧美一区二区三区| 欧美一区二区精品小视频在线| 美女高潮喷水抽搐中文字幕| 真人一进一出gif抽搐免费| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| 午夜福利在线观看吧| 一本精品99久久精品77| 久久久久久九九精品二区国产| 精品免费久久久久久久清纯| 国产 一区 欧美 日韩| 神马国产精品三级电影在线观看| 久久亚洲真实| 99精品在免费线老司机午夜| 美女免费视频网站| 这个男人来自地球电影免费观看| 免费在线观看日本一区| 午夜两性在线视频| 99在线人妻在线中文字幕| 真人一进一出gif抽搐免费| 人妻久久中文字幕网| 国产视频一区二区在线看| 国产91精品成人一区二区三区| 精品一区二区三区视频在线 | 熟女人妻精品中文字幕| 香蕉久久夜色| 美女午夜性视频免费| 久久久久久久午夜电影| 女人被狂操c到高潮| 国产激情欧美一区二区| 午夜福利高清视频| 成年女人永久免费观看视频| 日本撒尿小便嘘嘘汇集6| 两性午夜刺激爽爽歪歪视频在线观看| avwww免费| 美女cb高潮喷水在线观看 | 亚洲精品粉嫩美女一区| aaaaa片日本免费| 国产午夜精品论理片| 少妇人妻一区二区三区视频| 免费搜索国产男女视频| 免费观看精品视频网站| 亚洲国产精品sss在线观看| 母亲3免费完整高清在线观看| 91老司机精品| 午夜福利18| 嫁个100分男人电影在线观看| 国产欧美日韩精品一区二区| 国内揄拍国产精品人妻在线| 欧美成狂野欧美在线观看| 国产亚洲精品av在线| av欧美777| bbb黄色大片| 国产精品国产高清国产av| 国产精品爽爽va在线观看网站| 亚洲av成人精品一区久久| 亚洲乱码一区二区免费版| 久久久色成人| 变态另类成人亚洲欧美熟女| 91av网站免费观看| 中文字幕高清在线视频| 丰满的人妻完整版| 欧美丝袜亚洲另类 | 国产精品一及| 成年免费大片在线观看| 我的老师免费观看完整版| 最近最新免费中文字幕在线| 亚洲国产日韩欧美精品在线观看 | 亚洲精品粉嫩美女一区| 99久久综合精品五月天人人| 97碰自拍视频| 精品欧美国产一区二区三| 国产精品一区二区精品视频观看| 叶爱在线成人免费视频播放| 亚洲avbb在线观看| 久久国产精品影院| 国产伦人伦偷精品视频| 看黄色毛片网站| 亚洲精品美女久久久久99蜜臀| 亚洲男人的天堂狠狠| 久久精品人妻少妇| 淫秽高清视频在线观看| 观看免费一级毛片| 日本免费一区二区三区高清不卡| 亚洲,欧美精品.| 久久精品国产99精品国产亚洲性色| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 成人18禁在线播放| 男人舔女人的私密视频| 国产精品久久久久久人妻精品电影| 可以在线观看毛片的网站| 亚洲国产看品久久| 亚洲人成电影免费在线| 亚洲欧美精品综合一区二区三区| 嫩草影院精品99| 亚洲精品乱码久久久v下载方式 | www日本黄色视频网| 久9热在线精品视频| 精品不卡国产一区二区三区| 18禁美女被吸乳视频| 人妻丰满熟妇av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲真实伦在线观看| 最近最新中文字幕大全免费视频| 国产在线精品亚洲第一网站| 美女午夜性视频免费| 熟女人妻精品中文字幕| 成人高潮视频无遮挡免费网站| 亚洲熟妇熟女久久| 99久久成人亚洲精品观看| 日本五十路高清| 成年人黄色毛片网站| 国产精品国产高清国产av| 90打野战视频偷拍视频| 观看免费一级毛片| 国产一区二区在线观看日韩 | 久久精品91无色码中文字幕| 国产成人精品久久二区二区免费| 1024香蕉在线观看| 免费看日本二区| 免费在线观看影片大全网站| 男女视频在线观看网站免费| 少妇的逼水好多| 91麻豆精品激情在线观看国产| 国产精品亚洲av一区麻豆| 欧美一区二区精品小视频在线| 免费看十八禁软件| 老司机深夜福利视频在线观看| 十八禁网站免费在线| 色av中文字幕| 女人高潮潮喷娇喘18禁视频| 久9热在线精品视频| 在线a可以看的网站| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 色噜噜av男人的天堂激情| 精品一区二区三区四区五区乱码| 小说图片视频综合网站| 免费大片18禁| 99精品在免费线老司机午夜| 精品国产三级普通话版| 婷婷丁香在线五月| 亚洲欧美日韩无卡精品| 首页视频小说图片口味搜索| 2021天堂中文幕一二区在线观| 可以在线观看的亚洲视频| 一区二区三区激情视频| 欧美xxxx黑人xx丫x性爽| 日韩欧美在线乱码| 看免费av毛片| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 又紧又爽又黄一区二区| 成人午夜高清在线视频| 国产亚洲精品综合一区在线观看| 97超视频在线观看视频| 精品99又大又爽又粗少妇毛片 | 变态另类丝袜制服| 国产 一区 欧美 日韩| 亚洲片人在线观看| 亚洲一区高清亚洲精品| 免费在线观看影片大全网站| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| 少妇裸体淫交视频免费看高清| 亚洲av中文字字幕乱码综合| 长腿黑丝高跟| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 国产单亲对白刺激| 日韩欧美国产一区二区入口| 午夜精品久久久久久毛片777| 成人特级黄色片久久久久久久| 观看美女的网站| 欧美高清成人免费视频www| 久久久久久国产a免费观看| 欧美一级a爱片免费观看看| 国产91精品成人一区二区三区| 日韩欧美在线二视频| 欧美在线黄色| 国产成人av教育| 亚洲中文字幕一区二区三区有码在线看 | 欧美极品一区二区三区四区| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 亚洲国产欧美网| 色综合站精品国产| 亚洲 欧美一区二区三区| 波多野结衣高清作品| 精品久久久久久久久久免费视频| 国产精品免费一区二区三区在线| 国产又黄又爽又无遮挡在线| 国产精品久久久久久精品电影| 欧美日韩亚洲国产一区二区在线观看| 巨乳人妻的诱惑在线观看| 亚洲一区高清亚洲精品| 国产野战对白在线观看| 成年免费大片在线观看| 无遮挡黄片免费观看| 欧美一级毛片孕妇| 一a级毛片在线观看| 精品日产1卡2卡| 亚洲成人中文字幕在线播放| 狂野欧美激情性xxxx| 欧美日本亚洲视频在线播放| 99久久国产精品久久久| 久久精品人妻少妇| 热99re8久久精品国产| 伊人久久大香线蕉亚洲五| 可以在线观看的亚洲视频| 12—13女人毛片做爰片一| 午夜精品在线福利| 国产高清视频在线播放一区| 国产成人aa在线观看| 国产男靠女视频免费网站| 99久国产av精品| 亚洲 国产 在线| 很黄的视频免费| 国产三级中文精品| 在线观看免费午夜福利视频| 12—13女人毛片做爰片一| 午夜精品久久久久久毛片777| 午夜视频精品福利| 女警被强在线播放| 18禁观看日本| 国产成人aa在线观看| 好男人电影高清在线观看| 青草久久国产| 国产一区二区三区视频了| 男人的好看免费观看在线视频| 99视频精品全部免费 在线 | 99久久精品一区二区三区| 国产av在哪里看| 免费人成视频x8x8入口观看| 在线观看免费午夜福利视频| 999久久久国产精品视频| www.www免费av| 日本黄色视频三级网站网址| 久久久久国内视频| 亚洲午夜精品一区,二区,三区| 亚洲av免费在线观看| 国产成人啪精品午夜网站| 免费看a级黄色片| 国产成人一区二区三区免费视频网站| 国产精品98久久久久久宅男小说| 国产精华一区二区三区| 日本免费一区二区三区高清不卡| 在线免费观看不下载黄p国产 | 久久精品综合一区二区三区| 中文字幕高清在线视频| 国产精品女同一区二区软件 | 日本五十路高清| 午夜激情欧美在线| 国产男靠女视频免费网站| 色视频www国产| 啦啦啦免费观看视频1| 99热这里只有精品一区 | 日韩欧美国产在线观看| 草草在线视频免费看| svipshipincom国产片| 久久亚洲精品不卡| 夜夜躁狠狠躁天天躁| 国产av不卡久久| 国产亚洲精品av在线| 日日干狠狠操夜夜爽| 最近最新中文字幕大全电影3| 久久精品aⅴ一区二区三区四区| 精品久久久久久成人av| 日韩 欧美 亚洲 中文字幕| 性欧美人与动物交配| 伊人久久大香线蕉亚洲五| 国产精品亚洲av一区麻豆| 久久这里只有精品中国| 国产aⅴ精品一区二区三区波| 久久精品综合一区二区三区| 熟妇人妻久久中文字幕3abv| 中亚洲国语对白在线视频| 精品不卡国产一区二区三区| 国产亚洲精品久久久com| 午夜视频精品福利| 首页视频小说图片口味搜索| 精华霜和精华液先用哪个| 最新在线观看一区二区三区| 久久久久久九九精品二区国产| 操出白浆在线播放| 在线视频色国产色| 日本黄色片子视频| 欧美又色又爽又黄视频| 国产精品美女特级片免费视频播放器 | 美女午夜性视频免费| 高清在线国产一区| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 午夜a级毛片| 亚洲av成人精品一区久久| 九色成人免费人妻av| 国产成人精品久久二区二区免费| 一个人免费在线观看电影 | 色精品久久人妻99蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 午夜成年电影在线免费观看| 国产成人系列免费观看| 黄色丝袜av网址大全| 国产精品一区二区三区四区久久| 99久久国产精品久久久| 中文字幕人成人乱码亚洲影| 一本精品99久久精品77| 国产精品日韩av在线免费观看| 一区二区三区国产精品乱码| 一夜夜www| 国产精品久久久久久久电影 | 亚洲男人的天堂狠狠| 哪里可以看免费的av片| 亚洲午夜理论影院| 男女之事视频高清在线观看| 成人特级黄色片久久久久久久| 日韩有码中文字幕| 欧美av亚洲av综合av国产av| 黄色丝袜av网址大全| 色综合站精品国产| 手机成人av网站| 亚洲欧美日韩高清专用| 国产精品久久久av美女十八| 国产亚洲欧美98| 国产免费av片在线观看野外av| 国产成人av激情在线播放| av在线天堂中文字幕| 在线观看免费视频日本深夜| 99热6这里只有精品| 久久婷婷人人爽人人干人人爱| 看免费av毛片| 欧美性猛交╳xxx乱大交人| 国产精品99久久99久久久不卡|