The Inverse Optimization of an Optical Lithographic Source with a Hybrid Genetic Algorithm

: As an effective resolution enhancement technology, source optimization (SO) is considered key for signiﬁcantly improving the image quality of optical lithography at advanced nodes. To solve the problem of unsatisfactory SO performance, it is necessary to combine it with optimization algorithms. In this study, an SO method based on a hybrid genetic algorithm is proposed to achieve an acceptable source shape in the imaging process for optical lithography. To overcome the problems of local optima and the small search scope, an update strategy that uses particle swarm optimization and the tabu list method from the tabu search algorithm are utilized to enhance the optimization performance. Meanwhile, different feature patterns were employed as the input of the optimization model. These simulation results show that the proposed SO method exhibits dominant optimization performance for SO in optical lithography.


Introduction
Optical lithography is a critical part of the manufacturing of very-large-scale integrated circuits (VLSIs), and it can be used to implement pattern transfer from a mask to a wafer with a reduction in the magnification by using the optical projection method [1][2][3]. The imaging process in optical lithography is shown in Figure 1. The rays emitted by an illumination source pass through a group of optical mirrors to form Kohler illumination. These rays irradiate the mask to generate diffracted light, which carries the feature information. An objective lens collects the diffracted light, which finally creates an image on the surface of a wafer. Because of the low-pass filtering effect of the objective lens, the high-frequency information of the feature pattern is lost, which causes aerial image to be distorted. Therefore, improving the imaging quality in optical lithography is essential in current advanced technology nodes.
According to the Rayleigh criterion R = k 1 λ/N A, the lithographic resolution (R) can be absolutely improved to shrink the illumination wavelength (λ) and increase the numerical aperture (N A) of the objective lens, where k 1 is the process factor [4,5]. In advanced nodes, the loss of high-frequency information is becoming more and more serious, and this causes an obvious optical proximity effect (OPE) [6,7]. These two approaches can no longer guarantee that the imaging quality requirements of continuously evolving technology will be met. Several resolution enhancement technologies (RETs) for the attainment of high-quality aerial image have been proposed [8][9][10]. In these methods, pixelated source optimization (SO), which is an inverse optimization method, can efficiently promote the performance in optical lithographic imaging [11][12][13]. This can be achieved with mathematical and algorithmic approaches, which can inversely optimize the intensity distribution of an optical lithographic source according to aerial image or resistance pattern. mathematical and algorithmic approaches, which can inversely optimize the intensity distribution of an optical lithographic source according to aerial image or resistance pattern. The enhancement of optimization performance with pixelated SO methods has become a research hotspot. Freeform illuminations, which can be created with a set of micromirror arrays, have been developed to provide more degrees of freedom for regulating the illumination shape in optical lithography [11,14,15]. In a pixelated source model, changes in intensity of the pixels can indicate the mirror flip angle, where a pixel represents a unit mirror.
To improve the imaging quality in optical lithography, various types of pixelated SO methods have been proposed; of these, gradient-based methods, such as gradient descent (GD) [16], conjugate gradient (CG) [17,18], and steepest descent (SD) [19], have become the most significant. As local optimization algorithms, gradient-based methods have quick convergence, robustness, and the ability to search for the optimal solution [16,20,21]. These methods are not limited by the size of the optimized matrix, but only depend on the complexity of the gradient function. However, with the decrease in the critical dimensions of feature patterns toward smaller process nodes, the complexity of photo-resistance models is continuously expanded. It is significantly complicated to derive the gradient function of a photo-resistance model. Furthermore, the optimal solutions of complex models have an increased probability of falling into local optima when applying these methods to process the optimization of lithographic sources in advanced applications. To improving the imaging performance in optical lithography, it is not enough to address the optimization efficiency.
The genetic algorithm (GA), which is a typical example of a heuristic approach, has been widely employed as a key tool for achieving the optimization of complex models in Figure 1. Schematic representation of the optical lithographic imaging process r s represents the effective radius of the illumination source, r p is the radius of the pupil, and θ max is the maximal incident angle of the rays in the objective exit.
Freeform illuminations, which can be created with a set of micromirror arrays, have been developed to provide more degrees of freedom for regulating the illumination shape in optical lithography [11,14,15]. In a pixelated source model, changes in intensity of the pixels can indicate the mirror flip angle, where a pixel represents a unit mirror.
To improve the imaging quality in optical lithography, various types of pixelated SO methods have been proposed; of these, gradient-based methods, such as gradient descent (GD) [16], conjugate gradient (CG) [17,18], and steepest descent (SD) [19], have become the most significant. As local optimization algorithms, gradient-based methods have quick convergence, robustness, and the ability to search for the optimal solution [16,20,21]. These methods are not limited by the size of the optimized matrix, but only depend on the complexity of the gradient function. However, with the decrease in the critical dimensions of feature patterns toward smaller process nodes, the complexity of photoresistance models is continuously expanded. It is significantly complicated to derive the gradient function of a photo-resistance model. Furthermore, the optimal solutions of complex models have an increased probability of falling into local optima when applying these methods to process the optimization of lithographic sources in advanced applications. To improving the imaging performance in optical lithography, it is not enough to address the optimization efficiency.
The genetic algorithm (GA), which is a typical example of a heuristic approach, has been widely employed as a key tool for achieving the optimization of complex models in several fields, such as production scheduling, image processing, and data mining [22,23]. N. Sortrakul et al. utilized a GA to solve a comprehensive optimization model for production scheduling and preventive maintenance planning [24]. Zhang et al. used a new adaptive operator to achieve a nonlinear adaptive adjustment for image segmentation via a GA [25]. Tzung et al. found a suitable membership function for a mining problem through a GAbased framework and applied the best-set membership function in fuzzy association rule mining [26]. As this is a global optimization method, it is more convenient to find the optimal condition of the merit function for a complex model because of its gradient-free features. Moreover, due to its inherent property of being based on a probabilistic rule, the more flexible search function of a GA can be implemented to decrease the influence of multiple parameters of a complex model to improve the optimization efficiency. Several studies have indicated that GAs have displayed considerable capacities for optimizing the imaging performance in optical lithography. Fuhner et al. proposed that it is significantly effective to improve the imaging quality by using a GA to optimize the illumination geometries in optical lithography [27]. Li et al. applied a GA for optical proximity correction to compensate for the image distortion induced by the optical proximity effect [28]. Hong et al. used a modified GA method to improve the process windows in inverse lithography for various common objectives, including the depth of focus, normalized image log slope, and process variability band [29]. Nevertheless, the multiple parameters of a complex model potentially cause the search range of an optimization model based on a GA to be reduced, thus causing premature convergence [30,31]. Due to the scalability of GAs, it is possible to combine a GA with other algorithms to enhance the optimization capabilities and improve the convergence efficiency [24,32,33]. To achieve the optimal source shape in the imaging process of optical lithography for advanced applications, a hybrid GA can be considered as an effective approach for complex imaging models.
In this study, a hybrid algorithm created by combining a GA with two decision methods is proposed to find the optimal illumination mode according to different feature patterns in order to improve the imaging performance in optical lithography. These decision methods are the iterative search approaches of the tabu search algorithm (TS) and particle swarm optimization algorithm (PSO). Because of the complex model of the optical lithographic imaging process, the convergence efficiency of GA-based SO methods is limited and can be excessively reduced. To overcome the search limitations, an update strategy using PSO can expand the search scope. The tabu list function from the TS can improve the optimization capabilities within a reasonable range. Furthermore, the combination of a crossover operation and mutation operation in the GA can ensure the diversification of the population for the SO model. The simulation results prove that the imaging performance in optical lithography can be significantly improved by utilizing an SO method with the proposed hybrid GA approach.

