|
1.IntroductionIn attempts to push the limit of optical lithography, several resolution enhancement techniques (RETs) have been widely applied to improve the lithographic performance.1–3 Inverse lithography technology (ILT),4,5 as an active approach of RETs, is considered as an effective and economically viable way to meet various challenges in current and future technology nodes.6,7 To put ILT into practice, various methods have been proposed by academic as well as industrial communities on mask manufacturing rule constraints,8,9 pattern grouping strategy to accelerate computing,10 GPU-based hardware accelerated techniques,11 and mathematical solution methods for the inverse problem.4,5 In this paper, we focus on improvement of mathematical solution methods, where the computational efficiency is one of the most noteworthy issues. Currently, as critical dimension (CD) shrinks, the pattern density of integrated circuits gets much denser and lithographic process variations, such as lens-wafer defocus and exposure dose variation, become more pronounced. It, therefore, requires finer mask models and through process window compensated approaches, which in return significantly increase the computational burden. Inverse lithography technique treats the mask design as an inverse mathematical problem that aims at synthesizing an input mask to deliver a desired output pattern on the wafer. Liu and Zakhor pioneered a mask design method based on the branch and bound algorithm and simulated annealing.12,13 But this is a time-consuming process. In order to reduce the computational complexity, iterative methods were proposed to solve the inverse problem via an optimization process,14,15 and they were further classified as linear, quadratic, and nonlinear optimization problems.16,17 Meanwhile, Poonawala and Milanfar designed the model-based optical proximity correction system and introduced the steepest descent algorithm for the optimization framework.18,19 Subsequently, the optimization framework was further generalized for phase-shifting masks20,21 and partially coherent imaging systems,22–26 and the optimization algorithm was improved with an active set method21 and with an augmented Lagrangian method.27 Most recently, Lv et al. further improved the computational efficiency of the mask synthesis by using the conjugate gradient and an optimal iterative step,28 and by introducing a mask filtering technique to enhance the regularity of the synthesized mask pattern.29,30 On the other hand, with ever-shrinking feature size, the printed feature becomes increasingly sensitive to the fluctuation of the fabrication process, which limits the yield in semiconductor industry. In order to synthesize masks that are robust to process variations, the average wafer performance is optimized via minimizing the expectation of the difference between the desired pattern and the output pattern with respect to process fluctuations.31–36 This approach takes into account process variations explicitly, and is well understood and easily accomplished. However, it further increases the computational complexity. It is noted that all these methods discretize mask as a raster image constituted by grids (pixels) and perform mask synthesis on a certain fine grid throughout the mask synthesizing process. Since the computational complexity of such iterative methods is ,14,17,18 where is total grid numbers used in the discretization of mask, the mask grid size strongly impacts the computational cost. Naturally, one would like to employ a coarser grid to model mask and, thus, to cut down the runtime. Unfortunately, such straightforward strategy usually leads to inaccurate simulation results and a poorly performing mask that fails to optimally utilize the imaging capability of the lithography scanner. There has been a tradeoff between the runtime and mask quality. In this paper, a cascadic multigrid (CMG) algorithm37,38 is introduced to improve the computational efficiency of robust inverse mask synthesis, while maintaining the model accuracy and mask pattern quality. The fundamental idea of CMG algorithms is to perform more iterations on coarser grids and fewer iterations on fine ones, so that the overall runtime may be dramatically reduced. The CMG algorithm is a one-way multigrid method that requires no coarse grid correction and is, thus, easy to implement. Since the computational complexity of iterative methods is , the CMG algorithm significantly reduces the runtime compared to the conventional methods, which accomplish mask synthesis on a fixed fine grid. More details of the CMG algorithm will be presented in Sec. 3.3. For theoretical convergence results and complexity analysis on CMG algorithms, interested readers may refer to Refs. 3738.39. to 40. Moreover, we have recently developed an analytical circle-sampling technique for fast imaging simulation of partially coherent systems by decomposing the double-impulse response (DIR) function into analytical kernels.41 The DIR function,42 whose Fourier transform is the familiar transmission cross-coefficient (TCC), is projected onto the circle-sampling function (CSF)41 space and converted into a much smaller projected matrix. Singular value decomposition (SVD) is performed to the smaller projected matrix, and then we obtain analytical optical kernels with its eigenvectors and the CSFs. We have demonstrated that this method avoids directly performing SVD to the large TCC matrix; the simulation speed is, thus, dramatically improved. In particular, the derived optical kernels have analytical forms; as a result, the grid size of the kernels can be set to any desired value without re-preforming decomposition. Also, most recently we developed a metric called edge distance error (EDE) to guide mask synthesis.29,30 Compared to the commonly used metric pattern error,16–28 the metric EDE has a dimension of length and is independent of the simulation grid size. The analytical circle-sampling technique and EDE are both independent of grid size. Such nature is well suited for the CMG algorithm, which frequently changes the mask grid size during its synthesizing process. In this paper, we will apply the analytical circle-sampling technique to the forward lithography imaging modeling and employ EDE to guide mask synthesis. The remainder of this paper is organized as follows. Section 2 details the forward lithography imaging model using an analytical circle-sampling technique. Section 3 introduces the robust inverse lithography problem and presents CMG algorithm in detail. Section 4 provides simulation results to demonstrate the validity and efficiency of the proposed method. Finally, we draw some conclusions in Sec. 5. 2.Forward Lithography Imaging ModelWe apply an analytical circle-sampling technique to the forward lithography imaging simulation. The lithography imaging process is often decomposed into two parts, namely the optical image formation and the resist development. The optical image generated by a partially coherent imaging system can be expressed by a bilinear transform in the spatial domain as43 where represents the spatial coordinate (), is the mask transmittance function, denotes complex conjugation, and is called the DIR function whose Fourier transform is the familiar TCC.42 Here, is the mutual intensity function that describes the coherence of the illumination source, and is the point spread function (PSF) of the optical system that may be obtained by an inverse Fourier transform of the pupil function as where is the pupil function without any aberrations and is the normalized pupil plane coordinate.Since the DIR function is Hermitian and band-limited, it is possible to approximate the DIR function as a finite CSF series as41 where is the total number of CSFs, () is the ’th CSF, and are the expansion coefficients. It is noted that the projected matrix , i.e., [], is also Hermitian and positive definite, so that we can perform SVD to further decompose the projected matrix into eigenvalues and eigenvectors, as follows: where is the ’th eigenvector with eigenvectors in total, is its Hermitian conjugate, and is the corresponding eigenvalue. Substituting Eq. (5) into Eq. (4), we can analytically express the DIR function as where () is called the ’th analytical optical kernel, which is calculated as where is the ’th element of . Figure 1 illustrates the flow chart of the DIR function decomposition into analytical kernels.Substituting Eq. (6) into Eq. (1), the optical image is given by Subsequently, the optical image goes through the resist development to form the printed image on the wafer. The resist effect is often approximated by a constant threshold resist model using the following logarithmic Sigmoid function:18 where is the steepness of the Sigmoid function and is the threshold level of the resist.Putting Eqs. (8) and (9) together, the lithography imaging equation can be formulated as where the operator implements the forward mapping from the input mask to the output pattern on the wafer.3.Inverse Mask SynthesisIdeally, the output pattern on the wafer is desired to match the input design intent that serves as the start point in most mask synthesis routines. However, the optical imaging system typically acts as a low-pass spatial frequency filter and cannot deliver the high-frequency components of the intended pattern on wafer.7 The band-limited system causes the output pattern to be a distorted version of the input intent. So, the objective of inverse lithography is to synthesize a mask pattern to precompensate the effect of frequency loss so as to deliver a wafer image closer to the desired pattern. 3.1.Evaluation MetricsIn order to evaluate the performance of the synthesized mask pattern, we introduce some evaluation metrics. An EDE is introduced as a metric to evaluate the difference between the output pattern of the input mask and the desired pattern. Figure 2 depicts the pixel-based representation of a mask pattern and its output pattern on the wafer, where the red dots are discrete sampling elements (grids or pixels) of the patterns, denotes the absolute difference area between the desired pattern contour and the output pattern contour, is the perimeter of the desired pattern contour, and and are the lengths of the element along the and directions, respectively. EDE is defined as Assuming that the element size is small enough in Fig. 2, the absolute difference area can be approximated by multiplying the total number of elements in shadow and the element area as Noticing that the value of the grid in the output pattern is either 0 or 1, the number is approximately equal to the pattern error as where is the norm, is the desired pattern, and is usually called pattern error, which is often employed as a metric to evaluate the difference between the output pattern of the input mask and the desired pattern. However, pattern error is a dimensionless quantity and highly depends on mask feature and simulation grid size. Therefore, pattern error is inappropriate to act as a metric, especially when the simulation grid size is constantly changed. Substituting Eq. (13) into Eq. (12) and then into Eq. (11), we have the expression of EDE asIt is noted that EDE has the dimension of length and, thus, is independent of the simulation grid size. More detailed physical description of EDE can be found in Ref. 30. In addition, in order to quantify the manufacturability of the synthesized mask, mask quadratic error metric and mask complexity metric are adopted.18,19 Since we focus on the binary mask in this paper, the mask quadratic error metric and the mask complexity metric are expressed, respectively, as where , , and is the simulation area, or the total number of grids in mask .3.2.Formulation of Robust Inverse LithographyFrom Sec. 3.1, there are, overall, three evaluation metrics that need to be optimized, i.e., the EDE, the mask quadratic error metric, and the mask complexity metric. In the literature, they are usually combined with certain proportions and to be stated as one single optimization objective as Here, is also called cost function. Therefore, the inverse lithography problem under the nominal condition is defined as However, the formulation Eq. (18) of inverse lithography problem is under the assumption of an ideal imaging system without any process variations. However, in practice, process variations are significant and should not be ignored. In this paper, we take into account the two most important variations in lithography process, i.e., lens-wafer defocus and exposure dose variation, and employ a statistical strategy to optimize the average wafer performance via minimizing the expectation of the difference between the desired pattern and the output pattern with respect to process fluctuations. Specifically, exposure dose variation can be accounted for by varying the resist threshold in Eq. (9); the effect of defocus on the image intensity distribution is effectively an even-type lens aberration, leading to a defocused PSF function .32 where , NA is the numerical aperture, and is the incident light wavelength. Therefore, the expectation of the difference between the desired pattern and the output pattern with respect to process fluctuations is expressed by where are the discretized values of defocus with a density distribution ; are the discretized values of equivalent exposure dose with a density distribution , and and are the corresponding discretization number. is called statistical EDE, and it takes process variations into account explicitly. In this paper, the statistical EDE is used to evaluate the average wafer performance with respect to process fluctuations, namely, to evaluate the robustness of mask.Now, the cost function in the robust inverse lithography is expressed as The robust inverse lithography problem is, thus, defined as 3.3.Cascadic Multigrid AlgorithmIn this section, we develop a CMG algorithm to solve the robust inverse lithography problem. For the setup of CMG algorithm, we define a sequence of spaces with the corresponding space grid sizeSince the grid is usually square, and we use only a parameter to describe the space grid size for simplicity, that means . denotes the finest grid space currently used for discretization of a mask. First, we solve the robust inverse lithography problem, i.e., Eq. (22), on the coarsest space and obtain a mask . The obtained mask on the space is then embedded in the next space via an interpolation operation, and then the interpolated mask servers as the initial value on the space . With the initial value , it calculates a mask on the space . The computation repeats in this manner until a mask on the finest space is obtained. The pseudo-code of CMG algorithm is described in Table 1. Table 1The pseudo-code of cascadic multigrid (CMG) algorithm.
The CMG algorithm starts mask synthesis on the coarsest space . Since the coarsest mask contains much less optimization variables (grids) than the finest mask on space , it spends much less runtime to find a solution and is more likely to change each variable. From the perspective of signal processing, the coarse mask can be considered as the low-frequency portions of mask. The low-frequency portions are first optimized and then further corrected by adding high-frequency details on the refined spaces. By contrast, conventional methods perform mask synthesis directly on the finest space , where it is expensive to change each variable and impossible to separate and treat efficiently the low-frequency information content of the mask. Therefore, the CMG algorithm can achieve a mask pattern by taking less runtime with better solution performance, compared to the conventional methods. This is the fundamental rationale behind the CMG algorithm, and simulation results below demonstrate its advantages. 3.4.Conjugate Gradient MethodIn this section, we present a solution procedure of solving Eq. (22) on a certain space with a conjugate gradient (CG) method. First of all, we calculate the gradient of the cost function with respect to mask as Detailed expression of Eq. (25) is given in the Appendix. In CG method, the optimization direction in the ’th iteration is defined by Here, is a matrix representing the gradient of cost function at mask , and is a factor that depends on different CG methods. When , the evolution velocity is the negative gradient, and the CG method reduces to the steepest descent method. In this paper, we employ the Polak-Ribière-Polyak (PRP) CG method.28 In practical computations, the PRP CG method is generally believed to be one of the most efficient CG methods. It makes use of the velocity information in the previous iteration and can automatically adjust to avoid jamming. It essentially performs a restart when a bad direction occurs. In PRP CG method, is defined as Now, the optimization procedure is described as follows:
In the above procedure, the iteration is terminated when or , where is defined as the minimum value of the norm of velocity and is the prescribed upper limit of the number of iterations. The termination criterion means that the iteration stops when the gradient is zero or rather small. 4.SimulationsSimulations were performed on a partially coherent imaging system with a quasar source illumination (). The wavelength in the simulations was set at 193 nm and NA was 1.35. The resist effect was approximated by a Sigmoid function with and . The mask quadratic error regularization factor and the mask complexity regularization factor were both set at 0.01. The parameter and of the initial mask in Eq. (29) were 0.90 and 0.05, respectively. The value of was set at 0.3 in termination criterion . Instead of computing the step length in Eq. (34) accurately, we directly set at a small constant 0.3. All the simulations were carried out with in-house MATLAB® codes on an HPZ820 Workstation ( GHz Intel Xeons with 8Cores/16Threads) using a Windows 7 (64 bit) operating system. 4.1.Results of Robust Inverse LithographyFigure 3 depicts the synthesized mask patterns by nominal and robust ILT. Figure 3(a) shows the desired pattern consisting of pixels with a grid size of 2.5 nm, which is commonly encountered in logic circuit design. From Fig. 3, the synthesized mask patterns by robust ILT are quite different from that by nominal ILT. It is particularly interesting to note that subresolution assist features (SRAFs) are generated in Fig. 3(c), but not in Fig. 3(b). The reason might be that defocus from the nominal condition results in much worsened image contrast; thus SRAFs tend to be produced in order to enhance the image quality. Moreover, we compare the exposure-defocus (E-D) trees of the synthesized mask patterns by setting the EDE to within of the CD target (45 nm in this case) in Fig. 4. As expected, the synthesized mask patterns by robust ILT have wider E-D windows. This is because the robust ILT further takes the process variations into consideration. As shown in Figs. 3 and 4, the robust ILT tends to generate SRAFs and has the capability of achieving wider E-D windows compared to that synthesized for the nominal condition only. 4.2.Results of Cascadic Multigrid AlgorithmIn this case, we set three level spaces in total, , , and , with the corresponding grid size of 15, 5, and 2.5 nm. A bilinear method is used to interpolate from a coarse grid space to a finer one. Figure 5 depicts the simulated images of the CMG algorithm. Since the grid size on the space is very coarse, the synthesized mask pattern, Fig. 5(b), on the space reveals a general pattern and consists of numerous gray-level SRAFs. These general pattern and gray-level SRAFs are difficult to manufacture in practice and require further correction on the fine spaces. Compared to the mask pattern synthesized by the conventional method [Fig. 3(c)], the mask pattern synthesized by CMG algorithm [Fig. 5(f)] is mostly the same but possesses more SRAFs. Figure 6(a) illustrates the convergence comparison for the statistical EDE. The CMG algorithm takes 132 iterations to find an optimal mask pattern, Fig. 5(b), with a statistical EDE of 4.20 nm, and this mask pattern is then interpolated to serve as an initial value on the space , as shown in Fig. 5(c). It is noted that the convergence history jumps to 5.34 nm here [shown in the blue circle in Fig. 6(a)]; this is because the synthesized coarse mask pattern [Fig. 5(b)] is not exactly accurate and, thus, required further correction. The CMG algorithm takes further 87 iterations on space and reaches a mask pattern, Fig. 5(d), with a statistical EDE of 4.44 nm. Finally, the CMG algorithm achieves a mask pattern, Fig. 5(f), after 20 iterations on the finest space . From Fig. 6(a), it is observed that the CMG algorithm achieves a smaller statistical EDE of 4.46 nm compared to 4.51 nm by the conventional method, which directly synthesizes mask on the finest space . Moreover, we also depict the convergence of pattern error in the CMG algorithm. It is observed that pattern error changed significantly when mask is interpolated from one grid space to another, and it is, thus, not intuitive to study on the convergence. Since pattern error is a dimensionless quantity and highly depends on mask feature and simulation parameters (such as simulation grid resolution), it is not suited for the CMG algorithm. The runtime of the conventional method and the CMG algorithm on each space is shown in Table 2. Since the optimization variables (grids) are relatively less on the space , it synthesizes an optimal mask pattern, Fig. 5(b), after 8.67 s compared to 620.03 s by directly synthesizing mask on the finest space . Overall, the CMG algorithm takes 239 iterations using 135.35 s compared to 193 iterations using 620.03 s by the conventional method; in other words, the CMG algorithm is faster than the conventional method although taking more iterations. Moreover, we compare the E-D windows for the synthesized mask pattern by the conventional method and the CMG algorithm in Fig. 6(b). It is observed that the synthesized mask pattern by the CMG algorithm has slightly improved the E-D window. The reason might be that CMG algorithm generates more SRAFs. Table 2Runtime of the conventional method and CMG algorithm on each space for a desired pattern Fig. 5(a).
Since there are rather less optimization variables grids on the space than grids on the space , it takes much less runtime and is more likely to change each variable to find a better solution. As expected, the synthesized result on the coarsest space contains a general pattern with numerous SRAFs. From the perspective of signal processing, these patterns can be considered low-frequency portions in mask, and these patterns are further corrected by adding high-frequency details on the fine space. The conventional method accomplishes mask synthesis directly on the fine space ; it, therefore, is expensive to change each variable and cannot handle the low-frequency portions in mask explicitly. As a result, the CMG algorithm has the capacity of achieving a mask pattern with smaller statistical EDE and wider E-D window compared to the conventional method. We also performed simulations for a more complicated pattern in Fig. 7. Figure 7(a) shows the desired pattern, which is a metal layer of the benchmark clock gate circuit pattern within a simulation area of .44 As expected, the synthesized mask pattern, Fig. 7(b), by nominal ILT is quite different from that synthesized by robust ILT, Figs. 7(c) and 7(d), and the robust ILT has significantly improved the E-D window (see Fig. 8). On the other hand, compared to the synthesized mask pattern by the conventional method, CMG algorithm reaches a mostly same pattern but with more SRAFs; the E-D window has, thus, slightly improved. The runtime of the conventional method and the CMG algorithm on each space is shown in Table 3. It is revealed that the CMG algorithm is overall faster than the conventional method. This set of simulations further demonstrates that the CMG algorithm has big improvement in computation efficiency, along with a slightly better lithographic performance (wider E-D window). Table 3Runtime of the conventional method and CMG algorithm on each space for a desired pattern Fig. 7(a).
5.ConclusionsIn this paper, we report a CMG method for robust inverse mask synthesis. This method synthesizes mask hierarchically. It first computes an initial mask pattern on a coarse grid and, therefore, spends less runtime and is more likely to change each variable to find a better solution. This initial mask is then interpolated to a fine grid and to be further corrected, which requires a small number of iterations on the fine grid. Overall, the CMG algorithm is able to run more than four times faster than conventional mask synthesis methods working on a fixed fine grid. In addition, we apply an analytical circle-sampling technique to the forward lithography imaging modeling and employ a novel EDE metric to guide mask synthesis. These two techniques are both independent of grid size and are well suited for the CMG algorithm. It is expected that the proposed method will provide a useful mask synthesis strategy in practical optical lithography. AcknowledgmentsThis work was funded by the National Natural Science Foundation of China (Grant No. 91023032, 51005091, 51121002), the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20120142110019), the National Science and Technology Major Project of China (Grant No. 2012ZX02701001), and the National Instrument Development Specific Project of China (Grant No. 2011YQ160002. It was also supported in part by the UGC Areas of Excellence project Theory, Modeling, and Simulation of Emerging Electronics, and by the State Key Lab of Digital Manufacturing Equipment and Technology under Project DMETKF2013003. AppendicesAppendix:Derivation of Eq. (25)The gradient of the cost function with respect to mask is expressed in Eq. (25). In the following, we present the detailed expression of each term in Eq. (25), i.e., , , and . In the present content, and denote the spatial coordinate (). According to Ref. 28, with a fixed defocus and dose variation , the gradient of EDE with respect to mask is where () is the ’th optical kernels with the defocus and is the corresponding resist image of input mask under the defocus and dose variation . So the gradient of the statistical EDE with respect to mask isThe gradient of the mask quadratic error with respect to mask is The gradient of the mask complexity with respect to mask is where sign[·] is the Signum function, and is an operator of the first derivative asReferencesA. K. Wong, Resolution Enhancement Techniques in Optical Lithography, SPIE Press, Bellingham, Washington
(2001). Google Scholar
L. W. Liebmannet al.,
“TCAD development for lithography resolution enhancement,”
IBM J. Res. Dev., 45
(5), 651
–665
(2001). http://dx.doi.org/10.1147/rd.455.0651 IBMJAE 0018-8646 Google Scholar
F. Schellenberg,
“Resolution enhancement technology: the past, the present, and extensions for the future,”
Proc. SPIE, 5377 1
–20
(2004). http://dx.doi.org/10.1117/12.548923 PSISDG 0277-786X Google Scholar
L. PangY. LiuD. Abrams,
“Inverse lithography technology (ILT): what is the impact to the photomask industry?,”
Proc. SPIE, 6283 62830X
(2006). http://dx.doi.org/10.1117/12.681857 PSISDG 0277-786X Google Scholar
L. PangY. LiuD. Abrams,
“Inverse lithography technology (ILT): a natural solution for model-based SRAF at 45 nm and 32 nm,”
Proc. SPIE, 6607 660739
(2007). http://dx.doi.org/10.1117/12.729028 PSISDG 0277-786X Google Scholar
D. O. S. Melvilleet al.,
“Computational lithography: exhausting the resolution limits of 193-nm projection lithography systems,”
J. Vac. Sci. Technol. B, 29
(6), 06FH04
(2011). http://dx.doi.org/10.1116/1.3662090 JVTBD9 0734-211X Google Scholar
M. L. Rieger,
“Communication theory in optical lithography,”
J. Micro/Nanolith. MEMS MOEMS, 11
(1), 013003
(2012). http://dx.doi.org/10.1117/1.JMM.11.1.013003 JMMMGF 1932-5134 Google Scholar
A. Villaretet al.,
“Inverse lithography technique for advanced CMOS nodes,”
Proc. SPIE, 8683 86830E
(2013). http://dx.doi.org/10.1117/12.2010111 PSISDG 0277-786X Google Scholar
T. Cecilet al., , “Enhancing fullchip ILT mask synthesis capability for IC manufacturability,” Proc. SPIE 7973, 79731C (2011). Google Scholar S. S. Yehet al.,
“Study of the pattern aware OPC,”
Proc. SPIE, 9052 90521Y
(2014). http://dx.doi.org/10.1117/12.2046059 Google Scholar
I. Torunogluet al.,
“A GPU-based full-chip inverse lithography solution for random patterns,”
Proc. SPIE, 7641 764115
(2010). http://dx.doi.org/10.1117/12.846638 Google Scholar
Y. LiuA. Zakhor,
“Optimal binary image design for optical lithography,”
Proc. SPIE, 1264 401
–412
(1990). http://dx.doi.org/10.1117/12.20216 PSISDG 0277-786X Google Scholar
Y. LiuA. Zakhor,
“Binary and phase shifting mask design for optical lithography,”
IEEE Trans. Semiconduct. Manuf., 5
(2), 138
–152
(1992). http://dx.doi.org/10.1109/66.136275 ITSMED 0894-6507 Google Scholar
Y. C. PatiT. Kailath,
“Phase-shifting masks for microlithography: automated design and mask requirements,”
J. Opt. Soc. Am. A, 11
(9), 2438
–2452
(1994). http://dx.doi.org/10.1364/JOSAA.11.002438 JOAOD6 0740-3232 Google Scholar
S. SherifB. SalehR. D. Leone,
“Binary image synthesis using mixed linear integer programming,”
IEEE Trans. Image Process., 4
(9), 1252
–1257
(1995). http://dx.doi.org/10.1109/83.413169 IIPRE4 1057-7149 Google Scholar
Y. Granik,
“Solving inverse problems of optical microlithography,”
Proc. SPIE, 5754 506
–526
(2005). http://dx.doi.org/10.1117/12.600141 PSISDG 0277-786X Google Scholar
Y. Granik,
“Fast pixel-based mask optimization for inverse lithography,”
J. Micro/Nanolith. MEMS MOEMS, 5
(4), 043002
(2006). http://dx.doi.org/10.1117/1.2399537 JMMMGF 1932-5134 Google Scholar
A. PoonawalaP. Milanfar,
“Mask design for optical microlithography—an inverse imaging problem,”
IEEE Trans. Image Process., 16
(3), 774
–788
(2007). http://dx.doi.org/10.1109/TIP.2006.891332 IIPRE4 1057-7149 Google Scholar
A. PoonawalaP. Milanfar,
“A pixel-based regularization approach to inverse lithography,”
Microelectron. Eng., 84
(12), 2837
–2852
(2007). http://dx.doi.org/10.1016/j.mee.2007.02.005 MIENEF 0167-9317 Google Scholar
X. MaG. R. Arce,
“Generalized inverse lithography methods for phase shifting mask design,”
Opt. Express, 15
(23), 15066
–15079
(2007). http://dx.doi.org/10.1364/OE.15.015066 OPEXFF 1094-4087 Google Scholar
S. H. ChanA. K. WongE. Y. Lam,
“Initialization for robust inverse synthesis of phase-shifting masks in optical projection lithography,”
Opt. Express, 16
(19), 14746
–14760
(2008). http://dx.doi.org/10.1364/OE.16.014746 OPEXFF 1094-4087 Google Scholar
X. MaG. R. Arce,
“Binary mask optimization for inverse lithography with partially coherent illumination,”
J. Opt. Soc. Am. A, 25
(12), 2960
–2970
(2008). http://dx.doi.org/10.1364/JOSAA.25.002960 JOAOD6 0740-3232 Google Scholar
E. Y. LamA. K. Wong,
“Computation lithography: virtual reality and virtual virtuality,”
Opt. Express, 17
(15), 12259
–12268
(2009). http://dx.doi.org/10.1364/OE.17.012259 OPEXFF 1094-4087 Google Scholar
J. C. YuP. Yu,
“Impacts of cost functions on inverse lithography patterning,”
Opt. Express, 18
(22), 23331
–23342
(2010). http://dx.doi.org/10.1364/OE.18.023331 OPEXFF 1094-4087 Google Scholar
X. MaG. R. Arce,
“Pixel-based OPC optimization based on conjugate gradients,”
Opt. Express, 19
(3), 2165
–2180
(2011). http://dx.doi.org/10.1364/OE.19.002165 OPEXFF 1094-4087 Google Scholar
J. LiY. ShenE. Y. Lam,
“Hotspot-aware fast source and mask optimization,”
Opt. Express, 20
(19), 21792
–21804
(2012). http://dx.doi.org/10.1364/OE.20.021792 OPEXFF 1094-4087 Google Scholar
J. LiS. Y. LiuE. Y. Lam,
“Efficient source and mask optimization with augmented Lagrangian methods in optical lithography,”
Opt. Express, 21
(7), 8076
–8090
(2013). http://dx.doi.org/10.1364/OE.21.008076 OPEXFF 1094-4087 Google Scholar
W. Lvet al.,
“Level-set-based inverse lithography for mask synthesis using the conjugate gradient and an optimal time step,”
J. Vac. Sci. Technol. B, 31
(4), 041605
(2013). http://dx.doi.org/10.1116/1.4813781 JVTBD9 0734-211X Google Scholar
W. LvQ. XiaS. Y. Liu,
“Pixel-based inverse lithography using a mask filtering technique,”
Proc. SPIE, 8683 868325
(2013). http://dx.doi.org/10.1117/12.2011469 PSISDG 0277-786X Google Scholar
W. LvQ. XiaS. Y. Liu,
“Mask-filtering-based inverse lithography,”
J. Micro/Nanolith. MEMS MOEMS, 12
(4), 043003
(2013). http://dx.doi.org/10.1117/1.JMM.12.4.043003 JMMMGF 1932-5134 Google Scholar
N. JiaE. Y. Lam,
“Machine learning for inverse lithography: using stochastic gradient descent for robust photomask synthesis,”
J. Opt., 12
(4), 045601
(2010). http://dx.doi.org/10.1088/2040-8978/12/4/045601 JOOPDB 0150-536X Google Scholar
Y. Shenet al.,
“Robust level-set-based inverse lithography,”
Opt. Express, 19
(6), 5511
–5521
(2011). http://dx.doi.org/10.1364/OE.19.005511 OPEXFF 1094-4087 Google Scholar
N. JiaE. Y. Lam,
“Pixelated source mask optimization for process robustness in optical lithography,”
Opt. Express, 19
(20), 19384
–19398
(2011). http://dx.doi.org/10.1364/OE.19.019384 OPEXFF 1094-4087 Google Scholar
S. K. Choyet al.,
“A robust computational algorithm for inverse photomask synthesis in optical projection lithography,”
SIAM J. Imaging Sci., 5
(2), 625
–651
(2012). http://dx.doi.org/10.1137/110830356 1936-4954 Google Scholar
S. LiX. WangY. Bu,
“Robust pixel-based source and mask optimization for inverse lithography,”
Opt. Laser Technol., 45 285
–293
(2013). http://dx.doi.org/10.1016/j.optlastec.2012.06.033 OLTCAS 0030-3992 Google Scholar
S. Y. Liuet al.,
“Convolution-variation separation method for efficient modeling of optical lithography,”
Opt. Lett., 38
(13), 2168
–2170
(2013). http://dx.doi.org/10.1364/OL.38.002168 OPLEDP 0146-9592 Google Scholar
F. A. BornemannP. Deuflhard,
“The cascadic multigrid method for elliptic problems,”
Numer. Math., 75
(2), 135
–152
(1996). http://dx.doi.org/10.1007/s002110050234 NUMMA7 0029-599X Google Scholar
D. BraessP. DeuflhardK. Lipnikov,
“A subspace cascadic multigrid method for mortar elements,”
Computing, 69
(3), 205
–225
(2002). http://dx.doi.org/10.1007/s00607-002-1460-2 0010-485X Google Scholar
W. L. BriggsV. E. HensonS. F. McCormick, A Multigrid Tutorial, 2nd ed.SIAM Press, Philadelphia, PA
(2000). Google Scholar
L. ReichelA. Shyshkov,
“Cascadic multilevel methods for ill-posed problems,”
J. Comput. Appl. Math., 233
(5), 1314
–1325
(2010). http://dx.doi.org/10.1016/j.cam.2009.03.019 JCAMDI 0377-0427 Google Scholar
P. Gonget al.,
“Fast aerial image simulations for partially coherent systems by transmission cross coefficient decomposition with analytical kernels,”
J. Vac. Sci. Technol. B, 30
(6), 06FG03
(2012). http://dx.doi.org/10.1116/1.4767442 JVTBD9 0734-211X Google Scholar
B. E. A. Saleh,
“Optical bilinear transformations: general properties,”
Opt. Acta, 26
(6), 777
–799
(1979). http://dx.doi.org/10.1080/713820055 OPACAT 0030-3909 Google Scholar
A. K. Wong, Optical Imaging in Projection Microlithography, SPIE Press, Bellingham, Washington
(2005). Google Scholar
“The NanGate 45nm Open Cell Library, v1.3,”
(2013) http://www.si2.org/openeda.si2.org/projects/nangatelib July ). 2013). Google Scholar
BiographyWen Lv is currently a PhD candidate at Huazhong University of Science and Technology under the guidance of Professor Shiyuan Liu. He received his BS degree from the School of Mechanical Science and Engineering of the same university in 2011. His research involves various issues in optical lithography, including inverse lithography, fast optical image simulation, and mask writing technique. He is a student member of SPIE and IEEE. Edmund Y. Lam received his BS (with distinction), MS, and PhD degrees in electrical engineering from Stanford University. He was a senior engineer at KLA-Tencor Corporation in San Jose, California, before joining the Electrical and Electronic Engineering Department at the University of Hong Kong, where he is currently an associate professor. He is also the founding director of its Imaging Systems Laboratory. His research interests include computational optics and imaging, particularly their applications in the semiconductor manufacturing process. He is a fellow of SPIE and OSA. Haiqing Wei received his BS in physics from Peking University, China, MS in physics, and PhD in electrical and computer engineering from McGill University, Canada. He has been a research scientist and engineer for over a decade, working in various subjects of fiber-optic communications, optical imaging, nanophotonics, and computational lithography. He has over 20 patent applications in these subject areas and has published over 30 scientific articles in peer-reviewed journals. Since April 2012, he has been appointed as a guest professor at Huazhong University of Science and Technology. Shiyuan Liu is a professor of mechanical engineering at Huazhong University of Science and Technology, leading his Nanoscale and Optical Metrology Group with research interest in metrology and instrumentation for nanomanufacturing. He also actively works in the area of optical lithography, including partially coherent imaging theory, wavefront aberration metrology, optical proximity correction, source mask optimization, and inverse lithography technology. He received his PhD in mechanical engineering from Huazhong University of Science and Technology in 1998. He is a member of SPIE, OSA, AVS, IEEE, and Chinese Society of Micro/Nano Technology. He holds 21 patents and has authored or coauthored more than 100 technical papers. |