# Thin Metal Superlens Imaging in Nanolithography

Authors: Jing Wang; Yunlong Sheng
Journal: International Journal of Optics
 (2019)

---

## 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.


---

## Body

## 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 [3–15], 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 [4–7].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εs and a thick metallic mask of permittivity εo 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 εd, 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 εp. The image plane is placed in a semi-infinite dielectric medium εq at a distance dq from the superlens. The incident plane wave of TM polarization was launched in dielectric medium εp and then encountered the object mask εo. 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 εd 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 εq. The evanescent waves will simply undergo an exponential delay in the semi-infinite dielectric medium εq. 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 εd with two metal/dielectric interfaces, εd/εo and εd/εs, respectively. The second cavity is the metal superlens layer with two interfaces, εs/εd and εs/εq, 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.Figure 1
Metal superlens imaging system for nanolithography.


### 2.1. Transmission Coefficient of Superlens Cavity

The superlens cavity of permittivityεs is between two metal/dielectric interfaces of εs/εd and εs/εq. Its transmission coefficient can be calculated by considering resonance of multiple reflected fields in the superlens cavity as [1, 5](1)τs=estdstsq1-es2rsdrsq,where the Fresnel reflection and transmission coefficients from mediumi to mediumj, with subindicesi, j =o, d, s,andq, are rij=εjkzi-εikzj/εjkzi+εikzj and tij=2εjkzi/εjkzi+εikzj with kzi=εik02-kx2; the propagation factor ei=expikzidi, describing the phase change of the propagating waves with kx2<εik02 along the distance di and the exponent decay in amplitude of the evanescent waves kx2>εik02 over di, respectively.According to (1) and for a silver superlens with variable thickness, the amplitude of transmission coefficient τs as a function of normalized wavevector kx0=kx/εdk0 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 kx01 slightly larger than kx0=1 and a relative broad and low peak is at a high spatial frequency kx02. When increasing the thickness only one peak will appear at a spatial frequency between spatial frequency kx01 and kx02. 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).Figure 2
For Ag superlens of thickness ranging from 5 nm to 120 nm sandwiched by two dielectric mediumsεd=εq=2.34, (a) the amplitude of transmission coefficient τs<1 for propagating waves with kx0<1, while for evanescent waves with kx0>1, it can be amplified. (b) Mode effective index of DMD waveguide.
(a) (b)Figure 3
The amplitude of transmission efficient of superlens and dielectric and transfer function of the system as a function of normalized wavevector. (a)dd = 40nm, ds = 30 nm, (b) dd = 10nm, and ds = 15nm.
(a) (b)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 kx01 and a SRSP mode associated with a broad peak at higher spatial frequencies kx02. 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 [4–7], then the very thin superlens is usually inappropriate for imaging due to the existence of sharp peak of τs resulting from excitation of LRSP mode. The LRSP mode can be cutoff [5] if rds+es2rsq≈0 is met, which means the reflection of metal layer is around 0.


### 2.2. Transmission Coefficient of Dielectric Cavity

