## Abstract

Brewster effect has attracted extensive interests through microwave to optical regime. However, previous work mainly focused on the cases with single linear polarization. We demonstrate in this paper an equal Brewster’s angle for both TE (transverse-electric) and TM (transverse-magnetic) waves. Tunable Brewster effect is furthermore achieved. For practical implementation of this effect, we propose an anisotropic metagrating composed of modified split ring resonators (MSRR) as a proof-of-concept in X-band. Measurement results of the manufactured prototypes validate our method. The proposed metagrating featuring Brewster’s angle manipulation potentially paves the way for many applications throughout the spectrum, e.g., angular selectivity.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Well-studied since Sir David Brewster and Malus, traditional Brewster effect has various applications in both optical and microwave regimes: polarizer [1], Brewster angle microscopy (BAM) [2], dielectric properties quantifying [3], and extraordinary optical transmission (EOT) [4,5]. These classical applications are yet limited to transverse-magnetic (TM, or *p*-polarized) electromagnetic (EM) waves. For example, in the most common laser cavity, where the Brewster window acts as a polarizer, the angular selectivity is always entangled with TM polarization. This restriction is caused by the fact that the classic Brewster angle for transverse-electric (TE, or *s*-polarized) waves only exists when the waves impinge on the surface of magnetic media [1], but the magnetic response in natural media is usually tenuous, especially in the optical regime. Therefore, it is of practical significance to find a medium with an identical Brewster’s angle that is not sensitive to the polarization. Another possible application of such a medium is an angle-filter promoting the signal-noise-ratio (SNR) in the receivers of the oriented wireless communication, especially for unknown polarization cases. More importantly, it would inspire many novel applications spanning the entire spectrum, e.g., information encryption with geometric angles, and efficient energy harvesting and conversion [6].

The confines abovementioned have been overcome with the advances in metamaterials and metasurfaces. The generalized Brewster effect—realizations of Brewster’s angle for TE waves—were recently reported in these artificially composite media. For the visible frequencies, Ryosuke et al. realized non-reflection only for TE waves at a certain angle (Brewster’s angle), using stratified metal-dielectric metamaterial [7]. Paniagua-Domínguez et al. applied the first Kerker’s condition [8] and the Mie theory [9] to achieve scattering cancellation, realizing the generalized Brewster effect in all-dielectric metasurfaces (silicon nanodisks) [10]. Similarly, high-refractive index (HRI) nanophotonic structures were properly described by coupled electric and magnetic dipole (CEMD) analytical formulation, and were further numerically confirmed to achieve the generalized Brewster effect [11]. Besides, the optical angular selectivity was proposed by Y. Shen [12,13], where the broadband total transmission at a narrow-angle range was enabled by hetero-structured photonic crystals.

The reported meta-atoms in the optical or terahertz regimes, e.g., all-dielectric disks, spheres, cylinders, stacks, graphene, or other Mie scatters and gratings [10-18], feature geometrically simple structures and thus their electrical and magnetic polarizabilities can be described with analytical expressions, which in turn simplifies yet constrains the applicable range of the design method. Moreover, these three-dimensional structures are inconvenient to be miniaturized or integrated if applied to lower frequencies. Thus, the generalized Brewster effect with classical microwave fabrication techniques remains relatively elusive. Yasuhiro Tamayama experimentally observed TE Brewster effect in a bulk metamaterial [19], and theoretically proposed bi-anisotropic metafilms with 45° Brewster’s angle for only TE waves [20]. Recently, Guillaume Lavigne and Christophe Caloz extended the Brewster effect to arbitrary angles and polarizations in bi-anisotropic metasurfaces [21] and analyzed it with generalized sheet transition conditions (GSTC) [22,23]. Unfortunately, their proposed metasurface consisted of purely bi-anisotropic scattering particle tend to be difficult to implement [24]. Despite the provided possible solutions, onerous susceptibility derivation could not be avoided [25].

