Issue 
EPJ Appl. Metamat.
Volume 7, 2020
Frontiers in microwave, photonic, and mechanical metamaterials



Article Number  10  
Number of page(s)  9  
DOI  https://doi.org/10.1051/epjam/2020009  
Published online  19 January 2021 
https://doi.org/10.1051/epjam/2020009
Research article
Inverse design of EBG waveguides through scattering matrices
^{1}
CNRIREA, National Research Council of Italy – Institute for Electromagnetic Sensing of Environment, Via Diocleziano 328, Napoli 80124, Italy
^{2}
DIIES, Università Mediterranea of Reggio Calabria, Via Graziella, Feo di Vito, Reggio Calabria 89122, Italy
^{*} email: palmeri.r@irea.cnr.it
Received:
1
November
2020
Accepted:
1
December
2020
Published online: 19 January 2021
Inverse design procedures aim at determining optimal parameters for a given device in order to satisfy assigned specifications. In this contribution, the design of optimal EBG waveguides through inverse problems tools is addressed. In particular, an inversion tool based on the so called ‘scattering matrices’ is proposed and assessed to optimize the guiding effect for straight and bent waveguides.
Key words: Artificial materials / EBG / inverse design / inverse problems / optimization problem / scattering matrix
© R. Palmeri and T. Isernia, published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The evergrowing interest in artificial materials, including metamaterials/metasurfaces and photonic crystals, is related to their very intriguing properties and possibilities allowing very numerous and even exotic applications in all ranges of the electromagnetic spectrum including optics [1], microwaves [2] and millimeterwaves frequencies [3].
In this contest, devices based on electromagnetic bandgap (EBG) structures result very attractive and are exploited for a wide class of applications, ranging from radiation [4] and propagation [5] applications, to particle accelerators [6]. Amongst them, in this contribution we focus on EBG waveguides.
EBG structures [7] are periodic arrangement of inclusions usually adopted to improve the performance of devices or to realize adhoc functionalities for a given device. In particular, EBG structures have the fascinating property of suppressing field propagation within a frequency band (depending on the structure itself) while allowing it in other bands. If a defect is introduced in the periodic arrangement by removing some inclusions, such a defect will support an electromagnetic mode and hence it is possible to generate localized or guided fields [7].
Despite the very large number of analysis tools for this kind of structures, both at ‘macroscopic’ (like homogenization techniques) and ‘microscopic’ (as finiteelementmethods) level, as well as at a ‘mesoscale’ level [8], formalized synthesis and design tools are quite limited in literature. Basic approaches are founded on trialanderror cyclic processes, making extensive use of numerical simulations that are occasionally combined with geneticalgorithmsbased optimizations [9]. Obviously, because of the huge number of involved parameters and the corresponding high number of trials, this is a farfrom effective way to proceed. More rigorous design approaches are those based on inverse design, such as the inverse homogenization problems [10,11] which are mainly devoted to synthesizing effective electromagnetic properties of unit cells. Unfortunately, the arising dielectric profiles exhibit a complex distribution which is difficult to fabricate [10], so that fabricationoriented tools are needed. In this respect, topology optimization [12] based design tools have been widely adopted during the years. As the name suggest, this kind of techniques aims at modifying the topology of the device at hand, that is, the geometrical and e.m. characteristics of the object being designed, in such a way some assigned specifications are satisfied. Developed design tools based on topology optimization [13–17] or the like [18,19] are able to realize binary devices wherein just two materials (one of which is usually air) are involved. Other types of design tools based on group theory [9] and the emerging concept of deep neural network [20] have been also proposed.
Recently, the design of artificial materials through an inverse design approach has been introduced within the framework of inverse scattering problems for antennas [21] and cloaking [22,23] applications, showing very promising performance compared to homogenizationbased approaches, but also evidencing limitations as far as the computational burden is concerned. In fact, the adopted tools in [21] require very small mesh elements to correctly represent small scatterers and hence the field they scatter.
As an alternative to above techniques and with the aim of overcoming relates issues, we have recently proposed, for a completely different problem concerned with antennas (see below), a novel inverse design approach based on inverse scattering coadjuvated by the use of the so called Scattering Matrix Method (SMM) [24]. The SMM is a very accurate and fast analysis tool for structures made of multiple dielectric or metallic scatterers. In fact, it allows to exploit analytical (or offline numerical) results related to the ‘response’ of each single inclusion to have an accurate modelling of the overall device. Notably, it operates somehow at a ‘mesoscale’ level, as the modelling is not based on the microscopic level of numerical simulations acting on a dense grid, and avoids the macroscopic but possibly approximated way of thinking based on homogenization techniques. Notably, the basic SMM method can be eventually improved (see [8] for more details) by introducing aggregations of elementary scatterers and computing the scattering matrix of the larger scatterer which is obtained.
As we are going to summarize in the next section, the scattering matrix, through proper expansions of the incident field (acting as input) and scattered fields (acting as output), allows to take into account in a very effective fashion the ‘response’ of each single scatter at hand. Then, one is left with a system of equations governing (just) the interactions amongst the different scatterers, while interactions internal to the single scatterers constituting the device are taken into account through the corresponding scattering matrix. Of course, the number of equations and the dimension and structure of the scattering matrices depend on the number of scatterers as well as on their geometrical and electromagnetic characteristics.
With the aim of pursuing the design of EBG devices, the SMM formulation can be conveniently turned into an inverse problem. More in detail, the role of data and unknowns is inverted so that the scattering matrices of each inclusion can be considered the new unknowns of the problem, while expansion coefficients of the scattered fields must be determined in such a way to obey assigned field specifications, depending on the device at hand.
Very recently, we have proposed our Inverse SMM (ISMM) tool in the framework of innovative antennas design. In fact, we have successfully applied it in [25] for the same design problem considered in [21]. Note that in the antenna applications above, the specifications were derived starting from an optimal synthesis of the coefficients of a representation of the scattered and total field, and the pursuing of such a field is then the goal of the inverse scattering step, so that such a field substitute the field measurements which are the usual data of inverse scattering for reconstruction problems.
In this contribution, we consider a completely different problem with a completely different kind of design goal. As a matter of fact, EBG guiding structures are considered rather than antennas, and the actual goal is the optimization of the Voltage Standing Wave Ratio (VSWR) along the structure in such a way to avoid undesired reflections in case of bends, rather than the realization of a given field.
Notwithstanding the completely different problem and the different kind of goals, it is shown that scattering matrices, as well as the kind of (state and data) equation which are the core of inverse scattering procedure can still be used.
As well known, the electromagnetic response of a class of EBG devices can be understood by analysing the bands diagram which depends in turn on the kind of lattice, polarization, and the used defects (if any). On the other side, such a ‘basic’ way of thinking does not exploit a number of degrees of freedom which can be gained by changing (for example) the permittivity or dimensions of (some of) the inclusions. For example, which is the basic problem we deal with in the following as a first demonstration of potentialities, one could use a few degrees of freedom regarding a limited number of inclusions to improve the behaviour of an EBG waveguide. In fact, by starting from a basic structure, inverse design procedures are exploited to modify (part of) the lattice for an improving of the VSWR (i.e., saying it in other words, for avoiding waves coming back along the defect). A similar problem has been dealt with in [13–17]. In particular, we apply the ISMM to improve guiding performance of straight and bent EBG waveguides. Our procedure resembles the one in [26], wherein the multiple scattering theory (MST) is adopted. Both approaches are based on a harmonic expansion for fields and their formulation is general, so they can be formalized for different kind of inclusions; the main difference is that, opposite to the SMM, the MST exploits the Tmatrix to link the scattered field to the electromagnetic properties of the scatterer. As far as the design procedure is concerned, an important difference can be highlighted amongst our approach and the one based on the MST. In fact, in [26] a global optimization is performed for the solution of the inverse problem. As a consequence, to limit the computational burden, several apriori information must be used. In particular, in [26] the authors not only set arrangement and dimension of scatterers, but also their electromagnetic parameters, the remaining degrees of freedom being just the presence or not of the scatterers themselves. In our approach instead the inverse problem is solved through a conjugate gradient based algorithm that is able to handle both geometric and electromagnetic features, assumed that a favourable starting point can be provided for the optimization step.
The remainder of the paper is structured as follows. In Section 2 the SMM is briefly recalled while in Section 3 its inverse formulation is introduced. Section 4 is devoted to present and assess the procedure for designing optimized EBG waveguides. Conclusions follow.
For the sake of simplicity, 2D scalar fields are considered and the time harmonic factor e ^{ jωt } is assumed and dropped.
2 Basics of the Scattering Matrix Method (SMM)
Let consider the geometry depicted in Figure 1 wherein the ℓth cylinder with cross section C _{ℓ} is located at (r _{ℓ}, θ _{ℓ}) with respect to a global reference system Oxy. Let also E _{ inc } be a given incident field impinging on the set of scatterers. By exploiting the cylindrical harmonics expansion of fields [27], the field scattered by the set of inclusions in the regions of space outside the single scatterers reads as:(1)wherein N is the total number of scatterers, k _{ b } is the wavenumber of the background medium and is the mth order Hankel function of second kind. Notably, representation (1) is expressed in the reference system of the ℓth scatterer (see Fig. 1).
As soon as the observation point is a little bit outside an inclusion C _{ℓ}, the corresponding internal summation can be safely truncated^{1} to a finite number, say M _{ℓ}, which also depends on the single inclusion maximum dimensions [28].
Accordingly, the total field can be represented as:(2)
From equation (2) it is possible to note three contributions: the first one is the primary incident field; the second addendum is a summation of scattered fields from all inclusions but the ℓth one, that accordingly acts as a secondary incident fields on it; finally, the last term simply represents the field scattered by the ℓth cylinder as defined in equation (1).
In equation (2), the fields representations are still expressed in the reference system linked to the pertaining cylinder. In order to have a field representation in a unique global reference system (e.g., the ℓth one), the Graf's formula [29] for the Hankel function can be used. We finally obtain what follows (see [24] for more details):(3)
a _{ℓ,m } coefficients being the following compact expression:(4)wherein Q _{ℓ,m } encodes terms arising from the harmonic expansion of the primary source, while T _{ℓ,i,m,q } takes into account mutual coupling effects amongst scatterers [24], that is, it is related to the secondary incident field. For the sake of completeness, analytic expressions of T _{ℓ,i,m,q } and Q _{ℓ,m } taken from [24] are given in Appendix.
If b _{ℓ} is the vector encoding the field representation coefficients while a _{ℓ} is the vector encoding the source representation coefficients, and according to the fact that a finite number of parameters can be used to represent both fields and sources [28], for any fixed scatterer one can conveniently design a finite dimensional ( M _{ℓ} × M _{ℓ}) scattering matrix S _{ℓ} encoding all relevant information of the scatterer, such that:(5)
Hence, after multiplying equation (4) by the scattering matrix, we finally achieve the linear system in the unknown b _{ℓ}:(6)
that formally describe the scattering model at hand.
Once equation (6) is solved, the scattered field in the region of space outside scatterers can be finally evaluated through equation (1) [24].
Fig. 1
Geometry sketch for the problem at hand. 
3 The InverseSMM (ISMM)
The aim of an inverse design problem is to infer optimal properties on certain parameters, within a set of constraints, to accomplish assigned specifications at best. In the present case, the parameters to be optimized are the geometric and electromagnetic features of inclusions in such a way the scattering from the optimized set is able to satisfy some given specification on the scattered or total field.
If a nominal assigned field to be pursued is given, in system (6) both the scattering matrices S _{ℓ} and the coefficients b _{ℓ} are unknowns, and these latter must be determined such to fit the given assigned field as well, that is, such that:(7)
Accordingly, in the ISMM design problem, S _{ℓ} and b _{ℓ} must be optimally retrieved to contemporarily satisfy the model (6) and the data constraints (7). Notably, such an inverse problem is not trivial to solve because of its nonlinearity that could cause the optimization problem (based on a gradient scheme) to be trapped into local minima. However, as long as a good matching is gathered with both the model and the data, suboptimum solutions can be accepted.
A further issue to be taken into account concerns the existence of a solution for the problem, especially in terms of a final actual device. More in detail, even if scattering matrices S _{ℓ} can be found as solution of the inverse design problem, the existence of corresponding physical reliable scatterers is not ensured. For this reason, apriori information on the scatterers are required to gather a feasible solution of the problem. This allow, for instance, to look for diagonal matrices if cylindrical scatterers with circular symmetric crosssection are dealt with. Also, one can look for a single S matrix if the scatterers are identical. Moreover, we can directly search for the parameters actually determining the entries of the matrices if the electromagnetic features of the inclusions are supposed to be of a given kind, namely if we are dealing, for instance, with homogeneous dielectric or metallic inclusions, single layer or stratified inclusions, and so on. In this case, in fact, analytic expressions for the scattering coefficients are available [25].
At the end of the day, thanks to the fact one is reducing as much as possible the number of unknowns of the problem and by taking advantage from convenient starting points derived from a general physical understanding or homogenization based suggestions, the minimization procedure will ensure actual solutions of the design problem, without resorting to global optimization algorithms, with the inherent advantages in terms of computational burden.
4 An ‘ad hoc’ procedure for optimal EBG waveguides design
The SMM and ISMM presented above are general formulation for the scattering model dealt with, namely they still apply whatever the arrangement of the scatterers, their electromagnetic (e.g., dielectric, metallic, stratified) and geometric (cylinders with circular or elliptic crosssection) features, whatever the background medium and the primary source^{2}. Thereinafter, in order to deal with EBG devices, we consider circular lossless dielectric scatterers with radius ρ _{ℓ} and constant relative permittivity ϵ _{ℓ}, embedded in air and arranged in a periodic (triangular or square) lattice with spacing d.
On the basis of ρ _{ℓ}, ϵ _{ℓ},d values and lattice type, the EBG structure presents a gap within a definite frequency band, that represents a forbidden region for fields to propagate. When a line of inclusions is removed from the periodical lattice, a new guiding effect can be obtained. As a matter of fact, the line defect turns into a waveguide for an electromagnetic wave with a frequency within the band gap of the (nondefective) structure. An example of this mechanism is shown in Figure 2, from which it is possible to see that an impinging field at f = 288.46 THz cannot propagate within the full structure because it belongs to the bandgap frequencies, while it propagates when the line defect is realized.
However, the defective EBG structure is able to confine and somehow guide the field, but an unsatisfactory VSWR is obtained. In fact, wave reflections due to the discontinuity in the space produce oscillations of the field inside the guiding region. Notably, this behaviour occurs whatever the kind of waveguide (straight, bent).
Hence, we can exploit the ISMM framework to optimize (some of) the inclusions in order to improve VSWR and get a better design. In particular, we fix the arrangement of the scatterers and their dimension, which is supposed to be identical; then, we look for the optimal permittivity value of inclusions belonging to the socalled ‘design region’ (DR) allowing a minimum ripple for the field amplitude within the waveguide and, in particular, in the socalled ‘shaping region’ (Λ_{shaping }).
Hence, differently from [25], the inverse design problem can be finally formulated as follows:(7.1)
Subject to:(7.2) (7.3)wherein:(8)constraint (7.2) allows to fit the scattering model, while constraint (7.3) enforces a symmetry in the permittivity of the structure.
In the following subparagraphs an assessment of the inverse design procedure is performed for straight and bent waveguides. The considered primary source is a plane wave from the left side of the structure. As we are interested in perspective to reconfigurable devices, the permittivity of the inclusions are the unknowns of our (fixed dimensions) elementary bricks. Of course, a completely similar approach is possible by looking for dimensions (with a fixed permittivity) or both.
As far as the starting point of the optimization problem is concerned, in the specific problem dealt with herein (which can be seen as the problem of optimizing the field which is already propagating within waveguides), the standard defective structure is an obvious (and effective) choice. Also note that constraints on the admissible range values can be considered in the optimization problem to take into account actual feasibility issues.
Fig. 2
Electromagnetic field propagation mechanism for EBG devices. (a) Bandmap in logarithmic scale. Amplitude of electric field with λ = 1.04 μm for (b) nondefective and (c) defective EBG structure. 
4.1 Straight waveguide
As a first example, let us consider the geometry shown in Figure 3a. The full EBG structure is composed of N=168 rods with radius ρ =0.075 μm and permittivity ϵ = 12.25 arranged on a triangular lattice with spacing d = 0.45 μm. The central row composed of 21 rods is removed to realize a straight waveguide. The amplitude of the total electric field on domain for the defective structure is shown in Figure 3b. As previously outlined, the field amplitude achieved within the channel is oscillating, corresponding to a Ripple value equal to 0.3821. To improve such a parameter, the design problem (7) is solved by considering rods on the edge of the channel as DR, see Figure 3a. The optimized structure and the corresponding amplitude field are shown in Figures 3c and 3d, respectively. As it can be seen, a proper adjustment of some inclusions' permittivity allows to reduce the Ripple to 0.0173 (with a comparable value of the average field) and hence to improve the guiding effect.
Fig. 3
Optimization of a straight waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (b) the corresponding electric field amplitude distribution. (c) and (d) the same as (a) and (b) but for the optimized EBG device through the ISMM [ϵ _{ opt }=8.4,10.8,13.2,12.4,10.1,11.4,13.4,10.7,10.2,13.7,12.7,9.6,12.3,16.4,10.9,9.4,16, 16.7,9.3,11.5,18.1]. 
4.2 60° bend waveguide
The second example concerns the optimization of a double 60° bend waveguide. In this case, the total number of scatterers is N = 315, while the structure's geometry as well as the cylinders' characteristics are unchanged. By removing 17 rods as in Figure 4a the amplitude field distribution of Figure 4c is achieved, for which the corresponding Ripple, evaluated for horizontal and diagonal channels separately, is summarized in Table 1. By optimizing the device by using the DR highlighted in Figure 4a, we obtain the device in Figure 4b and a reduced Ripple within the channel, see Figure 4d and Table 1.
Fig. 4
Optimization of a double 60° bend waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (c) the corresponding electric field amplitude distribution. (b) and (d) the same as (a) and (c) but for the optimized EBG device through the ISMM [ϵ _{ opt }=83,77.7,79.2,78.3,79,79.5,82.4,79.9,80,79.7,79.5,82.7,80.1,79.7,78.8,79.3,79,87.1]. 
Quantitative comparison of Ripple for the double 60° bend waveguide example.
4.3 90° bend waveguide
As last example we deal with a double 90° bent waveguide. The rods of the EBG structure (N = 225) have the same parameters as before but for the lattice that has a square arrangement in this case. By removing 20 inclusions as shown in Figure 5a we realize a waveguide with sharp bends, for which the field distribution of Figure 5b is obtained. A similar device is considered in [16] but with rods arranged in a triangular lattice, which is unusual for this kind of sharp bends. Conversely, in [14] the same authors investigate a single 90° bent waveguide but just the corner area is optimized.
In the present case the DR still corresponds to the rods on the channel's edges. As can be apprised from Figure 5d, that is the field pertaining to the optimized structure of Figure 5c, the design through the ISMM allows to reduce the Ripple of the guided field. See Table 2 for a quantitative comparison.
Fig. 5
Optimization of a double 90° bend waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (b) the corresponding electric field amplitude distribution. (c) and (d) the same as (a) and (b) but for the optimized EBG device through the ISMM [ϵ _{ opt }=12.2,14.3,12.2,12.5,11.9,12.5,17.7,11.3,12.4,7.4,8.3,13.7,14.4,14.2,12.9,13.2,12.8,12.9, 11.8,21.6]. 
Quantitative comparison of Ripple for the double 90° bend waveguide example.
5 Conclusions
In this paper the inverse design of artificialmaterialsbased devices has been addressed. In particular, electromagnetic bandgap (EBG) devices have been considered. Differently from approaches relying on fullwave methods or trialanderror procedures, inverse design based techniques represent powerful tools to obtain optimal devices realizing desired properties in a more effective fashion. In this contribution, an inverse design tool for EBG devices based on the scattering matrix method (ISMM) has been proposed. More in detail, the SMM has been adopted to model EBG devices (which are a set of scatterers actually) and hence to develop an ad hoc procedure to determine optimal parameters describing the inclusions composing the device in order to accomplish desired specifications. The proposed ISMM represents a general framework and makes use of simple geometry structures with multiple materials that could be possibly tuned to accomplish with different specification (i.e., with reconfigurable devices). The proposed approach has been successfully assessed by optimizing the guided electromagnetic field within straight and bent waveguides. In particular, by optimizing the permittivity value of inclusions located on the waveguide's edges, our tool is able to control oscillations of the field while keeping an adequate transmitted power.
As we are interested in perspective to reconfigurable devices, we have focused herein on fixed spacings and fixed dimensions of the inclusions. Interestingly, the model we have used can be both immediately used for the case of fixed permittivitiesvariable dimensions as well as extended to the case of variable spacings (by acting on the T matrices rather than on the S matrices).
Acknowledgments
This work was supported in part by the Italian Ministry of University and Research under the PRIN research project “CYBERPHYSICAL ELECTROMAGNETIC VISION: ContextAware Electromagnetic Sensing and Smart Reaction,” prot. 2017HZJXSZ.
Appendix
The present Appendix is devoted to report the analytic expressions of entries of T and Q matrices involved in the SMM model of equation (6).
Accordingly, in the chosen coordinates system, they read:(A1)
wherein represents the mutual position between the chosen ℓth coordinate system (O ^{ℓ} x ^{ℓ} y ^{ℓ}) and the coordinate system of the generic ith inclusion, while, depending on the primary source, a different Q _{ℓ,m } formula must be derived.
In particular, if a plane wave E _{ inc } (x, y) = e ^{−jk b (x sin θ inc −y cos θ inc ) } is considered (θ _{ inc } being the angle with respect to the yaxis) the pertaining expression for Q _{ℓ,m } is the first one in (A2), while if a line current I _{ s } located at (r _{ s }, θ _{ s }) is supposed as primary source, that is, , the second relation in (A2) holds, being the mutual position between the chosen ℓth coordinate system and the source position. Other kind of sources can be obviously accommodated as well.(A2)
References
 A. Silva, F. Monticone, G. Castaldi, V. Galdi, A. Alù, N. Engheta, Performing mathematical operations with metamaterials, Science 343, 160 (2014) [Google Scholar]
 M. Moccia, G. Castaldi, G. D'Alterio, M. Feo, R. Vitiello, V. Galdi, Transformationopticsbased design of a metamaterial radome for extending the scanning angle of a phasedarray antenna, IEEE J. Multiscale Multiphys. Comput. Tech. 2, 159 (2017) [Google Scholar]
 H. Tao, W.J. Padilla, X. Zhang, R.D. Averitt, Recent progress in electromagnetic metamaterial devices for terahertz applications, IEEE J. Sel. Top. Quantum Electron. 17, 92 (2010) [Google Scholar]
 F. Frezza, L. Pajewski, E. Piuzzi, C. Ponti, G. Schettini, Radiation enhancement properties of an Xband woodpile EBG and its application to a planar antenna, Int. J. Antennas Propag. 2014, 1 (2014) [Google Scholar]
 A.R. Weily, K.P. Esselle, T.S. Bird, B.C. Sanders, Experimental woodpile EBG waveguides, bends and power dividers at microwave frequencies, Electron. Lett. 42, 3 (2006) [Google Scholar]
 G. Torrisi et al., Numerical study of photoniccrystalbased dielectric accelerators, in Proceedings of 10th International Particle Accelerator Conference (IPAC), Melbourne, VIC, Australia, 2019, JACOW, Geneva, Switzerland, pp. 3653 [Google Scholar]
 F. Yang, Y. RahmatSamii, Electromagnetic Band Gap Structures in Antenna Engineering, ser. Cambridge books online (Cambridge University Press, 2008) [Google Scholar]
 L. Crocco, F. Cuomo, T. Isernia, Improved scattering matrix method for the analysis of twodimensional PBG devices, Microw. Opt. Technol. Lett. 48, 2564 (2006) [Google Scholar]
 C.M. Reinke, M. Teofilo, M.F. Su, M.B. Sinclair, I. ElKady, Grouptheory approach to tailored electromagnetic properties of metamaterials: an inverseproblem solution, Phys. Rev. E 83, 066603 (2011) [Google Scholar]
 Y. ElKahlout, G. Kiziltas, Inverse synthesis of electromagnetic materials using homogenization based topology optimization, Prog. Electromagn. Res. 115, 343 (2011) [Google Scholar]
 D.R. Smith, D.C. Vier, T. Koschny, C.M. Soukoulis, Electromagnetic parameter retrieval from inhomogeneous metamaterials, Phys. Rev. E 71, 036617 (2005) [Google Scholar]
 M.P. Bendsoe, O. Sigmund, Topology optimization: theory, methods, and applications (Springer Science & Business Media, 2013) [Google Scholar]
 L.H. Frandsen, A. Harpøth, P.I. Borel, M. Kristensen, J.S. Jensen, O. Sigmund, Broadband photonic crystal waveguide 60 bend obtained utilizing topology optimization, Opt. Express 12, 5916 (2004) [Google Scholar]
 J.S. Jensen, O. Sigmund, Systematic design of photonic crystal structures using topology optimization: lowloss waveguide bends, Appl. Phys. Lett. 84, 2022 (2004) [Google Scholar]
 J.S. Jensen, O. Sigmund, Topology optimization of photonic crystal structures: a highbandwidth lowloss Tjunction waveguide, JOSA B 22, 1191 (2005) [Google Scholar]
 J.S. Jensen, O. Sigmund, L.H. Frandsen, P.I. Borel, A. Harpoth, M. Kristensen, Topology design and fabrication of an efficient double 90/spl deg/photonic crystal waveguide bend, IEEE Photon. Technol. Lett. 17, 1202 (2005) [Google Scholar]
 M. Otomori, J. Andkjaer, O. Sigmund, K. Izui, S. Nishiwaki, Inverse design of dielectric materials by topology optimization, Prog. Electromagn. Res. PIER 127, 93 (2012) [Google Scholar]
 F. Callewaert, V. Velev, P. Kumar, A.V. Sahakian, K. Aydin, Inversedesigned broadband alldielectric electromagnetic metadevices, Sci. Rep. 8, 1358 (2018) [Google Scholar]
 J. Han, J. Huang, J. Wu, J. Yang, Inverse designed tunable fourchannel wavelength demultiplexer, Opt. Commun. 465, 125606 (2020) [Google Scholar]
 M.H. Tahersima, K. Kojima, T. KoikeAkino, D. Jha, B. Wang, C. Lin, K. Parsons, Deep neural network inverse design of integrated photonic power splitters, Sci. Rep. 9, 1368 (2019) [Google Scholar]
 D. Felbacq, G. Tayeb, D. Maystre, Scattering by a random set of parallel cylinders, J. Opt. Soc. Am. A 11, 2526 (1994) [Google Scholar]
 R. Palmeri, M.T. Bevacqua, A.F. Morabito, T. Isernia, Design of artificialmaterialbased antennas using inverse scattering techniques, IEEE Trans. Antennas Propag. 66, 7076 (2018) [Google Scholar]
 R. Palmeri, M.T. Bevacqua, T. Iserni, Design of invisibility devices through artificial materials: further possible tools from the inverse scattering perspective, 2019 Thirteenth International Congress on Artificial Materials for Novel Wave Phenomena (Metamaterials), IEEE, 2019, pp. X–290 [Google Scholar]
 R. Palmeri, T. Isernia, Volumetric invisibility cloaks design through spectral coverage optimization, IEEE Access 7, 30860 (2019) [Google Scholar]
 R. Palmeri, T. Isernia, Inverse design of artificial materials based lens antennas through the scattering matrix method, Electronics 9, 559 (2020) [Google Scholar]
 A. Hakansson, J. SánchezDehesa, L. Sanchis, Inverse design of photonic crystal devices, IEEE J. Sel. Areas Commun. 23, 1365 (2005) [Google Scholar]
 D.S. Jones, Acoustic and electromagnetic waves (Oxford Science Publications, 1986) [Google Scholar]
 O.M. Bucci, T. Isernia, Electromagnetic inverse scattering: retrievable information and measurement strategies, Radio Sci. 32, 2123 (1997) [Google Scholar]
 M. Abramovitz, I. Stegun, Handbook of Mathematical Functions (Dover, New York, 1970) [Google Scholar]
