Arxiv:1911.02884V1 [Physics.Optics] 7 Nov 2019
Total Page:16
File Type:pdf, Size:1020Kb
Boundary conditions for the effective-medium description of subwavelength multilayered structures Maxim A. Gorlach1, ∗ and Mikhail Lapine2 1ITMO University, Saint Petersburg 197101, Russia 2School of Mathematical and Physical Sciences, University of Technology Sydney, NSW 2007, Australia Nanostructures with one-dimensional periodicity, such as multilayered structures, are currently in the focus of active research in the context of hyperbolic metamaterials and photonic topological structures. An efficient way to describe the materials with subwavelength periodicity is based on the concept of effective material parameters, which can be rigorously derived incorporating both local and nonlocal responses. However, to provide any predictions relevant for applications, effective material parameters have to be supplemented by appropriate boundary conditions. In this work, we provide a comprehensive treatment of spatially dispersive bulk properties of multilayered metamaterials as well as derive boundary conditions for the averaged fields. We demonstrate that local bianisotropic model does not capture all the features related to second-order nonlocal effects in the bulk of metamaterial. As we prove, while the bulk response of multilayers does not depend on the unit cell choice, effective boundary conditions are strongly sensitive to the sequence of layers and multilayer termination. The developed theory provides a clear interpretation of the recent experiments on the reflectance of all-dielectric deeply subwavelength multilayers suggesting further avenues to experimentally probe electromagnetic nonlocality in metamaterials. I. INTRODUCTION However, it is not just bulk nonlocality that deter- mines the properties of multilayers, since boundary ef- fects should also be taken into account [13{15]. There- Electromagnetic properties of multilayered structures fore, to make the nonlocal description self-consistent, one and their effective material parameters have attracted re- has to supplement bulk effective permittivity by appro- search interest since early days [1]. Current fabrication priate boundary conditions. Since homogenized descrip- techniques based on sputtering or atomic layer deposi- tion includes averaged fields, boundary conditions should tion allow one to fabricate multilayers with deeply sub- also be formulated for the averaged fields, and therefore wavelength thickness and low roughness truly approach- it is not obvious a priori that the continuity of tangential ing the regime of metamaterial [2 and 3]. Nevertheless, components of electric and magnetic fields will still hold. a series of studies has warned against application of the Furthermore, previous studies contain a clear indica- standard frequency-dependent (i.e. local) effective mate- tion that the effective boundary conditions should be rial parameters to the multilayered metamaterials even in modified. For instance, it has recently been suggested the subwavelength regime [4{7]. In particular, multilay- theoretically [16] and verified experimentally [2] that the ers can give rise to tri-refringence phenomenon [6], while reflectance of all-dielectric multilayered structure with the standard techniques of effective parameters retrieval the layers of subwavelength thickness deviates from the can yield unphysical results for multilayers [8]. predictions of the local effective medium model also de- Such peculiar behavior has recently been attributed to pending on the structure termination. Shortly after- the nonlocal (or spatially dispersive) electromagnetic re- wards, quite a few proposals have been put forward on sponse of multilayers, which manifests itself through the how to interpret this peculiar feature [14, 15, 17, and dependence of polarization on electric field in the neigh- 18]. Most importantly, this feature can not be explained boring regions of space. Such electromagnetic nonlocality by bulk nonlocality alone, and therefore one has to work is described via the dependence of effective permittivity out the correct form of boundary conditions in order to tensor" ^(!; k) on wave vector k which is considered as capture this phenomenon. a variable independent of ! in order to capture linear The rest of the article is organized as follows. In Sec- response of the structure to the arbitrary excitation [9]. tion II, we provide a general analysis of the bulk non- arXiv:1911.02884v1 [physics.optics] 7 Nov 2019 To evaluate effective nonlocal permittivity tensor, it has local response of multilayered structure with arbitrary been proposed to use current-driven homogenization ap- number of layers in the unit cell and derive an explicit proach based on the analysis of the structure response to expression for the nonlocal permittivity tensor" ^(!; k) in- the external distributed currents [10 and 11]. The effec- corporating spatial dispersion corrections up to the sec- tive susceptibility is then found as a matrix which relates ond order. Using this explicit expression and the link polarization averaged over the unit cell to the averaged between local and nonlocal descriptions [10 and 19], we electric field. Recently, this strategy has been applied analyze in Sec. III, whether it is possible to formally de- to the multilayered metamaterial composed of isotropic scribe all second-order spatial dispersion effects in the layers of two types, and an explicit though cumbersome bulk of multilayered structure in terms of local permit- expression for the effective permittivity tensor has been tivity and permeability tensors. The answer appears to derived for this particular case [7 and 12]. be negative. In Section IV we derive the effective bound- 2 ary conditions for multilayered structure incorporating spatial dispersion corrections up to the second order. As 2 3 we show, the form of boundary conditions appears to be 1 iky y ikx x X i(kx+nb) x sensitive to the termination of a multilayered structure, D(r) = e 4D0 e + Dn e 5 ; which provides an elegant interpretation of the depen- n6=0;n=−∞ dence of reflectance on the sequence of layers. To illus- (3) trate our results, in Sec. V we perform calculations for where b = 2π=a is the period of reciprocal lattice, and the a specific metamaterial with two layers in the unit cell. dependence of the fields on x coordinate is determined Finally, in Sec. VI we discuss the obtained results con- by the Bloch theorem. The amplitudes of Fourier coef- cluding with the outlook for future applications. ficients En and Dn should be chosen in order to satisfy inhomogeneous wave equation with external current: r(divE) − ∆ E = q2 D + J ; (4) II. DERIVATION OF THE BULK NONLOCAL RESPONSE where q = !=c, J(r) ≡ 4π i q j(r)=c and CGS system of units is employed. In this section, we provide a general derivation of the Besides that, using the structure periodicity, we ex- bulk response of multilayered structure assuming peri- pand permittivity "(x) and inverse permittivity ζ(x) ≡ odic permittivity modulation along x axis, "(x), with 1="(x) into the Fourier series: the period equal to a as illustrated in Fig. 1. Analo- 1 X gously to Refs. [7 and 12], we use nonlocal homogeniza- "(x) = " + " einbx ; (5) tion approach assuming excitation of the structure by the 0 n n6=0;n=−∞ external distributed currents 1 X ζ(x) = ζ + ζ einbx ; (6) ikx x+iky y 0 n j(r) = j0 e (1) n6=0;n=−∞ oscillating with the frequency !. For simplicity, we as- where the coefficients "n and ζn are related to each other sume that ! a=c 1 and ka 1, keeping the terms up due to the fact that "(x) ζ(x) ≡ 1. We aim to calculate to the second order in ! a=c or ka. To obtain the effec- the effective permittivity tensor of a multilayered struc- tive permittivity expanded in terms of small parameter, ture defined as we adopt the technique similar to Refs. [20 and 21] and expand the fields into Floquet harmonics as follows: D0 =" ^(!; k) E0 : (7) 2 1 3 Note that the amplitudes E0 and D0 precisely corre- iky y ikx x X i(kx+nb) x spond to the averaged fields used in metamaterial ho- E(r) = e 4E0 e + En e 5 ; n6=0;n=−∞ mogenization procedure [10], while the obtained" ^(!; k) (2) corresponds to effective permittivity definition in any pe- riodic medium [9]. It is also assumed that the external current J has only zeroth Floquet harmonic. TE TM Combining the expansions Eqs. (2), (3) with Eq. (4), we get a linear system of equations for En and Dn coeffi- E H cients. The entire system splits into two independent sets of equations which correspond to TE and TM-polarized H E waves propagating in the structure. We analyze TE case k k below, whereas more involved case of TM polarization is z y examined in Appendix A. In case of TE polarization, wave equation (4) and the a material equation D(r) = "(x) E(r) yield: 2 2 2 kx + ky E0z = q D0z + Jz ; (8) ε()x 2 2 2 (kx + nb) + ky Enz = q Dnz (n 6= 0) ; (9) x X D0z = "0 E0z + "−n Enz ; (10) FIG. 1. The scheme of multilayered structure with period a n6=0 X along x axis and arbitrary permittivity profile "(x) with two Dnz = "n E0z + "−m E(n+m)z ; (n 6= 0) (11) possible polarizations of plane wave sketched. For simplicity, m6=−n we assume that the wave vector k defined by the external distributed currents exciting the structure lies in Oxy plane. where the summation is performed from −∞ to 1. As a first step, we examine the quasistatic case when q = 3 k = 0. Inspecting Eq. (9), we recover that Enz = 0. nonzero off-diagonal components "xy and "yx. Their sym- Then Eq. (10) yields D0z = "0 E0z and the zz component metric part ("xy + "yx)=2 appears to be proportional to of effective permittivity is simply equal to the average the product kx ky, and it is this term which has been in- permittivity "0. terpreted in Ref. [7] as the rotation of the optical axis Hence, in the limit qa 1 and ka 1 we can consider of metamaterial induced by spatial dispersion.