Simultaneous Generation of Two Pairs of Stokes and Terahertz Waves from Coupled Optical Parametric Oscillations with Quasi-Phase-Matching

We present a theoretical investigation on simultaneous generation of two pairs of Stokes and terahertz (THz) waves from coupled optical parametric oscillations (OPOs) with a quasi-phase-matching (QPM) scheme. The two pairs of Stokes and THz waves are generated by stimulated polariton scattering (SPS) from periodically-inverted GaP. By analyzing the QPM conditions of coupled OPOs we find that the two THz waves with any frequency below the transverse optical (TO) mode frequency of GaP can be simultaneously generated with a suitable pump wavelength. We calculate the photon flux densities of the two THz waves by solving the coupled wave equations. The calculation results indicate that the two THz waves can be efficiently generated with high pump intensities, particularly in lower THz frequency band.


Introduction
In recent years, it has been demonstrated that coupled optical parametric oscillations (OPOs) can simultaneously generate two pairs of signal and idler waves from periodically-inverted KTiOPO 4 (KTP) plates by a single pump wavelength [1,2].The wavelengths for each pair of signal and idler waves satisfy quasi-phase-matching (QPM) conditions.The wavelengths of the signal and idler waves are near infrared when pump wavelength is 532 nm.Based on the research above, we consider that the frequencies of the two idler waves generated from coupled OPOs with periodically-inverted GaP by stimulated polariton scattering (SPS) can extend to terahertz (THz) frequencies.The frequency separation between the two THz waves can be tuned by selecting the thickness of each GaP plate.The generated two THz waves are useful for imaging and food inspection [3,4].
SPS which was used in MgO: LiNbO 3 , LiTaO 3 , and KTiOPO 4 , has proved to be an efficient scheme to generate THz waves [5][6][7][8][9][10][11][12].SPS can also be used in zinc-blende structure crystals, since the crystals have a high nonlinear coefficient and a low THz wave absorption coefficient [13,14].GaP which has an infrared-and Raman-active transverse optical (TO) mode with frequency of 367 cm −1 is an attractive material for THz wave generation via SPS [15] since GaP has a high second-order nonlinear coefficient (d 36 = 70.6 pm/V at 1064 nm) [16] and a low THz wave absorption coefficient below TO mode frequency [15].Moreover, collinear configuration can be realized in GaP by using cross-Reststrahlen band dispersion compensation phase-matching [17].
In this work, we theoretically study simultaneous generation of two pairs of Stokes and THz waves from coupled OPOs with periodically-inverted GaP by SPS.The QPM conditions of the coupled OPOs are analyzed.We calculate the photon flux densities of the two THz waves by solving the coupled wave equations.

