# High Resolution through Graded-Index Microoptics

Authors: Victor V. Kotlyar; Alexey A. Kovalev; Anton G. Nalimov; Sergey S. Stafeev
Journal: Advances in Optical Technologies
 (2012)

---

## Abstract

By solving Helmholtz equations, relationships to describe propagating modes in an arbitrary graded-index planar waveguide are derived. We show that in the quadratic- and secant-index waveguides a minimal mode width is 0.4λ/n, where λ is the wavelength in free space and n is the refractive index on the fiber axis. By modeling in FullWAVE, we show that the high-resolution imaging can be achieved with half-pitch graded-index Mikaelian microlenses (ML) and Maxwell’s “fisheye” lenses. It is shown that using a 2D ML, the point source can be imaged near the lens surface as a light spot with the full width at half maximum (FWHM) of 0.12λ. This value is close to the diffraction limit for silicon (n=3.47) in 2D media FWHM=0.44λ/n=0.127λ. We also show that half-pitch ML is able to resolve at half-maximum two close point sources separated by a 0.3λ distance.


---

## Body

## 1. Introduction

Recent advances in microoptics and nanophotonics have made possible the focusing of coherent laser light into a subwavelength spot or the superresolution imaging of a point source of light. The subwavelength focusing beyond the diffraction limit of0.5λ/n, where λ  is the wavelength in free space andnis the material refractive index at the focus, can be performed using a superlens [1]. In 2D case, instead of conventional diffraction limit 0.5λ/n one must use 0.44λ/n. This value can be obtained after replacing the Airy disk 2J1(x)/x by sinc-function sin(x)/x.The superlens is a 2D planar plate made up of the metamaterial that comprises alternating metallic and dielectric layers. The electric permittivities of the layers are selected so that an effective refractive index of the composite material be equal ton=-1. Experiments on the superresolution through superlenses were reported in [2, 3]. In the experiments, a superresolution of 0.4λ was achieved [2]. A similar experiment conducted in [4] with a subwavelength silver layer operating as a superlens has shown that two lines separated by a 145 nm distance can be resolved when illuminated by UV light of wavelength 365 nm, thus producing a superresolution of 0.4λ.A far-field hyperlens reported in [5] was able to resolve two lines of width 35 nm spaced 150 nm apart for a 365 nm wavelength, again achieving a superresolution of 0.4λ. Note, however, that a hyperlens modelled in the form of a grating [6] was shown to achieve a superresolution of 0.05λ at the imaging plane found 1.5λ apart from the surface. Apparently, the absorption and scattering of light by metamaterial that occurs in real experiments was disregarded in modelling. This argument was indirectly verified by results reported in [7], in which the laser light was focused with a zone plate made up of a gold film of thickness 100 nm. The focal spot size at half maximum was estimated to be 0.35λ,whereas the experimentally measured size of the focal spot at half maximum was found to equal λ.Multilayer and anisotropic nanostructures allowing one to achieve a subwavelength resolution were analyzed in [8–10]. For example, parameters of a 1D eight-layer Ag/SiO2  structure of thickness 400 nm to focus light from a 0.4λ source into the same-size focal spot were studied in [8]. An anisotropic 2D nanostructure characterized by the dielectric permittivity tensor components  εx=0.01-i0.01  and εz=-100, and a 400 nm thickness on the z-axis (for λ=700 nm) was proposed in [9]. The modelling has shown that such a structure is able to resolve two lines of width 3 nm placed 23 nm apart, providing a superresolution of 0.03λ. A nanostructure composed of two different anisotropic layers to resolve two narrow slits placed 50 nm apart when illuminated by a 1550 nm wavelength was studied in [10] by the same authors.Candidates for achieving the superresolution can be found among the photonic crystals. Modelling conducted in [11] has shown that a 2D photonic-crystal slab with permittivity ε=12 composed of a triangular array of circular holes of radius r=0.4a (a is the hole array period) can be used as an imaging lens for wavelength λ=a/0.3. In this case, a point source is imaged as a focal spot of size 0.3λ, whereas two point sources placed 0.5λ apart are resolved by the 20% criterion.In recent experiments with a 2D photonic-crystal slab used as a superlens, a point source of size0.4λ was imaged [12]. The latest publications propose an improved variant of superlenses using a nanoshell [13] or a graded-index boundary of the negative-refraction material [14]. In [15] an anisotropic layer was experimentally demonstrated to enhance and transform the evanescent surface waves into propagating light modes.From the previous survey it follows that although theoretically superlenses allow to achieve arbitrary high resolution, experimentally the values0.3λ-0.4λ have been obtained [2, 4, 5, 7, 11, 12].Such values of resolution can be obtained with help of gradient-index optics as well. Gradient optics works as near-field optics: gradient lens is placed near the object and the image emerges in vicinity of the exit surface of the lens. Therefore, lens material affects the resolution. The limiting resolution, which can be obtained with gradient lens, decreasesn times in comparison with conventional refraction lens and is 0.51λ/n in 3D and 0.44λ/n in 2D cases. If one will use refraction lens with immersion, then index of immersion liquids does not exceed 1.5, although gradient lenses can be made of silicon with index 3.47.In this work, we numerically demonstrate that widely known in optics graded-index Mikaelian lenses and Maxwell “fisheye” lenses may also be considered as candidates for subwavelength focusing. For a 2D Mikaelian lens, we show that a point source is imaged near the lens surface as a focal spot of sizeFWHM=0.12λ (full width at half maximum). This value is close to the diffraction limit for silicon (n=3.47) in 2D media FWHM=0.44λ/n=0.127λ. This value is smaller than values earlier reported in [7, 8, 11, 12]. We also show that half-pitch Mikaelian lens can resolve at half-maximum two point sources placed 0.3λ apart, which is smaller than reported in [2, 4, 5, 11].Analytical relationships for modes propagating in graded-index planar waveguides were derived in [16, 17]. A general constraint of the above works has been the assumption on the existence of one [16] or two [17] turning points of the refractive index profile of the planar waveguide. We have derived extended analytical relationships for mode solutions in the graded-index planar waveguide that have no constraints on the number of turning points of the refractive index function.


