Abstract:
Techniques for model calibration and alignment of measurement contours of printed layout features with simulation contours obtained with a model are disclosed. With various implementations of the invention, contour point errors are determined. Based on the contour point errors and a cost function, values of alignment parameters may be determined. The values of alignment parameters may be used to realign the measurement contours for model calibration. The alignment may be conducted concurrently with model calibration.

Description:
RELATED APPLICATIONS 
       [0001]    This application claims priority to U.S. Provisional Patent Application No. 61/304,807, entitled “Contour Self-Alignment For OPC Model Calibration,” filed on Feb. 16, 2010, and naming Ir Kusnadi et al. as inventors, which application is incorporated entirely herein by reference. 
     
    
     FIELD OF THE INVENTION 
       [0002]    The present invention is directed to the calibration of models used to simulate lithographic processes. Various implementations of the invention may be particularly beneficial to the calibration of lithographic processes simulation models for use in optical proximity correction, such as model-based optical proximity correction. 
       BACKGROUND OF THE INVENTION 
       [0003]    Electronic circuits, such as integrated microcircuits, are used in a variety of products, from automobiles to microwaves to personal computers. Designing and fabricating microcircuit devices typically involves many steps, sometimes referred to as the “design flow.” The particular steps of a design flow often are dependent upon the type of microcircuit, its complexity, the design team, and the microcircuit fabricator or foundry that will manufacture the microcircuit. Typically, software and hardware “tools” verify the design at various stages of the design flow by running software simulators and/or hardware emulators. These steps aid in the discovery of errors in the design, and allow the designers and engineers to correct or otherwise improve the design. These various microcircuits are often referred to as integrated circuits (IC&#39;s). 
         [0004]    Several steps are common to most design flows. Initially, the specification for a new circuit is transformed into a logical design, sometimes referred to as a register transfer level (RTL) description of the circuit. With this logical design, the circuit is described in terms of both the exchange of signals between hardware registers and the logical operations that are performed on those signals. The logical design typically employs a Hardware Design Language (HDL), such as the Very high speed integrated circuit Hardware Design Language (VHDL). The logic of the circuit is then analyzed, to confirm that it will accurately perform the functions desired for the circuit. This analysis is sometimes referred to as “functional verification.” 
         [0005]    After the accuracy of the logical design is confirmed, it is converted into a device design by synthesis software. The device design, which is typically in the form of a schematic or netlist, describes the specific electronic devices (such as transistors, resistors, and capacitors) that will be used in the circuit, along with their interconnections. This device design generally corresponds to the level of representation displayed in conventional circuit diagrams. The relationships between the electronic devices are then analyzed, to confirm that the circuit described by the device design will correctly perform the desired functions. This analysis is sometimes referred to as “formal verification.” Additionally, preliminary timing estimates for portions of the circuit are often made at this stage, using an assumed characteristic speed for each device, and incorporated into the verification process. 
         [0006]    Once the components and their interconnections are established, the design is again transformed, this time into a physical design that describes specific geometric elements. This type of design often is referred to as a “layout” design. The geometric elements, which typically are polygons, define the shapes that will be created in various layers of material to manufacture the circuit. Typically, a designer will select groups of geometric elements representing circuit device components (e.g., contacts, channels, gates, etc.) and place them in a design area. These groups of geometric elements may be custom designed, selected from a library of previously-created designs, or some combination of both. Lines are then routed between the geometric elements, which will form the wiring used to interconnect the electronic devices. Layout tools (often referred to as “place and route” tools), such as Mentor Graphics&#39; IC Station or Cadence&#39;s Virtuoso, are commonly used for both of these tasks. 
         [0007]    IC layout descriptions can be provided in many different formats. The Graphic Data System II (GDSII) format is a popular format for transferring and archiving two-dimensional graphical IC layout data. Among other features, it contains a hierarchy of structures, each structure containing layout elements (e.g., polygons, paths or poly-lines, circles and textboxes). Other formats include an open source format named Open Access, Milkyway by Synopsys, Inc., EDDM by Mentor Graphics, Inc., and the more recent Open Artwork System Interchange Standard (OASIS) proposed by Semiconductor Equipment and Materials International (SEMI). These various industry formats are used to define the geometrical information in IC layout designs that are employed to manufacture integrated circuits. Once the microcircuit device design is finalized, the layout portion of the design can be used by fabrication tools to manufacturer the device using a photolithographic process. 
         [0008]    There are many different fabrication processes for manufacturing a circuit, but most processes include a series of steps that deposit layers of different materials on a substrate, expose specific portions of each layer to radiation, and then etch the exposed (or non-exposed) portions of the layer away. For example, a simple semiconductor device component could be manufactured by the following steps. First, a positive type epitaxial layer is grown on a silicon substrate through chemical vapor deposition. Next, a nitride layer is deposited over the epitaxial layer. Then specific areas of the nitride layer are exposed to radiation, and the exposed areas are etched away, leaving behind exposed areas on the epitaxial layer, (i.e., areas no longer covered by the nitride layer). The exposed areas then are subjected to a diffusion or ion implantation process, causing dopants, for example phosphorus, to enter the exposed epitaxial layer and form charged wells. This process of depositing layers of material on the substrate or subsequent material layers, and then exposing specific patterns to radiation, etching, and dopants or other diffusion materials, is repeated a number of times, allowing the different physical layers of the circuit to be manufactured. 
         [0009]    Each time that a layer of material is exposed to radiation, a mask must be created to expose only the desired areas to the radiation, and to protect the other areas from exposure. The mask is created from circuit layout data. That is, the geometric elements described in layout design data define the relative locations or areas of the circuit device that will be exposed to radiation through the mask. A mask or reticle writing tool is used to create the mask based upon the layout design data, after which the mask can be used in a photolithographic process. The image created in the mask is often referred to as the intended or target image, while the image created on the substrate, by employing the mask in the photolithographic process is referred to as the printed image. 
         [0010]    As designers and manufacturers continue to increase the number of circuit components in a given area and/or shrink the size of circuit components, the shapes reproduced on the substrate (and thus the shapes in the mask) become smaller and are placed closer together. This reduction in feature size increases the difficulty of faithfully reproducing the image intended by the layout design onto the substrate. Adding to the difficulty associated with increasingly smaller feature size is the diffractive effects of light. As light illuminates the mask, the transmitted light diffracts at different angles in different regions of the mask. These effects often result in defects where the intended image is not accurately “printed” onto the substrate during the photolithographic process, creating flaws in the manufactured device. 
         [0011]    To address this problem, one or more resolution enhancement techniques are often employed to improve the resolution of the image that the mask forms on the substrate during the photolithographic process. Examples of various resolution enhancement techniques are discussed in “Resolution Enhancement Technology: The Past, the Present, and Extensions for the Future,” Frank M. Schellenberg, Optical Microlithography XVII, edited by Bruce W. Smith, Proceedings of SPIE Vol. 5377, which article is incorporated entirely herein by reference. One of these techniques, radiation amplitude control, is often facilitated by modifying the layout design data employed to create the lithographic mask. One way to implement this technique, for example, is to adjust the edges of the geometric elements in the layout design so that the mask created from the modified layout data will control the radiation amplitude in a desired way during a lithographic process. The process of modifying the layout design data in this manner is often referred to as “optical proximity correction” or “optical process correction” (OPC). 
         [0012]    As previously noted, a layout design is made up of a variety of geometric elements, which typically are polygons. In some conventional “rule-based” optical proximity correction processes, geometric features that will increase the fidelity of the photolithographic process are automatically be added to geometric elements in the design. For example, some optical proximity correction processes may create serifs at one or more corners of a polygonal geometric element. 
         [0013]    Still other types of optical proximity correction processes, sometimes referred to as “model-based” optical proximity correction processes, fragment the edges of the geometric elements in a design. More particularly, the individual edges of each geometric element are divided into smaller sections, often referred to as edge segments or edge fragments. The size of the fragments and which particular edges are to be fragmented are dependent upon parameters of the optical proximity correction process. The relative positions of the edge fragments relative are changed, and the new edge fragment locations are employed in a simulation of the actual lithographic process that will be used to manufacture the integrated circuit. The printed image produced by the lithographic process simulation is then compared with intended or “target” printed image, to see how closely the simulated printed image matches the target printed image. Based upon the degree of correspondence between the simulated printed image and the target printed image, the relative positions of the edge fragments are modified, and the lithographic process is simulated again using the new edge fragment locations. When the simulated printed image cannot be substantially improved by further displacement of the edge segments, it is often said that the optical proximity correction process has converged. This process of simulation, modification, and simulation is repeated until the simulated printed image sufficiently corresponds to the intended printed image, or until the optical proximity correction process has converged. 
         [0014]    Layout designs can be very large. For example, one layout data file for a single layer of a field programmable gate array may be approximately 58 gigabytes. Accordingly, performing even a single iteration of an optical proximity correction process on a design is computationally intensive. Repeating a model-based optical proximity correction process until the simulated printed image matches the intended printed image, or until the optical proximity correction process has converged, only adds to the time required to finalize the layout design. Often, it can take as many as eight or more iterations for an optical proximity correction process to converge. Due to the number of required iterations of optical proximity correction and the complexity and size of modern layout designs, the time required to perform optical proximity correction is often measured in days. Even where advanced computer processing techniques are employed, performing model-based optical proximity correction may still take days. 
         [0015]    One significant factor in model-based optical proximity correction is the model used to simulate the actual photolithographic manufacturing process. If the model is not accurate, then the corrected design data may still not print the desired image during the actual photolithographic manufacturing process. Accordingly, proper calibration of photolithographic manufacturing process simulation models is important. 
         [0016]    At least two types of techniques for calibrating photolithographic manufacturing process simulation models conventionally have been used. Critical dimension (CD) based model calibration techniques rely on printing simple test structures, such as lines or rectangles of varying sizes, followed by measuring their widths and/or spaces at prescribed locations. CD measurements are typically collected in the form of a table consisting of gauges and their CD values, to which photolithographic manufacturing process simulation models then are calibrated by data fitting. 
         [0017]    A second calibration technique, which uses contours of manufactured structures to augment CD-based model calibration, has been reported for many years now. Contour-based calibration uses contours extracted from top-down scanning electron microscope (SEM) images of printed features on a physical wafer. The contours of the actual manufactured structures then are compared with simulated printed images generated using the model. With contour-based calibration techniques, the calibration process is typically iterative, comparing SEM measured contours as input with simulated print images simulated using photolithographic manufacturing process simulation models that are improved with each iteration. In a hybrid calibration technique, the calibration uses a cost function that includes both CD and contour data. 
         [0018]    Limitations as well as benefits of using either CD-based or contour-based model calibration have been enumerated. One of the primary metrology challenges in contour-based calibration is to properly import and overlay the measured physical contours onto the corresponding calibration features for the comparison. Because contour-based calibration typically includes root-mean-square (RMS) differences as a metric for comparing contour shapes, it is important that the physical contours must be correctly aligned to be properly compared with the simulated printed image. Thus, contour-based calibration is sensitive to input misalignment of the measured contour data. 
         [0019]    Some work has been made on accurately aligning measured contour data for model calibration. During post-metrology/pre-calibration steps, alignment is initially done by pattern matching, and then refined using edge placement error (EPE) criterion, which has improved model calibration quality. When a model for simulating a printed image is not yet available, alignment using the original layout design data (i.e., the target image) is perfectly acceptable, especially for structures with multiple features in the SEM image field of view (FOV). However, design features are not always printable on a wafer, and hence aligning SEM measured contours to the original layout design data is not always straightforward. 
       BRIEF SUMMARY OF THE INVENTION 
       [0020]    Aspects of the invention relate to the calibration cost function and evaluation of a measured contour alignment using the contour of the simulated printed image as a reference. With various implementations of the invention, a contour of a simulated image generated using a given initial (e.g., uncalibrated) model is used as a reference for alignment of measured physical contours. By applying a fitness test to these two sets of contours, evaluation of the fitting of the simulation contour and at the same time the alignment of the measurement contour. This makes it possible to apply correction to the model fitness from within the cost function once misalignment info is known as calibration progresses. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0021]      FIG. 1A  illustrates a rectangular mask feature and a simulated printed image produced by the rectangular mask feature. 
           [0022]      FIG. 1B  illustrates a modified mask feature and a simulated printed image produced by the modified mask feature. 
           [0023]      FIG. 2A  illustrates edge fragmentation of a mask feature. 
           [0024]      FIG. 2B  illustrates target images and simulated printed images for three edge fragments. 
           [0025]      FIG. 2C  illustrates edge fragment position adjustments for improving the simulated printed image of an edge. 
           [0026]      FIG. 2D  illustrates final positions of the edge fragments and a corresponding simulated printed image. 
           [0027]      FIG. 3  illustrates an example of bidirectional contour point error determination. 
           [0028]      FIG. 5  illustrates an example of a contour point error that is Euclidean distance from a point Mi in a measurement contour to the closest point Si in a simulation contour. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     Optical Proximity Correction 
       [0029]    In a photolithographic process, as explained above, electromagnetic radiation is transmitted through selectively transparent areas of a mask. The radiation passing through these transparent areas then irradiates desired portions of a layer of photoresistive material on a semiconductor substrate. The mask in turn is created from layout design data describing the geometric features that should be manufactured on the semiconductor substrate in order to create the desired circuit. For example, if a transistor should have a rectangular gate region, then the layout design data will include a rectangle defining that gate region. This rectangle in the layout design data is then implemented in a mask for creating the rectangular gate region. 
         [0030]    During a photolithographic process, however, optical effects will prevent the shapes defined by the mask from being faithfully imaged onto the substrate. Diffractive effects, for example, may distort the image produced by a mask. Moreover, these distortions become more pronounced as the images produced by the mask become smaller relative to the wavelength of radiation used in the photolithographic process. Thus, the rectangular mask feature  101  illustrated in  FIG. 1  may produce only the image  103 . As seen in this figure, the image  103  is substantially narrower in the corners (e.g., corner  105 ) than the ideal rectangular shape intended by the mask feature  101 . Likewise, the image  103  may have areas (e.g.,  107 ) that extend beyond the ideal rectangular shape intended by the mask feature  101 . Often the intended shape or feature is referred to as the target shape, or the target image, and typically corresponds to the mask feature  101 . Additionally, the image created by employing the mask in a photolithographic process is often referred to as the printed image. 
         [0031]    To correct for these optical distortions, many circuit designers will attempt to modify the layout design data to enhance the resolution of the images that will be produced by the resulting mask during the photolithographic process. Thus, some designers will employ an optical proximity correction (OPC) process on the layout design data, in an effort to better control the amplitude of the radiation transmitted by the mask at specific locations. In a conventional optical proximity correction process, the edges of the geometric elements in the design are fragmented. For example, as shown in  FIG. 2A , an edge of the geometric element  201  used to create the mask feature  101  may be fragmented into edge fragments  201 A- 201 F. The size of the edge fragments in a given layout design depends upon the optical proximity correction process parameters, often referred to as the optical proximity correction recipe. The “recipe” specifies the size of the edge fragments. Accordingly, not all edges within a layout design are fragmented in every optical proximity correction process. 
         [0032]    The optical proximity correction process also simulates the printed image. That is, the photolithographic process is simulated in order to produce a simulated printed image, such as the example image  103  shown in  FIG. 1 . This simulated image is compared to the target image. Typically, this comparison is done at each edge fragment. For example, as shown in  FIG. 2B , the target image is a distance d 1  away from the simulated printed image at the edge fragment  201 A, the target image is a distance d 2  away from the simulated printed image at the edge fragment  201 C, while the target image intersects the simulated printed image at the edge fragment  201 B. The distances between the target image and the simulated printed image are often referred to as the edge placement error (EPE). Accordingly, in most conventional optical proximity correction process each edge fragment or unfragmented edge has an associated edge placement error. 
         [0033]    Next, the edge fragments are individually moved in order to improve the resolution of the simulated printed image for the resulting mask. For example, as shown in  FIG. 2C , the edge fragment  201 A is displaced in a direction away from the geometric element  201 , in an effort to widen the corresponding portion of the image that would be produced by the resulting mask. Similarly, the edge fragment  201 C is displaced in a direction toward from the geometric element  201 , in an effort to narrow the corresponding portion of the image that would produced by the resulting mask. Next, the image that would be produced by a mask using the displaced edge fragments is simulated, and the new simulated image is compared with the target image, and the edge placement errors for each edge fragment are computed. 
         [0034]    This process of moving the edge fragments, simulating the image that would be produced using the moved edge fragments, and comparing the simulated image to the target image may be repeated a number of times. Each cycle of moving edge fragments and comparing the new simulated image to target image is referred to as an iteration of the optical proximity correction process. Typically, edge fragments moved during an given iteration, and the distance the edge fragments are displaced is determined based upon the edge placement error. For example, an optical proximity correction process may move the edge fragments some factor of the edge placement error away from the simulated printed image. Additionally, each edge fragment could be displaced the same distance during a given iteration. The specific parameters than control edge movement is dependant upon the tool used to implement the optical proximity correction process and the optical proximity correction process recipe. 
         [0035]    Typically, these steps will be repeated until the simulated image is sufficiently similar to the target image (e.g., both d 1  and d 2  are smaller than a threshold value), or until it is determined that the displacements of the edge fragments already have converged on locations where no further movement of the edge fragments will improve the simulated image, as shown in  FIG. 2D . Once the final positions of the edge fragments are determined in the layout design data, as shown in  FIG. 2D , a modified mask feature  101 ′ can be created from the corrected layout design data. As shown in  FIG. 1B , the image  103 ′ produced by the modified mask feature  101 ′ should more closely correspond to the target image. 
       Contour Alignment 
       [0036]    A cost function to compare simulation to measured contours is used in contour-based lithographic model calibration. In the model building, calibration is performed by varying the model&#39;s parameters. This produces different simulation contours and hence different model fitness to be used to reduce the contour difference. 
         [0037]    The simulation to measurement contour difference is measured within a bounded neighborhood which typically is a layout area representing the SEM image FOV. This is done in bidirectional fashion from measurement to simulation and vice versa to capture extra print or misprint in either set of contours such as shown in  FIG. 1 . The metric is a weighted RMS of contour point errors each of which represents the distance of a point in one contour to a closest point on the other contour. 
         [0038]      FIG. 3  illustrates the bidirectional contour distance calculation: (a) distance from measurement to simulation, (b) distance from simulation to measurement. Area S is an example where an extra simulation occurs resulting in a larger overall distance in (b) since the closest measurement points belong to different design feature. 
         [0039]    For simplicity, the following description refers to a set of contours covered by single SEM image FOV and error is depicted as single direction from the measurement to simulation. Assuming a perfect alignment, the contour distance from measurement to simulation is given by the following equation (10): 
         [0000]    
       
         
           
             
               
                 
                   
                     E 
                     RMS 
                   
                   = 
                   
                     
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             1 
                           
                           N 
                         
                          
                         
                             
                         
                          
                         
                           
                             w 
                             i 
                           
                            
                           
                             D 
                             i 
                             2 
                           
                         
                       
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             1 
                           
                           N 
                         
                          
                         
                             
                         
                          
                         
                           w 
                           
                             i 
                             | 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where N is the number of points in the measurement contour, w i  the weight, and D i  is the error defined as a Euclidean distance from point i of the measurement contour to the closest point in simulation contour. 
         [0040]      FIG. 5 . An error is Euclidean distance from a point Mi in measurement contour to the closest point Si in simulation contour. Error of point M i  is given by: 
         [0000]        D   i =√{square root over (( S   ix   −M   ix ) 2 +( S   iy   −M   iy ) 2 )}{square root over (( S   ix   −M   ix ) 2 +( S   iy   −M   iy ) 2 )}=√{square root over ( D   ix   2   +D   iy   2 )}
 
         [0041]    Let the measurement contours be misaligned by (T i , 0), then the error at point M i  can be rewritten as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           E 
                           i 
                           2 
                         
                         = 
                           
                          
                         
                           
                             
                               ( 
                               
                                 
                                   S 
                                   ix 
                                 
                                 - 
                                 
                                   M 
                                   ix 
                                 
                                 + 
                                 
                                   T 
                                   x 
                                 
                               
                               ) 
                             
                             2 
                           
                           + 
                           
                             
                               ( 
                               
                                 
                                   S 
                                   iy 
                                 
                                 - 
                                 
                                   M 
                                   iy 
                                 
                                 + 
                                 0 
                               
                               ) 
                             
                             2 
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           
                             
                               ( 
                               
                                 
                                   D 
                                   ix 
                                 
                                 + 
                                 
                                   T 
                                   x 
                                 
                               
                               ) 
                             
                             2 
                           
                           + 
                           
                             D 
                             iy 
                             2 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   30 
                   ) 
                 
               
             
           
         
       
     
         [0042]    Taking into account the misalignment, the contour distance becomes: 
         [0000]    
       
         
           
             
               
                 
                   
                     E 
                     RMS 
                   
                   = 
                   
                     
                       
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               E 
                               i 
                               2 
                             
                           
                         
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             w 
                             i 
                           
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   40 
                   ) 
                 
               
             
           
         
       
     
         [0043]    In this case, the misalignment is equivalent to solving T x  such that E RMS  is minimum that is when the following is minimum: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           E 
                           SSQ 
                         
                         = 
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               E 
                               i 
                               2 
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               ( 
                               
                                 
                                   T 
                                   x 
                                   2 
                                 
                                 + 
                                 
                                   2 
                                    
                                   
                                     D 
                                     ix 
                                   
                                    
                                   
                                     T 
                                     x 
                                   
                                 
                                 + 
                                 
                                   D 
                                   ix 
                                   2 
                                 
                                 + 
                                 
                                   D 
                                   iy 
                                   2 
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                         
                           
                             
                               T 
                               x 
                               2 
                             
                              
                             
                               
                                 ∑ 
                                 
                                   i 
                                   = 
                                   1 
                                 
                                 N 
                               
                                
                               
                                   
                               
                                
                               
                                 w 
                                 i 
                               
                             
                           
                           + 
                           
                             2 
                              
                             
                               T 
                               x 
                             
                              
                             
                               
                                 ∑ 
                                 
                                   i 
                                   = 
                                   1 
                                 
                                 N 
                               
                                
                               
                                   
                               
                                
                               
                                 
                                   w 
                                   i 
                                 
                                  
                                 
                                   D 
                                   ix 
                                 
                               
                             
                           
                           + 
                           
                             
                               ∑ 
                               
                                 i 
                                 = 
                                 1 
                               
                               N 
                             
                              
                             
                                 
                             
                              
                             
                               
                                 w 
                                 i 
                               
                                
                               
                                 ( 
                                 
                                   
                                     D 
                                     ix 
                                     2 
                                   
                                   + 
                                   
                                     D 
                                     iy 
                                     2 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   50 
                   ) 
                 
               
             
           
         
       
     
         [0044]    or, when the following condition is met: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       ∂ 
                       
                         E 
                         SSQ 
                       
                     
                     
                       ∂ 
                       
                         T 
                         x 
                       
                     
                   
                   = 
                   
                     
                       
                         2 
                          
                         
                           T 
                           x 
                         
                          
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             w 
                             i 
                           
                         
                       
                       + 
                       
                         2 
                          
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               D 
                               ix 
                             
                           
                         
                       
                     
                     = 
                     0 
                   
                 
               
               
                 
                   ( 
                   60 
                   ) 
                 
               
             
           
         
       
     
         [0045]    This leads a solution for T x    
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                     x 
                   
                   = 
                   
                     
                       - 
                       
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               D 
                               ix 
                             
                           
                         
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             w 
                             i 
                           
                         
                       
                     
                     = 
                     
                       - 
                       
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             
                               w 
                               i 
                             
                              
                             
                               ( 
                               
                                 
                                   S 
                                   ix 
                                 
                                 - 
                                 
                                   M 
                                   ix 
                                 
                               
                               ) 
                             
                           
                         
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             N 
                           
                            
                           
                               
                           
                            
                           
                             w 
                             i 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   70 
                   ) 
                 
               
             
           
         
       
     
         [0046]    Similarly for misalignment of (0, T y)    
         [0000]    
       
         
           
             
               
                 
                   
                     T 
                     y 
                   
                   = 
                   
                     - 
                     
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             1 
                           
                           N 
                         
                          
                         
                             
                         
                          
                         
                           
                             w 
                             i 
                           
                            
                           
                             ( 
                             
                               
                                 S 
                                 iy 
                               
                               - 
                               
                                 M 
                                 iy 
                               
                             
                             ) 
                           
                         
                       
                       
                         
                           ∑ 
                           
                             i 
                             = 
                             1 
                           
                           N 
                         
                          
                         
                             
                         
                          
                         
                           w 
                           i 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   80 
                   ) 
                 
               
             
           
         
       
     
         [0047]    Algorithmically the (mis)alignment info (in this case T x  and T y ) is obtainable at the same time as calculating contour distance E RMS  if Eq. (80) and thus it serves as supplement to model fitness that is being reported. There are several possibilities as to how the info can be utilized. One option is to use the alignment info simply to supplement the fitness reporting. This is basically saying that the reported fitness for a given model includes some errors due to misalignment. Another option is to use the info to correct the fitness (1) by recalculating fitness each time misalignment occurs, or (2) by using it as a basis for calculating fitness in the next iteration, or (3) by calculating it up front to correct the fitness for rest of iterations. This is stating that the cost function self aligns the input contours as the calibration progresses. 
       Conclusion 
       [0048]    While the invention has been described with respect to specific examples including presently preferred modes of carrying out the invention, 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 invention 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 invention may be implemented using any desired combination of electronic design automation processes.