## 1.

## Introduction

Optimization techniques play an important role in the improvement of the pattern fidelity and depth of focus (DOF) of current optical lithography systems. Source optimization is mainly aimed at altering the source shape,^{1}^{,}^{2} whereas mask optimization is aimed at modulating the amplitude of the electric field to compensate for the optical proximity effect.^{3}^{,}^{4} In 2002, Rosenbluth et al. proposed the first source mask optimization (SMO) method that exploits the synergy between the source and mask to achieve a higher resolution.^{5} Since then, a number of SMO methods have been proposed in the literature.^{6}7.8.^{–}^{9} Most methods are based on a scalar imaging model that is no longer accurate for a numerical aperture $(\mathrm{NA})>0.6$.^{10} In high-NA immersion lithography systems, the vector nature of the electromagnetic field must be taken into account. Recently, we proposed a pixelated SMO based on a vector imaging model that significantly improved the simulation precision for lithography at the 45-nm node and beyond.^{11}^{,}^{12}

Previous SMO methods fixed the NA and fell short in considering the mutual impact of the NA with respect to the source and mask. Prior work has demonstrated that a larger NA could realize a higher resolution, but the DOF would decrease because of the relation $\mathrm{DOF}={k}_{2}\xb7\lambda /{(\mathrm{NA})}^{2}$,^{13} where $\lambda $ is the wavelength and ${k}_{2}$ is the process factor. Hence, it is highly important to pursue the optimal NA during the SMO procedure so as to realize the designated resolution and achieve high image fidelity within a larger DOF. In addition, current pixelated SMO methods dramatically increase the complexity and fabrication cost of the optimized source and mask patterns; thus, they suffer from an inherent disadvantage in manufacturing.^{14}^{,}^{15}

To overcome these limitations, this paper proposes a parametric source-mask-NA co-optimization (SMNO) method to improve the pattern fidelity within a large DOF and the complexity of the source and mask patterns. To our knowledge, this paper is the first to solve for the parametric SMNO problem based on a vector imaging model. First, the vector imaging model described in Refs. 11 and 16 is used to formulate the SMNO framework, which significantly improves the simulation precision for the 45-nm node in immersion lithography. Then, an analytic model of the parametric source, mask, and NA is built. Since the derivative of the arch function exists, it is used to approximate the parametric source pattern, phase-shifting mask (PSM) pattern, and NA values such that the gradient-based algorithm can be applied to the SMNO problem. In particular, the source is modeled by its partial coherent factor and opening angle. The mask is represented by the main feature and the serif. During the SMNO process, all parameters are simultaneously optimized by using the steepest-descent algorithm. In order to validate the proposed SMNO algorithm, simulations based on a quasar source and a two-dimensional mask are presented as an example. The simulations show that, in comparison with the SMO method, the proposed SMNO method can achieve the optimal combination of source, mask, and NA to achieve superior imaging performance over a wider DOF. In addition, the parametric source can maintain an extremely simplistic distribution to avoid customizing the diffractive optical elements or installing an expensive pixel source generator. The parametric mask greatly reduces the complexity, which is beneficial for mask fabrication and lowers the cost.

The remainder of the paper is organized as follows. Section 2 briefly summarizes the vector imaging model used in this paper. The analytic and parametric source, mask, and NA are modeled in Sec. 3. The SMNO algorithm is proposed and described at length in Sec. 4. The simulations are provided in Sec. 5. Finally, conclusions are drawn in Sec. 6.

## 2.

## Vector Imaging Model for Immersion Lithography

In this paper, we choose the vector imaging model described in Refs. 11 and 16 as the basis for developing the SMNO framework. This vector imaging model can provide accurate simulation results even when $\mathrm{NA}>1$ in immersion lithography. The accuracy of the model has been proven by comparison with PROLITH.^{11}

In the following part, we use $(x,y)$ and $(f,g)$ to represent the coordinate systems in the spatial and frequency domains. $(\alpha ,\beta ,\gamma )$ is the direction cosine of the light propagating through $(f,g)$. The direction cosines on the source side, mask side, and image side are denoted as $({\alpha}_{s},{\beta}_{s},{\gamma}_{s})$, $(\alpha ,\beta ,\gamma )$, and $({\alpha}_{i},{\beta}_{i},{\gamma}_{i})$, respectively.

A schematic of the lithography imaging process is illustrated in Fig. 1. A source point $S({\alpha}_{s},{\beta}_{s})$ in the source plane emits the polarized light ${\mathbf{E}}_{\text{in}}=[{E}_{x};{E}_{y}]$ propagating in the direction $\widehat{k}$ and incident on the mask $M(x,y)$, where $\widehat{k}=({\alpha}_{s},{\beta}_{s},{\gamma}_{s})$ is the direction cosine with ${\gamma}_{s}=\sqrt{1-{\alpha}_{s}^{2}-{\beta}_{s}^{2}}$. The electric field on the exit pupil ${\mathbf{E}}_{\mathrm{ext}}$ can be expressed as

## (1)

$${\mathbf{E}}_{\mathrm{ext}}({\alpha}_{i},{\beta}_{i})={Q}_{\mathrm{NA}}({\alpha}_{i},{\beta}_{i})\xb7C({\alpha}_{i},{\beta}_{i})\xb7\mathrm{\Psi}({\alpha}_{i},{\beta}_{i})\phantom{\rule{0ex}{0ex}}\xb7\mathbb{F}\{M(x,y)\}\xb7{\mathbf{E}}_{\text{in}}\xb7{e}^{i\delta},$$## (2)