## 2. Solution of the Helmholtz Equation for a 2D Graded-Index Waveguide

Figure1 gives a schematic representation of the problem. We consider the propagation of the TE-wave in a 2D graded-index medium with the refractive index n=n(x), with the electric field vector directed along the y-axis.Figure 1
Schematic representation of the TE-wave propagating in a 2D graded-index waveguide.The electric field amplitudeEy(x,z) satisfies the following Helmholtz equation [18]:
(1)[∂2∂x2+∂2∂z2+k2n2(x)]Ey(x,z)=0,
where k is the wave number. The expansion of the electromagnetic wave amplitude in terms of the transverse modes of the graded-index medium is
(2)Ey(x,z)=∑n=0∞Cn(x)exp⁡(iβnz),
where βn=kzn is the propagation constant of the nth mode. From (2), it follows that the light field has an axial period T, so that βn=2πn/T. For example, in a graded-index medium with the quadratic index profile, the modes are described by the Hermite-Gauss functions that form a countable basis [19].Substituting (2) in (1) gives an equation of the amplitude of propagating modes in the graded-index medium as
(3)d2Cn(x)dx2+pn(x)Cn(x)=0,
where
(4)pn(x)=k2n2(x)-βn2.By changing variables we have(5)Cn(x)=Cn(0)exp⁡[∫0xfn(ξ)dξ].
Equation (3) can be reduced to a nonlinear differential Whittaker equation [20] for the function fn(x) as follows:
(6)dfn(x)dx+fn2(x)+pn(x)=0.Equation (6) can be solved by the expansion into the Taylor series of the functions fn(x) and pn(x) as follows:
(7)fn(x)=∑m=0∞Cm(n)xm,pn(x)=∑m=0∞pm(n)xm,
where Cm(n) and pm(n) are the unknown and known expansion coefficients of the corresponding functions. Substituting (7) into (6) yields recurrent relations (m>0) for the unknown series terms in (7) as follows:
(8)Cm(n)=-m-1(pm-1(n)+∑s=0m-1Cs(n)Cm-1-s(n)),
where n is integer and C0(n) are indefinite constants. Then, the amplitudes of the propagating modes in a graded-index waveguide are explicitly given by
(9)Cn(x)=Cn(0)exp⁡(∑m=0∞Cm(n)xm+1m+1),
where Cm(n) are derived from the recurrent relationships (8). The coefficients pm(n) entering in (7) are derived from
(10)pm(n)=k2m!dmn2(x)dxm|x=0,m>0,p0(n)=k2n2(0)-βn2,m=0.Finally, the amplitude of the TE-wave propagating in the 2D graded-index waveguide is given by(11)Ey(x,z)=∑n=0∞Cn(0)exp⁡[iβnz+∑m=0∞Cm(n)xm+1m+1].For each mode of (11), there are two indefinite constants, Cn(0) and C0(n). Since the modes of (9) are neither orthogonal nor normalized, so in order for the field of (11) to be expanded in terms of the said modes (with the aim of finding the coefficients Cn(0)), we need to truncate both series in  (11) to finite sums, then solving sets of linear algebraic equations.The constantsC0(n) need to be selected in a special way for each mode. By way of illustration, consider one particular case. Assume the propagating modes in a quadratic-index medium: n2(x)=n02-α2x2. In this case, pn(x)=p0(n)+p1(n)x+p2(n)x2, where p0(n)=k2n02-βn2, p1(n)=0, p2(n)=-k2α2 and pm(n)=0 at m>2. Let C0(n)=0, then, C2(n)=-(C0(n)C1(n))=0 and C1(n)=-p0(n)=βn2-k2n02. For the remaining coefficients to equal zero, Cm(n)=0, at m>2, it will suffice to put the third coefficient to equal zero: C3(n)=-p2(n)/3-(C1(n))2/3=0. Whence, we obtain the following condition on the medium parameter α: k2α2=(k2n02-βn2)2. Assume that βn=kn0/2, then we obtain α=kn02/2. Thus, we can infer that the mode of the quadratic-index waveguide, n2(x)=n02(1-k2n02x2/4), is described by the Gaussian exponential function as follows:
(12)Ey(x,z)=C(0)exp⁡(ikn0z2-k2n02x24).


