CC BY
5
Research article
UDC 519.61; 519.62; 621.372.82
DOI: 10.22363/2658-4670-2020-28-3-252-273
Asymptotic method for constructing a model of adiabatic guided modes of smoothly irregular integrated optical waveguides
Anton L. Sevastianov
Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya St., Moscow, 117198, Russian Federation
(received: June 10, 2020; accepted: September 14, 2020)
The paper considers a class of smoothly irregular integrated optical multilayer waveguides, whose properties determine the characteristic features of guided propagation of monochromatic polarized light. An asymptotic approach to the description of such electromagnetic radiation is proposed, in which the solutions of Maxwell's equations are expressed in terms of the solutions of a system of four ordinary differential equations and two algebraic equations for six components of the electromagnetic field in the zero approximation. The gradient of the phase front of the adiabatic guided mode satisfies the eikonal equation with respect to the effective refractive index of the waveguide for the given mode.
The multilayer structure of waveguides allows one more stage of reducing the model to a homogeneous system of linear algebraic equations, the nontrivial solvability condition of which specifies the relationship between the gradient of the radiation phase front and the gradients of interfaces between thin homogeneous layers.
In the final part of the work, eigenvalue and eigenvector problems (differential and algebraic), describing adiabatic guided modes are formulated. The formulation of the problem of describing the single-mode propagation of adiabatic guided modes is also given, emphasizing the adiabatic nature of the described approximate solution of Maxwell's equations.
Key words and phrases: smoothly irregular integrated optical multilayer waveguides, eigenvalue and eigenvector problem, single-mode propagation of adiabatic guided modes
1. Introduction
Fundamental results in the theory of regular waveguides were obtained for closed (metallic) waveguides by A. N. Tikhonov and A. A. Samarskii [1], and for open (dielectric) waveguides by A. G. Sveshnikov [2] and V. V. Shevchenko [3]. Among the irregular waveguides, one can distinguish transversely irregular and longitudinally irregular waveguides. For transversely irregular waveguides, the equations and the corresponding solutions allow the separation of variables [4].
© Sevastianov A.L., 2020
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
Here the incomplete Galerkin method developed by A. G. Sveshnikov [2], [5], [6] received the greatest recognition.
For closed longitudinally irregular waveguides, B. Z. Katsenelenbaum developed the method of cross sections [7], which was generalized for open longitudinally irregular waveguides by V. V. Shevchenko [8]. These models do not describe depolarization and hybridization of guided modes in irregular sections of waveguides. A. A. Egorov, L. A. Sevastyanov and A. L. Sevastyanov developed the foundations of the theory of smoothly irregular 3D dielectric and, in particular, integrated optical waveguides [9], [10], which was successfully applied to a number of three-dimensional integrated optical waveguides and smoothly irregular 3D waveguide devices based on them [11]—[13]. The mathematical basis of the model of adiabatic guided modes (AGMs) is the asymptotic method and the method of coupled modes. The asymptotic method for solving a boundary value problem for a system of differential equations with respect to a small parameter 5 allows it to be reduced to a system of ordinary differential equations with special boundary conditions, the method of solving which is known. The coupling of two second-order equations for modes of two different polarizations when solving the original system of equations by the asymptotic method, manifests itself in the first approximation as a weak (of the order of 5) coupling of two linear oscillators. It reflects the violation of the structure regularity caused by a change in the phase constant of smoothly irregular dielectric waveguides.
In this paper, we consider an approach to the construction of a model of propagation of electromagnetic radiation in integrated optical smoothly irregular waveguide structures. Traditionally, such models are described using Maxwell's equations. The paper considers only monochromatic radiation, depending on time as exp(iut). Such time dependence of the solution allows considering a model of steady-state guided propagation of electromagnetic radiation.
2. Basic concepts and notations
Guided propagation of monochromatic polarized electromagnetic radiation in integrated optical waveguides is described by Maxwell's equations. The electromagnetic field is described using complex amplitudes. A material medium is considered, consisting of dielectric subdomains that fill the entire three-dimensional space. The permittivities of the subdomains are different and real, and the permeability is everywhere equal to the permeability of vacuum. It follows that in the absence of foreign currents and charges, the induced currents and charges are equal to zero.
In the absence of foreign charges and currents, the scalar Maxwell's equations follow from the vector ones, and the boundary conditions for the normal components follow from the boundary conditions for the tangential components [14]. The constitutive equations are assumed to be linear. Thus, the electromagnetic field in a space filled with dielectrics in the Gaussian system of units is described by equations
1 SB 1 SD
rotE = ——, rotH = ——, D = eE, B = ^H, (1)
where E, H are the electric and magnetic field strength vectors, D is the electric displacement vector, B is the magnetic flux density vector, c is the velocity of electromagnetic waves in vacuum.
At the interface between dielectric media 1 and 2, the tangential components of electric and magnetic field strengths satisfy the following boundary (matching) conditions:
H |1 = H |2, Er |1 = Et|2 . (2)
The asymptotic boundary conditions for guided modes at infinity
||E| -> 0, ||H|| -> 0, (3)
ensure uniqueness of the solution of the problem (1)-(3).
In equations (1), £ is the medium permittivity, ^ is the medium permeability. Let us denote by n = -^/Jle the refractive index of the medium (hereinafter — of a dielectric layer of the considered multilayer dielectric structure).
3. The considered class of objects
The object of our consideration is the guided propagation of monochromatic electromagnetic radiation of the optical range in thin-film integrated optical structures. Such structures are complex waveguide structures formed by the deposition of additional waveguide layers of various (smoothly irregular) geometric configurations on the base waveguide. As a base waveguide, we consider a regular planar three-layer waveguide filling the entire three-dimensional space (open waveguide): a substrate layer (substrate) is located in the lower half-space, then a guiding layer of constant thickness is located, and a cladding layer is located in the upper half-space. In this case, the guiding layer (core) is optically denser compared to the substrate and the cladding. By a thin-film waveguide, we mean a waveguide whose core thickness is comparable to the wavelength of propagating radiation.
Integrated optical waveguide structures are formed by introducing into the three-layer planar dielectric waveguide additional layers of variable thickness h(y, z). The additional waveguide layers are specified by the interface between the additional and waveguide layers x = h(y, z) and directly affect phase distribution tp(y, z) and effective refractive index neff (y, z) of the composite waveguide system. The phase velocity is smaller and the effective refractive index neff(y, z) is greater in the locations where the total thickness of the waveguide layers is greater.
In particular, it can be a few-mode integrated optical waveguide implemented in the form of a three-layer dielectric "plate". On this structure, a local but smoothly irregular optical inhomogeneity is deposited, i.e., an additional dielectric layer of variable thickness. A thin-film waveguide generalized lens (TFWGL) (see Figure 1) is an example [11], [12].
Figure 1. Schematic of Luneberg TFWGL with an additional waveguide layer having
cylindrical symmetry
In such integrated optical waveguide, a normal mode of the planar waveguide (allowed by the structure1, see Appendix) travelling from infinity (or from a point source in the plane of the planar waveguide) is deformed when meeting the localized irregularity (optical inhomogeneity). However, it preserves the structure of a 'transverse resonance', i.e., the structure of a standing wave in the 'transverse' (vertical) direction. Gorelyshev, Neishtadt et al. [20], [21] formulate this conditions as conservation of an adiabatic invariant.
After passing the irregularity region, the deformed "adiabatically invariant" guided (quasi)mode restores the initial (vertical) parameters of a normal mode of a regular waveguide (or transforms into a superposition of normal modes).
For convenience, let us define the Cartesian coordinate system so that all spatial subdomains corresponding to infinite dielectric layers be bounded by planes parallel to the yOz-plane and surfaces asymptotically parallel to the yOz-plane, so that hereafter e = e(x), ^ = 1.
We will call smoothly irregular the optical structures satisfying the inequalities specified by the geometry of the additional waveguide layer:
dh dh
dy dz
In weakly inhomogeneous 3D media the electromagnetic radiation propagation is described by locally planar waves or adiabatic approximations to the solutions of the Maxwell's equations, obtained using the asymptotic method [22]. By analogy with locally plane and locally spherical 3D waves [22],
1 From the theory of planar regular waveguides [15]—[19] it is known that an electromagnetic wave propagates through a regular waveguide in the form of a normal guided mode.
[23], we seek the guided propagation of electromagnetic monochromatic polarized radiation in a smoothly irregular integrated optical waveguide in the form of modified locally normal guided modes of reference waveguides1 .
The adiabatic approximation of the solution of Maxwell's equation obtained in this way will preserve the adiabatic invariants that reflect the guided character of light propagation (the so called transverse resonance condition [15], [26].
This work is devoted to the search for a model of adiabatic guided propagation of monochromatic electromagnetic radiation in smoothly irregular integrated optical structures. The technique for finding it is based on an asymptotic approach in the form of locally normal guided modes of a locally planar reference waveguide. The adiabatic approximation of the solution of Maxwell's equations obtained in this way will preserve adiabatic invariants reflecting the waveguide nature (the so-called transverse resonance condition) of light propagation.
4. Basic equations of the adiabatic guided mode model
Let us recall the earlier assumptions made for the considered integrated optical waveguides and the electromagnetic radiation propagating along them.
1. Electromagnetic radiation is optical and monochromatic with a fixed wavelength A e [380; 780], nm.
2. The thickness of the guiding layer (core) of the base thin-film waveguide is comparable to the length of the propagating monochromatic electromagnetic radiation d ~ A.
3. The surface of additional waveguide layer (x = h(y, z)) satisfies the
limiting conditions
dh dh
dy dz
« 1.
4. The integrated optical waveguide is a material medium consisting of dielectric subdomains that fill the entire space.
5. The permittivities of the subdomains are different and real, and the permeability is everywhere equal to the magnetic permeability of vacuum.
6. There are no external currents and charges. It follows from this that in the absence of external currents and charges, the induced currents and charges are equal to zero.
7. A Cartesian coordinate system is introduced as follows: the interfaces between the dielectric media of the basic three-layer waveguide are parallel to the yOz-plane. In this case, the subdomains of space corresponding to the cladding and substrate layers are semi-infinite, the additional waveguide layers are asymptotically parallel to the yOz-plane, so that £ = e(x).
In Cartesian coordinates associated with the geometry of the substrate (or a three-layer planar dielectric waveguide), Maxwell's equations are written in the form
xThe notion of reference waveguides (dielectric planar) is presented in papers by Kat-senelenbaum and Shevchenko [7], [8], [24], [25].
dHz dHy dEz dEy I1
dy dz c dt , dy dz c dt
dHx dHz edEy dEx dEz VdHy
dz dx c dt , dz dx c dt
dHy dHx edEz dEy dEx ^ dHz
dx dy c dt , dx dy c dt
(4)
To construct the model of adiabatic guided modes (AGMs) we represent the solutions of (4) in the form of locally normal guided modes of a locally planar reference waveguide (see [7], [8], [27]), which in the method of asymptotic expansion take the form
t^î \ \—^ Es (x, y, z) r. , . >
E(x, y, z,t) = 2_^ ( ■ exp {iut — ikQ<p(y, z)} ,
s=n ( 10
(5)
H(x, y, z,t) = ^ exp {iut — ik0ip(y, z)} . (6)
s=n
(—iu
In the notation Es (x, y, z, t), Hs (x, y, z, t) the separation of x with a semicolon means the following assumption
and
dÊs (x,y,z) dÊs(x,y,z) 1 dÊs(x,y,z)
dy , dz u dx
dHs (x,y,z) dHs(x,y,z) 1 dHs (x,y,z)
dy , dz u dx
J = X,y,z
J = x,y,z,
(7)
(8)
where || || is the Euclidean norm and u is the circular frequency of the propagating monochromatic electromagnetic radiation.
Using the approach of the asymptotic expansion method with respect to the dimensional small parameter u-1 [23], [28], [29], we substitute expressions (5), (6) into the system of equations (4) and equate the coefficients at equal powers of the small parameter u-1. As a result, with the relations (7) and (8) taken into account, in the zero approximation of the method of asymptotic expansion with respect to small parameter we arrive at a system of homogeneous equations:
tfn Qy^n + Q^^n = "i^oneEn^
dHZ
dV
!Tz1±n
ikn „ Hn „ — ikneEï, ox
(9) (10)
dx 0 dy
+ i ko^H* = i ko£EzQ, (11)
-iko eq + iko ^ E0 = ikooHH0, (12)
-k°%EX° -it = lk°llH°' (13)
8Ey 8 D
+ l k°dlE° = ' k°11H" - (14)
By means of simple transformations, we reduce it to the form of interest to us. Namely, from the relation (12) we get the expression
T* = ^1 RZ Tjy
j0 = dyeH0 + dz eH0,
which we substitute into equations (13)-(14) that take the form
dEo dD ( dD1TJZ , dD1
= ' «0 * + )-' «-M, (15)
— ^ -' ^o m + £ m )- (16)
dx
From relation (9) we get the expression
Uo = dyti0 dz ti 0, and substitute it into equations (10)-(11), which take the form
-T- = <* + *>). (17)
^ = '^ -'%($£* ^)- (18)
Ultimately, Maxwell's equation reduced in the zero approximation to two algebraic equations and four differential equations of the first order.
At any fixed values of (y,z), equations (15)-(16), (17)-(18) take the form of a system of ordinary differential equations of the first order. Hereinafter, we deal with the zero approximation of the asymptotic expansion with respect to the small parameter; therefore, the index of the order of smallness is omitted:
dH* + ikodi (dDEZ - d^ \ + .ko£Ey = ^
dHV + ikod^ (d£EZ E .^EZ = ^
dEZ + ^^ y -^H* ]- ikofiHy = 0,
dEy + ikockp ( chpRy -
a h v h' )+' t°»H'=o.
In addition, for any fixed values of (y, z) the algebraic equations are valid
1 (9v rTy df rrA rrX _ 1 (ddf
rpx ___Try _TTZ TTX ___rpz _Tp
E° =e{TzH° -~dy H° ), H° =jA~dy E°
Remark 1. In the proposed form (5)-(6) of the desired solutions of
Maxwell's equations the quantities — and —- have the meaning of phase
constants of guided propagation of radiation in the directions Oy and Oz in all layers of the waveguide simultaneously, i.e., they determine the effective refractive index of the waveguide under study for the given adiabatic guided mode:
(%(y^))2 +{jfz(y-z))2 = < (V, *)• (!9)
From the analysis carried out, we conclude that with the accepted assumptions the zero approximation to the guided solution of Maxwell's equations is given by the following relations:
E(x\eq (x;y^)\
> = <i> } exp y,z)},
H( x,y,z,t)\ \H0 (x;y,z)\
with
dE0=--0 («
dx
dEo „. (
2
= lko )H° +lko ^ H0
dHo _ ,7. („. (2\ W (d{p\ f9^
(20)
0,
and
Fx = d^1uz +9V1 ttV rrx = 9^1 p s d(P 1 p y (21)
Eo =-^~eH° + Tz~eH°, Ho = d^vE° -Tz JE* (21)
as well as
(^j(y,z)) +(^(y,z)) =nff(y,z)-
For thin-film multilayer waveguide, consisting of optically homogeneous layers, from (2) the conditions of the electromagnetic field matching at the interfaces between the media follow:
n x E + nx E+ = 0,
n x H +ñx H+ = 0. From (3) the asymptotic conditions follow
¡7° E° H° H°
Jy, ^Z , y, z
0.
(22)
(23)
(24)
The system of equations (20), (24) for any fixed (y, z) defines a problem of
finding eigenvalues (Vp) (y,z) and eigenfunctions (EJy ,EJZ ,HJy ,HJZ y (y,z), v j
normalized to unity:
\E3y\ dx = 1,
I ->|2
\H¡\ dx = 1.
5. Algebraic model of adiabatic guided modes
In the case of a multilayer integrated optical waveguide consisting of homogeneous dielectric layers (possible, with complex permittivities) the relations (20), (22)-(24) are valid with generally non-horizontal interfaces between the layers (see Figure 2).
Figure 2. Structure of a multilayer thin-film waveguide
In this case, in each inner layer the four-dimensional system of ordinary differential equations with constant coefficients has a four-parameter system of fundamental solutions. Then the general solutions of equations (20) in the layers take the form
U(x; y, z) = C1l1 e1+x + C2l2e^-x + C3I3e^+x + C4$4e*-x,
Ey (X) =
-I(A2 -A4^ + Çv exP {-ko V¥ + Ï -^£X}
dp dp dy dz
exp
{ko V
dtf2 + dy *
d>f2 dz
Hex} (A1 + A3) (jl£
dy2 dy
dp dp dy dz
-(A2 + A4) (ye - ¥) exp {-ko V^ + - >1£X} _ Hz (x) = d^&p
dy dz
I(AX - A4)£ exp {ko V^ + ^ - V£x} V^^
dp dp '' dy dz
dp>2 dp>2
Hy(x) = (A2 + ^4) exp {-k0+ — V£X
where
dp2 dp2
Ez(x) = (A1 + ^3 ) exp {ko \\ + - fi£x } ,
T
u = (Hz(x; y, z); Ez(x; y, z); Hy(x; y, z); Ey(x; y, z))
and j± = ±ko.
The rest components of the electromagnetic field are calculated using formulas (21), and all the six components enter the explicit form of relations (22), (23).
In the substrate and cladding layers due to asymptotic conditions the two-dimensional systems of ordinary differential equations with constant coefficients have two-parametric systems of fundamental solutions.
The solution unique for all layers satisfies the matching conditions at the interfaces between the layers, i.e., determines a system of linear algebraic equations for indefinite coefficients, which specify the expression of particular solutions in the layers in terms of the systems of fundamental solutions.
At an arbitrary smooth interface between two dielectric media, described by the equation F(x,y,z) = x — h(y,z) = 0 (see Figure 3), conditions (22) and (23) take the form
[nxÊ] = 1Ev Tz F dz. —Ez — E9h, x d z:
n x H~] = {Hv Tz Hdh- ~Hz W —hz R 9Z
5 — Ed-hT
y xd y,
7 — H9^ T
y xdy,
(25)
X
= (i, -dhjdy, -dhjdz)1
■KM
Figure 3. Equation of a normal to the interface between the layers
in
In the expressions (25) and (26) the denominator y1+ (^) + (ff)
the expressions for the normal was omitted, since it is nonzero and coincides in both sides of equations (25), (26). It is worth noting that only two of three components of the obtained vectors (25) and (26) are linearly independent. Therefore, for writing the boundary conditions we will use the following expressions:
Az + Ax Qz
Ai — Ai —
y x d y
x=h(y,z)
i . 2 A 2 ^ h
= I Az +
= I A 2 — A 2 — y d y
x=h(y,z)
x=h(y,z)
x=h(y,z)
where A = {E, H}. Moreover, for a planar boundary F(x, y,z) = x - const the above expressions simplify to the following form:
Ai \ = A2 \
y I=const y i=const
Ai
'i=const
= A2 \
'i=const
The above relations should be completed with the asymptotic condition (24).
Thus, in each inner layer of a Allayer waveguide the four-dimensional general solution is parametrized by four indefinite coefficients, in the substrate and cladding layers at the expense of asymptotic conditions the number of coefficients is reduced by two in each layer. Therefore, there is a total of 4(k — 2) + 2 + 2 = 4(k — 1) coefficients. To each of ( k — 1) interfaces four
equations correspond, making the total of 4(k — 1). Thus, we have 4(k — 1) linear homogeneous algebraic equations for 4(k — 1) unknown coefficients.
Returning to the calculation of the electromagnetic field in the multilayer waveguide, we seek a nontrivial solution to the homogeneous system of linear algebraic equations with respect to the coefficient of expansion over the fundamental system of solutions in each dielectric layer:
M(h(y, z), Vh(y, z), <p(y, z), Vp(y, z)) A-0. (27)
To find nontrivial fields E, H at an arbitrary point (y, z) the condition of solvability should be satisfied for the system of homogeneous linear algebraic equations
det (M(h(y, z), Vh(y, z), <p(y, z), Vp(y, z))) - 0. (28)
Thus the system of homogeneous linear algebraic equations (27), nontrivially solvable under the condition (28), is an algebraic model of adiabatic guided modes in a smoothly irregular multilayer integrated optical waveguide. The roots of equation (28) are a set of eigenvalues and the solution of the system of equations (27) after substitution of each particular root is the corresponding eigenvector normalized by the condition Mil = 1.
6. Results
Thus, in the course of several stages of sequential reduction, we have formulated a number of problems of modeling the adiabatic guided modes in a smoothly irregular integrated optical multilayer waveguide.
Problem 1. The problem of finding eigenvalues and eigenfunctions within the AGM model
E(x,y,z,t)\ (E0 (x;y,z)\
. } = { - , , } exp Put — гp{y, z)} , (29)
H(x,y,z,t)) [H0(x;y,z)j
is formulated as follows:
9E0 _ (dV\ Try
£-d^ = —lk»(-dy) (Tz)H°—tko(£l1 dH0 ,, ((dV\2\ W ,7. fd(P\ f9(P
^ = - ^(fz) ) F° —< k0 (fz)(t]
(30)
Eigenfunctions (E0, E0, H0, H0)T (x; (Vip)2,) correspond to eigenvalues
The accompanying components of the electromagnetic field are calculated by the formulas:
Wx = 9 V1 TTZ+d V1 Try TTX = 9 P1 J?z 9 <P 1 p y (32)
Problem 2. The problem of finding eigenvalues of eigenfunctions within the AGM model
{È( X, y,z,t)} = { Êo( x; y, z)} exp [tut — i <p( y, z)} , [H(x, y, z, t)) [H0(x; y, z) J
is formulated as follows:
2 2 ^ + k2 (en — (Vif)2) E0 = 0, + k2 (e„ — (V<p)2) H0 - 0. (33)
Eigenfunctions (E(0, H(0, )T (x;(Vip)2) correspond to eigenvalues (Vip)2.
The accompanying components of the electromagnetic field are calculated by the formulas
^^ (f)(S) -0 — < (^(^ )
^"o m^T) *+'I-*m m
(34)
^^ (^—(M) —'«o m m *
Problem 3. Following the ideology of the cross section method, the steady-state regime of guided propagation of electromagnetic radiation within the frameworks of the AGM model, i.e., the solution of Maxwell's equations rotH =
ikeE, rotE = — ik^H with asymptotic conditions (E,H) (x) -> 0 is
sought in the form of a sum
(E, H)T(x,y,z) = Y, C$j(y, z) (E, H)T(x, fy),
where (Ê,H) (x,f3j) are solutions to the equations (30)-(31), (32) with additional condition:
V (Ê, H) T (x, V,z) = -y£ i PjC$j (y, z) (Ê, H)T ( x, fy),
j
where (3(y, z) = V p(y, z).
Problem 4. Following the ideology of the method of cross sections, the steady-state regime of guided propagation of electromagnetic radiation within the frameworks of the AGM model, i.e., the solution of Maxwell's equations rot(rotH) = k2e/iE, rot(rotE) = k02with asymptotic conditions
—*■ —*■ T —*■
(E, H) (x)-> 0 is sought in the form of a sum
(Ê,H)T (x,y,z) = Y, C$j (y, z) (E,H)T (x, Pj),
3
where (Ê, H) (x, fij) are solutions to the equations (33), (34) with additional condition:
V (Ê, H) T (x, y,z) = -y£ tfij c$j (y, z) (Ê, H)T (x, fij ), where fi(y, z) = Vp(y, z).
The single-mode regime is specified by ansatz
(E,H)^ (x,y,z) = C$j (y,z) (E,H)T (xjj)
with the additional condition
V (Ê, H)T_ (x, y, z) = -0jC$j(y, z) (Ê, H)T(x, fij),
where fi(y,z) = Vp(y,z).
Remark 2. Similar to the method of cross sections [24], [25], [30], the solution of the last equation has the form
(y,z) = fi^)exp yko
^ Vo, 0
which together with relations (29), (30)-(31), (32) concludes the description of the adiabatic character of the model under consideration.
(fiy (y,z)dy + fiz (y,z)dz)
7. Conclusion
The development of methods for the rigorous and approximate analysis of smoothly irregular integrated optical waveguides requires the development of new mathematical models of the corresponding objects, as well as the use of new methods for studying the problems arising in this case. The fundamental problem of electrodynamics of smoothly irregular waveguide three-dimensional (3D) structures is the development of stable methods and algorithms for solving the corresponding Maxwell's equations.
The paper considers an approach to the formulation of the problem of propagation of electromagnetic radiation in integrated optical smoothly irregular waveguide structures. Traditionally, similar problems are formulated based on Maxwell's equations. In this paper, only monochromatic radiation is considered, which is reflected in the dependence of solution on the frequency of propagating radiation. This type of time dependence of the solution allows considering a steady-state electrodynamic problem for electromagnetic radiation.
The problem of finding the eigenvectors (guided modes) and eigenvalues is considered using the model of adiabatic guided modes in the framework of the zero approximation of the asymptotic expansion for a planar regular three-layer optical waveguide. Considering that the permittivity and permeability are piecewise constant functions, the problem is solved in each subdomain with constant values of e, ^ with subsequent matching of solutions at the interfaces between the dielectric media. In each layer e, ^ have constant values, and the construction of the entire fields E, H requires setting and solving the problem of finding p(y, z).
The paper considers a class of smoothly irregular integrated-optical multilayer waveguides, the properties of which determine the characteristic features of waveguide propagation of monochromatic polarized light in them. An asymptotic approach to the description of this type of electromagnetic radiation is proposed, reducing the solutions of the system of Maxwell's equations to a form, which is expressed in terms of the solutions of a system of four ordinary differential equations and two algebraic equations for six components of the electromagnetic field in the zero approximation.
The multilayer structure of waveguides allows one more stage of reducing the model to a homogeneous system of linear algebraic equations, the nontrivial solvability condition of which specifies the relationship between the gradient of the phase front of radiation and the gradients of interfaces between thin homogeneous layers.
In the final part of the work, eigenvalues and eigenvector problems (differential and algebraic) describing adiabatic guided modes are formulated. The problem of describing the single-mode propagation of adiabatic guided modes is also formulated with emphasis on the adiabatic nature of the described approximate solution of Maxwell's equations.
Acknowledgments
The publication has been prepared with the support of the Russian Foundation for Basic Research (RFBR) according to the research project No 19-01-00645.
Appendix. Normal modes of a regular planar optical
waveguide
"Plane" guided modes
The monographic literature [15]-[19] widely describes normal modes of a regular planar dielectric waveguide propagating along the axis Oz from —to to to (and/or back). In the vertical direction (along the Ox axis) they have the structure of standing waves, while in the horizontal direction along the waveguide they are travelling waves and are considered (not quite correctly) steady-state ('invariant') in the direction transverse with respect to the propagation direction (along the axis Oy).
The vertical distribution of the electromagnetic field of a TE mode expressed in terms of the 'leading' transverse component of the electric field Ey is given by the equation
^Jk (x) + (£^ — p2 )Ey (X) = 0. (35)
Two other components of the electromagnetic field of the TE mode are expressed in terms of the leading one by the formulas
Hx = ——Ey, Hz = tt---j-^. (36)
These three relations can be derived from Maxwell's equations in the form (4).
For TM modes analogous relations that follow from Maxwell's equations have the form
^ () (*) + (W — fy2) Hv(*) = 0, (37)
dx \e dx
P 1 dHy Ex — —Hy, Ez — — —---—. (38)
S %Ko S qjx
Equations (35) and (37) in multilayer waveguides composed of uniform regular dielectric layers, in each subdomain of the real axis (at the intersection of each layer with the vertical axis) take the form of second-order linear ordinary differential equations with constant coefficients. Therefore, the most common form of solutions to these equations is obtained using a fundamental system of solutions. Consequently, the distributions of the corresponding leading components of the electromagnetic field are written in the form of expansions of general solutions in terms of (sin(\x), cos(xx)) or (exp [±i(\x)}) and the expansion coefficients of the particular solution are determined from the boundary conditions (22)-(23).
These solutions, in particular, can be obtained according to the following algorithm.
For guided TE modes we get a homogeneous system of linear algebraic equations (SLAE) with matrix M^(fi) with respect to variables A+,A+ ,A+,A-
A+ exp {jia-} - exp {ix{a-} + A- exp {—ix-a-} ,
C C
A+ exp {7ja-} = j1 exp {ix{a1} — A- exp {—ix{a-}) , ik0 k0
exp {ix{a2} + A- exp {—ix{a2} - A- exp {—^ca2} ,
i c
J1 exp {ix10,2 } — A- exp {—ix10-2 }) - — A- exp {—lJc0,2 } .
K0 IK0
The homogeneous SLAE (ME) A - 0 is nontrivially solvable if and only if its determinant is zero,
det (ME) - 0. (39)
For guided TM modes a system of homogeneous linear algebraic equations is obtained with the matrix M^(fi) for unknowns B+, B+ ,B+ ,B-, the solutions of which yield the values of the unknown amplitude coefficients
B+ exp {^Sa-} - B+ exp {ix{ai} + B- exp {—ix{a-} , c c
ls B+ exp {jia1} - (B+ exp {ix{a1} — B- exp {—ix{ai}) ,
ikoss ko£1
B+ exp {ix{a2} + B- exp {-ix{02} = B- exp {-ji02} ,
j j
-j^ (B+ exp {ix{o.2} - B- exp {-ix{02}) = exp a2} .
The homogeneous SLAE (MH) B - 0 is nontrivially solvable if its determinant equals zero
det (MH) - 0. (40)
Here A and B are the coefficients of expansion of TE and TM modes in terms of the fundamental system of solutions, respectively.
Equations (39) and (40) are equivalent to the dispersion relations in the trigonometric form
Xfd = arctg (^j + arctg (^ j+mn (41)
for guided TE modes and
i F.^ m\
+ arctg
\£c XT/
for guided TM modes
Xfd = arctg (^ j + arctg (^ j + mn (42)
If the expressions for the electromagnetic field strength are presented in complex form with the described dependence on the rest coordinates and time taken into account, i.e., in the form
(Ey, Hx,HZ)T (x, y, z, t) = (A1 ,A2, A3)T(x, y, z) exp [iut — ieE(x, y, z)} , (Hy,EX,EZ) (x,y,z,t) = (Bi,B2,B3)T(x,y,z) exp [iut — ieH(x,y,z)}
with real-valued amplitude A(B) and phase (pH), then the phase remains constant in time along x, defines a travelling wave along z, and is constant along y. In other words, the phase front of the described solution to Maxwell's equations is "planar" (i.e., linear in the yOz-plane) and defines a "plane" ( i.e., linear in the yOz-plane) wave. There exist "forward" and "backward" waves travelling in opposite directions along the z-axis. The can be identified with the "plane" (in the yOz-plane) waves emitted by infinitely remote sources.
"Cylindrical" guided modes
Now let us proceed to the guided modes of a regular planar dielectric waveguide, excited by a source linear along the Ox-axis and point-like in the yOz-plane, localized at the point (0,zo). The structure of the modes along the Ox-axis completely coincides with that of "plane" guided and leaky TE and TM modes. Let us analyze the structure of cylindrical guided modes propagating in the yOz-plane
^J (x,y,z,t)= (^J (x) exp [iut — ik0/3r} ,
where r2 = y2 + (z — z0)2; y = r sin 0, z = z0 + r cos Q. In the yOz-plane circular fronts propagate from the origin of polar coordinates y = r sin d,
z = z0 + r cos Q.
Thus, in both cases the solutions for the normal modes are written as
E\ exp {iut — ik0 tp(y,z)}
h) (x,y,z,t)= UJ (x)—tws—,
ft
where:
— (p(y, z) = fiz for the modes from an infinitely remote source and
— ip(y, z) = fi^y2 + (z — z0)2 for the modes from a localized source.
References
[1] A. A. Samarskii and A. N. Tikhonov, "Representation of the field in
a waveguide as the sum of the TE and TM fields [O predstavlenii polya v volnovode v vide summy polej TE i TM]," Zhurnal tekhnicheskoy fiziki, vol. 18, no. 7, pp. 959-970, 1948, in Russian.
[2] A. G. Sveshnikov, "A substantiation of a method for computing the propagation of electromagnetic oscillations in irregular waveguides," USSR Computational Mathematics and Mathematical Physics, vol. 3, no. 2, pp. 413-429, 1963. DOI: 10.1016/0041-5553(63)90027-2.
[3] V. V. Shevchenko, "Spectral decomposition in eigen- and associated functions of a nonselfadjoint problem of Sturm-Liouville type on the entire axis [O spektral'nom razlozhenii po sobstvennym i prisoedinennym funkciyam odnoj nesamosopryazhennoj zadachi tipa SHturma-Liuvillya na vsej osi]," Differ. Uravn., vol. 15, no. 11, pp. 2004-2020, 1979, in Russian.
[4] L. N. Deryugin, A. N. Marchuk, and V. E. Sotin, "Properties of planar asymmetrical dielectric waveguides on a substrate of dielectric [Svojstva planarnogo asimmetrichnogo dielektricheskogo volnovoda na podlozhke iz dielektrika]," Izv. Vyssh. Uchebn. Zaved., Ser. Radioelektron., vol. 10, no. 2, pp. 134-141, 1967, in Russian.
[5] A. G. Sveshnikov, "The incomplete Galerkin method [Nepolnyj metod Galerkina]," Dokl. Akad. Nauk SSSR, vol. 236, no. 5, pp. 1076-1079, 1977, in Russian.
[6] A. A. Bykov, A. G. Sveshnikov, and M. K. Trubetskov, "Reduced Galerkin's method application to calculations of eigenwaves in open waveguides [Primenenie nepolnogo metoda Galerkina dlya rascheta sobstvennyh voln otkrytyh volnovodov]," Matem. Mod,., vol. 3, no. 7, pp. 111-123, 1991, in Russian.
[7] B. Z. Katsenelenbaum, Theory of Irregular Waveguides with Slowly Varying Parameters [Teoriya neregulyarnyh volnovodov s medlenno menyayushchimisya parametrami]. Moscow: Akad. Nauk SSSR, 1961, in Russian.
[8] V. V. Shevchenko, Continuous Transitions in Open Waveguides [Plavnye perekhody v otkrytyh volnovodah]. Moscow: Nauka, 1969, in Russian.
[9] L. A. Sevastianov and A. A. Egorov, "Theoretical analysis of the waveguide propagation of electromagnetic waves in dielectric smoothly-irregular integrated structures," Optics and Spectroscopy, vol. 105, no. 4, pp. 576-584, 2008. DOI: 10.1134/S0030400X08100123.
[10] A. A. Egorov and L. A. Sevastianov, "Structure of modes of a smoothly irregular integrated optical four-layer three-dimensional waveguide," Quantum Electronics, vol. 39, no. 6, pp. 566-574, 2009. DOI: 10.1070/ QE2009v039n06ABEH013966.
[11] A. A. Egorov, K. P. Lovetskiy, A. L. Sevastianov, and L. A. Sevastianov, "Simulation of guided modes (eigenmodes) and synthesis of a thin-film generalised waveguide Luneburg lens in the zero-order vector approximation," Quantum Electronics, vol. 40, no. 9, pp. 830-836, 2010. DOI: 10.1070/QE2010V040N09ABEH014332.
[12] A. A. Egorov, A. L. Sevast'yanov, and L. A. Sevast'yanov, "Stable computer modeling of thin-film generalized waveguide Luneburg lens," Quantum Electronics, vol. 44, no. 2, pp. 167-173, 2014. DOI: 10.1070/ QE2014v044n02ABEH015303.
[13] A. A. Egorov, A. L. Sevastyanov, E. A. Ayryan, and L. A. Sevastyanov, "Stable computer modeling of thin-film generalized waveguide Luneb-urg lens [Ustojchivoe komp'yuternoe modelirovanie tonkoplenochnoj obobshchennoj volnovodnoj linzy Lyuneberga]," Matem. Mod., vol. 26, no. 11, pp. 37-44, 2014, in Russian.
[14] A. S. Il'inskii, V. V. Kravtsov, and A. G. Sveshnikov, Mathematical Models of Electrodynamics [Matematicheskie modeli elektrodinamiki]. Moscow: Vyssh. Shkola, 1991, in Russian.
[15] M. J. Adams, An Introduction to Optical Waveguides. New York: Wiley, 1981.
[16] T. Tamir, "Guided-Wave Optoelectronics," in Integrated Optics, T. Tamir, Ed. Berlin: Springer-Verlag, 1990.
[17] A. W. Snyder and J. D. Love, Optical Waveguide Theory. New York: Chapman and Hall, 1983.
[18] D. Markuze, Theory of Dielectric Optical Waveguides. New York: Academic Press, 1974.
[19] M. Barnoski, Introduction to Integrated Optics. New York: Plenum Press, 1974.
[20] A. I. Neishtadt, "On the accuracy of conservation of the adiabatic invariant," Journal of Applied Mathematics and Mechanics, vol. 45, no. 1, pp. 58-63, 1981. DOI: 10.1016/0021-8928(81)90010-1.
[21] I. V. Gorelyshev and A. I. Neishtadt, "On the adiabatic perturbation theory for systems with impacts," Journal of Applied Mathematics and Mechanics, vol. 70, no. 1, pp. 4-17, 2006. DOI: 10. 1016/j . jappmathmech.2006.03.015.
[22] V. M. Babich and V. S. Buldyrev, Asymptotic Methods in Short-Wavelength Diffraction Theory, Alpha Science Series on Wave Phenomena. Harrow, UK: Alpha Science International, 2009.
[23] Y. A. Kravtsov and Y. I. Orlov, Geometrical Optics of Inhomogeneous Media. Berlin: Springer-Verlag, 1990.
[24] B. Z. Katsenelenbaum, "Irregular waveguides with slowly varying parameters [Neregulyarnye volnovody s medlenno menyayushchimisya parametrami]," Doklady Akademii Nauk SSSR, vol. 102, no. 4, p. 711, 1955, in Russian.
[25] B. Z. Katsenelenbaum, "A contribution to the general theory of non-regular wave guides [K obshchej teorii neregulyarnyh volnovodov]," Dokl. Akad. Nauk SSSR, vol. 116, no. 2, pp. 203-206, 1957, in Russian.
[26] A. I. Neishtadt, "Propagation of rays in smoothly irregular waveguides and perturbation theory of Hamiltonian systems [Rasprostranenie luchej v plavno neregulyarnyh volnovodah i teoriya vozmushchenij gamil'tonovyh sistem]," Izv. vuzov. Radiofizika, vol. 25, no. 2, pp. 218226, 1982, in Russian.
[27] A. L. Sevastyanov, "Single-mode waveguide spread of light in a smooth irregular integral optical waveguide [Komp'yuternoe modelirovanie polej napravlyaemyh mod tonkoplenochnoj obobshchennoj volnovodnoj linzy Lyuneberga]," in Russian, Ph.D. dissertation, Peoples' Friendship University of Russia, Moscow, 2010.
[28] S. Solimeno, B. Crosignani, and P. DiPorto, Guiding, Diffraction, and Confinement of Optical Radiation. Orlando, FL: Academic Press, 1986.
[29] M. Kline and I. W. Kay, Electromagnenic Theory and Geometricak Optics. New York: Wiley (Inter-science), 1965.
[30] M. V. Fedoryuk, "A justification of the method of transverse sections for an acoustic wave guide with nonhomogeneous content," Mathematical Physics, vol. 13, no. 1, pp. 162-173, 1973. DOI: 10.1016/0041-5553(74) 90012-3.
For citation:
A. L. Sevastianov, Asymptotic method for constructing a model of adiabatic guided modes of smoothly irregular integrated optical waveguides, Discrete and Continuous Models and Applied Computational Science 28 (3) (2020) 252-273. DOI: 10.22363/2658-4670-2020-28-3-252-273.
Information about the authors:
Sevastianov, Anton L. — Candidate of Physical and Mathematical Sciences, assistant professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9550927, ORCID: https://orcid.org/0000-0002-0280-485X, Scopus Author ID: 50462435500)
УДК 519.61; 519.62; 621.372.82
DOI: 10.22363/2658-4670-2020-28-3-252-273
Асимптотический метод построения модели адиабатических волноводных мод плавно-нерегулярных интегрально-оптических
волноводов
А. Л. Севастьянов
Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, 117198, Россия
В работе рассмотрен класс плавно нерегулярных интегрально-оптических многослойных волноводов, свойства которых определяют характерные черты волноводного распространения в них монохроматического поляризованного света. Предложен асимптотический подход к описанию данного вида электромагнитного излучения, в результате которого решения системы уравнений Максвелла редуцируется к такому виду, который выражается через решения системы четырёх обыкновенных дифференциальных уравнений и двух алгебраических уравнений для шести компонент электромагнитного поля в нулевом приближении. Градиент фазового фронта адиабатической волноводной моды удовлетворяет уравнению эйконала относительно эффективного показателя преломления волновода относительно данной моды.
Многослойная структура волноводов позволяет произвести ещё один этап редукции системы уравнений модели к однородной системе линейных алгебраических уравнений, условие нетривиальной разрешимости которой задаёт связь градиента фазового фронта излучения с градиентами поверхностей раздела между тонкими однородными слоями.
В завершающей части работы сформулированы задачи (дифференциальная и алгебраическая) на собственные значения и собственные векторы для описания адиабатических волноводных мод. Приведена также формулировка задачи описания одномодового режима распространения адиабатических волноводных мод, подчёркивающая адиабатический характер описываемого приближенного решения уравнений Максвелла.
Ключевые слова: плавно нерегулярные интегрально-оптические многослойные волноводы, задачи на собственные значения и собственные векторы, одномодовый режим распространения адиабатических волноводных мод
