Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 473079, 5 pages
Research Article

Impedance Operator Description of a Metasurface with Electric and Magnetic Dipoles

Université de Montpellier, Laboratoire Charles Coulomb UMR CNRS-UM 5221, 34095 Montpellier, France

Received 5 August 2015; Accepted 4 November 2015

Academic Editor: Ivan D. Rukhlenko

Copyright © 2015 Didier Felbacq. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


A metasurface made of a collection of nanoresonators characterized by an electric dipole and a magnetic dipole was studied in the regime where the wavelength is large with respect to the size of the resonators. An effective description in terms of an impedance operator was derived.

1. Introduction

Metasurfaces are the bidimensional analogue of metamaterials [1]. They are made of basic cells containing resonant elements disposed periodically or not on a surface. The cells can contain one or several elements (nanoparticles, nanoantennas, etc.) whose aim is to produce a dephasing of some incident field, whose wavelength is large with respect to the size of the cells. The collection of cells is then able to produce a collective effect resulting in refracting and transmitting properties far from that of a homogeneous material [2]. For instance, it was shown in [3] that they could be used to design flat achromatic lenses. Metasurfaces are thus a wide generalization of subwavelength diffraction gratings, where the basic cell can be quite complicated. The periodic structure of gratings can be relaxed so as to obtain interesting effect on the polarization of the incident field [4]. In this context, we study the field diffracted by a periodic set of linear nanoresonators, electromagnetically characterized by their scattering matrix . We are interested in the regime where the wavelength is much larger than the size of the nanoresonators. We proceed to an asymptotic analysis related to homogenization theory [5, 6]. The field diffracted by the structure is derived and it is shown that it is characterized by an impedance operator. Our results extend to electric and magnetic dipoles results in [7] where only electric dipoles were considered.

2. Setting of the Problem

The structure under study is made of an infinite number of resonators invariant along , periodically disposed at points , where is the period and . Each scatterer at position is characterized in the frequency domain by a scattering matrix as well as by electromagnetic parameters and . We assume that the wavelength in vacuum is much larger than the size of the resonators, which are assumed to be contained in a cylinder of diameter . We therefore consider a linearly polarized incident field , where and , is an angle of incidence, and . According to the polarization and can be either the relative permittivity or relative permeability.

The stationary scattering problem considered therefore reads the following: find a field (i.e., the space of locally square integrable functions with a locally square integrable weak derivative) satisfying in the sense of Schwartz distribution in and such that satisfies the outgoing wave conditions:(i)For .(ii)For .

The functions and are equal, respectively, to and inside the scatterers and to 1 outside. The following holds; see [8].

Proposition 1. Apart possibly from a discrete set of wavenumbers , the field exists and is unique.

Our point is to provide a simplified expression of the scattering problem by means of an impedance operator.

Let denote the translation along of amplitude ; that is, . For later purpose, we define , , , and , for . For , we impose . In the following, we denote .

Let and let us define for the field of Hilbert spaces . It obviously holds the following.

Lemma 2. The commutator of and vanishes: .

Applying Floquet-Bloch analysis [9], we obtain the following.

Proposition 3. The operator has a direct integral decomposition where with domain .

3. Multiple Scattering Approach

The incident field has the expansion [10] . For one scatterer alone, the incident field gives rise to a field where . (resp. ) is a Bessel (resp., Hankel) function. For the infinite set of scatterers, this gives a diffracted field that reads Multiple scattering theory [10] allows to write that for where and . The matrix is given by Here Proposition 3 was used through the introduction of a Bloch phase . In this expression ; that is,where (note that ). The following series [11] indexed by appear:And the entries of the matrix are .

