# Propagation Properties of Airy Beam through Periodic Slab System with Negative Index Materials

Authors: Long Jin; Xingqiang Zhang
Journal: International Journal of Optics
 (2018)

---

## Abstract

Based on light transfer matrix and electric field vector equation, the evolution of Airy beam propagating in periodic slab system with three negative index materials (NIMs) and its transmission mechanism are investigated. The intensity profiles on emergent surface of periodic slab system and side view of Airy beam propagating in each right handed material (RHM) and double negative material (DNM) unit including lossless and losses DNMs are discussed. It is revealed that the self-recovery Airy beam can be achieved in long distance by using lossless periodic slab system as long as the negative refractive indexnl=-nr and each unit length L=Z. As to losses slab system contained DNMs, the smaller the collision frequencies are, the better the Airy beam quality is formed. It is expected that the proposed manner of beam transmission and corresponding conclusions can be useful for extension applications of optical control, especially for optical communication and optical encryption technique.


---

## Body

## 1. Introduction

Airy beam, which makes up one of the important nondiffracting waves, was first described by Berry and Balazs in 1979 by solving Airy solution of the paraxial wave equation [1]. Recently, Siviloglou et al. observed the 1D and 2D Airy phenomena through experiment in 2007 [2]; they utilized spatial light modulator combined with Fourier transform technique to generate the accelerating finite energy Airy beam. From then on, the particular properties of this beam, approximate nondiffraction, parabolic path of the central lobe, self-accelerating in the transverse, and self-recovery, have been an attractive topic for both theoretical and applied study in modern optics [3]. At the same time, a considerable number of frontier and interdisciplinary applications, such as optical microparticles manipulation [4], formation of Airy beam lasers [5], generation of optical bullets [6], routing of surface plasmon polaritons [7], generation of curved plasma channels [8], and creating curved electron beams without any external field [9], have been took advantage of the Airy function and Airy beam technology.As far as the width and depth are concerned, optical film and related technology had been significantly developed with the progress of science and technology. The preparation process of optical film is using different refractive index materials to coat a periodic media structure [10]. These optical components possessing the properties of spectral transmittance, spectral reflection, spectral absorption, and altering the state of polarization and/or phase of light are widely used in manufacturing high reflectance film, antireflective film, interference filter, and so on [11]. Meanwhile, metamaterials so-called negative index materials (NIMs) have mushroomed in many aspects of material physics and chemistry science [12]. Double negative materials (DNMs) with both permittivity ε and permeability μ are known as left handed materials (LHMs), whose extraordinary electromagnetic (EM) properties were first presented by Veselago [13]. Beyond that, there are single negative materials (SNMs) among which only one of the EM parameters is negative; these materials include the mu negative media (MNGM, ε>0, μ<0) and the epsilon negative media (ENGM, ε<0, μ>0) [14, 15]. There exist abnormal negative refraction effect [16], inverse Cerenkov radiation, reverse Doppler effect [17], and negative radiation pressure [18] in NIMs because the electromagnetic wave group velocity and phase velocity of these materials are oriented in the opposite direction; based on these principles, optical film contained NIMs has broad application prospects in microwave communication, large capacity storage medium, and superresolution imaging.In the letter, we mainly research the evolution of Airy beam propagation in optical film with DNMs, but what we want to emphasize is that our optical film model is 1D periodic slab system, with its period number being very small and each unit is measured in metre in order to conveniently observe beam transmission mechanism. The rest of the paper is organized as follows: firstly, the general propagation formulae of the Airy beam propagating through the paraxial optical system are calculated in space domain by transfer matrix and the Collins integral. Furthermore, the intensity profiles on emergent surface of periodic slab system and side view of Airy beam transmission in each RHM and DNM unit including lossless and losses DNMs are discussed. Finally, some important conclusions are summarized. It is expected that the proposed manner of wave transmission and the corresponding conclusions can be useful for extension applications of optical control, especially for optical communication and optical encryption technique.


## 2. Physical Model and Formulae

