Abstract

Superlens imaging system in nanolithography can be regarded as a cascade of two F-P cavities, i.e., a superlens cavity and a dielectric cavity between superlens and introduced mask of high loss, and the transfer function of system is obtained by considering multiple reflections inside the two cavities. For the range of wavevector of interest, the typical high peak of transmission coefficient of superlens coincides with a local minimum of transmission coefficient of dielectric cavity. The peak of transfer function of system corresponds to the peak of transmission coefficient of dielectric cavity. Thin superlens imaging system in nanolithography is analyzed based on transfer function, which can be flattened by simply tuning transmission coefficient of dielectric cavity and superlens cavity. The results are further validated by Finite Element Method (FEM) simulations.

1. Introduction

The metal planar superlens, with a negative permittivity at optical frequency and a positive permeability, was proposed by Pendry in 2000 [1] as an alternative to negative index media (NIM) by Veselago in 1967 [2] to break the diffraction limit. The principle of superlens is to compensate the exponential decay of the evanescent field away from the object by amplifying the evanescent waves through surface plasmons (SP) resonances. Since then, many research efforts had been devoted into superlens imaging [315], including theoretical models, numerical simulations, experimental demonstrations, and applications.

The silver superlens was demonstrated experimentally with the resolution of one-sixth of the illumination wavelength by Zhang et al. in 2005 [3]. In their experiment, 35 nm thick silver slab was used because it gives the optimum transfer function. Thinner sliver slabs show higher but narrow enhancement bands, and thicker slabs show smaller enhancements in reference to zero-order transmission [4]. In order to eliminate the sharp peak of transfer function, Sheng et al. proposed designing the metallic superlens close to the cutoff condition of the long-range SP mode to balance the amplification by the SP resonance and the flatness of the transfer function [5, 6]. Moore et al. suggested a performance window for superlens with total thickness range from 120 nm to 140 nm in order to get a flatter transfer function [7]. Those were reported by considering transfer functions of the imaging system with three-layer structure, a superlens sandwiched by two semi-infinite dielectrics [47].

Normally, when dealing with the imaging system, a perfectly absorbing, thin screen with slit(s) has been widely adopted. The transmission of the idealized screen will be 1 in the slit or 0 otherwise. While in real nanolithography applications, there exists a metallic object mask. The metallic mask itself exists mainly for constructing a perfect object function in the object plane just behind the mask for the imaging system. To obtain a good enough object function, which means as close as the ideal case, i.e., 1 in the slit and 0 otherwise, a thick mask with intrinsic high loss, such as Chrome (Cr), is usually applied. Due to the introduction of mask, a dielectric cavity is naturally formed between the mask and superlens. Blaikie et al. approximately considered the neglected recursive reflections in the dielectric cavity between the superlens and the mask and modified transfer function by an improved transfer-matrix model [8]. Sheng et al. further took possible SPR by mask itself in some cases into account and gave a more general model with its transfer function, which is optimized by genetic algorithm (GA) [9].

For the superlens imaging system in nanolithography, we regarded it as two cascaded F-P cavities, a superlens cavity and a dielectric cavity between superlens and introduced mask of high loss, and the transfer function of system is obtained by considering multiple reflections inside the two cavities [9, 10]. We studied the transfer function of system and the transmission coefficient of two cavities and then revealed some relations among them. It is found that the peak shown in transmission coefficient of superlens always corresponds to the local minimum of transmission coefficient of dielectric cavity. Moreover, the peak of transfer function of system coincides with transmission coefficient of dielectric cavity instead of that of superlens. Based on these, we propose that, by simply tuning the transmission coefficient of superlens cavity and dielectric cavity, a very thin metal superlens imaging system can show well-balanced transfer function and thus produces improved image, which is validated by FEM simulations.

2. Metal Superlens Imaging System

Consider an imaging system for nanolithography with a metal superlens of permittivity and a thick metallic mask of permittivity with high intrinsic loss to construct an object function, as shown in Figure 1. For the purpose of the image analysis, herein a slit of nanosize is perforated through the object mask. A dielectric cavity, with permittivity , is formed between the superlens and the mask. The illuminating light of TM polarization is incident normally to the object mask from the semi-infinite dielectric medium . The image plane is placed in a semi-infinite dielectric medium at a distance from the superlens. The incident plane wave of TM polarization was launched in dielectric medium and then encountered the object mask . In the case where the lithographic metal mask layer is thick and of high intrinsic loss, the incident waves will be dissipated in the mask medium and can only pass through the slit. At the exit of the slit, the waves, labelled as source A, will be scattered in the medium launching the Surface Plasmon Polaritons (SPPs) along interfaces of the superlens. Both the homogeneous propagating waves and inhomogeneous SPP waves will reach the image plane and contribute to the image in dielectric medium . The evanescent waves will simply undergo an exponential delay in the semi-infinite dielectric medium . Therefore, the imaging system for source A shown in Figure 1 can be regarded as two cascaded Fabry-Perot (F-P) cavities. The first cavity is the dielectric layer with two metal/dielectric interfaces, and , respectively. The second cavity is the metal superlens layer with two interfaces, and , respectively. The transfer function of the imaging system can be computed as the product of transmission coefficients of the two cavities, which are obtained by considering multiple reflections inside the two cavities. We analyze the mathematical expressions of the transmission coefficients corresponding to the two cavities in order to get a flat transfer function of the imaging system.