## 3. General Solution for a Secant-IndexWaveguide

If there is a waveguide with the refractive index defined by a secant function on the transverse coordinate(13)n(x)=n0ch-1(kn0x2),
a particular solution of (1) is given by [21] as follows:
(14)E1y(x,z)=exp⁡(ikn0z2)ch-1(kn0x2).In this case, the Helmholtz equation takes the following form:(15)[∂2∂x2+∂2∂z2+k2n02ch-2(kn0x2)]Ey(x,z)=0.We will seek the general solution of (15) in the following form:
(16)E1y(x,z)=A(x)exp⁡(iγz).Substituting (16) into (15), we obtain
(17)d2A(x)dx2+g(x)A(x)=0,
where
(18)g(x)=k2n02ch-2(kn0x2)-γ2.It has been known [20] that (17) has a general solution
(19)A(x)=A1(x)[C1+C2∫0xA1-2(ξ)dξ],
where A1(x) is a particular solution of (17), and C1, C2 are indefinite constants. In our case, the solution in (14) can be chosen as the particular solution, that is,
(20)A1(x)=ch-1(kn0x2),γ=kn02.Then, the general solution for the secant-index waveguide is given by(21)E1y(x,z)=exp⁡(ikn0z2)ch-1(kn0x2)×{C1+C2kn02[12sh(kn0x2)+kn0x2]},
where
(22)C1=|E1y(x=0,z)|,C2=|dE1y(x=0,z)dx|.From (14), the mode width at half maximum in the secant-index waveguide is
(23)FWHM=ln⁡(3+22)πn02≈0.397λn0,
where λ is the wavelength of light in free space and n0 is the refractive index on the waveguide axis.


## 4. Partial Solution for a Quadratic-Index Waveguide

In Section2 of this paper, we have shown that the propagating mode of a quadratic-index waveguide with definite parameters can be described by the Gaussian exponential function (12). In this section, we shall demonstrate that this remains true of any quadratic-index medium with arbitrary parameters as follows:
(24)n2(x)=n02(1-w2x2),
where w is an arbitrary constant. Then, (1) reads as
(25)[∂2∂x2+∂2∂z2+k2n02(1-w2x2)]Ey(x,z)=0,
and the solution of (25) will be sought for in the following form:
(26)E2y(x,z)=E0exp⁡(ipz-q2x2),
where
(27)q2=wkn02,p=kn01-wkn0.Thus, a particular mode solution of (25) is given by
(28)E2y(x,z)=E0exp⁡(ikn0z1-wkn0-wkn02x2).Note that atw=kn0/2 the solution in (28) coincides with that in (12). At w=kn0/2, it follows from (28) that the Gaussian mode width (diameter) at half maximum is
(29)FWHM=ln⁡4λπn0≈0.38λn0.Comparison of (29) and (23) suggests that the both modes (the secant and the Gaussian mode) have nearly the same width. The effective width of the quadratic-index waveguide is derived from the condition n(x0)=1, being given by
(30)2x0=2n02-1πn0λ.Atn0=1.5, from (30), we find that 2x0≈0.48λ. Thus, the effective width of a glass quadratic-index planar waveguide that can only propagate the Gaussian mode of (28) equals nearly half the free-space wavelength. Note that presently such half-wave waveguides are actively used in applications [22, 23].


