Novel Design Framework for Dual-Band Frequency Selective Surfaces Using Multi-Variant Differential Evolution

Frequency Selective Surfaces (FSSs) have become increasingly popular during the last years due to their combined characteristics, which meet, in general, the requirements of the next-generation wireless communication networks. In this work, a cross-platform design framework for FSS structures is presented and evaluated by utilizing a recently introduced evolutionary optimization algorithm, namely, the Multi-Variant Differential Evolution (MVDE). To the best of the authors knowledge, this is the first time that the MVDE algorithm is applied to a design problem in Electromagnetics. The proposed design framework is described in detail and the utilized evolutionary algorithm is assessed in terms of its performance by applying several benchmark functions. In this context, the MVDE is comparatively evaluated against other popular evolutionary algorithms. Moreover, it is applied to the design and optimization of two different representative examples of FSS structures based on three use cases of unit cell geometry. Optimization results indicate the efficacy of the proposed framework by quantifying the performance of the designed FSS structures in terms of several system metrics. The optimized FSS structures exhibit dual-band operation and quite acceptable results in the ISM frequency bands of 2.45 GHz and 5.8 GHz.


Introduction
Fifth-generation (5G) mobile communication systems have been started to deploy worldwide almost for two years. These systems, both from the operator and the end user perspective, require several frequency bands with Multiple-Input Multiple-Output (MIMO) architecture to operate properly [1][2][3]. Within this complex environment, various characteristics, such as spatial filtering [4], isolation [5], and decoupling [1], are of importance in the current and the next mobile communication networks, as several heterogeneous services with diverse features and multiple-frequency requirements are combined. Frequency Selective Surfaces (FSS) are electromagnetic (EM) devices that can integrate the above characteristics [6]. In the forthcoming years, Next-Generation Wireless Networks (NGWN) will face several challenges in their deployment, including ultra-low power con-and 31% for a 4-unit supercell, and RF-to-DC conversion efficiency of 27% and 29% for an 8-unit supercell, both at the frequency bands of GSM-1800 and Wi-Fi 2.45 GHz, accordingly.
Evolutionary optimization algorithms [73,74] are one of the three main categories of population-based meta-heuristic algorithms; the other two are physics-based [75] and swarm-intelligence algorithms [76][77][78]. Evolutionary algorithms (EAs) adopt mechanisms that are usually inspired by the biological evolution of natural environments, such as mutation, crossover, and selection. Each solution to the optimization problem acts as an individual candidate. The result of the cost function for each of the available solutions determines the activation of the previously mentioned mechanisms. In the course of iterations, the evolution of the population is taking place (based on these mechanisms) towards the convergence of the optimal solution. Various meta-heuristic optimization algorithms, such as the Genetic Algorithm (GA) [79,80], the Differential Evolution (DE) [81], the Self-Adaptive Differential Evolution [82], the Particle Swarm Optimization (PSO) [83], the Grey Wolf Optimizer (GWO) [84], and the Wind Driven Optimization (WDO) [85], have been considerably exploited by their application to the field of Electromagnetics [85][86][87][88][89][90].
Evolutionary optimization algorithms have been also applied to the design of FSS structures. In [91], an optimization process based on GAs [79,80] for the design of dualband FSS structures was presented. The shape of the unit cell was optimized, and the polarization performance was assessed. A novel methodology for reconfigurable FSSs was introduced in [92] by utilizing the GAs. The optimization technique was applied to obtain the optimal choice of the switches to produce the desired frequency response. The authors of [93] discussed the synthesis of an FSS structure with a low cross-polarization level by utilizing the parallel binary-coded micro-genetic algorithm [94]. They optimized the pattern and the periodicity of the FSS along the x-and y-directions to minimize the cross-polarization level for a wide frequency band and a considerably large range of angles. The same algorithm was also applied to optimize FSS structures that operate as spatial filters [95] or as metamaterials in multi-frequency bands [96]. The combination of the genetic algorithm and the geometry-refinement technique [97] was introduced in [98] to design and optimize a multiband single-layer FSS. The derived result exhibited broader characteristics in terms of bandwidth operation.
Moreover, the DE [81] as well as the Covariance Matrix Adaptation Evolution Strategy (CMA-ES) [99] have been utilized in the design of FSS structures. In [100], the authors applied the optimization technique of differential evolution strategy to design an FSS structure. They also performed a comparison between the DE and the GA. From their results, they derived that the DE strategy exhibited comparative advantages against the GA. The same optimization technique (DE strategy) was also utilized in [101] to design a double-square loop as a unit cell in an FSS structure. Additionally, the authors of [102] applied the CMA-ES to efficaciously design and optimize the metallic shape of a unit cell in an FSS structure. Finally, FSS structures have been designed by the use of Swarm Intelligence (SI) algorithms. One of the most widely known algorithms, i.e., the PSO [83], has been applied to synthesize FSS structures [103][104][105]. Furthermore, the Multi-Objective Lazy Ant Colony Optimization (MOLACO), which is an extended version of Ant Colony Optimization (ACO) technique [106], was utilized in the design of a 3D FSS design [107]. Recently, the authors of [108] designed a dual-band single-layer FSS structure by utilizing the Harris Hawks Optimization (HHO) technique [109].
Various optimization tools and frameworks, mostly open-source, with applications in Electromagnetics, have been introduced in the literature over the past years. Among all, the authors of [110] presented two optimization packages for designing various technical devices and systems: the Agros environment and theĀrton optimization toolbox. The former represents an environment of systems consisting of Partial Differential Equations (PDEs); it is a multi-physical Finite Element Method (FEM) solver that supports various 3D problems in Engineering, including Electromagnetics. The latter is basically an optimization toolbox for Python. It provides an efficient programming environment for a large variety of optimization methods. An open-source C++ library using FEM techniques with high-order meshing is described in [111]. MFEM (Modular Finite Element Methods) supports various algorithms, solvers, as well as hardware platforms that are utilized by the scientific community. Finally, OpenEMS [112] is a free and open-source electromagnetic solver based on the Finite Difference Time Domain (FDTD) method. It supports various features, such as 3D mesh, multi-threading, remote simulations capability, and MATLAB interface. If we consider the main features of various optimization tools and frameworks presented in the literature, we can conclude that their common advantage is the provision of open-source distribution. Most of them are based on a programming language (Python, C++, etc.) in order to operate properly. However, their main limitation is the lack of support in cross platforms.
In this work, a new cross-platform framework for designing and optimizing FSS structures as applied to various problems in Electromagnetics is analyzed and validated. The proposed design framework combines a recently introduced evolutionary optimization algorithm, namely, the Multi-Variant Differential Evolution (MVDE) algorithm [113], and a commercial high-frequency electromagnetic solver (HFSS; © 2020 ANSYS, Inc., Canonsburg, PA, USA). The novelty of our work lies in the fact that a complete optimization framework for designing FSS structures, and consecutively, various EM models, such as antennas, antenna arrays, microwave filters, etc., is provided. To the best of the authors' knowledge, this is the first time that the MVDE algorithm is applied to solve an optimization problem in the field of Electromagnetics. Furthermore, a contributing characteristic of this work originates from the utilization of a pair of U-slots apertures in a corresponding set of patches as a unit cell in the optimization process of an FSS structure. Moreover, to the best of the authors' knowledge, this is the first dual-band FSS absorber, which is suitable for RF energy harvesting applications, that is operating in the ISM Wi-Fi 2.45 GHz and Wi-Fi 5.8 GHz frequency bands. The design of an FSS structure, especially when the unit cell holds a geometry with a medium of high complexity, is a problem in Electromagnetics that can be classified as an optimization one. Taking into consideration the requirements of the next-generation communication networks in multi-frequency band operation, the design of dual-or multi-band FSS structures is a straightforward optimization process since the number of parameters (decision variables in an optimization problem) is considerably large.
The remainder of this work is structured as follows. Section 2 includes the applied methods and techniques for the presentation and the evaluation of the introduced design framework. Specifically, in Section 2.1, the problem of the optimization process is defined; in Section 2.2, the proposed design framework is analytically described; in Section 2.3, a brief description of the utilized MVDE algorithm is included; and in Section 2.4, the applied algorithm is assessed in terms of its performance against 4 widely-known evolutionary algorithms. Section 3 presents and discusses the optimization results obtained by the utilization of the proposed framework to design and optimize two representative examples of FSSs structures that exhibit dual-band frequency operation. In detail, Section 3.1 describes the parameters that are considered to set up the optimization process for designing both the unit cell and the FSS structure, whereas Sections 3.2 and 3.3 include the main results based on several system metrics obtained by the optimization process. Finally, Section 4 concludes the findings and outlines the future steps of this work. Figure 1 pictures the generalized approach of a single-layer FSS structure. Without losing the generalization of the problem definition, a single-layer FSS design (a similar problem definition can be expressed for dual-or even triple-layer FSS structures) consists of a unit cell (usually of metal) which is periodically repeated along a two-dimensional lattice. Beneath the lattice of the unit cells, a grounded substrate layer (usually of dielectric) is placed. Based on the electromagnetic theory, when an incident field triggers the FSS structure, the latter one performs as an absorber [13], a reflector [14], or a transmitter [12], based on the frequency of the incident field. The shape and the size of the unit cell is usually the main challenging and demanding task in the optimization process [114] as they determine the resonance of the FSS structure in the desired frequency band and, consecutively, its performance as an aftereffect of its functional characteristics.