2.1. Transmission Coefficient of Superlens Cavity

The superlens cavity of permittivity is between two metal/dielectric interfaces of and . Its transmission coefficient can be calculated by considering resonance of multiple reflected fields in the superlens cavity as [1, 5]where the Fresnel reflection and transmission coefficients from medium i to medium j, with subindices i, j = o, d, s, and q, are and with ; the propagation factor , describing the phase change of the propagating waves with along the distance and the exponent decay in amplitude of the evanescent waves over , respectively.

According to (1) and for a silver superlens with variable thickness, the amplitude of transmission coefficient as a function of normalized wavevector is depicted in Figure 2(a). Typically for a thin superlens of less than 30 nm in thickness a sharp and high peak is shown at a low spatial frequency slightly larger than and a relative broad and low peak is at a high spatial frequency . When increasing the thickness only one peak will appear at a spatial frequency between spatial frequency and . For large thickness of superlens, the amplitude of transmission coefficient becomes low and flat due to high loss in the thick metal superlens, as shown in Figure 3(a).

According to Maxwell’s equations and the boundary conditions in the Dielectric-Metal-Dielectric (DMD) waveguide structure, the dispersion relation can be obtained and solved numerically [5]. The effective indices of waveguide modes are plotted for different thickness of Ag superlens in Figure 2(b). Two modes, long-range surface plasmon (LRSP) mode and short-range surface plasmon (SRSP) mode, are supported for DMD waveguide of a thin superlens. The two peaks correspond to such two excited SP modes: a LRSP mode associated with a narrow peak located at lower spatial frequencies and a SRSP mode associated with a broad peak at higher spatial frequencies . When the thickness of the metal superlens increases, location of the peak at low frequency (LRSP mode) tends to shift to higher frequency, and location of the peak at high frequency (SRSP mode) tends to lower frequency. At a certain thickness of the superlens, only one mode will be supported, which accounts for only one peak of the transmission coefficient.

If the imaging system is considered as a superlens between two semi-infinite dielectrics without presence of the metal lithographic mask, as that in [47], then the very thin superlens is usually inappropriate for imaging due to the existence of sharp peak of resulting from excitation of LRSP mode. The LRSP mode can be cutoff [5] if is met, which means the reflection of metal layer is around 0.

2.2. Transmission Coefficient of Dielectric Cavity

The dielectric slab of is sandwiched by the metal superlens and the metal mask, constituting a dielectric F-P cavity. Its transmission coefficient can be calculated [9, 10] by taking into account the recursive reflections in the cavity asFor a thick mask with high loss used in the superlens image system for nanolithography, ; then (2) can be further written as For evanescent waves with , is real. The denominator in (3) can be evaluated, based on and , asThus, when the amplitude of transmission coefficient of the superlens shows a very high peak associated with the excitation of LRSP mode at a low spatial frequency slightly larger than , the amplitude of transmission coefficient of the dielectric cavity will fall to a local minimum. Meanwhile, will reach a maximum when the absolute value of denominator of (3) approaches 0; i.e.,Besides, it is interesting to notice that if in (3), then the transmission coefficient , implying an exponential decay of the evanescent waves in the dielectric cavity, just as that in the free space without any resonances, and is independent of the cavity thickness.

2.3. Transfer Function of Imaging System

As the imaging system for source A is a cascade of two F-P cavities, the total transfer function is calculated byand when condition (5) is satisfied, and will both reach their maximums, according to (3) and (6), at the same spatial frequency . This reveals that the peak of the total transfer function for source A in the imaging system corresponds to the peak of transmission coefficient of dielectric cavity at instead of that of superlens at . Actually, the value of at is always balanced as the transmission coefficient of the superlens can be suppressed by the local minimum of at . The transfer function for source A is defined as the transmission coefficients of the system described in (6) as a function of the spatial spectral frequency. According to SPP waveguide theory, the LRSP mode is cutoff when the propagation constant is purely real, which occurs when the nature of the mode changes from attenuating () to growing () with the propagation along the metal layer [5, 6]. If the cutoff condition for LRSP mode, i.e., , is fulfilled, the transfer function . In this case, a relative flat transfer function can be expected when considering together the well-controlled, proper amplification of SPP by superlens and an exponential decay of amplitude of evanescent waves in dielectric cavity.