In this paper, we generalize angularly tunable Brewster effect for both TE and TM polarized electromagnetic waves impinging on our metagrating. Two main advances are made in our work in comparison with the previous research in microwave and optical regime. For one thing, the aforementioned polarization restriction is overcome by a properly designed thin metasurface with identical Brewster’s angle for TE and TM waves. For the other, the angular manipulation can be realized through changing either the grating space or the grating tilt angle. These two key phenomena can be explained by the anisotropic EM characteristic of the proposed metagrating. Therefore, we will first analytically demonstrate the derivation of the generalized Brewster effect conditions in anisotropic media. In accordance to these conditions, we then adopt the modified split ring resonator (MSRR) structure [26] with simplified EM parameters. Based on its simulated scattering parameters, the permittivity and permeability tensors will be retrieved to confirm the Brewster effect by eliminating the possibility of Fabry-Pérot-type of total transmission. The circularly polarized waves are also applied in the simulation and experiment to substantiate polarization arbitrariness. In the last section, the manipulation of the Brewster’s angle is demonstrated. The proposed two facile manipulating methods make it possible to achieve applications involving angular selectivity. Experimental results of the manufactured prototypes further validate our design method.

## 2. Brewster effect in anisotropic media

Considering the most general media, one can describe their EM characteristics with four dyadic tensors: $\overline{\overline{\epsilon}}$, $\overline{\overline{\mu}}$, $\overline{\overline{\xi}}$, $\overline{\overline{\zeta}}$. The constitutive relations should be written as

where*c*

_{0},

*ε*

_{0},

*μ*

_{0}denote the speed of light, permittivity, and permeability in the vacuum. In a time-harmonic field, Maxwell’s equations can be applied to get the wave equation for the electric field

**k**must satisfy Eq. (3) with a non-trivial electric field, which leads to

The Eq. (4) is one of the dispersion equations for any plane wave in the most general form. However, many hidden difficulties make these equations unpractical to apply. Chief among them is that the plane wave is hard to decompose to TE and TM waves in the most bi-anisotropic cases. In addition, a worse scenario lies in four undetermined dyadic tensors with 36 unknown parameters. To simplify this problem for further design, we assume that

*θ*), Brewster effect happens, so there should be no reflection. Regarding the

_{B}*xOz*plane as the incident plane, the dispersion equation for the TE polarized wave takes the form as the follows,

*k*= |

_{x}**k**|sin

*θ*and

_{t}*k*= |

_{z}**k**|cos

*θ*, and

_{t}*θ*is the refraction angle. Meanwhile, the no-reflection assumption can be used to simplify the EM boundary conditions, by which one gets

_{t}## 3. Practical realization

#### 3.1 Unit cell design of the metagrating

In the microwave regime, any non-magnetic dielectric generates the non-unity relative permittivity (*ε _{xx}*,

*ε*,

_{yy}*ε*≠ 1), which guarantees the prevalence of the Brewster effect for TM polarized waves. As for the TE Brewster’s angle, Eq. (8) requires at least non-unity relative permeability (

_{zz}*μ*≠ 1). Thus, only in the media with permittivity and permeability both non-unity and properly designed, can the Brewster’s angle exist for arbitrary polarized waves.

_{xx}μ_{zz}Looking back to the advance in metamaterial and metasurface, the classic artificially composed medium that fits this condition is called double negative material (DNG) proposed by D. R. Smith [27,28]. We similarly take the SRR structure in [19,20] and [27,28] into account (as shown in Fig. 1(a)). Unfortunately, the traditional SRR unavoidably behaves as a bi-anisotropic particle with its electric and magnetic moments shown as [26]

According to Eq. (10), an external magnetic field along the*x*axis would generate the electric moment along the

*y*axis. As shown in Fig. 1(c), this phenomenon can be verified by the normalized surface current distribution extracted from the CST Microwave Studio simulation. The currents on the inner and outer rings are in the opposite directions, accumulating the net electric charges in the way shown in Fig. 1(c) and producing the electric moment

**p**

*sketched in Fig. 1(a).*

_{m}Therefore, we then replace one of the rings by another one on the opposite side of the substrate with the split orientating to the opposite direction [26]. Our modified SRR (MSRR) unit cell and its simulated current distribution are shown in Figs. 1(b) and 1(d). The EM behaviors of the MSRR should be similar to those of traditional SRR, e.g., the negligible non-diagonal components of the electric permittivity and the magnetic permeability due to high polarization isolations, i.e., *ε _{ij}*,

*μ*with

_{ij}*i*≠

*j*. They behave disparately, however, with regard to the bi-anisotropic effect. In fact, as shown in Fig. 1(d), the currents on both rings afflux in the same direction. The oppositely-placed splits and exactly equal dimensions enable the cancellation of electric moments along the

*y*axis caused by the external magnetic field (

**p**

*= −*

_{mf}**p**

*), in that the electric charges on the front ring are just opposite to those on the back. Hence, the cross-coupling susceptibility*