Optimization Process
In this work, a new cross-platform design framework for FSSs based on evolutionary optimization algorithms is presented. The proposed framework combines the advantages of a commercial high-frequency electromagnetic solver (HFSS; © 2020 ANSYS, Inc., Canonsburg, PA, USA) and the robustness of an evolutionary optimization algorithm, i.e., the MVDE [113], which was recently introduced in the literature and is briefly described in Section 2.3. The proposed framework is a generalized approach for obtaining feasible solutions to the design of antenna structures based on optimization algorithms. Figure 2 displays the general concept and the main steps of the proposed design framework that can be applied in several antenna optimization problems and can be combined with various evolutionary optimization algorithms. We should point out that the proposed framework can be utilized with several different meta-heuristic algorithms, and with other optimization problems as well if certain modifications are applied.
Analyzing the proposed design framework that is displayed in Figure 2, the first step includes the initialization of the selected optimization technique, i.e., the MVDE algorithm in our case. Prior to the initialization of the algorithm (i.e., the definition of (a) the population number, (b) the number of the independent trials, (c) the dimensionality of the optimization problem, and (d) the control parameters), the decision variables that fully described the EM structure geometry must be defined, as well as their boundaries to the given optimization problem. Moreover, the definition of the EM model (i.e., in our case the FSS structure models, both unit-cell and total structure), as well as the systems metrics to be optimized, are also required. Based on the optimization problem and the model of the EM structure, various system metrics can be utilized, such as the magnitude of the reflection coefficient (S 11 ), the realized gain, the characteristic impedance, or the efficiency of the proposed EM structure [86,88,90,114]. In this work, we have selected the magnitude of the reflection coefficient as the system metric to be optimized. However, for each of the given solutions, the characteristic impedance of the FSS structure, as well as the realized gain at the desired frequencies are also computed and evaluated. The second step in the optimization process creates a set of values for the decision variables of the optimization problem [89], and, at each iteration, the population of the selected algorithm is initialized. Additionally, at each population member, a design solution based on the EM model (e.g., FSS structure) is created and parsed to the high-frequency EM solver, which computes the selected system metrics at the frequencies of interest. The results are reported back to the optimization algorithm to compute the objective function of the problem [86]. These results are stored temporarily, until the threshold criteria are applied. Once the computation of the objective function is completed for every member of the entire population, the results are populated to the optimization algorithm. Based on the selected algorithm (i.e., MVDE), the best values of the position and the objective function for the members of the population are stored [86], and the specific mechanisms for parent selection and crossover mutation are taking place. The current iteration is completed by updating the position of the population, and by storing the results for the next iteration.
The optimization process is concluded when the stopping criteria are met, i.e., when the number of maximum iterations or maximum function evaluations is achieved or the system metrics values are satisfied [86,88,89]. As a result, the feasible solution of the optimization problem is extracted and the final design of the EM structure is obtained. Based on the final design, a set of system metrics (optimized and/or not) is computed to obtain the final results of the proposed EM structure geometry.

