Abstract
Bi-anisotropic optical metamaterials are playing an increasingly important role in current wave-functional metamaterials and topological photonics due to their extra degree of freedom in addition to the permittivity and permeability. In this work, we derived the closed-form expressions for effective constitutive parameters of 2-dimensional (2D) bi-anisotropic metamaterials whose chirality tensors can possess both diagonal and off-diagonal components in the long-wavelength limit based on the Mie theory. Our formulas can be regarded as an extension of the Maxwell-Garnet formula to 2D bi-anisotriopic metamaterials and are verified through full wave numerical simulations. These closed-form formulas will benefit the design and analysis of the optical properties of 2D bi-anisotriopic metamaterials.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Recently, a great deal of attention has been paid to bi-anisotropic metamaterials due to their intriguing properties in wave manipulations, such as negative refraction [1–4], strong optical activity [5–7] and circular dichroism [8–10]. Bi-anisotropic metamaterials have also been widely noticed as a fruitful platform for realizing nontrivial topological phenomena in photonics [11–14], due to the inversion symmetry breaking and strong spin-orbital coupling effects.
In the framework of effective medium theory (EMT) [15], the scattering properties of metamaterials as well as the propagating properties of light inside metamaterials are sufficiently described by the effective constitutive parameters. As a result, acquiring these effective constitutive parameters correctly and efficiently is the main concern in the EMT. There have been lots of models and methods employed for determining the effective constitutive parameters, either analytically or numerically [16–23]. For the cases of achiral inclusions, the well-known Maxwell-Garnet (MG) formula is widely used for obtaining the effective constitutive parameters in the long-wavelength limit due to its accuracy and efficiency [24]. Although the MG formula has already been extended to three-dimensional (3D) chiral metamaterials composed of isotropic spherical chiral particles fora longtime[25,26], the more nontrivial extension tobi-anisotropic cylinders in two-dimensional (2D) system (the system is translation invariant along out-of-plane direction) is still lacking. Unlike the 3D isotropic chiral metamaterials, the 2D bi-anisotropic metamaterials is inherently anisotropic due to the asymmetry between the in-plane and out-of-plane directions even if the inclusions are isotropic. This structural anisotropy as well as the material anisotropy of the inclusion itself bring us great difficulty to obtain a close-form formula when applying the EMT and conventionally the effective constitutive parameters are only obtainable by numerical retrieval [27–33].
In this work, we derived for the first time the analytic formulas of the effective constitutive parameters for 2D bi-anisotropic metamaterials whose chirality tensors have both diagonal and off-diagonal terms, in the long-wavelength limit through the use of Mie theory. Two sets of formulas were derived, one is based on the Mie coefficients of the inclusions which are valid even for inclusions having complex underlying structures, the other is based on the constitutive parameters of the homogeneous inclusions. The latter set of formulas are directly reduced to the traditional MG formula for 2D systems in the limit of vanishing chirality. Consequently, the formulas can be regarded as an extension of the MG formula to the 2D bi-anisotriopic metamaterials. The validities of our analytical formulas are confirmed by the full wave numerical simulations. Our work will benefit the investigation of wave manipulating properties of the bi-anisotropic metamaterials, and will also be conductive to the design of metamaterials for practical uses.
2. Helmhotz equation for the bi-anisotropic medium
The geometrical sketch of the bianisotropic metamaterial considered in this study are schematically shown in Fig. 1, where identical bi-anisotropic inclusions are arranged into a regular or random lattice in the xy plane (in-plane) embedded in air with permittivity ε0 and permeability µ0. The inclusions are invariant along z (out-of-plane) direction, and only the in-plane propagating electromagnetic waves (the out-of-plane wavenumber kz = 0) are considered throughout this paper. Therefore, the electromagnetic fields are independent of z.
The constitutive relations of the bi-anisotropic medium are given by [34]
where c is the light speed in vacuum, and ,, are the relative permittivity, permeability and chirality tensors. For the sake of simplicity, in the following we will set ϵ0 = µ0 = c = 1. We assume that the medium is reciprocal and isotropic in the xy plane. Therefore the permittivity and permeability tensors can be written as where εt,εz (µt,µz) are the in-plane and out-of-plane relative permittivities (permeabilities). The chirality tensor considered here has both diagonal and off-diagonal terms and is written asWe will show that the electromagnetic wave equation inside this medium can be written into a Helmholtz form. Therefore, the fields can be expanded in terms of the vector cylindrical wave functions (VCWFs) [35].
According to the Maxwell’s equations and note that the system is invariant along z direction, namely, ∂z f = 0, where f is any field component, we can obtain that
whereNote that ∇×E = ik0B, ∇×H = −ik0D, where k0 = ω/c is the wavenumber in vacuum, and combine the first two subequations in Eq. (4), we obtain that
We can see that the electromagnetic fields are divergent-less only when κk = 0. Substituting Eq. (6) into Eq. (4), the last two subequations become
whereTherefore,
where we have used ∇×D = ∇ (∇·D) − ∇2D = −∇2D, andWe can formulate Eq. (9) into a Helmholtz equation as
whose basic solutions are the VCWFs.3. Mie coefficients of the bi-anisotropic cylinder
3.1. κk = 0
We first consider that the chirality tensor has no off-diagonal terms, i.e., κk = 0. This kind of bi-anisotropic metamaterial can be realized using a lattice of metallic helical structures [1, 4] and were used for realizing Weyl metamaterials [13, 14]. In the absence of κk, the Helmholtz equation Eq. (11) can be greatly simplified. Note that the in-plane fields are decoupled from the out-of-plane fields, we only need to consider the wave equation for the out-of-plane fields. Using the constitutive relations in Eq. (1), we formulated the eigenvalue problem as
The two wavenumbers inside the medium are determined by the eigenvalues of the 2 by 2 matrix in Eq. (12) which are
and the corresponding polarization settings are determined by the corresponding eigenvectorsConsider that an electromagnetic wave is scattered by a bi-anisotropic cylinder with radius r0. Based on the idea of Mie theory [35], we expanded the electric field inside the cylinder in series of the VCWFs as
where r = (r,ϕ) is the cylindrical coordinate with the origin located at the cylinder’s center, are expansion coefficients to be determined, and are the VCWFs with explicit expressions given by [35, 36]In Eq. (16), ρ = kr, and are the Bessel function and Hankel function of the first kind, and . The curls and divergences of VCWFs satisfy the following relations [35]
Note that the electromagnetic fields are divergent-less for κk = 0, according to Eq. (17), they are expanded only using and , see Eq. (15). The magnetic field inside the cylinder can be obtained according to
which givesAccording to the polarization settings defined in Eq. (14), and combing the expressions of the electric (Eq. (15)) and magnetic (Eq. (19)) fields, we obtained that
Assuming the background is air, the incident (Ei, Hi) and scattered (Es, Hs) fields can be expressed as
where pn, qn are beam shape coefficients and an, bn are scattering coefficients. They are related by the Mie coefficients as [35]Applying the boundary conditions that the tangential electromagnetic field components should be continuous at r = r0,
where x0 = k0r0, x± = k±r0 are the dimensionless size parameters, the Mie coefficients for the bi-anisotropic cylinder are obtained by solving Eqs. (20), (22) and (23). In the long-wavelength limit where x0 ≪ 1, the zeroth and first order Mie coefficients are3.2. κt = κz = 0
Next we consider that the chirality tensor possesses no diagonal terms, i.e., κt = κz = 0. This kind of bi-anisotropic metamaterial canbe realized using alattice of split-ring resonators [2, 10] or Ω-particles [37–39]. The matrix in Eq. (11) becomes diagonal. It is easy to obtain the two wavenumbers inside the medium as
and the corresponding polarization settings are Dx = Dy = Bz = 0 and Dz = Bx = By = 0, respectively. Therefore, the electric displacement and magnetic induced fields can be expanded as where cn, dn, en, fn are expansion coefficients, According to Eq. (7) and using the properties of VCWFs, see Eq. (17), we obtained that and the electromagnetic fields inside the medium asThe electromagnetic fields are non longer divergent-less because of the nonzero κk. Also imposing the boundary conditions,
where x1,2 = k1,2r0, the Mie coefficients are obtained by solving Eqs. (22), (27) and (29). In the long-wavelength limit, the zeroth and first order Mie coefficients areCombining Eqs. (24) and (30), we can conclude that the general expressions of Mie coefficients for bi-anisotropic cylinders with chirality tensor possessing both diagonal and off-diagonal terms are
4. Closed-form experssions for the effective constitutive parameters
Consider a metamaterial composed of identical bi-anisotropic inclusions arranged into a regular or random lattice embedded in air. The lattice constant fulfills k0a ≪ 1. Therefore, the local effective medium theory, such as the coherent potential approximation (CPA) method, can be applied due to the weak multiple scattering. Following the idea of CPA method [17], the zeroth and first order Mie coefficients of a single bi-anisotropic inclusion should equal to those of an effective cylinder with the effective constitutive parameters and occupying a unit cell volume, namely
where Ai, Bi, Ci, Di are Mie coefficients of the inclusion, and are the Mie coefficients of the effective cylinder calculated by substituting the effective constitutive parameters into Eq. (31). According to Eqs. (31) and (32), we can obtain the effective constitutive parameters as where with Ω being the volume of a unit cell. If the inclusion is cylindrical and homogeneous, substituting Eq. (31) into Eq. (33), we obtain that where is the filling ratio. For κk = 0, and consider the inclusion to be isotropic (its permittivity, permeability and chirality are scalars), Eqs. (33) and (34) are just reduced to Eqs. (20) and (22) in our previous work [40], which were derived using the multiple scattering technique. Eq. (34) is reduced to the well-known MG formula in the absence of chirality (κt = κz = κk = 0). Therefore, we can regard Eq. (34) as the extended MG formula for bi-anistropic metamateials. We note that although Eq. (33) is less efficient than Eq. (34) in calculating the effective constitutive parameters, it is more fundamental since it does not require the inclusion to be homogeneous. For the inclusion whose Mie coefficients satisfy A1 = A−1, B1 = B−1, C±1 = D∓1, Eq. (33) can be applied, even if the inclusion has complex underlying structures, e.g., the inclusion is made of helices [40]. The Mie coefficients of inclusion having complex structures need to be retrieved using the scattering fields in the far field region. However, the retrieval of the Mie coefficients only needs to be done once if the metamaterials with different filling ratios and lattice symmetries are composed of the same inclusions. This is a great advantage comparing with conventional parameter retrieval methods, which requires the repetition of the retrieval procedures when the filling ratio or lattice symmetry changes.Note that Eqs. (33) and (34) are valid only when the long-wavelength approximation holds, which requires that the incident wavelength is about 10 times larger than the lattice constant of the metamaterial. Therefore, the frequency usually should not be greater than c/(10a). When the frequency is too high, the effective constitutive parameters can only be obtained by the numerical retrieval, because for this case no analytic expressions can be derived due to the non-local effect and boundary sensitivity. Eqs. (33) and (34) are valid irrespective of the material parameters of the inclusions and the lattice structures. When the filling ratio is not too large, i.e., p < 0.5, Eqs. (33) and (34) are always accurate. However, if the filling ratio is too high, i.e., p > 0.5, the multiple scattering effect between cylinders will become important. Then the multiple scattering technique should be used and the high filling ratio corrections to the formulas should be taken into account [41].
5. Numerical testing
In this section, we will test the validity of our formula Eq. (34) using the full wave numerical simulations. Consider a plane wave obliquely incident on a metamaterial slab which is infinite along y direction and has a finite width d = 50a along x direction, where a is the lattice constant of the square lattice that the bi-anisotropic cylindrical inclusions are arranged into, see the inset in Fig. 2(a). For comparison, we also considered the same plane wave illuminates on an effective medium slab possessing the same width d and effective constitutive parameters calculated by Eq. (34). The validity of the effective parameters can be tested by checking the consistence between the spatially averaged lattice fields, such as the electric field or displacement field: , , where the spatial integrals are over the whole unit cell, and the fields in the effective medium Ee, De, since these two are equal in the long-wavelength limit [42,43].
In Fig. 2, we showed the spatially averaged lattice fields in the metamaterials having the square lattice structure and the corresponding fields in effective medium by circles and lines, respectively. All the fields are calculated using the finite-element method package COMSOL [44]. We can see that for different incident angles, the spatially averaged lattice fields (circles) both inside and outside the slab are in accordance with the fields in the effective medium (lines). Such consistence also exists for other field components and polarized incident waves. This indicates that our formula Eq. (34) can indeed determine the effective constitutive parameters correctly for this kind of bianisotropic metamaterials.
In Fig. 3, we also compared the spatially averaged lattice fields in the metamaterials having a triangular lattice structure and the corresponding fields in effective medium slab. Now 60 layers are used. The excellent agreement is also seen. In fact, because our formulas is obtained based on the local EMT, they are valid irrespective of the lattice symmetry.
The validity of Eq. (34) can also be tested according to the effective refractive indices, which are calculated either according to the slopes of the band dispersion of the metamaterial possessing a lattice structure, i.e., n/c = dk/dω in the limit of ω →0, k →0, or by n = k/k0 which are obtained by substituting the effective constitutive parameters into Eq. (13) or (25). The results by these two methods should be equal if the effective constitutive parameters are correct.
The refractive indices calculated by these two methods for κk = 0 and κt = κz = 0 are shown in Figs. 4(a) and 4(b), respectively. We first used Eq. (34) to calculate the effective constitutive parameters and then substituted them into Eqs. (13) and (25) to obtain the refractive indices analytically. The band dispersion for the square lattice is obtained using COMSOL. It is seen that the numerical results based on the band dispersion analysis (circles) match the analytic results using the formulas very well, again confirming the validity of Eq. (34).
The cases of strong resonances were also investigated. In Fig. 5, we showed the refractive indices calculated by these two methods when the cylinder is around the surface plasmon resonance. Resonances with very high quality factors occurr at κt ≈ 0.9 and κk ≈ 0.7 for Figs. 5(a) and 5(b), respectively. We can see that even very close to the resonance peak where the refractive index goes to infinity, the formulas can still produce the effective refractive indices accurately, indicating that our formulas are also valid for strong resonance cases. Note that the unplotted ranges for n− (κt < 0.9) and n2(κk < 0.7) are corresponding to the photonic band gaps for n− and n2, respectively.
6. Conclusion
In summary, by deriving the expressions of Mie coefficients and using the CPA method, we derived the closed-form expressions for the effective constitutive parameters of 2D bi-anisotropic metamaterials in the long-wavelength limit. Our formulas work for reciprocal metamaterials whose chirality tensors can have both diagonal and off-diagonal terms and satisfy κxy = −κyx = κk. The formulas can be regarded as the extended MG formulas for bi-anisotropic metamaterials since they are directly reduced to the traditional MG formulas when the chirality is vanishing. In principle, our method for deriving the closed-form formulas can be extended to the case that metamaterials are composed of non-reciprocal inclusions, e.g., the permittivity tensor possesses pure imaginary off-diagonal components, although the Mie coefficients will become more complicated. However, for κxy = κyx, our method fails because the wave equation inside the medium is no longer the Helmholtz form and the fields cannot be expanded in series of VCWFs.
Funding
Natural Science Foundation of SZU (2019009); National Natural Science Foundation of China (NSFC) (11734012 and 11574218).
References
1. J. B. Pendry, “A chiral route to negative refraction,” Science 306(5700), 1353–1355 (2004). [CrossRef] [PubMed]
2. 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]
3. S. Zhang, Y. S. Park, J. Li, X. Lu, W. Zhang, and X. Zhang, “Negative refractive index in chiral metamaterials,” Phys. Rev. Lett. 102(2), 023901 (2009). [CrossRef] [PubMed]
4. C. Wu, H. Li, Z. Wei, X. Yu, and C. T. Chan, “Theory and experimental realization of negative refraction in a metallic helix array,” Phys. Rev. Lett. 105(24), 247401 (2010). [CrossRef]
5. E. Plum, V. A. Fedotov, A. S. Schwanecke, and N. I. Zheludev, “Giant optical gyrotropy due to electromagnetic coupling,” Appl. Phys. Lett. 90(22), 223113 (2007). [CrossRef]
6. H. Liu, D. A. Genov, D. M. Wu, Y. M. Liu, Z. W. Liu, C. Sun, S. N. Zhu, and X. Zhang, “Magnetic plasmon hybridization and optical activity at optical frequencies in metallic nanostructures,” Phys. Rev. B 76(7), 073101 (2007). [CrossRef]
7. T. Q. Li, H. Liu, T. Li, S. M. Wang, F. M. Wang, R. X. Wu, P. Chen, S. N. Zhu, and X. Zhang, “Magnetic resonance hybridization and optical activity of microwaves in a chiral metamaterial,” Appl. Phys. Lett. 92(13), 131111(2008). [CrossRef]
8. S. L. Prosvirnin and N. I. Zheludev, “Polarization effects in the diffraction of light by a planar chiral structure,” Phys. Rev. E 71(3), 037603 (2005). [CrossRef]
9. V. A. Fedotov, P. L. Mladyonov, S. L. prosvirnin, A. V. Rogacheva, Y. Chan, and N. I. Zheludev, “Asymmetric propagation of electromagnetic waves through a planar chiral structure,” Phys. Rev. Lett. 97(16), 167401 (2006). [CrossRef] [PubMed]
10. E. Plum, X. X. Liu, V. A. Fedotov, Y. Chen, D. O. Tsai, and N. I. Zheludev, “Metamaterials: optical activity without chirality,” Phys. Rev. Lett. 102(11), 113902 (2009). [CrossRef] [PubMed]
11. A. B. Khanikaev, S. H. Mousavi, W.-K. Tse, M. Kargarian, A. H. MacDonald, and G. Shevets, “Photonic topological insulators,” Nat. Mater. 12(3), 233–239 (2013). [CrossRef]
12. W.-J. Chen, S.-J. Jiang, X.-D Chen, B. Zhu, L. Zhou, J.-W. Dong, and C. T. Chan, “Experimental realization of photonic topological insulator in a uniaxial metacrystal waveguide,” Nat. Commun. 5, 5782 (2014). [CrossRef] [PubMed]
13. W. Gao, M. Lawrence, B. Yang, F. Liu, F. Fang, B. Beri, J. Li, and S. Zhang, “Topological photonic phase in chiral hyperbolic metamaterials,” Phys. Rev. Lett. 114(3), 037402 (2015). [CrossRef] [PubMed]
14. B. Yang, Q. Guo, B. Tremain, R. Liu, L. E. Barr, Q. Yan, W. Gao, H. Liu, Y. Xiang, J. Chen, C. Fang, A. Hibbins, L. Lu, and S. Zhang, “Ideal Weyl points and helicoid surface states in artificial photonic crystal structures,” Science 359(6379), 1013–1016 (2018). [CrossRef] [PubMed]
15. T. C. Choy, Effective Medium Theory: Principles and Applications (Oxford University, 1999).
16. D. R. Smith, D. C. Vier, T. Koschny, and C. M. Soukoulis, “Electromagnetic parameter retrieval from inhomogeneous metamaterials,” Phys. Rev. E 71(3), 036617 (2005). [CrossRef]
17. Y. Wu, J. Li, Z.-Q. Zhang, and C. T. Chan, “Effective medium theory for magnetodielectric composites: Beyond the long-wavelength limit,” Phys. Rev. B 74(8), 085111 (2006). [CrossRef]
18. B. A. Slovick, Z. G. Yu, and S. Krishnamurthy, “Generalized effective-medium theory for metamaterials,” Phys. Rev. B 89(15), 155118 (2014). [CrossRef]
19. J. M. MacLaren, S. Crampin, D. D. Vvedensky, and J. B. Pendry, “Layer Korringa-Kohn-Rostoker technique for surface and interface electronic properties,” Phys. Rev. B 40(18), 12164–12175 (1989). [CrossRef]
20. A. Modinos, N. Stefanou, and V. Yannopapas, “Applications of the layer-KKR method to photonic crystals,” Opt. Express 8(3), 197–202 (2001). [CrossRef] [PubMed]
21. J. M. MacLaren, X.-G. Zhang, W. H. Butler, and X. Wang, “Layer KKR approach to Bloch-wave transmission and reflection: Application to spin-dependent tunneling,” Phys. Rev. B 59(8), 5470–5478 (1999). [CrossRef]
22. G.W. Milton, The Theory of Composites (Cambridge University, 2002). [CrossRef]
23. A. Priou, A. Sihvola, and S. Tretyakow, Advances in Complex Electromagnetic Materials (Springer Science & Business Media, 2012).
24. J. C. M. Garnet, “Colours in metal glasses and in metallic films,” Philos. Trans. R. Soc. London A 203, 385–420 (1904). [CrossRef]
25. A. H. Sihvola and I. V. Lindell, “Chiral Maxwell-Garnett mixing formula,” Electron. Lett. 26(2), 118–119 (1990). [CrossRef]
26. A. Lakhtakia, V. K. Varadan, and V. V. Varadan, “On the Maxwell-Garnett model of chiral composites,” J. Mater. Res. 8(4), 917–922 (1993). [CrossRef]
27. C. Menzel, C. Rockstuhl, T. Paul, and F. Lederer, “Retrieving effective parameters for quasiplanar chiral metamaterials,” Appl. Phys. Lett. 93(23), 233106 (2008). [CrossRef]
28. B. Wang, J. Zhou, T. Koschny, M. Kafesaki, and C. M. Soukoulis, “Chiral metamaterials: simulations and experiments,” J. Opt. A: Pure Appl. Opt. 11(11), 114003 (2009). [CrossRef]
29. R. Zhao, T. Koschny, and C. M. Soukoulis, “Chiral metamaterials: retrieval of the effective parameters with and without substrate,” Opt. Express 18(14), 14553–14567 (2010). [CrossRef] [PubMed]
30. X. Chen, B.-I. Wu, J. A. Kong, and T. M. Grzegorczyk, “Retrieval of the effective constitutive parameters of bianisotropic metamaterials,” Phys. Rev. E 71(4), 046610 (2005). [CrossRef]
31. Z. Li, K. Aydin, and E. Ozbay, “Determination of the effective constitutive parameters of bianisotropic metamaterials from reflection and transmission coefficients,” Phys. Rev. E 79(2), 026610 (2009). [CrossRef]
32. X. Jing, R. Xia, W. Wang, Y. Tian, and Z. Hong, “Determination of the effective constitutive parameters of bianisotropic planar metamaterials in the terahertz region,” J. Opt. Soc. Am. A 33(5), 954–961 (2016). [CrossRef]
33. J. Zhao, X. Jing, W. Wang, Y. Tian, D. Zhu, and G. Shi, “Steady method to retrieve effective electromagnetic parameters of bianisotropic metamaterials at one incident direction in the terahertz,” Opt. Laser Technol. 95, 56–62 (2017). [CrossRef]
34. J.A. Kong, Electromagnetic Wave Theory (Cambridge University, 2008).
35. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley and Sons, 1983).
36. N. Wang, J. Chen, S. Liu, and Zhifang Lin, “Dynamical and phase-diagram study on stable optical pulling force in Bessel beams,” Phys. Rev. A 87(6), 063812 (2013). [CrossRef]
37. A. N. Serdyukov, I. V. Semchenko, S. A. Tretyakov, and A. Sihvola, Electromagnetics of Bi-Anisotropic Materials: Theory and Applications (Gordon and Breach Science, 2001).
38. M. M. I. Saadoun and N. Engheta, “A reciprocal phase shifter using novel pseudochiral or ω medium,” Microw. Opt. Technol. Lett. 5(4), 184–188 (1992). [CrossRef]
39. S. A. Tretyakov, C. R. Simovski, and M. Hudlicka, “Bianisotropic route to the realization and matching of backward-wave metamaterial slabs,” Phys. Rev. B 75(15), 153104 (2007). [CrossRef]
40. N. Wang, S. Wang, Z.-Q. Zhang, and C. T. Chan, “Closed-form expressions for effective constitutive parameters: Electrostrictive and magnetostrictive tensors for bianisotropic metamaterials and their use in optical force density calculations,” Phys. Rev. B 98(4), 045426 (2018). [CrossRef]
41. W. Sun, S. B. Wang, Jack Ng, Lei Zhou, and C. T. Chan, “Analytic derivation of electrostrictive tensors and their application to optical force density calculations,” Phys. Rev. B 91(23), 235439 (2015). [CrossRef]
42. N. Wang, S. Wang, and J. Ng, “Electromagnetic stress tensor for an amorphous metamaterial medium,” Phys. Rev. A 97(3), 033839 (2018). [CrossRef]
43. M. G. Silveirinha, “Metamaterial homogenization approach with application to the characterization of microstructured composites with negative parameters,” Phys. Rev. B 75(11), 115104 (2007). [CrossRef]
44. www.comsol.com.