3. Imaging by Thin Superlens Imaging System

We first consider an Ag superlens imaging system with typical parameters: the metallic mask is chrome with permittivity and thickness of = 50nm and the silver slab with permittivity with incident wavelength of 365 nm and thickness = 30 nm; the dielectric layer between superlens and mask has the permittivity = 2.34 and thickness of = 40 nm, and the dielectric = 2.34, as shown in Figure 1. From Figure 3(a), the sharp peak of corresponds to the local minimum of at , which is the effective index of LRSP mode excited in DMD waveguide. The broad peak of at is shown because of the excitation of SRSP mode. The value of transfer function at is balanced by the multiplication of and , which indicates that we cannot judge the imaging system only by transmission coefficient of superlens. The peak of transfer function locates at , which is the same as that of transmission coefficient of dielectric cavity . We change the thickness of Ag to = 15 nm and thickness of dielectric =10 nm, which form a thinner superlens imaging system. The transmission coefficient of 15 nm thick superlens has a shaper and higher peak, as shown in Figure 3(b), than that of superlens of 30 nm thickness in Figure 3(a). Still, the peak of corresponds to the local minimum of at , and thus the value of at is balanced in the same manner. The peak of transfer function and transmission coefficient of the dielectric cavity is obtained around , as depicted in Figure 3(b). The very thin superlens imaging system can show a flatter transfer function over a broader range of spatial frequency by simply tuning the transmission coefficient of superlens cavity and dielectric cavity.

For further validations, 2D Finite Element Method (FEM) by commercial software package COMSOL Multiphysics 4.4 was employed to simulate the imaging system in (x, y) plane, as shown in Figure 4. A TM (, , ) plane wave with wavelength of 365 nm was launched along y direction from the input port on the top boundary in the FEM simulations. A perfectly matched layer (PML) was applied at the bottom boundary and periodical boundary conditions (PBC) were imposed on the left and right sides. A chrome mask with two slits of 40 nm was used. The centre-to-centre distance of two slits was 120 nm. The lateral length along x direction of each layer was 1 um. To clearly show the Hz fields after the mask of high loss, the incident waves before and in the mask and the mask itself are removed in Figures 4(a) and 4(c). The image of the two spaced slits is recorded at the imaging plane placed at a distance of 10 nm away from the superlens, as shown in Figures 4(b) and 4(d). Obviously, the image in Figure 4(d) by thin superlens system with transfer function in Figure 3(b) is better than the one in Figure 4(b) with the transfer function shown in Figure 3(a), because the side lobes of the image in Figure 4(d) were largely suppressed compared with that in Figure 4(b). However, the amplitude of the unwanted centre peak is still high, as shown in Figures 4(b) and 4(d).

The imaging performance of thin superlens imaging system can be further improved by approaching the LRSP mode cutoff condition as we demonstrated. Change the dielectric = 2.34 to = 4.08, and keep all others unchanged as those in Figure 4(c). The transfer function and transmission coefficient and are plotted in Figure 5(a), from which we can see that the high peak of associated with LRSP mode is mainly removed, and approximately follows an exponential decay for evanescent waves, and thus the transfer function is balanced to show flatness. The Hz field by a FEM simulation in this superlens system is shown in Figure 5(b) and the image is notably improved as shown in Figure 5(c), in which both the side lobes and the centre peak are greatly suppressed compared with those in Figures 4(b) and 4(d).

4. Conclusions

We investigated the imaging performance of the metal planar superlens for nanolithography, considering the superlens imaging system as a cascade of two F-P cavities: a superlens cavity and a dielectric cavity between superlens and introduced mask of high loss. By analyzing the transmission coefficient of superlens and dielectric cavity and the transfer function of whole system, we found that the peak of the transmission coefficient of the superlens always coincides with the local minimum of that of the dielectric cavity, and the peak of the whole transfer function of system actually corresponds to that of dielectric cavity instead of superlens. Then, we show that very thin superlens imaging system, which was usually believed to be improper for superlens imaging, has an improved imaging performance by simply tuning the transmission coefficient of superlens cavity and dielectric cavity. All the results are analyzed based on transfer functions and further confirmed by FEM simulations.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors acknowledge the support from the Natural Science and Engineering Research Council (NSERC) of Canada, Discovery Grant, to this research program.