SYSTEM AND METHOD FOR PRODUCT YIELD PREDICTION USING DEVICE AND PROCESS NEIGHBORHOOD CHARACTERIZATION VEHICLE
This application is a Continuation in Part of U.S. Patent Application 09/442, 699, filed November 18, 1999. This application claims the benefit of U.S. Provisional Application
60/166,307, filed November 18, 1999, and U.S. Provisional Application 60/166,308, filed November 18, 1999.
Background of the Invention The present invention pertains to fabrication of integrated circuits and more particularly to systems and methods for improving fabrication yields.
Silicon based integrated circuit technology, which has evolved into sub-micron line widths, is now able to produce chips containing millions of circuit elements. The process is extremely complex, requiring a large number of process steps to produce multi-level patterns of semiconductor, metals, and insulator types of materials. The circuits are interconnected by metal lines created in multilevel geometries through a large number of very small via holes. Every process step produces three dimensional statistical variations in the geometry and the materials properties of the final configuration. Such statistical variations, which include systematic and/or random defects, can result in both yield and performance degradation of the product. Yield and performance detractors can be found to vary across a chip, across the wafer and from wafer to wafer.
The initial design simulations for an integrated circuit chip, along with expert knowledge of the process capabilities generate a standard cell library that defines the standard device logic, memory and analog unit cells, and the design rules that define the limits and the desired dimensions of the multi-layer film and active device structures. This information is used to generate a mask set for the production of an integrated circuit product. A set of manufacturing process specifications is also generated which describes in detail the multitude of processes associated with each mask level. The mask generated for each process level defines the two dimensions parallel to the Si substrate, i.e. the planar dimensions of each processed layer. The manufacturing process specifications then determine the materials and their properties as well as the third dimension normal to the Si substrate, e.g. diffusion depths, metal thickness, and the thickness of thermally grown and deposited oxides.
For a new chip design there may be a number of iterations before an acceptable product process is defined for stable production. These iterations can include both changes to the mask
set and the manufacturing specifications. The classic s-shaped "learning curve" is a generally accepted concept that models the manufacturing cycle for the release of such high technology type products. The initial flat section of the curve represents the initial trials of the design and process, and generally is considered to represent essentially a very low and inconsistent yield output regime. In this initial stage, some changes to the manufacturing process specifications can be made in order to stabilize the process well enough to obtain a finite but consistent yield result. The so called "ramp up," section of the manufacturing cycle is the section where yield of the product is consistent and is increasing rapidly. The end of the "learning curve" is relatively flat where the product yield is flat and stable. At this stage, the cost of the product is essentially determined by the yield, because all the manufacturing costs are relatively fixed. It is well known that the manufacturing cost of the first two sections of this learning cycle are extremely high because of the amortized cost of multi-billion dollar manufacturing facilities, as well as the cost of highly skilled personnel. Thus, a profit greater than zero must be realized at some point of the "ramp up" cycle, and the projected business profit generally occurs at the beginning of the fixed yield cycle.
For the past thirty years, the integrated circuit technology has been increasing the density of circuits at an exponential rate. This has been accomplished by decreasing the characteristic 'line width" to sub-micron dimensions. Because of this, the economic requirements for the introduction of new products, as well as the maintenance of existing products, have now reached a level of great concern, because they represent very significant cost factors for the industry.
In general, the initial design phases of integrated circuits are optimized with respect to yield and performance factors through the use of elaborate simulation programs rather than through extensive process variations. The driving force for the use of simulation programs rather than the use of process variations is the much higher relative cost to manipulate the process steps.
Elegant simulation programs have also been written in an attempt to correlate classes of observed defects with the design rules and the process steps in order to provide statistically based yield models for integrated circuit products. For example, attempts have been made to predict the yield distributions in a chip, given the data from the mask set. Although these programs have contributed to the knowledge base of the integrated circuit technologies, it has been difficult for such programs to have a direct effect on the yield and performance of a new or established product. This is because it is extremely difficult for a program to represent a class of integrated circuit products when there is an extremely wide variation in the resulting designs. A case in point is the variation in large assemblies of random logic with respect to the distribution of
(multi-level) line lengths, shapes, and via holes in the interconnection scheme of a particular design. It is therefore difficult for such simulation models to take into account this unpredictable variability. One concludes that the art, with respect to simulation programs, has been useful in the creation of the initial product design, but has not proven to be very effective at optimizing the learning curve or enhancing the performance of the original design.
Integrated circuit chips all undergo an extensive test procedure at the wafer level. Such product testers, which are extremely expensive, are designed primarily to test for functionality. Only the nominal performance of the chip can be measured using these results since the probe contact configuration, as well as the limiting capability of the measuring circuits, prevent an accurate measurement of nanosecond switching speeds that can occur during the normal operation the chip. Further, product testers are not able to measure a large number of worst case line and via paths, or worst case gate logic fan-in and fan-out situations in order to evaluate the critical factors of the product design. Therefore, product testers cannot determine the cause of large yield variations when they occur, nor can they provide sufficient information that will lead to the improvement of the existing yield or performance of an integrated circuit product.
The prior art has addressed this problem through the creation of process monitoring circuits that are incorporated within the product chip and utilize the scribe line area, (or the area between the bonding pads) within the chip reticle. Such test configurations are commonly referred to as a "Scribe Line Monitor" (SLM). Early versions of such monitors attempted to extrapolate AC behavior through DC tests made on the test configurations through the use of simulation models. The more recent art has developed AC test methods using internal circuits within the SLM, e.g. ring oscillators and multiplexing functions that can generate limited performance and yield tests of the representative elements. In these cases it is found that the density of circuits in the SLM are not adequate to represent the behavior of a large assembly of dense circuits found on the product chip because of certain optical effects of the masks combined with the photolithography process. U.S. Patent No. 5,703,381, Iwasa, et al., attempted to alleviate this problem by making connections to test transistors that were shared within the product configuration. Other SLM designs included circuits that contain some combinations of line lengths and via holes, chosen to represent some worst case situations that affect circuit delays. Some designs include a number of inverter gates in series for the purpose of measuring the average switching time of the logic elements. U.S. Patent No. (6,124,143), Sugasawara, has also included representations of lines and via holes on more than one level.
The attempts to predict yields through simulation programs, or on chip testing methods , have been made with varying degrees of success.. Thus, there is a need for an improved system and method for integrated circuit product yield prediction.
SUMMARY OF THE INVENTION
A characteπzation vehicle includes a single die, including a device neighborhood testing module for allowing measurement of vaπations in electπcal performance of an active device due to a first plurality of process vaπations withm the active device, and a second plurality of process vaπations affecting an area suπoundmg the active device The device neighborhood testing module includes at least one active test device and an array of dummy devices, within which the active test device is located
BRIEF DESCRIPTION OF THE DRAWINGS
FIGURE 1 is a block diagram depicting the steps performed by a preferred embodiment of the system of the present invention.
FIGURE 2 is a block diagram depicting additional steps performed by the system of the present invention to effect a feedback loop. FIGURE 3 is a plan view of a Device Neighborhood/Process Neighborhood (DN/PN) characteπzation vehicle (CV) according to the present invention.
FIGURE 4 is a plan view of the DN CV shown in FIG. 3. FIGURE 5 is a plan view of the PN CV shown in FIG. 3. FIGURE 6A shows a device under test in the DN CV of FIG 4, surrounded by dummy devices.
FIGURE 6B is an enlarged view of a detail of FIG 6A FIGURE 7 is a diagram of typical DOE (design of expeπments) vaπables for gate length and active width in the DN CV of FIG. 4.
FIGURE 8 is a diagram of a typical DOE vaπables for the gate length and poly spacing for poly gate devices in the DN CV of FIG 4
FIGURES 9A and 9B are diagrams showing a typical (DOE) for the active area including the active width and length, and the active area spacing for the DN CV of FIG. 4.
FIGURE 9C is a representation of the structural dimensions referred to in 9A and 9B for the DN CV of FIG. 4.
FIGURES 10A-10C show special structures for the measurements referenced in Table 3 in the DN CV of FIG. 4.
FIGURE 11 is a plan view of a typical PN module which includes line width structures in PN CV of FIG. 5. FIGURES 12A and 12B show representations of bridge and comb structures in the PN
CV of FIG. 5, referenced in Table 4.
FIGURES 13 A and 13B show structures for the measurement of critical dimension and sheet resistance in the PN CV of FIG. 5, referenced in Table 5.
FIGURE 14 is a plan view of a Kelvin structure used to measure critical dimensions in the PN CV of FIG. 5.
FIGURES 15A-15C show structures used to measure junction leakage in the PN CV of FIG. 5 and referenced in Table 6.
FIGURE 16 shows structures used to study the alignment of contact and vias in the PN CV of FIG. 5. FIGURE 17 shows a structure used to measure the impact of gate to contact spacing in the TESTCHIP module shown in FIG. 5.
FIGURE 18 is a detailed schematic diagram of a basic device neighborhood cell.
FIG. 19 is a diagram showing definitions of "shortable area" in a characterization vehicle.
FIG. 20 is a diagram showing a test pattern for analyzing yield of T-shaped endings. FIG. 21 is a diagram of a nest structure for extracting defect size distributions.
FIG. 22 is a logarithmic diagram of failures plotted against a parameter relating to number of lines shorted, line spacing, and width.
DETAILED DESCRIPTION
U.S. Patent Application 09/442, 699, filed November 18, 1999. U.S. Provisional Patent Application 60/166,307, filed November 18, 1999, and U.S. Provisional Patent Application 60/166,308, filed November 18, 1999 are all hereby incorporated by reference herein in their entireties.
A key element of the exemplary embodiment of the invention is a separate chip or Device Neighborhood/Process Neighborhood Characterization Vehicle (DN/PN CV) that can comprehensively represent any given manufacturing process for various families of integrated
circuits. The DN/PN CV provides a means to make meaningful measurements of the critical factors in the manufacturing process which determine the yield and performance of the product.
Further, the DN/PN CV provides a means to significantly shorten the time dependence of the first two sections of the "learning curve" described above. The DN/PN CV can be used to increase the yield and performance of an integrated circuit product in the case where it has progressed to the final stage of its product cycle. In general, the DN/PN CV is able to be used in the product cycle at any stage of its development in order to provide yield and/or performance enhancements. Therefore, the DN/PN CV provides a new approach with regard to the solution of the most critical problems encountered in the development and the manufacturing of sub- micron integrated circuit technologies.
The DN module allows testing of the impact on device performance a first plurality of process variations affecting the device and a second plurality of process variations affecting other devices suπounding the device to be analyzed. Analogously, the PN module allows testing of the impact on structure electrical characteristics of a first plurality of process variations affecting the structure, and the simultaneous effect of a second plurality of process variations of the other surrounding structures. (The structures may include, for example, combs, snakes and the like.)
Referring now to Figure 1, there is shown a block diagram depicting the steps performed by a system, generally designated 10, for predicting integrated circuit yields in accordance with the present invention. The system 10 utilizes at least one type of characterization vehicle 12.
The characterization vehicle 12 preferably is in the form of software containing information required to build an integrated circuit structure which incorporates at least one specific feature representative of at least one type of feature to be incorporated into the final product. For example, the characterization vehicle 12 might define a test vehicle which can characterize selected lithographic layers for probing the health and manufacturability of the metal interconnection module of the process flow under consideration. The structures need to be large enough and similar enough to the actual product or type of products running in the fabrication process to enable a reliable capture or fingerprint of the various maladies that are likely to affect the product during the manufacturing process. . The characterization vehicle 12 defines features which match one or more attributes of the proposed product layout. For example, the characterization vehicle 12 might define a CV module having a layout which includes features which are representative of the proposed product layout (e.g. examples of line size, spacing and periodicity; line bends and runs; etc.) in order to determine the maladies likely afflicting those specific design types and causing yield loss.
The CV 12 might also define one or more active regions and neighboring features of the proposed design in order to explore impact of layout neighborhood on device performance and process parameters; model device parameters as a function of layout attributes; and determine which device correlate best with product performance. In addition to providing information for assessing and diagnosing yield problems likely to be seen by the product(s) under manufacture, the characterization vehicle is designed to produce yield models 16 which can be used for accurate yield prediction. These yield models 16 can be used for purposes including, but not limited to, product planning, prioritizing yield improvement activities across the entire process, and modifying the original design of the product itself to make it more manufacturable.
The majority of the test structures in the characterization vehicle 12 contemplated in the invention are designed for electrical testing. To this end, the reliability of detecting faults and defects in the modules evaluated by each characterization vehicle is very high. Inspection equipment cannot deliver or promise this high degree of reliability. Furthermore, the speed and volume of data collection is very fast and large respectively since electrical testing is fast and cheap. In this way, statistically valid diagnosis and/or yield models can be realized.
The characterization vehicle 12 is preferably in the form of a GDS 2 layout on a tape or disc which is then used to produce a reticle set. The reticle set is used during the selected portions of the fabrication cycle 14 to produce the yield model 16. Thus the yield model 16 is preferably constructed from data measured from at least a portion of a wafer which has undergone the selected fabrication process steps using the reticle set defined by the characterization vehicle 12.
The yield model 16 not only embodies the layout as defined by the characterization vehicle, it also includes artifacts introduced by the fabrication process operations themselves. The yield model 16 may also include prototype architecture and layout patterns as well as features which facilitate the gathering of electrical test data and testing prototype sections at operating speeds which enhances the accuracy and reliability of yield predictions.
An extraction engine 18 is a tool for extracting layout attributes from a proposed product layout 20 and plugging this information into the yield model 16 to obtain a product yield prediction 22. Such layout attributes might include, for example, via redundancy, critical area, net length distribution, and line width/space distribution. Then, given layout attributes from the proposed product layout 20 and data from yield models 16 which have been fabricated based upon information from the characterization vehicles 12, product yield 22 is predicted. Using the system and method of the present invention, the predictable product yield obtainable can be that
.
associated with each defined attribute, functional block, or layer, or the resultant yield prediction for the entire product layout.
Referring now to Figure 2, there is shown a block diagram of the system for predicting integrated circuit yields 10 in accordance with the present invention additionally comprising a feedback loop, generally designated 24, for extracting design attributes 26 from product layout
20 by means of extraction engine 28. In accordance with this feature of the present invention, the characterization vehicle 12 is developed using attributes of the product layout 20. In this case, attributes of the product layout are extracted, making sure that the range of attributes are spanned in the characterization vehicle 12. For example, the product layout is analyzed to determine line space distribution, width distribution, density distribution, the number of island patterns, in effect developing a subset of the entire set of design rules of the fabrication process, which subset is applicable to the particular product layout under consideration. With respect to patterns, the product layout analysis would determine the most common pattern, the second most common pattern, and so forth. These would be extracted by the extraction engine 28 yielding design attributes 26 encompassing all of these patterns for inclusion into the characterization vehicle 12.
With respect to densities, if the analysis of the product layout reveals that the density of a first metal is from 10% to 50%, then the characterization vehicle would include the entire range of 10% to 50% for the first metal.
Figure 3 shows a layout of the DN/PN CV test vehicle that includes the device neighborhood and process neighborhood test engineering groups (DN and PN CVs), and other special structure modules. The device neighborhood (DN) subsection incorporates 200 or more discreet NMOS and PMOS transistor structures. The purpose of the DN CV is to quantify the impact of within die pattern variations on the performance of the transistor. Placing the test devices within an aπay of dummy transistors derives the local pattern dependencies. The principal variants in the DN CV are polysilicon and active area density. In addition, special structures are included to investigate the impact of contact and gate spacing.
The PN CV portion of the design is typically comprised of 17 or more sub-modules. The area containing these modules represents a size comparable to a large portion of a product chip.. The purpose of these modules is to measure the impact of within die pattern variations on systematic yield loss associated with bridging, junction leakage, CD line width variation and contact and via formation. Like the DN CV, the PN CV accomplishes pattern-dependent characterization by varying the local process environment around the test structures. In this manner, the PN CV can effectively simulate the within-die pattern dependent process variations
expected during normal manufactuπng. By accounting for the effects of the process neighborhood, both the DN and PN CVs offer a much more comprehensive vehicle for debugging and ramping processes into production. Separate from the CVs, a special structure module provides a vehicle to characterize well and channel profiles as well as N+ and P+ sheet resistances, while other structures focuss on the effects of contact to gate spacing.
2.0 DN/PN CV Floorplan and Description
In this section, the floorplan and basic description of each module of the DN/PN CV is described. Detailed descriptions of each module discussed in this section can be found in section 3.
As is shown in Figure 3, the DN/PN CV test chip is composed of four major modules, the DN, PN, CVs MOSCAP and TESTCHIP... The DN module measures the variation in electrical device performance as a function of local process environment. The DN module is split into two types sub-modules: NMOS and PMOS and three process effects sub-modules: poly variant, active_variant, and special_structure.
Figure 4 is a plan view of the DN CV. The poly variant submodule quantifies the relationship between polysilicon density and electrical device performance. Both macro_ and micro_effects are studied. Similar to the poly variant module, the active_variant structures study the influence of active_area density on the electrical performance of the transistor. Finally, in the special_structures module, the effects of contact_to_poly, contact_to_contact and poly_to_poly spacing are studied along with the influence of contact size.
The PN test vehicle, shown in Figure 5, is split into four types of effects modules: bridging_structures, CD (critical dιmension),_and_sheet_resιstance, junction_leakage, and contact_and_via. The bridging_structures devices consist of five modules which investigate poly_to_poly bridging on field oxide, poly_to_active bridging (N and P) and poly_to_poly bridging due to stringers at the field oxide edge (N and P). The CD_and_sheet_resistance structures measure CD control (e.g. line width), and sheet resistance of active (N+ and P+), gate (N+ and P+), and metal 1, metal 2. metal 3 and metal 4 interconnection layers. In total, this part of the chip encompasses 4 sub-modules and 62 individual test structures. The junction_leakage
section is composed of 2 modules with a total of 24 devices, which quantify the area leakage of N+ and P+ active as well as the perimeter leakage for active area bounded by gate and field. Finally, the contact_and_via submodule is directed at establishing guidelines for all contact and via formation for the four_layer metal process. This section consists of six sub-modules containing a total of 90 test structures.
The MOSCAP and TESTCHTP modules, also shown in Figure 5, contain specialty devices that measure the active sheet resistance and well_ and channel_profiles as well as structures to determine the yield impact of contact to poly gate spacing.
3.0 Module and Structure Description
In this section, the composition and structure of the DN PN CV is discussed. The DN and PN CVs, MOSCAP and TESTCHTP modules are addressed in Sections 3.1 through 3.4, respectively. Section 4 provides specific device locations, pad definitions and measurement parameters.
3.1 DN CV
The purpose of the DN CV is to measure the effects of the local process environment on the final electrical performance of an MOS transistor. This is accomplished by placing the device under test in a "sea" of dummy devices which can be tailored to provide specific values for parameters such as poly and active area density. An example of such a layout is provided in Figure 6. In this figure, a large area of filled with dummy structures. The device under test is located at the center of the array of dummy devices. In this manner, the local process environment is effectively controlled. The DN module is composed of two type modules, NMOS and PMOS, which are identical in composition and consist of three effects sub-modules, poly variant, active_variant, and special_structure. Each sub-module is described in detail in Sections 3.1.1 through 3.1.3 below.
3.1.1 Poly_Variant
Polysilicon etch loading can potentially impact final circuit performance by creating unanticipated and un-testable variation in the gate length within the die due to both macro and
micro-density effects. The poly variant module incorporates 54 transistor structures in a design of experiments (DOE) that investigates the impact of micro_and macro-polysilicon density on gate length as a function of gate length and width. The DOE for gate length and active width is summarized in Figure 7. Six sets of transistors incorporate these gate/active_width variations and form another designed experiment in which the macroscopic polysilicon density is varied. Macroscopic variations are produced by tailoring the length and width of the gate on the dummy transistors within each test structure, while the micro-density effects are controlled by varying the number of gate fingers which are in close proximity. The list of approximate polysilicon densities and multi-gate variations is provided in Table 1.
Table 1
In addition to the macroscopic loading effects, micro-loading during the gate etch is investigated by the use of multi-gate structures. In this experiment, the active width is kept constant at 2 μm while the gate length, gate spacing and number of gates is varied according to Figure 8. The effects of macro-loading are also investigated by again varying the local polysilicon density. The polysilicon density variations for the multi-finger structures is shown in Table 2.Table 2 lists Target polysilicon densities for polysilicon multi-gateDOE in the PN-CV module.
Table 2
3.1.2 Active Variant
Active etch loading can potentially impact final circuit performance by creating unanticipated and un-quantifiable vaπations in the shallow_trench_ιsolatιon profiles. These vaπations can manifest themselves in systematic yield loss due to poly-bπdgmg, excessive leakage, etc. The active vaπant module incorporates 45 transistor structures in a design of expeπments (DOE) that investigates the impact of active density on device performance as a function of gate length and width. The DOE for gate length and active width is summaπzed in Figures 9 A and 9B. The active density is controlled by tailoπng the X and Y spacing on the dummy transistors within each test structure. The DOE of active spacmgs. Sx and Sy, is given in Figure 9C. For this set of test structures the gate length of the dummy transistors is fixed.
3.1.3 Specιal_Structures
In addition to polysilicon and active area density, the DN CV investigates the effects of contact- to-poly, contact-to-contact and poly-to-poly spacmgs, as well as the impact of contact size. This
DOE is contained with the specιal_structures module The specifics of the vaπations are listed m Table 3 and Figures lOA-lOC. Table 3 provides a Summary of a DOE for cπtical dimensions of the special structure and devices m the PN-CV module
Table 3
FIG. 18 is a plan view showing a basic DN cell 1800 DN cell 1800 has a structure 1802 at its center that can be used in a device and Kelvin sheet resistance measurement, to measure a sheet resistance of the gate of the device 1804. Structure 1802 includes a stπpe 1806 which is an extension of the gate 1805 Stπpe 1806 has four contacts 1808 at its ends Cuπent is injected into two of the contacts 1808, and the voltage drop is measured across the other two contacts The device 1804 can be tested for functionality and cuπent can also be injected through the gate
line 1805, to measure voltage drop at the ends of the gate line. This will give a measure proportional to w, which in combination with a Van der Pauw measurement, can give sheet resistance.
3.2 PN CV
The purpose of the PN CV is to vary the environment around a given test structure to stimulate within_die pattern dependent variation. Like the DN CV, the local environment is controlled over a large area, comparable to a large area on a product chip.. An example of such a layout is provided in Figure 11. In this figure, a structure to measure line- width is surrounded by a comb-in-comb which dictates the local etch environment. By varying the line-width and space of the comb, the local process environment is effectively controlled and CD can be coπelated with process layer density. In this manner, the test structure is a more accurate predictor for results on actual product.
The PN CV module is composed of four major sub-modules: bridging_structures, CD_and_sheet_resistance, junction_leakage, and contact_and_via. Each sub-module is described in detail in Sections 3.2.1 through 3.2.4 below.
3.2.1 Bridging Structures
The bridging _structures module, composed entirely of snakes and combs, is incorporated to measure the effects of polysilicon and active area pitch on systematic yield loss due to polysilicon bridging. The module is divided into five sub-modules: SNK_CMB_FLD, SNK_CMB_PAA, SNK_CMB_NAA, SNK_CMB_PAA_FLD, and SNK_CMB_NAA_FLD.
The SNK_CMB_FLD module looks at poly-to-poly bridging on field oxide. SNK_CMB_PAA and SNK_CMB_NAA investigate yield loss mechanisms due to polysilicon_to_active bridging across the sidewall spacer during the silicidation process. SNK_CMB_PAA_FLD and SNK_CMB_NAA_FLD focus on potential problems due to polysilicon-to-polysilicon bridging along a field edge caused by shallow trench isolation issues. Each of these sub-modules represents a DOE in which the dominant line width and space is varied. A summary of each of the modules is provided in Table 4. Visual references for the parameter definitions are provided in Figures 12A and 12B
Table 4
3.2.2 CD and sheet resistance
The CD_and_sheet_resistance module is incorporated to measure the effects of etch loading on the CD of a given layer. Sheet resistance is included to provide necessary input for the line- width measurement structure. The CD_and_sheet _resistance is composed of four sub-modules:
KELVΓN_NGC_MI, KELVΓN_MI_PGC,
KELVIN_M2_M3 and KELVIN_M3_M4. Line-width structures are included for each of the interconnect layers. These include N+ poly (NGC), P+ poly (PGC), Ml, M2, M3 and M4. Van der Pauw structures are included for each of the above layers as well as for N+ active (NAA) and P+ active (PAA). Table 5 provides the specifics for each of the DOEs, including a summary of CD and Sheet Resistance DOE incorporated in the PN-CV.
Table 5
Figures 13 A and 13B provide a visual definition of the parameters investigated. The structure used to measure CD is shown in Figure 14. In this structure, a known current is forced from contact II to 12 and the coπesponding voltage drop is measured at contacts VI and V2. The equation which governs the voltage drop is
where DSheet is the measured sheet resistance, W is the line width, L is the spacing between the contact pads 2 and 3, V23 is the measured voltage and 114 is the measured current. Within the CD module L is set to the same value for all test structures and DSheet is determined using the van der Pauw structures. Thus equation (1) can be re-arranged as
Since all of the values on the right_hand_side of equation (2) are known, accurate determination of the line width can be made.
Extraction of the sheet resistance is accomplished using the Van der Pauw structure illustrated in Figures 13A and 13B. In this structure a current is forced between contacts 1 and 2 while the voltage is measured between pads 3 and 4. Because the structure is symmetric, the sheet resistance can be determined using the following equation
R v;. ^.. = ^(2)7 /.,* (3)
Special attention must be paid to the cuπent used during the measurement as wafer heating (metal and silicon) and minority carrier injection (silicon) can skew the measurement results. Currents for these structures are outlined in Section 4.
3.2.3 Junction_Leakage
The junction_leakage module is implemented within the PN CV to resolve the measured reverse_bias leakage cuπent into its primary components of area, perimeter and corner for both field_ and gate_bcmnded diodes. The module also provides data on the influence of line width and pitch on the leakage current. The module consists of two identical sub-modules, NAA_JCT_LKG and PAA_JCT_LKG that incorporate N+/P and P+/N diodes, respectively. Table 6 lists the design attributes of the diodes contained within the junction eakage module. The parameters are defined schematically in Figures 15A-15C. Table 6 includes a summary of junction leakage DOE in the PN-CV module.
Table 6
In the leakage current analysis, area leakage is extracted using a rectangular diode with a low perimeter_to_area ratio (approximately 6%). Consequently, the measured leakage current is dominated by the diode area and is used to approximate the area leakage component. Perimeter leakage is extracted using the known area leakage cuπent in combination with a series of diodes that incorporate varying area_to_perimeter ratios. In the perimeter analysis, the calculated area leakage current is subtracted from the leakage current measured in the perimeter structure. The residual amount of current, which can be attributed to perimeter leakage, is divided by the diode perimeter value. The resulting perimeter leakage, in units of A/μm, extracted from the diode series is compared. If the values are approximately equal, the perimeter component value is assumed coπect. If the values differ significantly, a more thorough analysis must be performed in which the leakage cuπent is expressed as a function of area, perimeter and corner for each diode and the equations are solved simultaneously. A detailed breakdown of the area (A), perimeter (P) and number of corners (C) contained in each diode is provided in Section 4.
3.2.4 Contact and Via
The contact_and_via module is included in the PN CV to study systematic yield issues related to each of the contact_and _via etch steps included in the CMOSIS full_flow process. The module consists of six sub-modules: CONT_NAA, CONT_PAA, CONT_NGC, VIAl, VIA2 and VIA3.
Each sub-module includes 15 contact or via structures for a total of 90. The test structures are split between two designed experiments. The first measures the effects of line_width, redundancy and spacing on via yield. This DOE is given in Table 7, including DOE for contact and via module in the PN-CV. In this Table DR is the minimum design rule and WIDE denotes some value significantly larger that the minimum dimension.
Table 7
The parameter space for this designed experiment is shown in Table 8, including Summary of a DOE for via mis-alignment structures in the DN-CV module. The dimensions referred to in
Table 8 are defined in Figure 16.
Table 8
3.3 MOSCAP
The purpose of the MOSCAP module is to provide sheet resistance and capacitance measurements of individual implants not addressed in the PN CV. MOSCAP is composed of a single module containing five sheet resistance and 12 capacitor structures. The characteristics of the test structures are summarized in Table 9, including DOE for sheet resistance of source drain and capacity over implants in the DN-CV module.
Table 9
Sheet resistance is measured using the van der Pauw technique described in Section 3.2.2. Specifically, the van der Pauw structures in the MOSCAP measure the resistance values for the deep source/drain and shallow extension implants. The capacitors are incorporated to measure the capacitance from gate_to_well and gate_to_well+channel in a standard_oxide_thickness_device as well as the capacitance from gate to well+channel for the thicker_oxide, high_voltage transistors.
3.4 TESTCHIP
The purpose of the TESTCHIP module is to measure the impact of gate_to_contact spacing. The module incorporates 16 devices (8 NMOS/8 PMOS) in a single column of the CV. The structure used for this characterization incorporates asymetric contact spacing as is illustrated in Figure
17. Table 10 gives the specific range of values for gate-to-contact spacing investigated in the experimental design.
Table 10
The structure utilizes a modified comb_in_comb arrangement. Each comb structure incorporates 80 gate fingers of 57.6 μm length providing a total critical perimeter of 4608 μm. Device test is performed by placing a voltage across the gate and active contacts and measuring the current between the contacts.
4.0 Testing Procedures
In this section, testing procedures for each module in the DN/PN CV are presented. The pad aπangement and numbering, shown in Figure 18, is the same for all sub-modules. The individual measurements necessary to test the DN_ and PN CVs, MOSCAP and TESTCHIP modules are addressed in Sections 4.1 through 4.4, respectively.
4.1 DN CV Measurements
In the DN CV a single test configuration is required for all devices. Measurements are made using Toshiba standard MOS transistor current_voltage recipes which should vary the GATE and DRAIN voltages (VG and VD) while measuring the DRAIN cuπent (ID).
For 13 of the 15 pad frames, pads <1> and <17> represent a common contact for all transistors to the WELL and SOURCE, respectively. The remaining pads are arranged in GATE_DRADSf pairs. GATE contacts occupy pads <2> through <16>, while DRAIN contacts occupy pads <18> through <32>. The GATEJDRATN pairs are aπanged as
<GATE_PAD_NUMBER>_<GATE_PAD_NUMBER+16>.
The two exceptions to the standard aπangement are Column 8 and Column 15 (refer to the DN CV map provided in Figure 1 for column numbering scheme). Column 8 contains both NMOS and PMOS transistors. As a result, pads <1> and <17> provide connections to the PMOS WELL and SOURCE, respectively, while pads <10> and <26> provide contact to the NMOS WELL and SOURCE. Since Column 15 is unfilled, the WELL and SOURCE contacts are located at pads <3> and <19>, respectively. Pads <1>,<2>,<17> and <18> are not used.
4.2 PN CV Measurements
The test configurations for each of the 17 PN CV modules is in Sections 4.2.1 through 4.2.4. For each measurement, the pads involved as well as the measurement conditions on each pad (voltages, cuπents, etc.) are listed.
4.2.1 snake/comb configuration SNK CMB FLD:
SNK_CMB_NAA, SNK_CMB_PAA, SNK_CMB_NAA_FLD, SNK_CMB_PAA_FLD:
4.2.2 junction leakage configuration NAA JCT LKG:
M&AUlHWia ■ ,KM-Bffai-.tJΛl
PAA JCT LKG:
4.2.3 contact/via chain configuration
CONT_NAA, CONT .PAA, CONT_NGC, VIA2, VIA3:
VIA1:
4.2.4 Kelvin configuration
KELVIN NGC Ml:
Measurement
KELVIN PGC Ml:
KELVIN M3 M4:
4.3 MOSCAP Measurements
In the MOSCAP one test configuration is required. This configuration is shown in below.
4.4 TESTCHIP Measurements
In the TESTCHIP one test configuration is required. This configuration is shown in below.
Test Name ■y.aMW.T.'.'HM est 1 «-*ι> tø ix}, <n* 1 β> 10 V::ι>, measure c«*f *Λt Jlowwig through <n , n=l ihrough 16,
By including both the device neighborhood and process neighborhood modules on a single die, device measurements can be correlated with the process neighborhood measurements, so that simulation models used to infer device performance can be verified.
Other characterization vehicles for via, device, suicides, poly, el al, are often designed and utilized. However, the procedure and techniques for designing them are the same. For purposes of illustration, the example metal characterization vehicle will be carried through on extraction engines and yield models.
Reference is again made to FIGS. 1 and 2. The extraction engine 18 has two main purposes: (1) it is used in determining the range of levels (e.g. linewidth, linespace, density) to use when designing a characterization vehicle. (2) It is used to extract the attributes of a product layout which are then subsequently used in the yield models to predict yield. (1) has already been described above with reference to how the line width, space and density of the snake, comb and Kelvin structures were chosen in the example characterization vehicle. Thus, most of the following discussion focuses on (2).
Since there are nearly infinite numbers of attributes that can be extracted from the product layout, it is impossible to list or extract all of them for each product. Thus, a procedure is required to guide which attributes should be extracted. Usually, the characterization vehicle drives which attributes to extract. The process consists of: 1. List all structures in the characterization vehicle
2. Classify each structure into groups or families such that all structures in the family form an experiment over a particular attribute. For example, in the metal characterization vehicle discussed above, a table of family classifications might be:
Family Attributes Explored
Nest structures Basic defectivity over a few linewidths and spaces
Snakes and Combs Yield over wide range of linewidths and spaces including very large widths next to small spaces and very large spaces next to small widths.
Kelvin_CD CD variation across density, linewidth, and linespace.
+ van der Pauws
Border structures Effect of different OPC schemes on yield.
3. For each family, determine which attributes must be extracted from the product layout. The exact attributes to choose are driven from which attributes are explored. For example, if a particular family explores yield over different ranges of space, then either a histogram of spaces or the shortable area for each space must be extracted. For the above example, the required list of attributes might be:
Family Attributes Explored Attributes to Extract from
Product Layout
(A) Nest structures Basic defectivity over a few Critical area curves, linewidths and spaces.
(B) Snakes and combs Yield over wide range of Shortable area and/or instance counts linewidths and spaces for each line width and space including ... explored in the characterization vehicle.
(C ) Kelvin_CD and CD variation across density, Histograms of pattern density, van der Pauws linewidth, and space linewidth, and linespace
(D) Border structures Effect of different OPC For each OPC scheme selected to schemes on yield use on product layout, the shortable area or instance count.
4. Use the attributes extracted in the appropriate yield models as previously described.
For other characterization vehicles, the families and required attributes will obviously be different. However, the procedure and implementation is similar to the example described above.
As stated above, the yield model 16 is preferably constructed from data measured from at least a portion of a wafer which has undergone the selected fabrication process steps using the reticle set defined by the characterization vehicle 12. In the preferred embodiment, the yield is modeled as a product of random and systematic components:
The methods and techniques for determining Ysj. and Yrj are as follows.
SYSTEMATIC YIELD MODELING
Since there are so many types of systematic yield loss mechanisms and they vary from fab to fab, it is not practicable to list every possible systematic yield model. However, the following describes two very general techniques and gives an example of their use especially within the context of characterization vehicles and the methodology described herein.
AREA BASED MODELS The area based model can be written as:
Where q is a design factor explored in the characterization vehicle such as line width, line space, length, ratio of width/space, density, etc. Y0(q) is the yield of a structure with design factor q from the characterization vehicle. A0(q) is the shortable area of this structure and A(q) is the shortable area of all instances of type q on the product layout. Yr(q) is the predicted yield of this structure assuming random defects were the only yield loss mechanism. The procedure for calculating this quantity is described below in connection with random yield modeling.
The definition of shortable area is best illustrated with the example shown in Figure 19. This type of test structure can be used to determine if the fab is capable of yielding wide lines that have a bend .with a spacing of s. In this sample test structure, a short is measured by applying a voltage between terminal (1) and (2) and measuring the current flowing from terminal (1) to (2). If this current is larger than a specified threshold (usually 1-lOOnA), a short is detected. The shortable area is defined to be the area where if a bridging occurs, a short will be measured. In the example of Figure 19, the shortable area is approximately x*s). The A(q) term is the shortable area of all occurrences of the exact or nearly exact patten (i.e. a large line with a spacing of s and a bend of 45 degrees) shown in Figure 19 in a product layout. The Yr(q) term is extracted by predicting the random yield limit of this particular structure using the critical area method described below.
It is important to realize that the effectiveness of this model is only as good as the number of structures and size of structures placed on the characterization vehicle. For example, if the angled bend test structure shown in Figure 19 were never put on the characterization vehicle or was not placed frequently enough to get a meaningful yield number, then there would be no hope of modeling the yield loss of wide line bends on the product layout. While it is difficult to define exactly how many of how big the test structure should be on the characterization vehicle, practical experience has shown that the total shortable area of each test structure on the characterization vehicle should ideally be such that A(q)/Ao(q)<10.
The above discussion has concentrated on shorts since they generally tend to dominate over open yield loss mechanisms. However, open yield loss mechanisms can be modeled equally well with this yield model so long as shortable area is replaced by open causing area.
INSTANCE BASED YIELD MODEL The general form of the instance based yield model is:
Where Yo(q) and Yr(q) are exactly the same as in the area based yield model. Ni(q) is the number of times the unit cell pattern or very similar unit cell pattern to the test pattern on the characterization vehicle appears on the product layout. No(q) is the number of times the unit cell pattern appears on the characterization vehicle.
For example, Figure 20 shows a simple test pattern for examining the yield of T-shaped endings at the ends of lines near a space of s. This test pattern is measured by applying a voltage across terminals (1) and (2) and measuring the shorting current. If this pattern was repeated 25 times somewhere on the characterization vehicle, then No(q) would be 25x5=125 since there are five unit cells per each test structure.
If the number of times this unit cell occurs with a spacing of s near it is extracted from the product layout, the systematic yield of this type of structure can be predicted. For example, if there are five structures with 500 unit cells in each structure then No(q)=2500. If Ni(q) from some product was 10,000 and a yield of the test structures on the characterization vehicle of 98.20% was measured. Using the techniques described below, Yr(q) can be estimated as 99.67%. Using these numbers in the equation:
0.9820 25CD
», = =92.84% 0.9967J
RANDOM YIELD MODELING
The random component can be written as:
Where CA(x) is the critical area of defect size x and DSD(x) is the defective size distribution, as also described in"Modeling of Lithography Related Yield Losses for CAD of NSLI Circuits", W. Maly, IEEE Trans, on CAD, July 1985, ppl61-177, which is incorporated by reference as if fully set forth herein. Xo is the smallest defect size which can be confidently observed or measured.
This is usually set at the minimum line space design rule. The critical area is the area where if a defect of size x landed, a short would occur. For very small x, the critical area is near 0 while very large defect sizes have a critical area approaching the entire area of the chip. Additional description of critical area and extraction techniques can be found in P. K. Nag and W. Maly, "Yield Estimation of VLSI Circuits," Techcon90, Oct. 16-18, 1990. San Jose; P. K. Nag and W.
Maly, "Hierarchical Extraction of Critical Area for Shorts in Very Large ICs," in Proceedings of The IEEE International Workshop on Detect and Fault Tolerance in VLSI Systems, IEEE Computer Society Press 1995, pp. 10-18; I. Bubel, W. Maly, T. Waas, P. K. Nag, H. Hartmann, D. Schmitt-Landsiedel and S. Griep, "AFFCCA: A Tool for Critical Area Analysis with Circular Defects and Lithography Deformed Layout," in Proceedings of The IEEE International
Workshop on Detect and Fault Tolerance in VLSI Systems, IEEE Computer Society Press 1995, pp. 19-27; C. Ouyang and W. Maly, "Efficient Extraction of Critical Area in Large V1SI ICs," Proc. IEEE International Symposium on Semiconductor Manufacturing, 1996, pp. 301-304; C. Ouyang, W. Pleskacz, and W. Maly, "Extraction of Critical Area for Opens in Large VLSI Circuits," Proc. IEEE International Workshop on Defect and Fault Tolerance of VLSI Systems,
1996, pp. 21-29, all of which references are incorporated in this detailed description as if fully set forth herein.
The defect size distribution represents the defect density of defects of size x. There are many proposed models for defect size distributions (see, for example, "Yield Models - Comparative Study", W. Maly, Defect and Fault Tolerance in VLSI Systems, Ed. by C. Stapper, et al, Plenum Press, New York, 1990; and "Modeling of Integrated Circuit Defect Sensitivities", C.H. Stapper, IBM J. Res. Develop., Vol. 27, No. 6, November, 1983, both of which are incorporated by reference as if fully set forth herein), but for purposes of illustrations, the most common distribution:
D x k ∑SD(x) = =2=
will be used where Do represents the total number of defects/cm2 greater than x0 observed. P is a unitless value which represents the rate at which defects decay over size. Typically, p is between 2 and 4. K is a normalization factor such that
The following two sections describe techniques for extracting defect size distributions from characterization vehicles.
THE NEST STRUCTURE TECHNIQUE The nest structure is designed for extracting defect size distributions. It is composed of N lines of width w and space s as shown in Figure 21. This structure is tested by measuring the shorting current between lines 1 and 2, 2 and 3, 3 and 4, ..., and N-l and N. Any current above a given spec limit is deemed a short. In addition, opens can be testing by measuring the resistance of lines 1, 2, 3, ...., N-l, and N. Any resistance above a certain spec limit is deemed to be an open line. By examining how many lines are shorted together the defect size distribution can be determined.
If only two lines are shorted then the defect size must be greater than s and no larger than 3w + 2s. Any defects smaller than s will not cause a short at all while defects larger than 3w+2s are guaranteed to cause a short of at least 3 lines. For each number of lines shorted, an interval of sizes can be created:
It should be noted that the intervals overlap; thus, a defect size distribution cannot be directly computed. This restriction only places a limit on p extraction. Thus, in order to estimate p, a p estimate is computed from the distribution from all the even number lines and then from all the odd number lines. Finally, the two values are averaged together to estimate p. To extract p, the In (number of faults for x lines shorted) vs log ([x-l]s + [x-2]w) is plotted. It can be shown that the slope of this line is -p. The Do term is extracted by counting the number of failures at each grouping of lines and dividing by the area of the structure. However, for very large Do, this estimate will be too optimistic. Additional information on extracing defect size distribution from structures similar to the test structures can be found, for example, in "Extraction of Defect Size Distribution in an IC Layer Using Test Structure Data", J. Khare, W. Maly and M.E. Thomas, IEEE Transactions on Semiconductor Manufacturing, pp. 354-368, Vol. 7, No. 3, August, 1994, which is incorporated by reference as if fully set forth herein.
As an example, consider the following data taken from 1 wafer of 100 dies:
If the structure size is 1cm2 then the Do would be 98 + 11 +4 + 2 +1 = 133 / (100 * 1) =
1.33 defects/cm . Also, the plot of log (number of failures) vs log ([x-l]s + [x-2]w) (see Figure 22) shows that p=2.05.
THE COMB STRUCTURE TECHNIQUE Assuming a comb of width = space = s, it can be shown that the yield of this structure can be written as:
ln[| ln(r>|] = ln[ - f DSD(x)
* (1 -p) * ln(*)
*b
Thus, from the slope of the plot of In[ |ln(Y) | ] vs. ln(s), p can be estimated. The Do extraction technique is the same technique as mentioned above.
YIELD IMPACT AND ASSESSMENT
Once a sufficient number of characterization vehicles has been run and yield estimates are made for each characterization vehicle, the results are placed in a spread sheet to enable prioritization of yield activities. Tables XIV through XVI are examples of information contained in such a spread sheet. It has been divided into sections of metal yield, poly and active area (AA) yield (Table XIV), contact and via yield (Table XV), and device yield (Table XVI). The columns on the left indicate systematic yield loss mechanisms while the columns on the right indicate random yield loss mechanisms. Although the exact type of systematic failure mechanisms vary from product to product, and technology by technology, examples are shown in Tables XIV through XVI. Usually, targets are ascribed to each module listed in the spread sheet. The further a module yield is away from a target, the more emphasis and resources are devoted to fixing the problem. For example, if the target was set artificially at 95 percent for each module in the example shown in Tables XIV through XVI, then clearly (M2 H> M3 ) vias (75.12%) followed by similar vias (Mi =<> M2 ) (81.92%) , Mi shorts (82.25%), and contacts to poly (87.22%) are below target and, with vias (M2 =» M3) needing the most amount of work and contacts to poly needing the least amount of work.
Within each module, it is also possible to tell where the greatest yield loss is situated. That is, is it one particular systematic mechanism being the yield down or is it merely a random defectivity problem, or is it some combination of the two? For example, as shown in Table XV, via (M2 =o M3) yield loss is clearly dominated by a systematic problem affecting vias connected to long metal runners on the M3 level (77.40%). Vias from (Mi =*>M ) are affected by the same problems (91.52%) in addition to a random defectivity problem (92.49%). Solving vias (Mi .© M2 ) yield problems would require fixing both of these problems.
As shown in Table XIV, Mi yield loss is also dominated by a random defectivity issue (85.23%) in addition to a systematic problem affecting wide lines near small spaces (96.66%).
Fixing both of these problems would be required for improving Metal 1. Similar conclusions can be made for other modules in the spread sheet.
For the worst yielding modules, frequent running of further characterization vehicles for this module would be required. Usually, splits will be done on these characterization vehicles to try and improve and validate those improvements in module yield. For the modules which are within target, routine monitoring of short flow characterization vehicles would still be required to validate that there has been no down turn or other movement in module yield. However, these characterization vehicles can be run less frequently than for those modules with known problems.
Not fϊirnisched at time of publication
Not flirnisched at time of publication
TABLE XIV
TABLE XV
99.71%
99.53%
TABLE XVI
Although the invention has been described in terms of exemplary embodiments, it is not limited thereto. Rather, the appended claim should be construed broadly, to include other variants and embodiments of the invention which may be made by those skilled in the art without departing from the scope and range of equivalents of the invention.