The schematic diagram of 1D periodic slab system with NIMs is shown in Figure1; the medium layer with black colour stands for RHM and the light green is NIM. Here we choose isotropic DNM with its relative permittivity and permeability given by the Drude model [19]:(1)εrω=1-ωpe2ωω+iγe,μrω=1-ωpm2ωω+iγm,where ωpe and ωpm are the electric plasma frequency and the magnetic plasma frequency. γe, γm are the electric and magnetic collision frequencies, respectively. The frequency v=ω/(2π) is measured in 102 THz, and ω is angular frequency. DNM’s refractive index nl=-εr(ω)·μr(ω); the length of each DNM unit is L. With regard to RHM, the refractive index and each unit length are nr and Z, respectively.Figure 1
Light path schematic diagram of periodic slab system with DNMs.Considering an Airy beam that propagates fromO point along the direction of z, the electric field at z=0 (O point) can be described as [20, 21](2)Ex1,y1=Aix1x0exp⁡a0x1x0×Aiy1y0expb0y1y0,where Ai(·) denotes the Airy function, x0 and y0 are beams horizontal size and vertical size, respectively, a0 and b0 represent exponential modulation factors, and E(·) is electric field intensity on incident surface. The propagation of beam through paraxial periodic slab system with NIMs can be described by the Collins integral [22]:(3)E2x,y,z=-iλBeikz∬s1E1x1,y1,0eik/2BAx12+y12+Dx2+y2-2x1x+y1ydx1dy1,where λ=1.55 μm is the wavelength of Airy beam, k=2π/λ is the wave number. A, B, C, and D are the transfer matrix elements of the total paraxial optical systems T:(4)T=ABCD.Now, let us discuss some useful matrices of optical medium. When beam is traveling x distance in DNM or RHM, the transfer matrix can be written as(5)Mx=1x01.While beam propagates from RHM into DNM slab, the transfer matrix form in RHM-DNM interface is given by(6)Mnr,nl=100nrnl.Likewise, DNM-RHM interface transfer matrix is(7)Mnl,nr=100nlnr.By cross-producing these matrices in accordance with the order of Airy beam passing through periodic slab system, we can obtain the general propagation formula of the Airy beam propagating through the paraxial optical systems in the space domain:(8)Ex=1AAixAx0-B24A2k2x04+ia0BAkx02·exp-ikCx22A·expa0xAx0-a0B22A2k2x04-iB312A3k3x06+ia02B2Akx02+iBx2A2kx03,Ey=1AAiyAy0-B24A2k2y04+ib0BAky02·exp-ikCy22A·expb0yAy0-b0B22A2k2y04-iB312A3k3y06+ib02B2Aky02+iBy2A2ky03,Ex,y=Ex×Ey.And the light intensity is easily obtained by using complex conjugate of E(x,y,z):(9)I=n2cμ0Ex,y,z2∝Ex,y,z2.


## 3. Results and Discussion

