Spectra and field distribution of photonic-crystal structure with inclusions of metal nanoparticles I.A. Glukhov1, S.G. Moiseev1,2 1 Ulyanovsk State University, 42 Lev Tolstoy Str., 432017, Ulyanovsk, Russia 2 Kotelnikov Institute of Radio Engineering and Electronics of the Russian Academy of Sciences, Ulyanovsk Branch, 48/2 Goncharov Str., 432011 Ulyanovsk, Russia Abstract Transmittance and reflectance spectra as well as field distribution in 1D photonic-crystal structure with embedded dielectric layer and monolayer of metal nanoparticles are characterized. The influence of plasmonic monolayer location on defect modes of photonic-crystal structure is demonstrated with respect to domains of field confinement in the cavity area. Keywords: nanoplasmonics; photonic-crystal structure; defect mode; field localization 1. Introduction In recent years 1D photonic-crystal structures (PCS) created on the basis of different materials are of a special interest to researchers. Owing to periodic modulation of refractive index, photonic spectrum of these structures has a band gap, in which incident radiation is practically totally reflected. This property is critical for practical use as it enables to control optical radiation in data-transmission systems and in laser technology. Particularly remarkable is Fabry-Perot microresonator-like structure composed of two Bragg reflectors with defect layer there between. Defect layer in such-like structure plays the role of optical microcavity (microresonator) on which electromagnetic radiation can be localized. This can add to material-radiation interaction effects. Varying geometrical and physical properties of the structure it is possible to control spectral characteristics of PCS [1, 2] that enables to improve considerably their functionality. For example, through breakdown of the structure periodicity or using materials with controlled properties (non-linear, resonant, magnetogirotropic) photonic spectrum of PCS can be modified considerably. Metallic-dielectric nanocomposite media are advanced materials to be used as microcavity of photonic-crystal resonator. In the field of plasmonic resonance vigorous dispersion of optical properties of these materials is observed [3, 4]. This paper describes the case of ultrathin resonance structure as a monolayer of metal nanoparticles, plasmonic frequency of which coincides with defect mode frequency of PCS. 2. PCS material parameters and transfer matrixes In order to calculate reflectivity and transmission of plane-layered structure with embedded monolayer of nanoparticles we employ T-matrix technique. A special case is interface, optical qualities of which are determined by Fresnel reflection and transmission coefficients [5]. Since array of nanoparticles situated in the same plane interacts with electromagnetic wave like plane interface, it can be also treated as an interface with its own reflection and transmission coefficients. We assume that there are N interfaces in the layered medium, and they are formed by N-1 interfacial boundaries and a single layer of nanoparticles. A space between interfaces is packed by media with different refraction indexes ni  i  0... N  . Semi-infinite media are those that have n0 and nN refraction indexes. Let a harmonic wave is incident on a layered structure in z-direction. To describe its propagation in PCS we introduce the following notation for electric field components inside structure: Ei ( zi ) to the left of i number interface; Ei ( zi ) to the right of i number interface; E f for the propagating forward wave; Eb for the propagating backward wave. According to the introduced notation, complex amplitudes of counter-propagating waves on m interface in the layer with reflection index nm1 are equal to E f ( zm ) and Eb ( zm ) . At the same interface but in the layer with reflection index nm they are equal to E f ( zm ) and Eb ( zm ) . Relationship of these fields on m-interface (to the left and to the right of it) can be expressed as matrix equation:  E f ( zm )   E f ( zm )     = I m 1,m    , (1)  E b ( zm )   Eb ( zm )  1 1  rm ,m 1  I m 1,m = r , (2) tm 1,m  m 1,m tm 1,m tm ,m 1  rm 1,m rm ,m 1  3rd International conference “Information Technology and Nanotechnology 2017” 43 Computer Optics and Nanophotonics / I.A. Glukhov, S.G. Moiseev where ri , j , ti , j are complex reflection and transmission coefficients of the interface deviding media with refrection indexes ni and n j when the lightwave is incident from the medium with refraction index ni . In case of plane interface ri , j and ti , j are Fresnel coefficients [5]. Relationship of the fields on two interfaces of m and m  1 numbers confining homogeneous layer of m number is via transfer matrix Fˆm :  E f ( zm )  ˆ  E f ( zm )     = Fm    , (3)  E b ( zm )   Eb ( z m )   exp( i m ) 0  Fˆm =  , (4) 0 exp(i m )  where  m = knm Lm is phase thickness of the layer; k   / с is the wave number. Applying expressions (1) – (4) to the entire PCS we obtain relation for the amplitudes to the left of the first interface and to the right of the last (with number N) interface:  E f ( z1 )  ˆ  E f ( zN )     = G   , (5)  Eb ( z1 )   Eb ( zN )  Gˆ = I 0,1 F1 I1,2 F2 FN 1 I N 1, N . (6) Note that in semi-infinite medium with refraction index nN there exists only transmitted wave, therefore we shall assume Eb ( z N ) = 0 in (5). Reflectance and transmittance of PCS are calculated from the formulas 2 2 1 Gˆ T= , R = 21 . (7) Gˆ11 Gˆ11 To calculate field distribution in PCS we can use expressions similar to (5), in the left part of which column elements are substituted with amplitudes of local fields in the corresponding points. 3. Reflection and transmission spectra of nanoparticle monolayer For the analysis of the properties of PCS with embedded nanocomposite monolayer film of nanoparticles placed into dielectric matrix, it is necessary to know amplitude coefficients of monolayer reflection and refraction. Analytical calculation of these coefficients is not a trivial task; therefore we use numerical technique – FEM implemented in COMSOL Multiphysics software. For simplicity we consider the case of the ordered monolayer film in which nanoparticles are located at the sites of square lattice lying in the palne (xy). Taking into account structure symmetry, we took fourth of the structure’s unit cell (Fig. 1). Such domain contains fourth of nanoparticle and is of the size equal to ½ period of the structure along the direction of x and y coordinate axes. In order to obtain the entire monolayer it is necessary to apply reflection operations to the domain shown in Fig. 1 (to complete the unit to the full with a spherical shape particle), and then to apply transmission operations along x and y coordinate axes. The structure thus obtained will be 2D array of nanoparticles in (xy) plane. Boundary conditions are selected in such a way that the model fits the case of normal (to the plane of monolayer film) incidence of the light wave. Incident polarization is oriented parallel to one of the crystal axes of monolayer film. Fig. 1. Modeling domain in Comsol Multiphysics software. 3rd International conference “Information Technology and Nanotechnology 2017” 44 Computer Optics and Nanophotonics / I.A. Glukhov, S.G. Moiseev As a medium in which nanoparticle monolayer is weighted, we use material with dielectric constant  m  2.25 . To calculate dielectric constant of metal nanoparticles we use relation of Drude theory [5]:  p2  p ( ) =  0  , (8)  2  i where  p is plasmonic frequency,  0 is lattice contribution,  is relaxation parameter. For definiteness, as nanoparticle material we use silver, for which  p = 1.36 1016 c-1,  0 = 5 ,  = 3 1013 c-1. Modeling outcomes for optical properties of silver nanoparticle monolayer are shown in Fig. 2. It is obvious that in the domain of surface plasmonic resonance of nanoparticles (resonant wavelength falls within 435 nm) monolayer film reflection and transmission spectra are subject to strong changes. Amplitudes of the observed values in resonance region depend on the surface concentration of nanoparticles: with decrease in the average distance among nanoparticles resonance becomes more pronounced; the width of resonance increases and frequency shift of resonance towards short wavelength region is observed. Thus, spectral characteristics of the monolayer of metal nanoparticles depend on internal geometrical parameters which make it possible to control to a certain extent its influence on spectrum of PCS. 4. Analysis of the properties of photonic-crystal structure with monolayer of nanoparticles Let us consider PCS, in which between two dielectric reflectors there is a defect layer consisting of dielectric matrix and a monolayer film of nanoparticles. Transfer matrix of such structure can be expressed as: N = I 0,1 F1 I d 1, d Fd 1 Fs Fd 2 I d , d 1 FN 1 I N 1, N = M a DM b , (9) where Fs is transfer matrix of nanoparticle monolayer; Fd 1 and Fd 2 are transfer matrixes of the layers that edge monolayer of nanoparticles; M a and M b are transfer matrixes that describe dielectric reflectors containing a and b binary layers respectively. Binary layers of dielectric reflectors consist of two layers of isotropic dielectric material of real transmissivities  j and thickness L j ( j = 1, 2). For modeling optical properties of the structure we assumed 1 = 6.25 and  2 = 2.25. Thickness of layers of the structure meets the requirements L1,2 = 0 /4 1,2 ; thickness of defect layer is equal to Ld = 0 /  d , where 0 – wavelength in vacuum calculated for the central frequency of the photonic band gap. Presence of defect layer in PCS leads to the occurrence in the photonic band gap of narrow spectral transmission band with the peak value of transmission index which is close to 1. 1 0.8 0.6 0.4 0.2 0 360 400 440 480 520 wavelength, nm Fig. 2. Reflection, transmission and absorption spectra for monolayer of silver nanoparticles. Period of structure is 10 nm. Figures 3 and 4 show field distribution (square amplitude of dielectric field intensity) and transmission and reflection spectra of PCS of M 5 DM 5 -type depending on location of the monolayer of plasmonic nanoparticles. Location of the monolayer is given by relations Ld1 = Ld / 2   , Ld 2 = Ld / 2   , where Ld is thickness of the central layer, which is divided by a monolayer of nanoparticles into two domains of the thickness Ld 1 and Ld 2 respectively. The domain of the thickness Ld 1 is located from incidence of external electromagnetic wave. Thus,  is a value of monolayer film displacement from the center of PCS. Parameters of dielectric layers of PCS correlate such that in the center of PCS, i. e. for   0 , field amplitude is close to zero (standing-wave nod), and in displacement by   71 nm the amplitude reaches a maximum value (antinode of standing wave). The given relations indicate that if monolayer film of nanoparticles is located in the center (   0 ), field distribution and transmission and reflection spectra are practically comparable with the case of monolayer-free PCS. Such peculiarity is explained practically by total absence of electrodynamic interaction between monolayer film and lightwave in standing-wave node. Monolayer film is really in the domain of electromagnetic shadow. If nanoparticles are in the domain of strong field (   71 nm), when intensity 3rd International conference “Information Technology and Nanotechnology 2017” 45 Computer Optics and Nanophotonics / I.A. Glukhov, S.G. Moiseev of electromagnetic scattering by monolayer is the highest, defect mode transformation is observed. Transformation is displayed as decrease in its amplitude and splitting of spectrum curves. Fig. 3. The distribution of electromagnetic field amplitude (in relative units) throughout PCS as a function of the value of displacement of monolayer of nanoparticles from the center of the structure. R, T 1 (а) 0.8 0.6 0.4 0.2 0 400 500 600 wavelength, nm R, T 1 (b) 0.8 0.6 0.4 0.2 0 400 500 600 wavelength, nm Fig. 4. Ttransmission (dashed line) and reflection (continuous line) spectra of PCS when monolayer displacement is (а)   71 nm, (b)   0 . 5. Conclusions We have demonstrated that control of transmission and reflection coefficients corresponding to defect mode in photonic band gap of PCS is capable owing to the use of nanoparticle monolayer with plasmonic resonance. It is shown that defect mode amplitude is heavily dependent on the location of monolayer. Dependence of spectral characteristics of the layered structure on the location of plasmonic monolayer is attributed to the inhomogeneity of electromagnetic field distribution in the optical microcavity placed between distributed Bragg reflectors. Acknowledgements This work was supported by the Ministry of Education and Science of the Russian Federation (State Contracts Nos. 3.5698.2017/P220 and 3.8388.2017/ITR) and the Russian Foundation for Basic Research (Projects Nos. 15-07-08111 and 17-02-01382). 3rd International conference “Information Technology and Nanotechnology 2017” 46 Computer Optics and Nanophotonics / I.A. Glukhov, S.G. Moiseev References [1] Vorobev LE, Ivchenko EL, Firsov DA, Shalygin VA. Optical properties of nanostructures. Saint Petersburg: “Nauka” Publisher 2001. [2] Gaponenko SV, Rosanov NN, Ivchenko EL. Optics of nanostructures. Saint Petersburg: “Nedra” 2005. [3] Moiseev SG, Ostatochnikov VA. Defect modes of one-dimensional photonic-crystal structure with a resonance nanocomposite layer. Quantum Electron 2016; 46(8): 743–748. [4] Vetrov SY, Avdeeva AY, Timofeev IV. Spectral properties of a one-dimensional photonic crystal with a resonant defect nanocomposite layer. JETP 2011; 113(5): 755–761. [5] Born M, Wolf E. Principles of Optics. Cambridge: Cambridge University Press 1999. 3rd International conference “Information Technology and Nanotechnology 2017” 47