Metamaterial-Based Sub-Microwave Electromagnetic Field Energy Harvesting System

This paper presents the comprehensive analysis of the sub-microwave, radio frequency band resonant metastructures’ electromagnetic properties with a particular emphasis on the possibility of their application in energy harvesting systems. Selected structures based on representative topologies of metamaterials have been implemented in the simulation environment. The models have been analyzed and their substitute average electromagnetic parameters (absorption, reflection, transmission and homogenized permeability coefficients) have been determined. On the basis of simulation research, prototypes of electromagnetic field two-dimensional absorbers have been manufactured and verified experimentally in the proposed test system. The absorber has been implemented as a component of the low-cost energy harvesting system with a high-frequency rectifier and a voltage multiplier, obtaining usable DC energy from the electromagnetic field in certain frequency bands. The energy efficiency of the system has been determined and the potential application in energy harvesting technology has been assessed.


Introduction
The dynamic development of semiconductor integrated circuit technology, the growing share of wireless powered devices on the industrial and consumer market (in particular in the field of IoT technology, wearable devices and wireless measurement systems) as well as the current pro-ecological policy resulted in the dissemination of electronic devices for which the average demand for electric power does not exceed a single milliwatt. The existing trends in powering these devices are headed towards two directions-improving the efficiency of conventional electricity sources and searching for alternative, renewable sources. The unconventional approach to providing power for devices with low power consumption has resulted in the demand for a new type of technology in which energy is obtained directly from the device's surroundings, known as energy harvesting. The available energy harvesting (EH) solutions are mostly based on commonly used photovoltaic sources.
However, today there is another potential source of energy-the ubiquitous electromagnetic radiation with a frequency band from hundreds of kHz to tens of GHz [1], the so-called "electro-smog". The sources of this radiation are very diverse-from intentionally installed radio and telecommunications transmitters to sources whose emission is an undesirable side effect. The available literature presents various methods of recovering electric energy from the electromagnetic field, mostly based on antenna systems [1,2]. The undoubted disadvantages of antennas, however, are their low geometric compactness (small active surface in relation to the occupied volume) and relatively low absorption (the ability to absorb electromagnetic field energy) [3][4][5][6].
Therefore, it became necessary to search for other types of absorption systems with a much higher absorption efficiency and better application parameters, i.e., greater geometric compactness, the ability to shape frequency characteristics and easy integration on any working surface. A potential solution to the above problems is absorption systems based on resonant metamaterials [7][8][9][10][11][12][13][14][15][16][17], i.e., special electromagnetic structures exhibiting unusual circuit-field properties in certain frequency bands. The miniaturization potential, the possibility of adjusting the geometry of metamaterials to the substrate on which they are to be installed (e.g., a car body, a roof or walls of a building, housing of electronic appliances), much higher absorption, additional beneficial properties and extensive research into the technology [5,8,9,[14][15][16][17][18][19][20][21][22][23][24][25][26][27] (mostly in the microwave spectrum) are important factors supporting the development of effective solutions for obtaining energy from the electromagnetic field using metamaterial absorbers. Table 1 shows a brief summary of selected parameters of EH systems based on rectennas and resonant metamaterials presented in the technical literature in order to compare their electromagnetic energy absorption capabilities. The above-mentioned comparison suggests that the absorption systems based on resonance metamaterials have the efficiency of EM field energy absorption by at least several dozen percent higher than in the case of commonly used patch antennas and as much as several times greater than the conventional dipole-type antenna arrays. At the same time, it is worth noting that the efficiency of electromagnetic field energy absorption by antenna systems is inversely proportional to the bandwidth of the antenna and drastically decreases in the case of broadband antennas due to the reduction in the system quality factor Q. Metamaterial absorbers make it possible to expand the absorption bandwidth with a much smaller relative decrease in the absorption efficiency.

Electromagnetic Smog-Based Energy Harvesting
The available power of electromagnetic smog in modern urban and industrial environments allows it to effectively power many electronic devices whose power demand is not high [21,28,29,37], as shown in Figure 1. In order to verify the energy potential and determine the most advantageous frequency bands for EH purposes, research of the average (for the period of 10 min) spectral power density available in a residential area (urban environment) and an industrialized district of Gdańsk city (Poland) has been performed. Measurements have been made using the D-DOT type wideband field probe and verified with a set of omnidirectional antennas with specific properties at a 50 Ω load. The results are presented in Figures 2 and 3. Based on the supplementary research results, it can be concluded that there are narrow frequency bands, mainly used for telecommunication purposes, where significant levels of power density are available. However, due to the significant dependence of the frequency bands with the greatest availability of power density on the location of the measuring point, the designed energy absorbers should be oriented to a specific application. Hence, for the purposes of this study, several frequency ranges have been arbitrarily designated for which further research has been carried out. In the industrial environment with a significant concentration of various types of telecommunications transmitters and sources of radiated electromagnetic interference (in the frequency range below approx. 30 MHz), a significant local increase in the available power density is visible. Due to the similarity of the characteristics of the given environments in the frequency range above 500 MHz, the scale of Figure 3 has been adjusted accordingly.
This paper presents the results of studies on metamaterial absorber solutions designed to work in the sub-microwave, radio frequency spectrum for which there is a significant potential for energy harvesting, as indicated by environmental measurements. The proposed concept in terms of the energy absorption efficiency, the possibility of the selective formation of frequency characteristics and the application potential significantly exceeds EH systems based on conventional antenna arrays.

Selected Metamaterial Structures
Due to the manufacturing technology and macroscopic properties, metamaterial types can be divided into resonant and non-resonant. Resonant metamaterials are mainly metallic structures placed in a dielectric environment in which characteristic dimension L meets the relationship L < λ/4 [38,39], where λ is the incident wavelength. This condition guarantees that their properties do not only result from the influence of the incident electromagnetic wave but the influence of the secondary wave generated by the resonance of unit cell currents emulating the atomic resonances of the crystal lattice. In the case of non-resonant metamaterials, their characteristic dimension L is much smaller and it is assumed [38,40,41] that it meets the relationship L < (10 −4 to 10 −6 ) λ. The properties of such materials are similar; however, they arise due to the presence of other (microscopic) interactions. This paper focuses on the resonant metamaterials only, due to their higher efficiency [38,40,42] in the considered sub-microwave frequency band, as well as the ease of prototyping such structures for experimental purposes.
The paper presents the result of research on the most representative structure, well described in the literature, with a relatively narrow frequency band, but with a high absorption capability [24,25], i.e., the split ring resonator (SRR) cell proposed originally by [43,44]. The selection of the resonator unit cell topology was based on a review of the literature, and the decisive factor was the possibility of shaping a narrowband, selective absorption characteristic in the low RF range while maintaining the high geometric compactness of the structures. For comparison, the electric component unit cells (ELC) are characterized by much higher resonance frequencies for the same dimensions as for SRR [7,36], while spiral-type [40] cells, which (for comparable cell sizes) exhibit a lower resonance frequency, have a much larger surface filling which makes it impossible to cover transparent surfaces (e.g., windows) with metamaterial matrices. SRR cells represent a compromise between low resonance frequency and high geometry compactness. In the literature, there are examples of the use of SRR cells for EH purposes [16,30,31], but with an emphasis on microwave frequencies (and higher). In this work, models with an absorption band adapted to sub-microwave radio frequencies have been developed, and the geometric and simplified equivalent circuit models are presented in Figure 4. The equivalent circuit model complexity results from the adopted geometrical dimensions of the resonator unit cell, conductive and dielectric materials' parameters and the required accuracy of model representation. It is possible to design multi-mode cells with many resonance frequencies or with more diffused frequency characteristics. The number of resonator modes strictly depends on its electrical dimensions and the coupling between the conducting paths due to the presence of several resonances between the elements of the metamaterial cell (as shown in the equivalent diagram in Figure 4b) with shifted natural frequency points. In this work, the goal was to construct, as far as possible, single-mode metastructure matrices in order to narrow the selective absorption band.  [16,[30][31], but with an emphasis on microwave frequencies (and higher). In this work, models with an absorption band adapted to sub-microwave radio frequencies have been developed, and the geometric and simplified equivalent circuit models are presented in Figure 4. The equivalent circuit model complexity results from the adopted geometrical dimensions of the resonator unit cell, conductive and dielectric materials' parameters and the required accuracy of model representation. It is possible to design multi-mode cells with many resonance frequencies or with more diffused frequency characteristics. The number of resonator modes strictly depends on its electrical dimensions and the coupling between the conducting paths due to the presence of several resonances between the elements of the metamaterial cell (as shown in the equivalent diagram in Figure 4b) with shifted natural frequency points. In this work, the goal was to construct, as far as possible, single-mode metastructure matrices in order to narrow the selective absorption band.

Metamaterial Structure Simulation Research
The purpose of the selected resonant metamaterials' simulation tests is determining their electromagnetic properties in the frequency domain, as well as the suitability of the presented structures for the process of high frequency electromagnetic field energy absorption. Resonators have been modeled in the form of copper tracks with a thickness of 50 μm on a FR-4 (glass-epoxy) dielectric laminate with a thickness of 1.6 mm. The losses of both conductive and dielectric materials have been taken into account in the simulation model. The initial dimensions of the investigated resonator unit cell have been determined on the basis of analytical relationships [5,42], in order to adjust them to an arbitrarily selected frequency of approx. 250 MHz (VHF). The geometric dimensions specified in Figure 4a are presented in Table 2.  Figure 4. Magnetic split ring resonator unit cell (a) geometric model, and (b) simplified resonance circuit model with the optimal incident wave polarization condition: C og -outer ring gap capacitance, C ig -inner ring gap capacitance, C io -inner to outer ring traces capacitance, L i -inner ring selfinductance, L o -outer ring self-inductance, M io -inner to outer ring mutual inductance.

Metamaterial Structure Simulation Research
The purpose of the selected resonant metamaterials' simulation tests is determining their electromagnetic properties in the frequency domain, as well as the suitability of the presented structures for the process of high frequency electromagnetic field energy absorption. Resonators have been modeled in the form of copper tracks with a thickness of 50 µm on a FR-4 (glass-epoxy) dielectric laminate with a thickness of 1.6 mm. The losses of both conductive and dielectric materials have been taken into account in the simulation model. The initial dimensions of the investigated resonator unit cell have been determined on the basis of analytical relationships [5,42], in order to adjust them to an arbitrarily selected frequency of approx. 250 MHz (VHF). The geometric dimensions specified in Figure 4a are presented in Table 2. Geometric models of elementary SRR cells implemented in the CST Studio [45] computing environment and created finite element meshes are shown in Figure 5. As can be seen, metastructure unit cells tuned to the lower radio frequency range may show a small percentage of the coverage of the dielectric substrate with an opaque conducting material (due to the use of narrow path thickness w and a small inter-ring distance d). Therefore, it is possible to apply such structures on transparent functional surfaces, such as windows, without a significant reduction in visibility. In accordance with the purpose of the research, two-dimensional matrices consisting of unit cells of the SRR resonator have been created ( Figure 6). It should be noted that the resonant frequency of the resonator arrays in comparison to a single unit cell (consisting of the same unit cells) is different due to the interaction (mainly capacitive) between the individual matrix cells. For instance, a matrix composed of cells with a natural frequency of about 250 MHz shows a resonance around 350 MHz. This phenomenon can be minimized by increasing the distance between unit cells within one matrix; however, due to the aim to maximize the matrix density, and thus increase the efficiency of electromagnetic field energy absorption, this procedure has been neglected in the further prototyping process. It is also possible to tune matrices' natural frequency by changing the geometrical dimensions of the unit cells (Figure 4a).
For simulation purposes, two representative resonant frequencies have been selected, i.e., 350 MHz (for planar SSR matrix) and 850 MHz (for transverse SSR matrix). The scaled dimensions of the SRR matrices have been calculated and are summarized in Table 3. The simulation models presented in this paper do not take into account the connections between the matrix unit cells, which has a certain impact on the correlation of the parameters of simulation models and real prototypes of metastructure matrices. However, due to the use of some design steps when arranging the strip connections on the actual prototypes of the matrices related to the minimization of parasitic parameters, these discrepancies have been reduced. A series of full-wave parameter calculations as a function of the incident wave frequency have been performed. In each case, the systems were excited by a plane wave with appropriate polarization and a constant power density of 100 mW/cm 2 , e.g., the 350 MHz SSR unit cell was stimulated with a power of 1.78 W. The generation of the discrete finite element mesh has been carried out by adopting first-degree tetrahedral elements. The maximum size of the finite elements for the resonator models has been chosen experimentally, based on the criterion of the minimum increase in computational accuracy with each subsequent step of refining the model mesh. In the case of the space surrounding the resonators, the criterion of minimizing the number of nodes has been used due to the lack of significant influence of this area on the calculation results. In both cases, the necessary condition (resulting from the maximum ratio of the single finite element characteristic size γ to the smallest incident wavelength λ) in Equation (1) has been applied.
The models use boundary conditions appropriate for the propagation path of the electromagnetic wave in the waveguide. The direction of wave propagation has been determined on the basis of the polarization condition of the metamaterial cell incident wave.
On the basis of simulations, frequency characteristics of quantitative parameters characterizing selected structures, i.e., absorbance A, reflectance R, and transmittance T, have been plotted (Figure 7). At the same time, thermal losses in copper and dielectric material at the resonance point of a single SRR cell were determined at the level of approx. 3% of the absorbed power.
In Figure 7c, the presence of multiple absorbance peaks in the SRR two-dimensional transverse matrix can be noticed. This phenomenon results from the existence of a mutual shift of the resonance frequencies of individual stripes constituting the absorption structure. The inner strips of the structure interact more noticeably with the adjacent ones than the outer strips of the matrix. This fact makes it possible to shape the matrix absorption bandwidth at the cost of a decrease in the absorption efficiency (Q factor) of the system. Illustrative views of the electric field distribution around the elementary SRR cell at the resonance point and out of resonance are shown in Figure 8.  For the SRR in resonance, there is a significant concentration of electromagnetic field components in the inter-ring region and at rings' split points. Hence, the highest efficiency in draining electrical energy from the resonator can be expected when the leads are connected precisely at the split point [24].
The characteristics of relative complex permeability coefficients of selected structures (Figure 9) for the representative frequency bands (in which the resonance of metamaterial cells is visible) have been determined using the appropriate simplified transformation algorithm using Matlab [46]. In the literature, there are various approaches to the metamaterials' permeability and permittivity extraction algorithms [19,[47][48][49][50] that differ in the level of complexity and physical interpretation. However, the discussion on the physical interpretation of the µ and ε parameters is not the subject of this application-oriented paper.
The key output parameters of the computer simulation are the coefficients A, R and T, which were determined on the basis of the full-wave power balance in a waveguide.  In order to optimize the impedance matching of the load and metamaterial absorber, the output impedance characteristics of matrices around the resonant point have also been determined. An example of the transverse matrix output impedance Z out characteristics is shown in Figure 10.

Experimental Studies of Metamaterial Structures
For the final verification of metamaterial structures' quantitative parameters and to assess their suitability for use as energy absorbers for electromagnetic field energy harvesting, the prototypes of metamaterial matrices have been designed, manufactured and tested. The aim of the experimental tests was to determine the absorption and transmission coefficients of the metamaterial matrices and compare them to conventional antenna systems. The following assumptions have been made:

•
The experimental research aims to compare the absorption properties of metamaterial matrices with conventional antenna arrays and verifies the suitability of the metasurfaces for energy harvesting purposes. • External interferences and electromagnetic signals present in the test space have a negligible influence on the measurement results due to the significant level of power emitted by the transmitting antenna.

•
The influence of electromagnetic wave reflections in the test space on the measurement results has been ignored.

Materials and Methods
Due to the wide frequency band for which the measurements have been performed (200 MHz to 1500 MHz), two topologies of directional, high-gain transceiving antennas, i.e., the Yagi-Uda (350 MHz to 550 MHz band, with G = 13 dBi directional gain) and the broadband log-periodic antenna LPDA (700 MHz to 2700 MHz, with G = 12 dBi directional gain), were used.
In order to determine the absorption characteristics of antenna systems, a symmetrical arrangement resulting directly from the reciprocity principle was used. This system requires a uniform structure of the transmitting and receiving antenna.
The source antennas and energy receivers (in the form of receiving antennas and metamaterial absorbers) were placed with a constant distance of x = 2 m, which is approximately minimum distance ensuring the location of EM energy receivers in the far zone of the transmitting antenna (in order to ensure the specified, most optimal polarization of the wave incident on the metamaterial matrix, as shown in Figure 4), determined on the basis of the phase criterion [3] and at a constant height of h = 1.5 m on insulating support to eliminate the ground influence on the path of wave propagation. The tests were performed in a laboratory space equipped with shielding systems limiting the influence of external interferences on the measurement results.
About 10 W of power was fed to the transmitting antenna from the broadband RF amplifier. The Rohde & Schwarz SMA100B RF signal generator was used as the signal source. Based on VSWR measurements of the transmitting system, the percentage value of the power transmitted by the antenna was determined as not less than 96.8% of the supplied power. With the use of a broadband D-DOT (SFE3-5G by Montena Technology) probe, the additional measurements of the available surface power density at the location of the receiving systems were performed and determined at the average level of S ref = 82.2 µW/cm 2 (for Yagi-Uda as a transmitter), and S ref = 74.3 µW/cm 2 (for LPDA). The Rohde & Schwarz FSVA3004 spectrum analyzer was used as an energy receiver and power level measuring system.
The drawing of the test system for measuring the absorption properties of classical antennas is presented in Figure 11. The drawing of the test system for measuring the metamaterial matrices' absorbed power characteristics is shown in Figure 12. The use of wideband antennas with an impedance of 50 Ω as energy receivers guaranteed the input of the signal receiver impedance matching. H-1000 coaxial cables with a minimized attenuation coefficient were used as the connecting lines. In the case of metamaterial absorbers, the value of the load impedance, matched most acceptably to all tested matrices, can be defined as the geometric mean value of the output impedance of individual structures. This value was estimated at Z L = 353 Ω on the basis of simulation results. Following the above, in order to connect the matrices with the load, it is unacceptable to use coaxial cables with an impedance of 50 Ω or 75 Ω. Therefore, a cable with the closest impedance value was used, i.e., a symmetrical flat antenna cable with an impedance of 300 Ω. Additionally, in order to load the transmission line properly, an instrumentation amplifier (adapter) with an input impedance of 1 MΩ and an output impedance of 50 Ω was used, the input of which was terminated with a non-inductive 300 Ω resistor. The load mismatch attenuation coefficient was calculated as 0.7%, which is a satisfactory value.
The two-dimensional flat matrices' effective area was assumed as 256 cm 2 , which corresponds to a square with a side length of 16 cm. Flat matrices were manufactured in three variants of resonance frequencies, i.e., 332 MHz, 425 MHz, and 505 MHz. The transverse matrix was designed for a resonant frequency of approx. 816 MHz and made as a four-layer package, with the dimensions of each layer around 20 cm long by 4 cm wide and a distance between the layers of 1.5 cm. According to the observation of the highest concentration of the electromagnetic field in the split points of the SRR, unit cells were connected to the matrix at this point. Individual cells within one matrix were connected in series and parallel, while the outputs from the matrices were made in the form of high-frequency coaxial connectors. Connections between unit cells were designed as wide strips of copper to minimize the parasitic inductances. The printed circuit boards were made as double-sided and the connections on the bottom layer were shifted with respect to the connections on the top layer of the board in order to minimize parasitic capacitances. In order to increase the thickness of the printed paths to the level of about 50 µm, tin electroplating was used. The views of the metamaterial matrices prototypes are shown in Figure 13.

Experimental Results of the System Absorption Efficiency
In order to compare the efficiency of electromagnetic field energy absorption by conventional antenna systems and metamaterial absorbers in the radio frequency range, the frequency characteristics of the receiver power have been measured in the test system as shown in Figures 11 and 12. The experimental results have been presented in the form of relative power absorption efficiency characteristics of the considered systems, which have been defined as (Equation (2)): where: P A -power absorbed by the system [µW], A eff -effective surface area of the absorption system [cm 2 ], S ref -reference available surface power density at the location of the absorbing system [µW/cm 2 ]. The effective absorption areas A eff of the Yagi-Uda and LPDA antennas have also been determined based on the approximate relationship (Equation (3)), i.e., 2618 cm 2 and 560 cm 2 , respectively. The effective surface of metamaterial absorbers used for calculations was the real surface of the prepared matrices.  Figure 14 presents the frequency characteristics of the electromagnetic field energy absorption efficiency of conventional antenna systems used as energy receivers along with the maximum value of absorbed power and the determined maximum efficiency. Figure 15 presents the frequency characteristics of the efficiency of electromagnetic field energy absorption in the presented metamaterial matrices along with the maximum value of absorbed power and determined maximum efficiency at the resonance point.   Characteristic parameters of the tested absorption systems are shown in Table 4. The proposed topology of resonators is characterized by very narrow resonance bands, which makes it possible to design systems with high absorption for frequency ranges not useful for the powered device (e.g., electromagnetic interference) without the absorption of useful signals (e.g., communication). To shape the selectivity of the absorbers, a layered combination of two-dimensional flat matrices has been proposed, where for each layer it is possible to use an absorber with a different resonance frequency to cover the desired absorption band. The measurement results of this type of solution are shown in Figure 16. As predicted, the frequency response is a superposition of the three individual metamaterial matrices. The peak value of the power absorbed decreased to the value of 15.7 mW, which in terms of percentage absorption efficiency is 75.1%. The visible decrease in the efficiency of energy absorption due to the interactions between individual matrices has been found. Another possibility of shaping the frequency characteristics of metamaterial absorbers is to use the multi-mode structure topologies [39], e.g., based on fractal geometries. However, due to the resonant quality factor Q reduction in these types of solutions, their absorption efficiency reaches lower values.

Metamaterials in the Energy Harvesting System
In order to verify the suitability of metamaterial absorbers for energy harvesting purposes, a prototype of a system for converting energy from a high-frequency electromagnetic field to usable DC has been developed. Due to the low output voltages of the metamaterial absorbers, a symmetrical rectifier with a voltage multiplier has been used. The presented EH system has been tested in the test configuration shown in Figure 17 to determine the frequency characteristics, energy conversion efficiency and variability of output parameters as a function of load. The rectifier circuit has been built with the use of high-frequency Schottky diodes BAT15-04R and ceramic capacitors with capacitance values C 1 = 1 nF and C 2 = 10 nF, reduced ESR, and increased tolerance of 0.2%. In order to determine the characteristics of the output voltage and the power generated by the system as a function of the load, the values of the load resistance R o were changed in the range from 2 kΩ to 16 kΩ for the first resonant frequency (Figure 16) of the absorber. Based on the load characteristics, the optimal value of the load resistance has been determined at the level of 5.6 kΩ. For the optimal value of the load resistance of the system, the voltage and the output power value of the EH system as a function of the frequency of the input signal have been measured. The measurement results are presented in Figures 18 and 19.  The maximum value of the obtained power has been determined at the level of approx. 10.9 mW at a voltage of 9.1 V. Assuming the available power density from the transmitting antenna, the maximum value of energy conversion efficiency of the entire EH system has been calculated at the level of approx. 56% for the frequency of 332 MHz. The high efficiency drop of the EH system is caused by the adopted low-cost rectifier with a voltage multiplier circuit. A further increase in the system efficiency is possible through the use of low-loss switching power electronic elements, made, for example, of gallium nitride.

Conclusions
Based on the conducted research, it has been proven that high-frequency electromagnetic field energy absorbers based on resonance metamaterials constitute a competitive alternative to conventional antenna systems (with comparable sizes) as transducers in EH systems achieving efficiency exceeding 90% (more than five times higher than the conventional dipole-type antenna arrays). It should also be noted that extending the absorption band in the case of antenna systems comes at the significant expense of the maximum absorption efficiency. Metamaterial absorbers are more resistant to this phenomenon.
Effective shaping of the metamaterials' energy absorption frequency bands can be achieved by the superposition of many matrices with different resonance frequencies. This type of solution allows the formation of the absorption band of the EH system selectively. The absorption band can also be formed through the use of multi-mode metamaterial absorbers, e.g., fractal-based structures. The proposed EH system with a metamaterial absorber and a symmetrical rectifier with a voltage multiplier meets the desired expectations. The efficiency of energy transmission to the receiver or potential energy storage, taking into account losses at each stage of energy conversion, reaches the value of approximately 56% and is at a satisfactory level, comparable or greater than the effectiveness of solutions available in the literature. It is possible to further increase the efficiency of the system by using a rectifier and a voltage multiplier with lower switching losses.
The practical application of meta-surface structures is easier compared to classic antenna systems due to the greater flexibility in the geometry selection of the two-dimensional matrices. Many types of surfaces can be covered with metamaterial absorbers, including housings of electronic devices. Metamaterial absorbers, due to the possibility of selectively shaping the frequency absorption characteristics, will not attenuate the telecommunications signals necessary for the proper operation of devices. In particular, in the case of narrowstrip resonators such as the SRR (especially for the radio and microwave frequency range), it is possible to coat transparent utility surfaces such as glass car windows or residential windows with metastructures without a significant reduction in visibility. The proposed application enables powering electronic systems with low-power demand using only the energy harvesting technique or in cooperation with other energy sources. It should be noted, however, that due to the availability of sufficient electromagnetic field power only in the environments with a high level of electromagnetic field intensity, the use of this type of technology is so far limited mainly to urban and industrial environments.
Resonant metamaterials can also be used as high-efficiency couplers and receivers in the microwave wireless power transfer technique with the additional possibility of directing the electromagnetic radiation beam (metamaterial lenses) or highly effective and selective electromagnetic shielding systems for electromagnetic compatibility or military purposes.

Conflicts of Interest:
The author declares no conflict of interest.