To delve into the evolution of Airy beam passes through periodic slab systemABm, where A and B are RHM and DNM, respectively, and m=2 is period number of this system, we carry out Airy beam transmission numerical calculation and display some simulation findings. Here we assume x0=y0=0.2 mm, a0=b0=0.1, each DNM is lossless, and related parameters are set as γe=γm=0, v=2×102 THz, ωpe=ωpm=2πν×2, 2πν×2.55, and 2πν×3.1, resulting in nl=-1.00, −1.55, and −2.10, respectively, nr=1.55 fixed, L=Z=0.25 m. Firstly, we depict transverse normalized intensity distribution of incident Airy beam and three kinds of emerging beams on output surface. It is easily found that the incident beam distribution and emerging beam transverse profiles are completely the same when DNM’s plasma frequencies ωpe=ωpm=2πν×2.55 (i.e., abs(nl)=nr=1.55) from Figures 2(a) and 2(b), and the normalized central main lobe intensity slightly decreases when ωpe=ωpm=2πν×3.1, which was discerned from the colour bar in Figure 2(c); while ωpe=ωpm=2πν×2 and emerging central maximum intensity is attenuated significantly, the quantity value is 0.57 compared to that of ωpe=ωpm=2πν×2.55.Figure 2
Transverse normalized intensity distribution of incident Airy beam and three kinds of emerging beams on output surface. (a) Incident Airy beam distribution, (b) emerging Airy beam distribution, DNM’sωpe=ωpm=2πν×2.55, (c) emerging Airy beam distribution, DNM’s ωpe=ωpm=2πν×3.1, and (d) emerging Airy beam distribution, DNM’s ωpe=ωpm=2πν×2.
(a)
(b)
(c)
(d)The intensity evolution of Airy beam propagating in periodic slab systemAB2 is deeply investigated in order to show how wave transfers from RHM to DNM; the simulation results are depicted in Figure 3; green arrow lines show the material boundaries; this reveals that the maximum intensity decreases and intensity distribution generates self-bending as the propagation distance increases while the Airy beam propagates in the RHM; when passing through the RHM-DNM interface, the beam undergoes negative refraction; in other words, both incident beam and refracted beam are on the same side of the interface normal; subsequently, the beam experiences inverse parabolic trajectory in DNM, which indicates the self-bending in RHM can be compensated and refocused again when wave passes through the DNM. In Figure 3(a), as DNM’s nl=-nr=-1.55, the transverse pattern could go back to its original feature at 1st end of DNM and operate the stable and excellent beam quality; in consequence, we draw a conclusion that the intensity profiles of incident surface and emergent surface are completely the same on condition that L=Z no matter what the value of period number m is; this case is called Airy beam which can just be compensated by each DNM, and related principle diagram is presented in Figure 3(b). The calculation result in Figure 3(c) reveals that the self-bending in RHM cannot be completely compensated by DNM when nl=-2.10<-nr and L=Z=2.5 m; this effect accumulates accompanied by increasing of periodic slab system period number m and severely impacts the Airy beam’s quality of imaging and transmission. We also depict corresponding principle diagram in Figure 3(d) to describe beam transmission characteristics under this circumstance. Now an interesting question naturally arises: How long is the DNM needed to just compensate self-bending when nl=-2.10?Figure 3
Side view of Airy beam propagating in each RHM and DNM unit successively with three kinds of DNMs and corresponding principle diagrams. (a) DNM’sωpe=ωpm=2πν×2.55, (b) principle diagram where beam can just be compensated by each DNM, (c) DNM’s ωpe=ωpm=2πν×3.1, (d) principle diagram where beam cannot be completely compensated by each DNM, (e) Airy beam transmission in quasiperiodic slab system, Z=2.5, L=105/31, (f) DNM’s ωpe=ωpm=2πν×2, (g) principle diagram where beam is overcompensated by each DNM, and (h) Airy beam transmission in quasiperiodic slab system, Z=2.5, L=50/31.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)When Airy beam passes through the 1st RHM and DNM, the transfer matrixT can be described as(10)T=ABCD=100nl×1L01×100nrnl×1Z01×1001nr.Substituting Z=2.5 m, nr=1.55, and nl=-2.10 into (10), the analytical transfer matrix element B can be formulated as(11)B=5031-1021×L.Hence, the necessary distance for just compensating self-bending is L=105/31 m in DNM. The side intensity evolution of Airy beam propagating in quasiperiodic slab system with Z=2.5 m and L=105/31 m alternately is plotted in Figure 3(e); it is clearly found that the output Airy beam gets back to its original feature at the end of each DNM, which verifies our theoretical analysis. On the contrary, the length of DNM to just compensate self-bending when nl=-1.00 is smaller than 2.5 m. By substituting Z=2.5 m, nr=1.55, and nl=-1.00 into (10), analytical transfer matrix element B can also be written as(12)B=5031-Land the quantity value is L=50/31 m; therefore, the beam intensity self-bends again in the rest of the of 1st DNM because L=2.5 m fixed; when passing through the RHM-DNM interface, the beam undergoes negative refraction as well and then experiences inverse parabolic trajectory even in the next RHM in Figure 3(f); Pendry once explained this overcompensation phenomenon in detail through photoelectromagnetic theory [23], and the principle diagram is depicted in Figure 3(g).As we all known, the nondiffracting beams could be extended to infinite space and have infinite energy; they are generally unobtainable in the laboratory. However, from the above discussion, we can come to the conclusion that the self-recovery and self-reproducing Airy beam can be achieved in long distance by using periodic slab system with DNMs as long asnl=-nr,L=Z; the result is very important to form stable and excellent beam quality. When nl≠-nr, the larger nl requires smaller L to just compensate self-bending in RHMs; otherwise, the conclusion is also the opposite. L≠Z manifests that the slab system is quasiperiodic, which increases the difficulty of materials production process. Moreover, the quantitative relation between L and nl of these quasiperiodic slab systems remains to explore and perfect further.In practice, losses in DNMs are inevitable; in this subsection, emerging intensity profiles with four different collision frequenciesγe are displayed in Figure 4. It is clearly acquired that the larger the collision frequencies γe are, the more serious the Airy beam deformation are formed. When γe=20 THz, only main lobe of Airy beam existed; other side lobes are vanished or significantly weakened from Figure 4(d), which indicates that the beam is seriously distorted. Similarly, side views of Airy beam propagating in each RHM and DNM unit successively changing with two kinds of collision frequencies γe are exhibited in Figure 5. We can easily find that the beam propagation is rarely affected when γe=0.2 THz from Figure 5(a); however, as γe gradually increases, the beam’s attenuation and diffraction become more and more obvious; when γe=20 THz, the interaction and diffraction between main lobe and side lobes lead to a distorted Airy beam; at this point, the bit error rate of communication system would greatly increase as to optical encryption. Gris-Sánchez et al. have designed and made a periodic optical fiber that guides an approximate Airy pattern as one of the beam guided modes recently; the fiber’s attenuation was 11.0 dB/km at 1.55 μm wavelength [24]. Here we estimate maximum collision frequencies value allowed by using their experimental result. Beam loss equation can be defined as [25](13)αλ=-10Z+L2lgIoutIin,where α(λ) is loss coefficient, Iin; Iout denote incident beam intensity and emerging beam intensity, respectively. By substituting Z=L=2.5 m and α(λ)=11.0 dB/km into (13), we obtain Iout/Iin=0.975. When DNM’s collision frequencies γe and γm equal 0.2 THz, the corresponding Iout/Iin=0.980≈0.975; hence, we had better choose γe and γm no more than 0.2 THz (0.1% of the frequency v) in order to achieve accurate transmission of optical communications.Figure 4
Emerging beam intensity profiles with four different collision frequenciesγe. (a) γe=0.2 THz, (b) γe=2.0 THz, (c) γe=10 THz, and (d) γe=20 THz.
(a)
(b)
(c)
(d)Figure 5
Side view of Airy beam propagating in each RHM and DNM unit changing with two different collision frequenciesγe. (a) γe=0.2 THz; (b) γe=20 THz.
(a)
(b)


## 4. Conclusion

In summary, we have studied the evolution of Airy beam propagating in periodic slab system with different negative index materials and its transmission mechanism. The emerging intensity profiles on output surface of slab system and side view of Airy beam transmission in each RHM and DNM unit including lossless and losses DNMs have been discussed. It is revealed that the self-recovery Airy beam can be achieved in long distance by using lossless periodic slab system with negative index materials as long asnl=-nr, L=Z; the conclusion is very important to forming stable and excellent beam quality. As to losses periodic slab system contained DNMs, the larger the collision frequencies γe and γe are, the more serious the Airy beam deformation is formed; hence, we had better choose γe and γm no more than 0.1% of the frequency v in order to achieve accurate transmission of optical communications. It is expected that the proposed manner of wave transmission and the corresponding conclusions can be useful for the extension applications of optical control, especially for optical communication and optical encryption technique.



---

*Source: 9478483-2018-01-16.xml*