Generation of Higher Terahertz Harmonics in Nonlinear Paraelectrics under Focusing in a Wide Temperature Range ()

Volodymyr Grimalsky^{1}, Jesus Escobedo-Alatorre^{1}, Christian Castrejon-Martinez^{2}, Yered Gomez-Badillo^{1}

^{1}Center for Investigations on Engineering and Applied Science (CIICAp), Institute for Investigations on Basic and Applied Science (IICBA), Autonomous University of State Morelos (UAEM), Cuernavaca, Mor., Mexico.

^{2}Tecnológico Nacional de México/Instituto Tecnológico de Zacatepec, Zacatepec, Mor., Mexico.

**DOI: **10.4236/jemaa.2023.154004
PDF
HTML XML
103
Downloads
305
Views
Citations

It is theoretically investigated the generation of higher harmonics of two-dimensional and three-dimensional terahertz electromagnetic beams in nonlinear crystals. The attention is paid to crystalline paraelectrics like SrTiO_{3} under the temperatures 60 - 200 K, these crystals possess the cubic nonlinearity. The bias electric field is applied to provide the dominating quadratic nonlinearity. The initial focusing of the beams not only increases the efficiency of generation of higher harmonics, but alto makes possible to select maxima of different higher harmonics at some distances from the input. At lower temperatures the nonlinearity behaves at smaller input amplitudes, whereas at higher temperatures the harmonic generation can be observed at higher frequencies up to 1.5 THz. In three-dimensional beams the peak amplitudes of higher harmonics can be bigger than in two-dimensional beams, but the ratios of these peak values to the maximum values of the focused first harmonic are smaller than in two-dimensional beams.

Keywords

Terahertz Wave Beams, Nonlinear Crystalline Paraelectrics, Different Temperatures, Generation of Harmonics, Initial Focusing

Share and Cite:

Grimalsky, V. , Escobedo-Alatorre, J. , Castrejon-Martinez, C. and Gomez-Badillo, Y. (2023) Generation of Higher Terahertz Harmonics in Nonlinear Paraelectrics under Focusing in a Wide Temperature Range. *Journal of Electromagnetic Analysis and Applications*, **15**, 43-58. doi: 10.4236/jemaa.2023.154004.

1. Introduction

Now the assimilation of the terahertz (THz) range 0.1 - 30 THz takes place [1] - [17] . As the nonlinear crystals, the ferroelectrics in the non-polar phase are utilized as the nonlinear dielectrics in the lower part of THz range, so-called paraelectrics like SrTiO_{3}, KTaO_{3} [18] - [40] . The crystalline SrTiO_{3} possesses high cubic electrodynamic nonlinearity and low losses in the lower part of THz range 0.1 - 2.5 THz at moderately low temperatures *T* = 60 - 200 K. There exists the frequency dispersion in THz range when the frequency is near the soft mode frequency *i.e.* the lowest frequency of oscillations of the optical type of the crystalline lattice [19] [21] [22] [23] [24] [27] [28] [40] [41] [42] . In the crystalline SrTiO_{3} the soft mode frequency decreases with the decrease of the temperature whereas the dielectric nonlinearity increases there [19] . In the absence of the frequency dispersion in the microwave range, the nonlinearity results in the generation of higher harmonics and forming the shock electromagnetic (EM) waves [29] [30] . The waveguide dispersion can be used there [31] .

Because there is a problem of excitation of powerful THz radiation, the interesting and important phenomenon is generation of higher harmonics [33] [39] of input EM waves at relatively low frequencies. The transverse bias electric field should be applied to provide the induced quadratic nonlinearity [29] [33] . Earlier the generation of higher THz harmonics was investigated in 2D geometry, or for the beams that depend on a single transverse coordinate, at the temperature *T* ≈ 80 K [33] . Only the quadratic nonlinearity was taken into account there. Note that the frequencies of higher harmonics are limited by the soft mode frequency. At higher temperatures the soft mode frequency increases, whereas the nonlinearity decreases. Therefore, for the practical needs it is rather better to investigate the generation of higher harmonics in the wide temperature interval. In THz range the dissipation in the paraelectrics increases compared with the microwave range. A possible way to compensate losses is using the initial focusing of the EM wave beams [33] [34] [36] [37] [39] . Thus, also 3D THz beams should be considered to produce better initial focusing. Under maximum concentration of EM energy in 3D case, the role of the cubic nonlinearity may be essential compared with 2D case, as well as the self-action due the generation of the zeroth harmonic. The cubic nonlinearity results in both the influence on generation of higher harmonics and in the self-action.

The present paper is devoted to the theoretical investigations of generation of higher harmonics of THz EM waves in paraelectric crystals like SrTiO_{3} in the wide temperature range *T* = 60 - 200 K. The bias electric field is applied to provide the dominating quadratic nonlinearity. The terms due to the cubic nonlinearity are taken into account, too. Both the focusing of 2D, or plane, beams and 3D, or cylindrical, ones are considered. At lower temperatures the dielectric nonlinearity is higher, but the possible frequency range is lower due to decreasing the soft mode frequency. At higher temperatures, it is possible to increase the frequency range. It is demonstrated that the cubic nonlinearity is not important under focusing of 2D beams, but it is important under extreme focusing of 3D beams.