Partially Coherent Imaging Model
The forward imaging process in optical lithography can be described by a partially coherent imaging model (PCI), which is equivalent to the superposition of a series of coherent images. Assuming that a point source on the optical axis is utilized to complete completely coherent imaging (CCI), this imaging process can be expressed as follows [34]: where (x, y) are the spatial coordinates in the imaging plane, ( f , g) are the spectral coordinates of the pupil, and ( f , g ) are the spectral coordinates of the mask spectrum. I (x, y) represents the intensity distribution of an image obtained via CCI. H( f , g), which is a pupil function, can be regarded as a low-pass filter, which severely restricts the higher-order diffraction information of the feature pattern from participating in the imaging process. In addition, it can be used to simulate the effects of defocus, aberration, and apodization by modifying the phase information of the pupil. M( f , g ) denotes the frequency spectrum of the feature pattern in the mask, which can be acquired with a Fourier transform operation.
To facilitate the imaging process of optical lithography, the illumination source can be depicted as discrete matrices with a pixelated representation, where a pixel represents a point source. For the off-axis point source, in the process of PCI, the imaging process can be implemented by shifting the pupil according to the position of this point source. The intensity distribution of an image obtained via PCI can be found by accumulating the images obtained by each point source. Hence, the intensity distribution of an aerial image for optical lithography can be calculated as follows: In Equation (2), I(x,y) represents the intensity distribution of an aerial image for optical lithography. J(f,g) is a discrete illumination source. In order to concisely express the imaging process of PCI, Equation (2) can be transformed into a format for matrix multiplication-a source matrix and mask matrix-and the result of CCI can be described by using the illumination cross-coefficient (ICC).
To evaluate the imaging performance in optical lithography, it is vital to describe the resistance patterns on the surface of a wafer. In this study, the formation process of a resistance pattern can be approximately explained by a sigmoid model, which is given by: In Equation (4), R represents the contour layout of the resistance pattern formed due to the photo-resistance effect. sig(·) is a sigmoid function that can approximately achieve a binary distribution. The parameters α and t r denote the steepness index and the threshold of photo-resistance, respectively.