The dielectric slab ofεd 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 as(2)τd=ed1-ed2rds+es2rsq/1-es2rsdrsqrdo+eo2rop/1-eo2roprod.For a thick mask with high loss used in the superlens image system for nanolithography, rop=0; then (2) can be further written as (3)τd=11/ed-edrds+es2rsq/1-es2rsdrsqrdo=11/ed-rds+es2rsqrdoed/estdstsqτs.For evanescent waves with kx2>εdk02, ed=expiεdk02-kx2dd=exp-kx2-εdk02dd is real. The denominator in (3) can be evaluated, based on a-b≤a±b≤a+b and ab=ab, as(4)rds+es2rsqrdoedestdstsqτs-1ed≤1ed-rds+es2rsqrdoedestdstsqτs≤rds+es2rsqrdoedestdstsqτs+1ed.Thus, when the amplitude of transmission coefficient of the superlens τs shows a very high peak associated with the excitation of LRSP mode at a low spatial frequency kx01 slightly larger than kx0=1, the amplitude of transmission coefficient of the dielectric cavity τd will fall to a local minimum. Meanwhile, τd will reach a maximum when the absolute value of denominator of (3) approaches 0; i.e.,(5)1ed-rds+es2rsqrdoedestdstsqτs→0.Besides, it is interesting to notice that if rds+es2rsq≈0 in (3), then the transmission coefficient τd≈ed, 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τA is calculated by(6)τA=τd·τs=11/edτs-rds+es2rsqrdoed/estdstsq,and when condition (5) is satisfied, τA and τd will both reach their maximums, according to (3) and (6), at the same spatial frequency kx03. This reveals that the peak of the total transfer function τA for source A in the imaging system corresponds to the peak of transmission coefficient of dielectric cavity τd at kx03 instead of that of superlens τs at kx01. Actually, the value of τA at kx01 is always balanced as the transmission coefficient of the superlens τs can be suppressed by the local minimum of τd at kx01. 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 (Imβ>0) to growing (Imβ<0) with the propagation along the metal layer [5, 6]. If the cutoff condition for LRSP mode, i.e., rds+es2rsq≈0, is fulfilled, the transfer function τA=τd·τs≈ed·τs. 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.


## 2.1. Transmission Coefficient of Superlens Cavity

The superlens cavity of permittivityεs is between two metal/dielectric interfaces of εs/εd and εs/εq. Its transmission coefficient can be calculated by considering resonance of multiple reflected fields in the superlens cavity as [1, 5](1)τs=estdstsq1-es2rsdrsq,where the Fresnel reflection and transmission coefficients from mediumi to mediumj, with subindicesi, j =o, d, s,andq, are rij=εjkzi-εikzj/εjkzi+εikzj and tij=2εjkzi/εjkzi+εikzj with kzi=εik02-kx2; the propagation factor ei=expikzidi, describing the phase change of the propagating waves with kx2<εik02 along the distance di and the exponent decay in amplitude of the evanescent waves kx2>εik02 over di, respectively.According to (1) and for a silver superlens with variable thickness, the amplitude of transmission coefficient τs as a function of normalized wavevector kx0=kx/εdk0 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 kx01 slightly larger than kx0=1 and a relative broad and low peak is at a high spatial frequency kx02. When increasing the thickness only one peak will appear at a spatial frequency between spatial frequency kx01 and kx02. 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).Figure 2
For Ag superlens of thickness ranging from 5 nm to 120 nm sandwiched by two dielectric mediumsεd=εq=2.34, (a) the amplitude of transmission coefficient τs<1 for propagating waves with kx0<1, while for evanescent waves with kx0>1, it can be amplified. (b) Mode effective index of DMD waveguide.
(a) (b)Figure 3
The amplitude of transmission efficient of superlens and dielectric and transfer function of the system as a function of normalized wavevector. (a)dd = 40nm, ds = 30 nm, (b) dd = 10nm, and ds = 15nm.
(a) (b)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 kx01 and a SRSP mode associated with a broad peak at higher spatial frequencies kx02. 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 [4–7], then the very thin superlens is usually inappropriate for imaging due to the existence of sharp peak of τs resulting from excitation of LRSP mode. The LRSP mode can be cutoff [5] if rds+es2rsq≈0 is met, which means the reflection of metal layer is around 0.


## 2.2. Transmission Coefficient of Dielectric Cavity