_{mb}*χ*in Eq. (10) and Eq. (11) would be minimized to zero. To put it another way, the magnetic moment is purely contributed by the external magnetic field, while the electric moment totally comes from the external electric field.

^{em}#### 3.2 Brewster effect for arbitrary polarized EM waves

After cancelling the cross-coupling susceptibility, we can reasonably apply Eq. (8) in our design of metagrating and verification for the generalized Brewster effect. The proposed metagrating is sketched in Fig. 2.

When the TE polarized wave impinges on the metagrating, the electric and magnetic field can be illustrated as the red arrows in Fig. 2. The electric field along the *y* axis *E _{y}* and the existence of the substrate are able to produce non-unity

*ε*, while the non-unity

_{yy}*μ*can be realized by magnetic field along the

_{xx}*x*axis

*H*. Meanwhile, the periodical repetition of the unit cells provides more freedom for the parameters in Eq. (8). In this case, Brewster’s angle for the TE polarized wave can be achieved by our metagrating, and it is also true for the TM case with similar analysis. Besides, it is clear that the dimension of the ring and the space between the gratings could influence the EM response of the metagrating, and then make it possible to manipulate the Brewster’s angle.

_{x}Here, we choose F4B as the substrate with relative permittivity *ε _{r}* = 2.2 and its thickness

*t*= 0.8 mm. In addition, we set the grating space

_{s}*P*= 8 mm, unit cell dimension

_{x}*P*=

_{y}*d*= 20 mm, the outer radius of the ring

*R*= 9 mm, the inner radius

_{out}*R*= 3 mm, and the split width

_{in}*s*= 3.5 mm.

As shown in Fig. 3(a), when the TE or TM wave in the free space impinges on the metagrating, the reflection index can be retrieved from the simulated scattering parameters (*S _{ij}*,

*i, j*= 1, 2), which can be written as

*z*is the generalized equivalent impedance of the metagrating and derives from [29]

Based on the simulated scattering parameters, the retrieved reflection indices for TE and TM polarized waves are shown in Figs. 3(b) and 3(c), respectively. The Brewster’s angles are further outlined by the white dash lines. It is interesting to note that the critical angle, above which the waves would be totally reflected, can also be witnessed in Fig. 3. Yet we care much more about the case when TE and TM polarized waves have the same Brewster’s angle, as marked by the white crosses in Figs. 3(b) and 3(c). It can be witnessed that TE and TM polarized waves share 45° Brewster’s angle at 10.3 GHz.

However, the Brewster effect should be verified by retrieving the EM parameters of the metagrating and excluding the Fabry-Pérot type of total transmission. Retrieving from the simulated scattering parameters, one can usually calculate the equivalent EM parameters of a structure, but it should be noticed that the traditional retrieval method is not accurate here due to the oblique incidence—the retrieved EM parameters inherently vary with the incidence angle [30]. Hence, we apply a modified retrieval method in which scattering parameters in the normal incidence case are taken as assisted conditions. As shown in Fig. 3(a), one can apply the EM boundary conditions here and readily write the scattering parameters with the reflection index, which reads,

where*R*can be represented by Eq. (12) and

*Q*= exp(

*jn*

_{TE}(

*θ*)

_{t}*k*

_{0}

*d*cos

*θ*).

_{t}*n*

_{TE}(

*θ*) is the refraction index at the angle of

_{t}*θ*. To conveniently utilize the available retrieval formulation developed in the normal incidence case, we define the generalized refraction index

_{t}*n*

_{TE}

*’*as follows,In this form, it can readily be retrieved by Eq. (14)-(15).

The Snell’s law can be further applied as

Considering the TE polarized wave case,*n*

_{TE}(

*θ*) also conforms to

_{t}*θ*= 0),

_{t}*ε*can be solved from Eq. (18) and Eq. (20), which reads,It should be noted that to determine anisotropic components of these electromagnetic parameters is not straightforward since the number of unknown parameters (six in total) is larger than the conditions even with two polarizations applied (four in total). Nevertheless, since the electric field of TE waves is always along the

_{yy}*y*axis as shown in Fig. 2,

*ε*should remain the same value despite of the variation of incident angle. It is also true for

_{yy}*μ*in the TM case. Such conditions help us find a convenient and feasible approach to determine all the anisotropic components uniquely. Therefore,

_{yy}*ε*calculated by Eq. (21) can be then used to determine the other two EM parameters through Eq. (17)-(20) with retrieved