Pattern Error Minimization Problem
In a complicated optimization model, it is essential to construct an efficient merit function that can distinctly describe the convergence tendency of an iteration process. For the proposed hybrid optimization approach based on a heuristic algorithm, the optimization process has no complex gradient derivation procedures for the merit function. In this study, the pattern error (PE) is regarded as the merit function, and it is defined as the absolute sum of the unit differences between the mask and resistance patterns, where f PE represents the value of the PE for each iteration in the optimization model, and M is the discrete matrix formation of the mask. It is assumed that the mask is a 2D matrix with a size of N × N, which has the form of the binary distribution. The value of each pixel is set to 0 or 1. Let R and M be represented by the real matrices, which are denoted by R ∈ R N×N and M ∈ R N×N , respectively. Equation (5) shows that the minimum value of the PE can be achieved by updating the intensity distribution of the illumination source to cause the resistance pattern to be closer to the mask layout throughout the iterative process.

The Optimization Strategy for the Illumination Source
In imaging systems for optical lithography, an extended source is employed as the illumination source, and is deemed to be partially coherent illumination. According to their different partial coherence factors and different intensity distributions, the illumination modes can be divided into circular illumination, annular illumination, dipolar illumination, Appl. Sci. 2023, 13, 5708 5 of 12 quadrupole illumination, etc. In this study, the annular illumination mode is utilized as the initial source in the proposed SO model, as shown in Figure 2. The partially coherent factor (σ) of this annular source is determined by the source's effective radius (r s ) and the pupil's radius (r p ), σ = r s /r p . For annular illumination in the optimization model, the maximum and minimum partially coherent factors are, respectively, σ out = r s,max /r p and σ in = r s,min /r p , where r s,max and r s,min are the maximum and minimum radiuses of the illumination source, respectively. out the iterative process.

