1. Introduction
Conventional optical single-material glass lenses are unable to focus all the spectral components of light into the same convergence point, even in ideal circumstances where the effects of diffraction are negligible. The reason for this limitation is the frequency dispersion of the glass refractive index, which causes wavelengths associated with different colors to be refracted differently [
1M. Born and M. Wolf, Principles of Optics, 7th (expanded) ed. (Cambridge University Press, 1999).
]. Hence the image produced by a glass lens may be distorted, and in such a case the optical system is said to suffer from chromatic aberrations. The material dispersion is manifested in the form of beautiful rainbows when white light is separated into its spectral components by a glass prism [
2. Newton, “A Letter of Mr. Isaac Newton, Professor of the Mathematicks in the University of Cambridge; Containing His New Theory about Light and Colors: Sent by the Author to the Publisher from Cambridge, Febr. 6. 1671/72; In Order to be Communicated to the R. Society,” Philos. Trans. R. Soc. Lond. 6(69-80), 3075–3087 (1671). [CrossRef]
], and is rooted in fundamental physical restrictions, stemming from the causality and passivity of the dielectric response. Causality and passivity determine that the index of refraction of any conventional low loss dielectric material is a strictly increasing function of frequency [
3L. D. Landau, E. Lifshitz, and L. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Elsevier Butterworth-Heinemann, 2004).
].
Since the optical path length of a given ray in an optical system,
, is written in terms of the indices of refraction (
) of the involved materials, the correction of chromatic aberrations is a non-trivial problem. Indeed, the material dispersions are combined additively, and since
, it follows that
is a strictly increasing function of frequency. Nevertheless, several strategies to minimize the effect of chromatic aberrations are well documented in the literature and are typically based on the combination of materials with different positive dispersion (e.g [
4T. T. Smith, “The color correction of an achromatic doublet,” J. Opt. Soc. Am. 10(1), 39–62 (1925). [CrossRef]
–
7M. Herzberger and N. R. McClure, “The design of superachromatic lenses,” Appl. Opt. 2(6), 553–560 (1963). [CrossRef]
].). This is possible because the system can be designed so that the trajectories of the rays inside the lens change with frequency in such a manner that the profile of
at the exit surface is invariant (apart from the sum of an irrelevant constant). Lenses with reduced aberrations are known as achromatic doublets.
The definition of
suggests that the effects of material dispersion in glass (
) may be easier to suppress if a material with anomalous dispersion (
) is available. However, the aforementioned restrictions stemming from Kramers-Kronig relations for causal and passive dielectric media [
3L. D. Landau, E. Lifshitz, and L. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Elsevier Butterworth-Heinemann, 2004).
], indicate that a regime of anomalous dispersion implies very significant loss, and hence this solution seems to be impractical. Is it possible to overcome this limitation?
Structured functional materials with extended electromagnetic responses have been on the spotlight in recent years [
8J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85(18), 3966–3969 (2000). [CrossRef] [PubMed]
–
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
]. Typically, such metamaterials consist of periodic arrangements of metallic or dielectric particles embedded in a host material. The key feature of metamaterials is that their electromagnetic response is mainly determined by the geometry of the inclusions, rather from the chemical properties of its constituents, allowing for a greater control of their effective response. In particular, it has been shown in a recent work [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
] that surprisingly the restrictions on
coming from Kramers-Kronig relations are less severe for the class of metamaterials with spatial dispersion. Different from conventional dielectrics, spatially dispersive materials are characterized by the fact that the polarization vector responds in a nonlocal manner to the macroscopic electric field [
3L. D. Landau, E. Lifshitz, and L. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Elsevier Butterworth-Heinemann, 2004).
]. It was shown in Ref [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
], that a metamaterial formed by an array of nanowires (a crossed wire mesh) enables a broadband low-loss anomalous dispersion regime, such that
. Here, we theoretically demonstrate that such a metamaterial may permit reducing significantly the chromatic aberration of a conventional thin glass lens.
To begin with, we consider a thin compound lens formed by two materials. The refracting surfaces at the air interfaces have radius of curvature
and
. We use the convention that
and
are positive for convex surfaces (seen from the air region), and assume that the
i-th (
i = 1,2) refracting material is associated with a material with index
[see
Fig. 1(a)
].
Fig. 1 (a) Illustration of the chromatic aberration of a conventional thin biconvex glass lens with parameters . (b) Biconvex optical metamaterial lens that corrects the chromatic aberration for all the colors of light; the compound lens is formed by a thin plano-convex glass lens coated with a thin plano-convex double wire medium. Panels (c) and (d) show cuts of the “double wire medium” slab along the planes xoz and xoy, respectively. The slab has thickness .
Since the radii of conventional optical lenses have physical sizes that correspond to tens or even hundreds of wavelengths of the visible spectrum, geometrical optics may be used to describe how light propagates in such optical systems. For a thin compound lens standing in air, the focal length satisfies the well-known Lensmaker’s equation [
1M. Born and M. Wolf, Principles of Optics, 7th (expanded) ed. (Cambridge University Press, 1999).
]:
For a thin lens, the distances
and
, from the object and focal plane to the lens, respectively, satisfy
Obviously, a lens formed by a single material (let’s say glass), has focal length such that
where
is the refractive index of the glass. Since the permittivity
of glass is an increasing function of frequency in the optical domain [
3L. D. Landau, E. Lifshitz, and L. Pitaevskii, Electrodynamics of Continuous Media, 2nd ed. (Elsevier Butterworth-Heinemann, 2004).
], it is manifest from
Eq. (1) that the focal length
of the lens is a decreasing function of frequency. Therefore, the focusing provided by the biconvex lens of
Fig. 1(a) is imperfect, as the colors of light associated with longer wavelengths (“red” light) are less refracted by the optical system, converging to a longitudinal point farther than the components of light associated with shorter wavelengths (“violet” light). The change
in the focal length caused by a change
in the refractive index of the glass is the chromatic aberration.
For a compound lens we may suppose that to a first approximation both
and
vary linearly with frequency so that
and
where
and that
is some reference frequency at which the lens is designed (and for which
and
). Then, using
Eq. (1), we easily find that the optical power of the bi-layer lens can be made independent of frequency provided [
4T. T. Smith, “The color correction of an achromatic doublet,” J. Opt. Soc. Am. 10(1), 39–62 (1925). [CrossRef]
]:
From the above formula it is manifest that the correction of the chromatic aberrations can be achieved by considering two materials with positive dispersion (when
and
are both positive) and such that the radii of curvature of two refracting surfaces have opposite signs. This is the conventional solution that is the basis of achromatic doublets and that has been used for decades [
4T. T. Smith, “The color correction of an achromatic doublet,” J. Opt. Soc. Am. 10(1), 39–62 (1925). [CrossRef]
].
Let us now consider that the thin lens is instead formed by two convex refracting surfaces, i.e.
and
are positive. In such a scenario, in order to satisfy the condition
(2) a material with anomalous dispersion is required, because
and
must have opposite signs. As far as we know, this solution has not been seriously considered before because regimes of anomalous dispersion in conventional materials imply very significant loss. Here, we show that such a solution is within the realm of reality if the material with anomalous dispersion is taken as a metamaterial with engineered dispersion. For this case we put
and
being
the effective index of refraction of the metamaterial. Hereafter, a lens that satisfies
Eq. (2) with
will be referred to as
compensated biconvex lens.
2. Metamaterial with broadband anomalous dispersion
An interesting possibility to realize the required material with anomalous dispersion is the nanowire metamaterial investigated in Ref. [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
]. This material is formed at the nanoscopic level by two nonconnected arrays of metallic nanowires, such that each array of parallel wires is arranged in a square lattice with lattice constant
and tilted by
with respect to the
x axis. The two arrays of wires are mutually orthogonal and lie in planes parallel to the
xoz plane; the distance between adjacent perpendicular wires is
[
Figs. 1(c) and
1(d)]. For simplicity of modeling, and without loss of generality, it will be assumed that the wires stand in air.
The nanowire metamaterial is characterized by a strong spatial dispersion, i.e, the permittivity seen by a wave with propagation factor
depends on the wave vector
[
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
]. For propagation in the
xoy plane with
and assuming that the electric field is polarized along the
z direction, the effective permittivity of this material is [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
,
13M. G. Silveirinha and C. A. Fernandes, “Nonresonant structured material with extreme effective parameters,” Phys. Rev. B 78(3), 033108 (2008). [CrossRef]
]:
where
is the plasma wave number,
and
are the radius and the complex permittivity of the metallic wires, respectively, and
is the volume fraction of each set of wires. The effective index of refraction,
of the metamaterial is found by solving the dispersion equation
with respect to
. This yields:
In the particular case of perfectly conducting wires (
), the index of refraction reduces to
which is clearly a decreasing function of frequency [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
], since
is a parameter that depends merely on the geometry of the metamaterial. This property still holds for realistic metals at optical frequencies. This is illustrated in
Fig. 2(a)
, where we plot
as a function of frequency for Al nanowires modeled by a Drude dispersion model with parameters consistent with experimental data reported in the literature [
14M. A. Ordal, R. J. Bell, R. W. Alexander, Jr, L. L. Long, and M. R. Querry, “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W,” Appl. Opt. 24(24), 4493–4499 (1985). [CrossRef] [PubMed]
]. It is seen that the metamaterial is characterized by low-loss broadband anomalous dispersion in the entire visible domain.
Fig. 2 (a) refractive indices of dense flint glass SF10 (black curve) and of the nanowire metamaterial (blue curves). Panels (b), (c) and (d): normalized
in the vicinity of a metamaterial prism at: (b)
.(c)
. (d)
. (e) reversed rainbow obtained by blending the different light wavelengths [panels (b), (c) and (d)]. (f) transmission angle
as a function of the wavelength
for a Gaussian beam that illuminates the prism along the normal direction. The blue dashed curve was obtained using the theoretical formula
and the solid black curve was calculated using a full wave FDFD-SD [
18J. T. Costa and M. G. Silveirinha, “Macroscopic electromagnetic response of arbitrarily shaped spatially dispersive bodies formed by metallic wires,” arXiv:1205.1760v1 [physics.comp-ph], (2012).
] simulation based on the effective medium model
It is curious to mention that notwithstanding that at the “microscopic level” all the constituents (host and inclusions) of the metamaterial satisfy the effective medium is characterized by due to the complex electromagnetic interactions between its different elements. In particular, the optical path length of a wave that transverses a metamaterial slab decreases with frequency, whereas a naïve application of the formula (with standing either for the inclusions or for the host) would suggest the opposite. This apparent contradiction is explained by the fact that the formula is not valid when the distance between the different materials is small on the scale of the wavelength, because in such conditions the wave envelope is described by an effective index of refraction, which is not a simple average of those of the constituent materials. In some sense, in a nanowire metamaterial the sum of many positives () can yield a negative optical path length ().
To numerically solve the Maxwell’s equations in a scenario wherein electromagnetic waves interact with spatially dispersive bodies (using effective medium theory) such as the nanowire metamaterial, we use the finite differences frequency domain (FDFD) method [
15K. Yasumoto, Electromagnetic Theory and Applications for Photonic Crystals (Taylor and Francis, 2006).
]. The FDFD method that we use to characterize spatially dispersive materials (FDFD-SD) is based on ideas analogous to those described in Refs. [
16J. M. McMahon, S. K. Gray, and G. C. Schatz, “Nonlocal optical response of metal nanostructures with arbitrary shape,” Phys. Rev. Lett. 103(9), 097403 (2009). [CrossRef] [PubMed]
,
17Y. Zhao, P. A. Belov, and Y. Hao, “Modelling of wave propagation in wire media using spatially dispersive finite-difference time-domain method: numerical aspects,” IEEE Trans. Antenn. Propag. 55(6), 1506–1513 (2007). [CrossRef]
], and specific details are reported in Ref [
18J. T. Costa and M. G. Silveirinha, “Macroscopic electromagnetic response of arbitrarily shaped spatially dispersive bodies formed by metallic wires,” arXiv:1205.1760v1 [physics.comp-ph], (2012).
]. In order to confirm that the effective refractive index of the double wire medium is a decreasing function of frequency in the optical domain, we investigated the refraction of a cylindrical Gaussian beam by a metamaterial prism with
and aperture
. The metamaterial has lattice constant
and the nanowires have radius
. The beam waist of the incoming wave was taken equal to
.
Figures 2(b),
2(c) and
2(d) show density plots of
in the vicinity of the prism at
and
respectively, and support that shorter wavelengths are less refracted than longer wavelengths, confirming that the effective refractive index
of the metamaterial prism is a decreasing function of frequency in the optical domain. The spectral electromagnetic fields associated with
and
[
Figs. 2(b),
2(c) and
2(d)] can be blended and represented in a RGB color scale − taking into account the relative intensity of the fields for each wavelength—and this results in a reversed rainbow [
Fig. 2(e)]. In
Fig. 2(f) the angle of transmission,
is depicted as a function of frequency calculated with
(i) our FDFD-SD code (black solid curve) [
18J. T. Costa and M. G. Silveirinha, “Macroscopic electromagnetic response of arbitrarily shaped spatially dispersive bodies formed by metallic wires,” arXiv:1205.1760v1 [physics.comp-ph], (2012).
]
(ii) the theoretical formula
(blue dashed curve) [
12M. G. Silveirinha, “Anomalous refraction of light colors by a metamaterial prism,” Phys. Rev. Lett. 102(19), 193903 (2009). [CrossRef] [PubMed]
]. It is seen that the results concur very well.
3. Achromatic biconvex metamaterial lens
To illustrate how the considered metamaterial can, indeed, nearly eliminate the chromatic aberration, we consider the design of a compensated biconvex lens such that
and
,
. The nanowire material has the same parameters as in the previous example, except that the nanowires diameter is increased to
. We assume that the dielectric is a dense flint glass SF10. In the visible spectrum, the refractive index
of this glass is described to a first approximation by the so-called Cauchy’s equation
with
and
[
19F. A. Jenkins and H. E. White, Fundamentals of Optics, 4th ed. (McGraw-Hill, 1981).
]. At the central frequency of the visible spectrum (
), we can estimate
,
,
, and
. Hence, from
Eqs. (1) and
(2) it follows that the radius of curvature of the first and second interfaces is
and
, respectively. The central thickness of the two layers is taken
and
.
To have a reference against which we can compare the performance of the compensated biconvex lens, we consider as well an ordinary single material plano-convex lens (
) made exclusively of glass (
), and with the same optical power as the compensated lens (i.e.
is invariant). This requires that the radius of curvature
of the convex surface is taken equal to:
. The central thickness of the single-material lens is
.
Figure 3(a)
shows the electric field profile calculated along the mid-plane perpendicular to the single-material lens under the illumination of a cylindrical Gaussian beam. The dashed and solid curves in
Fig. 3(a) were obtained using a commercial full-wave electromagnetic simulator [
20] and our FDFD code, respectively. The FDFD code predicts that the chromatic aberration
(defined with respect to the central frequency
) at
and
, is
and
(
Fig. 3(a)), respectively.

