Terahertz (THz) generation by a GaP ridge waveguide with a collinear modal phase-matching scheme based on cascaded difference frequency generation (DFG) processes is theoretically analyzed. The cascaded Stokes interaction processes and the cascaded anti-Stokes interaction processes are investigated from coupled wave equations. THz intensities and quantum conversion efficiency are calculated. Compared with non-cascaded DFG processes, THz intensities from 11-order cascaded DFG processes are increased to 5.48. The quantum conversion efficiency of 177.9% in cascaded processes can be realized, exceeding the Manley-Rowe limit.
I. INTRODUCTION
The terahertz (THz) radiation, which is generally referred to as the frequency from 0.1 to 10 THz, has recently drawn much attention due to its tremendous potential applications, such as imaging, material detection, environmental monitoring, communication, astronomy and national defense security
[1
-
4]
. For such applications, a high-power, widely tunable, and compact source of THz-wave is required. Due to the interest in exploiting this region there have been many schemes proposed on source technologies over the last twenty years or so
[5
-
10]
. Among many electronic and optical methods for the coherent THz-wave generation, difference frequency generation (DFG)
[11
-
14]
is of importance because it offers the advantages of relative compactness, narrow linewidth, wide tuning range, high-power output and room-temperature working environment. In DFG, two optical pump beams, with their frequencies separated by a few THz, interact through a χ
(2)
process to generate a THz beam. Unfortunately, the quantum conversion efficiency of the DFG is extremely low as the THz-wave is intensely absorbed by the nonlinear optical crystal. To improve the low quantum conversion efficiency and overcome the Manley-Rowe limit, cascaded DFG in which more than one THz photon is generated from the depletion of a single pump photon is a promising method. Theoretical descriptions and experimental demonstrations of an enhancement output of THz wave via cascaded DFG processes have been reported recently. Liu
et al
.
[15]
proposed a scheme for monochromatic THz generation via cascading enhanced Cherenkov-type DFG in a sandwich-like waveguide. It is predicted that THz power can be boosted by nearly 8-fold with a 400 MW/cm pump in a 40-mm-long Si-LiNbO
3
-Si waveguide. Lee
et al
.
[16]
experimentally observed the fourth-Stokes order in cascaded Stimulated Polariton Scattering utilizing Mg:LiNbO
3
. Saito
et al
.
[17]
described a scheme for efficient THz generation using a cascaded optical parametric oscillator using a GaP sheet cavity. By choosing an appropriate pump wavelength and cavity design, the cascading process contributes to efficient THz-wave generation, resulting in a high output peak power of 1.8 MW and a high photon conversion efficiency of 1.086 at 3 THz.
The conversion efficiency of the cascaded DFG process is primarily determined by the effective interaction length and the absorption in the nonlinear optical crystal at THz frequencies. Usually, collinear phase matching is a preferred scheme to maximize the effective interaction length. THz wave generations by GaP rib waveguide via collinear modal phase-matched DFG have been observed
[18]
. While GaP has a large optical nonlinearity for THz DFG (50 pm/V at 1.55 μm
[19]
), its absorption coefficient in the THz range is relative large and increases rapidly with increasing frequency. Despite the fact that the THz wave is intensely absorbed by the optical crystal in collinear phase matching configuration, co-propagating the THz wave with the pump and signal waves in the same direction can achieve long- range amplification of the THz wave in a nonlinear crystal until pump depletion
[20]
.
In this paper, we present the theoretical analysis of THz generation by GaP ridge waveguide with a collinear modal phase-matching scheme based on cascaded DFG processes. We investigate the cascaded Stokes interaction processes and the cascaded anti-Stokes interaction processes. THz intensities and quantum conversion efficiency are calculated from coupled wave equations.
II. THEORETICAL MODEL
Figure 1
shows a schematic diagram of THz wave generation by collinear modal phase-matching cascaded DFG. TM-like guided THz wave in the GaP ridge waveguide was generated through type-I phase matching when the electric fields of both pump and signal waves are along [
1
10]. THz wave (
ω
T
) is generated via interactions between the incident pump (
ω
p
) and signal (
ω
s
) waves in the first-order DFG process, which consumes the higher frequency pump photon and amplifies the lower frequency signal photon. The amplified signal wave also acts as a higher frequency pump wave, which amplifies the THz wave and generates a new lower frequency cascaded signal (
ω
cs
) wave in the second-order DFG process. Simultaneously, anti-Stokes interactions will also occur that consume the THz photon and pump photon, resulting in a higher frequency anti-Stokes signal (
ω
cp
) wave. The cascaded Stokes processes and anti-Stokes processes can be continued to any high order as long as the phase-matching conditions are satisfied. The intensity of the THz wave is determined by a trade-off between the Stokes processes and the anti-Stokes processes.
Schematic diagram of the cascaded DFG to generate THz radiation by GaP ridge waveguide. DFG with Type-I phase matching condition, the electric field vectors of pump and signal wave Ep and Es are along [110], and that of THz wave ET is along [001] .
THz wave is guided by a GaP ridge waveguide with dimensions
t
,
h
, and
w
, as shown in
Fig. 1
. Single mode operation of the THz wave produced in the ridge waveguide is realized by satisfying the following condition
[21]
:
The collinear modal phase-matching condition in the ridge waveguide is expressed as
where
ni
(
i
=
p
,
s
and
T.eff
) correspond to the refractive index of the pump, signal and the TM-like guided THz wave, respectively.
λi
(
i
=
p
,
s
and
T
) correspond to the wavelength of the pump, signal and the TM-like guided THz wave, respectively. The mode effective index
nT.eff
of the THz wave is calculated using the effective-index method
[22]
.
The coupled wave equations of cascaded DFG can be derived from common nonlinear optical three-wave interaction equations, shown as
where
ωn
and
ωT
denote the frequency of pump and THz wave, respectively.
En
and
ET
denote the electric field amplitude of pump and THz wave, respectively.
αn
and
αT
denote the absorption coefficient of pump and THz wave in the optical crystal, respectively.Δ
kn
indicates the wave vector mismatch in the cascaded DFG process,
κn
and
κT
are the coupling coefficients.
deff
is the effective nonlinear coefficient, and the
deff
is 50 pm/V at 1.55 μm
[19]
.
c
is the speed of light in vacuum,
ε
0
is the vacuum dielectric constant,
I
is the power density,
nn
is the refractive index. The generation and consumption of THz photons are accomplished during the interaction between the
n
-order and (
n
+1)-order Stokes waves, as shown in Eq. (3). The second item in the right side of the equal sign in Eq. (4) shows the Stokes processes where THz photons and
n
-order Stokes photons are generated, and the third item in the right side of the equal sign in Eq. (4) shows the anti-Stokes processes where THz photons and (
n
+1)-order Stokes photons are consumed.
The theoretical values of refractive index are calculated using a wavelength-independent Sellmeier equation for GaP in the IR
[23]
and THz
[24]
range, respectively. The Sellmeier equation for GaP of Madarasz
et al
.
[23]
in the IR range can be written as
The Sellmeier equation for GaP in the THz range
[24]
can be written as
where
ε
(
v
) is the complex dielectric constant,
ν
is the wave-number,
ε
∞ is the high-frequency dielectric constant,
νTO
,
ρ
and
γ
are the eigenfrequency, oscillator strength and damping coefficient of the 367 cm
−1
polariton mode in GaP, and
ε
∞= 9.07,
ρ
= 1.945,
νTO
= 367 cm
−1
,
γ
=9.0 cm
−1
.
III. CALCULATIONS
Here, in simulating the cascaded DFG dynamics, pump wave
ω
p
and signal wave
ω
s
are supposed to be 193.55 and 192.55 THz, respectively. THz frequency
ω
T
is taken to be 1 THz. We set the dimensions of GaP ridge waveguide
t
,
h
, and
w
to 120, 200 and 80 μm, respectively. The mode effective index
nT.eff
for 1 THz is 3.1405 to realize collinear modal phase-matching DFG in GaP ridge waveguide. The wave vector mismatch Δ
k
and coherence length in cascaded DFG processes is shown in
Fig. 2
. In the cascaded Stokes processes, wave vector mismatch is less than 3.14 cm
−1
during 11-order cascaded processes. In the case of cascaded anti-Stokes processes, wave vector mismatch is less than 3.14 cm
−1
during 10-order cascaded processes. In the following calculations, 11-order cascaded Stokes and 10-order cascaded anti-Stokes processes are taken into account as the coherence length is larger than 1 cm. As shown in
Fig. 3
, THz intensities in GaP ridge waveguide based on cascaded DFG with cascading orders 1, 3, 5, 10 and 11 versus crystal length are calculated according to Eqs. (3) and (4). The intensity of both pump and signal wave are 20 MW/mm
2
. The absorption coefficient at 1 THz is 2.5 cm
−1
[25]
. From
Fig. 3
we find that THz intensities without cascading processes are extremely low. THz intensities with cascading order 3, 5, 10 and 11 are enhanced. THz intensity of 0.37 MW/mm
2
can be obtained with 11-order cascaded Stokes processes. Compared with non-cascaded DFG processes, THz intensities from 11-order cascaded DFG processes are increased to 5.48. In non-cascaded DFG processes, at best, a single THz photon is generated from each pump photon. The cascaded processes can enhance the THz output, simply by generating several THz photons from each pump photon.
Wave vector mismatch and coherence length of cascaded DFG. Assuming ωp =193.55 THz, ωs =192.55 THz, ωT = 1 THz.
THz intensities by GaP ridge waveguide based on cascaded DFG with cascading orders 1, 3, 5, 10 and 11.
As the Stokes processes generate THz photons and the anti-Stokes processes consume THz photons, THz intensities depend on the Stokes processes and the anti-Stokes processes.
Figure 4
shows the maximum intensities of the optical waves during the cascaded Stokes processes and anti-Stokes processes. In this figure we assume that the optical waves at interval of 1 THz with frequencies from 183.55 to 193.55 THz interact in the Stokes and anti-Stokes processes. The initial pump and signal waves are 188.55 and 187.55 THz, respectively, with a power density of 20 MW/mm
2
. From the figure we find that the power densities of optical waves in the Stokes processes is higher than that of optical waves in the anti-Stokes processes, which indicates that the Stokes processes are stronger than the anti-Stokes processes. In the cascaded Stokes processes,
n
-order Stokes photons are consumed and (
n
+1)-order Stokes and THz photons are amplified. The cascaded Stokes processes continue to amplify the (
n
+1)-order Stokes and THz photons as long as the phasematching conditions are satisfied. Actually, the anti-Stokes processes take place only if the Stokes processes generate THz photons. The anti-Stokes processes consume high-order Stokes and THz photons, and the processes will stop if the high-order Stokes and THz photons are exhausted.
The maximum intensity of the optical waves during the cascaded Stokes processes and anti-Stokes processes. Assuming the initial pump and signal waves are 188.55 and 187.55 THz with a power density of 20 MW/mm2, respectively.
Figure 5
shows the relationship between the maximum THz intensities and the pump wave frequencies. In this figure we assume that the optical waves at interval of 1 THz with frequencies from 183.55 to 193.55 THz interact in the cascaded Stokes and anti-Stokes processes. The frequency of the pump wave is 1 THz larger than that of the signal wave. Both the pump and signal intensities are 20 MW/mm
2
. From the figure we find that THz intensities are higher as the pump frequencies locate in the high-frequency area. The high THz intensities originate from the interaction of the high-order Stokes processes as the pump frequencies locate in the high-frequency area, which indicates that the Stokes processes are stronger than the anti-Stokes processes. As the pump frequency equals 193.55 THz, THz wave with a maximum intensity of 0.3686 MW/mm
2
can be obtained. In the high-frequency area where high-order Stokes processes interact, optimal crystal lengths are longer considering cascading, which is consistent with the principle of cascaded nonlinear processes.
The relationship between the maximum THz intensities and pump frequencies. Both the pump and signal intensities are 20 MW/mm2.
Pump intensity is directly related to the quantum conversion efficiency in a cascaded DFG processes. The maximum THz intensity and quantum conversion efficiency are calculated when the original pump intensities are changed from 1 MW/mm
2
to 20 MW/mm
2
, as shown in
Fig. 6
. In the calculations, pump wave and signal wave are supposed to be 193.55 and 192.55 THz, respectively.
Fig. 6
demonstrates that the maximum THz intensity and quantum conversion efficiency significantly increase with the pump intensity. THz wave with a maximum intensity of 0.3686 MW/mm
2
can be generated as pump intensity equals to 20 MW/mm
2
, corresponding to the quantum conversion efficiency of 177.9%. The quantum conversion efficiency of 177.9% in cascaded processes exceeds the Manley-Rowe limit.
The maximum THz intensity and quantum conversion efficiency versus pump intensity. Assuming pump and signal frequencies are 193.55 and 192.55 THz, respectively.
THz generation with the frequency of 1 THz with a collinear modal phase-matching scheme based on cascaded DFG processes is theoretically analyzed above. As for THz waves with frequencies lower than 1 THz, THz generation can be effectively enhanced based on cascaded DFG processes if the collinear modal phase-matching scheme is satisfied. As for THz waves with frequencies larger than several THz, the output of the THz wave is seriously affected by the absorption by the GaP crystal as the absorption coefficients of THz wave in GaP crystal rapidly increase with the increase of frequency.
IV. CONCLUSION
THz generation by GaP ridge waveguide with a collinear modal phase-matching scheme based on cascaded DFG processes is theoretically analyzed. The cascaded DFG processes comprise the Stokes interaction processes and the cascaded anti-Stokes interaction processes. The calculation results indicate that the Stokes processes are stronger than the anti-Stokes processes. Compared with non-cascaded DFG processes, THz intensities from 11-order cascaded DFG processes are increased to 5.48. THz wave with a maximum intensity of 0.3686 MW/mm
2
can be generated as pump intensity is 20 MW/mm
2
, corresponding to the quantum conversion efficiency of 177.9%. The quantum conversion efficiency of 177.9% exceeds the Manley-Rowe limit, providing us an efficient way to enhance the output of THz waves.
Color versions of one or more of the figures in this paper are available online.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant Nos. 61201101 and 61205003) and Young Backbone Teachers in University of Henan Province (Grant No. 2014GGJS-065).
Koenig S.
,
Lopez-Diaz D.
,
Antes J.
,
Boes F.
,
Henneberger R.
,
Leuther A.
,
Tessmann A.
,
Schmogrow R.
,
Hillerkuss D.
,
Palmer R.
,
Zwick T.
,
Koos C.
,
Freude W.
,
Ambacher O.
,
Leuthold J.
,
Kallfass I.
(2013)
“Wireless sub-THz communication system with high data rate,”
Nature Photon.
7
977 -
981
DOI : 10.1038/nphoton.2013.275
Watts C. M.
,
Shrekenhamer D.
,
Montoya J.
,
Lipworth G.
,
Hunt J.
,
Sleasman T.
,
Krishna S.
,
Smith D. R.
,
Padilla W. J.
(2014)
“Terahertz compressive imaging with metamaterial spatial light modulators,”
Nature Photon.
8
605 -
609
DOI : 10.1038/nphoton.2014.139
Köhler R.
,
Tredicucci A.
,
Beltram F.
,
Beere H. E.
,
Linfield E. H.
,
Davies A. G.
,
Ritchie D. A.
,
Iotti R. C.
,
Rossi F.
(2002)
“Terahertz semiconductor-heterostructure laser,”
Nature
417
156 -
159
DOI : 10.1038/417156a
Yeh K. L.
,
Hoffmann M. C.
,
Hebling J.
,
Nelson K. A.
(2007)
“Generation of 10 μJ ultrashort terahertz pulses by optical rectification,”
Appl. Phys. Lett.
90
171121 -
DOI : 10.1063/1.2734374
Carr G. L.
,
Martin M. C.
,
McKinney W. R.
,
Jordan K.
,
Neil G. R.
,
Williams G. P.
(2002)
“High-power terahertz radiation from relativistic electrons,”
Nature
420
153 -
156
DOI : 10.1038/nature01175
Ding Y. J.
(2004)
“Quasi-single-cycle terahertz pulses based on broadband-phase-matched difference-frequency generation in second-order nonlinear medium: high output powers and conversion efficiencies,”
IEEE J. Select. Topics Quantum Electron.
10
1171 -
1179
DOI : 10.1109/JSTQE.2004.837205
Knap W.
,
Lusakowski J.
,
Parenty T.
,
Bollaert S.
,
Cappy A.
,
Popov V. V.
,
Shur M. S.
(2004)
“Terahertz emission by plasma waves in 60 nm gate high electron mobility transistors,”
Appl. Phys. Lett.
84
2331 -
2333
DOI : 10.1063/1.1689401
Ding Y. J.
(2014)
“Progress in terahertz sources based on difference-frequency generation,”
J. Opt. Soc. Am. B
31
2696 -
2711
DOI : 10.1364/JOSAB.31.002696
Majkić A.
,
Zgonik M.
,
Petelin A.
,
Jazbinšek M.
,
Ruiz B.
,
Medrano C.
,
Günter P.
(2014)
“Terahertz source at 9.4 THz based on a dual-wavelength infrared laser and quasi-phase matching in organic crystals OH1,”
Appl. Phys. Lett.
105
141115 -
DOI : 10.1063/1.4897639
Dolasinski B.
,
Powers P. E.
,
Haus J. W.
,
Cooney A.
(2015)
“Tunable narrow band difference frequency THz wave generation in DAST via dual seed PPLN OPG,”
Opt. Express
23
3669 -
3680
DOI : 10.1364/OE.23.003669
Saito K.
,
Tanabe T.
,
Oyama Y.
(2014)
“Design of a GaP/Si composite waveguide for CW terahertz wave generation via difference frequency mixing,”
Appl. Opt.
53
3587 -
3592
DOI : 10.1364/AO.53.003587
Liu P.
,
Xu D.
,
Yu H.
,
Zhang H.
,
Li Z.
,
Zhong K.
,
Wang Y.
,
Yao J.
(2013)
“Coupled-mode theory for Cherenkov-type guided-wave terahertz generation via cascaded difference frequency generation,”
IEEE J. Lightwave Technol.
31
2508 -
2514
DOI : 10.1109/JLT.2013.2268876
Lee A. J.
,
Pask H. M.
(2015)
“Cascaded stimulated polariton scattering in a Mg:LiNbO3 terahertz laser,”
Opt. Express
23
8687 -
8698
DOI : 10.1364/OE.23.008687
Saito K.
,
Tanabe T.
,
Oyama Y.
(2015)
“Cascaded terahertz-wave generation efficiency in excess of the Manley-Rowe limit using a cavity phase-matched optical parametric oscillator,”
J. Opt. Soc. Am. B
32
617 -
621
DOI : 10.1364/JOSAB.32.000617
Saito K.
,
Tanabe T.
,
Oyama Y.
,
Suto K.
,
Nishizawa J.
(2009)
“Terahertz-wave generation by GaP rib waveguides via collinear phase-matched difference-frequency mixing of near- infrared lasers,”
J. Appl. Phys.
105
063102-1 -
063102-2
Shoji I.
,
Kondo T.
,
Kitamoto A.
,
Shirane M.
,
Ito R.
(1997)
“Absolute scale of second-order nonlinear-optical coefficients,”
J. Opt. Soc. Am. B
14
2268 -
2294
DOI : 10.1364/JOSAB.14.002268
Wang T. D.
,
Huang Y. C.
,
Chuang M. Y.
,
Lin Y. H.
,
Lee C. H.
,
Lin Y. Y.
,
Lin F. Y.
,
Kitaeva G. K.
(2013)
“Long-range parametric amplification of THz wave with absorption loss exceeding parametric gain,”
Opt. Express
21
2452 -
2462
DOI : 10.1364/OE.21.002452
Xia J.
,
Yu J.
,
Li Y.
,
Chen S.
(2004)
“Single-mode condition for silicon rib waveguides with large cross sections,”
Opt. Eng.
43
1953 -
1954
DOI : 10.1117/1.1776188
Yeh P.
,
Yariv A.
2007
Photonics: Optical Electronics in Modern Communication
Oxford University Press
New Delhi
Bass M.
,
Enoch J. M.
,
Van Stryland E. W.
,
Wolfe W. L.
2001
Handbook of Optics
2nd ed
McGraw-Hill
New York, USA
Palik E. D.
1998
Handbook of Optical Constants of Solids
Academic Press
San Diego, USA
Tanabe T.
,
Suto K.
,
Nishizawa J.
,
Saito K.
,
Kimura T.
(2003)
“Tunable terahertz wave generation in the 3- to 7-THz region from GaP,”
Appl. Phys. Lett.
83
237 -
239
DOI : 10.1063/1.1592889