The Optimization Strategy for the Illumination Source
In imaging systems for optical lithography, an extended source is employed as the illumination source, and is deemed to be partially coherent illumination. According to their different partial coherence factors and different intensity distributions, the illumination modes can be divided into circular illumination, annular illumination, dipolar illumination, quadrupole illumination, etc. In this study, the annular illumination mode is utilized as the initial source in the proposed SO model, as shown in Figure 2. The partially coherent factor ( σ ) of this annular source is determined by the source's effective radius (  Because of a special property of symmetry with respect to the optical axis, an optical lithographic source can be decomposed into four identical parts according to the axis quadrants, which are J 1 , J 2 , J 3 , and J 4 in Figure 2a. It is essential to solely optimize one of the four regions of an annular source to reduce the computational complexity of the optimization model. The region J 1 shown in Figure 2b is hypothetically employed as the part to be optimized. Assuming that the size of the pixelated source is N J × N J , the coordinate range of the X-Y axis in the first quadrant is 0 to (N J − 1)/2. An indexing operation that is implemented with Equation (6) can be utilized to record the positions of valid pixels in J 1 . P = Ψ{J 1 (x, y) == 1}, x = 0, 1, 2, · · · , N J − 1)/2, y = 0, 1, 2, · · · , N J − 1)/2.
In Equation (6), P is an index result that contains the position information of the valid pixel points in J 1 . Assuming that the number of valid pixel points is hypothetically m, P can be represented by a real matrix, which is denoted by P ∈ R m×2 . Ψ{·} is a query function that can achieve the indexing operation with the criterion of J 1 (x, y) == 1. Then, the 1-D vector matrix J 1 shown in Figure 2c is generated as the objective to be optimized, and it is explained by J 1 ∈ R 1×m . In the proposed hybrid algorithm, a GA is employed as the main body of the optimization framework to update the intensity value of a pixel point in J 1 . The initial population is necessarily created as the initial solution group, and it can be defined as a real matrix denoted by J ∈ R n×m . After several iterations, the initial population matrix can be continuously updated to select the optimal condition J 1 according to the value of the merit function f PE , where J 1 ∈ R 1×m . Then, the intensity distribution of one quarter of the source can be attained by realigning the matrix J 1 via the recorded position of J 1 . The integral source shape, which is shown in Figure 2e, can be generated by using a mirror operation to flip the quarter of the source displayed in Figure 2d. The procedure of optimizing the source is described in Figure 3. query function that can achieve the indexing operation with the criterion of 1 ( , ) 1 J x y == . Then, the 1-D vector matrix 1 J shown in Figure 2c is generated as the objective to be optimized, and it is explained by . In the proposed hybrid algorithm, a GA is employed as the main body of the optimization framework to update the intensity value of a pixel point in 1 J . The initial population is necessarily created as the initial solution group, and it can be defined as a real matrix denoted by '' n m × ∈ J . After several iterations, the initial population matrix can be continuously updated to select the optimal condition . Then, the intensity distribution of one quarter of the source can be attained by realigning the matrix 1 '' J via the recorded position of 1 J . The integral source shape, which is shown in Figure   2e, can be generated by using a mirror operation to flip the quarter of the source displayed in Figure 2d. The procedure of optimizing the source is described in Figure 3. Because of the complexity of imaging models for optical lithography, it is difficult to utilize a traditional GA to attain the optimal source. Although GAs have a great convergence effect in iterative processes for simple optimization models, premature convergence can cause the merit function to fall into a local optimum in a complicated model. To improve the performance of the optimization model based on a GA, optimization strategies-the method of updating the population in PSO and the search method in TS-were inserted into the GA to overcome the local optimal conditions. In the proposed hybrid algorithm, the GA executes mutation and crossover operations to revise the current population in each iteration. The updated population can be calculated with the following formulas in the mutation process. Because of the complexity of imaging models for optical lithography, it is difficult to utilize a traditional GA to attain the optimal source. Although GAs have a great convergence effect in iterative processes for simple optimization models, premature convergence can cause the merit function to fall into a local optimum in a complicated model. To improve the performance of the optimization model based on a GA, optimization strategies-the method of updating the population in PSO and the search method in TS-were inserted into the GA to overcome the local optimal conditions. In the proposed hybrid algorithm, the GA executes mutation and crossover operations to revise the current population in each iteration. The updated population can be calculated with the following formulas in the mutation process.
In Equation (7), represents a random increment that can be generated with a random function rand() in different iteration loops. i and τ are the current iteration number and the total number of iterations, respectively. The increase and decrease in individuals with respect to the magnification in P can be, respectively, implemented with Equations (8) and (9). P(j, k) represents the units of individuals in the jth row and kth column, where j = 1, 2, 3, · · · , n and k = 1, 2, 3, · · · , m. In addition, it is necessary for the crossover operation of the GA to be employed to enlarge the search range in order to enhance the convergence abilities. The node to be crossed in P can be found with where Ψ represents a mathematical function that implements a rounding operation to generate the nearest integer. The exchanged part of the individual can be decided according to the position (η obj , η node ). Hence, the population P can be updated with the mutation and crossover operations in the procedure of the GA. The optimal individual in the current iteration can be selected from P, and it is recorded on the tabu list. It is significant to establish the list in order to enhance the search capabilities and break away from local optima. The tabu list can store the number of individuals, and it is calculated with the following formula.
To ensure the global conditions in the complex SO model, PSO is employed to enlarge the update scope of the variates. The strategy of using a weight coefficient based on the hyperbolic tangent is proposed to manipulate the speed of updating the variates, which is also varied in a different loop, and this is denoted by the following: In Equation (13), W (i) represents the weight coefficient in the ith iteration, T [·] is the hyperbolic tangent function, and ξ is the constant. W max and W min are, respectively, the maximum and minimum weight coefficients, which are set before the iteration starts. Equations (14) and (15) indicate the strategy of updating the individual J l in PSO. V l (i + 1) represents the search speed of the variate for the lth individual in the ith iteration. Q current,l and Q global are the current optimal value and the global optimal value, respectively. C 1 and C 2 are both learning factors and are constants.