$${Q}_{\mathrm{NA}}({\alpha}_{i},{\beta}_{i})=\{\begin{array}{ll}1,& \sqrt{{\alpha}_{i}^{2}+{\beta}_{i}^{2}}\le \frac{\mathrm{NA}}{{n}_{i}}\\ 0,& \text{otherwise}\end{array},$$The electric field at the image plane can be written as

## (3)

$${\mathbf{E}}_{\mathrm{im}}({x}_{i},{y}_{i},{z}_{i})=\frac{2\pi}{{n}_{i}\xb7R}\xb7{\mathbb{F}}^{-1}\{{\mathbf{E}}_{\mathrm{ext}}({\alpha}_{i},{\beta}_{i})\},$$The complete aerial image intensity ${I}_{\mathrm{im}}$ of the partial coherent source $S$ can be obtained by

## (4)

$${I}_{\mathrm{im}}({x}_{i},{y}_{i},{z}_{i})=\iint S({\alpha}_{s},{\beta}_{s})\xb7[{\mathbf{E}}_{\mathrm{im}}({x}_{i},{y}_{i},{z}_{i})\phantom{\rule{0ex}{0ex}}\xb7{\mathbf{E}}_{\mathrm{im}}^{*}({x}_{i},{y}_{i},{z}_{i})]\mathrm{d}{\alpha}_{s}\mathrm{d}{\beta}_{s}.$$In order to obtain unit intensity in the aerial image,^{17} the normalized aerial image intensity ${I}_{\mathrm{AI}}$ can be written as

## (5)

$${I}_{\mathrm{AI}}({x}_{i},{y}_{i},{z}_{i})=\frac{{I}_{\mathrm{im}}({x}_{i},{y}_{i},{z}_{i})}{{I}_{\text{clear}}({x}_{i},{y}_{i},{z}_{i})},$$The aerial image represents the distribution on the wafer plane of the optical intensity that will cause the exposure of the resist. Usually, the exposure dose is described by the aerial image threshold value ${r}_{\mathrm{t}}$ in a constant threshold resist model.^{18} The resist can be developed when the aerial image intensity ${I}_{\mathrm{AI}}$ is larger than ${r}_{\mathrm{t}}$. For numerical consideration, we employ a sigmoid function instead of a hard threshold to calculate the resist image.^{19} Then, the exposure resist image $Z$ can be expressed as

## 3.

## Analytic Model of the Parametric Source, Mask, and NA

In the aerial image equation, the source, mask, and NA are all binarized and have been formulated using a rectangular function. We often avoid using the rectangular function in optimization because it is not differentiable. In this section, we use a differentiable function instead of the rectangular function to describe the source, mask, and NA as the use of a differentiable function is more convenient for optimization in inverse lithography. Let us introduce the name “arch function” for the following rational approximation of a rectangular function:

where $b$ is the steepness. As shown in Fig. 2, the arch function is approximately equivalent to the rectangular function when $b$ is sufficiently large.## 3.1.

### Source Model

Let $S({\alpha}_{s},{\beta}_{s})\in {\mathbb{R}}^{{N}_{s}\times {N}_{s}}$ denote the source with all entries equal to 0 or 1, where ${N}_{S}$ is the dimension of the source. The parametric source can be described by a partial coherent factor, such as the outer sigma ${\sigma}_{\text{out}}$, the sigma width ${\sigma}_{\text{width}}$, and the opening angle $\theta $. For annular illumination, the source shape is described by the rectangular function.

## (8)

$${S}_{A}({\alpha}_{s},{\beta}_{s};{\sigma}_{\text{out}},{\sigma}_{\text{width}})=\{\begin{array}{ll}1,& ({\sigma}_{\text{out}}-{\sigma}_{\text{width}})\le {r}_{s}\le {\sigma}_{\text{out}}\\ 0,& \text{otherwise}\end{array},$$In order to calculate the derivative of the source parameter, the source can be modeled by a two-dimensional arch function. The annular illumination ${S}_{A}$ can be formulated as

## (9)

$${S}_{A}({\alpha}_{s},{\beta}_{s};{\sigma}_{\text{out}},{\sigma}_{\text{width}})=\frac{1}{{({r}_{s}/{\sigma}_{\text{out}})}^{2{b}_{s}}+1}-\frac{1}{{[{r}_{s}/({\sigma}_{\text{out}}-{\sigma}_{\text{width}})]}^{2{b}_{s}}+1},$$For quasar illumination, the source is also restricted by the opening angle $\theta $. The quasar illumination ${S}_{Q}$ can be modeled by the product of ${S}_{A}$ and a radial function ${K}_{R}$.

## (10)

$${S}_{Q}({\alpha}_{s},{\beta}_{s};{\sigma}_{\text{out}},{\sigma}_{\text{width}};{\theta}_{1},{\theta}_{2})\phantom{\rule{0ex}{0ex}}={K}_{R}({\alpha}_{s},{\beta}_{s};{\theta}_{1},{\theta}_{2})\xb7{S}_{A}({\alpha}_{s},{\beta}_{s};{\sigma}_{\text{out}},{\sigma}_{\text{width}}),$$## (11)

$${K}_{R}({\alpha}_{s},{\beta}_{s};{\theta}_{1},{\theta}_{2})=1-\frac{1}{{\left(\frac{2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180}{180-{\theta}_{1}}\right)}^{2{b}_{s}}+1}\phantom{\rule{0ex}{0ex}}+\frac{1}{{\left(\frac{2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180}{{\theta}_{2}}\right)}^{2{b}_{s}}+1},$$## 3.2.

### Mask Model

Let $M(x,y)\in {\mathbb{R}}^{N\times N}$ be the mask with all entries equal to ${t}_{m}\xb7\mathrm{exp}(i{\phi}_{m})$, where $N$ is the dimension of the mask, and ${t}_{m}$ and ${\phi}_{m}$ are the feature transmittance and phase shift of the attenuated PSM (AttPSM), respectively. For a 6% AttPSM, ${t}_{m}=\sqrt{0.06}$ and ${\phi}_{m}=\pi $ in the phase-shift region, whereas ${t}_{m}=1$ and ${\phi}_{m}=0$ in the non-phase-shift region.

The parametric mask is composed of the main feature and the serif, where the serif is placed at the corner of the main feature. The main feature is controlled by the feature width and feature height, while the serif is modeled by the serif size and the serif offset. The serif offset is the distance that the serif deviates from the main feature. The schematic of the parametric mask is shown in Fig. 4.

For the binary mask with all entries equal to 0 or 1, the $l$’th main feature can be expressed by using the arch function.

## (12)

$${m}_{F}^{l}(x,y;{m}_{w}^{l},{m}_{h}^{l})=1-\frac{1}{{\left[\frac{2(x-{x}_{0}^{l})}{{m}_{w}^{l}}\right]}^{2{b}_{m}}+1}\xb7\frac{1}{{\left[\frac{2(y-{y}_{0}^{l})}{{m}_{h}^{l}}\right]}^{2{b}_{m}}+1},$$The $j$’th serif can be expressed as

## (13)

$${m}_{S}^{j}(x,y;{m}_{\mathrm{ss}}^{j},{m}_{\mathrm{so}}^{j})\phantom{\rule{0ex}{0ex}}=1-\frac{1}{{\left[\frac{2(x-{x}_{0}^{j}-{m}_{\mathrm{so}}^{j})}{{m}_{\mathrm{ss}}^{j}}\right]}^{2{b}_{m}}+1}\xb7\frac{1}{{\left[\frac{2(y-{y}_{0}^{j}-{m}_{\mathrm{so}}^{j})}{{m}_{\mathrm{ss}}^{j}}\right]}^{2{b}_{m}}+1},$$The function of the whole binary mask ${m}_{0}$ can then be derived from the summation of the main feature and the serif.

## (14)

$${m}_{0}(x,y;{m}_{w},{m}_{h},{m}_{\mathrm{ss}},{m}_{\mathrm{so}})\phantom{\rule{0ex}{0ex}}=\sum _{\chi =1}^{l}{m}_{F}^{\chi}(x,y;{m}_{w}^{\chi},{m}_{h}^{\chi})+\sum _{\zeta =1}^{j}{m}_{S}^{\zeta}(x,y;{m}_{\mathrm{ss}}^{\zeta},{m}_{\mathrm{so}}^{\zeta})-l-j+1,\phantom{\rule{0ex}{0ex}}$$Because the mask distribution is binarized, we convert ${m}_{0}$ into binary when simulating the aerial and resist image.

It is noted that we use ${m}_{0}$ in the process of calculating the derivative of the main feature and serif parameters, and use ${m}_{0}^{\prime}$ instead of ${m}_{0}$ when performing the aerial and resist image simulation.

The AttPSM can be modeled by the product of the attenuated layer ${M}_{\mathrm{att}}({m}_{0})$ and the phase shift layer $\mathrm{\Phi}({m}_{0})$.

## (16)

$$M(x,y;{m}_{w},{m}_{h},{m}_{\mathrm{ss}},{m}_{\mathrm{so}})={M}_{\mathrm{att}}({m}_{0})\xb7\mathrm{\Phi}({m}_{0}),$$Thus, each independent main feature ${M}_{F}^{l}$ and serif ${M}_{S}^{j}$ of the AttPSM can be written as

## 3.3.

### NA Model

NA can also be modeled by the differentiable arch function to replace the circle function of Eq. (2). Then, ${Q}_{\mathrm{NA}}({\alpha}_{i},{\beta}_{i})$ can be written as

## (19)

$${Q}_{\mathrm{NA}}({\alpha}_{i},{\beta}_{i})=\frac{1}{{\left(\frac{\sqrt{{\alpha}_{i}^{2}+{\beta}_{i}^{2}}}{\mathrm{NA}}\right)}^{2{b}_{\mathrm{na}}}+1},$$Just as with the binarization process of the mask in Eq. (15), NA should also be binarized as ${Q}_{\mathrm{NA}}^{\prime}$ while performing the aerial and resist image simulation.

## (20)

$${Q}_{\mathrm{NA}}^{\prime}=\{\begin{array}{ll}1,& {Q}_{\mathrm{NA}}>0.5\\ 0,& {Q}_{\mathrm{NA}}\le 0.5\end{array}.$$${Q}_{\mathrm{NA}}$ in Eq. (19) is only used to calculate the derivative of NA.

## 4.

## SMNO Based on the Vector Imaging Model

## 4.1.

### Cost Function

Given a binary target pattern $\overline{Z}\in {\mathbb{R}}^{N\times N}$ with all entries equal to 0 or 1, the pattern error (PE) ${F}_{Z}^{\prime}$ can be expressed as the difference between $\overline{Z}$ and the resist image $Z$.

## (21)

$${F}_{Z}^{\prime}={\Vert Z-\overline{Z}\Vert}_{2}^{2}=\sum _{\tau =1}^{{N}^{2}}{({z}_{\tau}-{\overline{z}}_{\tau})}^{2},$$In order to improve the convergence of the optimization, we add the difference between the aerial image ${I}_{\mathrm{AI}}$ and the target pattern $\overline{Z}$ in the cost function.

## (22)

$${F}_{\mathrm{AI}}^{\prime}={\Vert {I}_{\mathrm{AI}}-c\xb7\overline{Z}\Vert}_{2}^{2}=\sum _{\tau =1}^{{N}^{2}}{({I}_{\mathrm{AI},\tau}-c\xb7{\overline{z}}_{\tau})}^{2},$$Thus, we combine Eqs. (21) and (22) to form the cost function.

where ${\omega}_{g}\in [\mathrm{0,1}]$ is the weight of ${F}_{Z}^{\prime}$.In order to maintain high pattern fidelity over a large range of DOF, the final cost function $D$ should be adjusted by adding the off-focus term. Thus,

## (24)

$$D={\omega}_{\mathrm{foc}}\xb7{F}_{\mathrm{foc}}+(1-{\omega}_{\mathrm{foc}})\xb7{F}_{\mathrm{defoc}},$$## 4.2.

### SMNO Algorithm

The SMNO algorithm can be formulated as the search for the optimal source, mask, and NA to minimize the cost function $D$, such that

In order to find the best combination of the source, mask, and NA, the steepest-descent method is used to implement the proposed SMNO algorithm. Then, the $k$’th iterative parameter set ${\{{P}_{n}\}}^{(k)}$, which includes all parameters of the source, mask, and NA, can be calculated by

where $P$ represents the parameters of the source, mask, and NA. $\mathrm{\Lambda}$ is a vector that represents the step length of the parameters. The step length should be assigned before the optimization for each individual parameter $P$. The steepest-descent direction is denoted by $\mathbf{d}$, which can be derived byThe partial derivative of the cost function $D$ with respect to each parameter $P$ of the source, mask, and NA can be written as

## (28)

$$\frac{\partial D}{\partial P}={\omega}_{\mathrm{foc}}\xb7\frac{\partial {F}_{\mathrm{foc}}}{\partial P}+(1+{\omega}_{\mathrm{foc}})\xb7\frac{\partial {F}_{\mathrm{defoc}}}{\partial P}.$$Then, for the derivative of the cost function $F$ at the focal and defocus positions,

## (29)

$$\frac{\partial F}{\partial P}={\omega}_{g}\xb7\frac{\partial {F}_{Z}^{\prime}}{\partial P}+(1-{\omega}_{g})\xb7\frac{\partial {F}_{\mathrm{AI}}^{\prime}}{\partial P},$$## (30)

$$\frac{\partial {F}_{Z}^{\prime}}{\partial P}=\sum _{\tau =1}^{{N}^{2}}(2\{\frac{1}{\mathrm{exp}[-a({I}_{\mathrm{AI}}-{r}_{\mathrm{t}})]+1}-{\overline{z}}_{m}\}\phantom{\rule{0ex}{0ex}}\xb7\frac{-1}{{\{\mathrm{exp}[-a({I}_{\mathrm{AI}}-{r}_{\mathrm{t}})]+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\mathrm{exp}[-a({I}_{\mathrm{AI}}-{r}_{\mathrm{t}})]\xb7(-a)\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial P}),$$## (31)

$$\frac{\partial {F}_{\mathrm{AI}}^{\prime}}{\partial P}=\sum _{\tau =1}^{{N}^{2}}[2({I}_{\mathrm{AI},\tau}-c\xb7{\overline{z}}_{\tau})\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial P}].$$From Eqs. (28)–(31), we find that $\partial D/\partial P$ can be derived once $\partial {I}_{\mathrm{AI}}/\partial P$ is calculated. Therefore,

## (32)

$$\frac{\partial D}{\partial P}=\frac{\partial D}{\partial {I}_{\mathrm{AI}}}\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial P}.$$It is noted that all the parameters of source, mask, and NA are constrained by boundaries in lithography. In order to reduce the bound-constrained optimization problem to an unconstrained optimization problem, we adopt the following parametric transformation to convert $P$ into $\mathrm{\Omega}$.

## (33)

$$\mathrm{\Omega}=\frac{-1}{{b}_{p}}\xb7\mathrm{ln}(\frac{{P}_{\mathrm{max}}-{P}_{\mathrm{min}}}{P-{P}_{\mathrm{min}}}-1),$$## (34)

$$P=\frac{{P}_{\mathrm{max}}-{P}_{\mathrm{min}}}{\mathrm{exp}(-{b}_{p}\xb7\mathrm{\Omega})+1}+{P}_{\mathrm{min}}.$$Then, the derivative of the cost function $D$ with respect to the parameter $P$ can be converted to the derivative of $D$ with respect to $\mathrm{\Omega}$.

## (35)

$$\frac{\partial D}{\partial \mathrm{\Omega}}=\frac{\partial D}{\partial P}\xb7\frac{\partial P}{\partial \mathrm{\Omega}},$$## (36)

$$\frac{\partial P}{\partial \mathrm{\Omega}}=\frac{{b}_{p}\xb7({P}_{\mathrm{max}}-{P}_{\mathrm{min}})\xb7\mathrm{exp}(-{b}_{p}\xb7\mathrm{\Omega})}{{[\mathrm{exp}(-{b}_{p}\xb7\mathrm{\Omega})+1]}^{2}}.$$Therefore, the derivative $\partial {I}_{\mathrm{AI}}/\partial P$ can also be calculated by $\partial {I}_{\mathrm{AI}}/\partial \mathrm{\Omega}$. In the following part, we will calculate the derivatives of the parameters of the source, mask, and NA.

## 4.2.1.

#### Derivative of the source parameter

The derivative of the cost function $D$ with respect to the source parameter $\mathrm{\Omega}({s}_{v})$ is

## (37)

$$\frac{\partial D}{\partial \mathrm{\Omega}({s}_{v})}=\frac{\partial D}{\partial {I}_{\mathrm{AI}}}\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}({s}_{v})},$$## (38)

$$\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}({s}_{v})}=\frac{{I}_{\text{clear}}\xb7\partial {I}_{\mathrm{im}}/\partial \mathrm{\Omega}({s}_{v})-{I}_{\mathrm{im}}\xb7\partial {I}_{\text{clear}}/\partial \mathrm{\Omega}({s}_{v})}{{({I}_{\text{clear}})}^{2}}.$$We use ${I}_{\mathrm{im}}^{\prime}$ to represent ${I}_{\mathrm{im}}$ and ${I}_{\text{clear}}$; then, $\partial {I}_{\mathrm{im}}^{\prime}/\partial \mathrm{\Omega}({s}_{v})$ can be derived from Eq. (4).

