Different Geometries of Superheterodyne Amplification of Electromagnetic Beams in Waveguides Nitride-Dielectric

Show more

1. Introduction

The electromagnetic radiation of terahertz range (THz) is widely used in spectroscopy, medicine, scanning, and environmental science [1]. Now the sources of the THz radiation are both electron tubes (gyrotrons, backward wave tubes) and lasers. These sources are not well compatible with the traditional millimeter wave devices. An alternative way can be realized with the amplification of electromagnetic waves (EMW) under three-wave resonant interactions in semiconductor films where the negative differential conductivity (NDC) occurs. This method is based on the transfer of amplification of a space charge wave (SCW) at the microwave range or the lower part of THz range [2] to EMW at higher frequencies [3] [4]. EMW are either in the upper part of THz range or in the optical range, including infrared and visible ones.

Generally for the practical needs in the active media the transfer of the linear amplification of a low-frequency wave to high-frequency one is important [4] - [9], *i.e.* the superheterodyne amplification. A preference of the superheterodyne amplification is the using of relatively low levels of the electromagnetic pump wave at the intermediate frequency, which is utilized for the resonant wave coupling.

In this paper, different geometries of the superheterodyne amplification of EMW are investigated. The superheterodyne amplification of EM waves can realize in waveguides based on nitride *n-*GaN, *n-*InN films where the three-wave resonant interaction of two EMW with SCW amplified due to NDC can be realized. The superheterodyne amplification takes place both in the collinear geometry, where the propagating EM modes and the SCW propagate in the same direction, and in the anti-collinear one, where the pump EMW at the intermediate frequency propagates in the opposite direction. This amplification is due to the transfer of amplification from SCW to the signal EMW at the highest frequency and is effective, namely, the values of 20 - 40 dB can be reached at the lengths of 40 - 100 µm. An influence of the finite widths of EM waves in the plane of the film has been taken into account.

2. Model and Equations

It is considered *n-*GaN or *n-*InN film of a submicron thickness placed onto a dielectric substrate, see Figure 1. The nitride film is at 0 <*x *<2*l*, a dielectric substrate is at *x *< 0. Above at *x *> 2*l* there is either vacuum or a dielectric. Below the symmetric waveguide is considered, so above there is the same dielectric as at *x *< 0.

The nitride film is the waveguide for EM waves when the condition is satisfied: *ε*_{2EM}>*ε*_{1EM},*ε*_{3EM}. Here *ε*_{1,2,3EM} are corresponding permittivities in the optical range or upper part of THz one, which differ from their values in the microwave range due to the frequency dispersion [10] [11] [12] [13]. In such a dielectric waveguide it is possible to realize the resonant three-wave interaction of the following waves localized along *OX* axis: forward EMW at the frequency *ω*_{1} and the longitudinal wave number *k*_{1}, backward EMW *ω*_{2},*k _{2}*, and SCW

$\begin{array}{l}{\omega}_{3}={\omega}_{1}-{\omega}_{2},\\ {k}_{3}={k}_{1}-{k}_{2}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{for}\text{\hspace{0.17em}}\text{the}\text{\hspace{0.17em}}\text{collinear}\text{\hspace{0.17em}}\text{geometry};\\ {k}_{3}={k}_{1}+\left|{k}_{2}\right|\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{for}\text{\hspace{0.17em}}\text{anti-collinear}\text{\hspace{0.17em}}\text{one}.\end{array}$ (1)