2. Basic Equations

The nonlinear propagation of EM wave beams with the dominating electric field component *E _{y}* =

$\begin{array}{l}\frac{{\partial}^{2}P}{\partial {t}^{2}}+\gamma \frac{\partial P}{\partial t}+{\omega}_{T}^{2}\left(1+\frac{{P}^{2}}{{P}_{n}^{2}}\right)P=\epsilon \left(0\right){\omega}_{T}^{2}E;\\ \frac{{\partial}^{2}E}{\partial {z}^{2}}+{\Delta}_{\perp}E=\frac{1}{{c}^{2}}\frac{{\partial}^{2}P}{\partial {t}^{2}};\\ D\equiv {\epsilon}_{0}\left(E+P\right)\approx {\epsilon}_{0}P;\\ {\Delta}_{\perp}E\equiv \{\begin{array}{l}\frac{{\partial}^{2}E}{\partial {x}^{2}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\text{2D}\text{\hspace{0.17em}}\text{case},\\ \frac{1}{\rho}\frac{\partial}{\partial \rho}\left(\rho \frac{\partial E}{\partial \rho}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\text{3D}\text{\hspace{0.17em}}\text{case};\end{array}\\ \rho \equiv {\left({x}^{2}+{y}^{2}\right)}^{1/2}.\end{array}$ (1)

Below both 3D beams and 2D ones are considered. Here
${\omega}_{T}$ is the soft mode frequency, which is in THz range for SrTiO_{3},
$\gamma $ is the lattice dissipation. At the temperature *T* ≈ 80 K it is
${\omega}_{T}\approx 6\times {10}^{12}\text{\hspace{0.17em}}{\text{s}}^{-1}$ , and the static linear dielectric permittivity is
$\epsilon \left(0\right)\equiv \epsilon \left(\omega =0\right)=1.8\times {10}^{3}$ [19] [26] [33] . In paraelectric crystalline SrTiO_{3} the permittivity increases with the decrease of temperature. In Equations (1) the parameter *P _{n}* determines the value of the cubic nonlinearity of the polarization. It is connected with the characteristic magnitude of the electric field

Figure 1. The temperature dependencies of the static linear permittivity *ε*(0), curve 1; of the soft mode frequency *ω _{T}* (10

For 3D, or cylindrical, EM beams the validity of the paraxial approximation is assumed and checked [37] , namely the possible focusing is moderate and the longitudinal component of the electric field is small $\left|{E}_{z}\right|\ll \left|{E}_{y}\right|$ .

In the linear case Equations (1) result in the well-known expressions of the complex dielectric permittivity *ε*(*ω*) º *ε'*(*ω*) + *iε**''*(*ω*) and to the dispersion relation for the plane EM wave *k* = *k*(ω) [19] [26] [33] :

$\begin{array}{l}\epsilon \left(\omega \right)\approx \frac{\epsilon \left(\omega =0\right)\cdot {\omega}_{T}^{2}}{{\omega}_{T}^{2}-{\omega}^{2}+i\gamma \omega},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\omega >0;\\ k\left(\omega \right)\equiv {k}^{\prime}+i{k}^{\u2033}=\frac{\omega}{c}\epsilon {\left(\omega \right)}^{1/2};\text{\hspace{0.17em}}\text{\hspace{0.17em}}E~\mathrm{exp}\left(i\left(\omega t-kz\right)\right).\end{array}$ (2)

In the nonlinear stationary case Equations (1) result in the known formula for the nonlinear permittivity [19] [26] :

$\epsilon \left(E\right)=\frac{\epsilon \left(\omega =0\right)}{1+\frac{{P}^{2}}{{P}_{n}^{2}}}\approx \frac{\epsilon \left(\omega =0\right)}{1+\frac{{E}^{2}}{{E}_{n}^{2}}}\approx \epsilon \left(\omega =0\right)\left(1-\frac{{E}^{2}}{{E}_{n}^{2}}\right).$ (3)

Equations (3) for the static nonlinear permittivity are equivalent in the case of moderate electric fields |*E*| < *E _{n}*, that is considered below.

3. Equations for Slowly Varying Amplitudes

Here the nonlinear EM wave propagation is investigated in the crystalline SrTiO_{3}, when the bias electric field *E _{s}* is applied [33] :

$\begin{array}{l}E={E}_{s}+\stackrel{\u02dc}{E},\text{\hspace{0.17em}}\text{\hspace{0.17em}}P={P}_{s}+\stackrel{\u02dc}{P};\\ {P}_{s}\cdot \left(1+\frac{{P}_{s}^{2}}{{P}_{n}^{2}}\right)=\epsilon \left(0\right){\omega}_{T}^{2}{E}_{s}.\end{array}$ (4)

In the absence of the bias field only the odd harmonics are excited, and the efficiency is low, as our simulations have been demonstrated. In the presence of the bias field the induced quadratic nonlinearity becomes dominating, but the cubic nonlinearity should be preserved generally in the equations. Generally all harmonics, even and odd ones, are essential. The equations for the variable components of the polarization $\stackrel{\u02dc}{P}$ and the electric field $\stackrel{\u02dc}{E}$ are:

$\begin{array}{l}\frac{{\partial}^{2}\stackrel{\u02dc}{P}}{\partial {t}^{2}}+\gamma \frac{\partial \stackrel{\u02dc}{P}}{\partial t}+{\stackrel{\u02dc}{\omega}}_{T}^{2}\stackrel{\u02dc}{P}+{\omega}_{T}^{2}\left(\frac{3{P}_{s}{\stackrel{\u02dc}{P}}^{2}}{{P}_{n}^{2}}+\frac{{\stackrel{\u02dc}{P}}^{3}}{{P}_{n}^{2}}\right)=\epsilon \left(0\right){\omega}_{T}^{2}\stackrel{\u02dc}{E};\\ \Delta \stackrel{\u02dc}{E}\equiv \frac{{\partial}^{2}\stackrel{\u02dc}{E}}{\partial {z}^{2}}+{\Delta}_{\perp}\stackrel{\u02dc}{E}=\frac{1}{{c}^{2}}\frac{{\partial}^{2}\stackrel{\u02dc}{P}}{\partial {t}^{2}};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\u02dc}{\omega}}_{T}^{2}\equiv {\omega}_{T}^{2}\left(1+\frac{3{P}_{s}{}^{2}}{{P}_{n}^{2}}\right).\end{array}$ (5)

From Equations (5) it is possible to obtain a single equation for the polarization [33] :

$\begin{array}{l}\Delta \left(\frac{{\partial}^{2}\stackrel{\u02dc}{P}}{\partial {t}^{2}}+\gamma \frac{\partial \stackrel{\u02dc}{P}}{\partial t}+{\stackrel{\u02dc}{\omega}}_{T}^{2}\stackrel{\u02dc}{P}\right)-\frac{{\omega}_{T}^{2}}{{c}^{2}}\epsilon \left(0\right)\frac{{\partial}^{2}\stackrel{\u02dc}{P}}{\partial {t}^{2}}\\ =-{\omega}_{T}^{2}\left(\frac{3{P}_{s}}{{P}_{n}^{2}}\Delta \left({\stackrel{\u02dc}{P}}^{2}\right)+\frac{1}{{P}_{n}^{2}}\Delta \left({\stackrel{\u02dc}{P}}^{3}\right)\right).\end{array}$ (6)

The nonlinearity is considered as moderate here, and the method of slowly varying amplitudes is applied [33] [43] [44] [45] .

Below the propagation of THz EM beams along *OZ* axis is considered. The solution of Equation (6) is searched as the set of harmonics including the zeroth harmonic:

$\begin{array}{l}\stackrel{\u02dc}{P}=\frac{1}{2}{\displaystyle \underset{j=1,2,3,\cdots}{\sum}{B}_{j}\left(z,{\stackrel{\to}{r}}_{\perp},t\right){\text{e}}^{i{\phi}_{j}}}+c.c.+{B}_{0}\left(z,{\stackrel{\to}{r}}_{\perp},t\right);\\ {\phi}_{j}\equiv {\omega}_{j}t-{k}^{\prime}\left({\omega}_{j}\right)z;\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\omega}_{j}=j\cdot {\omega}_{1};\\ k\left(\omega \right)\equiv {k}^{\prime}+i{k}^{\u2033}=\frac{\omega}{c}\epsilon {\left(\omega \right)}^{1/2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\epsilon \left(\omega \right)\approx \frac{\epsilon \left(\omega =0\right)\cdot {\omega}_{T}^{2}}{{\stackrel{\u02dc}{\omega}}_{T}^{2}-{\omega}^{2}+i\gamma \omega}.\end{array}$ (7)

Here *B _{j}*(

The frequency dependencies of the real parts of the permittivity, the real and imaginary parts of wave numbers of EM waves, and the ratios of the imaginary part to the real one of the permittivity are given in Figure 2. The bias field *E _{s}* = 0.3

In Equations (7) the zeroth harmonic term is taken into account that is due to the induced quadratic nonlinearity. From Equation (6) it is possible to obtain the following expression for *B _{0}*:

${\stackrel{\u02dc}{\omega}}_{T}^{2}{B}_{0}+\frac{3{P}_{s}}{2{P}_{n}^{2}}{\omega}_{T}^{2}{\displaystyle \underset{l=1,2,3,\cdots}{\sum}{\left|{B}_{l}\right|}^{2}}=0.$ (8)

In another words, the term with *B _{0}* < 0 results in the effective decreasing of the bias polarization

${P}_{s}\to {P}_{s}\cdot \left(1-\frac{3{\omega}_{T}^{2}}{2{P}_{n}^{2}{\stackrel{\u02dc}{\omega}}_{T}^{2}}{\displaystyle \underset{l=1,2,3,\cdots}{\sum}{\left|{B}_{l}\right|}^{2}}\right).$ (9)

With using the standard procedure described in details in [33] , the following nonlinear parabolic equations for harmonics have been derived:

$\begin{array}{l}\frac{\partial {B}_{j}}{\partial z}+\frac{i}{2{{k}^{\prime}}_{j}}{\Delta}_{\perp}{B}_{j}+{\Gamma}_{j}{B}_{j}\\ =\frac{i{{k}^{\prime}}_{j}{\omega}_{T}^{2}}{\left({\stackrel{\u02dc}{\omega}}_{T}^{2}-{\omega}_{j}^{2}+i\gamma {\omega}_{j}\right){P}_{n}^{2}}\{\frac{3{P}_{s}}{4}[{\displaystyle \underset{m}{\sum}{B}_{m}{B}_{j-m}}\mathrm{exp}\left(i\left({\phi}_{m}+{\phi}_{j-m}-{\phi}_{j}\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+2{\displaystyle \underset{m}{\sum}{B}_{m}^{\ast}{B}_{j+m}}\mathrm{exp}\left(i\left(-{\phi}_{m}+{\phi}_{j+m}-{\phi}_{j}\right)\right)+4{B}_{0}{B}_{j}]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{1}{8}[{\displaystyle \underset{m,p}{\sum}{B}_{j-m-p}{B}_{m}{B}_{p}}\mathrm{exp}\left(i\left({\phi}_{m}+{\phi}_{p}+{\phi}_{j-m-p}-{\phi}_{j}\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+3{\displaystyle \underset{m,p}{\sum}{B}_{m+p-j}^{\ast}{B}_{m}{B}_{p}}\mathrm{exp}\left(i\left({\phi}_{m}+{\phi}_{p}-{\phi}_{m+p-j}-{\phi}_{j}\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+3{\displaystyle \underset{m,p}{\sum}{B}_{m+p+j}{B}_{m}^{\ast}{B}_{p}^{\ast}}\mathrm{exp}\left(i\left({\phi}_{m}+{\phi}_{p}-{\phi}_{m+p-j}-{\phi}_{j}\right)\right)]\}.\end{array}$ (10)

(a)(b)

Figure 2. The dependencies on frequency *ω* of the parameters of THz waves in the crystalline SrTiO_{3}. Part (a) is the dependencies of the real parts of the linear permittivity *ε**'*, solid curves, and the wave numbers of EM wave *k**'*, dot curves. Part (b) is the dependencies of the imaginary part of the wave number *k**''*, dot curves, and tan(*δ*) º |*ε**''*|/*ε**'*, solid curves. The curves 1, 2, 3, 4 are at the temperatures *T* = 60 K, 77 K, 150 K, and 200 K correspondingly. The bias electric field is *E _{s}* = 0.3

For simulations it is better to rewrite Equations (10) with using the slowly varying amplitudes for harmonics *A _{j}* of the electric field

$\begin{array}{l}{\left({{k}^{\prime}}_{j}\right)}^{2}{A}_{j}\mathrm{exp}\left(-ij{{k}^{\prime}}_{1}z\right)=\frac{{\omega}_{j}^{2}}{{c}^{2}}{B}_{j}\mathrm{exp}\left(-i{{k}^{\prime}}_{j}z\right);\\ \text{or}\text{\hspace{0.17em}}\text{\hspace{0.05em}}{B}_{j}=\epsilon \left(0\right){\alpha}_{j}{A}_{j}\mathrm{exp}\left(i\left({{k}^{\prime}}_{j}-j{{k}^{\prime}}_{1}\right)z\right);\\ {\alpha}_{j}\equiv \frac{1}{\epsilon \left(0\right)}{\left(\frac{{{k}^{\prime}}_{j}c}{{\omega}_{j}}\right)}^{2}\approx 1.\end{array}$ (11)

From Equation (11) the following set of equation is written down:

$\begin{array}{l}\frac{\partial {A}_{j}}{\partial z}+\frac{i}{2{{k}^{\prime}}_{j}}{\Delta}_{\perp}{A}_{j}+i\left({k}_{j}-j{{k}^{\prime}}_{1}\right){A}_{j}\\ =i\frac{{{k}^{\prime}}_{j}{\omega}_{T}^{2}{\epsilon}^{2}\left(0\right)}{\left({\stackrel{\u02dc}{\omega}}_{T}^{2}-{\omega}_{j}^{2}+i\gamma {\omega}_{j}\right){P}_{n}^{2}{\alpha}_{j}}\{\frac{3{P}_{s}}{4\epsilon \left(0\right)}[{\displaystyle \underset{m}{\sum}{\alpha}_{j-m}{\alpha}_{m}}{A}_{j-m}{A}_{m}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+2{\displaystyle \underset{m}{\sum}{\alpha}_{m-j}{\alpha}_{m}}{A}_{m-j}^{\ast}{A}_{m}+4{A}_{0}{\alpha}_{j}{A}_{j}]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+\frac{1}{8}[{\displaystyle \underset{m,p}{\sum}{\alpha}_{j-m-p}}{\alpha}_{m}{\alpha}_{p}{A}_{j-m-p}{A}_{m}{A}_{p}+3{\displaystyle \underset{v,p}{\sum}{\alpha}_{m+p-j}}{\alpha}_{m}{\alpha}_{p}{A}_{m+p-j}^{\ast}{A}_{m}{A}_{p}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}+3{\displaystyle \underset{v,p}{\sum}{\alpha}_{j+m-p}}{\alpha}_{m}{\alpha}_{p}{A}_{j+m+p}{A}_{m}^{\ast}{A}_{p}^{\ast}]\};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{P}_{n}/\epsilon \left(0\right)\equiv {E}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{P}_{s}/\epsilon \left(0\right)\approx {E}_{s}.\end{array}$ (12)

In Equations (12) the expression for the zeroth harmonic *A _{0}* is

${A}_{0}\equiv \frac{{B}_{0}}{\epsilon \left(0\right)}=-\frac{3{P}_{s}\epsilon \left(0\right){\omega}_{T}^{2}}{2{P}_{n}^{2}{\stackrel{\u02dc}{\omega}}_{T}^{2}}{\displaystyle \underset{l=1,2,3,\cdots}{\sum}{\alpha}_{l}^{2}{\left|{A}_{l}\right|}^{2}}.$ (13)

It is assumed that the first harmonic only is at the input of the system, *ω _{1}* <

To increase the efficiency of generation of higher harmonics, the initial focusing by the circular antenna is applied, see Figure 3. The boundary condition for the first harmonic is in 2D case

${A}_{1}\left(z=0,x,t\right)={A}_{10}{F}_{x}\left(x\right)M\left(x\right).$ (14)

In 3D case there is the radial distance *ρ* instead *x*. Here *A _{10}* is the maximum value of the amplitude at the input

The factor *M*(*x*) is due to a possible initial focusing of the pulse due to the excitation of the circular antenna, see Figure 3. It is considered the symmetrical case with respect to *x* = 0.

Figure 3. Possible initial focusing of THz beams due to the circular antenna. The symmetric case is considered with respect to *x* = 0. In 3D axially symmetric case the transverse coordinate is *ρ*.

At the input of the nonlinear crystal *z* = 0 the EM wave gets the phase shift
$-k\Delta \left(x\right)=-k\left(R-{\left({z}_{c}^{2}+{x}^{2}\right)}^{1/2}\right)$ , where
${z}_{c}$ is the distance from the input to the focus point
${z}_{c}={\left({R}^{2}-{\left({L}_{x}/2\right)}^{2}\right)}^{1/2}$ , *k* is the wave number; thus the phase multiplier is
$M\left(x\right)=\mathrm{exp}\left(-ik\Delta \left(x\right)\right)$ . In 3D case *x* is replaced by the radial coordinate *ρ*.

4. Simulations of Generation of Harmonics

The simulations of generation of higher harmonics in the nonlinear crystal SrTiO_{3} have been realized under different temperatures, both in 2D and 3D geometries. It has been obtained that the pure cubic nonlinearity without the bias electric field is not effective for this generation. There is the optimum value of the bias electric field *E _{s}* ≈ (0.2 - 0.4)

The numerical simulations have been provided by the implicit finite difference schemes [46] . Several iterations have been applied that demonstrate good convergence. The simulations have demonstrated a possibility of generation of higher harmonics with the numbers *n* ≥ 5, *i.e.* the frequency multiplication, or frequency up-conversion, in the whole temperature interval *T* = 60 - 200 K.

The typical results of simulations are presented in Figures 4-8. It is seen that in the focused beams the maxima of different higher harmonics are realized at different distances from the input of the crystal.

At lower temperatures the nonlinearity is higher, so the input amplitudes of the first harmonics can quite lower, of about 1 kV/cm. But the highest frequencies that can be obtained under the up-conversion are of about 0.5 THz. In turn, at higher temperatures it is possible to excite the frequencies up to *f* º *ω*/2π = 1.5 THz.

From our simulations in is seen that the concentration of the energy near the focus is naturally higher in 3D beams. For this reason, in 2D beams the influence

Table 1. The used parameters of crystalline SrTiO_{3}.

(a) (b) (c) (d) (e) (f)

Figure 4. The generation of higher harmonics under the temperature *T* = 60 K in 2D beam. Part (a) is the longitudinal dependences of |*A _{j}*(

(a) (b) (c) (d) (e) (f)

Figure 5. The generation of higher harmonics under the temperature *T* = 60 K in 3D beam. Part (a) is the longitudinal dependences of |*A _{j}*(

(a)(b)

Figure 6. The generation of higher harmonics under the temperature *T* = 77 K. The longitudinal dependences of |*A _{j}*(

(a)(b)

Figure 7. The generation of higher harmonics under the temperature *T* = 150 K. The longitudinal dependences of |*A _{j}*(

(a)(b)

Figure 8. The generation of higher harmonics under the temperature *T* = 200 K. The longitudinal dependences of |*A _{j}*(

of the cubic nonlinearity and the zeroth harmonic is not essential. In 3D beams this influence should be taken into account generally, because near the maximum of the focusing the value of the amplitude of the first harmonic becomes comparable with the value of the bias electric field. The using of 3D beams makes possible to increase of the maxima of higher harmonics compared with 2D beams, but the relative values of these maxima are smaller in 3D beams than in 2D ones.

Because the value of the temperature of the crystal is the important parameter to realize the frequency multiplication, it is better to use the input pulses of durations < 1 μs to avoid the heating of crystals.

5. Conclusions

In the lower part of the terahertz range in the nonlinear crystalline paraelectrics like SrTiO_{3}, it is possible to observe the generation of higher harmonics and thus to realize the frequency up-conversion. The crystalline paralectric SrTiO_{3} possesses the cubic dielectric nonlinearity, so to increase the efficiency of harmonic generation, it is rather better to apply the bias electric field and to get the induced quadratic nonlinearity. In THz range, the frequency dispersion takes place, so the number of harmonics is large but finite, of about 10 - 20. The initial focusing of the beams by the circular antenna results not only in higher efficiency of generation, but also makes possible to select the maxima of different harmonics at specified distances from the input.

The generation of higher harmonics can be realized in the wide temperature range 60 - 200 K. At lower temperatures, the possible frequencies are smaller, of about 0.5 THz, but the nonlinearity is higher. At higher temperatures, the frequency up-conversion can be observed up till the frequencies of about 1.5 THz.

Both plane and cylindrical focused beams can be used for generation of harmonics. In cylindrical beams, the absolute maxima of harmonics can be higher than in plane ones, but the ratios of these maxima to ones of the first harmonic are lower in the cylindrical beams. The influence of the cubic nonlinearity and the zeroth harmonic is essential in the cylindrical beams due to high values of the focused first harmonic.

Acknowledgements

The authors are grateful to SEP-CONAHCyT (Mexico) for a partial support of our work.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

[1] |
Biswas, A., Banerjee, A., Acharyya, A., Inokawa, H. and Roy, J.N. (2020) Emerging Trends in Terahertz Solid-State Physics and Devices: Sources, Detectors, Advanced Materials, and Light-Matter Interactions. Springer, New York.
https://doi.org/10.1007/978-981-15-3235-1 |

[2] |
Carpintero, G., Garcıa Muñoz, L.E., Hartnagel, H.L., Preu, S. and Räisänen, A.V. (2015) Semiconductor Terahertz Technology Devices and Systems at Room Temperature Operation. Wiley, New York. https://doi.org/10.1002/9781118920411 |

[3] |
Bründermann, E., Hübers, H.-W. and Kimmitt, M.F. (2012) Terahertz Techniques. Springer, New York. https://doi.org/10.1007/978-3-642-02592-1 |

[4] | Song, H.-J. and Nagatsuma, T. (2015) Handbook of Terahertz Technologies. Devices and Applications. CRC Press, Taylor & Francis Group, Boca Raton. |

[5] |
Rieh, J.-S. (2021) Introduction to Terahertz Electronics. Springer, New York.
https://doi.org/10.1007/978-3-030-51842-4 |

[6] |
Choudhury, B., Rakesh, A.M. and Jha, M. (2016) Active Terahertz Metamaterial for Biomedical Applications. Springer, New York.
https://doi.org/10.1007/978-981-287-793-2 |

[7] |
Ganichev, S.D. and Prettl, W. (2006) Intense Terahertz Excitation of Semiconductors. Oxford University Press, Oxford.
https://doi.org/10.1093/acprof:oso/9780198528302.001.0001 |

[8] |
Kolejáka, P., Postava, K., Miáica, M., Kužel, P., Kadlec, F. and Pištora, J. (2018) Experimental Gouy Phase Shift Compensation in Terahertz Time-Domain Spectroscopy. Photonics and Nanostructures—Fundamentals and Applications, 31, 129-133.
https://doi.org/10.1016/j.photonics.2018.06.011 |

[9] | Dexheimer, S.L. (2008) Terahertz Spectroscopy Principles and Applications. CRC Press, Boca Raton. |

[10] |
O’Sullivan, C. and Murphy, J.A. (2012) Field Guide to Terahertz Sources, Detectors, and Optics. SPIE Press, Bellingham. https://doi.org/10.1117/3.952851 |

[11] |
Perenzoni, M. and Paul, D.J. (2014) Physics and Applications of Terahertz Radiation. Springer, New York. https://doi.org/10.1007/978-94-007-3837-9 |

[12] |
Siegel, P.H. (2002) Terahertz Technology. IEEE Transactions on Microwave Theory and Techniques, 50, 910-928. https://doi.org/10.1109/22.989974 |

[13] |
Mittleman, D. (2003) Sensing with Terahertz Radiation. Springer, New York.
https://doi.org/10.1007/978-3-540-45601-8 |

[14] |
Atakaramians, S., Shahraam Afshar, V., Monro, T.M. and Abbott, D. (2013) Terahertz Dielectric Waveguides. Advances in Optics and Photonics, 5, 169-215.
https://doi.org/10.1364/AOP.5.000169 |

[15] |
Woolard, D.L., Loerop, W.R. and Shur, M.S. (2003) Terahertz Sensing Technology, 2 Vols. World Scientific Publ., Singapore. https://doi.org/10.1142/5244 |

[16] | Lee, Y.-S. (2009) Principles of Terahertz Science and Technology. Springer, New York. |

[17] | Robertson, W.M. (1995) Optoelectronic Techniques for Microwave and Millimeter-Wave Engineering. Artech Publ., Norwood. |

[18] |
Gevorgian, S. (2009) Ferroelectrics in Microwave Devices, Circuits and Systems. Springer, New York. https://doi.org/10.1007/978-1-84882-507-9 |

[19] | Vendik, O.G. (1979) Ferroelectrics in Microwave Technology. Sov. Radio, Moscow. (In Russian) |

[20] |
Vendik, O.G. and Zubko, S.P. (1997) Phenomenological Description of the Permittivity of Strontium Titanate as a Function of Applied Electric Field and Temperature. Technical Physics, 42, 278-281. https://doi.org/10.1134/1.1258678 |

[21] | Blinc, R. and Žekš, B. (1974) Soft Modes in Ferroelectrics and Antiferroelectrics. North-Holland, Amsterdam. |

[22] | Lines, M.E. and Glass, A.M. (1977) Principles and Applications of Ferroelectric and Related Materials. Clarendon Press, Oxford. |

[23] | Rabe, K.M., Ahn, C.H. and Triscone, J.-M. (2007) Physics of Ferroelectrics. A Modern Perspective. Springer, New York. |

[24] |
Strukov, B.A. and Levanyuk, A.P. (1998) Ferroelectric Phenomena in Crystals. Springer, New York. https://doi.org/10.1007/978-3-642-60293-1 |

[25] | Maugin, G.A. (1988) Continuum Mechanics of Electromagnetic Solids. North-Holland, Amsterdam. |

[26] |
Buzin, I.M., Ivanov, I.V., Belokopytov, G.V., Sychev, V.M. and Chuprakov, V.F. (1981) Low-Temperature Ferroelectrics: Dielectric Permittivity, Losses, and Parametric Interactions at Ultrahigh Frequencies. Izvestiya VUZ, Fizika, 24, 6-28. (English Transl. Soviet Physics Journal) https://doi.org/10.1007/BF00941340 |

[27] | Rez, I.S. and Poplavko, Yu.M. (1989) Dielectrics. Basic Properties and Applications in Electronics. Radio and Svyaz’, Moscow. (In Russian) |

[28] | Poplavko, Yu.M., Pereverzeva, L.P., Voronov, S.O. and Yakimenko, Yu.I. (2007) Physical Material Science. Vol. 2. Dielectrics. KPI Publ., Kiev. (In Ukrainian) |

[29] | Gassanov, L.G., Koshevaya, S.V., Narytnik, T.N. and Omel’yanenko, M.Yu. (1978) Parametric and Nonlinear Interaction of Electromagnetic Waves in Paraelectrics. Izv VUZ Radioelektronika, 21, 56-63. (English Transl. Radioelectronics and Communications Systems) |

[30] | Gassanov, L.G., Koshevaya, S.V. and Omel’yanenko, M.Yu. (1980) On Frequency Multiplication in Paraelectrics. Radiotekhnika i Elektronika, 25, 1238-1243. (English Transl. Radio Engineering and Electronic Physics) |

[31] | Koshevaya, S.V., Kononov, M.V. and Omel’yanenko, M.Yu. (1985) An Influence of Dispersion in Waveguiding Systems with Cubically Nonlinear Dielectrics. Izv VUZ Radioelektronika, 28, 53-56. (English Transl. Radioelectronics and Communications Systems) |

[32] |
Koshevaya, S.V., Grimalsky, V.V., Kotsarenko, Y.N. and Tecpoyotl-T, M. (2016) Modulation Instability of Transversely Limited Electromagnetic Waves of Terahertz Range in Strontium Titanate Paraelectric. Radioelectronics and Communications Systems, 59, 489-495. https://doi.org/10.3103/S0735272716110029 |

[33] |
Grimalsky, V., Koshevaya, S., Escobedo-Alatorre, J. and Tecpoyotl-Torres, M. (2016) Nonlinear Terahertz Electromagnetic Waves in SrTiO_{3} Crystals under Focusing. Journal of Electromagnetic Analysis and Applications (JEMAA), 8, 226-239.
https://doi.org/10.4236/jemaa.2016.810021 |

[34] |
Grimalsky, V., Koshevaya, S., Escobedo-Alatorre, J. and Jatirian-Foltides, E. (2017) Formation of Short Terahertz Electromagnetic Pulses in Nonlinear Paraelectrics. Proceedings 30th International Conference on Microelectronics MIEL-2017, Nis, 9-11 October 2017, 91-94. https://doi.org/10.1109/MIEL.2017.8190076 |

[35] |
Grimalsky, V., Koshevaya, S., Escobedo-Alatorre, J. and Jatirian-Foltides, E. (2019) Stimulated Brillouin Scattering of Terahertz Electromagnetic Pulses in Paraelectrics. Applied Physics B, 125, 15-22. https://doi.org/10.1007/s00340-018-7125-4 |

[36] |
Grimalsky, V.V., Rapoport, Y.G., Boardman, A.D. and Koshevaya, S.V. (2018) Nonlinear Focusing of Picosecond Baseband Pulses in Paraelectric Crystals. Optical and Quantum Electronics, 50, 102-114. https://doi.org/10.1007/s11082-018-1369-4 |

[37] |
Grimalsky, V.V., Rapoport, Yu.G., Koshevaya, S.V., Escobedo-Alatorre, J. and Tecpoyotl-Torres, M. (2021) Nonlinear Focusing of Picosecond Baseband Pulses in Paraelectric Crystals in a Wide Temperature Range. Optical and Quantum Electronics, 53, Paper 484. https://doi.org/10.1007/s11082-021-03104-6 |

[38] |
Rapoport, Yu.G., Grimalsky, V.V., Koshevaya, S.V. and Melendez-Isidoro, D.L. (2015) Modulation Instability of Terahertz Electromagnetic Pulses in SrTiO_{3} Paraelectric. Proceedings 2015 IEEE 35th International Conference on Electronics and Nanotechnology (ELNANO), Kyiv, 21-23 April 2015, 131-134.
https://doi.org/10.1109/ELNANO.2015.7146851 |

[39] |
Zamudio-Lara, A., Koshevaya, S.V., Grimalsky, V.V. and Ya?ez-Cortes, F. (2015) Frequency Multiplication of Terahertz Radiation in the Crystals of Strontium Titanate Paraelectric. Radioelectronics and Communications Systems (Izvestiya VUZ, Radioelektronika), 58, 411-416. https://doi.org/10.3103/S0735272715090034 |

[40] |
Kozina, M., Fechner, M., Marsik, P., et al. (2019) Terahertz-Driven Phonon Upconversion in SrTiO_{3}. Nature Physics, 15, 387-392.
https://doi.org/10.1038/s41567-018-0408-1 |

[41] |
Kamarás, K., Barth, K.L., Keilmann, F., Henn, R., Reedyk, M., Thomsen, C., Cardona, M., Kircher, J., Richards, P.L. and Stehle, J.L. (1995) The Low Temperature Infrared Optical Functions of SrTiO_{3} Determined by Refectance Spectroscopy and Spectroscopic Ellipsometry. Journal of Applied Physics, 78, 1235-1240.
https://doi.org/10.1063/1.360364 |

[42] |
Yashchyshyn, Y., Godziszewski, K., Bajurko, P., Modelski, J., Szafran, M., Bobryk, E., Pawlikowska, E., Tarapata, G., Weremczuk, J. and Jachowicz, R. (2013) Tunable Ferroelectric Ceramic-Polymer Composites for Sub-THz Applications. Proceedings 43rd European Microwave Conference, Nuremberg, 7-10 October, 2013, 676-679.
https://ieeexplore.ieee.org/document/6686746 https://doi.org/10.1109/APMC.2013.6695095 |

[43] | Bloembergen, N. (1965) Nonlinear Optics. W. A. Benjamin Inc., New York. |

[44] |
Kivshar, Y.S. and Agrawal, G.P. (2003) Optical Solitons. From Fibers to Photonic Crystals. Academic Press, New York.
https://doi.org/10.1016/B978-012410590-4/50012-7 |

[45] |
Weiland, J. and Wilhelmsson, H. (1977) Coherent Non-Linear Interaction of Waves in Plasmas. Pergamon Press, London.
https://doi.org/10.1007/978-1-4757-1571-2_29 |

[46] |
Samarskii, A.A. (2001) The Theory of Difference Schemes. Marcel Dekker Inc., New York. https://doi.org/10.1201/9780203908518 |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2024 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.