Phase-Matching Characteristics
In the SPS processes with cross-Reststrahlen band dispersion compensation phase-matching, pump and Stokes waves are in the near-infrared transmission window of GaP crystal, and THz wave is in the far-infrared transmission window, on the other side of the crystal's Reststrahlen band.Since the refractive indices of the three waves in GaP are approximately equal, the collinear phase-matching condition can be realized.Figure 1 shows THz frequencies ν T versus pump wavelengths λ p with collinear phase-matching condition in GaP.The theoretical values of the refractive index are calculated using a Sellmeier equation for GaP in the infrared range [17] and in the THz range [15], respectively.From the figure we find that pump wavelengths λ p in the range of 0.69-1.03µm fulfill collinear phase-matching, corresponding to the THz frequencies ν T in the range of 0.15-8.8THz. Figure 2 shows phase mismatching ∆k among collinear SPS processes in GaP when λ p = 0.9 µm.
where k p , k s , and k T are the wave vectors of pump, Stokes, and THz waves respectively.From the figure we find that ∆k is zero at 5.65 THz, where phase-matching can realize.At area I where frequencies are slightly smaller than 5.65 THz, ∆k is slightly larger than zero, whereas at area II where frequencies are slightly larger than 5.65 THz, ∆k is slightly smaller than zero.At area I the phase mismatch ∆k in optical parametric oscillation (OPO) can be compensated by grating vector k Λ of periodically-inverted GaP.k Λ = 2π/Λ and Λ is the QPM period.Simultaneously, at area II, the phase mismatch ∆k in the other OPO can also be compensated by the same grating vector k Λ .As a result, the two phase mismatch of the two OPOs can be simultaneously compensated by the same grating vector k Λ .The two OPOs will couple with each other with a same pump wavelength and a same QPM period.In order to achieve efficient conversion of coupled OPOs from pump wave to THz wave, a precise phase-matching condition must be satisfied.Figure 3 shows the QPM scheme of coupled OPOs with periodically-inverted GaP.In order to simultaneously satisfy the two QPM conditions, accurate value of grating vector k Λ is required.The two QPM conditions of the coupled OPOs are where k s1 and k s2 are the wave vectors of the two Stokes waves respectively, k T1 and k T2 are the wave vectors of the two THz waves respectively.The energy conservation condition has to be fulfilled, where λ s1 and λ s2 are the wavelengths of the two Stokes waves respectively, λ T1 and λ T2 are the wavelengths of the two THz waves respectively.If the two QPM conditions of the coupled OPOs are simultaneously satisfied, the two THz waves can be simultaneously generated with a single pump wave.Figure 4 shows THz frequencies and Stokes wavelengths versus QPM period Λ with coupled OPOs with different pump wavelengths.ν T1 and ν T2 are frequencies of the two THz waves,

THz Wave Photon Flux Density
The generated THz intensities by SPS processes depend on the nonlinear optical susceptibility and THz absorption coefficients.The nonlinear optical susceptibility in GaP in the THz frequency range is governed by the superposition of electronic and ionic contributions.Faust showed that the ionic and electronic contributions are of opposite sign, leading to a cancellation of both contributions below the TO mode frequency of GaP [18].Sussman deduced coupled wave equations for SPS processes [19].
this work we calculate THz intensity based on the theories of Faust and Sussman.The analytical expression of THz parametric gain coefficient g T in SPS processes under the QPM condition in the international system of units can be written as [19,20] ) where ω TO , S, and Γ denote eigenfrequency, oscillator strength of the polariton modes, and the bandwidth of the TO mode in GaP crystal, respectively.ω s and ω T are the angular frequencies of Stoke and THz waves, respectively.ε 0 and ε ∞ are vacuum dielectric constant and high-frequency dielectric constant, respectively.I p is the power density of pump wave, g 0 is the low-loss parametric gain.n p , n s , and n T are the refractive indices of the pump, Stokes and THz waves, respectively.α T is material absorption coefficient in THz region.d e is the electronic second-order nonlinear coefficient, and d e f f is the bulk value of the second-order nonlinear coefficient involving electronic and ionic contributions.C 1 is a coupling constant, and C 1 is −0.53 for GaP [18].THz wave absorption, without phase mismatch and pump depletion, the coupled wave equations can be solved to give the THz photon flux density φ T with a general solution [21,22], given by where L is the crystal length.The initial THz photon flux density φ T is assumed to be zero, and φ s (0) is the initial Stokes wave photon flux density.Figure 6 shows THz wave photon flux density φ T /φ s (0) from coupled OPOs by SPS processes.φ T1 and φ T2 are the photon flux densities of ν T1 and ν T2 , respectively, φ s1 (0) and φ s2 (0) are the initial photon flux densities of λ s1 and λ s2 , respectively.In the calculations pump intensity I p is 500 MW/cm 2 since the damage threshold of GaP is 650 MW/cm 2 [23].
From Figure 6a we find that φ T /φ s (0) rapidly increases with crystal length L as ν T1 = 1.13 THz and ν T2 = 0.95 THz.As shown in Figure 5, d e f f is large and α T is small as THz wave frequencies are around 1 THz.The large d e f f and small α T enhance the rapid increase of φ T /φ s (0) with crystal length L.
As the coefficient α T of 1.13 THz is larger than that of 0.95 THz, φ T1 /φ s1 (0) is smaller than φ T2 /φ s2 (0).φ T2 /φ s2 (0) with value of 0.016 can be reached as crystal length L is 5 cm.From Figure 6b we find that φ T /φ s (0) rapidly increases with crystal length L and then keeps invariant as ν T1 = 6.69 THz and ν T2 = 6.45 THz.Compared with Figure 6a, d e f f becomes smaller and α T becomes larger, which induces an intensive decrease of φ T1 /φ s1 (0) and φ T2 /φ s2 (0) in Figure 6b.From Figure 6 we conclude that high-power THz waves can be generated in lower THz frequency band.The intensity of THz wave can be boosted by injecting an initial Stokes wave as a seed light.As shown in Equation ( 6), low-loss parametric gain g 0 is proportional to the pump intensity I p .A pump wave with a pump intensity on the order of GW/cm 2 can enlarge the gain coefficient g T by several orders of magnitude [24,25], which can enhance the intensity of THz wave. Figure 7 shows the THz wave photon flux density φ T /φ s (0) versus pump intensity I p .φ T1 /φ s1 (0) and φ T2 /φ s2 (0) monotonously increase with the increase of pump intensity I p .φ T2 /φ s2 (0), with a value of 0.026, can be reached as pump intensity I p is 650 MW/cm 2 .The intense THz wave can be generated using a pump wave with a high intensity.