Figure 1. Geometry of the problem. The region 0<*x *<2*l* is a nitride film, *x *<0 is a dielectric substrate, *x *>2*l* is either vacuum or dielectric. The interacting waves are two EM waves (*ω*_{1},*k*_{1}), (*ω*_{2},*k _{2}*) and a space charge wave (

The frequencies of EMW and SCW are of about:
${\omega}_{1,2}~{10}^{14}\text{\hspace{0.17em}}\text{-}\text{\hspace{0.17em}}4\times {10}^{15}{\text{s}}^{-1}$,
${\omega}_{3}\approx 2{\omega}_{1}\times \left({v}_{0}{\epsilon}_{2EM}^{1/2}/c\right)~{10}^{12}\text{\hspace{0.17em}}\text{-}\text{\hspace{0.17em}}4\times {10}^{12}{\text{s}}^{-1}$ (
${f}_{3}={\omega}_{3}/2\pi ~150\text{\hspace{0.17em}}-\text{\hspace{0.17em}}700\text{\hspace{0.17em}}\text{GHz}$ ) in the anti-collinear geometry. In the collinear geometry the SCW frequency is 3 - 5 times lower, it is
${\omega}_{3}={v}_{0}\cdot \left({k}_{1}-{k}_{2}\right)\equiv \left({\epsilon}_{{\omega}_{1}}^{1/2}-{\epsilon}_{{\omega}_{2}}^{1/2}\right)\cdot {\omega}_{1}\cdot {v}_{0}/c$. Here
${v}_{0}\approx 2.4\times {10}^{7}\text{cm}/\text{s}$ is the velocity of SCW in *n-*GaN [10] [11],
${\epsilon}_{{\omega}_{1,2}}$ * _{ }*are effective permittivities of EM waveguide modes at the frequencies

The dynamics of SCW is described by the equations of motion of the electron fluid jointly with the Poisson equation for the electric field. At the frequency range *f*≤ 400 GHz the simplest diffusion-drift equation can be applied:

$\begin{array}{l}\frac{\partial n}{\partial t}+div\left(v\left(E\right)n-D\nabla n\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}v=\mu \left(\left|E\right|\right)E;\\ div\left({\epsilon}_{0}\epsilon \left(x\right)\nabla \stackrel{\u02dc}{\phi}\right)=-e\left(n-{n}_{0}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}E=-\nabla \stackrel{\u02dc}{\phi}+{E}_{0}\end{array}$ (2)

Here *n* is the electron concentration,
$\stackrel{\u02dc}{\phi}$ is the potential of the variable electric field, *v* is the electron drift velocity, *n*_{0} is the equilibrium electron concentration, which is equal to the donor one; *D* is the diffusion coefficient, *μ*(*E*) is the electron mobility that depends on the value of the electric field, *E*_{0} is a bias constant electric field.

The data for nitrides GaN, InN are taken from [10] [11]. The coordinate frame is aligned along the crystalline axes. The lower indices 1, 2, 3 are related to the substrate, film, and the region over the film. The corresponding dielectric permittivities in the microwave range are *ε*_{1},*ε*_{2},*ε*_{3}.

For EMW the Maxwell equations are:

$\begin{array}{l}\nabla \times H={\epsilon}_{0}\epsilon \left(x\right)\frac{\partial E}{\partial t}+j,\text{\hspace{0.17em}}\text{\hspace{0.17em}}j=e\left({n}_{0}+\stackrel{\u02dc}{n}\right)v;\\ \nabla \times E=-\left(1/{\epsilon}_{0}{c}^{2}\right)\frac{\partial H}{\partial t},\\ B={\mu}_{0}H,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}^{*}\frac{\partial v}{\partial t}\approx eE.\end{array}$ (3)

Here *m ^{*}* is the effective electron mass,
$\stackrel{\u02dc}{n}$ is the variable electron concentration of SCW;

The influence on EMW on SCW is due to taking into account the Lorentz force in Equation (2), *i.e.* due to the substitution there

${E}_{z}\to {{E}^{\prime}}_{z}={E}_{z}-{\mu}_{0}{v}_{y}{H}_{x}.$ (4)

3. Linear Amplification of Space Charge Waves in Nitride Films

Here the linear amplification of SCW due to NDC is considered. At the surfaces of the film the following boundary conditions are used for the density of the electric current *j* [2]:

${j}_{x}\left(x=0\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{j}_{x}\left(x=2l\right)=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}j=e\left(v\left(E\right)n-D\nabla n\right).$ (5)

Equation (2) have been linearized, and the solutions for the perturbation of the electron concentration
$\stackrel{\u02dc}{n}$ and for the variable part of the electric potential
$\stackrel{\u02dc}{\phi}$ are searched as the travelling wave
$\stackrel{\u02dc}{n}$,
$\stackrel{\u02dc}{\phi}~\mathrm{exp}\left(i\omega t-i{k}_{z}z-i{k}_{y}y\right)$. Now an attention is paid to the transversely non-uniform case *k _{y}*¹ 0. The dispersion equation for SCW
${k}_{z}=k\left(\omega ,{k}_{y}\right)$ has been got from the substitution of the solutions within each partial region into the boundary conditions (5) and standard electric boundary conditions. In the case of spatial amplification of SCW, when a frequency

