1. Experimental Study of Long-Wavelength Mid-Infrared Laser Sources Based on DFG
Early DFG schemes are mainly based on bulky Ti: Sapphire femtosecond pulsed laser systems to generate dual-wavelength synchronous pulses, restricting the average power of the MIR below the milliwatt level due to the thermal load
[95,96,97,98,99,100,101,102][1][2][3][4][5][6][7][8]. Modern DFG schemes that can achieve relatively high powers are mainly based on dual-wavelength OPO, OPA, and ultrafast fiber laser systems
[103,104,105,106,107,108,109,110,111,112,113,114,115,116,117][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23].
A MIR femtosecond pulsed laser with a wide tuning range (5–17 μm) and an average power of 69 mW at the wavelength of 6 μm was achieved by Beutler et al. based on the DFG in AgGaSe
2 crystals utilizing a synchronously pumped OPO ultrafast laser source in 2015
[106][12]. A DFG-based MIR laser with an average power of 4.3 mW and a 10.5–16.5 μm wavelength-tunable range (as shown in
Figure 1) was built by employing a two-color femtosecond OPO system in 2012 by Harald Giessen’s research group of the Fourth Institute of Physics at the University of Stuttgart, Germany
[107][13]. Later in 2016, a highly stable MIR 350-fs pulsed laser system with nJ-level pulse energy and an ultra-wide tunability ranging from 1.33–2.0 μm and 2.13–20 μm has been achieved by Giessen’s group by combining OPO, OPA, and DFG technologies
[108][14]. More recently, Wilson et al. proposed a DFG scheme for MIR generation tunable from 3.5 μm to 9 μm in 2019. They utilized a few-cycle Ti: Sapphire laser to pump an OPA to generate a total power of up to 6 mJ of signal (1300–1450 nm) and idler (1730–2000 nm) lights in the NIR. Then, the two synchronous pulses were used for the DFG process, and a maximum of 120 μJ at 5.3 μm and up to 80 μJ at 8.9 μm center wavelengths with 70–90 fs of pulse width were obtained through their setup (as shown in
Figure 2)
[109][15].
Figure 1. A 10.5–16.5 μm tunable femtosecond MIR pulsed laser system. Here, a 1 mm long MgO-doped periodically poled lithium niobate (PPLN) crystal is used for the dual-signal-wavelength OPO working at room temperature, which is synchronously pumped by a 1040 nm, 7.4 W, and 530 fs Yb: KGW oscillator at 42 MHz repetition rate. The dual-wavelength output beams of the OPO are focused on the GaSe crystal for DFG-based MIR generation.
Figure 2. A DFG scheme for MIR generation tunable from 3.5 μm to 9 μm
[109][15]. An intense, few-cycle source in the long-wave infrared., Wilson, D.J.; Summers, A.M.; Zigo, S.; Davis, B.; Robatjazi, S.J.; Powell, J.A.; Rolles, D.; Rudenko, A.; Trallero-Herrero, C.A., COPYRIGHT 2019. (
a) a 20 mJ Ti: Sapphire 800 nm laser is used to pump an optical parametric amplifier (OPA) with 18 mJ, 26 fs pulses at 1 kHz. A total of 6 mJ signal and idler pulses are generated from the OPA, and up to 4.2 mJ are used to generate LWIR (MIR) pulses via DFG in a AgGaS
2 crystal (AGS) through type-I phase-matching. The residual pump and signal pulses are filtered by anti-reflective coated Zinc Selenide (ZS) and Germanium (Ge) windows. The Fresnel reflection of the Ge window is filtered with an additional Ge window so residual MIR power can be recorded on a pyroelectric detector (PyD) for power monitoring. A periscope sends the remaining pulse energy (>95%) into an experimental chamber with a ZS window, and then the MIR pulse is back-focused in Xenon (Xe) gas by a −25 mm spherical concave mirror (F). A pair of electrostatic lenses guide the Xe ions onto a micro-channel plate (MCP) detector for time-of-flight measurements. A small portion of the remaining 2 mJ OPA beam (pink) is used as a gating field in the XFROG. To perform electric field characterization, Mirror M is first removed from the periscope, and then the MIR pulse is coupled into the XFROG. (
b) A typical time-of-flight spectrum for Xe
+ when ionized with 8.9 μm pulse. (
c) A measured XFROG spectrogram when the MIR laser pulse is working at 8.9 μm.
Ultrafast fiber lasers have been adopted as the pump sources in some DFG experiments due to their high surface-area-to-volume ratios, strong heat dissipation performances, high repetition rates and average powers, as well as all-fiber compactness. The output pulses from a fiber laser are usually divided into two parts for MIR generation. One part is used as pump pulses and the other produces a new spectral component on the longer wavelength side through the nonlinear fiber effect as wavelength-tunable signal pulses
[110,111,112,113,114,115,116,117][16][17][18][19][20][21][22][23].
In 2013, Gambetta et al. reported a MIR pulsed laser with a tunability of 8–14 μm and a maximum average power of 4 mW in GaSe crystals (as shown in
Figure 3)
[111][17]. In the experiment, the Er: fiber laser oscillator (1.55 μm, 50 fs, and 550 mW) was used as the pump pulse, while the self-frequency shifted (SFS) soliton pulse (1.76–1.93 μm) generated from the nonlinear Raman fiber was serving as the signal pulse. The synchronism of the two pulses was adjusted through the optical delay line, and the tunable output of the MIR laser pulses was generated in a 1 mm-thick GaSe crystal. Zhou et al. developed a Yb: fiber laser system with an average power of 14.5 W, a repetition frequency of 30 MHz, and a pulse width of 165 fs. The output pulses of the system are divided into two parts in the optical path through the half-wave plate (HWP) and the polarization beam splitter (PBS). A portion of the pulses acted directly as pump pulses of the DFG. The other part of the pulses passed through a large-mode field fiber (4 cm), then the 125 fs signal pulses with the wavelength-tunable range of 1030–1215 nm were generated through the self-phase modulation effect. In the GaSe crystal, MIR pulses with a highest average power of 5.4 mW and a tunability between 7.4 and 18 μm (as shown in
Figure 4) were finally obtained
[114][20].
Figure 3. The 8–14 μm tunable MIR laser pulse generation. The experimental setup used for the mid-IR comb generation is based on a commercial dual-branch Er: fiber laser oscillator. The main output provides 1.55 μm, 550 mW, and 50 fs pulses, which are used as pump pulses. The second output produces tunable SFS solitons (as short as 84 fs) in the range 1.76–1.93 μm with an optical power varying from 100 to 250 mW by coupling pulse into a nonlinear Raman fiber, which are used as signal pulses.
Figure 4. Schematic diagram of the ultrafast MIR sources with a tunable range of 7–18 μm. A home-built, 14.5 W, 165 fs, and 30 MHz Yb: fiber system is used as the front-end experimental setup for the MIR generation. The laser output is split into two paths by a half-wave plate (HWP) together with a polarization beam splitter (PBS). One path of light is coupled into a piece of large-mode-area (LMA) fiber for SPM-enabled spectra broadening, and then a bandpass filter is used for targeted wavelength selection to generate tunable femtosecond signal pulses. The other optical path serves as the pump for subsequent DFG in a GaSe crystal.
In 2019, Ma et al. proposed an all-polarization-maintaining fiber amplification method (except the optical delay line) for MIR generation
[117][23]. The dual-wavelength pulses are split by an optical fiber coupler and enter different nonlinear amplification systems. Finally, the peak wavelengths of the two pulses are separated at least by 50 nm. The MIR generation from the system can be tuned from 7 to 10.5 µm. In 2018, Sotor et al. reported a MIR source ranging from 6 to 9 µm based on an all-fiber system
[115][21]. The pulses from the oscillator were split and travelled through different fiber amplification models. However, a fiber-coupled delay line was used to control the synchronism of the pulses, making the front-end system truly all-fiberized (as shown in
Figure 5). The maximum output power of the DFG was 7.4 mW with a repetition rate of 48 MHz at 7.5 µm, and the power was kept at 1 mW over the entire tuning range. More recently, Sotor et al. updated their system by adding a repetition rate stabilization model and a pulse-synchronism controlling model
[116][22], which realized the stabilization of the repetition rate and the temporal pulse overlap in the MIR generation. An average power of 5 mW was obtained in the tunable range of 6.5–9 µm at 125 MHz (as shown in
Figure 6).
Figure 5. MIR source ranging from 6 to 9 µm based on an all-fiber system. The pulses from the Er-doped fiber oscillator (OSC) are equally split into two branches. One part is amplified in an Er: fiber amplifier (EDFA-1) to approximately 120 mW and enters a polarization-maintaining (PM) highly nonlinear fiber (HNLF). The spectrum of the input pulse is then shifted towards longer wavelengths via Raman-induced soliton self-frequency shift (SSFS) in the HNLF. Finally, the output pulses from the HNLF are amplified to 150–235 mW in a Tm: fiber amplifier (TDFA), serving as the signal pulses for the DFG process. The second part is amplified to 175 mW in EDFA-2. Both arms are combined together into a common single-mode fiber with a 1550/2000 nm wavelength division multiplexer (WDM). The temporal overlap between the pump and signal pulses is adjusted by a fiber-coupled variable delay line (VDL). The pump and signal beams are focused on the nonlinear crystal together for MIR generation. The generated MIR (idler) beam is collimated with a 75 mm barium fluoride (BaF2) lens. The residual pump and signal lights are blocked by a long-pass filter (F) with a 4.5 µm cut-on wavelength, meaning that only MIR pulses can travel through it.
Figure 6. All-fiber MIR source with a tunable range of 6.5–9 µm. The system is composed of four main parts: (1) a mode-locked (ML) Er: fiber seed laser, (2) a 1.56 µm chirped pulse amplification (CPA) stage, (3) a stage for achieving a soliton shift for switching to 2 µm region and subsequent amplification of the ~2 µm red-shifted part of the spectrum through a CPA setup, and (4) difference frequency generation.
In recent years, the intrapulse difference frequency generation method (IDFG, also known as optical rectification) has become a new method to generate MIR laser source owning to its simplicity and phase stability
[118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41]. The coherent femtosecond-level pulses with a wavelength range covering the whole MIR region are successively realized through IDFG. Only one beam of femtosecond pulses is needed to travel through the nonlinear crystal to obtain MIR pulses, making it structurally simpler than general DFG. However, IDFG involves a variety of nonlinear effects, and the mechanism is relatively unclear. The demand for relevant pump pulse width (usually few-cycle pulses) is also challenging to meet. Furthermore, tunability seems hard to realize since a very broad spectrum is generated via IDFG. In addition, for output wavelength longer than 10 μm, especially over 15 μm, the optical-optical conversion efficiency (MIR power/pump power) drops by tens of times compared to the short wavelength (see
Table 1) regardless of the traditional DFG technology or the emerging IDFG technology.
Table 1. Studies on the generation of ultrafast MIR sources above 8 μm based on DFG and IDFG
[118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41].
Donna Strickland’s ultrafast laser group has been working on long-wavelength MIR pulse generation for tens of years, aiming on output wavelengths of over 16 μm based on the dual-wavelength ultrashort pulsed laser system and the DFG method
[136,137,138,139,140][42][43][44][45][46]. In 2012, an 18 μm pulsed laser with an average power of 30 μW under subpicosecond level was achieved with the Yb: fiber amplification system, and the wavelength is tunable between 16 and 20 μm
[138][44]. However, the experimental system is composed of bulk-optic elements, which is large, inefficient, and unstable. In 2018, they switched the whole bulk-optic system into a partially fiberized structure by replacing the traditional spatial grating with the chirped fiber Bragg grating as the wavelength selector, and two target wavelengths (1025 nm and 1085 nm) are selected. The dual-wavelength seed source was then amplified in the subsequent two-stage amplification system. As a result, the compact two-color femtosecond chirped pulse amplification (CPA) system with high average power (as shown in
Figure 7) was achieved
[139][45]. Furthermore, based on this system, a 17.4 μm MIR pulsed laser output with an average power of 2.5 mW was achieved by DFG under the condition of tight-focusing
[140][46]. Compared to the previous work in 2012, not only was the MIR power increased by about 80 times, but the stability was also improved
[138][44].
Figure 7. The dual-wavelength CPA system built by Donna Strickland’s group. A femtosecond laser system delivers 200 fs, 400 mW average power, and 65 MHz repetition rate laser pulses. The laser is then coupled into a photonic crystal fiber to achieve supercontinuum generation ranging from 900 to 1100 nm. A chirped fiber Bragg grating (CFBG) blocking the transmission from 1025 to 1085 nm is used to select the dual-wavelength seed source. Finally, the dual-wavelength seed source enters the CPA system to achieve final amplification with a total average power of 2.3 W (1.9 W at 1025 nm and 0.4 W at 1085 nm). The full width at half-maximum (FWHM) durations of the final compressed output pulses of each color are 900 fs and 600 fs, respectively.
In conclusion, the efficiency and stability of the current DFG laser systems can be vigorously improved if the following drawbacks can be addressed:
First of all, the components of the most current systems are partially fiberized, and structures of bulk-optic components are still adopted in the amplification and pulse compression modules, resulting in low conversion efficiency, significant loss, and strong instability.
Secondly, the increase of the total power of the amplification system is limited by the drawbacks in the process of seed spectrum selection and amplification. The effective spectral component that can be amplified in the seed spectrum selection process only accounts for less than one-tenth of the total power, and the other part of the energy in the supercontinuum is sacrificed. In addition, there are serious gain competitions in the amplification process when the dual-wavelength seeds are amplified simultaneously.
Thirdly, the impact of the accumulation of nonlinear phase shift on the final output pulse width is not considered in the whole CPA system. Compared with the initial seed pulse, the compressed pulse is wider and owns a prominent pedestal. The pulse quality still has a large room for improvement.
Last but not least, the generated MIR wavelength range can only cover 16–20 μm, which is a limited wavelength-tunable range. In addition, the power is also limited to several milliwatts. These two factors severely restrict the application field of the MIR laser source.
2. Theoretical Study of the DFG for Mid-Infrared Generation
In terms of theoretical research, there are few relevant theoretical calculations on transient DFG processes, and the transient model is usually complex. Therefore, the average power of the idler frequency light is always estimated with the steady-state DFG formula. However, since the formula is established under steady-state conditions, it is defective in the states of ultrashort pulses. To this end, the steady-state DFG formula should be modified. At present, the representative DFG formula of correction calculation of MIR power is as follows
[141][47]:
In Equation (1),
Tp,
Ts, and
Ti are Fresnel transmission coefficients of pump, signal, and idler pulse in nonlinear crystals, respectively.
deff is the effective second-order nonlinear coefficient.
Pp and
Ps are the average powers of the pump and signal light.
L is the interaction length of the pump, signal, and idler light in the nonlinear optical crystal.
ε0 is the dielectric constant in the vacuum.
c is the speed of light in the vacuum.
ni,
np, and
ns are the refractive index of nonlinear optical crystals towards the idler, pump, and signal light at the phase-matching angle, respectively.
λi is the wavelength of the idler pulse in a vacuum,
f is the pulse repetition rate.
t is the average value of the pump, signal, and idler pulse when the pulse widths of the three are assumed to be equal.
g is the correction coefficient after considering the walk-off effect. The formula is modified compared to the DFG calculation formula in the continuous-wave state by taking the influence of peak powers into consideration. However, it is still a rough model because the value of the correction coefficient
g (0 <
g < 1) is not well defined. Moreover, the formula does not consider the thermal effect of the crystal, and there is a significant deviation between the calculation results and the actual experiments
[138,141][44][47]. Therefore, the computational modified model of the average power of DFG should be reconstructed from the transient equation of the DFG process combined with an amount of experimental data.
Moreover, it is mentioned above that the physical process of OPA and DFG is basically equivalent, hence the generation of idler frequency light also happens in the realization of OPA. This raises new questions: can a critical condition be found during the DFG process by changing the laser power injected into nonlinear crystals to shift the DFG into the OPA process? If so, will the power of the idler frequency behave like an exponential nonlinear growth trend? Efforts have been made to explain the relevant problems through theoretical calculation. For the repetition rate of 30 MHz and the pulse width of 200 fs, the evolution curve of the idler pulse energy was obtained by varying with the pump and signal pulse energy, predicting that there are linear growth, nonlinear amplification, and saturation zones in the DFG process (as shown in
Figure 8)
[142][48]. It can be seen that when the signal light power is kept to a certain value, the pump power plays a decisive role in the improvement of the MIR (idler light) power, and the increase of the signal power can easily make the MIR power reach saturation. Therefore, when certain critical conditions of nonlinear amplification are met and the DFG process is transformed into an OPA process, the MIR power will be nonlinearly increased, and the optical-optical conversion efficiency will be greatly improved (>1%).
Figure 8. Energy evolution curve of idler pulse energy obtained through theoretical calculations. (a) Evolution of the idler energy versus the input pump energy when the input signal energy is fixed at 100 pJ, 1 nJ, and 10 nJ. (b) Evolution of the idler energy versus the input pump energy when the input signal energy is fixed at 1 nJ, 10 nJ, and 100 nJ. The thickness of GaSe crystal for MIR generation is 2 mm.
3. Feasibility of Mid-Infrared Power Scaling by Nonlinear-Amplification DFG
HRere we presentgarding the current problems for MIR ultrashort pulse generation above 8 μm. They are as follows:
(1) For current MIR lasers based on DFG or IDFG schemes, when the output wavelength is longer than 10 μm or even longer than 15 μm, the power will drop rapidly to the order of several milliwatts, and the optical-optical conversion efficiency is relatively low.
(2) Most of the research groups separate the relationship between DFG and OPA, and do not consider the nonlinear amplification and saturation states that appear in the DFG process in theoretical and experimental research.
From the author’s point of view, there is no doubt that IDFG could generate broadband MIR sources with relatively higher powers. However, the problem is that the output MIR spectrum is not easily controllable and tunable. Except for picking up crystals with high nonlinearity and transmission, there is much work to be done in the front-end pumping system due to the complexity of overcoming dispersion and nonlinearity problems for obtaining the few-cycle pulses. Compared to IDFG, DFG is still a cost-effective scheme since fiber lasers under hundreds of femtosecond level as pump sources are commercially available and relatively economic. To solve the problem of low power of MIR for the DFG scheme, the following solutions with nonlinear-amplification DFG are proposed:
(1) Make reasonable adjustments to the focused beam size on the crystals. The focused spot size greatly affects the phase matching of the DFG process. In order to achieve high-power MIR pulsed laser under DFG scheme, an appropriate laser focusing lens should be selected to optimize the spot size focused on the nonlinear crystal. A small focused spot size will lead to a higher MIR power, but the divergence angle will be too large for the MIR beam to be collected and serious thermal effects and damaging problems will occur. On the other hand, the divergence angle will become small, and the MIR beam will be easy to collect when the focused spot size is large, but the intensity threshold of DFG in the nonlinear amplification region will not be reached. Therefore, the focused spot size should be optimized to make the focused light intensity higher than the threshold for nonlinear amplification and less than the crystal damage threshold to achieve high power MIR pulses.
(2) Adjust the power/pulse energy ratio entering the crystal. As predicted in
[142][48], there are linear growth, nonlinear amplification and saturation regions in the DFG process. The MIR power is proportional to the product of the pump and signal light powers within a certain range. However, a large number of experiments show that it is not feasible to increase the signal power to the same level as the pump power since the wavelengths of the signal are usually not located in the gain curve of the materials.
We propose thaIt i
t is better to deliberately create a power gap between the pump and signal lights, which leads to the transformation of the linear growth state to the nonlinear amplification state and results in a conversion efficiency of >1%. Assuming a pump laser source with 10 W average power and a signal light with 100 mW average power, if the nonlinear amplification state is reached and kept, an average power of no less than 100 mW can be maintained throughout the whole MIR range.
4. Methods for Measuring the Power of DFG-Based Mid-Infrared
As shown in Figure 9, a large amount of thermal radiation is released from the back direction of the long-pass filter due to the strong absorption of pump and signal lights by the filter. Therefore, when the MIR power is measured by a thermal power meter behind the filter, the thermal radiation released from the filter will interfere with the measurement results, making the measured results much larger than the actual power. Under the condition of our experiment, only one-third to one-fifth of the total measured power shown on the display of the power meter is the real power of the coherent mid-infrared source generated by DFG.
Figure 9. Heat effect during the DFG process.
In view of this, we have explored three measurement methods for measuring the power of DFG-based MIR will be tested here, and their pros and cons are also weighed. Before all the measurement work, we tweak and adjust a all components to achieve the highest MIR power will be tweaked and adjusted.
Method 1: Nonlinear Crystal Rotating Method. Firstly, the total power of the MIR and the heat radiation by the filter is measured behind the filter. Then, the residual power (i.e., heat radiation) is measured when the output MIR disappears completely due to phase mismatch by slightly rotating the crystal by about 3°. Consequently, the power difference between the two measurements can be considered the actual MIR power. This method could achieve multiple measurements in a short time, since it is easy to restore the crystal to the position of the optimal phase-matching angle. The tweaked angle should not be more than 3° because the refraction of the crystal may affect the optical path.
Method 2: Pulse Synchronization Control method. As shown in Figure 7, the synchronism of the two pulses can be adjusted by the optical delay line placed behind the grating pulse compressors. When the two pulses are synchronized, the MIR power is the highest. We fFirst measure the power in this state. On the contrary, if the optical delay line is adjusted so that one pulse is far away from the other, the MIR will disappear. Then the MIR power could be determined by measuring the power difference between the two states. However, because the adjustment of the delay line is too fine, it is difficult to tweak back to the initial optimized position after one-time measurement, which means it will be hard to measure the power several times in a short time.
Method 3: Seed Laser Blocking Method. We first measure the power after the filter as we did in method 1, and then use an opaque plate to block the output of the femtosecond seed laser. When the seed laser is blocked, only ASE exists in the dual-wavelength CPA system. In this state, MIR will not be generated anymore. Therefore, the MIR power can be determined by measuring the difference between the two states. Multiple measurements could also be achieved in a short time. However, because the total power of the ASE is not equal to that of the total amplified output power (the power of the ASE is usually lower than that of the power in the normal amplification state), the measurement results are relatively larger than the actual one.
The measurement results by the three methods are shown in
Table 2 during
our experiment
[140][46], which is the average value obtained by five-times repeating the measurement for each method. It can be seen from the table that the power value and standard deviation measured by method 3 are large due to the thermal effect of the filter, while the results measured by methods 1 and 2 are relatively close. Compared to method 1, a relatively lower result and larger deviation are obtained by method 2 due to the difficulty of tweaking the delay line back. Consequently, method 1 is strongly recommended in the power measurement of DFG-based MIR due to its advantage of simplicity and time-saving.
Table 2. Measured results by three different methods.
5. Methods for Obtaining a High Beam Quality of DFG-Based Mid-Infrared
The beam quality of DFG-based MIR is also important since there are some application fields requiring high-quality MIR beams. The phase-matching condition largely affects the beam quality of DFG-based MIR, and the focused beam radius on the crystal and the crystal temperature are the main factors determining the phase-matching condition. Thus, there are mainly two ways to obtain a high beam quality of the MIR when applying the DFG method: controlling the focused beam size on the crystal or the crystal temperature.
With the tightly focused Gaussian beams, the optimum phase matching does not occur when
Δk=0 (
Δk=0 is the phase-mismatch vector) as in the plane wave approximation. As the focal spot size decreases to the dimensions of the MIR generated wavelengths, the axial phase matching is not optimized and the beam begins to propagate as conical emission, which will lower the beam quality. Consequently, a large-size focused beam with a radius of more than 100 μm on the crystal can ensure a better MIR beam quality
[140,143][46][49].
The crystal temperature is another parameter to determine the MIR beam quality, and there have been some mode-control DFG experiments in the quasi-phase-matched periodically poled lithium niobate (PPLN) nonlinear crystals with different temperatures
[144[50][51],
145], as shown in
Figure 10. Because the refractive index is highly related to the temperature in the PPLN crystal, the MIR beam profile is presented with different shapes in thermal imaging as the phase-mismatch vector varies with the temperature. Thus, it can be seen that temperature plays an important role in the MIR beam shaping under the condition of periodically poled crystals. However, there is almost no temperature-dependent research in the long-wavelength MIR region with other kinds of nonlinear crystals.
Figure 10. Experimental recordings and computations of the beam profile for three temperatures at normal incidence. Spatial mode control of radiation generated by DFG in periodically poled crystals.
Considering the currently reported works of the DFG-based MIR beam in the long-wavelength MIR region, thermal radiation is not a very big problem if the focused beam is large enough as most of the research group did not use a temperature-stabilized crystal oven. However, the thermal conductivity of the crystal is still an important parameter if the input peak intensity on the crystal reaches a relatively high value. According to reference
[103][9], the thermal conductivity of LiGaS
2 is 3 to 5 times higher than that of AgGaS
2 and AgGaSe
2 and ~2.5 times higher than the thermal conductivity of layered GaSe along its optic axis, and the thermal conductivity of LiGaS
2 is 10.05 W/(m·k)
[146][52]. However, the thermal conductivity of OP-GaP crystal is 110 W/(m·k) according to reference
[116[22][41],
135], which is the highest value
to the best of our knowledge. Thus, there is no doubt that OP-GaP is the preferable crystal for DFG-based MIR generation experiment if
we only take its thermal conductivity into consideration.