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

    Towards high-power mid-IR light source tunable from 3.8 to 4.5 μm by HBr-filled hollow-core silica fibres

    2022-03-19 09:27:12ZhiyueZhouZefengWangWeiHuangYulongCuiHaoLiMengWangXiaomingXiShoufeiGaoandYingyingWang
    Light: Science & Applications 2022年1期

    Zhiyue Zhou,Zefeng Wang?,Wei Huang,Yulong Cui,Hao Li,Meng Wang,Xiaoming Xi,Shoufei Gao and Yingying Wang

    1College of Advanced Interdisciplinary Studies,National University of Defense Technology,Changsha 410073,China

    2State Key Laboratory of Pulsed Power Laser Technology,Changsha 410073,China

    Abstract Fibre lasers operating at the mid-IR have attracted enormous interest due to the plethora of applications in defence,security,medicine,and so on.However,no continuous-wave (CW) fibre lasers beyond 4 μm based on rare-earthdoped fibres have been demonstrated thus far.Here,we report efficient mid-IR laser emission from HBr-filled silica hollow-core fibres(HCFs)for the first time.By pumping with a self-developed thulium-doped fibre amplifier seeded by several diode lasers over the range of 1940-1983 nm,narrow linewidth mid-IR emission from 3810 to 4496 nm has been achieved with a maximum laser power of about 500 mW and a slope efficiency of approximately 18%.To the best of our knowledge,the wavelength of 4496 nm with strong absorption in silica-based fibres is the longest emission wavelength from a CW fibre laser,and the span of 686 nm is also the largest tuning range achieved to date for any CW fibre laser.By further reducing the HCF transmission loss,increasing the pump power,improving the coupling efficiency,and optimizing the fibre length together with the pressure,the laser efficiency and output power are expected to increase significantly.This work opens new opportunities for broadly tunable high-power mid-IR fibre lasers,especially beyond 4 μm.

    Introduction

    Laser sources operating at the mid-infrared(mid-IR)are of considerable interest due to wide applications in biomedicine,spectroscopy,defence and manufacturing fields.Among these laser sources,fibre lasers have been favoured in recent years owing to their diffraction-limited beam quality,high conversion efficiency,long interaction length and compact system configuration1,2.However,due to the high phonon energy at approximately 1100 cm-1,silica-based fibres suffer from strong material absorption in the mid-IR range above 2.2 μm1.Extending the emission wavelength coverage of fibre lasers further into the mid-IR region,especially above 3 μm,is still an ongoing challenge1.For longer wavelengths,we usually use soft glass fibres of different host materials,including fluoride,telluride and chalcogenide glasses that have higher transparency in the mid-IR.A variety of transitions offered by rare-earth-doped soft glasses have been demonstrated for continuous-wave (CW) mid-IR fibre laser emission3-14.Figure 1 shows a summary of state-ofart CW mid-IR fibre output laser sources in terms of output power and slope efficiency as a function of lasing wavelengths.All the reported CW mid-IR lasers based on rare-earth-doped fibres operate below 4 μm thus far,and the output power decreases exponentially with the increasing emission wavelength,which is mainly due to the increase in the quantum defect between the pump and laser photon energies at longer wavelengths1.The solid lines with symbols show the tunable mid-IR fibre lasers.The tunable range of the fibre lasers is limited by the gain bandwidth of the doped ions.For mid-IR fibre lasers,the largest tunable range of 573 nm (2807-3380 nm) was obtained based on dysprosium-doped ZBLAN fibres with a maximum output power of 170 mW10.Recently,a pulsed laser operating at 5.38 μm was reported in terbium doped chalcogenide glass fibre15.Due to the high loss of the fibre,the laser threshold was just reached,and only emission spectrum was recorded without power values.In terms of mechanical strength,chemical durability,thermal conductivity,and optical nonlinearity,the performance of systems in the mid-IR based on these soft glasses is well behind that at shorter wavelengths using silica1.

    Fig.1 Summary of the state-of-art CW mid-IR fibre lasers in terms of output power and lasing wavelengths.

    Compared to silica-based solid core fibre lasers with difficulties in extending emission wavelength into the mid-IR spectral region,silica-based HCF gas lasers have opened new avenues for mid-IR light source.The mid-IR absorption loss of silica host materials in HCF can be reduced substantially owing to the small overlap of the core modes with the silica material,and the field in the silica material is at least an order of magnitude smaller than the peak field in the core,giving a higher damage threshold.The fast development of silica-based HCFs with low attenuation at mid-IR wavelengths16-19is a significant motivation to develop mid-IR gas-filled HCF lasers which combine the advantages of both traditional fibre lasers and gas lasers20-22.In addition,the HCF can provide an ideal circumstance for the interaction of gas and light with a much longer interaction distance than traditional gas cells23-25.To date,several mid-IR lasers based on gas-filled HCFs have been demonstrated,which can be divided into two categories according to the operating mechanism,namely stimulated Raman scattering (usually using H2and CH426-29) and population inversion realized by intrinsic absorption of gas molecules(usually using C2H2,CO,N2O,HCN and CO230-37).The laser threshold based on stimulated Raman scattering is about 5-6 orders of magnitude higher than that based on population inversion.Therefore,all the reported mid-IR Raman lasers based on gas-filled HCFs are pulsed26-29,while lasers based on population inversion are easier to realize CW mid-IR output32-37.Recently,a pulsed 4.6 μm laser based on N2O-filled Kagome HCFs was demonstrated30.However,limited by the transmission loss at the mid-IR,as well as the linewidth of the pump optical parametric oscillator,the laser slope efficiency is <3%.We recently reported a CW fibre laser at 4.3 μm by CO2-filled HCFs37,however,the possible output laser wavelength range is small (4.28-4.42 μm) due to the transition properties of CO2molecules.In many applications,such as multispecies trace gas detection,IR countermeasures and free-space optical communications,a mid-IR fibre laser covering much wider wavelength range over 3-5 μm is required.The HBr gas medium has sparse rotationalvibrational energy level and relaxation characteristics,enabling the achievement of broadly tunable mid-IR lasers.Up to now,although traditional HBr gas lasers operating at mid-IR spectral region have been demonstrated in gas cells38-42.The effective interaction length of the gas media and the laser beam is very short,and the systems are usually bulky and cumbersome,seriously limiting the applications of these lasers.Therefore,the silica-based HCF with low mid-IR loss and long interaction length is envisioned to provide an ideal waveguide for the HBr mid-IR emission lines.

    Here,we report the first CW HBr laser in a silica-based HCF.To demonstrate the widely wavelength tunable characteristics,a narrow linewidth 2 μm thulium-doped fibre amplifier (TDFA) seeded by a group of fine-tunable diode lasers is used to pump a 5 m-long anti-resonant HCF filled with low-pressure HBr gas.The wavelengths of the selected diode lasers cover the absorption lines of R(11),R(7),R(5),R(3),R(2) and R(0) for both isotopes H79Br and H81Br.A total of eleven narrow linewidth laser transitions,with five R-branch and six P-branch transitions,covering 3810-4496 nm are observed,which is the broadest tuning range and the 4496 nm output with strong absorption in silica-based fibres is the longestwavelength among CW fibre lasers to the best of our knowledge.When the HBr pressure is 5 mbar,the maximum laser output of 500 mW at approximately 4.2μm is achieved when pumped by the R(3)absorption line,with a slope efficiency of 18%.In addition,the output spectral components can be efficiently controlled under appropriate conditions and the laser exhibits excellent beam quality performance in the 4 μm CW region,with a measured beam quality factorM2of approximately 1.2.This work paves the way for compact widely tunable highpower mid-IR fibre laser sources in the future.

    Results

    Theoretical analysis of the energy level transitions of the HBr molecule

    The HBr molecule occurs in two isotopes,H79Br and H81Br,with almost equal natural abundances of 50.678 and 49.306%43,respectively.The corresponding isotope energy level mismatch between the two isotopes is approximately 50 GHz39.In addition,HBr is a diatomic molecule and consequently has only one vibrational normal mode,displayed in the inset of Fig.2a.The vibrational normal mode can be expressed by vibrational quantum numberv,which takes integer values from zero upwards.In each vibrational state,there is a series of rotational states due to molecular rotation,which are expressed by the rotational quantum numberJ.Owing to the small moments of inertia and strong molecular bonds,the HBr molecules have characteristically sparse rotational-vibrational energy levels38,39.The energy diagram of HBr energy levels for 4 μm emission is shown in Fig.2a.The HBr molecules can be excited from the ground state (v=0 vibrational state) to the upper level(v=2 vibrational state) through rotational-vibrational transitions.Due to the large vibrational level spacing,the population of the lower laser level(v=1 vibrational state)at room temperature is nearly zero40.Then,the HBr molecules can leave the upper level through radiation transition to the essentially emptyv=1 vibrational state,from which lasing occurs according to the selection rule ΔJ=+1 (or ΔJ=-1),referred to as R-branch (or P-branch) transitions.The population of thev=1 vibrational state is transferred back to the vibrational ground state via vibrational relaxation (non-radiative transition)induced by intermolecular collisions.

    Fig.2 Characteristics of the radiative transition of the HBr molecule.

    From the HITRAN database43,Fig.2b illustrates the first overtone absorption band of H79Br molecules at approximately 2 μm from the ground vibrational statev=0 to the excited vibrational statev=2.The number at each line represents the corresponding absorption lines of the R- or P-branch.The absorption intensity of different lines is directly proportional to the population of the rotational levels,which is governed by the Boltzmann distribution and the possibility of degeneracy44,indicating that the population rises to a maximum and then diminishes,as shown in Fig.2b.The corresponding emission band is illustrated in Fig.2c,which is similar to the absorption band.Table 1 summarizes the specific absorption and corresponding emission wavelengths of the R-branch.Owing to the increasing spacing between adjacent rotational states,a small tunable range of the pump,for example,from 1939 to 1983 nm,can generate an ultrabroad range of lasing wavelengths from 3760 to 4540 nm.In this experiment,the R(11),R(7),R(5),R(3),R(2) and R(0) absorption lines are chosen as pump wavelengths to achieve a tunable output lasing spectrum.Indeed,due to the strong absorption in the HCF,the relatively weak pump lines can also be absorbed by increasing the fibre length,so the potential lasing wavelength range can be broader.In addition,rotational relaxation usually occurs with a pulsed pump source with a long pulse duration due to inelastic collision between molecules,which will cause the population of the upper level to redistribute to other rotational levels and can also be used to achieve broadly tunable HBr lasers in the pulsed region41.

    Table.1 Calculated R-branch absorption lines and corresponding lasing lines of H79Br.

    Experimental layout and properties of the used HCFs and narrow linewidth pump source

    The experimental layout is shown in Fig.3a.High gain resulting from tight confinement of the pump light together with the active gas in the HCF permits to operate the laser in a single pass-configuration,which is actually the process of amplified spontaneous emission (ASE),without the use of any external resonator structure.The setup centres around a 5 m length of an in-house fabricated anti-resonant HCF surrounded by a single-ring 6-tube-cladding fused silica-based HCF structure,and a scanning electron micrograph (SEM) of the HCF crosssection is shown in Fig.3b.The core diameter of the HCF is approximately 80 μm,the average inner diameter of the cladding capillary tubes is approximately 38 μm,and the average wall thickness is approximately 700 nm.This design can guide light over a large wavelength range from 1.9 to 4.5 μm,covering both the pump and corresponding laser wavelengths.Figure 3c shows the simulated HCF confinement loss and material absorption loss using a finite-element solver with a perfectly matched layer45,plotted in red and green curves respectively while the blue curve implies the total simulated loss including both two losses.The black dots illustrate the transmission loss value from 0.2 to 0.4 dB/m obtained using the cut-back measurements (see “Methods”),which are much higher than the simulated loss in shorter wavelength below 3.9 μm while close to the simulated loss in longer wavelength.The higher measured loss may be contributed by the structural nonuniformity in the longitudinal and transverse directions.The variation tendency of the measured HCF loss is consistent with that of simulated loss.The pump light is collimated passing through the first lens and aligned by two parallel silver-coated mirrors.Then the pump light is focused into the input end of the HCF through the second lens with a coupling efficiency of~60%.At the output,an IR bandpass filter (>80% transmission at 4 μm and <0.1% transmission at 2 μm) is used to separate the generated mid-IR laser from the unabsorbed fraction of pump light that is transmitted by the HCF filled with HBr gas.The filtered output can be monitored using a thermal power metre or an optical spectrum analyzer (OSA).As displayed in Fig.3a,the pump system is a TDFA with a two-stage amplifier arrangement seeded by one of six diode lasers corresponding to the R(11),R(7),R(5),R(3),R(2) and R(0)absorption lines of HBr which are marked in blue dashed line box.All diode lasers are driven by four voltage pins,namely Vtec (the temperature-controlling voltage),Vbias(the bias voltage),Vcc (the supply voltage) and Gnd (the ground connection).Each diode laser can be precisely tuned in a few nanometres by adjusting the Vbias and Vtec,covering the absorption line of HBr.Figure 3d shows the measured spectrum of the pump system seeded by different absorption line wavelengths,which has a very weak ASE background.Even with increasing output power of the pump system seeded by the R(5)absorption line to the maximum output power of 8 W,the ASE background is not obviously enhanced,illustrating that the output power is mainly concentrated at the central wavelength.Figure 3e shows the wavelength of the R(5)absorption line-seeded pump system varying with the Vtec when the Vbias is 1.2 V.The wavelength can be accurately tuned from 1957.4 to 1959.8 nm with a good linear relationship of 0.85 nm/V,covering the R(5) absorption line of both the H79Br isotope at 1957.8 nm and H81Br isotope at 1958.1 nm.The pump system seeded by the other wavelength diode lasers has similar wavelength tuning results (see the Supplementary Information).

    Fig.3 Experimental layout and measured and calculated properties of the HCF and the homemade pump source.

    For a gas laser,the linewidth and wavelength stability of the pump source are very important due to the narrow absorption linewidth of molecules,which is usually several hundred MHz at low gas pressure21,determined by collisional broadening and Doppler broadening.By precisely tuning the pump wavelength across the R(3)absorption line based on wavelength tuning results and then measuring the power transmitted by the HCF filled with different HBr gas pressures at each wavelength,the absorption linewidth can be measured (see “Methods”),as shown in Fig.3f(the absorption linewidth of the other absorption lines is shown in the Supplementary Information).The discrete points are the measured data,and the curve is the corresponding smooth fitting curve.The three curves correspond to pressures of 2.2,4.8 and 6.5 mbar with 3.6 pm(270 MHz),4.8 pm(360 MHz)and 5.6 pm (420 MHz) absorption linewidths (full-width at half-maximum,FWHM),respectively.The spectral linewidth of the pump system is measured using a scanning Fabry-Perot (F-P) interferometer46(setup shown in the Supplementary Information).When seeded by 1971 nm diode laser,the measured pump linewidth (FWHM) Δνis approximately 23 MHz,which is much narrower than the absorption linewidth of HBr.As long as the pressure is high enough,the pump power can be efficiently absorbed by HBr gas.Additionally,the wavelength stability of the pump system can be measured by monitoring the fluctuation of the residual pump power at the output end when the wavelength is tuned to the centre of the R(3) absorption line,where the transmitted pump power should be the minimum,as shown in Fig.3f.The results for an approximately 30 min duration are shown in Fig.3g.The residual pump power remains stable at approximately 5 mW,showing a very good wavelength stability of <1 pm,which is estimated from Fig.3f,g.

    Measured widely step-tunable mid-IR spectra

    In the experiment,the measured mid-IR spectra are closely related to the HBr pressure and pump power.When pumped by any absorption line wavelength,both R-branch and P-branch laser lines should have been observed from the same upper level according to the Δj=±1 transitions,as explained in relation to Fig.2a.As the R-branch and P-branch laser lines both share a common upper level,they compete with each other32.However,the intermolecular collisions are enhanced with increasing HBr pressure,leading to the rotational relaxation.Then the accumulated population of the upper level will transfer to other rotational energy level,causing other laser line emissions,which will be explained in the next section.Figure 4a shows the evolution of the laser peak intensity of the R(3)and P(5)laser transitions with respect to the incident pump power when the H81Br isotope gas is pumped by the 1966 nm R(3) absorption line at lower pressure of 0.9 mbar.The concrete peak intensity ratio of R(3) to P(5) is displayed in Fig.4b.The peak intensity ratio of R(3) to P(5) increases from 0 to 1.8 with increasing incident pump power.Due to the larger Einstein A coefficient of the P-branch transition32,a single P-branch transition line might be observed first.Specifically,with increasing pump power,the population in the lower lasing level of the P(5) transition accumulates,causing a decrease in population inversion and reduced gain until it is exceeded by the gain provided to the competing line R(3) with a still empty lower lasing level,allowing this second line R(3)to lase as well.However,for a higher pressure of approximately 5 mbar which can be regarded as the transition pressure,the maximum 8 W power scaling of the pump system is not high enough to saturate the P-branch transition and the intermolecular collisions are not enhanced enough to cause the rotational relaxation.In this case,a pure spectrum with only the single laser line P(5) is observed while the R(3) laser transition has vanished owing to the lower pump power,as shown in Fig.4c.By choosing appropriate HBr pressure that cannot give rise to rotational relaxation and pump power that is not high enough to allow R-branch transition to lase,the output spectral components can be efficiently controlled with a pure P-branch transition spectrum.Figure 4d shows all measured broadband steptunable mid-IR emission spectra pumped with the six absorption line wavelengths in turn.A total of eleven laser transitions with five R-branch transitions and six P-branch transitions covering 3810-4496 nm are individually observed.The R(0) laser transition,which was expected at 4025 nm when pumped by the 1983 nm R(0)absorption line,is conspicuously absent due to the smallest emission cross section43.

    Fig.4 Measured widely step-tunable mid-IR emission spectra.

    Energy level relaxation analysis

    Relaxation caused by molecular collisions is a common process in gas lasers,leading to energy transfer between different states,including energy transfer between rotational states within the same vibrational state (R-R relaxation) and energy transfer between vibrational states(V-V relaxation)47.As the energy gap for R-R relaxation is smaller than that for V-V relaxation,the rate of R-R relaxation is two orders of magnitude larger than the rate of V-V relaxation48.The net population density transfer rate in R-R relaxation can be given by47:

    wherekBis the Boltzmann constant andEjis the state energy of theJenergy level.This formula also reveals that R-R relaxation will lead to a Boltzmann distribution for each rotational state of the upper vibrational state to achieve equilibrium,and the larger the density of possible collision partners is,the faster the equilibrium is reached.For the CW pump,there is sufficient time for the population to reach equilibrium.Therefore,regardless of which rotational state is pumped,before the signal laser is generated,the population in the upper vibrational state nearly obeys the Boltzmann distribution for each rotational state.Figure 5a shows the Boltzmann distribution in the upper vibrational state at 293 K (see Supplementary Information).The population in theJ=3 rotational state is dominant.Because the emission cross-sections of each P-branch emission line are relatively close,the gain of the signal emission line is mainly dependent on the population of the upper state.Therefore,in theory,the P(4) emission line is generated first under high pressure where the R-R relaxation is strong.Interestingly,in our experiment,the exact laser emission caused by R-R relaxation is related to the type of isotope.For H79Br,the P(5) emission line is dominant in the output spectra,while for H81Br,the P(4)emission line is dominant.The mechanism of the influence of the isotope on relaxation remains to be studied.

    Figure 5b,c plots the output spectra changing with the gas pressure when the pump wavelength is tuned to the R(7)absorption line.The R(7)and P(9)transitions occur at a low gas pressure of 1.3 mbar,as similarly explained in relation to Fig.4a.For a higher pressure of 5.7 mbar,disappearance of the R(7) signal occurs because the threshold of R(7) exceeds the maximum output power level of the pump system used in the present experiment,as the threshold increases with increasing pressure.With a further increase in gas pressure to 8.2 mbar,the P(4)or P(5) transition (corresponding to the H81Br or H79Br isotope,respectively) starts to occur owing to R-R relaxation.In addition,for R(2) or R(3) pumping,P(4)or P(5) corresponds not only to the targeted laser transition line but also to a relaxation transition line.In the experiment,only the corresponding targeted laser transition line without a relaxation transition line is observed under R(2) or R(3) pumping by changing the incident power.Figure 5e,f plots the output spectra changing with the incident pump power at a gas pressure of 8.2 mbar.With increasing incident pump power,the intensity of the P(4)or P(5)line gradually decreases,while the intensity of the P(9) line gradually increases.This is because for the P(9) transition,the upper state is the directly pumped state,but for the P(4)or P(5)transition,the population of the upper state comes from R-R relaxation (see Fig.5d).Thus,the pump rate for the P(4) or P(5) transition is mainly determined by the R-R relaxation rate,while the pump rate for the P(9)transition is mainly determined by the pump power.Although under small-signal conditions,the gain of the P(4)or P(5)transition is larger than that of the P(9) transition,the gain of the P(4) or P(5) transition more easily reaches saturation due to the unchanged R-R relaxation rate,which inhibits the growth of the P(4)or P(5) laser power.However,for the P(9) transition,the stimulated emission is strengthened with increasing pump power.Therefore,in the competition between the P(9)transition and the P(4) or P(5) transition,the P(9) transition dominates with increasing incident power.

    Fig.5 Energy level relaxation characteristics.

    Mid-IR laser power output properties

    Since the two isotopes,H79Br and H81Br,gas molecules filled in the HCF with various pressures can be pumped by different absorption line wavelengths,a series of laser output power characteristics are measured.A simplified theoretical model is built to qualitatively analyze the laser generation process (see the Supplementary Information)and compare with the experimentally obtained results.Figure 6a illustrates the measured maximum output laser power of the individual 4 μm wavelength when H79Br and H81Br gas molecules are pumped.Due to the relatively low gain of the R-branch transition at approximately 4-6 mbar pressure (illustrated in Fig.4c),only the P-branch transition laser occurs.In contrast,at approximately 0.9-1.3 mbar with higher incident pump power,the R-branch transition dominates over the whole output spectra,as explained in relation to Fig.4a.Since the wavelength tuning range of the 1966 nm diode laser can only cover the R(3) absorption line of H81Br isotope molecules and the wavelength tuning range of the 1971 nm diode laser can only cover the R(2) absorption line of H79Br isotope molecules,R(11),R(7),R(5),R(2),P(2),P(4),P(7),P(9) and P(13) lasing transitions are measured when H79Br isotope molecules are pumped,while R(11),R(7),R(5),R(3),P(2),P(5),P(7),P(9) and P(13) lasing transitions are measured when H81Br isotope molecules are pumped.In addition,owing to the smallest emission cross-section,the R(0) lasing transition is absent.The maximum output of 500 mW is realized for the P(5)lasing transition when pumped by the R(3) absorption line under 5 mbar pressure.At lower pressure approximately 0.9-1.3 mbar,the pump power can just be partially absorbed with lower gain,which causes the measured maximum output laser power of R-branch transition is lower than that of P-branch transition.The measured output laser power at different HBr gas pressures as a function of the absorbed pump power when pumped by the R(3) absorption line is plotted in Fig.6b (the output laser power when pumped by the other absorption lines is shown in the Supplementary Information).The 4 μm laser power increases with the absorbed pump power without saturation,reaching a maximum value of approximately 500 mW at 5 mbar.At a lower pressure of 0.9 mbar,the pump power is only partially absorbed,leading to a lower output.In addition,the lasing threshold,defined as the minimum absorbed pump power necessary to observe a 4 μm laser output,increases with the HBr pressure due to the enhanced intermolecular collisions.The inset in Fig.6b illustrates the 1 h stability of the laser.The decline trend of the output laser is mainly owing to the ambient air leaking into the gas cells and HCF with time,decreasing the HBr gas purity and the output laser power.Especially for lower HBr gas pressure in the HCF,the leaky air will severely impair the output laser performance.Furthermore,Fig.6c compares the measured and simulated maximum output laser power and residual pump power as a function of HBr pressure when pumped by the R(3)absorption line for an incident pump power of approximately 7.5 W.With increasing gas pressure,the increased molecular density in the HCF results in higher pump power absorption and gain,giving a higher laser output and a lower residual pump power.The output laser power reaches a maximum at an optimum pressure.However,beyond the optimum pressure,due to enhanced intermolecular collisions,the lifetime of the laser upper level is shortened,leading to reduced gain and low output power,although almost all pump power is absorbed.The differences between measured and simulated results occur because the actual pressure in HCFs remains imprecise and may be lower than the pressure showed in the vacuum gauge.Thus,higher molecular density in simulation model is considered compared with experimental data.The measured and simulated 4 μm laser power pumped by the R(3) absorption line as a function of incident pump power and absorbed pump power at the optimum pressure of 5 mbar is plotted in Fig.6d.The measured output laser power increases almost linearly with the incident pump power and absorbed pump power beyond the threshold,reaching a maximum of 500 mW with slope efficiencies of 8.8 and 18% with respect to the incident pump power and absorbed pump power,respectively.Largely due to the discrepancy of pressures,the theoretical simulations are found to be in general agreement with the experimentally obtained data.

    Fig.6 Measured mid-IR laser power properties.

    Measurement of the mid-IR laser linewidth and beam quality factor M2

    Both the linewidth and beam quality of the mid-IR laser are measured when the pump wavelength is tuned to the R(7) absorption line under 5.7 mbar pressure.Under such conditions,only the 4335 nm P(9) laser line is emitted,while the 3875 nm R(7) laser line is totally suppressed,as mentioned previously.Specifically,the linewidth of the molecular transition is determined by three broadening processes,namely,collisional broadening,Doppler broadening and negligible natural broadening.The total line shape is a Voigt profile,a convolution of a collisional line shape and a Doppler line shape.The line shape of the mid-IR laser is measured using a scanning F-P interferometer with a free spectral range (FSR) of 1.5 GHz (see “Methods”,setup shown in the Supplementary Information) and is shown in Fig.7a.The time interval between the two peaks ΔTis approximately 12.86 ms,and the FWHM of one of the peaks Δtis approximately 0.48 ms.Therefore,the linewidth of the mid-IR laser can be calculated by the formula(FSR/ΔT)×Δt.The calculated linewidth is approximately 56 MHz,which is spectrally narrow without additional linewidthlimiting elements owing to the nature of gas molecule transitions.In addition,the output laser beam quality is characterized in terms ofM2,which is determined by measuring the beam profiles of the 4 μm laser output at different longitudinal positions (see “Methods”).As illustrated in Fig.7b,the output laser beam is observed to be near diffraction limited with anM2of 1.19±0.12.The inset shows the mode profile of the mid-IR laser beam at the waist position.

    Fig.7 Measured linewidth and beam quality of the mid-IR laser.

    Discussion

    In summary,an optically pumped CW 4 μm HBr gas laser has been demonstrated in HCFs for the first time.By using a self-developed precisely tunable narrow linewidth TDFA as the pump source,eleven laser transitions covering 3810-4496 nm are individually observed in the HBrfilled HCF,with the longest emission wavelength and largest tuning range.The maximum output power of the mid-IR laser is 500 mW,with a slope efficiency of 18%with respect to the absorbed pump power when the H81Br isotope gas is pumped by the R(3)absorption line under 5 mbar pressure.The output mid-IR laser exhibits a narrow linewidth of approximately 56 MHz and a near diffractionlimited beam quality with anM2of approximately 1.2.By further optimizing the fibre length and gas pressure,the efficiency and output power of the 4 μm laser are expected to be scaled significantly.In addition,widely tunable HBr lasers could be obtained by using a tunable 2 μm fibre laser covering many absorption lines.

    Looking forward,several strategies are available to further improve the performance of mid-IR gas-filled HCF lasers.Since the pump light is usually coupled into HCFs by the method of spatial optical path coupling,the coupling efficiency is unstable and easily influenced by the external environment.One of the major development directions is the employment of all-fibre structure coupling with low loss between HCFs and solid-core fibres to replace the gas cell in the experimental setup and it includes directly fusion splice49,fibre tapered technology50,reverse-tapering method51and so on.Achieving high power output is another important direction.To date,the highest power reported in gas-filled HCF lasers is only at the watt level.Several key issues need to be resolved for higher power output,mainly including a suitable narrow linewidth high-power pump source and low-loss coupling of a high-power pump laser.Obtaining more abundant laser wavelengths is also an important direction in the future.Compared with solid-core rareearth-doped fibre lasers,gas gain media are more convenient to be replaced in gas fibre lasers,and there are more choices.If the HCF transmission bands are properly designed,then with suitable gases and pump sources,various laser wavelengths can be obtained,especially in the mid-IR region,which is not easy to achieve with traditional fibre lasers.The soft glass has been used to manufacture HCFs52,and far-IR gas-filled HCF lasers are expected to be realized.Additionally,gas-filled HCF lasers also have certain advantages in realizing laser output in the visible and ultraviolet bands.Especially in the ultraviolet band,the photon darkening effect is much weaker than that of solid-core fibres.The choice of gain media in the visible and ultraviolet bands is very rich,including common inert gases,various chlorides and metal vapours.

    If the all-fibre coupling structure is employed to develop robust gas-filled HCF lasers with improved power level in the future,such mid-IR light source can provide more practical applications in lots of fields,such as remote monitoring of polluted gases or trace gases,free space communication in the atmosphere,infrared directional countermeasures,far-infrared frequency conversion,and so on.

    Methods

    HCF transmission loss measurement using cut-back methods

    Cut-back methods are employed to experimentally measure the fibre transmission loss.Usually,a broadband light source is used to couple the light into the input end of the fibre,and the transmitted light at the output end of the fibre is recorded by an OSA before the cut-back.Then,a part of the fibre as far away from the fibre input end is cut off to avoid the effects of high order modes and the output spectrum is recorded again.Notably,the input to the fibre should be kept the same at all times and the fibre should be kept in a large loop(the loop diameter is always larger than 1 m in the experiment) to prevent the bending loss.In the experiment,we use an optical parametric oscillator (the wavelength can be individually tuned from 3 to 4.1 μm)to replace the broadband light source.A power metre is also used to replace the OSA to measure the output power at several wavelengths.The fibre attenuation is expressed as:

    in whichLis the cut HCF length.P1(λ) is the output power of the original long HCF before cut-back,whileP2(λ) is the output power of the remaining HCF after cut-back.

    Absorption line measurement

    Each seed diode laser source of the pump system has four pins(shown in the experimental layout of Fig.3a),in which Vcc is the fixed supply voltage (usually 5 V),Gnd represents the ground connection,Vtec is the temperature-controlling voltage (from 0.1 to 3 V),and Vbias is the bias voltage (from 0 to 1.4 V).Specifically,except that the wavelength tuning range of the 1966 nm diode laser can only cover the R(3) absorption line of the H81Br isotope molecule and the wavelength tuning range of the 1971 nm diode laser can only cover the R(2)absorption line of the H79Br isotope molecule,the wavelength tuning range of other four diode lasers can cover the corresponding absorption lines of both the H79Br and H81Br molecules.The output centre wavelength can be adjusted by Vtec and Vbias,and in the experiment,we usually set Vbias to 1.2 V.Figure 3e shows the centre wavelength of the R(5) absorption line as a function of Vtec with a good linear relationship.By adjusting the Vtec step by step to tune the pump wavelength across the absorption line and then measuring the output power at each single wavelength,the absorption linewidth can be measured and normalized.

    Mid-IR laser power and spectra measurement

    As the experimental setup in Fig.3a shows,the output 4 μm laser and the residual pump transmit through the output window (WG31050,Thorlabs),which has a transmission of 87%at both the 2 and 4 μm bands.An IR bandpass filter (FB4250-500 or FB4000-500,Thorlabs) is set at the output end of the system to filter out the residual pump.A thermal power metre or a spectrometer(OSA-207C,Thorlabs) is used to measure the mid-IR laser power or the spectra.

    Laser linewidth measurement

    The mid-IR laser beam is aligned with the scanning F-P interferometer installed in the standard optical adjusting frame through two mirrors.Then,a lens is placed in front of the interferometer so that the beam waist is in the centre of the interferometer.The scanning voltage provided by the control box is connected to the oscilloscope while driving the piezoelectric ceramic (PZT) in the interferometer.The signal obtained by the photodiode detector in the interferometer is also connected to the oscilloscope after being amplified by the control box(setup shown in the Supplementary Information).

    Beam quality factor M2 measurement

    TheM2reported here is measured using two planoconvex lenses of 40 and 100 mm focal lengths.The output end of the HCF is placed at the focal point of the first lens,collimating the mid-IR laser beam.The second lens is employed to focus the collimated beam.Then,a CCD camera installed on a translation stage can scan the beam profiles at different positions along the propagation axis.A criterion of 1/e2of maximal intensity is applied to define the beam size at all measurement points and the value ofM2can be calculated.

    Acknowledgements

    This work is supported by the Outstanding Youth Science Fund Project of Hunan Province Natural Science Foundation (2019JJ20023),National Natural Science Foundation of China (NSFC) (11974427,12004431),State Key Laboratory of Pulsed Power Laser Technology(SKL-2020-ZR05,SKL-2021-ZR01)and Postgraduate Scientific Research Innovation Project of Hunan Province(CX20190026,CX20200047).

    Author details

    1College of Advanced Interdisciplinary Studies,National University of Defense Technology,Changsha 410073,China.2State Key Laboratory of Pulsed Power Laser Technology,Changsha 410073,China.3Hunan Provincial Key Laboratory of High Energy Laser Technology,Changsha 410073,China.4Institute of Photonics Technology,Jinan University,Guangzhou 511443,China

    Author contributions

    Z.W.conceived the idea of the 4 μm laser in HBr-filled HCF.Z.Z.assembled the setup and performed the measurements.Y.C.assisted in measuring the absorption line and output spectra.W.H.and H.L.assisted in measuring beam quality factorM2.W.H.carried out the theoretical analysis of relaxation process and gave the simplified simulation model of the laser generation process.M.W.and X.X.contributed to the interpretation of the experimental results.S.G.and Y.W.designed and fabricated the HCF.S.G.and Y.W.simulated the HCF transmission loss.Z.Z.,Y.C.and W.H.prepared the manuscript.Z.W.supervised this work and led the scientific collaboration.

    Data availability

    The data that support the results in this paper and other findings of this study are available from the corresponding author upon a reasonable request.

    Competing interests

    The authors declare no conflict of interest.

    Supplementary informationThe online version contains supplementary material available at https://doi.org/10.1038/s41377-021-00703-6.

    黑人巨大精品欧美一区二区蜜桃 | 国产成人精品福利久久| 两个人免费观看高清视频 | 久久久久久伊人网av| 国产极品天堂在线| 国产永久视频网站| 男的添女的下面高潮视频| av网站免费在线观看视频| 寂寞人妻少妇视频99o| 草草在线视频免费看| 美女xxoo啪啪120秒动态图| 亚洲四区av| 内地一区二区视频在线| 黄色一级大片看看| 国产极品粉嫩免费观看在线 | 91久久精品国产一区二区成人| 久久久久国产精品人妻一区二区| 在线观看免费视频网站a站| 亚洲成色77777| 免费观看a级毛片全部| 国产亚洲最大av| 色婷婷久久久亚洲欧美| 亚洲三级黄色毛片| 日本-黄色视频高清免费观看| 又大又黄又爽视频免费| 王馨瑶露胸无遮挡在线观看| 亚洲精品456在线播放app| 亚洲欧美成人精品一区二区| 亚洲精品亚洲一区二区| 成人特级av手机在线观看| 女性被躁到高潮视频| 伊人久久国产一区二区| 国产av精品麻豆| 最黄视频免费看| 国产精品欧美亚洲77777| 久久午夜福利片| 99久久综合免费| 国产伦理片在线播放av一区| 内地一区二区视频在线| 纯流量卡能插随身wifi吗| 新久久久久国产一级毛片| 欧美一级a爱片免费观看看| 精品一区二区三卡| 精品亚洲成国产av| 国产在线男女| 免费大片黄手机在线观看| 日韩精品免费视频一区二区三区 | 男人添女人高潮全过程视频| 国产伦理片在线播放av一区| 夫妻性生交免费视频一级片| 2022亚洲国产成人精品| 国产欧美日韩综合在线一区二区 | 国产亚洲精品久久久com| 亚洲精品国产av蜜桃| 久久精品国产亚洲av涩爱| 久久久午夜欧美精品| 免费av中文字幕在线| 精品一区二区三区视频在线| 亚洲国产成人一精品久久久| 久久国内精品自在自线图片| 欧美老熟妇乱子伦牲交| 一级av片app| 另类精品久久| 内地一区二区视频在线| 美女中出高潮动态图| 久久精品国产亚洲av涩爱| 一级毛片黄色毛片免费观看视频| a级毛片在线看网站| 69精品国产乱码久久久| 国产成人91sexporn| 国国产精品蜜臀av免费| 最近中文字幕高清免费大全6| 插阴视频在线观看视频| 日韩一区二区三区影片| 91成人精品电影| 人妻一区二区av| av.在线天堂| 亚洲中文av在线| 在线观看av片永久免费下载| 18禁在线播放成人免费| 一级av片app| 国产成人freesex在线| 亚洲美女黄色视频免费看| 18禁在线播放成人免费| 内地一区二区视频在线| 亚洲欧美一区二区三区黑人 | 亚洲va在线va天堂va国产| 亚洲欧美日韩另类电影网站| 亚洲国产精品国产精品| 久久久久久久久久成人| 国产有黄有色有爽视频| 亚洲一区二区三区欧美精品| 搡女人真爽免费视频火全软件| 国国产精品蜜臀av免费| 99九九在线精品视频 | 亚洲精品日韩在线中文字幕| 九九久久精品国产亚洲av麻豆| 国产精品三级大全| 99九九在线精品视频 | 欧美日韩在线观看h| 免费不卡的大黄色大毛片视频在线观看| 国产精品女同一区二区软件| 中文字幕精品免费在线观看视频 | 亚洲欧洲精品一区二区精品久久久 | 免费观看性生交大片5| 亚洲av免费高清在线观看| 高清毛片免费看| 国产成人freesex在线| 精品人妻偷拍中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 色视频在线一区二区三区| 黄色视频在线播放观看不卡| 少妇熟女欧美另类| 三上悠亚av全集在线观看 | 蜜臀久久99精品久久宅男| 看非洲黑人一级黄片| 精品久久久久久久久av| 嫩草影院入口| 啦啦啦啦在线视频资源| 纯流量卡能插随身wifi吗| av国产精品久久久久影院| 国产精品一区www在线观看| 国产乱人偷精品视频| 在线播放无遮挡| 人体艺术视频欧美日本| 亚洲av中文av极速乱| 国产亚洲一区二区精品| 久久韩国三级中文字幕| videossex国产| 超碰97精品在线观看| 赤兔流量卡办理| 欧美xxⅹ黑人| 最新中文字幕久久久久| 又粗又硬又长又爽又黄的视频| 欧美性感艳星| 亚洲欧美日韩卡通动漫| 亚洲成人av在线免费| 十八禁网站网址无遮挡 | 欧美激情国产日韩精品一区| 免费高清在线观看视频在线观看| 久久99蜜桃精品久久| 欧美性感艳星| 久久久久精品性色| av专区在线播放| 国产高清不卡午夜福利| 最新中文字幕久久久久| 国产综合精华液| 久久99热6这里只有精品| av在线老鸭窝| 制服丝袜香蕉在线| 国产伦精品一区二区三区视频9| 99九九线精品视频在线观看视频| 日韩中字成人| 中文欧美无线码| 菩萨蛮人人尽说江南好唐韦庄| 国产精品免费大片| 婷婷色综合www| 久久精品夜色国产| 久久久久久人妻| 乱人伦中国视频| 纯流量卡能插随身wifi吗| 午夜免费观看性视频| 夜夜看夜夜爽夜夜摸| 日韩大片免费观看网站| 久久午夜福利片| 亚洲精品乱码久久久v下载方式| 在线 av 中文字幕| 国产精品一二三区在线看| h视频一区二区三区| 国产高清不卡午夜福利| 日韩制服骚丝袜av| 久久精品久久精品一区二区三区| 一二三四中文在线观看免费高清| 国产精品伦人一区二区| 日韩精品免费视频一区二区三区 | 欧美激情国产日韩精品一区| 精品人妻熟女毛片av久久网站| 人人妻人人澡人人爽人人夜夜| 国产在线免费精品| 人妻一区二区av| 日本-黄色视频高清免费观看| 国内少妇人妻偷人精品xxx网站| 七月丁香在线播放| av国产久精品久网站免费入址| 欧美日韩亚洲高清精品| 欧美+日韩+精品| 国产欧美日韩精品一区二区| 色婷婷久久久亚洲欧美| 亚洲国产成人一精品久久久| 精品亚洲乱码少妇综合久久| 少妇的逼水好多| 亚洲电影在线观看av| 久久久久精品久久久久真实原创| 男女边吃奶边做爰视频| 午夜视频国产福利| 麻豆成人av视频| 日本免费在线观看一区| 亚洲自偷自拍三级| 国产精品熟女久久久久浪| 伊人久久国产一区二区| 不卡视频在线观看欧美| 国产精品国产三级专区第一集| 大码成人一级视频| 日韩欧美精品免费久久| 成年av动漫网址| 精品人妻熟女毛片av久久网站| 国产亚洲5aaaaa淫片| 色哟哟·www| 麻豆成人午夜福利视频| 亚洲av成人精品一区久久| 国产av精品麻豆| 七月丁香在线播放| 亚洲国产最新在线播放| 熟妇人妻不卡中文字幕| 人妻制服诱惑在线中文字幕| 黄片无遮挡物在线观看| 伦理电影免费视频| 国产美女午夜福利| 最近中文字幕2019免费版| 国产高清国产精品国产三级| 亚洲第一av免费看| 亚洲天堂av无毛| 国产精品熟女久久久久浪| 久久久久久久久久成人| 97在线视频观看| 亚洲人成网站在线播| 中文天堂在线官网| 亚洲国产最新在线播放| 黄色配什么色好看| 一级二级三级毛片免费看| 高清毛片免费看| 午夜日本视频在线| 国产精品熟女久久久久浪| 亚洲欧美日韩东京热| 2018国产大陆天天弄谢| 久久久久久久大尺度免费视频| av福利片在线观看| 欧美日韩视频高清一区二区三区二| 一级爰片在线观看| 精品一区二区免费观看| 国产无遮挡羞羞视频在线观看| 精品熟女少妇av免费看| 国精品久久久久久国模美| 男女无遮挡免费网站观看| 婷婷色综合大香蕉| 少妇熟女欧美另类| 91久久精品电影网| 国产亚洲精品久久久com| 亚洲欧美日韩东京热| 国产精品免费大片| 精品亚洲成a人片在线观看| 只有这里有精品99| 国产精品久久久久成人av| 国产毛片在线视频| 男女国产视频网站| 韩国高清视频一区二区三区| h视频一区二区三区| 亚洲精华国产精华液的使用体验| 极品教师在线视频| 国产在线免费精品| 3wmmmm亚洲av在线观看| 啦啦啦中文免费视频观看日本| 爱豆传媒免费全集在线观看| 各种免费的搞黄视频| 久久99热这里只频精品6学生| 狂野欧美激情性bbbbbb| 精品一品国产午夜福利视频| 欧美精品一区二区大全| 九色成人免费人妻av| 国产精品99久久99久久久不卡 | 国产精品久久久久久久电影| 午夜av观看不卡| 午夜91福利影院| 精品视频人人做人人爽| 一区在线观看完整版| 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 日本黄大片高清| 国产精品一区二区在线不卡| 亚洲美女搞黄在线观看| 我的老师免费观看完整版| 国产亚洲91精品色在线| 丝袜脚勾引网站| 国产精品国产av在线观看| av播播在线观看一区| 久久韩国三级中文字幕| 丁香六月天网| 天天躁夜夜躁狠狠久久av| 久久午夜综合久久蜜桃| 中文字幕亚洲精品专区| 插逼视频在线观看| 一级片'在线观看视频| 午夜福利网站1000一区二区三区| 黑人猛操日本美女一级片| 韩国高清视频一区二区三区| 免费大片黄手机在线观看| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区| 在线 av 中文字幕| 一区二区三区乱码不卡18| av在线app专区| 99精国产麻豆久久婷婷| 国产伦理片在线播放av一区| 91精品一卡2卡3卡4卡| 亚洲欧美精品专区久久| av在线播放精品| 18禁动态无遮挡网站| 亚洲成色77777| 在线观看国产h片| 三上悠亚av全集在线观看 | 国产精品欧美亚洲77777| 中文乱码字字幕精品一区二区三区| kizo精华| 丰满人妻一区二区三区视频av| 秋霞在线观看毛片| h日本视频在线播放| 国产一区亚洲一区在线观看| 欧美日韩av久久| 国产在线视频一区二区| 亚洲怡红院男人天堂| 熟女人妻精品中文字幕| 国产一区二区三区综合在线观看 | 18禁动态无遮挡网站| 波野结衣二区三区在线| 少妇的逼好多水| 国产精品蜜桃在线观看| 国产免费一级a男人的天堂| av在线播放精品| 夜夜爽夜夜爽视频| 亚洲美女黄色视频免费看| av国产精品久久久久影院| 老熟女久久久| 纵有疾风起免费观看全集完整版| 校园人妻丝袜中文字幕| 精品亚洲成国产av| 伊人亚洲综合成人网| 99九九线精品视频在线观看视频| av在线app专区| 一级,二级,三级黄色视频| 日本av手机在线免费观看| av在线老鸭窝| 国产精品成人在线| 午夜福利影视在线免费观看| 18+在线观看网站| 久久99精品国语久久久| 免费观看的影片在线观看| 久久久欧美国产精品| 69精品国产乱码久久久| 婷婷色综合www| 蜜桃在线观看..| a级毛片在线看网站| 亚洲av福利一区| 男人和女人高潮做爰伦理| 国产乱人偷精品视频| 久久热精品热| 18禁在线播放成人免费| 久久久久久久久久久丰满| 91午夜精品亚洲一区二区三区| 免费观看无遮挡的男女| 久久国产精品男人的天堂亚洲 | 国产亚洲一区二区精品| 插逼视频在线观看| 99九九线精品视频在线观看视频| 极品少妇高潮喷水抽搐| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 国产精品不卡视频一区二区| 久久久久久伊人网av| 久久99热6这里只有精品| 熟女电影av网| 9色porny在线观看| 涩涩av久久男人的天堂| 国产精品偷伦视频观看了| 婷婷色综合www| 国产一区二区三区综合在线观看 | 各种免费的搞黄视频| 成人二区视频| 国产在线免费精品| 大陆偷拍与自拍| 最新中文字幕久久久久| 一级黄片播放器| 精品国产露脸久久av麻豆| 日韩强制内射视频| 麻豆精品久久久久久蜜桃| 热99国产精品久久久久久7| 国产69精品久久久久777片| .国产精品久久| 免费观看性生交大片5| 久热这里只有精品99| 国产一区二区在线观看av| 婷婷色av中文字幕| 高清黄色对白视频在线免费看 | 狂野欧美激情性bbbbbb| 精品少妇黑人巨大在线播放| 国产精品欧美亚洲77777| av黄色大香蕉| 少妇人妻 视频| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜制服| 国产黄片视频在线免费观看| 国产午夜精品一二区理论片| 亚洲av.av天堂| 国产淫片久久久久久久久| 少妇高潮的动态图| 一级毛片电影观看| 亚洲国产精品国产精品| 纵有疾风起免费观看全集完整版| 最近中文字幕高清免费大全6| 亚洲人与动物交配视频| av女优亚洲男人天堂| 久久久欧美国产精品| 成人国产麻豆网| 下体分泌物呈黄色| 国产在线免费精品| 2018国产大陆天天弄谢| 亚洲一级一片aⅴ在线观看| 久久国内精品自在自线图片| 亚洲无线观看免费| 女人精品久久久久毛片| 国产精品99久久99久久久不卡 | 在线精品无人区一区二区三| 26uuu在线亚洲综合色| 久久6这里有精品| 免费看不卡的av| 乱人伦中国视频| 国产一区二区三区av在线| a 毛片基地| 熟女电影av网| 日本免费在线观看一区| 久久久国产一区二区| av在线app专区| 亚洲精品日韩av片在线观看| 最近手机中文字幕大全| 婷婷色综合www| 亚洲在久久综合| 国产精品免费大片| 丝袜在线中文字幕| 国内少妇人妻偷人精品xxx网站| 国产黄色免费在线视频| 国产黄频视频在线观看| 男人爽女人下面视频在线观看| 人妻制服诱惑在线中文字幕| 亚洲精品国产色婷婷电影| av黄色大香蕉| a级毛片在线看网站| 国产成人精品久久久久久| 男的添女的下面高潮视频| 2021少妇久久久久久久久久久| 波野结衣二区三区在线| 成人特级av手机在线观看| 亚洲不卡免费看| 一区二区三区精品91| 97精品久久久久久久久久精品| 夫妻性生交免费视频一级片| 亚洲成人手机| 精品久久久久久久久av| 欧美老熟妇乱子伦牲交| 亚洲国产av新网站| 欧美丝袜亚洲另类| 久久热精品热| av不卡在线播放| 欧美日韩亚洲高清精品| 久久久精品免费免费高清| 99热网站在线观看| 丝袜在线中文字幕| 最新的欧美精品一区二区| 免费av中文字幕在线| 爱豆传媒免费全集在线观看| 精品99又大又爽又粗少妇毛片| 街头女战士在线观看网站| 王馨瑶露胸无遮挡在线观看| 免费看光身美女| 青青草视频在线视频观看| 麻豆成人av视频| 亚洲情色 制服丝袜| 在线播放无遮挡| 国产精品一区二区在线不卡| 午夜视频国产福利| 精品久久国产蜜桃| 在线观看人妻少妇| 日韩中字成人| 欧美丝袜亚洲另类| 大陆偷拍与自拍| 日韩成人av中文字幕在线观看| 十八禁高潮呻吟视频 | 卡戴珊不雅视频在线播放| 国产一区二区三区综合在线观看 | 日本免费在线观看一区| 99九九在线精品视频 | 天堂8中文在线网| 国产av一区二区精品久久| 国产成人精品婷婷| tube8黄色片| 69精品国产乱码久久久| 国产成人aa在线观看| 国产伦理片在线播放av一区| 黑丝袜美女国产一区| 综合色丁香网| 久久久欧美国产精品| 国产亚洲欧美精品永久| videos熟女内射| 我的女老师完整版在线观看| 久久久久久人妻| 亚洲欧洲国产日韩| 中文字幕制服av| 久久99蜜桃精品久久| 日本黄色片子视频| 国产精品成人在线| 欧美丝袜亚洲另类| 国产精品成人在线| 日本黄色片子视频| 少妇人妻一区二区三区视频| 韩国高清视频一区二区三区| 国产精品福利在线免费观看| 天堂俺去俺来也www色官网| 国产亚洲5aaaaa淫片| 人妻少妇偷人精品九色| 久久热精品热| 97精品久久久久久久久久精品| 男人舔奶头视频| 最近的中文字幕免费完整| 婷婷色麻豆天堂久久| 久久人人爽人人片av| 国产亚洲5aaaaa淫片| 国产无遮挡羞羞视频在线观看| 国产精品成人在线| 人妻系列 视频| 日产精品乱码卡一卡2卡三| 2022亚洲国产成人精品| 亚洲欧美清纯卡通| 免费av不卡在线播放| 亚洲,一卡二卡三卡| 午夜福利,免费看| 久久午夜综合久久蜜桃| 在线精品无人区一区二区三| 亚洲色图综合在线观看| h日本视频在线播放| 精品国产国语对白av| 国产av码专区亚洲av| 美女cb高潮喷水在线观看| 久久99一区二区三区| 下体分泌物呈黄色| 日韩人妻高清精品专区| 99热6这里只有精品| 黄片无遮挡物在线观看| 美女视频免费永久观看网站| 看非洲黑人一级黄片| 丝瓜视频免费看黄片| 一级毛片 在线播放| 国产一区有黄有色的免费视频| 各种免费的搞黄视频| 男女无遮挡免费网站观看| 欧美少妇被猛烈插入视频| 亚洲精品一二三| 国产黄片视频在线免费观看| 三级经典国产精品| 一区二区三区乱码不卡18| 亚洲av日韩在线播放| 午夜激情福利司机影院| 欧美日韩在线观看h| 在线亚洲精品国产二区图片欧美 | 久久99精品国语久久久| 男人狂女人下面高潮的视频| 五月伊人婷婷丁香| 日韩强制内射视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产av码专区亚洲av| 一本一本综合久久| 男女免费视频国产| 18禁在线无遮挡免费观看视频| 麻豆成人午夜福利视频| 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 岛国毛片在线播放| 91精品国产国语对白视频| 久久久久久久久久人人人人人人| 免费观看的影片在线观看| 边亲边吃奶的免费视频| 国产日韩欧美视频二区| av免费观看日本| 亚洲国产精品国产精品| 国产精品免费大片| 十八禁网站网址无遮挡 | 欧美变态另类bdsm刘玥| 狂野欧美激情性xxxx在线观看| 99re6热这里在线精品视频| 中文字幕av电影在线播放| 国产精品.久久久| 国产精品熟女久久久久浪| 汤姆久久久久久久影院中文字幕| 纯流量卡能插随身wifi吗| 成人国产av品久久久| 久久精品久久久久久噜噜老黄| 蜜桃久久精品国产亚洲av| 亚洲精品乱久久久久久| 亚洲精品一二三| 国产成人a∨麻豆精品| 久久久久久久久久久免费av| 五月天丁香电影| 黄片无遮挡物在线观看| 久久久久久久久久久丰满| 国产一区亚洲一区在线观看| 两个人免费观看高清视频 | 欧美+日韩+精品| 亚洲高清免费不卡视频| 两个人的视频大全免费| 丝袜在线中文字幕| 最近最新中文字幕免费大全7| 日韩欧美一区视频在线观看 | 色94色欧美一区二区| 亚洲欧美清纯卡通| 精品人妻偷拍中文字幕| 久热久热在线精品观看| 亚洲怡红院男人天堂| 日韩视频在线欧美|