In the regime where , the cylinder can be described by a scattering matrix (this corresponds to an electric dipole and a magnetic dipole) and the field by coefficients , , and [10]. Therefore, only series are involved: , , and . It holds thatIn the extreme limit () where the scatterers are very small as compared to the wavelength and the period, the scattering matrix reduces to a scalar matrix : the scatterers are thus dipoles with a dipole moment along and the only involved series is ; this situation was addressed in [7]. The multiple scattering relation (2) then becomeswhere the series can be written [7, 12, 13]:For the more general case of an electric dipole and a magnetic dipole the following asymptotic expressions hold in the limit [14].

Proposition 4. Consider the following:

4. Scattering Properties of the Meta Surface

Define the electric moment and the magnetic moment. We write , , and , .

We can now state the following.

Theorem 5. The total field has the following expression:where

Proof. We start with the following relation, obtained from Poisson formula: Upon applying the operator , using the fact that the series on the right-hand side is normally convergent (thanks to the term ) and using the relation , we obtain where . Therefore we get The result follows after some simple algebra.

A simple, but interesting corollary is as follows.

Corollary 6. As ,

We are not a priori in the homogenization regime where and hence there can be several reflected and transmitted orders. In expression (10), the propagative waves (i.e., the diffractive orders of the grating) correspond to the ’s that are real. They are labelled by the finite set . The evanescent waves are labelled by the infinite set .

5. Impedance Operator Formulation

Our point is now to replace the set of nanoresonators by a metasurface which is simply the line . This requires to specify the boundary conditions there in terms of an impedance operator.

So, consider the continuation of the field obtained by making in (10). The continued field is still denoted as . It is a singular distribution. To handle this situation, let us introduce the following fields of Sobolev spaces and the dual spaces Let be the pseudodifferential operator defined by , whereIt is straightforward to show the following.

Proposition 7. is continuous and invertible from to , for .

The inverse of is the admittance operator defined by .

Let us denote as an element of , representing the discontinuity of the field and its derivative through . The traces of the field and its derivative are By definition, it holds that . The Calderòn projectors [15] and are defined by and . The preceding shows the following.

Proposition 8. Consider the following:

Obviously, it holds that , , , and , as it should.

The transmission conditions on can be written as follows: This suggests to define the following pseudodifferential operators, acting on : Both and are bounded with respect to (see (15)); hence see the following.

Proposition 9. The pseudodifferential operators and are isomorphism of .

These conditions can be rewritten conveniently in the operator form.

Theorem 10. The traces of the field diffracted by the metasurface under the illumination of a plane wave satisfy the impedance conditions:where

The operator is the transfer matrix of the meta surface. The discontinuity of the (effective) field at is due to the existence of a magnetic dipole moment. In the homogenization limit of very large wavelengths, that is, larger than the wavelengths corresponding to magnetic resonances and larger than twice the period, there are only one transmitted and one reflected wave; the evanescent waves can be discarded and the magnetic resonances have no effect; we then have the following.

Proposition 11. For a wavelength and larger than the largest magnetic resonance, the propagative part of the field is given by where and . The transfer matrix becomes

The form found for the transfer matrix matches that obtained in another context in [12].

6. Conclusion