## (39)

$$\frac{\partial {I}_{\mathrm{im}}^{\prime}}{\partial \mathrm{\Omega}({s}_{v})}=\iint \frac{\partial S}{\partial \mathrm{\Omega}({s}_{v})}\xb7[\mathbf{E}({x}_{i},{y}_{i},{z}_{i})\xb7{\mathbf{E}}^{*}({x}_{i},{y}_{i},{z}_{i})]\mathrm{d}{\alpha}_{s}\mathrm{d}{\beta}_{s},$$## 4.2.2.

#### Derivative of the mask parameter

The derivative of the cost function $D$ with respect to the mask parameter $\mathrm{\Omega}({m}_{z})$ is

## (40)

$$\frac{\partial D}{\partial \mathrm{\Omega}({m}_{z})}=\frac{\partial D}{\partial {I}_{\mathrm{AI}}}\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}({m}_{z})},$$## (41)

$$\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}({m}_{z})}=\frac{\partial {I}_{\mathrm{im}}/\partial \mathrm{\Omega}({m}_{z})}{{I}_{\text{clear}}},$$## (42)

$$\frac{\partial {I}_{\mathrm{im}}}{\partial \mathrm{\Omega}({m}_{z})}=\iint S({\alpha}_{s},{\beta}_{s})\xb7[\frac{\partial {\mathbf{E}}_{\mathrm{im}}}{\partial \mathrm{\Omega}({m}_{z})}\xb7{\mathbf{E}}_{\mathrm{im}}^{*}+{\mathbf{E}}_{\mathrm{im}}\xb7\frac{\partial {\mathbf{E}}_{\mathrm{im}}^{*}}{\partial \mathrm{\Omega}({m}_{z})}]\mathrm{d}{\alpha}_{s}\mathrm{d}{\beta}_{s}.$$For ${\mathbf{E}}_{\mathrm{im}}$ and its complex conjugate ${\mathbf{E}}_{\mathrm{im}}^{*}$,

