Horizontal development bias in negative tone development of photoresist

Abstract

Aspects of the disclosed techniques relate to techniques for resist simulation in lithography. Local minimal light intensity values are determined for a plurality of sample points in boundary regions of an aerial image of a feature to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points. Based on the local minimal light intensity values, horizontal development bias values for the plurality of sample points are then determined. Finally, resist contour data of the feature are determined based at least on the horizontal development bias values.

Claims

1. A method, executed by at least one processor of a computer, comprising: determining local minimal light intensity values for a plurality of sample points in boundary regions of an aerial image of a feature of an electronic design to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points; determining horizontal development bias values for the plurality of sample points based on the local minimal light intensity values; and computing, as part of a process for verifying printability of the feature or for modifying the feature, resist contour data of the feature based at least on the horizontal development bias values.

2. The method recited in claim 1, wherein the local minimal light intensity values are determined according to $I_m\left(x,y\right)=\frac{1}{2*b}*\left(b^2-\left({\left(I_{-b}\right)}^2.Math.G_s\right)\left(x,y\right)\right),$ where $G_s\left(x,y\right)=\frac{1}{s^2}\exp \left(-\frac{x^2+y^2}{s^2}\right)$ is a Gaussian kernel with diffusion length s, I.sub.b(x, y) is a result of the minus neutralization with the neutralization constant b, and x, y are coordinates, and wherein the horizontal development bias values are determined according to $H\left(x,y\right)=c*\sqrt{a^2+{\left(I_m\left(x,y\right)\right)}^2}*\sqrt{{\left(\frac{I\left(x,y\right)}{x}\right)}^2+{\left(\frac{I\left(x,y\right)}{y}\right)}^2}.$ where c is a linear coefficient, a is a constant, and I(x,y) is light intensity at (x,y).