Simulations
In this section, several simulation results involving aerial images with optimized sources are provided to verify the superior performance in optical lithography after utilizing the pixelated SO method based on the proposed hybrid GA approach. An imaging model for optical lithography is established, and this was based on a 193 nm immersion lithography system in the 45 nm node. In this model, the numerical aperture of the projection objective was 1.25. The annular illumination mode was employed for the imaging source in the initial stage to provide the initial aerial image, where σ out = 0.98 and σ in = 0.65. Furthermore, this pixelated annular source was set as a 41 × 41 matrix. To generate the resistance pattern, a sigmoid function was utilized in the optimization model, where α = 85 and t r = 0.19. Figure 4 illustrates the initial source and mask employed to accomplish the imaging process for optical lithography. For Figure 4a-c, the value of the pixels in the white region was set to 1; otherwise, their value was 0. To ensure the continuous variation of the valid pixels and to reduce the dispersion degree of the annular source in the optimization process, the blurred source shown in Figure 4b could be generated by performing a fuzzy operation on the binarized source in Figure 4a. To verify the optimization performance of the pixelated SO based on the proposed method, two patterns were employed as input feature patterns to form an aerial image for optical lithography; these were the horizontal block pattern (#M1) and the periodic optical grating pattern (#M2). They could be represented as real matrices with the size N × N, where N = 521. pixels and to reduce the dispersion degree of the annular source in the optimization process, the blurred source shown in Figure 4b could be generated by performing a fuzzy operation on the binarized source in Figure 4a. To verify the optimization performance of the pixelated SO based on the proposed method, two patterns were employed as input feature patterns to form an aerial image for optical lithography; these were the horizontal block pattern (#M1) and the periodic optical grating pattern (#M2). They could be represented as real matrices with the size N N × , where 521 N = . In these simulations, the proposed pixelated-SO-based method was carried out in three runs to evaluate its optimization performance. In each simulation, there were 150 iterations to optimize the intensity distribution of the lithographic source for the horizontal block pattern (#M1) and the periodic optical grating pattern (#M2). The final lithographic source, which was defined as the optimal source, was used as the input condition to complete the lithographic imaging process and to attain the resistance pattern. These simulation results are shown in Figure 5. Rows 1 and 2 present the source optimization results and the resistance patterns from three runs. The intensity distributions of the optimal sources for (#M1) and (#M2) are, respectively, displayed in the first and third columns of Figure 5. Their resistance patterns are, respectively, shown in the second and fourth columns, and they were obtained by utilizing the optimal source. In these optimized sources, the intensity distributions for #M1 were closer to a vertical quadrupole illumination mode. This is because, in the vertical direction, the pattern had a similar periodic grating structure to that of the vertical bipolar illumination mode. Horizontally, the middle area lacked pattern information for the horizontal illumination mode. For #M2, a 45degree quadrupole illumination mode was more necessary for ensuring better image quality in the imaging process for optical lithography due to the periodic horizontal grating structure of the array. In these simulations, the proposed pixelated-SO-based method was carried out in three runs to evaluate its optimization performance. In each simulation, there were 150 iterations to optimize the intensity distribution of the lithographic source for the horizontal block pattern (#M1) and the periodic optical grating pattern (#M2). The final lithographic source, which was defined as the optimal source, was used as the input condition to complete the lithographic imaging process and to attain the resistance pattern. These simulation results are shown in Figure 5. Rows 1 and 2 present the source optimization results and the resistance patterns from three runs. The intensity distributions of the optimal sources for (#M1) and (#M2) are, respectively, displayed in the first and third columns of Figure 5. Their resistance patterns are, respectively, shown in the second and fourth columns, and they were obtained by utilizing the optimal source. In these optimized sources, the intensity distributions for #M1 were closer to a vertical quadrupole illumination mode. This is because, in the vertical direction, the pattern had a similar periodic grating structure to that of the vertical bipolar illumination mode. Horizontally, the middle area lacked pattern information for the horizontal illumination mode. For #M2, a 45-degree quadrupole illumination mode was more necessary for ensuring better image quality in the imaging process for optical lithography due to the periodic horizontal grating structure of the array. Figure 6 shows the simulation results when using two other methods with 150 iterations; these were a hybrid GA based on TS (TS-GA) and a method that combined the GA with the optimization strategy of PSO (PSO-GA). Figure 6a,e show the final source shapes for #M1 when utilizing TS-GA and PSO-GA, respectively. Similarly, Figure 6c,g were, respectively, achieved by using these two methods for #M2. Compared with the simulation results for the proposed method, the results shown in Figure 6 are not satisfactory. Although the tendencies of the source shapes were similar, there was a significant decrease in the intensity value for TS-GA and PSO-GA. Furthermore, it was obvious that the quality of the resistance patterns was not high enough when these optimal sources were employed as the lithographic imaging conditions. It was verified that the search process fell into a local optimum in the optimization process for these two methods. Furthermore, the convergence curves in the iterative process for these simulations are shown in Figure 7. Although the PEs obtained with PSO-GA decreased faster than those obtained with TS-GA for these two patterns before 25 iterations, those obtained with PSO-GA fell into a local optimum earlier than the others. This was because PSO had high efficiency early in the optimization process, which improved the overall optimization efficiency of PSO-GA. In comparison with these two simulation results, TS-GA had stronger convergence stability. On the contrary, the proposed method not only had a stable convergence efficiency, but also had a stronger optimization ability than that of the other two methods for PEs with #M1 and #M2. Hence, the proposed method had enhanced search capabilities and higher convergence efficiency. Appl. Sci. 2023, 13, x FOR PEER REVIEW 9 of 12  Figure 6 shows the simulation results when using two other methods with 150 iterations; these were a hybrid GA based on TS (TS-GA) and a method that combined the GA with the optimization strategy of PSO (PSO-GA). Figure 6a,e show the final source shapes for #M1 when utilizing TS-GA and PSO-GA, respectively. Similarly, Figure 6c,g were, respectively, achieved by using these two methods for #M2. Compared with the simulation results for the proposed method, the results shown in Figure 6 are not satisfactory. Although the tendencies of the source shapes were similar, there was a significant decrease in the intensity value for TS-GA and PSO-GA. Furthermore, it was obvious that the quality of the resistance patterns was not high enough when these optimal sources were employed as the lithographic imaging conditions. It was verified that the search process fell into a local optimum in the optimization process for these two methods. Furthermore, the convergence curves in the iterative process for these simulations are shown in Figure 7. Although the PEs obtained with PSO-GA decreased faster than those obtained with TS-GA for these two patterns before 25 iterations, those obtained with PSO-GA fell into a local optimum earlier than the others. This was because PSO had high efficiency early in the optimization process, which improved the overall optimization efficiency of PSO-GA. In comparison with these two simulation results, TS-GA had stronger convergence stability. On the contrary, the proposed method not only had a stable convergence efficiency, but also had a stronger optimization ability than that of the other two methods for PEs with #M1 and #M2. Hence, the proposed method had enhanced search capabilities and higher convergence efficiency.   (e-h) present the simulation results from TS-GA and PSO-GA, respectively; (a,c,e,g) are the source optimization results; (b,d,f,h) are the resistance patterns with the optimal source as the input condition in these simulations.

Conclusions
In this study, a high-efficiency hybrid GA method was proposed to generate the optimal intensity distribution of a source in the imaging process for optical lithography in order to achieve an acceptable image quality. To improve the optimization efficiency in the SO model, a GA was used as the main body in the optimization model. Meanwhile, to expand the search scope, an update strategy using PSO was employed to modify the values of variates in the population. In this complex optimization model, the tabu list method was utilized to record the individuals that appeared, making it possible to break away from local optima to avoid the premature convergence of the optimization model. Moreover, it was significantly important for the pattern errors (PEs) to be defined as a merit function for selecting the current optimal individual of the population. To reduce the computational complexity, a quarter region of the source was used as the optimization target, and it could be flipped to form the integral source with a mirror operation. Furthermore, to improve the convergence efficiency, two other methods were employed to optimize the lithographic source shape, namely, TS-GA and PSO-GA. The simulation results confirmed that the proposed SO method based on the hybrid GA had efficient optimization

Conclusions
In this study, a high-efficiency hybrid GA method was proposed to generate the optimal intensity distribution of a source in the imaging process for optical lithography in order to achieve an acceptable image quality. To improve the optimization efficiency in the SO model, a GA was used as the main body in the optimization model. Meanwhile, to expand the search scope, an update strategy using PSO was employed to modify the values of variates in the population. In this complex optimization model, the tabu list method was utilized to record the individuals that appeared, making it possible to break away from local optima to avoid the premature convergence of the optimization model. Moreover, it was significantly important for the pattern errors (PEs) to be defined as a merit function for selecting the current optimal individual of the population. To reduce the computational complexity, a quarter region of the source was used as the optimization target, and it could be flipped to form the integral source with a mirror operation. Furthermore, to improve the convergence efficiency, two other methods were employed to optimize the lithographic source shape, namely, TS-GA and PSO-GA. The simulation results confirmed that the proposed SO method based on the hybrid GA had efficient optimization performance when improving the image quality during the imaging process for optical lithography.