## (43)

$$\frac{\partial {\mathbf{E}}_{\mathrm{im}}^{*}}{\partial \mathrm{\Omega}({m}_{z})}={\left[\frac{\partial {\mathbf{E}}_{\mathrm{im}}}{\partial \mathrm{\Omega}({m}_{z})}\right]}^{*},$$## (44)

$$\frac{\partial {\mathbf{E}}_{\mathrm{im}}}{\partial \mathrm{\Omega}({m}_{z})}=\frac{2\pi}{{n}_{i}\xb7R}\xb7{\mathbb{F}}^{-1}\{{Q}_{\mathrm{NA}}({\alpha}_{i},{\beta}_{i})\xb7C({\alpha}_{i},{\beta}_{i})\xb7\mathrm{\Psi}({\alpha}_{i},{\beta}_{i})\phantom{\rule{0ex}{0ex}}\xb7\mathbb{F}\{\frac{\partial M}{\partial \mathrm{\Omega}({m}_{z})}\}\xb7{\mathbf{E}}_{\text{in}}\xb7{e}^{i\delta}\}.$$For each individual feature on the mask $M$,

## (45)

$$\frac{\partial M}{\partial \mathrm{\Omega}({m}_{z})}=\frac{\partial {M}_{V}}{\partial \mathrm{\Omega}({m}_{z})},$$## (46)