MVDE Algorithm Description
MVDE algorithm has been recently introduced in the literature [113]. It is a stochastic and population-based evolutionary algorithm that combines five different variants of the legacy DE algorithm [81] to balance the exploration and the exploitation phases of searching the computational space during the iterative process. The utilized DE variants in MVDE are DE/rand/1, DE/current-to-best/1, DE/rand/2, DE/best/2, and DE/best/1 [115].
The multi-variant mutation schemes are incorporated with a self-adaptive scaling factor based on the probabilistic cosine and logistic distributions. The introduction of the adaptive scaling factor in the MVDE algorithm is applied to enhance the performance of the legacy DE algorithm, i.e., to mitigate the probability of trapping in local minima and to avoid premature convergence. The pseudocode of the MVDE algorithm is outlined in Algorithm 1.
To briefly describe the optimization technique of the MVDE algorithm, let us consider the definition of N pop as the population number (k = 1...NPop), D as the number of decision variables (j = 1...D), MaxVar and MinVar as the upper and the lower boundary of the decision variables, and MAxIt as the maximum number of iterations (i = 1...MaxIt). At the first iteration (i = 1), the initial values of the position u k and the corresponding objective function OF k for each member of the population are selected. A ranking process determines the best values of the position and the corresponding objective function for each member of the population. For all the following iterations (i > 1), the crossover factor (CF) is computed by At the current iteration i, and for each member of the population, the self-adaptive scaling factor [SF] k,D is computed based on the cosine or the logistic distribution. It is expressed by where [rand] k,D is a uniformly distributed random number. The derived self-adaptive scaling factor is further computed by applying a binary matrix, which is given by Therefore, the expression of the self-adaptive scaling factor is defined as Moreover, the parents' selection is based on the expression The next step in the optimization technique of the MVDE algorithm is the computation of the position for each member of the population at the next iteration, based on one of the five different variants of the legacy DE algorithm. The position is expressed by Finally, based on the position [u] i+1 k,D , the objective function OF i+1 k of the optimization problem is computed for each member of the population. The best values of the position u best k and the corresponding objective function OF best k are updated and the process is iterated until the stopping criteria are met. Ultimately, the feasible solution of the EM structure is extracted and the final system metrics are computed.