The spatial increment
${{k}^{\u2033}}_{z}$ depends essentially on the value of *k _{y}*. From the linearized equations it is possible to write down:

${{k}^{\u2033}}_{z}\left(\omega ,{k}_{y}\right)\approx {k}^{\u2033}\left(\omega ,0\right)-g{k}_{y}^{2};\text{\hspace{0.17em}}\text{\hspace{0.17em}}g\approx 2l\frac{{\omega}_{M}}{\omega}\frac{{\epsilon}_{2}}{{\epsilon}_{1}+{\epsilon}_{3}}$ (6)

Here
${\omega}_{M}=e{n}_{0}\mu /{\epsilon}_{0}{\epsilon}_{2}$ is the Maxwellian relaxation frequency,*
$\mu =\mu \left({E}_{0}\right)={v}_{0}/{E}_{0}$ *;
${v}_{0}=\mu \left({E}_{0}\right){E}_{0}$ is the constant part of the drift velocity. Therefore, a relatively small transverse non-uniformity of the beam of SCW ≤5 µm results in the essential decrease of the amplification, because the value of the relaxation frequency
${\omega}_{M}~{10}^{13}{\text{s}}^{-\text{1}}$ is high.

The results for *n-*GaN film with different thicknesses 2*l* and electron concentrations *n*_{0} for the bias electric field
${E}_{0}=\text{15}0\text{\hspace{0.17em}}\text{kV}/\text{cm}$ are presented in Figure 2 for the case of transversely uniform SCW *k _{y}*

Figure 2. Dependence of spatial increments of amplification of SCW on frequency in different *n-*GaN films.

${n}_{0}=2\times {10}^{17}{\text{cm}}^{-3}$,
${\epsilon}_{1}={\epsilon}_{3}=3.9$ (SiO_{2} above and below the film); the curve 2 is for
$2l=0.5\text{\hspace{0.17em}}\mu \text{m}$,
${n}_{0}=2\times {10}^{17}{\text{cm}}^{-3}$,
${\epsilon}_{1}={\epsilon}_{3}=8.5\left(\text{AlN}\right)$ ; the curve 3 is for
$2l=0.7\text{\hspace{0.17em}}\mu \text{m}$,
${n}_{0}=2\times {10}^{17}{\text{cm}}^{-3}$,
${\epsilon}_{1}={\epsilon}_{3}=8.5\left(\text{AlN}\right)$ ; the curve 4 is for
$2l=0.5\text{\hspace{0.17em}}\mu \text{m}$,
${n}_{0}=3\times {10}^{17}{\text{cm}}^{-3}$,
${\epsilon}_{1}={\epsilon}_{3}=8.5\left(\text{AlN}\right)$. Analogous results have been obtained for the waveguides on the base of *n-*InN films, where the increments of amplification are even greater.

Some parameters that are used below to investigate the superheterodyne amplification are: the frequency range is *f *≤ 400 GHz, the spatial increment is of about
${{k}^{\u2033}}_{z}\approx \text{2}\times {10}^{4}{\text{cm}}^{-\text{1}}$ ^{ }for the frequency *f* ≈ 200 GHz.

4. The Equations for Three-Wave Interaction

The nonlinear interaction is due to the modulation of the permittivity in the optical range (at higher frequencies) by SCW and due to the ponderomotive action of EMW to SCW in the microwave range (at lower frequencies). In the case of moderate nonlinearity it is possible to describe this resonant interaction by means of slowly varying wave amplitudes [3] [4] [5] [8] [12] [14].

It is investigated the interaction of transverse electric (TE) EM modes: *E = E _{y}*, the magnetic field is
$H=\left({H}_{x},0,{H}_{z}\right)$. For the transverse magnetic modes the results are similar.

At the boundaries of the film the standard EM boundary conditions are applied. For SCW the condition of the absence of the surface space charge leads to the conditions (5).

The solutions of Equations (2) and (3) for SCW and for EM waves are searched in the form:

$\begin{array}{l}E\left(x,y,z,t\right)=\frac{1}{2}({F}_{1}\left(x\right){A}_{1}\left(z,t,y\right)\mathrm{exp}\left(i\left({\omega}_{1}t-{k}_{1}z\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{F}_{2}\left(x\right){A}_{2}\left(z,t,y\right)\mathrm{exp}\left(i\left({\omega}_{2}t-{k}_{2}z\right)\right))+c.c.;\\ \stackrel{\u02dc}{n}=-\frac{i{n}_{0}}{2}U\left(z,t,y\right){F}_{3}\left(x\right)\mathrm{exp}\left(i\left({\omega}_{3}t-{k}_{3}z\right)\right)+c.c.\end{array}$ (7)

Here
${A}_{1,2}\left(z,t,y\right),U\left(z,t,y\right)$ are slowly varying amplitudes for EMW and SCW, *F*_{1,2}(*x*),*F*_{3}(*x*) are linear transverse profiles of the waves. The pump wave is *A*_{2}; *A*_{1} is EM pulse under amplification, *U* is SCW, which is produced under the three-wave interaction.

When using the orthogonality of waveguide modes [13], it is possible to get the coupled equations for the slowly varying amplitudes:

$\begin{array}{l}\frac{\partial U}{\partial t}+\frac{\partial U}{\partial z}-{\Gamma}_{U}U\cdot \left(1-{\left|U\right|}^{2}\right)-{g}_{U}\frac{{\partial}^{2}U}{\partial {y}^{2}}=i{A}_{1}{A}_{2}^{\ast};\\ \frac{\partial {A}_{1}}{\partial z}+i{g}_{1}\frac{{\partial}^{2}{A}_{1}}{\partial {y}^{2}}+{\Gamma}_{1}{A}_{1}=-{A}_{2}U;\end{array}$ (8a)

$\text{collinear}\text{\hspace{0.17em}}\text{geometry}:\frac{\partial {A}_{2}}{\partial z}+i{g}_{2}\frac{{\partial}^{2}{A}_{2}}{\partial {y}^{2}}+{\Gamma}_{2}{A}_{2}={A}_{1}{U}^{*};$ (8b)

$\text{anti-collinear}\text{\hspace{0.17em}}\text{geometry}:\frac{\partial {A}_{2}}{\partial z}-i{g}_{2}\frac{{\partial}^{2}{A}_{2}}{\partial {y}^{2}}-{\Gamma}_{2}{A}_{2}=-{A}_{1}{U}^{*}.$ (8c)

Equation (8) are presented in the non-dimensional form. The unity for EMW amplitudes corresponds to the intensity of about 10 MW/cm^{2}. The temporal scale *t _{n}* is related to the motion of SCW,

Because the nonlinear interaction of the waveguide modes is considered, the efficiency of the interaction is determined by the overlap integral *S* [13] of the interacting EM modes and SCW. The overlap integral is *S* ≈ 1 for the anti-collinear interaction and *S* ≈ 0.4 for collinear one, because under the anti-collinear interaction it is possible to choose the same fundamental modes for EMW at the frequencies *ω*_{1,2},* _{ }*whereas the interacting modes should be different for the collinear interaction.

One can see that the transverse nonuniformity along *OY* axis of EMW leads to the wave diffraction [9] [14], whereas the nonuniformity of SCW appears like the diffusion. The mechanism of the wave coupling is as follows. SCW modulates the effective permittivity of EM waves; the influence of EM waves on SCW is ponderomotive, due to the term Equation (4).

The boundary conditions for the slowly varying amplitudes are:

$\begin{array}{l}U\left(z=0,t,y\right)=0;\\ {A}_{1}\left(z=0,t,y\right)={A}_{10}\Phi \left(t\right)\Psi \left(y\right);\\ \text{the}\text{\hspace{0.17em}}\text{collinear}\text{\hspace{0.17em}}\text{geometry}:{A}_{2}\left(z=0,t,y\right)={A}_{20}\Psi \left(y\right)\\ \text{the}\text{\hspace{0.17em}}\text{anti-collinear}\text{\hspace{0.17em}}\text{geometry}:{A}_{2}\left(z={L}_{z},t,y\right)={A}_{20}\Psi \left(y\right).\end{array}$ (9)

Here *A*_{10} is a maximum amplitude of the small input pulse at the frequency *ω*_{1}, Φ(*t*) is the temporal shape of the pulse, Ψ(*y*) is its transverse shape; *A*_{20} is the constant amplitude of the EM pump wave;
${A}_{10}\ll {A}_{20}$. The length of the system is *L _{z}*. The width of the film is

5. Superheterodyne Amplification of Electromagnetic Pulses

It is investigated the superheterodyne amplification of small input pulses of EMW at the carrier frequency *ω*_{1}, *i.e.* the transfer of amplification of SCW due to NDC to EMW at higher frequency *ω*_{1} in the presence of the pump EMW at the frequency *ω*_{2}.

The mechanism of the superheterodyne amplification is as follows [3] [4]. Because of the three-wave interaction, the mixing of two EMW results in the generation of SCW at lower frequency *ω*_{3}. Then SCW is amplified in a medium with NDC. In that turn, in the output of the system the amplified EM wave appears at the frequency *ω*_{1}. The lengths of the waveguides are *L _{z}*

The difference schemes are used to approximate Equationi (8). The splitting with respect to physical factors has been used to take into account the wave diffusion term with ¶^{2}*U/*¶*y*^{2} in the equation for *U* [15] [16]. The upwind difference approximations have been applied for EM amplitudes *A*_{1,2}; for the diffraction terms the implicit difference approximation has been used.

The typical results of simulations are presented in Figure 3. The part a) is for the anti-collinear geometry, the part b) is for collinear one. The waveguide SiO_{2} – *n-*GaN – SiO_{2} is considered, see Figure 2, curve 1. In the anti-collinear geometry the signal EM frequency is
${\omega}_{1}\approx 3\times {10}^{14}{\text{s}}^{-1}$, SCW frequency is
${\omega}_{3}\approx 1.5\times {10}^{12}{\text{s}}^{-1}$. In the collinear geometry they are
${\omega}_{1}\approx 2\times {10}^{15}{\text{s}}^{-1}$, SCW frequency is
${\omega}_{3}\approx 1.5\times {10}^{12}{\text{s}}^{-1}$. The input amplitude for the pump is *A*_{20}*=* 0.04, it corresponds to the input intensity of about 20 kW/cm^{2}. The length of the system is *L _{z}*

Both geometries, anti-collinear and collinear, are suitable for the superheterodyne amplification. But the anti-collinear geometry is preferable for the efficiency, because the overlap integral of the interacting modes is bigger. But the collinear geometry seems preferable for decreasing the corresponding resonant frequency of SCW, *i.e.* for the easier realization of amplification due to NDC. In another words, in the collinear geometry it is possible to amplify EM signals at higher frequencies.

It is possible to amplify EM pulses of durations *t*_{0}*≥* 10 ps. At smaller durations

Figure 3. The superheterodyne amplification of EM pulses under different input amplitudes *A*_{10}. The input pump amplitude is *A*_{20}*=* 0.04. Curve 1 is for *A*_{10}*=* 10^{−9}, 2 is for *A*_{10}*=* 10^{−8}, 3 is for *A*_{10}*=* 10^{−7}, curve 4 is for *A*_{10}*=* 10^{−6}. The wave diffraction and diffusion are taken into account. The dot lines are the input pulses, the solid lines are the output ones. Part (a) is the anti-collinear geometry and (b) is collinear one.

the broadening of the pulses under the nonlinear interaction occurs.

When the amplitudes of input EM pulses are quite big, as ones at the curves 4 in Figure 3, the amplitudes of SCW under amplification become comparable with the equilibrium concentration, and the saturation of amplification takes place. This results in some distortion of the output EM pulse, especially near its maximum.

6. Conclusions

The superheterodyne amplification in the optical range and in the upper part of the terahertz range can be realized in waveguides on the base of *n-*GaN, *n-*InN films in two different geometries of the resonant three-wave interaction between two electromagnetic waves with the space charge wave. The superheterodyne amplification is due to the transfer of the amplification of the space charge wave in the medium with the negative differential conductivity to the signal electromagnetic wave. It can be realized in the waveguides of 40 - 100 μm lengths. The input durations of electromagnetic pulses can be 10 ps and more.

In the first geometry the propagation of the signal and pump electromagnetic waves is opposite, whereas in the second geometry all interacting waves propagate in the same direction. In the first geometry, the efficiency of the three-wave resonant interaction is higher than the second geometry, but possible frequencies of input optical signals can be higher in the second geometry.

The finite transverse widths of the signal and pump electromagnetic waves in the plane of the film have been taken into account. They result in the diffraction of electromagnetic modes and in the analog of the diffusion of the space charge wave. To avoid the influence of the wave diffraction and diffusion, the transverse sizes of the electromagnetic modes should be no smaller than 5 μm.

The dynamic range of the superheterodyne amplification without the saturation is of about 7 orders and more. The nonlinear saturation results in some distortion of the amplified pulses.

Acknowledgements

The authors thank SEP-CONACyT, Mexico, for partial support of our work.

References

[1] Siegel, P.H. (2002) Terahertz Technology. IEEE Transactions on Microwave Theory and Techniques, 50, 910-928.

https://doi.org/10.1109/22.989974

[2] Barybin, A.A. (1986) Waves in Thin Film Semiconductor Structures with Hot Electrons. Nauka, Moscow.

[3] Koshevaya, S., Grimalsky, V., Escobedo-Alatorre, J. and Tecpoyotl-Torres, M. (2003) Superheterodyne Amplification of Sub-Millimeter Electromagnetic Waves in an n-GaAs Film. Microelectronics Journal, 34, 231-235.

https://doi.org/10.1016/S0026-2692(03)00003-X

[4] Kalapusha, A.L. and Kotsarenko, N.Y. (1983) Acoustoelectronic Parametric Amplification of EM Waves of Infrared and Optical Ranges in Planar Optical Waveguides. Radioelectronics and Communications Systems (Izvestiya VUZ, Radioelektronika), 26, 71-72.

[5] Weiland, J. and Wilhelmsson, H. (1977) Coherent Non-Linear Interaction of Waves in Plasmas. Pergamon Press, London.

[6] Stenflo, L. (1996) Theory of Nonlinear Plasma Surface Waves. Physica Scripta, T63, 59-62.

https://doi.org/10.1088/0031-8949/1996/T63/008

[7] Rabinovich, M.I. and Fainshtein, S.M. (1973) High-Frequency Instability of an Electromagnetic Wave in a Nonequilibrium Magnetized Plasma, Journal of Experimental and Theoretical Physics, 36, 884-886.

http://jetp.ac.ru/cgi-bin/dn/e_036_05_0884.pdf

[8] Rabinovich, M.I. and Trubetskov, D.I. (1989) Oscillations and Waves: In Linear and Nonlinear Systems (Mathematics and Its Applications). Springer, Berlin.

[9] Sukhorukov, A.P. (1988) Nonlinear Wave Interactions in Optics and Radiophysics (Russian Edition). Nauka, Moscow.

[10] Levinshtein, M.E., Rumyantsev, S.L and Shur, M.S. (2001) Properties of Advanced Semiconductor Materials: GaN, AlN, InN. Wiley-Interscience, New York.
www.ioffe.ru/SVA/NSM/Semicond/GaN/

[11] Jain, S.C., Willander, M., Narayan, J. and Van Overstraeten, R. (2000) III-Nitrides: Growth, Characterization, and Properties. Journal of Applied Physics, 87, 965-1006.

https://doi.org/10.1063/1.371971

[12] Chaika, G.E., Malnev, V.N. and Panfilov, M.I. (1996) Interaction of Light with Space Charge Waves. Proceedings of SPIE, 2795, 279-285.

https://doi.org/10.1117/12.239220

[13] Marcuse, D. (1991) Theory of Dielectric Optical Waveguides. Academic Press, New York.

[14] Kalinovich, A., Lobanov, V.E., Sukhorukov, A.P., and Zverev, D.M. (2013) Interaction of Pulsed Laser Beams in Quadratic Nonlinear Media. Physics of Wave Phenomena, 21, 5-9.

https://doi.org/10.3103/S1541308X13010020

[15] Marchuk, G.I. (1990) Splitting and Alternating Direction Methods. In: Ciarlet, P.G. and Lions, J.L., Eds., Handbook of Numerical Analysis, Elsevier, Amsterdam, 203-462.

https://doi.org/10.1016/S1570-8659(05)80035-3

[16] Samarskii, A.A. (2001) The Theory of Difference Schemes. CRC Press, Boca Raton.

https://doi.org/10.1201/9780203908518