$$\frac{\partial {M}_{V}}{\partial \mathrm{\Omega}({m}_{z})}=(1-t)\xb7\frac{\partial {m}_{V}}{\partial \mathrm{\Omega}({m}_{z})}\xb7\mathrm{exp}[i\pi \xb7(1-{m}_{V})]\phantom{\rule{0ex}{0ex}}-[(1-t)\xb7{m}_{V}+t]\xb7i\pi \phantom{\rule{0ex}{0ex}}\xb7\mathrm{exp}[i\pi \xb7(1-{m}_{V})]\xb7\frac{\partial {m}_{V}}{\partial \mathrm{\Omega}({m}_{z})},$$The derivatives of the mask parameters $\partial {m}_{V}/\partial \mathrm{\Omega}({m}_{z})$ are listed in Appendix B. By substituting $\partial {m}_{V}/\partial \mathrm{\Omega}({m}_{z})$ into Eq. (46), the derivative of the cost function $D$ with respect to the mask parameter $\mathrm{\Omega}({m}_{z})$ can be solved.

## 4.2.3.

#### Derivative of the NA

The derivative of the cost function $D$ with respect to $\mathrm{\Omega}(\mathrm{NA})$ is

## (47)

$$\frac{\partial D}{\partial \mathrm{\Omega}(\mathrm{NA})}=\frac{\partial D}{\partial {I}_{\mathrm{AI}}}\xb7\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}(\mathrm{NA})}.$$The term $\partial {I}_{\mathrm{AI}}/\partial \mathrm{\Omega}(\mathrm{NA})$ can be derived from Eq. (5).

## (48)

$$\frac{\partial {I}_{\mathrm{AI}}}{\partial \mathrm{\Omega}(\mathrm{NA})}=\frac{{I}_{\text{clear}}\xb7\partial {I}_{\mathrm{im}}/\partial \mathrm{\Omega}(\mathrm{NA})-{I}_{\mathrm{im}}\xb7\partial {I}_{\text{clear}}/\partial \mathrm{\Omega}(\mathrm{NA})}{{({I}_{\text{clear}})}^{2}},$$## (49)

$$\frac{\partial I}{\partial \mathrm{\Omega}(\mathrm{NA})}\phantom{\rule{0ex}{0ex}}=\iint S({\alpha}_{s},{\beta}_{s})\xb7[\frac{\partial \mathbf{E}}{\partial \mathrm{\Omega}(\mathrm{NA})}\xb7{\mathbf{E}}^{*}+\mathbf{E}\xb7\frac{\partial {\mathbf{E}}^{*}}{\partial \mathrm{\Omega}(\mathrm{NA})}]\mathrm{d}{\alpha}_{s}\mathrm{d}{\beta}_{s}.$$In Eq. (49), $I$ represents both ${I}_{\mathrm{im}}$ and ${I}_{\text{clear}}$, while $\mathbf{E}$ represents the electric field of the mask $M$ and the mask ${M}_{\text{clear}}$. $\partial \mathbf{E}/\partial \mathrm{\Omega}(\mathrm{NA})$ can be derived from Eq. (3):

## (50)

$$\frac{\partial \mathbf{E}}{\partial \mathrm{\Omega}(\mathrm{NA})}=\frac{2\pi}{{n}_{i}\xb7R}\xb7{\mathbb{F}}^{-1}\{\frac{\partial {Q}_{\mathrm{NA}}}{\partial \mathrm{\Omega}(\mathrm{NA})}\xb7C({\alpha}_{i},{\beta}_{i})\xb7\mathrm{\Psi}({\alpha}_{i},{\beta}_{i})\phantom{\rule{0ex}{0ex}}\xb7\mathbb{F}\{M(x,y)\}\xb7{\mathbf{E}}_{\text{in}}\xb7{e}^{i\delta}\}.$$The term $\partial {Q}_{\mathrm{NA}}/\partial \mathrm{\Omega}(\mathrm{NA})$ can be derived from Eq. (19).

## (51)