Fig. 3 (a) Profile of the normalized squared electric field near the focal region of the single-material glass lens. Solid curves: obtained using Microwave Studio [
20]. Dashed curves: obtained with the FDFD simulator. (b) normalized
(obtained using the FDFD-SD full wave simulator [
18J. T. Costa and M. G. Silveirinha, “Macroscopic electromagnetic response of arbitrarily shaped spatially dispersive bodies formed by metallic wires,” arXiv:1205.1760v1 [physics.comp-ph], (2012).
]) in the vicinity of the compensated biconvex metamaterial lens. (c) analogous to (a), but for the compensated biconvex metamaterial lens. (d) focal plane curve as a function of the wavelength. The yellow stars (FDFD-SD method based on the effective medium model [
18J. T. Costa and M. G. Silveirinha, “Macroscopic electromagnetic response of arbitrarily shaped spatially dispersive bodies formed by metallic wires,” arXiv:1205.1760v1 [physics.comp-ph], (2012).
]) and the yellow circles (Microwave Studio [
20]) represent the position of the foci of the compensated metamaterial lens. The red triangles and the diamond symbols (FDFD) represent the focal curve of a conventional achromatic doublet and of the single-material lens, respectively.
The achromatic biconvex metamaterial lens performs far better than the ordinary lens.
Figure 3(b) shows a density plot of
obtained with the FDFD-SD code at
in the vicinity of the biconvex lens. The focal spot created by the lens under the illumination of a cylindrical Gaussian beam is evident.
Figure 3(c) shows the electric field profile calculated along the mid-plane perpendicular to the biconvex lens. For this case, our FDFD-SD code predicts that the chromatic aberration at the wavelength
is drastically reduced to
, whereas the chromatic aberration associated with the wavelength
is reduced to
. Thus, the chromatic aberrations associated with the wavelengths at the edges of the visible spectrum are practically eliminated. This property is actually valid in the entire visible spectrum, as supported by
Fig. 3(d) that shows the focal curve, .i.e., the chromatic aberration
as a function of the wavelength. The yellow stars were calculated using the FDFD-SD simulations whereas the yellow circles were obtained using the full wave simulator taking into account all the minute details of the metamaterial [
20]. It can be seen that the compensated lens permits focusing all the colors into the same convergence point. The diamond symbols in
Fig. 3(d) represent the focal curve of the ordinary single-material glass lens of
Fig. 3(a), showing a significant chromatic aberration.
It should be noted that
Eq. (1) is based on the thin lens approximation, and thus even within the framework of geometrical optics it is not exact. To take this into account and have a more robust correction of the chromatic aberrations, we have slightly adjusted the parameters of the metamaterial in the FDFD-SD simulations of
Fig. 3 so that
and
. Moreover, in the CST simulations of
Fig. 3, the lattice constant of the nanowire material and the diameter of the wires were adjusted to
and
, in order to have a physical response more consistent with the homogenization model.
We have also investigated the performance of a standard doublet formed by a biconvex crown glass N-BK7 lens [
21OSA, Handbook of Optics, 2nd ed. (McGraw-Hill Professional, 1994), vol. II.
] and a concave-plane dense flint glass SF10 lens. In this scenario, the doublet is formed by two materials with positive dispersion (
and
are both positive). Note that this achromat is not of the same type as that represented in
Fig. 1, wherein the middle interface is planar. Indeed, we verified with both ray tracing and full wave simulations that the generic configuration of
Fig. 1 is less effective in the correction of the aberrations when both materials have positive dispersion than the solution (based on biconvex and concave-plane layers) reported here. In order that the optical power of the doublet is the same as that of the compensated lens, we chose the radii of curvature of the first and second convex surfaces (N-BK7)
and
, and the radius of curvature of the concave-plane surfaces (SF10)
and
.The central thicknesses of both glasses is
. The focal curve of this doublet is also represented in
Fig. 3(d) (red triangles), showing a reduction of the chromatic aberration comparable to that of the compensated lens.