Conclusions
Two pairs of Stokes and THz waves can be simultaneously generated from coupled OPOs with SPS processes in periodically-inverted GaP.The two QPM conditions of the coupled OPOs are simultaneously satisfied.The two THz waves with any frequency below the TO mode frequency of GaP can be generated with a suitable pump wavelength.THz waves can be efficiently generated with high pump intensities, particularly in a lower THz frequency band.

Figure 1 .
Figure 1.THz frequencies ν T versus pump wavelengths λ p with collinear phase-matching condition in GaP.

Figure 3 .
Figure 3. Schematic diagram of QPM scheme of the coupled OPOs with periodically-inverted GaP.Pump wave, Stokes waves, and THz waves are collinear.k p is the wave vector of pump wave and k s1 and k s2 are the wave vectors of the two Stokes waves respectively, k T1 and k T2 are the wave vectors of the two THz waves respectively, k Λ is the grating vector of periodically-inverted GaP.

Figure 5
shows calculated d e f f and α T versus THz frequencies.In the calculation d e is 70.6 pm/V with pump wavelength of 1.064 µm.From the figure we find that α T monotonously increases with frequencies, and d e f f smoothly decreases first and then rapidly increases.d e f f is smaller than 70.6 pm/V below the TO mode frequency of GaP due to the cancellation of ionic and electronic contributions.The minimum value of d e f f is 5.0 pm/V at 7.54 THz.When approaching TO mode frequency, d e f f is larger than 70.6 pm/V because polaritons at this area induce giant ionic nonlinearities.By comprehensive consideration of d e f f and α T we conclude that the lower frequencies of THz wave in the range of 0-7.54 THz, the higher intensities of THz wave generated from SPS processes.

Figure 5 .
Figure 5. Calculated d e f f and α T versus THz frequencies.d e is 70.6 pm/V with pump wavelength of 1.064 µm.
c is velocity of light.From the figure we find that two pairs of Stokes and THz waves can be generated by coupled OPOs with a single pump wave.As QPM period Λ changes the THz frequencies and Stokes wavelengths are tuned.The difference frequencies between the two THz waves become small as QPM period Λ increases.The minimum difference frequencies are smaller than 0.01 THz.As pump wavelengths change from 0.75 µm to 0.852, 0.95, and 1.025 µm, ν T1 and ν T2 decrease from 8.1 THz to 1.05 THz.From Figure4we conclude that the two THz waves with any frequency from 1.05 THz to 8.1 THz can be simultaneously generated by coupled OPOs with a suitable pump wavelength.