$$\frac{\partial {Q}_{\mathrm{NA}}}{\partial \mathrm{\Omega}(\mathrm{NA})}=\frac{2{b}_{\mathrm{na}}{\left(\frac{\sqrt{{\alpha}_{i}^{2}+{\beta}_{i}^{2}}}{\mathrm{NA}}\right)}^{2{b}_{\mathrm{na}}-1}}{{[{\left(\frac{\sqrt{{\alpha}_{i}^{2}+{\beta}_{i}^{2}}}{\mathrm{NA}}\right)}^{2{b}_{\mathrm{na}}}+1]}^{2}}\xb7\frac{\sqrt{{\alpha}_{i}^{2}+{\beta}_{i}^{2}}}{{(\mathrm{NA})}^{2}}\xb7\frac{\partial \mathrm{NA}}{\partial \mathrm{\Omega}(\mathrm{NA})}.$$Finally, the derivative of the cost function $D$ with respect to the NA $\partial D/\partial \mathrm{\Omega}(\mathrm{NA})$ can be solved.

The steepest-descent direction $\mathbf{d}$ can be calculated after all the derivatives have been solved. Finally, the parameters can be optimized by using the steepest-descent method iteratively until the termination condition is satisfied.

## 5.

## Implementation of Parametric SMNO

In order to demonstrate the validity of the proposed optimization method, we illustrate the simulation results for the SMNO in this section. For comparison consideration, a parametric SMO without NA optimization is also performed in the simulations.

The lithography system is an argon fluoride (ArF) immersion lithography system with variable NA. The reflection index of the immersion medium is $n=1.44$. The wavelength is $\lambda =193\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$. The reduction of the projector is $R=4$.

The simulations use the parametric source and mask. The source shape is quasar with $Y$ polarization. Four source parameters, ${\sigma}_{\text{out}}$, ${\sigma}_{\text{width}}$, ${\theta}_{1}$, and ${\theta}_{2}$, are optimized. The source array dimension ${N}_{S}$ is 41. A 45-nm AttPSM line-space pattern is used as a target in the simulation. There are three lines in the $600\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}\times 600\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$ wafer, while the mask array dimension $N$ is 600. As shown in Fig. 5(a), each line feature is 45 nm in width and 360 nm in height, while the space between the lines is 90 nm. The feature transmittance is ${t}_{m}=\sqrt{0.06}$, with a 180 deg phase shift to enhance the resolution. With the modeling in Sec. 3.2, we can change the feature width, feature height, and serif to compensate for the optical proximity effects. All the main features and serifs are numbered in Fig. 5(b). Because the mask pattern has fourfold symmetry, we need only optimize the mask parameters on the top-left quarter, which simplifies the optimization procedure. Thus, 10 mask parameters are optimized, including ${m}_{w}^{1}$ and ${m}_{h}^{1}$ of main feature 1, ${m}_{w}^{2}$ and ${m}_{h}^{2}$ of main feature 2, ${m}_{\mathrm{ss}}^{4}$ and ${m}_{\mathrm{so}}^{4}$ of serif 4, ${m}_{\mathrm{ss}}^{5}$ and ${m}_{\mathrm{so}}^{5}$ of serif 5, and ${m}_{\mathrm{ss}}^{6}$ and ${m}_{\mathrm{so}}^{6}$ of serif 6.

There are 15 parameters in optimization, 4 source parameters, 10 mask parameters, and NA. Nonoptimized initial parameter values are given to serve as a starting point for optimization. The initial parameter values, minimum value ${P}_{\mathrm{min}}$, and maximum value ${P}_{\mathrm{max}}$ are listed in Table 1. The step lengths $\mathrm{\Lambda}$ in Eq. (26) are 0.01, 0.1, and 0.02 for source, mask, and NA parameters, respectively. The steepness of the parameters are all equal to 40 for ${b}_{s}$, ${b}_{m}$, and ${b}_{\mathrm{na}}$. The constant ${b}_{p}$ in Eq. (33) is 0.05.

## Table 1

Parameter settings, and source mask optimization (SMO) and source-mask-numerical aperture co-optimization (SMNO) results.

Parameters | Initial | Pmin | Pmax | SMO result | SMNO result | |
---|---|---|---|---|---|---|

Source | σout | 0.78 | 0.6 | 1 | 0.631 | 0.734 |

σwidth | 0.15 | 0.1 | 0.3 | 0.106 | 0.128 | |

θ1 (deg) | 45 | 30 | 60 | 30.033 | 30.085 | |

θ2 (deg) | 45 | 30 | 60 | 30 | 30 | |

Mask | mw1 (nm) | 45 | 40 | 100 | 52.337 | 59.386 |

mh1 (nm) | 360 | 350 | 480 | 425.3 | 406.847 | |

mw2 (nm) | 45 | 40 | 100 | 51.098 | 56.567 | |

mh2 (nm) | 360 | 350 | 480 | 425.333 | 429.145 | |

mss4 (nm) | 20 | 0 | 50 | 28.4 | 33.5 | |

mso4 (nm) | 10 | −10 | 20 | 14.234 | 14.662 | |

mss5 (nm) | 20 | 0 | 50 | 22.285 | 24.904 | |

mso5 (nm) | 10 | −10 | 20 | 15.096 | 14.146 | |

mss6 (nm) | 20 | 0 | 50 | 27.691 | 25.484 | |

mso6 (nm) | 10 | −10 | 20 | 14.582 | 14.148 | |

NA | NA | 1.3 | 1 | 1.35 | — | 1.013 |

The steepness ${a}_{r}=40$ and the image threshold ${r}_{\mathrm{t}}=0.25$ when calculating the resist image. In order to increase the DOF, the resist images at the best focus position ($\delta =0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$) and the 150-nm defocus position ($\delta =150\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$) are simulated. We assign ${\omega}_{\mathrm{foc}}=0.3$ in Eq. (24) and ${\omega}_{g}=0.8$ in Eq. (23). The constant $c$ in Eq. (22) is 1. The optimization terminates after 100 iterations.

The simulations are performed using the optimization algorithm in Sec. 4. The simulation for SMO uses the same optimization condition with SMNO except for NA. Simulation results for the SMO and SMNO are listed in the last two rows of Table 1 and are also shown in Fig. 6. The results show that both source and mask parameters are optimized after SMO and SMNO. It is noted that the optimal NA for this kind of target pattern is 1.013, which can realize the designated resolution and maintain high pattern fidelity within a large DOF.