_{yy}*n*

_{TE}

*’*in the TE polarized case. They are,

*j*2

*n*(

*θ*)

_{t}*k*

_{0}

*d*cos

*θ*)| ≠ 0.

_{t}To substantiate that 45°-Brewster’s angle applies to arbitrary polarized EM waves, we take TE, TM, and circularly polarized waves as excitations, separately. The experiment setup photograph and schematic are shown in Figs. 4(a) and 4(b), respectively. In order to measure the return loss (*S*_{11}) of the metagrating, two standard gain horn antennas (XB-GH90-20N with the BJ-100 waveguide) are worked as the source (Tx) and the receiver (Rx) in our experiments. In addition, the zoom-in photographs of the experimental prototype of the metagrating (24 × 9 units, around 6.32*λ*_{0} × 6.18*λ*_{0} at 10.3 GHz) and two sides (marked as A and B separately) of its single strip are demonstrated in Figs. 4(c) and 4(d).

Adopting the experiment setup abovementioned, the simulated and measured return loss (*S*_{11}) are illustrated in Figs. 4(e) and 4(f), where good agreement can be found between the simulation and the experiment results. When the circularly polarized waves are implemented with cone horns, the measured axial ratio of transmitted waves is as low as 0.78 dB (after calibration) at 10.3 GHz and 45° incident angle, which is also well consistent with the simulated 0.11 dB.

#### 3.3 Brewster’s Angle Manipulation

As aforementioned, many structural parameters of the metagrating play a notable role in determining Brewster’s angle and its corresponding frequency. Specifically, the average radius of MSRR would mainly determine at which frequency the Brewster effect occurs. On the other hand, the slot width (*s*) and the grating space (*P _{x}*) can both be changed to manipulate the Brewster’s angle. Yet once the MSRR structure is printed, the dimension of the metagrating units will be fixed, so the only possible way to manipulate the Brewster’s angle is to change the grating space or its tilt angle. The metagratings with different grating space (as shown in Fig. 5(a)) are manufactured and applied to our experiments. In Figs. 5 (b) and 5(c), we demonstrate the simulated and measured TE and TM Brewster’s angles at different frequencies when the transvers space

*P*varies from 4 mm to 10 mm. Both of them cover a broad angle range (TE: 15~80°, TM: 20~80°) in the X-band, which proves the efficiency of this manipulating method.

_{x}The intrinsic reason for the aforementioned manipulation is to change the equivalent EM parameters of the meta-grating by tuning the space between the strips. From the microscopic perspective, the underlying mechanism is based on coupled electric/magnetic dipole theories. Accordingly, we could also mechanically achieve the similar effect by adjusting the tilt angle of strips (Θ) as shown in Fig. 5(d). The simulated variation of Brewster’s angle at different tilt angles for several frequencies is shown in Figs. 5(e) and 5(f). By tuning the tilt angle of the metagrating, the Brewster’s angle covers the range from 0° through 80°. The tunability of EM properties achieved by tilting the meta-atoms’ axes with respect to the metasurface plane has also been reported by Sanchez-Gil’s group [31]. Overall, these two facile manipulation methods for Brewster effect are indeed of practical application.

## 4. Conclusion

To summarize, we have both theoretically investigated and experimentally achieved the generalized Brewster effect in anisotropic media. To this aim, the analytical formulations of the Brewster’s angles for both TE and TM polarized electromagnetic waves have been firstly derived. Although these two formulations have been properly simplified from the most general forms, they are still suitable for many occasions throughout all frequency bands. Most importantly, their explicit forms make it possible to instruct the practical designing process, especially when more sophisticated structures, e.g., metal and dielectric combined, are required.

Based on the theoretical analysis, a metagrating composed of modified split ring resonator (MSRR) has been proposed, featuring the identical and tunable Brewster’s angle for arbitrary polarized electromagnetic waves. To justify its commitment to the mentioned simplification, the discrepancies between the MSRR and the traditional split ring resonator (SRR) have been circumspectly analyzed. As a case study, 45°-Brewster’s-angle metagrating has been verified by full wave simulation, experiments and numerical retrieval of electromagnetic parameters. Furthermore, the manipulation of the Brewster’s angle for both two linear polarized waves have also been proved efficient enough.

