• <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è)的思想
    久久精品国产亚洲网站| 熟妇人妻久久中文字幕3abv| 青春草国产在线视频 | 成人二区视频| 久久久国产成人免费| 欧美在线一区亚洲| 日日干狠狠操夜夜爽| 日本免费一区二区三区高清不卡| 欧美3d第一页| 狂野欧美白嫩少妇大欣赏| 亚洲精品久久国产高清桃花| 欧美在线一区亚洲| 夜夜爽天天搞| 免费搜索国产男女视频| 欧美zozozo另类| 在线观看午夜福利视频| 欧美成人a在线观看| 男女视频在线观看网站免费| 99久久无色码亚洲精品果冻| 一边摸一边抽搐一进一小说| 一级毛片我不卡| 好男人视频免费观看在线| 国产高清有码在线观看视频| av在线天堂中文字幕| 亚洲无线在线观看| 一级毛片久久久久久久久女| 天美传媒精品一区二区| 综合色av麻豆| 一级毛片电影观看 | 亚洲av中文av极速乱| 日韩一本色道免费dvd| 此物有八面人人有两片| 神马国产精品三级电影在线观看| 日本与韩国留学比较| 精品久久久久久久久亚洲| 自拍偷自拍亚洲精品老妇| 欧美又色又爽又黄视频| 国产三级在线视频| 国产亚洲5aaaaa淫片| 久久精品91蜜桃| 九草在线视频观看| 亚洲国产精品合色在线| 免费看光身美女| 色综合色国产| h日本视频在线播放| 精品不卡国产一区二区三区| 国产成人freesex在线| 少妇高潮的动态图| 尾随美女入室| 尾随美女入室| 国产探花极品一区二区| 免费av观看视频| 中国美白少妇内射xxxbb| 久久久久久久久中文| 亚洲欧美成人精品一区二区| 在线播放无遮挡| 人妻制服诱惑在线中文字幕| www.av在线官网国产| 美女脱内裤让男人舔精品视频 | 18禁黄网站禁片免费观看直播| 中文欧美无线码| 国产三级在线视频| 黄片wwwwww| 我的女老师完整版在线观看| 亚洲av免费在线观看| 男女边吃奶边做爰视频| 97在线视频观看| 在线观看美女被高潮喷水网站| 午夜亚洲福利在线播放| 毛片一级片免费看久久久久| 少妇猛男粗大的猛烈进出视频 | 欧美xxxx性猛交bbbb| 成人鲁丝片一二三区免费| 校园春色视频在线观看| 国产v大片淫在线免费观看| 一本精品99久久精品77| 成人午夜精彩视频在线观看| 色5月婷婷丁香| 波野结衣二区三区在线| 夜夜看夜夜爽夜夜摸| 天天躁夜夜躁狠狠久久av| 日韩人妻高清精品专区| 黄片无遮挡物在线观看| 国产精品久久久久久精品电影| 国产精品一区二区三区四区久久| 亚洲18禁久久av| 一级毛片久久久久久久久女| 一个人观看的视频www高清免费观看| 最近手机中文字幕大全| 热99re8久久精品国产| 如何舔出高潮| 成人综合一区亚洲| 久久亚洲国产成人精品v| av在线亚洲专区| 白带黄色成豆腐渣| av.在线天堂| 天堂中文最新版在线下载 | 免费观看精品视频网站| 国产在视频线在精品| 亚洲国产精品国产精品| 男女那种视频在线观看| 国产av麻豆久久久久久久| 国产伦理片在线播放av一区 | av黄色大香蕉| 国内久久婷婷六月综合欲色啪| 青春草视频在线免费观看| 99久久精品一区二区三区| 天堂av国产一区二区熟女人妻| 成人一区二区视频在线观看| 久久久久久久久久成人| 成年版毛片免费区| 日韩欧美 国产精品| 国产三级中文精品| 男女下面进入的视频免费午夜| 亚洲内射少妇av| 日本av手机在线免费观看| av卡一久久| 午夜激情福利司机影院| 国产v大片淫在线免费观看| 91精品国产九色| 国产成人精品一,二区 | 五月玫瑰六月丁香| 校园人妻丝袜中文字幕| 午夜福利视频1000在线观看| 国产日本99.免费观看| 如何舔出高潮| 亚洲欧美日韩高清在线视频| 成人永久免费在线观看视频| .国产精品久久| 天天躁夜夜躁狠狠久久av| 久久这里只有精品中国| 亚洲国产欧洲综合997久久,| 精品久久国产蜜桃| 狠狠狠狠99中文字幕| 国产精品嫩草影院av在线观看| 国产亚洲欧美98| 99久久成人亚洲精品观看| 深夜a级毛片| 边亲边吃奶的免费视频| 婷婷色综合大香蕉| 狂野欧美激情性xxxx在线观看| 国产老妇伦熟女老妇高清| 欧美丝袜亚洲另类| 免费在线观看成人毛片| 国产乱人偷精品视频| 久久这里有精品视频免费| 夫妻性生交免费视频一级片| 中文字幕av成人在线电影| 我的老师免费观看完整版| 国产一区二区三区av在线 | 免费看av在线观看网站| 最近最新中文字幕大全电影3| 久久久久久久久中文| 在线观看免费视频日本深夜| 日韩av不卡免费在线播放| 欧美日韩精品成人综合77777| 美女黄网站色视频| 岛国在线免费视频观看| av女优亚洲男人天堂| 尤物成人国产欧美一区二区三区| 国产av麻豆久久久久久久| 伦理电影大哥的女人| 久久精品人妻少妇| 亚洲欧洲国产日韩| 男女边吃奶边做爰视频| 一区二区三区免费毛片| 成人av在线播放网站| 高清在线视频一区二区三区 | 熟女人妻精品中文字幕| 欧美成人免费av一区二区三区| 黄色欧美视频在线观看| 一个人看的www免费观看视频| 国内精品美女久久久久久| 国产一区二区在线观看日韩| 国产高清视频在线观看网站| 嫩草影院新地址| 国产美女午夜福利| 亚洲在线自拍视频| 伦精品一区二区三区| 国产免费一级a男人的天堂| 国产一区二区在线观看日韩| 色哟哟·www| 97热精品久久久久久| 女人被狂操c到高潮| 美女被艹到高潮喷水动态| 亚洲无线观看免费| 身体一侧抽搐| av免费观看日本| 天堂√8在线中文| 一个人观看的视频www高清免费观看| 一级黄片播放器| 久久久a久久爽久久v久久| 久久久久久久久久成人| 久久欧美精品欧美久久欧美| 国产精品免费一区二区三区在线| 亚洲av熟女| a级毛色黄片| 亚洲三级黄色毛片| 婷婷亚洲欧美| 国产亚洲5aaaaa淫片| 国产男人的电影天堂91| av卡一久久| 一级av片app| 国产成人午夜福利电影在线观看| 在线天堂最新版资源| 亚洲欧美日韩无卡精品| 91久久精品国产一区二区三区| 午夜福利视频1000在线观看| 特大巨黑吊av在线直播| 大又大粗又爽又黄少妇毛片口| 一个人看的www免费观看视频| 国产熟女欧美一区二区| 欧美日韩精品成人综合77777| 欧美性感艳星| 最近的中文字幕免费完整| 国产探花在线观看一区二区| 永久网站在线| 成人国产麻豆网| 搡老妇女老女人老熟妇| 国产美女午夜福利| 日韩视频在线欧美| 亚洲中文字幕日韩| 国产一区二区在线观看日韩| 综合色丁香网| 国产69精品久久久久777片| 联通29元200g的流量卡| 日本成人三级电影网站| 色噜噜av男人的天堂激情| 男女边吃奶边做爰视频| 春色校园在线视频观看| 亚洲天堂国产精品一区在线| 黄色日韩在线| 国产真实伦视频高清在线观看| 噜噜噜噜噜久久久久久91| 亚洲国产精品成人久久小说 | 看非洲黑人一级黄片| 日韩强制内射视频| 最近2019中文字幕mv第一页| 欧美最黄视频在线播放免费| av.在线天堂| 国产精品,欧美在线| 亚洲七黄色美女视频| 桃色一区二区三区在线观看| 99精品在免费线老司机午夜| 久久久久久久久久久免费av| 日韩中字成人| 欧美在线一区亚洲| 国产一级毛片在线| 中文字幕久久专区| 不卡一级毛片| 亚洲国产欧洲综合997久久,| 久久久色成人| 国产日韩欧美在线精品| 天堂中文最新版在线下载 | 日韩av不卡免费在线播放| 国产精品人妻久久久影院| 亚洲国产精品成人久久小说 | 12—13女人毛片做爰片一| 亚洲美女视频黄频| 久久韩国三级中文字幕| 亚洲av电影不卡..在线观看| 日韩欧美在线乱码| 欧美+亚洲+日韩+国产| 天堂√8在线中文| 日韩亚洲欧美综合| 国产 一区 欧美 日韩| 久久鲁丝午夜福利片| 神马国产精品三级电影在线观看| 国产亚洲91精品色在线| 日韩亚洲欧美综合| av专区在线播放| 永久网站在线| 国产男人的电影天堂91| 1024手机看黄色片| 日韩视频在线欧美| 日韩欧美 国产精品| 变态另类丝袜制服| 亚洲欧美中文字幕日韩二区| 日韩大尺度精品在线看网址| 永久网站在线| 亚洲18禁久久av| 亚洲欧美日韩高清专用| 亚洲欧美日韩高清在线视频| 国产黄片美女视频| 国产高清激情床上av| 久久韩国三级中文字幕| 99久久人妻综合| 在线观看66精品国产| 国产高清激情床上av| 99久国产av精品国产电影| 又黄又爽又刺激的免费视频.| 亚洲一区二区三区色噜噜| 天堂√8在线中文| 亚洲成人久久爱视频| 国产午夜精品论理片| 欧美日韩在线观看h| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 国产成人一区二区在线| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 热99在线观看视频| 卡戴珊不雅视频在线播放| 亚洲图色成人| 女同久久另类99精品国产91| 97超视频在线观看视频| 亚洲四区av| 久久这里有精品视频免费| 久久热精品热| 97在线视频观看| 欧美3d第一页| 国产伦一二天堂av在线观看| 国产精品久久久久久亚洲av鲁大| 国产精品爽爽va在线观看网站| 精品无人区乱码1区二区| 国产精品人妻久久久久久| 国产精品久久久久久亚洲av鲁大| 日韩强制内射视频| 女人被狂操c到高潮| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 午夜免费激情av| 欧美一级a爱片免费观看看| 一区二区三区高清视频在线| 成年女人永久免费观看视频| 久久久久网色| 啦啦啦观看免费观看视频高清| 国产精品一区二区性色av| a级毛片免费高清观看在线播放| 久久久久性生活片| av在线亚洲专区| 搡女人真爽免费视频火全软件| av卡一久久| 亚洲欧美日韩高清在线视频| 久久久久久九九精品二区国产| 舔av片在线| 免费看美女性在线毛片视频| 精品久久久久久久末码| 人妻久久中文字幕网| 婷婷精品国产亚洲av| 精品午夜福利在线看| 99热6这里只有精品| 亚洲欧美成人精品一区二区| 国产不卡一卡二| 日本欧美国产在线视频| 伦精品一区二区三区| www.色视频.com| 国产人妻一区二区三区在| 久久久久免费精品人妻一区二区| 少妇人妻一区二区三区视频| 亚洲va在线va天堂va国产| 成人一区二区视频在线观看| 黄片无遮挡物在线观看| 18禁在线无遮挡免费观看视频| 午夜精品一区二区三区免费看| 亚洲激情五月婷婷啪啪| 国产黄片视频在线免费观看| 精品久久久久久久久久久久久| 国产精品综合久久久久久久免费| 综合色丁香网| 最近最新中文字幕大全电影3| 99久久久亚洲精品蜜臀av| 青春草亚洲视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲电影在线观看av| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 久久久久国产网址| 美女黄网站色视频| 秋霞在线观看毛片| 黄色日韩在线| 久久久精品欧美日韩精品| 国产精品久久电影中文字幕| 人妻夜夜爽99麻豆av| 久久久久久国产a免费观看| 桃色一区二区三区在线观看| 日本成人三级电影网站| 淫秽高清视频在线观看| 高清毛片免费看| 久久久久久九九精品二区国产| 一夜夜www| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 国产成人a∨麻豆精品| 青春草视频在线免费观看| 国产一区二区三区av在线 | 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 久久精品夜色国产| 中国美白少妇内射xxxbb| 搡女人真爽免费视频火全软件| 亚洲最大成人手机在线| 看黄色毛片网站| 久久中文看片网| 亚洲欧美中文字幕日韩二区| 内射极品少妇av片p| 看非洲黑人一级黄片| 一本一本综合久久| 国产免费男女视频| 久久精品国产清高在天天线| 麻豆精品久久久久久蜜桃| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 日本一本二区三区精品| 国产真实伦视频高清在线观看| 亚洲av免费高清在线观看| 91狼人影院| 国产精品国产三级国产av玫瑰| 麻豆av噜噜一区二区三区| 亚洲人成网站在线播放欧美日韩| or卡值多少钱| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱| 99riav亚洲国产免费| 亚洲性久久影院| 国内揄拍国产精品人妻在线| 热99re8久久精品国产| www.av在线官网国产| 国产精品一区二区三区四区久久| 成人亚洲精品av一区二区| 你懂的网址亚洲精品在线观看 | 最近视频中文字幕2019在线8| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 国产伦在线观看视频一区| 长腿黑丝高跟| 欧美另类亚洲清纯唯美| 寂寞人妻少妇视频99o| 黄色视频,在线免费观看| 国产精品麻豆人妻色哟哟久久 | 欧美丝袜亚洲另类| 99视频精品全部免费 在线| 搡老妇女老女人老熟妇| 成人欧美大片| 国产伦精品一区二区三区视频9| 日本撒尿小便嘘嘘汇集6| 寂寞人妻少妇视频99o| 日韩av在线大香蕉| 日韩成人伦理影院| 国产av不卡久久| 日本在线视频免费播放| 亚洲乱码一区二区免费版| 国产人妻一区二区三区在| 国产毛片a区久久久久| 久久久久久大精品| 99热这里只有是精品50| 在线观看一区二区三区| 秋霞在线观看毛片| 亚洲国产高清在线一区二区三| 岛国在线免费视频观看| 18+在线观看网站| 长腿黑丝高跟| 婷婷色av中文字幕| 中文字幕av成人在线电影| 午夜福利在线观看吧| 一个人看视频在线观看www免费| 国产成人福利小说| 国产大屁股一区二区在线视频| 在线观看美女被高潮喷水网站| 午夜精品在线福利| 亚洲国产精品合色在线| 免费大片18禁| 欧美潮喷喷水| 国产精品,欧美在线| 亚洲美女搞黄在线观看| 成人三级黄色视频| 日韩 亚洲 欧美在线| 天堂中文最新版在线下载 | 九九热线精品视视频播放| 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 久久久久性生活片| 一本久久精品| 久久久久久久久久成人| 国产高清激情床上av| 99国产精品一区二区蜜桃av| 99热6这里只有精品| 国产精品永久免费网站| 亚洲国产日韩欧美精品在线观看| 欧美最新免费一区二区三区| 亚洲成人精品中文字幕电影| 久久精品影院6| 久久午夜福利片| 成人欧美大片| 97人妻精品一区二区三区麻豆| 1024手机看黄色片| 插逼视频在线观看| 久久韩国三级中文字幕| 亚洲性久久影院| 久久久久久久午夜电影| 成年av动漫网址| 午夜a级毛片| 性欧美人与动物交配| 神马国产精品三级电影在线观看| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av| 69av精品久久久久久| 不卡视频在线观看欧美| 久久综合国产亚洲精品| 看黄色毛片网站| 成人特级av手机在线观看| 免费观看精品视频网站| 99热只有精品国产| 国产av一区在线观看免费| 国产精品无大码| 久久人人爽人人片av| 欧美+亚洲+日韩+国产| 亚洲欧美中文字幕日韩二区| 亚洲欧美精品自产自拍| 免费观看的影片在线观看| 在线观看免费视频日本深夜| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区三区在线臀色熟女| 亚洲丝袜综合中文字幕| 国产男人的电影天堂91| 天堂网av新在线| 国产成年人精品一区二区| 观看免费一级毛片| 日韩欧美三级三区| 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 亚洲国产精品久久男人天堂| 亚洲国产精品sss在线观看| 又爽又黄无遮挡网站| 亚洲自拍偷在线| 干丝袜人妻中文字幕| 国产午夜精品久久久久久一区二区三区| 国产中年淑女户外野战色| 国产极品天堂在线| 久久久久性生活片| 亚洲真实伦在线观看| 久久精品久久久久久久性| 一区二区三区高清视频在线| 国产淫片久久久久久久久| 亚洲国产高清在线一区二区三| 插阴视频在线观看视频| 国产老妇女一区| 国产黄色小视频在线观看| 亚洲天堂国产精品一区在线| 丰满的人妻完整版| 久久午夜亚洲精品久久| 亚洲精华国产精华液的使用体验 | 国产乱人偷精品视频| 国产成人91sexporn| 亚洲国产精品国产精品| 国产成人91sexporn| 天堂影院成人在线观看| 美女国产视频在线观看| 亚洲国产精品国产精品| 亚洲一级一片aⅴ在线观看| 丝袜喷水一区| 国产成人精品婷婷| 99热精品在线国产| 一边亲一边摸免费视频| 黄色欧美视频在线观看| 久久亚洲精品不卡| 久久久久久久久久久丰满| 丰满人妻一区二区三区视频av| 老女人水多毛片| 99热只有精品国产| 亚洲丝袜综合中文字幕| 欧美性猛交黑人性爽| 91精品一卡2卡3卡4卡| 免费大片18禁| 国产成人freesex在线| 日韩av在线大香蕉| 在线a可以看的网站| 能在线免费看毛片的网站| 亚洲中文字幕日韩| 久久九九热精品免费| 1000部很黄的大片| 能在线免费观看的黄片| 日日啪夜夜撸| 国产精品伦人一区二区| 成人国产麻豆网| 国产大屁股一区二区在线视频| av又黄又爽大尺度在线免费看 | 秋霞在线观看毛片| 久久人人精品亚洲av| 中文字幕制服av| 日韩制服骚丝袜av| 日本一本二区三区精品| 亚洲内射少妇av| 18+在线观看网站| 听说在线观看完整版免费高清| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲经典国产精华液单| 在线观看午夜福利视频| 国产精品一区二区性色av| 久久精品夜色国产| 成人三级黄色视频| 国产精品永久免费网站| 国产精品人妻久久久影院| 美女被艹到高潮喷水动态| 中文字幕熟女人妻在线| 久久久久国产网址| 久久综合国产亚洲精品| 国内精品一区二区在线观看| 一边摸一边抽搐一进一小说| 亚洲一区二区三区色噜噜| 午夜激情福利司机影院| 日韩欧美精品v在线| 亚洲国产精品久久男人天堂| 免费不卡的大黄色大毛片视频在线观看 | 直男gayav资源| 久久中文看片网| 天堂av国产一区二区熟女人妻| 久久综合国产亚洲精品| 免费av不卡在线播放| 久久午夜亚洲精品久久|