The optimized source, mask, and resist image at the best focus, 100-nm defocus, and 150-nm defocus positions are all shown in Fig. 6. The red solid line in the figure indicates the target pattern. The PEs are also marked on the top of each resist image. From Fig. 6, one can see that without optimization, the mask can only be printed on the wafer at the best focus position with a large PE of 22,656, whereas no pattern can be printed at the 100-nm and 150-nm defocus positions. When the parametric source and mask have been optimized by using SMO, the mask can be printed at both the best focus and the defocus positions. The resist image is highly faithful to the target pattern with a PE of 5800 at the best focus. However, the PE equals 36,168 at the 150-nm defocus position, which is too large to maintain pattern fidelity. When performing SMNO, we found that the resist image could maintain high pattern fidelity through a large DOF. The PE of SMNO at the best focus, 100-nm defocus, and 150-nm defocus positions are 8038, 7430, and 14,186, respectively. Although the PE of SMNO is larger than that of SMO by 2238 at the best focus position, the PE values at the defocus positions are significantly less. The PE of SMNO at the 150-nm defocus position is 656 and 21,982 smaller than the PE of SMO at the 100-nm and 150-nm defocus positions. Figure 7 shows a comparison of the PEs at various defocus positions, where we found that the PE of SMO increases drastically at a defocus position, but SMNO maintains a low level of PE over a large range of DOF. The comparison of optimization results between the SMO and SMNO demonstrates that SMNO could effectively improve the pattern fidelity and enlarge the DOF, which also reveals that it is necessary to include NA in the optimization.

In contrast to pixelated SMO, we may note that the source and mask always maintain an exceedingly low complexity after optimization. The source has been optimized by adjusting the partial coherent factor, while the mask has been optimized by changing the feature size and adding a serif. There is no need to customize the diffractive optical elements or purchase an expensive pixel source generator when using the parametric source. The mask can maintain a low complexity after optimization. The parametric source and mask can, therefore, effectively reduce the cost of fabrication and provide high stability in production.

Figure 8 presents the convergent curves of the cost functions for SMO and SMNO. The red-circle line represents the cost function of SMO, while the black-triangle line illustrates SMNO. By incorporating NA in optimization, the cost function could be reduced much further.

The above performance comparisons reveal that the proposed parametric SMNO method can effectively improve the pattern fidelity and enhance the robustness of the optical lithography systems. NA should be considered in the optimization when optimizing the source and mask.

## 6.

## Conclusion

This paper proposes a parametric SMNO method using a vector imaging model to improve the pattern fidelity and DOF. We develop an analytical approach to the parametric source, mask, and NA, which could be effectively applied in inverse lithography. The mathematical expressions of the derivatives of source, mask, and NA parameters are derived on the basis of the vector imaging model. The steepest-descent algorithm is used to optimize the source, mask, and NA iteratively. Simulation results show that the SMNO produces a better source, mask, and NA combination compared to SMO. Simulation results also reveal that it is necessary to optimize NA to achieve superior lithography performance. The optimized parametric source, mask, and NA are capable of maintaining high pattern fidelity within a large DOF. The parametric source shape and mask layout are in an extremely simple distribution, which could effectively reduce the cost of fabrication and provide high stability in high-volume production. In our future work, we will improve the parametric expression of the mask to model a more complex shape, such as an L-shape and T-shape, and apply the SMNO to full-chip mask optimization.

## Appendices

## Appendix A

According to Eq. (10), the derivatives of the source parameters $\partial {S}_{Q}/\partial \mathrm{\Omega}({s}_{v})$ are

## (52)