Performance Evaluation
The MVDE algorithm in the proposed design framework has been evaluated in terms of its performance by utilizing four widely known evolutionary algorithms, i.e., the GA [79,80], the Biogeography-Based Optimization (BBO) [116], the DE [81], and the CMA-ES [99]. The performance score for each of the aforementioned algorithms is assessed by utilizing ten commonly used benchmark functions: Ackley ( f 1 ), Griewank ( f 2 ), Rastrigin  Table 1 lists the obtained results of the performance evaluation for the MVDE algorithm in conjunction with four popular evolutionary algorithms by computing the average values of the cost function of 100 independent trials, having the number of the decision variables equal to 30 and 50. From the obtained results, we can infer that the utilized MVDE algorithm that is integrated into the proposed design framework outperforms in both cases of the solutions to the optimization problem (i.e., D = 30 and D = 50) and for the whole set of the selected benchmark functions. The MVDE algorithm exhibits its effectiveness to balance between the exploration and the exploitation in the search space, and its robustness to efficiently converge to the global minima in a variety of benchmark functions. Compute the crossover factor CF using (1) 6: for (k = 1; k + +; k ≤ NPop) do 7: if (rnd > 3 × CF) then 8: Compute the self-adaptive scaling factor [SF] k,D by applying the cosine distribution of (2) 9: else 10: Compute the self-adaptive scaling factor [SF] k,D by applying the logistic distribution of (2) 11: end if 12: Compute the binary matrix using (3) 13: Compute the self-adaptive scaling factor [SF] k,D using (4) 14: Compute the parent selection in each iteration using (5) 15: if (i < 0.2 × MaxIt) then 16: Compute the position [u] i+1 k,D using the DE/rand/1 variant of (6) 17: Update the best values of u best k , OF best k 28: end for 29: end for Table 1. Performance metric (average values of the cost function of 100 independent trials) of the MVDE algorithm [113] in conjunction with four popular evolutionary algorithms (GA [79,80], BBO [116], DE [81], and CMA-ES [99]) and for each of the selected benchmark functions (the best values are marked in bold). To quantify the derived conclusion based on the presented results, the Friedman nonparametric statistical test is applied as a meta-performance metric to assess the efficacy for each of the selected algorithms based on the mean ranking. Table 2 summarizes the obtained results. We can easily conclude that the MVDE algorithm achieves the best score of the mean ranking, having the second-best score by CMA-ES. We should also point out that the legacy DE algorithm achieves the last score of the mean ranking in the Friedman test, mostly due to the drawbacks of trapping in local minima and of early convergence. Table 2. Meta-performance metric (Friedman nonparametric statistical test) of the MVDE algorithm [113] in conjunction with four popular evolutionary algorithms (GA [79,80], BBO [116], DE [81], and CMA-ES [99]) based on the results of Table 1  From the presented results, we can conclude that the MVDE algorithm, which is utilized in the proposed design framework, exhibits considerable stability and robustness in a variety of benchmark functions, i.e., functions with many local minima, bowl-shaped, valley-shaped, steep ridges, or other characteristics. We should also point out that the 25th percentile and the 75th percentile of the given distributions for the MVDE algorithm achieve an indistinguishable variation (for any of the selected benchmark functions).