The dielectric slab ofεd 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 as(2)τd=ed1-ed2rds+es2rsq/1-es2rsdrsqrdo+eo2rop/1-eo2roprod.For a thick mask with high loss used in the superlens image system for nanolithography, rop=0; then (2) can be further written as (3)τd=11/ed-edrds+es2rsq/1-es2rsdrsqrdo=11/ed-rds+es2rsqrdoed/estdstsqτs.For evanescent waves with kx2>εdk02, ed=expiεdk02-kx2dd=exp-kx2-εdk02dd is real. The denominator in (3) can be evaluated, based on a-b≤a±b≤a+b and ab=ab, as(4)rds+es2rsqrdoedestdstsqτs-1ed≤1ed-rds+es2rsqrdoedestdstsqτs≤rds+es2rsqrdoedestdstsqτs+1ed.Thus, when the amplitude of transmission coefficient of the superlens τs shows a very high peak associated with the excitation of LRSP mode at a low spatial frequency kx01 slightly larger than kx0=1, the amplitude of transmission coefficient of the dielectric cavity τd will fall to a local minimum. Meanwhile, τd will reach a maximum when the absolute value of denominator of (3) approaches 0; i.e.,(5)1ed-rds+es2rsqrdoedestdstsqτs→0.Besides, it is interesting to notice that if rds+es2rsq≈0 in (3), then the transmission coefficient τd≈ed, 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τA is calculated by(6)τA=τd·τs=11/edτs-rds+es2rsqrdoed/estdstsq,and when condition (5) is satisfied, τA and τd will both reach their maximums, according to (3) and (6), at the same spatial frequency kx03. This reveals that the peak of the total transfer function τA for source A in the imaging system corresponds to the peak of transmission coefficient of dielectric cavity τd at kx03 instead of that of superlens τs at kx01. Actually, the value of τA at kx01 is always balanced as the transmission coefficient of the superlens τs can be suppressed by the local minimum of τd at kx01. 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 (Imβ>0) to growing (Imβ<0) with the propagation along the metal layer [5, 6]. If the cutoff condition for LRSP mode, i.e., rds+es2rsq≈0, is fulfilled, the transfer function τA=τd·τs≈ed·τs. 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εo=-8.55+i8.96 and thickness of do = 50nm and the silver slab with permittivity εs=-2.6+i0.25 with incident wavelength of 365 nm and thickness ds = 30 nm; the dielectric layer between superlens and mask has the permittivity εd = 2.34 and thickness of dd = 40 nm, and the dielectric εq= 2.34, as shown in Figure 1. From Figure 3(a), the sharp peak of τs corresponds to the local minimum of τd at kx01≈1.3, which is the effective index of LRSP mode excited in DMD waveguide. The broad peak of τs at kx02≈3.6 is shown because of the excitation of SRSP mode. The value of transfer function τA at kx01≈1.3 is balanced by the multiplication of τs and τd, which indicates that we cannot judge the imaging system only by transmission coefficient of superlens. The peak of transfer function τA locates at kx03≈1.6, which is the same as that of transmission coefficient of dielectric cavity τd. We change the thickness of Ag to ds = 15 nm and thickness of dielectric dd =10 nm, which form a thinner superlens imaging system. The transmission coefficient τs 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 τs corresponds to the local minimum of τd at kx01≈1.1, and thus the value of τA at kx01≈1.1 is balanced in the same manner. The peak of transfer function τA and transmission coefficient of the dielectric cavity τd is obtained around kx03≈1.4, 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 Figure4. A TM (Ex, Ey, Hz) 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).Figure 4
Simulation results of the superlens imaging system by FEM. (a) Hz field distributions of system withdd = 40nm and ds = 30 nm; (b) image of two-slit object in the imaging plane in (a); (c) Hz field distributions of system with dd = 10nm and ds = 15 nm; (d) image of two-slit object in the imaging plane in (c).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εd = 2.34 to εd = 4.08, and keep all others unchanged as those in Figure 4(c). The transfer function τA and transmission coefficient τd and τs are plotted in Figure 5(a), from which we can see that the high peak of τs associated with LRSP mode is mainly removed, and τd approximately follows an exponential decay for evanescent waves, and thus the transfer function τA 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).Figure 5
(a) Amplitude of transmission coefficient of superlens and dielectric and transfer function of the system as a function of normalized wavevector. (a) LRSP mode cut-off condition is approached withdd = 10nm, ds = 15nm, and εd = 4.08; (b) Hz field distributions in superlens imaging system; (c) image of two-slit object in the imaging plane in (b).


## 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.



---

*Source: 6513836-2019-05-02.xml*