KnE Energy | VII International Conference on Photonics and Information Optics (PhIO) | pages: 235–240


1. Introduction

The one dimensional array (or chain) of the waveguides in the linear approximation represents the system for discrete diffraction realization. In this paper we will consider the binary linear waveguide array. The unit cell contains two kind of the waveguide, i.e., A-type and B-type. In slowly varying envelopes of the electric field of the continue wave radiation in waveguide is governed by the system of equations for coupled waves [1,2]. If we would like to consider only discrete diffraction [2,3] then the nonlinear properties of the waveguides may be neglected. It cab be considered as the first approximation.

The system of the equations describing the electromagnetic wave propagation in the coupled wave approximation takes the following form


(1) where An and Bn are normalized amplitudes of the electric field in waveguide from n -th unit cell, ζ is the normalized coordinate [1]. We assume that the phase mismatch is zero.

2. Analytical solution of the base equations

To obtain the solution of the system of equations (1) we can use the generation function method. Let us introduce the following functions


Using the equations (1) we can find the equations for the generation functions PA and PB :



From (1) it follows equation ζ2PA+Ω2PA=0 , where Ω2=4cos2(y/2) . The solution of this equation takes the form


The expression for the generation function PB follows from the first equation of (2). If the initial conditions for the amplitudes An and Bn are known, we can define the initial conditions for the generation function


The initial conditions for PA and PB allows us to determine integration constants C1 and C2 . Thus the generation function can be written as




(4) where Ω=2cos(y/2) . Using the orthogonality condition


the solution of the initial system of equations (1) can be written as



3. Particular examples of the field distribution over waveguides

Let us consider the following initial condition that corresponds to strong focusing radiation at ζ=0 : An(0)=A0δn0 and Bn(0)=0 . In this case PA0=A0 and PB0=0 . By the use the expressions (3), (4) and (5) we can write


To compute the integrals in these expressions we may use the formula by Anger, which in this case is looking like





Substitution of the equation (6) into (5) results in the following expression (here we will use the term η=2ζ )


Hence, at n=0 we have A0(ζ)=A0J0(2ζ) . At n1 it follows that



The second equation of (1) may be used to obtain the amplitudes Bn(ζ) . By using the expression (8) for An(ζ) , one can write




the equation for Bn(ζ) can be rewritten as


It follows that



Thus the distribution of the field amplitudes over waveguides in array under considered initial conditions is presented by the expressions (8) and (9). These equations describe the discrete diffraction in binary waveguide array under the condition of the strong focusing at ζ=0 .

Let us consider the case where the both waveguides in unit cell n=0 are illuminated. Initial conditions are following An(0)=A0δn0 and Bn(0)=B0δn0 . In this case we have PA0=A0 and PB0=B0 . The generation function can be written as


Using the (5) we can write


Two incoming here the integrals have been found previously


The third integral can be defined by the similar way. It results in following expression


So can immediately write down expressions for the distributions of the electric field amplitudes over waveguides for the selected initial conditions


where n=±1,±2,... . The third terms in these expressions represent the interference phenomena in waveguide array.

4. Conclusion

This distribution of field strengths describes discrete diffraction in a binary array of waveguides. For more complex cases of the initial conditions the expressions for the electric fields in waveguides contain terms that account for interference fields in the neighboring waveguides. It should be noted that if the initial conditions are selected as An(0)=(-1)nA0 and Bn(0)=(-1)nB0 , then the diffraction is absent. The equations (1) show that under these conditions the fields in waveguides are invariants. However, the flat band [5,6,7] in the spectrum of the linear waves is absent, as the number of nodes in the unit cell is less than three. In the spectrum there are only two branches that meet dispersive waves along the chain.


We are grateful to Prof. I.R. Gabitov, Dr. E.V. Kazantseva, Dr. V. A. Patrikeev, Dr. N.V. Bikov for enlightening discussions. This research was supported by the Russian Basic Research Foundation (Grant No 15-02-02764).



S. Somekh, E. Garmire, A. Yariv, H.L. Garvin, and R.G. Hunsperger, “Channel Optical Waveguide Directional Couplers,” Appl. Phys. Letts. vol. 22, no 2, pp. 46-47, 1973


A.I. Maimistov, and I.G. Gabitov. “Nonlinear optical effects in artificial materials”, Eur. Phys. J. Special Topics. vol. 147, no. 1, pp. 265-286, 2007.


T. Pertsch, T. Zentgraf, U. Peschel, A. Brauer, and F. Lederer, “Anomalous Refraction and Diffraction in Discrete Optical Systems,” Phys.Rev.Lett. vol. 88, 093901 (4 pp), 2002


O V Borovkova, V E Lobanov, A P Sukhorukov, and A K Sukhorukova, “Controllable discrete diffraction in cascade-induced waveguides", Quantum electron, vol. 39, no. 11, pp. 1050–1054, 2009.


S. Flach, D. Leykam, J.D. Bodyfelt, P. Matthies, and A.S. Desyatnikov, “Detangling flat bands into Fano lattices”, Europhys. Lett. vol. 105, 30001 (6 pp), 2014.


A.I. Maimistov, and I.R. Gabitov, “Optical flat bands in 2D waveguide arrays with alternating sign of refraction index”, J.Phys. Conference Series, vol. 714, 012013 (4 pp), 2016.


A.I. Maimistov and V.A. Patrikeev, “Electromagnetic Wave Propagation in a Quasi-1D Rhombic Linear Optical Waveguide Array”, J.Phys. Conference Series, vol. 737, 012008 (4 pp). 2016, doi:10.1088/1742-6596/737/1/012008



  • Downloads 11
  • Views 65



ISSN: 2413-5453