Optimization Results and Discussion
The proposed design framework, by incorporating a recently introduced evolutionary algorithm, namely, the MVDE algorithm [113], and a commercial high-frequency electromagnetic solver (HFSS; © 2020 ANSYS, Inc., Canonsburg, PA, USA), has been applied in an iterative process to design, optimize, and obtain a feasible solution of an FSS structure suitable for RF EH applications. To this end, the absorptive characteristic of the FSS structure is exploited and the corresponding system metrics are utilized in the optimization process. Sections 3.1-3.3 present the optimization setup that has been applied in the process, the general concept, as well as the main steps in designing and optimizing an FSS structure with a complex unit cell.

Optimization Setup
The optimization setup of designing a unit cell, and consecutively an FSS structure based on the selected unit cell, includes the definition of the unit cell geometry, i.e., the type of the absorber and the decision variables that fully describe its geometry [114]. Figure 5 portrays the proposed geometry of a unit cell as an EM structure. The proposed unit cell that forms an absorber (Figure 5a) consists of a pair of symmetrical patches (two axes of symmetry occur), which harvest the EM radiation of the surrounding environment in the desired frequency bands. Each of the two patches in the proposed unit cell includes a U-slot within a predefined by the patch surface, to achieve the second resonance of the total EM structure (the first resonance occurs from the patch dimensions). The harvested energy is directed to the terminal port (which is the reference point where the energy is converted) of the EM structure through the microstrip lines, one for each of the corresponding patches. The described geometry of the proposed unit cell can be repeated in both of the axes that define the plane of the unit cell (in our case this is the XY plane as Figure 5a indicates).
Based on the absorber of Figure 5a, three use cases of the proposed unit cell are determined and included in the optimization process. These are portrayed in Figure 5b-d. The first use case of the proposed unit cell exhibits two axes of symmetry (axes X and Y of the Figure 5b), where the U-slot is centered along the X-axis. In the second use case, the Y-axis of symmetry is omitted, thus allowing to offset the U-slot within the patch of the absorber, by introducing the Ox s decision variable (Figure 5c). The third use case of the proposed unit cell differentiates the starting point of the Ox s decision variable that affects the offset of the U-slot (Figure 5d).
If we take into consideration the Figure 5b, we can conclude that the number of decision variables that are required to fully describe the geometry of the proposed unit cell is equal to 10. We should also clarify that D is equal to 10 only for the first use case of the proposed unit cell, the second and the third use cases require an extra decision variable for their geometry description; therefore, in these cases, D is equal to 11. Such a geometry is practically impossible to describe analytically using the electromagnetic theory; the use of an optimization method to obtain a feasible solution of a unit cell and consecutively, of a dual-band FSS design as an EM structure, is a straightforward process. In this work, we analyze the operation of the proposed design framework by combining an evolutionary optimization algorithm (MVDE) [113] and a commercial electromagnetic solver (HFSS).
To design and optimize the performance of various use cases of a unit cell and the corresponding FSS structure, we exploit the characteristic of absorptivity [13]. This characteristic results in EM structures with maximum values of the magnitude of absorptivity in the frequency bands of interest or, equivalently, with minimum values of the magnitude of the reflection coefficient [114]. As a result, the objective of the given optimization problem can be reworded as the minimization of the magnitude of the reflection coefficient at the frequency bands of interest for each of the designed EM structures (unit cell, FSS structure). By minimizing the reflection coefficient at the terminal port of the EM structure, the impedance matching between the designed EM structure and the reference point (in practical RF EH systems, the reference point is referred to the module that transforms the RF energy to DC voltage, i.e., the RF-to-DC rectifier in a rectenna system) occurs, and con-secutively, the impedance of the designed EM structure tends to reach the characteristic value of 50 Ω. In this work, we express the given optimization problem by minimizing the reflection coefficient as the system metric of the designed EM structures (objective function of the optimization problem); however, the system metrics of the impedance, the radiation pattern, and the realized gain are also considered [117,118] to extract the feasible solution of the optimization problem. During the initialization phase of the optimization process in the proposed design framework, the following parameters are applied to the MVDE algorithm: The number of function evaluations, which acts as a stopping criterion in the optimization process, is equal to 10,000, and the threshold criterion for an acceptable solution of the optimization problem is set to −10 dB. As a result, the objective function of the previously described optimization problem can be formulated as where • u is the position vector for each member of the population of the utilized MVDE algorithm, • S i 11 (i = {2.45 GHz, 5.8 GHz}) is the system metric (magnitude of the reflection coefficient) of the designed EM structure at the specific frequencies of interest, • T dB is the threshold criterion for an acceptable solution of the optimization problem provided by the members of the population (in our case T dB =−10 dB), • Ψ is a positive number (multiplying factor in the objective function) that is triggered when the obtained solution is above the threshold criterion, and • OF is the objective function of the optimization problem. Figure 6 illustrates the final design solution of the unit cell for each of the use cases that are obtained by the proposed design framework and the applied optimization process. Moreover, Table 3 lists the final values of the decision variables that fully describe the geometry for each of the presented unit cell use case of Figure 6.  In order to design an FSS structure suitable for RF EH applications, the performance results of the unit cell use cases should be assessed. Figure 7 displays the comparative results of the S 11 magnitude of the reflection coefficient (final values of the system metric extracted from the obtained EM model using the proposed design framework and the optimization process of Figure 2) for each of the unit cell use cases. From the presented results, we can conclude that the first use case of the unit cell design achieves the best performance in terms of its reflection coefficient. In detail, the magnitude of the reflection coefficient is −48.    Table 4 lists the extracted results of the impedance at the frequencies of interest. From the presented results, we can infer that the first use case of the designed and optimized unit cell use cases exhibits the best results against the other two solutions in terms of its matching performance to the characteristic impedance of 50 Ω.  Figure 9 depicts the comparative results of the radiation pattern at the XZ and YZ planes of interest for each of the unit cell use cases. From the depicted results, we can deduce that the first and the third use cases of the unit cell exhibit the best performance in terms of its beamwidth operation. For example, both the first and the third use cases achieve a half-power bandwidth (HPBW) larger than 105 deg for the XZ plane and at the frequency band of 2.45 GHz (the corresponding HPBW of the second use case is smaller than 100 deg). Consecutively, the HPBW for the first and the third use cases is larger than 90 deg for the XZ plane and the frequency band of 5.8 GHz, whereas the second use case achieves an HPBW of less than 55 deg for the same main plane and frequency band, respectively. Table 4. Obtained results of the impedance values for each of the unit cell use cases at the frequencies of interest (extracted results from Figure 8).  In conclusion, the first of the presented unit cell use cases in the design of an FSS structure by utilizing the proposed design framework exhibits comparative benefits in terms of the main system metrics that are taking into consideration in an EM structure suitable for RF EH applications. Therefore, the geometry of the first use case will be utilized as a unit cell in two different FSS structures presented in Section 3.3. Figure 10 illustrates two representative designed and optimized solutions of an FSS structure based on the design, optimization, and the comparative results of the first of the unit cell use cases presented in Section 3.2. The first one (Figure 10a) is an FSS structure of 3 × 2 unit cells, whereas the second solution (Figure 10b) consists of 5 × 3 unit cells. Such a design, operating in the ISM frequency bands of 2.45 GHz and 5.8 GHz, and having a relatively small size, makes it a strong candidate for a receiving module in practical indoor RF EH applications. In Section 2.4, the MVDE algorithm that is included in the proposed design framework is assessed in terms of its performance. To this end, four evolutionary optimization algorithms, namely, the GA [79,80], the BBO [116], the DE [81], and the CMA-ES [99], and 10 popular benchmark functions, are utilized to compare their performance. To further validate our proposed framework when designing an FSS structure, a separate optimization process is carried out by utilizing the legacy DE algorithm [81] (besides the optimization process where the proposed MVDE algorithm [113] is applied), in order to compare and evaluate their performances using the selected system metrics of the EM structure. We should also point out that, in the optimization process of the proposed design framework, the utilized algorithms (MVDE, DE) are applied using the same parameters, as described in Section 3.1. The result of the comparison between the MVDE and the legacy DE algorithm in the optimization process of the proposed framework is depicted in Figure 11. From the displayed curves of Figure 11 , we can conclude that the optimized 3 × 2 and 5 × 3 FSS structures resonate at the desired frequency bands of Wi-Fi 2.45 GHz and 5.8 GHz. Moreover, we can infer that the result for the magnitude of the reflection coefficient of the MVDE algorithm outperforms the result of the legacy DE algorithm both for the two representative solutions of FSS structures, i.e., 3 × 2 FSS design and 5 × 3 FSS design, at the frequency bands of interest. Table 5 summarizes the performance comparison between the two evolutionary algorithms by listing the minimum values of the reflection coefficient in the desired frequency bands. From the depicted results, we can conclude that the incorporation of the recently introduced MVDE algorithm to the optimization process of the proposed design framework obtains an improvement in the magnitude of the reflection coefficient at the frequency bands of interest. In detail, if we consider the 3 × 2 FSS design, there is an improvement of 4.14 dB and 6.59 dB for the ISM Wi-Fi frequency bands of 2.45 GHz and 5.8 GHz, accordingly. Additionally, for the 5 × 3 FSS design, the improvement raises to 5.24 dB and 10.27 dB, for the same frequency bands of interest, accordingly.  Figure 12 displays the magnitude of the S ij parameters (i = j, 3 × 2 FSS design: (i, j = 1, 2), 5 × 2 FSS design: (i, j = 1, 2, 3)) that is obtained by utilizing the aforementioned optimization process and the proposed design framework. From the displayed results, we can conclude that the mutual coupling (usually the system metric that express the isolation between two or more terminal ports is the S ij parameter) is below −20 dB between the two and the three terminal ports, for the use case of 3 × 2 FSS design and 5 × 3 FSS design, respectively. Therefore, the isolation in both cased between the reference points is kept at acceptable values for the whole frequency bands of interest.  Finally, Table 6 summarizes the comparative results of this work against selected published works from the literature. The selection was made based on the functional characteristics of the published FSS structures and the provided data. The parameters for the comparison are the PCB substrate, the type of the unit cell, the frequency band of operation, the derived layout of the FSS structure, the total occupied area of the FSS structure expressed in wavelengths with respect to the lowest resonant frequency of the FSS, the magnitude of the reflection coefficient, and the reference frequency of the S 11 parameter. From the listed results of Table 6, we can infer that most of the published work of FSS structures that are designed for RF energy harvesting applications resonate in a single frequency band. Only the work in [72] and our work exhibit dual-band operation, however in different frequency bands. Furthermore, the Wi-Fi 2.45 GHz frequency band seems to be dominant in most of the works in the literature. Moreover, various types of unit cells have been utilized for the implementation of FSS structures, with most of them having low to medium complexity. Additionally, the PCB substrate of FR-4 and Rogers RO4003C are the most applied layers for the design of FSS structures. However, the occupied area in terms of wavelength units is considerably large in most of the works. Our proposed FSS structure occupies the smallest PCB area with respect to the wavelength of the minimum resonant frequency (i.e., 2.45 GHz), with the FSS in [72] having a similar size. As for the magnitude of the reflection coefficient, our proposed FSS structure achieves quite satisfactory values among the selected published works, thus making it a promising candidate for RF energy harvesting applications.