$$\frac{\partial S}{\partial \mathrm{\Omega}({\sigma}_{\text{out}})}={K}_{R}\xb7(\frac{2{b}_{s}\xb7{({r}_{s}/{\sigma}_{\text{out}})}^{2{b}_{s}-1}}{{[{({r}_{s}/{\sigma}_{\text{out}})}^{2{b}_{s}}+1]}^{2}}\xb7\frac{{r}_{s}}{{({\sigma}_{\text{out}})}^{2}}\phantom{\rule{0ex}{0ex}}-\frac{2{b}_{s}\xb7{[{r}_{s}/({\sigma}_{\text{out}}-{\sigma}_{\text{width}})]}^{2{b}_{s}-1}}{{\{{[{r}_{s}/({\sigma}_{\text{out}}-{\sigma}_{\text{width}})]}^{2{b}_{s}}+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{{r}_{s}}{{({\sigma}_{\text{out}}-{\sigma}_{\text{width}})}^{2}})\xb7\frac{\partial {\sigma}_{\text{out}}}{\partial \mathrm{\Omega}({\sigma}_{\text{out}})},$$## (53)

$$\frac{\partial S}{\partial \mathrm{\Omega}({\sigma}_{\text{width}})}={K}_{R}\xb7(\frac{2{b}_{s}\xb7{[{r}_{s}/({\sigma}_{\text{out}}-{\sigma}_{\text{width}})]}^{2{b}_{s}-1}}{{\{{[{r}_{s}/({\sigma}_{\text{out}}-{\sigma}_{\text{width}})]}^{2{b}_{s}}+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{{r}_{s}}{{({\sigma}_{\text{out}}-{\sigma}_{\text{width}})}^{2}})\xb7\frac{\partial {\sigma}_{\text{width}}}{\partial \mathrm{\Omega}({\sigma}_{\text{width}})},$$## (54)

$$\frac{\partial S}{\partial \mathrm{\Omega}({\theta}_{1})}\phantom{\rule{0ex}{0ex}}={S}_{A}\xb7(\frac{2{b}_{s}\xb7{\left[\right(2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180)/(180-{\theta}_{1})]}^{2{b}_{s}-1}}{{\{{\left[\right(2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180)/(180-{\theta}_{1})]}^{2{b}_{s}}+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180}{{(180-{\theta}_{1})}^{2}})\xb7\frac{\partial {\theta}_{1}}{\partial \mathrm{\Omega}({\theta}_{1})},$$## (55)

$$\frac{\partial S}{\partial \mathrm{\Omega}({\theta}_{2})}\phantom{\rule{0ex}{0ex}}={S}_{A}\xb7(\frac{2{b}_{s}\xb7{\left[\right(2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180)/{\theta}_{2}]}^{2{b}_{s}-1}}{{\{{\left[\right(2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180)/{\theta}_{2}]}^{2{b}_{s}}+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{2\xb7\frac{180}{\pi}\xb7a\text{\hspace{0.17em}}\mathrm{cos}\frac{{\alpha}_{s}}{{r}_{s}}-180}{{({\theta}_{2})}^{2}})\xb7\frac{\partial {\theta}_{2}}{\partial \mathrm{\Omega}({\theta}_{2})}.$$## Appendix B

The derivatives of the mask with respect to the mask parameter $\partial {m}_{V}/\partial \mathrm{\Omega}({m}_{z})$ are

## (56)

$$\frac{\partial {m}_{F}^{i}}{\partial \mathrm{\Omega}({m}_{w}^{i})}=-\frac{2{b}_{m}{\left[\frac{2(x-{x}_{0}^{i})}{{m}_{w}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(x-{x}_{0}^{i})}{{m}_{w}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\xb7\frac{2(x-{x}_{0}^{i})}{{({m}_{w}^{i})}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{1}{{\left[\frac{2(y-{y}_{0}^{i})}{{m}_{w}^{i}}\right]}^{2{b}_{m}}+1}\xb7\frac{\partial {m}_{w}^{i}}{\partial \mathrm{\Omega}({m}_{w}^{i})},$$## (57)

$$\frac{\partial {m}_{F}^{i}}{\partial \mathrm{\Omega}({m}_{h}^{i})}=-\frac{2{b}_{m}{\left[\frac{2(y-{y}_{0}^{i})}{{m}_{h}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(y-{y}_{0}^{i})}{{m}_{h}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\xb7\frac{2(y-{y}_{0}^{i})}{{({m}_{h}^{i})}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{1}{{\left[\frac{2(x-{x}_{0}^{i})}{{m}_{h}^{i}}\right]}^{2{b}_{m}}+1}\xb7\frac{\partial {m}_{h}^{i}}{\partial \mathrm{\Omega}({m}_{h}^{i})},$$## (58)

$$\frac{\partial {m}_{S}^{i}}{\partial \mathrm{\Omega}({m}_{\mathrm{ss}}^{i})}=-(\frac{2{b}_{m}{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\xb7\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{({m}_{\mathrm{ss}}^{i})}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{1}{{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1}+\frac{2{b}_{m}{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\phantom{\rule{0ex}{0ex}}\xb7\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{({m}_{\mathrm{ss}}^{i})}^{2}}\xb7\frac{1}{{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1})\xb7\frac{\partial {m}_{\mathrm{ss}}^{i}}{\partial \mathrm{\Omega}({m}_{\mathrm{ss}}^{i})},\phantom{\rule{0ex}{0ex}}$$## (59)

$$\frac{\partial {m}_{S}^{i}}{\partial \mathrm{\Omega}({m}_{\mathrm{so}}^{i})}\phantom{\rule{0ex}{0ex}}=-(\frac{2{b}_{m}{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\xb7\frac{-2}{{m}_{\mathrm{ss}}^{i}}\xb7\frac{1}{{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1}\phantom{\rule{0ex}{0ex}}+\frac{2{b}_{m}{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{M}_{\mathrm{SS}}^{i}}\right]}^{2{b}_{m}-1}}{{\{{\left[\frac{2(y-{y}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1\}}^{2}}\xb7\frac{-2}{{m}_{\mathrm{ss}}^{i}}\xb7\frac{1}{{\left[\frac{2(x-{x}_{0}^{i}-{m}_{\mathrm{so}}^{i})}{{m}_{\mathrm{ss}}^{i}}\right]}^{2{b}_{m}}+1})\phantom{\rule{0ex}{0ex}}\xb7\frac{\partial {m}_{\mathrm{so}}^{i}}{\partial \mathrm{\Omega}({m}_{\mathrm{so}}^{i})}.$$## Acknowledgments

We thank the financial support by National Science and Technology Major Project, National Natural Science Foundation of China under No. 61204113, Key Laboratory of Photoelectronic Imaging Technology and System (Beijing Institute of Technology, Ministry of Education of China) under Grant No. 2013OEIOF05, and the project of Mentor Graphics and Beijing Institute of Technology Collaboration in RET. We also thank KLA-Tencor Corporation for providing academic use of PROLITH.

## References

## Biography

**Xuejia Guo** received his BS degree in optical information science and technology from University of Jinan in 2008 and his PhD degree in optical engineering from Beijing Institute of Technology in 2014. His current interests include optical manufacturing, lithography simulation, and optimization.

**Yanqiu Li** received her MS and PhD degrees in optics from Harbin Institute of Technology of China. She is currently a professor in the School of Optoelectronics at Beijing Institute of Technology, Beijing, China. She holds over 30 Chinese patents and has published numerous articles on lithographic science.

**Lisong Dong** received his BS degree in applied physics from Hefei University of Technology in 2008 and his PhD degree in optical engineering from Beijing Institute of Technology in 2014. His research interests include lithography simulation and resolution enhancement technology.

**Lihui Liu** received his BS and MS degrees in optics from Hebei University, respectively, in 1999 and 2002. He received his PhD from the Institute of Modern Optics in NanKai University in 2005. He is currently an associate professor of School of Optoelectronics at Beijing Institute of Technology, Beijing, China. His research is focused on optical sensing, lithography simulation, and optimization.

**Xu Ma** received his BS degree in electrical engineering from Tsinghua University, Beijing, China, in 2005, and MS and PhD degrees in electrical and computer engineering from the University of Delaware in 2007 and 2009, respectively. Until 2010, he was a postdoctoral scholar at the University of California, Berkeley. He is currently a professor in the School of Optoelectronics at Beijing Institute of Technology, China. His research interests include computational lithography and signal and image processing.

**Chunying Han** received his BS degree in physics from Changsha University of Science & Technology, China, in 2009. He is currently a PhD candidate directed by Professor Yanqiu Li at Beijing Institute of Technology, Beijing, China. His research interests include lithography simulation and computational lithography.