## 5. Modeling the Propagation of Light through Graded-Index Microoptics

First experiments on superresolution imaging with superlenses in the optical regime were reported several years ago [3, 5]. Thus, a superresolution of 0.4λ has been achieved in the experiment reported in [3, 5].Theoretically (disregarding the absorption of the material), any degree of superresolution can be achieved with a superlens. In [9], it was shown by modeling that a hyper-lens that would form a magnified subwavelength image in the near field can be implemented as a plane-parallel layer. A 400 nm thick anisotropic slab with dielectric permittivities εx=0.01-i0.01 and εz=-100 was shown to resolve two 3 nm slits (directed along the y-axis) spaced 23 nm apart in a metallic screen with dielectric permittivity ε=1-i104, illuminated by a 700 nm TM wave. The resulting superresolution achieved is 0.05λ.In the following subsections, using a well-known FDTD method, implemented in FullWAVE, we numerically show that the high resolution is also achievable with the aid of 2D graded-index microoptics.


### 5.1. High Resolution Imaging Using a Mikaelian Microlens (ML)

Because of diffraction in uniform space, two nearby point sources of light cease to be resolved at a distance much smaller than the wavelength of incident light. By way of illustration, Figure2(a) depicts profiles of five original coherent Gaussian light sources of width λ/200, spaced λ/50 apart in the plane z=5 nm. Shown in Figure 2(b) is the intensity profile of these sources obtained in the image plane z=30 nm for the wavelength λ=1550 nm. Figure 2 suggests that two nearby light sources cease to be resolved at a distance z approximately equal to the spacing between them (z=λ/50).Light intensity profile near five Gaussian light sources of widthλ/200, spaced λ/50 apart, at different distances: (a) z=5 nm, (b) z=30 nm for λ=1550 nm (intensity is plotted in relative units).
(a)
(b)Shown in Figure3(a) is a numerically simulated image of two 35 nm sources spaced 150 nm apart obtained with a half-pitch ML [24–26]. The refractive index of this 2D secant-graded microlens is
(31)n(x)=n0ch-1(πx2L),
where 2L is the lens length. The lens width 2R is found from the condition n(R)=1. The lens in Figure 3(a) has the axial refractive index n0=2.1, width 2R=1 μm, length 2L=1,144 μm, and operating at wavelength λ=365 nm. These parameters are similar to those described in the experiment in [5].(a) The snapshot of the amplitudeEy(x,z) in the ML. (b) Averaged E-vector intensities of the TE wave calculated directly on the lens “rear” side of the ML (the horizontal line indicates the half-maximum intensity plotted in relative units).
(a)
.
(b)The distance between these sources is0.41λ, and they cannot be resolved by refraction lens. But gradient-index lens (31), which is placed close (20 nm) to sources, works as immersion and allows to increase resolution n0 times, that is, to resolve the two sources separated by distance 0.44λ/n0.Figure3(a) shows the instantaneous distribution of the electric field amplitude for the TE wave (propagating from the bottom upwards). The sources are seen to be imaged on the opposite side of the lens. Figure 3(b) shows a time-averaged intensity profile of the electric field, I(x,z=z0)=|Ey(x,z=z0)|2, directly on the lens “rear” side, that is, at a distance of z0=2L from the lens “front” side. In the experiment, the sources under imaging were put 20 nm before the lens “front” side.From Figure3(b), it is seen that two sources spaced 150 nm apart (with center-to-center distance of 180 nm) can be confidently resolved. The high resolution achieved in Figure 3(b) is 0.41λ. From Figure 3(b), it can also be found that the point source imaged with the half-pitch ML has a width at half maximum of FWHM=100nm=0.27λ. This value is comparable with the diffraction limit FWHM=0.44λ/n=0.21λ in 2D media with index n=2.1. This value of high resolution, 0.41λ, is very close to that reported in [2, 3, 5] (0.40λ). Note that the source image width of 0.27λ well agrees with the minimal width of the propagating mode in a secant-graded waveguide: 0.4λ/n0=0.27λ at n0=1.5 (see (23)).Numerical aperture of the ML isNA=(n02-1)1/2/n0=0.88 for n0=2.1. Therefore, for such lens the focal spot width at half-maximum intensity is FWHM=0.44λ/(n0NA)=0.24λ. This value is slightly less than obtained by simulation: FWHM=0.27λ.To improve the resolution of half-pitch ML, as suggested by (23) and (29), the axial refractive index was changed to n0=3.47 (Si). The other simulation parameters were also changed (Figure 4(a)) to wavelength λ=1 μm, lens width 2R=6 μm, and lens length 2L=4.92 μm. The simulation step along the spatial axes in all examples considered was λ/100. The width of the Gaussian source at the original plane was λ/20. Shown in Figure 4(a) is an instantaneous pattern of the E-vector of the light wave in the half-pitch ML calculated for the instance when the wave has travelled a 220 μm distance from the source. Figure 4(b) depicts an averaged intensity of the TE-wave at the ML output plane (Figure 4(a)). The computation has shown that the intensity of Figure 4(b) corresponds to the source image width at half maximum of FWHM=0.12λ. However, the intensity (or light power density) provides no information as to the proportion of source power propagating along the z-axis. This information can be derived from the projection of Poynting vector onto the optical axis, calculated at the ML output (Figure 4(a)), with the point source found at the lens input, as shown in Figure 5.(a) Instantaneous pattern of the E-vector of the TE wave in the half-pitch ML from a point source found at the front plane (b). An averaged intensity pattern in the lens rear plane (with arbitrary units plotted on they-axis).
(a)
(b)Figure 5
Profile of thez-projection of Poynting vector onto the x-axis calculated at the ML output (Figure 4(a)) with the point source at the input (in relative units).The central-maximum width of the flow of energy propagating along thez-axis in Figure 5 is equal to that in Figure 4(b), being equal to FWHM=0.12λ. The value of the diffraction limit (in 2D case) that can be achieved when focusing light in the homogeneous medium is known to be equal to 0.44λ/n, where n is the refractive index of the homogeneous medium. In our case, n0=3.47, therefore the diffraction limit is FWHM=0.127λ. Numerical aperture for the silicon (n0=3.47) ML is NA=(n02-1)1/2/n0=0.96. Therefore, the focal spot width at half maximum of intensity is FWHM=0.44λ/(n0NA)=0.132λ for such lens. This value is greater than obtained by simulation FWHM=0.12λ.Presumably, this is because when focusing light at the two-medium interface, there is also a contribution to the focus from inhomogeneous surface waves that form interference and diffraction patterns with deep subwavelength features [9, 10]. Figure 5 displays the fact that the surface wave play a part when imaging the point source. In a certain region on the x-axis, the z-projection of Poynting vector is negative, which means that near the output surface of ML there are both outgoing and incoming light waves, with a surface wave propagating along the microlens surface. As light propagates further in free space from the lens output surface, the light spot size quickly increases, becoming equal to the diffraction limit of 0.5λ at a wavelength distance from the surface.Notice that the change of the sign of the Poynting vector projection similar to that shown in Figure5 was earlier reported in [8], being termed as the optical vortex and interpreted as resulting from the interference between the propagating wave and the enhanced surface plasmon. Note, however, that [8] handled a 1D multilayered structure (a 1D photonic crystal).Figure6(a) shows the instantaneous distribution of the E-vector of the light wave in the Mikaelian lens (for the same parameters as in Figure 4(a)) illuminated by two 50 nm sources with center-to-center distance 300 nm, put 10 nm away from the lens bottom plane. Shown in Figure 6(b) is the relative time-averaged distribution of the z-projection of the Poynting vector calculated at the 10 nm distance from the microlens output surface. It can be seen from Figure 6(b) that the resulting superresolution value is 0.3λ, which is smaller than that reported in [2, 4, 5, 11]. Although point source image has width FWHM=0.12λ (Figures 5 and 4(b)), only two point sources, separated by 0.3λ, are surely resolved (Figure 6(b)). This is because images of point sources interfere with each other.(a) Instantaneous distribution of the E-vector of the TE wave in the half-pitch ML of Figure4(a) illuminated by two 50 nm sources with center-to-center distance 300 nm, put 10 nm away from the input plane (bottom horizontal line). (b) Relative time-averaged distribution of the z-projection of Poynting vector calculated at 10 nm distance from the output plane (top horizontal line).
(a)
(b)Let us note that such gradient-index lens can be fabricated as a photonic-crystal lens [27].