Conclusions
In this work, a new cross-platform design framework for Frequency Selective Surfaces by exploiting the robustness of evolutionary optimization algorithms is introduced and evaluated. The main advantage of utilizing such a design framework is the capability of designing FSS structures with high complexity, i.e., FSS structures that are based on unit cells with complex shapes. Furthermore, the application of the design framework raises the advantage of developing complex FSS structures with multi-band and even, multifunctional operations. Moreover, the incorporation of modern meta-heuristic algorithms in the optimization process of the design framework can mitigate various optimization limitations that lead to poor convergence, such as trapping in local minima or premature convergence. However, the utilization of the proposed design framework holds some limitations. One of the most important limitations is the proper definition of the boundaries in the optimization problem. Furthermore, the choice of an optimization algorithm in the proposed design framework is of importance as there is no suitable algorithm for solving all optimization problems. Within this context, the MVDE evolutionary algorithm is combined along with a commercial high-frequency electromagnetic solver in the proposed framework to design and optimize two representative examples of FSS structures based on three different use cases of unit cell geometry. The optimized FSS structures exhibit satisfactory results and dual-band operation in the ISM frequency bands of 2.45 GHz and 5.8 GHz. Specific system metrics are computed to assess the performance of the derived FSS structures, i.e., the reflection coefficient, the mutual coupling, the realized gain, and the HPBW. From the obtained results, we can conclude that both of the FSS structures resonate at the desired frequency bands, exhibit acceptable values of isolation between their terminal ports, present satisfactory values of realized gain, and achieve high values of HPBW. Future work includes the integration of the proposed design framework with other promising and recently introduced evolutionary optimization algorithms, the extension of the design framework to include multi-objective optimization problems of Electromagnetics, and the assessment of the presented design framework by optimizing EM structures with even more complex geometries.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest:
The authors declare no conflict of interest.