The scattering of a plane wave by a grating of nanorods was described in the framework of metasurfaces by exhibiting an impedance condition. This takes into account both the electric and the magnetic dipoles characterizing each nanorod. This study can be generalized to higher multipoles but also to nonperiodic, for instance, quasiperiodic, structures [16] and to elementary scatterers deposited on an arbitrary smooth surface as well. A similar approach was used in [17] to study the coupling of a quantum emitter with the modes supported by the metasurface; see also [18]. The proposed formalism can be used in this context to analyze the role of resonances.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.


  1. C. L. Holloway, E. F. Kuester, J. A. Gordon, J. O'Hara, J. Booth, and D. R. Smith, “An overview of the theory and applications of metasurfaces: the two-dimensional equivalents of metamaterials,” IEEE Antennas and Propagation Magazine, vol. 54, no. 2, pp. 10–35, 2012. View at Publisher · View at Google Scholar · View at Scopus
  2. N. Yu and F. Capasso, “Flat optics with designer metasurfaces,” Nature Materials, vol. 13, no. 2, pp. 139–150, 2014. View at Publisher · View at Google Scholar · View at Scopus
  3. M. Khorasaninejad, F. Aieta, P. Kanhaiya et al., “Achromatic metasurface lens at telecommunication wavelengths,” Nano Letters, vol. 15, no. 8, pp. 5358–5362, 2015. View at Publisher · View at Google Scholar
  4. X. Ma, M. Pu, X. Li et al., “A planar chiral meta-surface for optical vortex generation and focusing,” Scientific Reports, vol. 5, Article ID 10365, 2015. View at Publisher · View at Google Scholar
  5. L. Tartar, The General Theory of Homogenization, Springer, New York, NY, USA, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  6. D. Felbacq and G. Bouchitté, “Homogenization of a set of parallel fibres,” Waves in Random Media, vol. 7, no. 2, pp. 245–256, 1997. View at Publisher · View at Google Scholar · View at MathSciNet
  7. D. Felbacq, “Layer homogenization of a 2D periodic array of scatterers,” Photonics and Nanostructures—Fundamentals and Applications, vol. 11, no. 4, pp. 436–441, 2013. View at Publisher · View at Google Scholar · View at Scopus
  8. A.-S. Bonnet-Bendhia and F. Starling, “Guided waves by electromagnetic gratings and non-uniqueness examples for the diffraction problem,” Mathematical Methods in the Applied Sciences, vol. 17, no. 5, pp. 305–338, 1994. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. M. Reed and B. Simon, Methods of Modern Mathematical Physics, vol. 4 of Analysis of Operators, Academic Press, New York, NY, USA, 1978. View at MathSciNet
  10. D. Felbacq, G. Tayeb, and D. Maystre, “Scattering by a random set of parallel cylinders,” Journal of the Optical Society of America A, vol. 11, no. 9, pp. 2526–2538, 1994. View at Publisher · View at Google Scholar · View at MathSciNet
  11. C. M. Linton, “Lattice sums for the Helmholtz equation,” SIAM Review, vol. 52, no. 4, pp. 630–674, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  12. F. Zolla, D. Felbacq, and G. Bouchitté, “Bloch vector dependence of the plasma frequency in metallic photonic crystals,” Physical Review E, vol. 74, no. 5, Article ID 056612, 2006. View at Publisher · View at Google Scholar · View at Scopus
  13. P. A. Belov, S. A. Tretyakov, and A. J. Viitanen, “Dispersion and reflection properties of artificial media formed by regular lattices of ideally conducting wires,” Journal of Electromagnetic Waves and Applications, vol. 16, no. 8, pp. 1153–1170, 2002. View at Publisher · View at Google Scholar · View at Scopus
  14. A. I. Cabuz, Electromagnetic metamaterials—from photonic crystals to negative index composites [Ph.D. thesis], Université Montpellier II, 2007,
  15. M. Cessenat, Mathematical Methods in Electromagnetism, World Scientific, Singapore, 1996. View at Publisher · View at Google Scholar · View at MathSciNet
  16. Z. V. Vardeny, A. Nahata, and A. Agrawal, “Optics of photonic quasicrystals,” Nature Photonics, vol. 7, no. 3, pp. 177–184, 2013. View at Publisher · View at Google Scholar · View at Scopus
  17. D. Felbacq, “Weak and strong coupling of a quantum emitter with a meta-surface,” Superlattices and Microstructures, vol. 78, pp. 79–87, 2015. View at Publisher · View at Google Scholar · View at Scopus
  18. A. Cazé, R. Pierrat, and R. Carminati, “Strong coupling to two-dimensional Anderson localized modes,” Physical Review Letters, vol. 111, no. 5, Article ID 053901, 2013. View at Publisher · View at Google Scholar · View at Scopus