CC BY
17
Mathematical modeling
Research article
UDC 519.6 PACS 07.05.Tp
DOI: 10.22363/2658-4670-2019-27-4-325-342
Leaky waves in planar dielectric waveguide
Dmitriy V. Divakov1, Alexandre A. Egorov2, Konstantin P. Lovetskiy1, Leonid A. Sevastianov1,3, Andrey S. Drevitskiy1
1 Department of Applied Probability and Informatics Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya St., Moscow 117198, Russian Federation 2 A. M. Prokhorov General Physics Institute Russian Academy of Sciences, Moscow 119991, Russian Federation 3 Bogoliubov Laboratory of Theoretical Physics Joint Institute for Nuclear Research 6, Joliot-Curie St., Dubna, Moscow region 141980, Russian Federation
(received: December 23, 2019; accepted: December 30, 2019)
A new analytical and numerical solution of the electrodynamic waveguide problem for leaky modes of a planar dielectric symmetric waveguide is proposed. The conditions of leaky modes, corresponding to the Gamow-Siegert model, were used as asymptotic boundary conditions. The resulting initial-boundary problem allows the separation of variables. The emerging problem of the eigen-modes of open three-layer waveguides is formulated as the Sturm-Liouville problem with the corresponding boundary and asymptotic conditions. In the case of guided and radiation modes, the Sturm-Liouville problem is self-adjoint and the corresponding eigenvalues are real quantities for dielectric media. The search for eigenvalues and eigenfunctions corresponding to the leaky modes involves a number of difficulties: the problem for leaky modes is not self-adjoint, so the eigenvalues are complex quantities. The problem of finding eigenvalues and eigenfunctions is associated with finding the complex roots of the nonlinear dispersion equation. To solve this problem, we used the method of minimizing the zero order. An analysis of the calculated distributions of the electric field strength of the first three leaky modes is given, showing the possibilities and advantages of our approach to the study of leaky modes.
Key words and phrases: integrated optics, waveguide, Sturm-Liouville problem, dispersion relation, leaky modes, computer simulation
© Divakov D. V., Egorov A. A., Lovetskiy K. P., Sevastianov L.A., Drevitskiy A. S., 2019
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
1. Introduction
In the books by Marcuse [1], [2], Adams [3], Snyder and Love [4], Tamir [5], and other authors the terms "leaky rays" and "leaky modes" appear when discussing the propagation of polarized light in fiber optical waveguides with the refractive index of the core smaller than that of the cladding, and in planar waveguides with plates of material optically denser than the waveguide layer itself. In this case, Marcuse writes that the outflow of light from such a waveguide is akin to tunneling through a potential barrier in quantum mechanics. The "leaky light", in contrast to the "emitted light", propagates for quite a long time along the axis of the optical fiber. Similarly, in a planar waveguide, the resulting electromagnetic radiation propagates for some time at some distance along the waveguide, in contrast to the emitted light.
At the initial stage of the study of "leaky" modes, T. Tamir et al. [6]-[12], A.W. Snyder et al. [13]-[18], as well as other research teams [19]-[25], investigated the dispersion equations of optical waveguides written in terms of transition matrices from the point of view of choosing one (two) roots of an analytical function. The studies were executed using the theory of residues in Cauchy integrals.
In the papers by V. Shevchenko [24], [25] the behavior of guided modes during the transition of their wave numbers beyond critical values was analyzed, their transformation into leaky modes was shown, and the choice of quadrants to which the wave numbers should shift when passing through critical values was justified. At frequencies below critical, the reflection from the waveguide walls ceases to be complete, so that the waveguide modes continue to propagate experiencing incomplete internal reflection, because of which some radiation from the waveguide occurs. Such (improper) waveguide waves with radiative damping are called leaky waves.
Open waveguides as radiating systems were first investigated by Hansen [19], who proposed an antenna structure implemented using leaky waves. However, there was no understanding of the physical mechanism of the resulting waves. After all, the leaky waveguide mode is characterized by a complex longitudinal wave number with constant attenuation due to radiation losses.
In this case, the longitudinal attenuation leads to an exponential increase in the wave amplitude in the transverse direction. This fact violates the usual radiation condition for guided modes, described by the solutions of self-adjoint problems for the Helmholtz equation. The behavior of the resulting waves that seems non-physical was clarified by Marcuvitz [23] and Oliner [6]-[9].
In the papers by Oliner et al. [6]-[12], a detailed study of the complex roots of the dispersion equation that do not correspond to the guided modes is given. The study begins with the assumption that exponentially damped waves correspond to such roots, the experimental observation of which was earlier reported in Refs. [20]-[22]. First, using ray technique, and then with the help of mode analysis, the authors analyzed the wave solutions corresponding to four different roots of the fourth-power dispersion equation. Two of these roots correspond to solutions that exponentially decrease in the direction of propagation and are located symmetrically with respect to the axes of coordinates and the origin of coordinates. The rest two roots are rejected. Many publications of that time have been devoted to the analysis of the relative position of the variety of roots [6]-[18].
In the first decades of the study of leaky waves, the method of steepest descent was used as the most common method for their numerical search. In this case, the trajectory of the fastest descent comes close to the leakage poles, they begin to make a significant contribution or even dominate in the general directional pattern of an open waveguide. The field distribution of the resulting leaky wave increases in the transverse direction. However, the field amplitude remains finite in a wedge-like region of space that allows leakage.
As shown by Marcuvitz [23], these complex poles can correspond to leaky modes. Although they do not make a direct contribution to the correct spectral solution and can therefore be characterized by non-physical growth towards infinity, they can nevertheless accurately describe the radiation field in limited spatial domains. In Ref. [26], e.g., it is noted that in most publications on leaky modes there are no plots of fields of various types for leaky modes calculated numerically (see, e.g., [4], [5], [27]-[35]). In this case, the authors of some publications (e.g., [34]) propose to replace the leaky modes with radiative ones in limited domains. Our studies have shown that this can lead, firstly, to a large error in the calculation of losses, and secondly, to an inaccurate calculation of the field profiles of leaky modes at distances exceeding several wavelengths (> 2) of the electromagnetic radiation used. The replacement of one wave with another sometimes used requires serious analysis in each specific case. As a consequence, there is an urgent need to develop new algorithms for calculating the fields of both radiative and leaky modes, surpassing the standard methods, e.g., the FDTD method, in count rate and not inferior to them in accuracy.
In quantum physics, such solutions of the stationary Schrodinger equation are called Gamow resonances [36], [37] or Siegert quasi-states [38]. In recent decades, some researchers (see [39], [40]) solve boundary-value problems for the Helmholtz equations with the asymptotic conditions of Siegert leaky waves, obtaining numerical results interpreted by them as leaky waves. We propose to obtain (using a numerical method) the solutions of boundary problems for wave equations with asymptotic conditions of Siegert leaky waves. The numerical solutions obtained using this approach coincide with the solutions of Refs. [27]-[29], [33], [35], [41], but additionally allow description of the phase fronts of leaky waves and "angular outflow cones".
In our opinion, a more rigorous justification of the model of leaky waves of open waveguide systems can be obtained by starting calculations not from the Helmholtz equation, as is traditionally done, but from the wave equation preceding the Helmholtz equation, and most importantly, more adequately reflecting the wave nature of leaky modes.
2. Statement of the problem of modeling leaky modes of symmetric waveguides
Consider (Figure 1) a symmetric three-layer planar waveguide consisting of a dielectric film having the height h with real refractive index nj, surrounded by a cladding layer with real refractive index nc <nj.
The propagation of radiation in such structures is described by the Maxwell equations, material equations [3]-[5], and boundary conditions that distinguish
the class of solutions interesting for the researcher — in the present case, the leaky modes [6]-[12], [26], [42].
Figure 1. Symmetric three-layer dielectric waveguide
The generally accepted model of the electromagnetic field in a planar (infinitely extended along the Oy-axis) are fields that are independent of the variable y. In this case, Maxwell's equations are considerably simplified, since dEa/dy = dHa/dy = 0 for any a = x,y,z, and they are divided into two independent subsystems — the subsystem for the so-called TE-modes and for the TM-modes. The subsystem for the TE-modes can be represented as a single wave equation for the master component Ey
da_
dx2 dz2
n2(x) d2 c2 dt2
Ev = 0,
(1)
with the boundary conditions
dEv
dx
d_Ev
dx
- ik0Pc (P) Ev
+ ik0Pc (P) Ev
and the initial conditions
0,
x=h+0
x=-0
(2)
0,
Ev (x,z,t)\ z=o =^0 W,
(3)
t=0
where pc (ft) = ^/n2 -ft2, and c is the electrodynamic constant, n (x) is the variable refractive index of the considered three-layer waveguide, defined below. The subsystem for the TE-mode also includes two equations for the connection of components Hx, Hz with the master component Ey.
As a model of leaky modes propagation, we will consider Eq. (1), i.e., in other words, we will consider the propagation of leaky modes in terms of a wave process. As asymptotic boundary conditions, we will consider the conditions of leaky modes corresponding to the Gamow-Siegert model [36]-
[40].
In the case under consideration, the function describing the refractive index depends only on x, which makes it possible to separate the variables in Eq. (1).
As a result, we obtain solutions corresponding to leaky modes propagating in the positive direction of the z-axis:
E+31 = A+ ■ exp lik0jn2, — Pj2x + ik0z — iutX , (4)
Eyj\ x<h = A% ' exP Yko Vn2 -ßi2x + ikoPjz — +
x>0 k )
+ B+fj ■ exp { —
{-ik0 V
n2 -ßj2x + ik0ßrjz — iut} , (5)
E+ \ ^ = B+ ■ exp {—ik0 Vrf2 — ß]x + ik0ßJz — iut} ,
(6)
where u is the frequency, and are the eigenvalues of the non-self-adjoint Sturm-Liouville problem with boundary conditions that extract the leaky modes [3], [34], [41]:
-X" + k2] n2 (x)X = k]ß2X, X' (0) + ik0 ^n2c — ß2X(0) = 0, ^X' (h) — ik0^n2c —ß2X(h) = 0,
(7)
The eigenfunctions of the problem (7) are defined as general solution of the ordinary differential equation subject to the boundary conditions of this problem, that is, they have the form
A ■ piko^n2-ß2{x-h)
J-Lc C- ,
x > h,
X(x) = <
Af ■e
Bc ■ e
ikoJn2f-ß2x + ^
-ik0^n2-ß2x
f ^ e x <0
ik0^/n2—ß2x
0 < x < h,
(8)
and the constants Ac j, Bcj are determined from the field joining conditions at the boundaries of the waveguide layer x = 0 and x = h, which with Eq. (8) taken into account constitute a homogeneous system of linear algebraic equations:
( 1
ik{) Pc
0 0
—eikoPfh
1
ikçj pf
—e-ik0pfh
ik0 pfe-lkoPfh 1
—ikçj pf
0 0
—1
ik() Pc )
(Ac\
Ai B,
\BC y
0 0 0
W
(9)
where pc = ^n2c — ß2, Pf = jJn2 — ß2.
The homogeneous system of Eqs. (9) has a nontrivial solution if and only if the determinant of the matrix of the system (9) is zero [43]. The equality to zero of the determinant of the matrix of the system (9) can be achieved
for some values of the spectral parameter ft, which, in turn, determine the eigenvalues of the problem (7).
In each subdomain x > h, 0 < x < h, and x < 0 the solution of the wave equation corresponding to the leaky modes is representable as a wave with a complex wave vector. In the case of a field corresponding to leaky modes running in the positive direction of the z-axis for x > h and x > h due to the symmetry of the waveguide, the wave vector is determined as
= (10)
and in the waveguide layer (0 < x < h) there are two waves with wave vectors
t% = ko H't) (11)
a
the modules of the wave vectors being equal to the corresponding wave
= k0nj.
finding solutions corresponding to the leaky modes as an eigenvalue problem for a differential operator with non-self-adjoint boundary conditions (7), which we will further solve numerically.
numbers: \k±| = k0nc and We formulate the problem o
3. Description of the algorithm for numerical solution of the leaky mode problem
The spectral problem for a differential operator with non-self-adjoint boundary conditions (7) is formulated numerically as a problem of approximate determination of complex solutions of the equation
det M(ft) = 0, (12)
where M(ft) denotes the coefficient matrix of Eqs. (9) [42]. In the case of a problem similar to (7), but with self-adjoint boundary conditions, any classical method of finding the real roots of the equation can be applied (see, e.g., [44]-[46]). Our problem (7) is not self-adjoint, therefore, the eigenvalues of this problem are generally complex and the standard methods for root search can no longer be applied.
The problem (12) can be reformulated as a problem of finding the minimum of a function of two variables as follows. The desired quantity ft = ft' + ift" is
a complex number. Any solution ft = ft' + ift" of Eq. (12) will be also a local minimum of the non-negative function
F(ft' ,ft") = |det M(ft' + ift" )|2. (13)
The eigenvalues corresponding to the leaky modes are localized in the first quadrant of the complex plane Re (ft) > 0, Im (ft) > 0; moreover 0 < Re (ft) < nc (see Refs. [34], [47], as well as our papers [42]).
To find all the local minima of function (11) in this region, it is proposed to introduce a mesh in the region 0 < Re (ft) < nc, 0 < Im (ft) < B, where B is a constant that defines the boundary of the search for eigenvalues [42]. The nodes of the introduced mesh are used as initial approximations for the numerical method of finding the minimum of the function of two variables (13). In Refs. [42] the Hook-Jeeves method [48] was used, but there are also other efficient numerical methods for zero-order multidimensional minimization [48].
4. Analysis of leaky modes in terms of inhomogeneous
plane waves
Consider the complex values
Pcj = ^n2 — ft2 = V'cj + VP"C3, Pfj = ^n2f -ftf = P'fj + ip"3,
and ftj = ft' + iftj in the solutions (4)-(6), explicitly distinguishing their real and imaginary parts, which will allow us to reformulate Eqs. (4)-(6) in terms of inhomogeneous waves, whose amplitude is also a function of coordinates x and z:
E+1 ^ = A+ ■ exp {—k0p'{,jx — k0 ftj z} ■ exp {ik0p'cjx + ik0ft'z — iut} , (14) E+1 x<h = A+f] ■ exp {—k(}p"f]x — ko ftz} ■ exp [ik(}p'fjx + ik(} ft'z — iut] +
x>0
+ B+- ■ exp {k0p"fjx — k0 ftj z} ■ exp {—ik0p'^jx + ik0 ft'jz — iut} , (15)
E+j | ^ = A+j ■ exp {k0p'{,jx — k0 ftj z} ■ exp {—ik0p'cjx + ik0ft'z — iut} . (16) Consider the expression (14) in the form of an inhomogeneous wave
E++1 ^ = Acj (x,z)■ exp {ik0p'cjx + ik0 ft'z — iut} , (17)
where ACj (x, z) is the amplitude of the inhomogeneous wave defined as
Acj (x, z) = A+J ■ exp {—k0p"x — k0 ft"z} . (18)
Consider in more detail the inhomogeneous wave in the form (17) with variable amplitude Acj (x,z). If such waveguide parameters exist for which the inhomogeneous wave (17) of some j-th mode propagates in the direction
Sj = (k0p;cj, k0ft')T, along which Acj (x, z) ^ 0, then such wave will decay
(of course, if such direction exists). If a direction Sj = (k0p'Cj, k0ft') exists, along which Acj (x, z) = Const, then the wave (17) will become homogeneous
and, correspondingly, if a direction Sj = (k0p'Cj, k0ft')T exists, along which
^ then the wave will have infinitely growing amplitude along such a direction.
Consider first the case Ac- (x,z) = Const. The representation (18) allows choosing a direction in the plane xOz, along which Acj (x,z) = Const. This direction is described by the equation —k0p"c-x — k0z = 0. The symmetric direction k0p"jx — k0fi" z = 0 will be responsible for the constant amplitude of the inhomogeneous wave, corresponding to Eq. (16). The domain above the line —k0p"jx — k0Pj z = 0 corresponds to the domain of the amplitude growth, as well as the domain below the line k0p"c-x — k0fi"z = 0 (see Figure 2).
Xi I -hp"«* -kjfr>0 -hp"Cjx ».x- koPJ - -kafc) -KP]z< 0
Kp"cJx .. k p"jx *** ** - -KP]z> 0 - kQP"z < 0 Z - - ^nP"x ,
Figure 2. Amplitude growth/attenuation regions and a constant amplitude line of an inhomogeneous plane wave on xOz plane
The region of existence of leaky modes corresponding to non-uniform waves with non-increasing amplitude is shown in Figure 2 (the cone between two dashed lines). Namely, if the wave vector of the non-uniform wave is located in the cone between two dashed lines, then this leaky mode has a non-increasing amplitude and can propagate in the positive direction of the z-axis.
Let us consider the inhomogeneous wave (14) in more detail. As shown in Appendix (see (22)), in the region k0p"jx + k0fi"z < 0 of a non-uniform wave the amplitude will increase indefinitely, therefore we will consider the wave (14) in the region k0p"- x + ko fij z > 0, see Figure 2. This exponential growth is real within a limited transverse distance surrounding the origin [49]. Using conservation of energy flux, one can show [49] that any mode that decreases exponentially as it propagates must increase exponentially transverse to the direction of propagation. However, it is evident that exponential growth of the field (and mode energy) that extends to infinity is unphysical since we have a finite energy source. A more detailed analysis of this problem is beyond the scope of our paper.
5. Numerical analysis of leaky modes of symmetric three-layer waveguides
Let us proceed to numerical analysis of the obtained representation of the leaky modes (14)-(16). Since the structure of the waveguide under
consideration is symmetric, it is sufficient to consider only Eq. (14). We give in Figure 3 the calculated values of the complex phase deceleration coefficient calculated for a waveguide with nc=1.47, nf =1.565, A = 0.55 ^m and h = 1.1A.
m
0.5 0.6 0.7 0.8 0.9 1
Ш(13)
1.1 1.2 1.3 1.4
Figure 3. Complex eigenvalues corresponding to the leaky modes of a symmetric waveguide
Let us analyze the distribution of the electric field strength of the first three leaky modes at a fixed point in time.
a) three-dimensional image
(b) two-dimensional projection
Figure 4. The real part Ey (x, z, t*) for the first leaky mode
As seen from Figure 4, the field is concentrated in a cone formed by lines of constant amplitude, and the maximum intensity is observed at the boundaries of the cone where the amplitude is maximal. Outside the cone, there is an area of infinite growth of the amplitude of inhomogeneous waves represented by Eqs. (14)—(16); in this area the inhomogeneous waves characterizing the leaky mode in cladding layers cannot exist.
Inside the guiding layer, on the contrary, the field represented by Eqs. (14)-(16) attenuates rather rapidly and becomes almost completely damped at a distance of several wavelengths. The fields in the coating layer and the substrate are inhomogeneous waves, whose amplitudes decay exponentially the stronger, the smaller the distance to the waveguide layer. Due
to the rapid attenuation of the field in the waveguide layer and the gradual removal (escape) of inhomogeneous waves, which characterize the behavior of the leaky mode in the cladding layers, one can observe a virtual "separation" of the leaky mode from the waveguide layer. A similar "separation" is also characteristic of higher leaky modes (see Figures 5,6).
1.5
1
0.5
0
0.5 W"
-1
a) three-dimensional image
0.5 1 1.5 2 2.5 z
(b) two-dimensional projection
Figure 5. The real part E, (x, z, t*) for the second leaky mode
(a) three-dimensional image (b) two-dimensional projection
Figure 6. The real part Ey (x,z,t*) for the third leaky mode
Further propagation of the emerging mode occurs outside the waveguide layer in the direction of the corresponding wave vectors in the upper and lower cladding layers. The leaky mode will propagate at an angle dj =
Re (p
arctg I -—-r- ) in the top layer and at an angle —9, in the substrate in
Re (Pj)
the form of two independent inhomogeneous waves.
We also note that experimental data on the propagation of leaky modes, which qualitatively agree with the results obtained by us, are given in [50] (see Figures 4(b), 5(b), 6(b)). In the experimental studies given in [50], leaky modes also propagate in a cone, and are also characterized by the presence of "separation" of the leaky mode from the waveguide layer.
It is commonly assumed that for the leaky modes the amplitude increases with the distance x from the waveguide along the vertical axis (at a fixed longitudinal distance zand in the absence of losses in the waveguide). However, as it propagates along the axis z, this mode decays due to permanent energy losses from the waveguide layer to the environment. Functionally, the fields of leaky modes (vertical profile) are identical to the fields of ordinary guided modes; however, since unlike normal guided (homogeneous) modes, the leaky modes are inhomogeneous waves. In this regard, the representation of leaky waves of planar waveguides using the solutions of the wave equation seems to be preferable for us, compared to the traditionally used representation using the solutions of the Helmholtz equation.
At the same time, some features were revealed that we plan to analyze in our further publications. It is important to emphasize that the region of existence of leaky modes corresponding to inhomogeneous waves with non-increasing amplitude is detected (the cone between two dashed lines in Figure 2). Moreover, if the wave vector of an inhomogeneous wave is located in the region of the cone between two dashed lines, then such a leaky mode has a non-increasing amplitude and can propagate for a sufficiently long distance in the waveguide without absorption.
6. Conclusion
As is well known, conventional guided modes that exist when the waveguide layer thickness is above the critical value are considered in the optical beam representation as plane waves propagating in a regular waveguide due to the total internal reflection of waves at the interfaces between the waveguide media. From this point of view, the leaky waves propagate due to the effect of disturbed total internal reflection: during each act of disturbed total reflection at the interfaces of the media forming the waveguide, some of the power of this guided mode is radiated, i.e. "flows out" into the space surrounding the waveguide.
It is important to emphasize that the number of leaky modes with a gradual leakage is limited, unlike the continuum of radiative modes. The resulting gradual leakage waves form a discrete spectrum and are plane inhomogeneous waves. On the contrary, the radiative modes form a continuum (their spectrum is continuous) and are plane homogeneous waves. As a result, the replacement of one kind of these waves with another kind requires serious analysis in each specific case. From this point of view, the methods developed by us are undoubtedly useful for theoretical and numerical studies of dielectric and, in particular, optical waveguides supporting leaky modes, for example, when
used as basic elements in the development of advanced sensors or various interface elements in integrated optical processors.
Acknowledgement
The publication funded by RFBR according to the research projects no. 1807-00567, no. 18-51-18005, and no. 19-01-00645.
Appendix
Consider a wave similar to a plane monochromatic wave, but having a com/ft' + ik" \
plex wave vector k = I 1 1 ), that can be represented as a non-uniform
\k z + k z J
plane wave:
U(x, z,t) = C- eik'x+ik'z-iut . e-k"xx-k"zz_ (19)
Lines of equal phase will be given by equations k'Xx + k'zz = Const. Lines of equal amplitude will be given by equations k"Xx + k"zz = Const. Lines of equal phase and equal amplitude will be orthogonal to each other in non-absorbing media due to the fact that k2x + k2z = k2, where k is the wave number corresponding to the medium in which the wave propagates.
In more detail: k'2x + k'2z — k"2x — k"2z + 2i k'xk"x + 2i k'zk"z = kz and equating the imaginary parts of the quantities in the right-hand and left-hand sides of this equality, we obtain a condition k'x k" x + k'zk" z = 0 that
represents the scalar product of vectors k' = Re (ft) and k" = Im (k), which is zero, therefore
k' ±k". (20)
Transition to phase-ray coordinates
We introduce new coordinates attached to the lines of constant phase and constant amplitude:
—— k xx + k gz I ^ —— k ^x ^r k ^z
By virtue of the previously established orthogonality property (20), the introduced coordinate system will be orthogonal. In the coordinates the form of the considered inhomogeneous plane wave (19) is considerably simplified:
U(£,ri,t) = C ■ e-r< .eii-iwt. (22)
In the new variables, according to (22), the amplitude of the inhomogeneous wave under consideration decreases along the positive direction of the axis Or],. For each fixed •q = q > 0 the non-uniform wave is characterized by the amplitude C ■ e-q and the behavior harmonic along £ (see Figure 7).
For each fixed £ = p, the amplitude of the non-uniform wave decreases in the positive direction of the axis O•q (see Figure 7).
>1
Figure 7. Profile of a leaky waveguide wave (mode) in phase-ray coordinates
In the half-plane •q < 0, the amplitude of the non-uniform wave will increase indefinitely, therefore from the physical point of view the half-plane •q < 0 corresponds the so-called shadow region of the non-uniform wave under consideration.
The proposed representation of a wave with a complex wave vector in the form of an inhomogeneous plane wave and the subsequent introduction of phase-ray coordinates clearly demonstrate the essential properties of such waves and will be further used to analyze leaky waves.
References
[1] D. Marcuse, Light Transmission Optics. New York: Van Nostrand Reinhold, 1972.
[2] D. Marcuse, Theory of Dielectric Optical Waveguides. New York: Academic, 1973.
[3] M. Adams, An Introduction to Optical Waveguides. Chichester: Wiley, 1981.
[4] A. Snyder and J. D. Love, Optical Waveguide Theory. London: Chapman and Hall, 1983.
[5] T. Tamir, Integrated Optics. Berlin: Springer-Verlag, 1982.
[6] L. O. Goldstone and A. A. Oliner, "Leaky-wave antennas i: Rectangular waveguides," IRE Transactions on Antennas and Propagation, vol. 7, no. 4, pp. 307-319, 1959, doi: 10.1109/TAP.1959.1144702.
[7] L. O. Goldstone and A. A. Oliner, "Leaky-wave antennas ii: Circular waveguides," IRE Transactions on Antennas and Propagation, vol. 9, no. 3, pp. 280-290, 1961, doi: 10.1109/TAP.1961.1144995.
[8] T. Tamir and A. A. Oliner, "Guided complex waves, part i: Fields at an interface," Proc. inst. Elec. Eng., vol. 110, no. 2, pp. 310-324, 1963, doi: 10.1049/piee.1963.0044.
[9] T. Tamir and A. A. Oliner, "Guided complex waves, part ii: Relation to radiation patterns," Proc. inst. Elec. Eng., vol. 110, no. 2, pp. 325-334, 1963, doi: 10.1049/piee.1963.0045.
[10] S.-T. Peng and A. Oliner, "Guidance and leakage properties of a class of open dielectric waveguides: Part i-mathematical formulations," IEEE Transactions on Microwave Theory and Techniques, vol. 29, no. 9, pp. 843855, 1981, doi: 10.1109/TMTT.1981.1130465.
[11] A. Oliner, S.-T. Peng, T.-I. Hsu, and A. Sanchez, "Guidance and leakage properties of a class of open dielectric waveguides: Part ii-new physical effects," IEEE Transactions on Microwave Theory and Techniques, vol. 29, no. 9, pp. 855-869, 1981, doi: 10.1109/TMTT.1981.1130466.
[12] F. Tamir and F. Y. Kou, "Varieties of leaky waves and their excitation along multilayered structures," IEEE Journal of Quantum Electronics, vol. 22, no. 4, pp. 544-551, 1986, doi: 10.1109/JQE.1986.1072991.
[13] R. Sammut and A. W. Snyder, "Leaky modes on a dielectric waveguide: Orthogonality and excitation," Applied optics, vol. 15, no. 4, pp. 10401044, 1976, doi: 10.1364/AO.15.001040.
[14] R. Sammut, C. Pask, and A. W. Snyder, "Excitation and power of the unbound modes within a circular dielectric waveguide," Proc. inst. Elec. Eng., vol. 122, no. 1, pp. 25-33, 1975, doi: 10.1049/piee.1975.0004.
[15] R. Sammut and A. Snyder, "Contribution of umbound modes to light absorption in visual photoreceptors," J. Opt. Soc. Am., vol. 64, no. 12, pp. 1171-1174, 1974, doi: 10.1364/JOSA.64.001711.
[16] A. W. Snyder, "Leaky-ray theory of optical waveguides of circular cross section," Appl. Phys., vol. 4, pp. 273-298, 1974, doi: 10.1007/BF00928381.
[17] A. W. Snyder and D. J. Mitchell, "Leaky rays on circular optical fibers," J. Opt. Soc. Am., vol. 69, no. 5, pp. 599-607, 1974, doi: 10.1364/JOSA.64.000599.
[18] A. W. Snyder and D. J. Mitchell, "Leaky mode analysis of circular optical waveguides," Optoelectronics, vol. 6, pp. 287-296, 1974, doi: 10.1007/BF01423378.
[19] W. W. Hansen, "Radiating electromagnetic waveguide," [U.S. Patent No. 2 402 622], 1940.
[20] S. Barone, "Leaky wave contributions to the field of a line source above a dielectric slab," Microwave Research Institute, Polytechnic Institute of Brooklyn, Report R-532-546, PIB-462, 1956.
[21] S. Barone and A. Hessel, "Leaky wave contributions to the field of a line source above a dielectric slab-part ii," Microwave Research Institute, Polytechnic Institute of Brooklyn, Report R-698-58, PIB-626, 1958.
[22] E. S. Cassedy and M. Cohn, "On the existence of leaky waves due to a line source above a grounded dielectric slab," IRE Trans. Microwave Theory Tech., vol. 9, no. 3, pp. 243-247, 1961, doi: 10.1109/TMTT.1961.1125314.
[23] N. Marcuvitz, "On field representations in terms of leaky modes or eigenmodes," IRE Trans. Antennas Propag., vol. 4, no. 3, pp. 192-194, 1956, doi: 10.1109/TAP.1956.1144410.
[24] V. V. Shevchenko, "On the behavior of wave numbers beyond the critical value for waves in dielectric waveguides (media with losses)," Radiophys. Quantum Electron., vol. 15, pp. 194-200, 1972, doi: 10.1007/BF02209117.
[25] V. V. Shevchenko, "The expansion of the fields of open waveguides in proper and improper modes," Radiophys. Quantum Electron., vol. 14, no. 8, p. 972, 1974, doi: 10.1007/BF01029499.
[26] A. A. Egorov, "Theoretical and numerical analysis of propagation and scattering of eigen- and non-eigenmodes of an irregular integrated-optical waveguide," Quantum Electronics, vol. 42, no. 4, pp. 337-344, 2012, doi: 10.1070/QE2012v042n04ABEH014809.
[27] K. Ogusu, M. Miyag, and S. Nishida, "Leaky te modes in an asymmetic three-layered slab waveguide," J. Opt. Soc. Am., vol. 70, no. 1, pp. 6872, 1980, doi: 10.1364/JOSA.70.000048.
[28] S. Yamaguchi, A. Shimojima, and T. Hosono, "Analysis of leaky modes supported by a slab waveguide," Electronics and Communications in Japan, vol. 73, no. 11, pp. 20-31, 1990, doi: 10.1002/ecjb.4420731103.
[29] E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, "Determination of guided and leaky modes in lossless and lossy planar multilayer optical waveguides: Reflection pole method and wavevector density method," J. Lightwave Technol., vol. 17, no. 5, pp. 929-941, 1999, doi: 10.1109/50.762914.
[30] J. Petracek and K. Singh, "Determination of leaky modes in planar multilayer waveguides," IEEE Photon. Technol. Lett., vol. 14, no. 6, pp. 810-812, 2002, doi: 10.1109/LPT.2002.1003101.
[31] A. G. Rzhanov and S. E. Grigas, "Numerical algorithm for waveguide and leaky modes determination in multilayer optical waveguides," Technical Physics, vol. 55, no. 11, pp. 1614-1618, 2010, doi: 10.1134/S1063784210110113.
[32] A. B. Manenkov, "Orthogonality conditions for leaky modes," Ra-diophysics and Quantum Electronics, vol. 48, pp. 348-360, 2005, doi: 10.1007/s11141-005-0076-8.
[33] A. A. Romanenko and A. B. Sotskii, "The solution of the dispersion relations for planar waveguides in the case of complex roots," Tech. Phys., vol. 43, no. 4, pp. 427-433, 1998, doi: 10.1134/1.1258999.
[34] E. I. Golant and K. M. Golant, "New method for calculating the spectra and radiation losses of leaky waves in multilayer optical waveguides," Technical Physics, vol. 51, no. 8, pp. 1060-1068, 2006, doi: 10.1134/S1063784206080160.
[35] A. B. Sotsky, L. Steingart, J. Jackson, and et al., "Prism excitation of leaky modes of thin films," Tech. Phys., vol. 58, no. 11, pp. 1651-1660, 2013, doi: 10.1134/S106378421311025X.
[36] G. Gamow, "Zur quantentheorie des atomkernes," Z. Phys., vol. 51, pp. 204-212, 1928, doi: 10.1007/BF01343196.
[37] O. Civitarese and M. Gadella, "Physical and mathematical aspects of gamow states," Phys. Rep., vol. 396, no. 2, pp. 41-113, 2004, doi: 10.1016/j.physrep.2004.03.001.
[38] A. J. F. Siegert, "On the derivation of the dispersion formula for nuclear reactions," Phys. Rev., vol. 56, pp. 750-752, 1939, doi: 10.1103/Phys-Rev.56.750.
[39] O. I. Tolstikhin, V. N. Ostrovsky, and H. Nakamura, "Siegert pseudo-states as a universal tool: Resonances, s matrix, green function," Phys. Rev. Lett., vol. 79, pp. 2026-2029, 1997, doi: 10.1103/Phys-RevLett.79.2026.
[40] O. I. Tolstikhin, V. N. Ostrovsky, and H. Nakamura, "Siegert pseudostate formulation of scattering theory: One-channel case," Phys. Rev. A, vol. 58, pp. 2077-2096, 1998, doi: 10.1103/PhysRevA.58.2077.
[41] D. Stowell and J. Tausch, "Guided and leaky modes of planar waveguides: Computation via high order finite elements and iterative methods," PIERS Online, vol. 6, no. 7, pp. 669-673, 2010, doi: 10.2529/PIERS091216124247.
[42] D. Divakov, A. Tiutiunnik, and A. Sevastianov, "Symbolic-numeric computation of the eigenvalues and eigenfunctions of the leaky modes in a regular homogeneous open waveguide," MATEC Web of Conferences, vol. 186, p. 01009, 2018, doi: 10.1051/matecconf/201818601009.
[43] R. E. Bellman, Introduction to matrix analysis. New York: McGraw-Hill, 1960.
[44] J. E. Dennis and R. B. Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Prentice Hall, 1983.
[45] R. Fletcher, Practical Methods of Optimization, 2nd. Wiley, 1987.
[46] J. Nocedal and S. Wright, Numerical Optimization, 2nd. Springer, 2006.
[47] J. Zhu and Y. Y. Lu, "Leaky modes of slab waveguides - asymptotic solutions," Journal of Lightwave Technology, vol. 24, no. 3, pp. 16191623, 2006, doi: 10.1109/JLT.2005.863275.
[48] C. T. Kelley, Iterative Methods for Optimization. SIAM, 1999.
[49] J. Hu and C. R. Menyuk, "Understanding leaky modes: Slab waveguide revisited," Advances in Optics and Photonics, vol. 1, no. 1, pp. 58-106, 2009, doi: 10.1364/AOP.1.000058.
[50] A. J. Martinez-Ros, J. L. Gómez-Tornero, F. J. Clemente-Fernandez, and J. Monzo-Cabrera, "Microwave near-field focusing properties of width-tapered microstrip leaky-wave antenna," IEEE Trans, on Antennas and Propagation, vol. 61, no. 6, pp. 2981-2990, 2013, doi: 10.1109/TAP.2013.2252138.
For citation:
D. V. Divakov, A. A. Egorov, K. P. Lovetskiy, L. A. Sevastianov, A. S. Drevitskiy, Leaky waves in planar dielectric waveguide, Discrete and Continuous Models and Applied Computational Science 27 (4) (2019) 325-342. DOI: 10.22363/2658-4670-2019-27-4-325-342.
Information about the authors:
Dmitriy V. Divakov — Candidate of Physical and Mathematical Sciences, assistant of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9550783, ORCID: https://orcid.org/0000-0001-5944-6187, ResearcherlD: N-3949-2013, Scopus Author ID: 55091105300) Alexandre A. Egorov — Doctor of Physical and Mathematical Sciences, Chief Researcher of Department of Oscillations of A. M. Prokhorov General Physics Institute (e-mail: [email protected], phone: +7(499)1350210, ORCID: https://orcid.org/0000-0002-1999-3810, ResearcherID: F-4681-2012, Scopus Author ID: 7202362326)
Konstantin P. Lovetskiy (Russian Federation) — Associate Professor, Ph.D., Associate Professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9522572, ORCID: https://orcid.org/0000-0002-3645-1060, ResearcherID: A-5725-2017, Scopus Author ID: 18634692900) Leonid A. Sevastianov (Russian Federation) — professor, Doctor of Physical and Mathematical Sciences, professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University); leading researcher of the Bogoliubov Laboratory of Theoretical Physics (e-mail: [email protected], phone: +7(495)9522572, ORCID: https://orcid.org/0000-0002-1856-4643, ResearcherID: B-8497-2016, Scopus Author ID: 8783969400) Andrey S. Drevitskiy — PhD student of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9550783, Scopus Author ID: 57198778786)
УДК 519.6 PACS 07.05.Tp
DOI: 10.22363/2658-4670-2019-27-4-325-342
Вытекающие моды в планарных диэлектрических
волноводах
Д. В. Диваков1, А. А. Егоров2, К. П. Ловецкий1, Л. А. Севастьянов1,3, А. С. Древицкий1
1 Кафедра прикладной информатики и теории вероятностей Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, 117198, Россия 2 Институт общей физики имени А.М. Прохорова Российской академии наук ул. Вавилова, д. 38, Москва, 119991, Россия 3 Лаборатория теоретической физики Объединённый институт ядерных исследований ул. Жолио-Кюри, д. 6, г. Дубна, Московская область, 141980, Россия
В работе предложено новое аналитическое и численное решение волновод-ной задачи для вытекающих мод планарного диэлектрического симметричного волновода. В качестве асимптотических граничных условий использовались граничные условия, соответствующие модели Гамова-Зигерта. Поставленная начально-краевая задача допускает разделение переменных. Возникающая в результате разделения переменных задача отыскания собственных мод открытых трёхслойных волноводов формулируется как задача Штурма-Лиувилля с соответствующими граничными и асимптотическими условиями. В случае направляемых и излучательных мод задача Штурма-Лиувилля является самосопряжённой, поэтому её собственные значения — действительные величины для диэлектрических сред. Поиск собственных значений и собственных функций, соответствующих вытекающим модам, сопряжён с рядом трудностей: задача на собственные значения и собственные функции не является самосопряжённой, поэтому собственные значения являются комплексными величинами, таким образом, задача нахождения собственных значений и собственных функций связана с нахождением комплексных корней нелинейного дисперсионного уравнения. В работе для решения этой задачи использовался метод минимизации нулевого порядка. В работе дан анализ рассчитанных распределений напряжённости электрического поля первых трёх вытекающих мод, показывающий возможности и преимущества предложенного подхода.
Ключевые слова: интегральная оптика, волновод, задача Штурма-Лиувилля, дисперсионное соотношение, вытекающие моды, компьютерное моделирование
