The Bruggeman formalism provides a well-established technique for estimating the effective constitutive parameters of homogenized composite materials (HCMs).12.–3 The scope of its applicability is not restricted to dilute composite materials and it is easy to implement numerically, both of which contribute to its enduring popularity.
The Bruggeman formalism was originally devised for isotropic dielectric HCMs, comprising two (or more) isotropic dielectric component materials distributed randomly as electrically small spherical particles.4 Generalizations of the Bruggeman formalism which accommodate anisotropic and bianisotropic HCMs have been developed.5 A rigorous basis for the Bruggeman formalism—for isotropic dielectric,6 anisotropic dielectric,7,8 and bianisotropic9,10 HCMs—is provided by the strong-property-fluctuation theory, whose lowest-order formulation is the Bruggeman formalism. (In Ref. 7, the formula of Polder and van Santen, for the effective relative permittivity of an isotropic dielectric HCM, yields the same results as the formula of Bruggeman.4)
Our focus in this letter is on HCMs arising from two isotropic dielectric component materials, labeled and . Their relative permittivities are and , while their volume fractions are and . Both component materials are assumed to be randomly distributed as electrically small ellipsoidal particles. For simplicity, all component particles have the same shape and orientation. The surface of each ellipsoid, relative to its centroid, may be represented by the vector
The ellipsoidal shape of the component particles results in the corresponding HCM being an orthorhombic biaxial dielectric material. That is, the Bruggeman estimate of the HCM relative permittivity dyadic has the form11 5 4) means that numerical methods are generally needed to extract the relative permittivity parameters from them.
An alternative formalism for the homogenization of the same composite material as in the foregoing paragraph is also referred to as the Bruggeman formalism.1213.14.–15 Let us write the estimate of the HCM’s relative permittivity dyadic provided by this alternative formalism as16 9) is a quadratic equation in whose solution may be explicitly expressed as 12) may be determined by appealing to the anisotropic dielectric generalization of the Hashin–Shtrikman bounds,17 for example.
Let us illustrate the differences between the estimates and by means of some representative numerical results. The two estimates are identical for the limiting case represented by (i.e., for isotropic dielectric HCMs), but differences emerge as the asphericity of the component particles intensifies. Suppose that the shape parameters describing the component ellipsoids have the form , , and . Thus, the asphericity of the ellipsoids is governed by the scalar parameter . We begin with the nondissipative scenario wherein and . Also, we fix . Plots of the relative permittivity parameters versus the asphericity parameter are presented in Fig. 1. The difference between and grows steadily as increases, reaches a maximum for , and then slowly shrinks as increases beyond . The difference between and follows a similar pattern. However, in the case of and , the difference increases uniformly as increases. The differences between and are generally greater for than for . In the former case the maximum difference is approximately 15%, whereas in the latter case it is approximately 5%.
We turn now to the effect of volume fraction. The calculations of Fig. 1 are repeated for Fig. 2 except that here the relative permittivity parameters are plotted versus the volume fraction , while the asphericity parameter is fixed at . The differences between the estimates of the two formalisms are clearly greatest at mid-range values of , and they are generally greater for than for .
Lastly, the effects of dissipation are considered. We repeated the calculations of Fig. 1 but with and . Here governs the degree of dissipation exhibited by component material . The real and imaginary parts of the relative permittivity parameters are plotted versus the dissipation parameter in Fig. 3. The differences between the real parts of the estimates and are largest when component material is nondissipative and they decrease uniformly as increases. In contrast, the differences between the imaginary parts of the estimates and increase as increases. These differences in the imaginary parts generally reach a maximum for mid-range values of and thereafter decrease as increases. For both the real and imaginary parts of the estimates and , generally larger differences arise for than for .
Thus, there are significant differences between the estimates and when ellipsoidal component particles are considered. These differences may be exacerbated: (i) if the component particles become more aspherical, (ii) at mid-range values of the volume fractions of the component materials, and (iii) if the HCM is dissipative. The differences between the two estimates may be further exacerbated if one of the component materials has a positive-valued relative permittivity which is less than unity (or a relative permittivity whose real part is positive-valued and less than unity). Relative permittivities in this range are associated with novel materials possessing engineered nanostructures; these artificial materials have been the subject of intense research lately.1819.–20 (The parameter regime wherein one of the component materials has a positive-valued relative permittivity while the other has a negative-valued relative permittivity (or likewise for the real parts of the relative permittivities in the case of dissipative HCMs) is avoided here because the Bruggeman formalism can deliver estimates in this regime which are not physically plausible.21)
The differences between the two formalisms stem from the differences between the depolarization dyadics and . The Bruggeman formalism conceptually employs an average-polarizability-density approach:22 Suppose the composite material has been homogenized into an HCM. Into this HCM, let additional aligned ellipsoidal particles of the two component materials be dispersed in such a way as to maintain the overall volume fractions of and . This dispersal cannot change the effective properties of the HCM—because the volume-fraction-average of the polarization density dyadics of two particles, one of each component material, embedded in the HCM with relative permittivity is exactly zero. In computing the polarizability density dyadic of each particle, it must therefore be assumed that the particle is surrounded by the HCM. This fact legitimizes the use of , which contains the anisotropic HCM’s effective constitutive properties via the scalar of Eq. (7). Indeed, Eq. (4) is directly derivable from the frequency-domain Maxwell postulates specialized for the actual HCM.67.8.9.–10
On the other hand, use of to compute the polarizability density dyadic of a particle implies that it is surrounded by an isotropic HCM, which is clearly incorrect. The alternative formalism that delivers is a heuristic extrapolation of the Bruggeman formalism for isotropic dielectric HCMs,12,13 and it lacks the rigorous basis that underpins the estimate of . The dispersal of aligned ellipsoidal particles of the two component materials without disturbing the overall volume fractions into the HCM with relative permittivity would change the effective properties of the HCM—because the volume-fraction-average of the polarization density dyadics of two particles, one of each component material, embedded in the HCM with relative permittivity would not be zero.
We have thus delineated the differences between the two formalisms and identified one of them as the correct Bruggeman formalism. We hope that this exposition will prevent confusion between the two formalisms from perpetuating.
L. Ward, The Optical Constants of Bulk Materials and Films, 2nd ed., Institute of Physics, Bristol, UK (2000).Google Scholar
A. Lakhtakia, Ed., Selected Papers on Linear Optical Composite Materials, SPIE Optical Engineering Press, Bellingham, WA (1996).Google Scholar
T. G. Mackay, “Effective constitutive parameters of linear nanocomposites in the long-wavelength regime,” J. Nanophoton. 5(1), 051001 (2011), http://dx.doi.org/10.1117/1.3626857.JNOACQ1934-2608Google Scholar
D. A. G. Bruggeman, “Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen, I. Dielektrizitätskonstanten und Leitfähigkeiten der Mischkörper aus isotropen Substanzen,” Ann. Phys. Lpz. 24, 636–679 (1935). (Reproduced in Ref. 2).Google Scholar
W. S. WeiglhoferA. LakhtakiaB. Michel, “Maxwell Garnett and Bruggeman formalisms for a particulate composite with bianisotropic host medium,” Microw. Opt. Technol. Lett. 15(4), 263–266 (1997), http://dx.doi.org/10.1002/(ISSN)1098-2760;MOTLEO0895-2477Google Scholar
L. TsangJ. A. Kong, “Scattering of electromagnetic waves from random media with strong permittivity fluctuations,” Radio. Sci. 16(3), 303–320 (1981), http://dx.doi.org/10.1029/RS016i003p00303. (Reproduced in Ref. 2).RASCAD0048-6604Google Scholar
Z. D. Genchev, “Anisotropic and gyrotropic version of Polder and van Santen’s mixing formula,” Waves Random Media 2, 99–110 (1992), http://dx.doi.org/10.1088/0959-7174/2/2/001.WRMEEV0959-7174Google Scholar
N. P. Zhuck, “Strong–fluctuation theory for a mean electromagnetic field in a statistically homogeneous random medium with arbitrary anisotropy of electrical and statistical properties,” Phys. Rev. B 50(21), 15636–15645 (1994), http://dx.doi.org/10.1103/PhysRevB.50.15636.PRBMDO0163-1829Google Scholar
T. G. MackayA. LakhtakiaW. S. Weiglhofer, “Strong–property–fluctuation theory for homogenization of bianisotropic composites: formulation,” Phys. Rev. E 62(5), 6052–6064 (2000), http://dx.doi.org/10.1103/PhysRevE.62.6052; PLEEE81063-651XGoogle Scholar
T. G. MackayA. LakhtakiaW. S. Weiglhofer, “Third-order implementation and convergence of the strong-property-fluctuation theory in electromagnetic homogenisation,” Phys. Rev. E 64(6 Pt 2), 066616 (2001), http://dx.doi.org/10.1103/PhysRevE.64.066616.PLEEE81063-651XGoogle Scholar
G. B. Smith, “Effective medium theory and angular dispersion of optical constants in films with oblique columnar structure,” Opt. Commun. 71(5), 279–284 (1989), http://dx.doi.org/10.1016/0030-4018(89)90008-4. (Reproduced in Ref. 2).OPCOB80030-4018Google Scholar
C. G. GranqvistD. Le BellacG. A. Niklasson, “Angular selective window coatings: effective medium theory and experimental data on sputter-deposited films,” Renew. Energ. 8(1–4), 530–539 (1996), http://dx.doi.org/10.1016/0960-1481(96)88913-0.RNENE30960-1481Google Scholar
D. SchmidtE. SchubertM. Schubert, “Optical properties of cobalt slanted columnar thin films passivated by atomic layer deposition,” Appl. Phys. Lett. 100, 011912 (2012), http://dx.doi.org/10.1063/1.3675549.APPLAB0003-6951Google Scholar
D. PolderJ. H. van Santen, “The effective permeability of mixtures of solids,” Physica 12(5), 257–271 (1946), http://dx.doi.org/10.1016/S0031-8914(46)80066-1. (Reproduced in Ref. 2).PHYSAG0031-8914Google Scholar
Z. HashinS. Shtrikman, “A variational approach to the theory of the effective magnetic permeability of multiphase materials,” J. Appl. Phys. 33, 3125–3131 (1962), http://dx.doi.org/10.1063/1.1728579. (Reproduced in Ref. 2).JAPIAU0021-8979Google Scholar
A. Alùet al., “Epsilon-near-zero metamaterials and electromagnetic sources: tailoring the radiation phase pattern,” Phys. Rev. B 75, 155410 (2007), http://dx.doi.org/10.1103/PhysRevB.75.155410.PRBMDO0163-1829Google Scholar
G. Lovatet al., “Combinations of low/high permittivity and/or permeability substrates for highly directive planar metamaterial antennas,” IET Microw. Antennas Propagat. 1, 177–183 (2007), http://dx.doi.org/10.1049/iet-map:20050353.1751-8725Google Scholar
M. N. Navarro-Cíaet al., “Enhanced lens by and near-zero metamaterial boosted by extraordinary optical transmission,” Phys. Rev. B 83, 115112 (2011), http://dx.doi.org/10.1103/PhysRevB.83.115112.PRBMDO0163-1829Google Scholar