All the simulated, measured and relating retrieved results are satisfactorily congruous with one another. The proof-of-concept prototype of our metagrating validates our design, which could illuminate further research on Brewster-effect-related topics. For instance, two-dimensional (2D) grating might also be able to achieve generalized Brewster effect with broader bandwidth. Once achieved, it will release the frequency-polarization entanglement in the traditional Brewster effect, as well as be easier to be miniaturized and integrated on a chip in possibly any frequency regime.

## Funding

National Natural Science Foundation of China (NSFC) (61671178, 61301013).

## References

**1. **E. Hecht, *Optics* (Addison-Wesley, 2002).

**2. **W. Daear, M. Mahadeo, and E. J. Prenner, “Applications of Brewster angle microscopy from biological materials to biological systems,” Biochim Biophys Acta Biomembr **1859**(10), 1749–1766 (2017). [CrossRef] [PubMed]

**3. **L. I. Thomson, G. R. Osinski, and W. H. Pollard, “Application of the Brewster angle to quantify the dielectric properties of ground ice formations,” J. Appl. Geophys. **99**, 12–17 (2013). [CrossRef]

**4. **A. Alù, G. D’Aguanno, N. Mattiucci, and M. J. Bloemer, “Plasmonic Brewster Angle: Broadband Extraordinary Transmission through Optical Gratings,” Phys. Rev. Lett. **106**(12), 123902 (2011). [CrossRef] [PubMed]

**5. **N. Aközbek, N. Mattiucci, D. de Ceglia, R. Trimm, A. Alù, G. D’Aguanno, M. A. Vincenti, M. Scalora, and M. J. Bloemer, “Experimental demonstration of plasmonic Brewster angle extraordinary transmission through extreme subwavelength slit arrays in the microwave,” Phys. Rev. B Condens. Matter Mater. Phys. **85**(20), 205430 (2012). [CrossRef]

**6. **P. Bermel, M. Ghebrebrhan, W. Chan, Y. X. Yeng, M. Araghchini, R. Hamam, C. H. Marton, K. F. Jensen, M. Soljačić, J. D. Joannopoulos, S. G. Johnson, and I. Celanovic, “Design and global optimization of high-efficiency thermophotovoltaic systems,” Opt. Express **18**(S3Suppl 3), A314–A334 (2010). [CrossRef] [PubMed]

**7. **R. Watanabe, M. Iwanaga, and T. Ishihara, “s‐polarization Brewster’s angle of stratified metal–dielectric metamaterial in optical regime,” Phys. Status Solidi, B Basic Res. **245**(12), 2696–2701 (2008). [CrossRef]

**8. **M. Kerker, D.-S. Wang, and C. L. Giles, “Electromagnetic scattering by magnetic spheres,” J. Opt. Soc. Am. **73**(6), 765–767 (1983). [CrossRef]

**9. **C. F. Bohren and D. R. Huffman, *Absorption and scattering of light by small particles* (John Wiley & Sons, 2008).

**10. **R. Paniagua-Domínguez, Y. F. Yu, A. E. Miroshnichenko, L. A. Krivitsky, Y. H. Fu, V. Valuckas, L. Gonzaga, Y. T. Toh, A. Y. Kay, B. Luk’yanchuk, and A. I. Kuznetsov, “Generalized Brewster effect in dielectric metasurfaces,” Nat. Commun. **7**(1), 10362 (2016). [CrossRef] [PubMed]

**11. **D. R. Abujetas, J. A. Sánchez-Gil, and J. J. Sáenz, “Generalized Brewster effect in high-refractive-index nanorod-based metasurfaces,” Opt. Express **26**(24), 31523–31541 (2018). [CrossRef] [PubMed]

**12. **Y. Shen, D. Ye, I. Celanovic, S. G. Johnson, J. D. Joannopoulos, and M. Soljačić, “Optical broadband angular selectivity,” Science **343**(6178), 1499–1501 (2014). [CrossRef] [PubMed]

**13. **Y. Shen, D. Ye, L. Wang, I. Celanovic, L. Ran, J. D. Joannopoulos, and M. Soljačić, “Metamaterial broadband angular selectivity,” Phys. Rev. B Condens. Matter Mater. Phys. **90**(12), 125422 (2014). [CrossRef]

**14. **V. Neder, Y. Ra’di, A. Alù, and A. Polman, “Combined Metagratings for Efficient Broad-Angle Scattering Metasurface,” ACS Photonics **6**(4), 1010-1017 (2019).

**15. **Y. Ra’di, D. L. Sounas, and A. Alù, “Metagratings: Beyond the Limits of Graded Metasurfaces for Wave Front Control,” Phys. Rev. Lett. **119**(6), 067404 (2017). [CrossRef] [PubMed]