### 5.2. High Resolution through a “Fisheye” Microlens

Alongside the ML discussed above, the high resolution can be achieved with other graded-index imaging optical elements for which the refractive index as a function of coordinates has been derived in an explicit analytical form. The possibility of obtaining ideal image of the point source with help of Maxwell’s “fisheye” lens is justified in [28]. One such optical element is represented by Maxwell’s “fisheye” [18], whose 2D refractive index is given in polar coordinates as
(32)n(r)=n0[1+(rR)2]-1,
where n0 is the refractive index at the circle center and R is the element radius. According to (32), the refractive index is halved at r=R. A disadvantage of this element is that the refractive index center-to-edge contrast cannot be larger than two, whereas in the ML the refractive index contrast is a function of the material, for example, for Si ranging from 3.47 to 1. However, with the fisheye having a circular symmetry, any source found on its surface will be perfectly imaged at the diametrically opposite point of its surface.Figure7(a) shows the profile of the E-vector of the TE wave in the 2D fisheye microlens that has two nearby point sources on its surface. The simulation was performed for the refractive index at the lens center n0=3.47, lens radius R=2.5 μm, wavelength λ=1 μm, and two 0.05λ sources separated by a 440 nm distance, which is 0.44λ. Shown in Figure 7(b) is a time-averaged intensity pattern in the image plane (top horizontal line in Figure 7(a)). The two sources are seen to be resolved (the resolution being 0.44λ by the 20% Rayleigh criterion). Thus, the lens superresolution of 0.44λ is insignificantly beyond the diffraction limit of 0.5λ.(a) An instantaneous distribution of the E-vector of the TE wave in the 2D fisheye microlens simulated in FullWAVE, when illuminated from bottom by two point sources (b). Relative time-averaged intensity distribution in the image plane.
(a)
(b)When the microlens of Figure7(a) is illuminated by a solitary source, the resulting intensity distribution is shown in Figure 8. The image size at half maximum is FWHM=0.3λ. From the comparison of Figures 8 and 4(b),  we can infer that the fisheye forms a (about two-times) wider image of the point source when compared with the ML at similar parameters. Note that the magnitudes of high resolution obtained, 0.4λ (Figure 7(b)) and 0.3λ (Figure 8), are comparable with those reported in [2, 4, 5, 8, 11, 12].Figure 8
Relative time-averaged intensity pattern of the electric field in the image plane of the eyefish microlens of Figure7(a) when illuminated by a single point source found on its surface.Although point source image has widthFWHM=0.3λ (Figure 8), two point sources must be separated by 0.4λ (Figure 7(b)) distance in order to resolve them by Rayleigh criterion. This is because images of point sources interfere with each other.