A truncation rule can be derived by considering the singular value decomposition of the radiation operator. For the present case, M _{ℓ} = 2k _{ b } ρ _{ℓ} + 1, ρ _{ℓ} being the radius of the ℓth inclusion [28].
Cite this article as: Roberta Palmeri, Tommaso Isernia, Inverse design of EBG waveguides through scattering matrices, EPJ Appl. Metamat. 7, 10 (2020)
All Tables
All Figures
Fig. 1
Geometry sketch for the problem at hand. 

In the text 
Fig. 2
Electromagnetic field propagation mechanism for EBG devices. (a) Bandmap in logarithmic scale. Amplitude of electric field with λ = 1.04 μm for (b) nondefective and (c) defective EBG structure. 

In the text 
Fig. 3
Optimization of a straight waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (b) the corresponding electric field amplitude distribution. (c) and (d) the same as (a) and (b) but for the optimized EBG device through the ISMM [ϵ _{ opt }=8.4,10.8,13.2,12.4,10.1,11.4,13.4,10.7,10.2,13.7,12.7,9.6,12.3,16.4,10.9,9.4,16, 16.7,9.3,11.5,18.1]. 

In the text 
Fig. 4
Optimization of a double 60° bend waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (c) the corresponding electric field amplitude distribution. (b) and (d) the same as (a) and (c) but for the optimized EBG device through the ISMM [ϵ _{ opt }=83,77.7,79.2,78.3,79,79.5,82.4,79.9,80,79.7,79.5,82.7,80.1,79.7,78.8,79.3,79,87.1]. 

In the text 
Fig. 5
Optimization of a double 90° bend waveguide. (a) Permittivity map of the defective EBG device where the DR is superimposed with light grey rectangle and (b) the corresponding electric field amplitude distribution. (c) and (d) the same as (a) and (b) but for the optimized EBG device through the ISMM [ϵ _{ opt }=12.2,14.3,12.2,12.5,11.9,12.5,17.7,11.3,12.4,7.4,8.3,13.7,14.4,14.2,12.9,13.2,12.8,12.9, 11.8,21.6]. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.