**16. **A. Cordaro, J. van de Groep, S. Raza, E. F. Pecora, F. Priolo, and M. L. Brongersma, “Antireflection High-Index Metasurfaces Combining Mie and Fabry-Pérot Resonances,” ACS Photonics **6**(2), 453–459 (2019). [CrossRef]

**17. **M. Nieto-Vesperinas, R. Gomez-Medina, and J. J. Saenz, “Angle-suppressed scattering and optical forces on submicrometer dielectric particles,” J. Opt. Soc. Am. A **28**(1), 54–60 (2011). [CrossRef] [PubMed]

**18. **X. Lin, Y. Shen, I. Kaminer, H. Chen, and M. Soljačić, “Transverse-electric Brewster effect enabled by nonmagnetic two-dimensional materials,” Phys. Rev. A (Coll. Park) **94**(2), 023836 (2016). [CrossRef]

**19. **Y. Tamayama, T. Nakanishi, K. Sugiyama, and M. Kitano, “Observation of Brewster’s effect for transverse-electric electromagnetic waves in metamaterials: Experiment and theory,” Phys. Rev. B Condens. Matter Mater. Phys. **73**(19), 193104 (2006). [CrossRef]

**20. **Y. Tamayama, “Brewster effect in metafilms composed of bi-anisotropic split-ring resonators,” Opt. Lett. **40**(7), 1382–1385 (2015). [CrossRef] [PubMed]

**21. **Guillaume Lavigne and Christophe Caloz, “Extending the Brewster effect to arbitrary angle and polarization using bianisotropic metasurfaces,” in 2018 IEEE International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting (IEEE, 2018), pp. 771–772.

**22. **E. F. Kuester, M. A. Mohamed, M. Piket-May, and C. L. Holloway, “Averaged transition conditions for electromagnetic fields at a metafilm,” IEEE Trans. Antenn. Propag. **51**(10), 2641–2651 (2003). [CrossRef]

**23. **C. L. Holloway, E. F. Kuester, and A. Dienstfrey, “Characterizing Metasurfaces/Metafilms: The Connection Between Surface Susceptibilities and Effective Material Properties,” IEEE Antennas Wirel. Propag. Lett. **10**, 1507–1511 (2011). [CrossRef]

**24. **M. Albooyeh, V. S. Asadchy, R. Alaee, S. M. Hashemi, M. Yazdi, M. S. Mirmoosa, C. Rockstuhl, C. R. Simovski, and S. A. Tretyakov, “Purely bianisotropic scatterers,” Phys. Rev. B **94**(24), 245428 (2016). [CrossRef]

**25. **G. Lavigne, K. Achouri, V. S. Asadchy, S. A. Tretyakov, and C. Caloz, “Susceptibility Derivation and Experimental Demonstration of Refracting Metasurfaces Without Spurious Diffraction,” IEEE Antennas Wirel. Propag. Lett. **66**(3), 1321–1330 (2018). [CrossRef]

**26. **R. Marqués, F. Medina, and R. Rafii-El-Idrissi, “Role of bianisotropy in negative permeability and left-handed metamaterials,” Phys. Rev. B Condens. Matter Mater. Phys. **65**(14), 144440 (2002). [CrossRef]

**27. **D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite Medium with Simultaneously Negative Permeability and Permittivity,” Phys. Rev. Lett. **84**(18), 4184–4187 (2000). [CrossRef] [PubMed]

**28. **R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**(5514), 77–79 (2001). [CrossRef] [PubMed]

**29. **D. R. Smith, D. C. Vier, T. Koschny, and C. M. Soukoulis, “Electromagnetic parameter retrieval from inhomogeneous metamaterials,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. **71**(33 Pt 2B), 036617 (2005). [CrossRef] [PubMed]

**30. **J. Qi, H. Kettunen, H. Wallen, and A. Sihvola, “Compensation of Fabry–Pérot Resonances in Homogenization of Dielectric Composites,” IEEE Antennas Wirel. Propag. Lett. **9**, 1057–1060 (2010). [CrossRef]

**31. **R.D. Abujetas, Á. Barreda, F. Moreno, J. J. Sáenz, A. Litman, J. Geffrin, and J. A. Sánchez-Gil, “Brewster quasi bound states in the continuum in all-dielectric metasurfaces from single magnetic-dipole resonance meta-atoms,” arXiv:1902.07148 (2019).