3. The method recited in claim 1, wherein the horizontal development bias values are determined according to
xRmin*[t.sub.devh where Imin represents the local minimal light intensity values, R.sub.max is dissolution rate of fully exposed resist, R.sub.min is dissolution rate of unexposed resist, t.sub.dev is horizontal development time, h is depth, and and r.sub.0 are constants.

4. A non-transitory processor-readable medium storing processor-executable instructions for causing one or more processors to perform a method, the method comprising: determining local minimal light intensity values for a plurality of sample points in boundary regions of an aerial image of a feature of an electronic design to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points; determining horizontal development bias values for the plurality of sample points based on the local minimal light intensity values; and computing, as part of a process for verifying printability of the feature or for modifying the feature, resist contour data of the feature based at least on the horizontal development bias values.

5. The non-transitory processor-readable medium recited in claim 4, wherein the local minimal light intensity values are determined according to $I_m\left(x,y\right)=\frac{1}{2*b}*\left(b^2-\left({\left(I_{-b}\right)}^2.Math.G_s\right)\left(x,y\right)\right),$ where $G_s\left(x,y\right)=\frac{1}{s^2}\exp \left(-\frac{x^2+y^2}{s^2}\right)$ is a Gaussian kernel with diffusion length s, I.sub.b(x, y) is a result of the minus neutralization with the neutralization constant b, and x, y are coordinates, and wherein the horizontal development bias values are determined according to $H\left(x,y\right)=c*\sqrt{a^2+{\left(I_m\left(x,y\right)\right)}^2}*\sqrt{{\left(\frac{I\left(x,y\right)}{x}\right)}^2+{\left(\frac{I\left(x,y\right)}{y}\right)}^2}$ where c is a linear coefficient, a is a constant, and I(x,y) is light intensity at (x,y).

6. The non-transitory processor-readable medium recited in claim 4, wherein the horizontal development bias values are determined according to
xRmin*[t.sub.devh where Imin represents the local minimal light intensity values, R.sub.max is dissolution rate of fully exposed resist, R.sub.min is dissolution rate of unexposed resist, t.sub.dev is horizontal development time, h is depth, and and r.sub.0 are constants.

7. A system, comprising: a local minimum light intensity determination unit configured to determine local minimal light intensity values for a plurality of sample points in boundary regions of an aerial image of a feature of an electronic design to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points; a horizontal development bias values determination unit configured to determine horizontal development bias values for the plurality of sample points based on the local minimal light intensity values; and a resist contour simulation unit configured to compute, as part of a process for verifying printability of the feature or for modifying the feature, resist contour data of the feature based at least on the horizontal development bias values.

8. The system recited in claim 7, wherein the local minimal light intensity values are determined according to $I_m\left(x,y\right)=\frac{1}{2*b}*\left(b^2-\left({\left(I_{-b}\right)}^2.Math.G_s\right)\left(x,y\right)\right),$ where $G_s\left(x,y\right)=\frac{1}{s^2}\exp \left(-\frac{x^2+y^2}{s^2}\right)$ is a Gaussian kernel with diffusion length s, I.sub.b(x, y) is a result of the minus neutralization with the neutralization constant b, and x, y are coordinates, and wherein the horizontal development bias values are determined according to $H\left(x,y\right)=c*\sqrt{a^2+{\left(I_m\left(x,y\right)\right)}^2}*\sqrt{{\left(\frac{I\left(x,y\right)}{x}\right)}^2+{\left(\frac{I\left(x,y\right)}{y}\right)}^2}$ where c is a linear coefficient, a is a constant, and I(x,y) is light intensity at (x,y).

9. The system recited in claim 7, wherein the horizontal development bias values are determined according to
xRmin*[t.sub.devh where Imin represents the local minimal light intensity values, R.sub.max is dissolution rate of fully exposed resist, R.sub.min is dissolution rate of unexposed resist, t.sub.dev is horizontal development time, h is depth, and and r.sub.0 are constants.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 illustrates an example of a computing system that may be used to implement various embodiments of the disclosed techniques.

(2) FIG. 2 illustrates an example of a multi-core processor unit that may be used to implement various embodiments of the disclosed techniques.

(3) FIG. 3 illustrates an example of a chemically amplified photoresist comprising a poly(4-hydroxystyrene) (pHOST) protected by an acid-labile tert-butoxycarbonyl (t-BOC) functional group and a triphenylsulfonium hexafluoroantimonate (TPS.SbF.sub.6).

(4) FIG. 4A illustrates an example of a photochemical reaction of a photo acid generator.

(5) FIG. 4B illustrates an example of an acid-catalyzed deprotection reaction.

(6) FIG. 5 illustrates one example for small r.sub.min resulting in a negative slope of resist profile and another example for large r.sub.min resulting in a positive slope of resist profile.

(7) FIG. 6 illustrates an example of a resist simulation tool that may be employed according to various embodiments of the disclosed techniques

(8) FIG. 7 illustrates a flowchart describing a process for resist simulation that may be employed by various embodiments of the disclosed techniques.

DETAILED DESCRIPTION OF THE DISCLOSED TECHNIQUES

General Considerations

(9) Various aspects of the present disclosed techniques relate to techniques for resist simulation in lithography. In the following description, numerous details are set forth for the purpose of explanation. However, one of ordinary skill in the art will realize that the disclosed techniques may be practiced without the use of these specific details. In other instances, well-known features have not been described in details to avoid obscuring the present disclosed techniques.

(10) Some of the techniques described herein can be implemented in software instructions stored on a computer-readable medium, software instructions executed on a computer, or some combination of both. Some of the disclosed techniques, for example, can be implemented as part of an electronic design automation (EDA) tool. Such methods can be executed on a single computer or on networked computers.

(11) Although the operations of the disclosed methods are described in a particular sequential order for convenient presentation, it should be understood that this manner of description encompasses rearrangements, unless a particular ordering is required by specific language set forth below. For example, operations described sequentially may in some cases be rearranged or performed concurrently. Moreover, for the sake of simplicity, the disclosed flow charts and block diagrams typically do not show the various ways in which particular methods can be used in conjunction with other methods. Additionally, the detailed description sometimes uses terms like determine and compute to describe the disclosed methods. Such terms are high-level abstractions of the actual operations that are performed. The actual operations that correspond to these terms will vary depending on the particular implementation and are readily discernible by one of ordinary skill in the art.

(12) Also, as used herein, the term design is intended to encompass data describing an entire integrated circuit device. This term also is intended to encompass a smaller group of data describing one or more components of an entire device, however, such as a portion of an integrated circuit device. Still further, the term design also is intended to encompass data describing more than one microdevice, such as data to be used to form multiple microdevices on a single wafer.

(13) Illustrative Operating Environment

(14) The execution of various electronic design automation processes according to embodiments of the disclosed techniques may be implemented using computer-executable software instructions executed by one or more programmable computing devices. Because these embodiments of the disclosed techniques may be implemented using software instructions, the components and operation of a generic programmable computer system on which various embodiments of the disclosed techniques may be employed will first be described. Further, because of the complexity of some electronic design automation processes and the large size of many circuit designs, various electronic design automation tools are configured to operate on a computing system capable of concurrently running multiple processing threads. The components and operation of a computer network having a host or master computer and one or more remote or servant computers therefore will be described with reference to FIG. 1. This operating environment is only one example of a suitable operating environment, however, and is not intended to suggest any limitation as to the scope of use or functionality of the disclosed techniques.

(15) In FIG. 1, the computer network 101 includes a master computer 103. In the illustrated example, the master computer 103 is a multi-processor computer that includes a plurality of input and output devices 105 and a memory 107. The input and output devices 105 may include any device for receiving input data from or providing output data to a user. The input devices may include, for example, a keyboard, microphone, scanner or pointing device for receiving input from a user. The output devices may then include a display monitor, speaker, printer or tactile feedback device. These devices and their connections are well known in the art, and thus will not be discussed at length here.

(16) The memory 107 may similarly be implemented using any combination of computer readable media that can be accessed by the master computer 103. The computer readable media may include, for example, microcircuit memory devices such as read-write memory (RAM), read-only memory (ROM), electronically erasable and programmable read-only memory (EEPROM) or flash memory microcircuit devices, CD-ROM disks, digital video disks (DVD), or other optical storage devices. The computer readable media may also include magnetic cassettes, magnetic tapes, magnetic disks or other magnetic storage devices, punched media, holographic storage devices, or any other medium that can be used to store desired information.

(17) As will be discussed in detail below, the master computer 103 runs a software application for performing one or more operations according to various examples of the disclosed techniques. Accordingly, the memory 107 stores software instructions 109A that, when executed, will implement a software application for performing one or more operations. The memory 107 also stores data 109B to be used with the software application. In the illustrated embodiment, the data 109B contains process data that the software application uses to perform the operations, at least some of which may be parallel.

(18) The master computer 103 also includes a plurality of processor units 111 and an interface device 113. The processor units 111 may be any type of processor device that can be programmed to execute the software instructions 109A, but will conventionally be a microprocessor device. For example, one or more of the processor units 111 may be a commercially generic programmable microprocessor, such as Intel Pentium or Xeon microprocessors, Advanced Micro Devices Athlon microprocessors or Motorola 68K/Coldfire microprocessors. Alternately or additionally, one or more of the processor units 111 may be a custom-manufactured processor, such as a microprocessor designed to optimally perform specific types of mathematical operations. The interface device 113, the processor units 111, the memory 107 and the input/output devices 105 are connected together by a bus 115.

(19) With some implementations of the disclosed techniques, the master computing device 103 may employ one or more processing units 111 having more than one processor core. Accordingly, FIG. 2 illustrates an example of a multi-core processor unit 111 that may be employed with various embodiments of the disclosed techniques. As seen in this figure, the processor unit 111 includes a plurality of processor cores 201. Each processor core 201 includes a computing engine 203 and a memory cache 205. As known to those of ordinary skill in the art, a computing engine contains logic devices for performing various computing functions, such as fetching software instructions and then performing the actions specified in the fetched instructions. These actions may include, for example, adding, subtracting, multiplying, and comparing numbers, performing logical operations such as AND, OR, NOR and XOR, and retrieving data. Each computing engine 203 may then use its corresponding memory cache 205 to quickly store and retrieve data and/or instructions for execution.

(20) Each processor core 201 is connected to an interconnect 207. The particular construction of the interconnect 207 may vary depending upon the architecture of the processor unit 111. With some processor cores 201, such as the Cell microprocessor created by Sony Corporation, Toshiba Corporation and IBM Corporation, the interconnect 207 may be implemented as an interconnect bus. With other processor units 111, however, such as the Opteron and Athlon dual-core processors available from Advanced Micro Devices of Sunnyvale, Calif., the interconnect 207 may be implemented as a system request interface device. In any case, the processor cores 201 communicate through the interconnect 207 with an input/output interface 209 and a memory controller 210. The input/output interface 209 provides a communication interface between the processor unit 111 and the bus 115. Similarly, the memory controller 210 controls the exchange of information between the processor unit 111 and the system memory 107. With some implementations of the disclosed techniques, the processor units 111 may include additional components, such as a high-level cache memory accessible shared by the processor cores 201.

(21) While FIG. 2 shows one illustration of a processor unit 111 that may be employed by some embodiments of the disclosed techniques, it should be appreciated that this illustration is representative only, and is not intended to be limiting. Also, with some implementations, a multi-core processor unit 111 can be used in lieu of multiple, separate processor units 111. For example, rather than employing six separate processor units 111, an alternate implementation of the disclosed techniques may employ a single processor unit 111 having six cores, two multi-core processor units each having three cores, a multi-core processor unit 111 with four cores together with two separate single-core processor units 111, etc.

(22) Returning now to FIG. 1, the interface device 113 allows the master computer 103 to communicate with the servant computers 117A, 117B, 117C . . . 117x through a communication interface. The communication interface may be any suitable type of interface including, for example, a conventional wired network connection or an optically transmissive wired network connection. The communication interface may also be a wireless connection, such as a wireless optical connection, a radio frequency connection, an infrared connection, or even an acoustic connection. The interface device 113 translates data and control signals from the master computer 103 and each of the servant computers 117 into network messages according to one or more communication protocols, such as the transmission control protocol (TCP), the user datagram protocol (UDP), and the Internet protocol (IP). These and other conventional communication protocols are well known in the art, and thus will not be discussed here in more detail.

(23) Each servant computer 117 may include a memory 119, a processor unit 121, an interface device 123, and, optionally, one more input/output devices 125 connected together by a system bus 127. As with the master computer 103, the optional input/output devices 125 for the servant computers 117 may include any conventional input or output devices, such as keyboards, pointing devices, microphones, display monitors, speakers, and printers. Similarly, the processor units 121 may be any type of conventional or custom-manufactured programmable processor device. For example, one or more of the processor units 121 may be commercially generic programmable microprocessors, such as Intel Pentium or Xeon microprocessors, Advanced Micro Devices Athlon microprocessors or Motorola 68K/Coldfire microprocessors. Alternately, one or more of the processor units 121 may be custom-manufactured processors, such as microprocessors designed to optimally perform specific types of mathematical operations. Still further, one or more of the processor units 121 may have more than one core, as described with reference to FIG. 2 above. For example, with some implementations of the disclosed techniques, one or more of the processor units 121 may be a Cell processor. The memory 119 then may be implemented using any combination of the computer readable media discussed above. Like the interface device 113, the interface devices 123 allow the servant computers 117 to communicate with the master computer 103 over the communication interface.

(24) In the illustrated example, the master computer 103 is a multi-processor unit computer with multiple processor units 111, while each servant computer 117 has a single processor unit 121. It should be noted, however, that alternate implementations of the disclosed techniques may employ a master computer having single processor unit 111. Further, one or more of the servant computers 117 may have multiple processor units 121, depending upon their intended use, as previously discussed. Also, while only a single interface device 113 or 123 is illustrated for both the master computer 103 and the servant computers, it should be noted that, with alternate embodiments of the disclosed techniques, either the computer 103, one or more of the servant computers 117, or some combination of both may use two or more different interface devices 113 or 123 for communicating over multiple communication interfaces.

(25) With various examples of the disclosed techniques, the master computer 103 may be connected to one or more external data storage devices. These external data storage devices may be implemented using any combination of computer readable media that can be accessed by the master computer 103. The computer readable media may include, for example, microcircuit memory devices such as read-write memory (RAM), read-only memory (ROM), electronically erasable and programmable read-only memory (EEPROM) or flash memory microcircuit devices, CD-ROM disks, digital video disks (DVD), or other optical storage devices. The computer readable media may also include magnetic cassettes, magnetic tapes, magnetic disks or other magnetic storage devices, punched media, holographic storage devices, or any other medium that can be used to store desired information. According to some implementations of the disclosed techniques, one or more of the servant computers 117 may alternately or additionally be connected to one or more external data storage devices. Typically, these external data storage devices will include data storage devices that also are connected to the master computer 103, but they also may be different from any data storage devices accessible by the master computer 103.

(26) It also should be appreciated that the description of the computer network illustrated in FIG. 1 and FIG. 2 is provided as an example only, and it not intended to suggest any limitation as to the scope of use or functionality of alternate embodiments of the disclosed techniques.

(27) Development Rates and Horizontal Development Bias

(28) There are a variety of photoresists in the market, with properties optimized for various exposure and processing conditions. Most photoresists in use today are positive resists, meaning that exposure to light causes bonds in the resist to break. This scission of bonds reduces local molecular weight and increases solubility. Portions exposed to light therefore wash away when developed, leaving only the unexposed portions on the wafer. This is referred to as positive tone development. In a negative tone development process, an organic solvent is used to wash away unexposed regions of the resist film. There are also negative resists, which cross-link and form less soluble regions when exposed to light, but these are not as common in modern lithographic processes.

(29) As noted previously, chemically amplified resists have been developed to increase sensitivity for technology nodes in deep UV. Upon exposure, the photo-acid generator (PAG) creates an acid that diffuse through the resist, removing protecting groups from the polymer chains. This deprotection reaction is a catalyzed reaction. The exposure to a single exposing photon, which could only break a single bond in the past, can break as many bonds as the catalysis (i.e., the acid) can contact. The extent and number of bonds that can be broken by the catalyst varies from resist to resist, and the distance over which the acid can diffuse to break bonds can be highly dependent on local development conditions, such as temperature and the base quencher concentration.

(30) A resist model usually includes modeling the development step. The negative tone development may be modeled as:

(31) $\begin{array}{cc}r=r_{\max }\frac{\left(a+1\right){\left(1-m\right)}^n}{a+{\left(1-m\right)}^n}+r_{\min }a=r_{\max }\frac{\left(n+1\right)}{\left(n-1\right)}{\left(1-m_{\mathrm{TH}}\right)}^{n}m=\frac{C_{\mathrm{UnBlockedPolymer}}}{C_{\mathrm{Blocked}\mathrm{Polymer}}^{\mathrm{Initial}}}& \left(1\right)\end{array}$
where m.sub.TH is the threshold inhibitor concentration, r.sub.max is the dissolution rate of unexposed resist, r.sub.min is the dissolution rate of fully exposed resist, and n is the dissolution selectivity parameter controlling the contrast of the resist. For the positive tone development, Eq. (1) still applies except:

(32) $m=\frac{C_{\mathrm{BlockedPolymer}}}{C_{\mathrm{Blocked}\mathrm{Polymer}}^{\mathrm{Initial}}}$
and r.sub.max is the dissolution rate of fully exposed resist, r.sub.min is the dissolution rate of unexposed resist.

(33) Compared to the positive tone development, the negative tone development has a r.sub.min large enough to have a non-negligible effect on the developed resist profile. The resist development path may be considered to include two parts approximately: first a vertical path and then a horizontal path. The effect of the horizontal path is herein referred to as horizontal development bias (or additional horizontal development bias). The r.sub.min difference may lead to different slope of resist profile. FIG. 5 illustrates an example (510) for small r.sub.min and an example (520) for large r.sub.min. The horizontal development bias for 520 is significantly larger than that for 510. As a result, the resist profile of 510 has a negative slope while the resist profile of 520 have a positive slope.

(34) As shown in FIG. 5, the horizontal development bias is z-dependent or depth-dependent. The z-dependent bias relates to horizontal development time:
t.sub.dev=t.sub.z+t.sub.x(2)
where

(35) $t_z={}_0^h\frac{z}{r\left(x,z\right)}{}_0^h\frac{z}{\min \left(R\max ,{r_0\left[I\left(x0\right)\right]}^{}\right)}=\frac{h}{\min \left(R\max ,{r_0\left[I\left(x0\right)\right]}^{}\right)}$$t_x=t_{\mathrm{dev}}$

(36) For the negative tone development, assume at the horizontal development path: r(x, z)=Rmin+r.sub.0[I(x, z)].sup.. Then the corresponding horizontal development bias is:

(37) $\begin{array}{cc}\begin{array}{c}x=\mathrm{Rmin}*t_x\mathrm{Rmin}*\left[t_{\mathrm{dev}}-\frac{h}{\min \left(\mathrm{Rmax},{r_0\left[I\left(x0\right)\right]}^{}\right)}\right]\\ \mathrm{Rmin}*\left[t_{\mathrm{dev}}-h*\max \left(1/\mathrm{Rmax},{\left[I\left(x0\right)\right]}^{.Math..Math.}/r_0\right)\right]\\ \mathrm{Rmin}*\left[t_{\mathrm{dev}}-h*\max \left(1/\mathrm{Rmax},{\mathrm{Imin}}^{.Math..Math.}/r_0\right)\right]\end{array}& \left(3\right)\end{array}$

(38) The horizontal development bias may be incorporated into the CM1 resist model discussed in the background section by adding a horizontal development bias term into the resist surface function. The horizontal development bias term may be expressed as:

(39) $\begin{array}{cc}H\left(x,y\right)=c*\sqrt{a^2+{\left(I_m\left(x,y\right)\right)}^2}*\sqrt{{\left(\frac{I\left(x,y\right)}{x}\right)}^2+{\left(\frac{I\left(x,y\right)}{y}\right)}^2}& \left(4\right)\end{array}$
where c is a linear coefficient, a is a constant, and the approximate minimal aerial image intensity is defined as

(40) $\begin{array}{cc}{I}_m\left(x,y\right)=\frac{1}{2*b}*\left(b^2-\left({\left(I_{-b}\right)}^2.Math.G_s\right)\left(x,y\right)\right)& \left(5\right)\end{array}$
where I.sub.b(x, y) is the result of the minus neutralization with the neutralization constant b, and

(41) 0$\begin{array}{cc}{G}_s\left(x,y\right)=\frac{1}{s^2}\exp \left(-\frac{x^2+y^2}{s^2}\right)& \left(6\right)\end{array}$
is a Gaussian kernel with diffusion length s. To see why I.sub.m(x,y) approximately represents the minimal aerial image within the distance s from the observation point (x; y), assume that the result of the convolution:
((I.sub.b).sup.2G.sub.s)(x,y)=(bI.sub.min(x,y)).sup.2(7)

(42) (this would be correct if the calibration with the Gaussian kernel exactly picked up the maximum within its diffusion length). Then, I.sub.m(x, y)=I.sub.min(x, y) assuming that the constant b is large compared with the minimal aerial image I.sub.min(x, y).

(43) The square root {square root over (a.sup.2+(I.sub.m(x, y)).sup.2)} is a smooth approximation to the maximum function max{a, I.sub.m(x, y)}. Thus, the total horizontal bias term is an approximation to the quantity H.sub.b(x, y)=|I(x,y)|*max{a, I.sub.min(x, y)}.

(44) Resist Simulation Tools and Methods

(45) FIG. 6 illustrates an example of a resist simulation tool according to various embodiments of the disclosed techniques. As seen in the figure, the resist simulation tool 600 includes three units: a local minimum light intensity determination unit 620, a horizontal development bias values determination unit 640 and a resist contour simulation unit 660. Some implementations of the resist simulation tool 600 may cooperate with (or incorporate) one or more of an input database 615 and an output database 695.

(46) As will be discussed in more detail below, the local minimum light intensity determination unit 620 determines local minimal light intensity values for a plurality of sample points in boundary regions of an aerial image of a feature to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points. The aerial image may be obtained by an optical simulator and stored in the input database 615. Based on the local minimal light intensity values, the horizontal development bias values determination unit 640 determines horizontal development bias values for the plurality of sample points. Based at least on the horizontal development bias values, the resist contour simulation unit 660 simulates resist contour data of the feature, which may be stored in the output database 695.

(47) As previously noted, various examples of the disclosed techniques may be implemented by a multiprocessor computing system, such as the computing system illustrated in FIGS. 1 and 2. Accordingly, one or more of the local minimum light intensity determination unit 620, the horizontal development bias values determination unit 640 and the resist contour simulation unit 660 may be implemented by executing programming instructions on one or more processors in a computing system such as the computing system illustrated in FIG. 1 and FIG. 2. Correspondingly, some other embodiments of the disclosed techniques may be implemented by software instructions, stored on a non-transitory computer-readable medium, for instructing one or more programmable computers/computer systems to perform the functions of one or more of the local minimum light intensity determination unit 620, the horizontal development bias values determination unit 640 and the resist contour simulation unit 660. As used herein, the term non-transitory computer-readable medium refers to computer-readable medium that are capable of storing data for future retrieval, and not propagating electro-magnetic waves. The non-transitory computer-readable medium may be, for example, a magnetic storage device, an optical storage device, a punched surface type device, or a solid state storage device.

(48) It also should be appreciated that, while the local minimum light intensity determination unit 620, the horizontal development bias values determination unit 640 and the resist contour simulation unit 660 are shown as separate units in FIG. 6, a single servant computer (or a single processor within a master computer) may be used to implement two or more of these units at different times, or components of two or more of these units at different times.

(49) With various examples of the disclosed techniques, the input database 615 and the output database 695 may be implemented using any suitable computer readable storage device. That is, either of the input database 615 and the output database 695 may be implemented using any combination of computer readable storage devices including, for example, microcircuit memory devices such as read-write memory (RAM), read-only memory (ROM), electronically erasable and programmable read-only memory (EEPROM) or flash memory microcircuit devices, CD-ROM disks, digital video disks (DVD), or other optical storage devices. The computer readable storage devices may also include magnetic cassettes, magnetic tapes, magnetic disks or other magnetic storage devices, punched media, holographic storage devices, or any other non-transitory storage medium that can be used to store desired information. While the input database 615 and the output database 695 are shown as separate units in FIG. 6, a single data storage medium may be used to implement some or all of these databases.

(50) FIG. 7 illustrates a flowchart showing a process for resist simulation that may be implemented according to various examples of the disclosed techniques. For ease of understanding, methods of resist simulation that may be employed according to various embodiments of the disclosed techniques will be described with reference to the resist simulation tool 600 illustrated in FIG. 6 and the flow chart 700 in FIG. 7. It should be appreciated, however, that alternate implementations of an resist simulation tool may be used to perform the method of resist simulation shown in the flow chart 700 according to various embodiments of the disclosed techniques. In addition, it should be appreciated that implementations of the resist simulation tool 600 may be employed to implement methods of resist simulation according to different embodiments of the disclosed techniques other than the one illustrated by the flow chart 700 in FIG. 7.

(51) In operation 720 of the flow chart 700, the local minimum light intensity determination unit 620 determines local minimal light intensity values for a plurality of sample points in boundary regions of an aerial image of a feature to be printed on a resist coating, wherein each of the local minimal light intensity values represents a minimum light intensity value for an area surrounding one of the plurality of sample points. The aerial image of the feature may be derived by an optical simulator, such as those found in the CALIBRE family of software tools available from Mentor Graphics Corporation, Wilsonville, Oreg. As noted previously, the disclosed techniques may be incorporated into the CM1 resist model. The plurality of sample points may be some of the grid points employed by the CM1 resist model that are in the boundary regions of the aerial image. The area may have a circle shape. The circle's radius may be adjusted when the resist model is calibrated. With some implementations of the disclosed techniques in which the resist simulation tool 600 employs, for example, the CM1 resist model, the local minimal light intensity values may be determined according to Eq. (5).

(52) In operation 740, the horizontal development bias values determination unit 640 determines horizontal development bias values for the plurality of sample points based on the local minimal light intensity values. Some embodiments of the disclosed techniques use Eq. (3) to compute the horizontal development bias values. Some other embodiments of the disclosed techniques use Eq. (4) for the operation 740.

(53) In operation 760, the resist contour simulation unit 660 computes resist contour data of the feature based at least on the horizontal development bias values. When the CM1 resist model is employed by the resist simulation tool 600, the horizontal development bias term defined by Eq. (4) is added into the resist surface function. The resist surface function may be calibrated against critical dimension data measured or calculated using a rigorous model. By using a threshold constant, the resist contour data may be derived from the resist surface function.

CONCLUSION

(54) While the disclosed techniques has been described with respect to specific examples including presently preferred modes of carrying out the disclosed techniques, those skilled in the art will appreciate that there are numerous variations and permutations of the above described systems and techniques that fall within the spirit and scope of the disclosed techniques as set forth in the appended claims. For example, while specific terminology has been employed above to refer to electronic design automation processes, it should be appreciated that various examples of the disclosed techniques may be implemented using any desired combination of electronic design automation processes.