## 5.1. High Resolution Imaging Using a Mikaelian Microlens (ML)

Because of diffraction in uniform space, two nearby point sources of light cease to be resolved at a distance much smaller than the wavelength of incident light. By way of illustration, Figure2(a) depicts profiles of five original coherent Gaussian light sources of width λ/200, spaced λ/50 apart in the plane z=5 nm. Shown in Figure 2(b) is the intensity profile of these sources obtained in the image plane z=30 nm for the wavelength λ=1550 nm. Figure 2 suggests that two nearby light sources cease to be resolved at a distance z approximately equal to the spacing between them (z=λ/50).Light intensity profile near five Gaussian light sources of widthλ/200, spaced λ/50 apart, at different distances: (a) z=5 nm, (b) z=30 nm for λ=1550 nm (intensity is plotted in relative units).
(a)
(b)Shown in Figure3(a) is a numerically simulated image of two 35 nm sources spaced 150 nm apart obtained with a half-pitch ML [24–26]. The refractive index of this 2D secant-graded microlens is
(31)n(x)=n0ch-1(πx2L),
where 2L is the lens length. The lens width 2R is found from the condition n(R)=1. The lens in Figure 3(a) has the axial refractive index n0=2.1, width 2R=1 μm, length 2L=1,144 μm, and operating at wavelength λ=365 nm. These parameters are similar to those described in the experiment in [5].(a) The snapshot of the amplitudeEy(x,z) in the ML. (b) Averaged E-vector intensities of the TE wave calculated directly on the lens “rear” side of the ML (the horizontal line indicates the half-maximum intensity plotted in relative units).
(a)
.
(b)The distance between these sources is0.41λ, and they cannot be resolved by refraction lens. But gradient-index lens (31), which is placed close (20 nm) to sources, works as immersion and allows to increase resolution n0 times, that is, to resolve the two sources separated by distance 0.44λ/n0.Figure3(a) shows the instantaneous distribution of the electric field amplitude for the TE wave (propagating from the bottom upwards). The sources are seen to be imaged on the opposite side of the lens. Figure 3(b) shows a time-averaged intensity profile of the electric field, I(x,z=z0)=|Ey(x,z=z0)|2, directly on the lens “rear” side, that is, at a distance of z0=2L from the lens “front” side. In the experiment, the sources under imaging were put 20 nm before the lens “front” side.From Figure3(b), it is seen that two sources spaced 150 nm apart (with center-to-center distance of 180 nm) can be confidently resolved. The high resolution achieved in Figure 3(b) is 0.41λ. From Figure 3(b), it can also be found that the point source imaged with the half-pitch ML has a width at half maximum of FWHM=100nm=0.27λ. This value is comparable with the diffraction limit FWHM=0.44λ/n=0.21λ in 2D media with index n=2.1. This value of high resolution, 0.41λ, is very close to that reported in [2, 3, 5] (0.40λ). Note that the source image width of 0.27λ well agrees with the minimal width of the propagating mode in a secant-graded waveguide: 0.4λ/n0=0.27λ at n0=1.5 (see (23)).Numerical aperture of the ML isNA=(n02-1)1/2/n0=0.88 for n0=2.1. Therefore, for such lens the focal spot width at half-maximum intensity is FWHM=0.44λ/(n0NA)=0.24λ. This value is slightly less than obtained by simulation: FWHM=0.27λ.To improve the resolution of half-pitch ML, as suggested by (23) and (29), the axial refractive index was changed to n0=3.47 (Si). The other simulation parameters were also changed (Figure 4(a)) to wavelength λ=1 μm, lens width 2R=6 μm, and lens length 2L=4.92 μm. The simulation step along the spatial axes in all examples considered was λ/100. The width of the Gaussian source at the original plane was λ/20. Shown in Figure 4(a) is an instantaneous pattern of the E-vector of the light wave in the half-pitch ML calculated for the instance when the wave has travelled a 220 μm distance from the source. Figure 4(b) depicts an averaged intensity of the TE-wave at the ML output plane (Figure 4(a)). The computation has shown that the intensity of Figure 4(b) corresponds to the source image width at half maximum of FWHM=0.12λ. However, the intensity (or light power density) provides no information as to the proportion of source power propagating along the z-axis. This information can be derived from the projection of Poynting vector onto the optical axis, calculated at the ML output (Figure 4(a)), with the point source found at the lens input, as shown in Figure 5.(a) Instantaneous pattern of the E-vector of the TE wave in the half-pitch ML from a point source found at the front plane (b). An averaged intensity pattern in the lens rear plane (with arbitrary units plotted on they-axis).
(a)
(b)Figure 5
Profile of thez-projection of Poynting vector onto the x-axis calculated at the ML output (Figure 4(a)) with the point source at the input (in relative units).The central-maximum width of the flow of energy propagating along thez-axis in Figure 5 is equal to that in Figure 4(b), being equal to FWHM=0.12λ. The value of the diffraction limit (in 2D case) that can be achieved when focusing light in the homogeneous medium is known to be equal to 0.44λ/n, where n is the refractive index of the homogeneous medium. In our case, n0=3.47, therefore the diffraction limit is FWHM=0.127λ. Numerical aperture for the silicon (n0=3.47) ML is NA=(n02-1)1/2/n0=0.96. Therefore, the focal spot width at half maximum of intensity is FWHM=0.44λ/(n0NA)=0.132λ for such lens. This value is greater than obtained by simulation FWHM=0.12λ.Presumably, this is because when focusing light at the two-medium interface, there is also a contribution to the focus from inhomogeneous surface waves that form interference and diffraction patterns with deep subwavelength features [9, 10]. Figure 5 displays the fact that the surface wave play a part when imaging the point source. In a certain region on the x-axis, the z-projection of Poynting vector is negative, which means that near the output surface of ML there are both outgoing and incoming light waves, with a surface wave propagating along the microlens surface. As light propagates further in free space from the lens output surface, the light spot size quickly increases, becoming equal to the diffraction limit of 0.5λ at a wavelength distance from the surface.Notice that the change of the sign of the Poynting vector projection similar to that shown in Figure5 was earlier reported in [8], being termed as the optical vortex and interpreted as resulting from the interference between the propagating wave and the enhanced surface plasmon. Note, however, that [8] handled a 1D multilayered structure (a 1D photonic crystal).Figure6(a) shows the instantaneous distribution of the E-vector of the light wave in the Mikaelian lens (for the same parameters as in Figure 4(a)) illuminated by two 50 nm sources with center-to-center distance 300 nm, put 10 nm away from the lens bottom plane. Shown in Figure 6(b) is the relative time-averaged distribution of the z-projection of the Poynting vector calculated at the 10 nm distance from the microlens output surface. It can be seen from Figure 6(b) that the resulting superresolution value is 0.3λ, which is smaller than that reported in [2, 4, 5, 11]. Although point source image has width FWHM=0.12λ (Figures 5 and 4(b)), only two point sources, separated by 0.3λ, are surely resolved (Figure 6(b)). This is because images of point sources interfere with each other.(a) Instantaneous distribution of the E-vector of the TE wave in the half-pitch ML of Figure4(a) illuminated by two 50 nm sources with center-to-center distance 300 nm, put 10 nm away from the input plane (bottom horizontal line). (b) Relative time-averaged distribution of the z-projection of Poynting vector calculated at 10 nm distance from the output plane (top horizontal line).
(a)
(b)Let us note that such gradient-index lens can be fabricated as a photonic-crystal lens [27].


