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

Jesus Escobedo-Alatorre, Volodymyr Grimalsky, Svetlana Koshevaya, Margarita Tecpoyotl-Torres

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

DOI: 10.4236/jemaa.2020.1211013                

Abstract

The superheterodyne amplification of electromagnetic waves is investigated when the resonant three-wave interaction of two electromagnetic waves with the space charge wave occurs in the waveguides nitride n-GaN, n-InN films-dielectric. The amplification of SCW waves due to the negative differential conductivity is investigated in nitride n-GaN, n-InN films at the frequencies f ≤ 400 GHz in the lower part of the terahertz (THz) range. The electromagnetic waves are either in the upper part of THz range or in the optical range. The superheterodyne amplification is considered in two geometries, the collinear one in which the three interacting waves travel in the same direction and the anti-collinear geometry where the second electromagnetic wave propagates in the opposite direction. The preferences and drawbacks of each geometry are pointed out. The finite width of space charge waves leads to decrease of increments of amplification.

Keywords

Superheterodyne Amplification, Negative Differential Conductivity, Space Charge Wave, Nitride Films, Terahertz Range

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 <2l, a dielectric substrate is at x < 0. Above at x > 2l 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 ω₁ and the longitudinal wave number k₁, backward EMW ω₂,k₂, and SCW ω₃,k₃. The resonant matching conditions are (see Figure 1):

$\begin{array}{l}{\omega }_3={\omega }_1-{\omega }_2,\\ k_3=k_1-k_2\text{for}\text{the}\text{collinear}\text{geometry};\\ k_3=k_1+\left|k_2\right|\text{for}\text{anti-collinear}\text{one}.\end{array}$ (1)

Figure 1. Geometry of the problem. The region 0<x <2l is a nitride film, x <0 is a dielectric substrate, x >2l is either vacuum or dielectric. The interacting waves are two EM waves (ω₁,k₁), (ω₂,k₂) and a space charge wave (ω₃,k₃). The waves are non-uniform along OY axis. The bias constant electric field is aligned along OZ axis, SCW propagates also along OZ. Part a) is the collinear geometry and b) is anti-collinear one.

The frequencies of EMW and SCW are of about: ${\omega }_{1,2}~{10}^{14}\text{-}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{-}4\times {10}^{12}{\text{s}}^{-1}$ ( $f_3={\omega }_3/2\pi ~150-700\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 ω_1,2. For the collinear geometry is rather better to choose the lowest (fundamental) EM mode at the frequency ω₁ and the third EM mode at ω₂.

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,v=\mu \left(\left|E\right|\right)E;\\ div\left({\epsilon }_0\epsilon \left(x\right)\nabla \tilde{\phi }\right)=-e\left(n-n_0\right),E=-\nabla \tilde{\phi }+E_0\end{array}$ (2)

Here n is the electron concentration, $\tilde{\phi }$ is the potential of the variable electric field, v is the electron drift velocity, n₀ 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₀ 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 ε₁,ε₂,ε₃.

For EMW the Maxwell equations are:

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

Here m^* is the effective electron mass, $\tilde{n}$ is the variable electron concentration of SCW; v is the high frequency electron velocity.

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,j_x\left(x=2l\right)=0;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 $\tilde{n}$ and for the variable part of the electric potential $\tilde{\phi }$ are searched as the travelling wave $\tilde{n}$, $\tilde{\phi }~\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 ω and a transverse wave number k_y are real (ω > 0), the longitudinal wave number is complex, $k_z={k^{\prime }}_z+i{k^{″}}_z$ ( ${k^{\prime }}_z>0$ ); within a certain frequency range there exists the imaginary part ${k^{″}}_z>0$. The case of NDC is under consideration: dv/dE<0.

The spatial increment ${k^{″}}_z$ depends essentially on the value of k_y. From the linearized equations it is possible to write down:

${k^{″}}_z\left(\omega ,k_y\right)\approx k^{″}\left(\omega ,0\right)-g{k}_y^2;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 2l and electron concentrations n₀ for the bias electric field $E_0=\text{15}0\text{kV}/\text{cm}$ are presented in Figure 2 for the case of transversely uniform SCW k_y= 0. The curve 1 is for $2l=0.5\mu \text{m}$,

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₂ above and below the film); the curve 2 is for $2l=0.5\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\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\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^{″}}_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)\exp \left(i\left({\omega }_{1}t-k_{1}z\right)\right)\\ +F_2\left(x\right)A_2\left(z,t,y\right)\exp \left(i\left({\omega }_{2}t-k_{2}z\right)\right))+c.c.;\\ \tilde{n}=-\frac{i{n}_0}{2}U\left(z,t,y\right)F_3\left(x\right)\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₃(x) are linear transverse profiles of the waves. The pump wave is A₂; A₁ 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{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{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². The temporal scale t_n is related to the motion of SCW, i.e. t_n= l_n/v₀, where l_n is the spatial scale, v₀ is the velocity of SCW. Here the dissipation coefficients for EM modes have been introduced Γ_1,2; Γ_U > 0 is the increment of spatial amplification of SCW that has been considered for the case k_y= 0. Also the saturation of amplification of SCW is taken into account when the maximum values of SCW are close to n₀, i.e. at |U|≤ 1.

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{collinear}\text{geometry}:A_2\left(z=0,t,y\right)=A_{20}\Psi \left(y\right)\\ \text{the}\text{anti-collinear}\text{geometry}:A_2\left(z=L_z,t,y\right)=A_{20}\Psi \left(y\right).\end{array}$ (9)

Here A₁₀ is a maximum amplitude of the small input pulse at the frequency ω₁, Φ(t) is the temporal shape of the pulse, Ψ(y) is its transverse shape; A₂₀ 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 L_y, it is quite big to avoid an influence of the boundaries along OY axis. The temporal shape of the input signal pulse A₁ is chosen as Gaussian-like, the transverse shape of EMW Ψ(y)is bell-like with the half-width of about y₀ ≥ 5 μm.

5. Superheterodyne Amplification of Electromagnetic Pulses

It is investigated the superheterodyne amplification of small input pulses of EMW at the carrier frequency ω₁, i.e. the transfer of amplification of SCW due to NDC to EMW at higher frequency ω₁ in the presence of the pump EMW at the frequency ω₂.

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 ω₃. 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 ω₁. The lengths of the waveguides are L_z= 40 - 100 μm.

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 ¶²U/¶y² 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₂ – n-GaN – SiO₂ 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₂₀= 0.04, it corresponds to the input intensity of about 20 kW/cm². The length of the system is L_z= 100 μm. The temporal scale is t_n= 4 ns for the waveguide on the base of n-GaN and t_n= 3 ns for n-InN. The duration of the input EM pulse is t₀= 40 ps for n-GaN and t₀= 40 ps for n-InN. The input width of EMW is y₀= 5 μm. At smaller widths an influence of wave diffraction and diffusion becomes important and reduces the efficiency of the amplification. It is seen that the superheterodyne amplification is characterized by the wide dynamic range, i.e. it is possible to amplify the input EM pulses with the intensities that differ 7 orders and more.

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₀≥ 10 ps. At smaller durations

Figure 3. The superheterodyne amplification of EM pulses under different input amplitudes A₁₀. The input pump amplitude is A₂₀= 0.04. Curve 1 is for A₁₀= 10⁻⁹, 2 is for A₁₀= 10⁻⁸, 3 is for A₁₀= 10⁻⁷, curve 4 is for A₁₀= 10⁻⁶. 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.

Conflicts of Interest

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

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