## 5.2. High Resolution through a “Fisheye” Microlens

Alongside the ML discussed above, the high resolution can be achieved with other graded-index imaging optical elements for which the refractive index as a function of coordinates has been derived in an explicit analytical form. The possibility of obtaining ideal image of the point source with help of Maxwell’s “fisheye” lens is justified in [28]. One such optical element is represented by Maxwell’s “fisheye” [18], whose 2D refractive index is given in polar coordinates as
(32)n(r)=n0[1+(rR)2]-1,
where n0 is the refractive index at the circle center and R is the element radius. According to (32), the refractive index is halved at r=R. A disadvantage of this element is that the refractive index center-to-edge contrast cannot be larger than two, whereas in the ML the refractive index contrast is a function of the material, for example, for Si ranging from 3.47 to 1. However, with the fisheye having a circular symmetry, any source found on its surface will be perfectly imaged at the diametrically opposite point of its surface.Figure7(a) shows the profile of the E-vector of the TE wave in the 2D fisheye microlens that has two nearby point sources on its surface. The simulation was performed for the refractive index at the lens center n0=3.47, lens radius R=2.5 μm, wavelength λ=1 μm, and two 0.05λ sources separated by a 440 nm distance, which is 0.44λ. Shown in Figure 7(b) is a time-averaged intensity pattern in the image plane (top horizontal line in Figure 7(a)). The two sources are seen to be resolved (the resolution being 0.44λ by the 20% Rayleigh criterion). Thus, the lens superresolution of 0.44λ is insignificantly beyond the diffraction limit of 0.5λ.(a) An instantaneous distribution of the E-vector of the TE wave in the 2D fisheye microlens simulated in FullWAVE, when illuminated from bottom by two point sources (b). Relative time-averaged intensity distribution in the image plane.
(a)
(b)When the microlens of Figure7(a) is illuminated by a solitary source, the resulting intensity distribution is shown in Figure 8. The image size at half maximum is FWHM=0.3λ. From the comparison of Figures 8 and 4(b),  we can infer that the fisheye forms a (about two-times) wider image of the point source when compared with the ML at similar parameters. Note that the magnitudes of high resolution obtained, 0.4λ (Figure 7(b)) and 0.3λ (Figure 8), are comparable with those reported in [2, 4, 5, 8, 11, 12].Figure 8
Relative time-averaged intensity pattern of the electric field in the image plane of the eyefish microlens of Figure7(a) when illuminated by a single point source found on its surface.Although point source image has widthFWHM=0.3λ (Figure 8), two point sources must be separated by 0.4λ (Figure 7(b)) distance in order to resolve them by Rayleigh criterion. This is because images of point sources interfere with each other.


## 6. Conclusion

Thus, in this work we have(i)
derived mode solutions of the Helmholtz equation for an arbitrary graded-index planar waveguide, with the mode amplitude represented as the exponential of the Taylor series whose coefficients are deduced from the recurrent relations ((8), (9), and (11));(ii)
shown that the minimal mode width in quadratic- and secant-index planar waveguides amounts to 0.4 of the wavelength in free space divided by the refractive index on the waveguide axis ((23) and (29));(iii)
shown by modeling in FullWAVE that graded-index half-pitch ML and fisheye lenses are capable of high-resolution imaging (Figures6 and 7);(iv)
shown that a point source is imaged through the 2D half-pitch ML as a near-surface light spot of sizeFWHM=0.12λ (λ is wavelength in free space), which is close to the diffraction limit for silicon (n=3.47) in 2D media FWHM=0.44λ/n=0.127λ and smaller than values reported in [8, 11, 12] (Figures 4 and 5);(v)
also shown that the half-pitch ML is able to resolve at half-maximum two nearby sources placed 0.3λ apart, which is beyond the diffraction limit in free space of 0.44λ for 2D case and smaller than values reported in [2, 4, 5, 11] (Figure 6).



---

*Source: 647165-2012-07-11.xml*