Abstract:
The present invention provides a device, method and storage medium, stored with software programs, which, when a memory LSI defect analysis apparatus is used as a monitoring device to estimate reductions in yield, can shorten the time needed for full manual interpretation of the obtained results, by automatically interpreting the analyzed results obtained, and calculating the period of distribution patterns and the mix rate of regular patterned defects. First in defect number calculation process  71 , the total defect number of bits is found; and in factor selection process  61 , the factor f is then selected. Next, in regular patterned decision process  62 , the value of expected value function T(f) for the selected f is found, and it is decided whether or not it includes regularly patterned defects; and if it is decided that regularly patterned defects are included (process  63 ), then in regular patterned defect mix rate function calculation process  73 , regular pattern defect mix rate function MR(f) is calculated from number of bits n, factor f, and the value of estimated value function T(f). If it is decided that it does not contain regularly patterned defects (process  65 ), in regular patterned defect mix rate function calculation process  66 , regular patterned defect mix rate function MR(f) is assumed to be zero; and finally, in MR(f) calculation completion confirmation process  67 , it is confirmed whether or not MR(f) has been found for every f, and if there are factors for which MR(f) has not been found, then it returns to process  61.

Description:
BACKGROUND OF THE INVENTION  
         [0001]    The present invention relates to a failure analysis apparatus for semiconductor application devices, and in particular relates to a failure analysis apparatus for memory LSI.  
           [0002]    Description of the Related Art  
           [0003]    Conventional LSI failure analysis devices will now be explained. As an example of a memory LSI failure analysis device for improving yield and making known the cause of failures, Japanese Patent Application Laid-open No. 07-072206 provides an expert system in which the expertise of a process engineer, a circuit engineer and a layout engineer are all incorporated in a personal computer.  
           [0004]    In addition, Japanese Patent Application Laid-Open Nos. 2000-200814, 2000-321333, 2000-21195, and 2000-187062 provide device structures, which distinguish whether or not a defect is a result of the design by analyzing the types of divisors of the intervals between defective bits and their frequencies.  
           [0005]    Moreover, in Japanese Patent Application Laid-open No. 11-306793, as a defective bits map analysis method, a defect analysis method is provided, in which defective bits map is subjected to wavelet transform, an X-direction/high pass filtered and Y-direction/low pass filtered information XH (i, j) is summed for each i in the Y-direction to create a histogram, defective addresses X are obtained from those addresses i not having their summed value equal to zero, and the number of defects is then obtained from the respective absolute values of those summed values and outputted.  
           [0006]    Since the speed at which the large capacity and the high density of memory LSI tend to grow is ever increasing, it has now become necessary to respond to the defect analysis of dynamic random access memory (DRAM) having 256 megabits, 1 gigabit or more. Added to that is the fact that wafer size has also been increasing so that in the future they will most certainly have large diameters of 300 mm, and in those cases, the number of chips and the number of defects to be analyzed together with increase drastically.  
           [0007]    In order to solve the above-mentioned problems, high-speed techniques have been provided, for example, in Japanese Patent Application Laid-Open No. 2000-21195, which discloses defect analysis using distributed processing; Japanese Patent Application Laid-Open No. 2000-187062, which discloses defect analysis using regional divisions; and Japanese Patent Application Laid-Open No. 2000-321333, which discloses defect analysis algorithms.  
           [0008]    Not only does an increase in the number of defects increase the load to be analyzed, but it also means that the amount of analyzed results output as a result of LSI defect analysis also becomes enormous, thus making it difficult to manually judge the cause of the defects.  
           [0009]    For example, considering cases where it is used as a monitoring device in order to estimate the reduction in yield from the production line, since it is impossible to fully confirm the analyzed results manually, the defect analysis device itself needs to interpret the analyzed results and give off an alarm.  
           [0010]    Accordingly, in future LSI defect analysis devices, not only does the analytical ability need to be improved, but also they must have the ability to automatically interpret the analyzed results.  
           [0011]    There are no conventional LSI defect analysis devices equipped with the ability to automatically interpret the analyzed results.  
           [0012]    Furthermore, the device disclosed in Japanese Patent Application Laid-Open No. 2000-200814 analyzes the address difference between defective bits and judges whether or not there is a pattern to the defect distribution; however, the result output by this device is the value of an expected value function. In cases where an LSI defect analysis device is used as the monitoring device mentioned above, it can be thought of as only being useful in setting off an alarm when the rate of patterned defects exceeds a preset threshold value. Nevertheless, the expected value function output from the device disclosed in Japanese Patent Application Laid-Open No. 2000-200814 has qualities whereby even if the rate of patterned defects remains the same, its value will increase according to factor f; therefore, it is not suitable for such usage.  
           [0013]    Accordingly, while taking into consideration the above problems and limitations, the objective of the present invention is to provide a device, method and storage medium, stored with software programs, which, when a memory LSI defect analysis device is used as a monitoring device to estimate reductions in yield, can shorten the time needed for full manual interpretation of the obtained results, by automatically interpreting the analyzed results obtained, and calculating the period of distribution patterns and the mix rate of regular patterned defects. Further objectives, characteristics, and benefits of the present invention will become apparent to those skilled in the art from the description of the following preferred embodiments.  
         SUMMARY OF THE INVENTION  
         [0014]    The present invention, which achieves above-mentioned objectives, comprises: a testing means, which tests LSI; a testing means, which tests said memory LSI; a data readout means, which reads data that is output form said testing means and holds it in memory region of an analytical computer; an address difference calculation means, which calculates address differences between two defect data; an address difference histogram production means, which produces an address difference histogram based on said address difference; an expected value function calculation means, which calculates expected value function T(f) for function f of an address difference based on said address difference histogram; and a regular patterned defect mix rate calculation means, which calculates mix rate of regular patterned defects included in a defect distribution from said expected value function.  
           [0015]    The present invention, further comprises a regular patterned defect mix rate function calculation means, which finds a regular patterned defect mix rate function for every factor f.  
           [0016]    The present invention, further comprises a means of compensating a baseline; and a regular patterned defect mix rate function calculation means, which finds a regular patterned defect mix rate function based on said baseline.  
           [0017]    According to another invention, a memory LSI defect analysis method is provided which comprises a first step of testing a memory LSI; a second step of reading data associated with defect bits which is obtained from results of testing in said first step, and holding the data in a memory section of an analytical calculator; a third step of calculating address differences between two pieces of defect data; a fourth step of producing an address difference histogram based on said address differences; a fifth step of calculating expected value function T(f) for factor f of an address difference based on said address difference histogram; and a sixth step of calculating a regular patterned defect mix rate of a defect distribution from said expected value function.  
           [0018]    According to still another invention, a computer program product comprising a computer usable medium having computer readable program code means embodied in said medium for performing electrical testing on memory LSI to be analyzed, by a computer having a memory LSI defect analysis apparatus that comprises a memory test system, said product comprising; a first computer readable program code means for reading data that is obtained as a result of said testing and holding the data in memory region of an analytic calculator; a second computer readable program code means for calculating address differences between two pieces of defect data; a third computer readable program code means for producing an address difference histogram based on said address difference; a fourth computer readable program code means for calculating expected value function T(f) for factor f of an address difference based on said address difference histogram; and a fifth computer readable program code means for calculating mix rate of regular patterned defects included in a defect distribution from said expected value function. 
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0019]    The above-mentioned and other objects, features and advantages of this invention will become more apparent by referencing the following detailed description of the invention taken in conjunction with the accompanying drawings, wherein:  
         [0020]    [0020]FIG. 1 is a block diagram illustrating the constitution of the first embodiment of the present invention;  
         [0021]    [0021]FIG. 2 is a block diagram illustrating the constitution of the second embodiment of the present invention;  
         [0022]    [0022]FIG. 3 is a block diagram illustrating the constitution of the third embodiment of the present invention;  
         [0023]    [0023]FIG. 4 is a flowchart illustrating the constitution of the first embodiment of the present invention;  
         [0024]    [0024]FIG. 5 is a flowchart showing the process flow of regular patterned defect mix rate calculation means  16  according to the first embodiment of the present invention;  
         [0025]    [0025]FIG. 6 is a flowchart showing the process flow of regular patterned defect mix rate function calculation means  21  according to the second embodiment of the present invention;  
         [0026]    [0026]FIG. 7 is a flowchart showing the process flow of regular patterned defect mix rate calculation means  32  according to the third embodiment of the present invention;  
         [0027]    [0027]FIG. 8 is a diagram for describing the first embodiment of the present invention, and is a diagram showing an example of defect data including regular patterned defects at period  10 ;  
         [0028]    [0028]FIG. 9 shows a graph of expected value function T(f) for defect data A in FIG. 8;  
         [0029]    [0029]FIG. 10 illustrates a graph of expected value function T(f) for defect data B in FIG. 8;  
         [0030]    [0030]FIG. 11 is a graph of expected value function T(f) for defect data C in FIG. 8;  
         [0031]    [0031]FIG. 12 shows a graph of expected value function T(f) for defect data D in FIG. 8;  
         [0032]    [0032]FIG. 13 is a graph of expected value function T(f) for defect data E in FIG. 8;  
         [0033]    [0033]FIG. 14 is a diagram showing a graph in which the relationship between the ratio of (T max −1)/(f max −1) and the ratio of the number of regular patterned defects over the total number of defects is shown;  
         [0034]    [0034]FIG. 15 is a diagram showing a graph in which the relationship between the ratio of (T max −1)/(f max −1) and the square of the ratio of the number of regular patterned defects over the total number of defects is shown;  
         [0035]    [0035]FIGS. 16A and 16B are a diagram describing baseline compensation;  
         [0036]    [0036]FIG. 17 shows an example of a graph of expected value function T(f); and  
         [0037]    [0037]FIG. 18 is a diagram showing a graph of regular patterned defect mix rate function MR(f) corresponding to the defect data used in FIG. 17.  
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS  
       [0038]    The preferred embodiments of the present invention will now be described. In a first embodiment of the present invention, which will be described in detail later, a memory LSI defect analysis device, which comprises a memory test system that performs electrical testing of memory LSI, is comprised of testing means  11  (see FIG. 1), which tests memory LSI; data readout means  12 , which reads out the data of defect bits that have been read out from the testing means and records it; address difference calculation means  13 , which calculates the address difference between two defect data; address difference histogram preparation means  14 , which creates a histogram of address differences based on the address differences; expected value function calculation means  15 , which calculates expected value function T(f) for factor f based on the address difference histogram; and, regular patterned defect mix rate calculation means  16 , which calculates the mix rate of regular patterned defects that are included in the defect distribution by the expected value function.  
         [0039]    In more detail, regular patterned defect mix rate calculation means  16 , as shown in FIG. 5,  
         [0040]    (a) finds maximum value T max  of expected value function T(f) for factor f,  
         [0041]    (b) finds value f max  of f at the time of maximum value T max ,  
         [0042]    (c) decides if maximum value T max  is larger than 1,  
         [0043]    (d) decides if larger than 1 then “regular patterned defect”, or if 1 or smaller then “irregular patterned defect”, and  
         [0044]    (e) when regular patterned defect, calculates mix rate with the following formula:  
         regular patterned defect mix rate=         mix                 rate     =             T                 max     -   1         f                 max     -   1         .                             
         [0045]    In another embodiment of the present invention, instead of regular patterned defect mix rate calculation means  16 , regular patterned defect mix rate function calculation means  21  may be provided, which calculates regular patterned defect mix rate function MR(f) from expected value function T(f).  
         [0046]    As shown in FIG. 6, regular patterned defect mix rate function calculation means  21   
         [0047]    (a) selects factor f and finds the value of expected value function T(f) for the factor f,  
         [0048]    (b) decides if T(f) is larger than 1 or not,  
         [0049]    (c) if T(f) is larger than 1, then decides the defect distribution has a regular pattern at period f and finds the value of regular patterned defect mix rate function MR(f) using the following formula:  
           MR        (   f   )       =           T        (   f   )       -   1       f   -   1           ,                         
 
         [0050]    (d) if T(f) is 1 or less, then decides the defect distribution does not have a regular pattern at period f and makes  
         regular patterned defect mix rate function  MR ( f )=0, and  
         [0051]    (f) confirms whether or not regular patterned defect mix rate function MR(f) has been found for all factors f, and if finished then ends one string of steps, if not returns to step (a).  
         [0052]    In another embodiment of the present invention, baseline compensation means  31 , which allows calculation of a more accurate mix rate, and regular patterned defect mix rate function calculation means  32 , which calculates a regular patterned defect mix rate function corresponding to a baseline compensated by the baseline compensation means  31  are provided.  
         [0053]    As shown in FIG. 7, regular patterned defect mix rate function calculation means  32   
         [0054]    (a) finds total number of defects n using defect calculation processing,  
         [0055]    (b) selects factor f,  
         [0056]    (c) finds the value of expected value function T(f) for that f, and decides whether or not the value is (n−f)/(n−1) or greater and that the number of defects n is equal to or greater than factor f,  
         [0057]    (d) when the conditions of T(f)&gt;(n−f)/(n−1) and n&gt;f have been met, decides defect distribution has a regular pattern at period f, and finds regular patterned defect mix rate function MR(f) using the following formula  
           MR        (   f   )       =             (     n   -   1     )          T        (   f   )         -   n   +   f       n        (     f   -   1     )             ,                         
 
         [0058]    (e) when the above-mentioned conditions are not met, decides defect distribution does not have a regular pattern at period f, and makes  
         regular patterned defect mix rate function  MR ( f )=0,  
         [0059]    (f) confirms whether or not regular patterned defect mix rate function MR(f) has been found for all of the factors f, and if finished then ends one string of steps, if not returns to step (b).  
         [0060]    In the present invention, data readout means  12 , address difference calculation means  13 , histogram preparation means  14 , expected value function calculation means  15 , regular patterned defect mix rate function calculation means  16 , regular patterned defect mix rate function calculation means  21 , baseline compensation means  31 , and regular patterned defect mix rate function calculation means  32  are able to execute their processes and functions using the program run in the computer which comprises the memory defect analysis device. In this case, in the computer that contains a storage medium readout device, the storage medium (e.g. floppy disk (FD), hard disk drive (HDD), magnetic tape (MT), compact disk ROM (CD-ROM), digital versatile disk (DVD), semiconductor memory, etc.), which records the program, reads it out to the computer via the storage medium readout device. Alternatively, the program can be downloaded to the computer via a communication means, then be loaded into the main memory and be executed in order to put into operation the memory defect analysis device of the present invention.  
         [0061]    The embodiments of the present invention will now be described below while referencing the Figures.  
         [0062]    First Embodiment:  
         [0063]    [0063]FIG. 1 is a diagram showing the structure of a first embodiment of the memory LSI defect analysis device of the present invention. FIG. 4 is a flowchart showing the process flow according to the first embodiment of the present invention.  
         [0064]    As shown in FIGS. 1 and 4, in the first embodiment of the present invention, testing means  11  performs electrical testing on the memory LSI to be tested, the results are then output as bitmap data to a memory device. The bitmap data is read out by the data readout means  12 , and the coordinate data of each defect bit is held in the memory of analytical computer (Step  41 ).  
         [0065]    In the address difference calculation means  13 , two defect bits, a and b, of the defect bits that have been read are selected, and the difference of their addresses d(a, b) is calculated (Step  42 ). The method of calculating this address difference varies depending on the type of analysis; however, when X address analysis is performed, the address difference can be shown as the absolute value of the difference of the X coordinates of defect bits a and b. In other words, if the X coordinates of a and b are given as a(x) and b(x), respectively, then  
           d ( a, b )=| a ( x )− b ( x )|  (1).  
         [0066]    On the other hand, when Y address analysis is performed, the address difference is calculated to be the absolute value of the difference of the Y coordinates of defect bits a and b. That is, if the Y coordinates of a and b are given as a(y) and b(y), respectively, then  
           d ( a, b )=| a ( y )− b ( y )|  (2).  
         [0067]    In the histogram preparation means  14 , summation is performed in order to prepare histogram H(d) for address difference d. Specifically, for defect bits a and b, H(d(a, b)) is incremented by 1 (Step  43 ).  
         [0068]    A completed address difference histogram H(d) can be obtained by performing the operations of address difference calculation means  13  and histogram preparation means  14  for each respective pair of defect bits (Steps  42  to  44 ).  
         [0069]    In the expected value function calculation means  15 , the calculation of expected value function T(f) is performed based on the address difference histogram H(d) (Step  45 ).  
         [0070]    It should be noted that the expected value function T(f) has been defined as equation (3) shown below (reference: Japanese Patent Application Laid-Open No. 2000-200814).  
               T        (   f   )       =       f                   ∑     m        (   f   )             (     N   -   ux     )               (   3   )                               
 
         [0071]    where,  
         [0072]    Σm(ƒ): the number of combinations of defect bits where address difference has ƒ as a factor,  
         [0073]    N: total number of combinations of defect bits, and  
         [0074]    ux: the number of combinations of defect bits where address difference becomes 0.  
         [0075]    When max(d) denotes the maximum address difference of defect bits, n denotes the number of defect bits, and j denotes the counter, the following equations are given:  
               N   =       n        (     n   -   1     )       2       ,   and               ux   =     H        (   0   )                                   
  ux=H (0)  
         [0076]    Therefore, expected value function T(f) can be found from address difference histogram H(d) using equation ( 4 ) shown below.  
           T ( f )=2 f/{n ( n− 1)−2 H( 0)}×Σ( fj )  (4)  
         [0077]    Here, the summing of ΣH(fj) is performed for j where j=1 to fj≦max(d).  
         [0078]    In the regular patterned defect mix rate function calculation means  16 , the mix rate of regular patterned defects is calculated based on expected value function T(f) (Step  46 ). It is possible to find the relationship between the expected value function T(f) needed for this, and the mix ratio of regular patterned defects in the following manner:  
         [0079]    The mix ratio of regular patterned defects is defined as the ‘ratio of the regular patterned defects to the total number of defects’; however, there exists the following relationship between the expected value function T(f) for defective data, which includes regular patterned defects with period λ, and the mix rate of regular patterned defects.  
         [0080]    [0080]FIGS. 9 through 13 show the results from calculating the expected value function T(f) for defect data A through E, respectively, which are shown in FIG. 8 and include regular patterned defects with period  10 .  
         [0081]    As shown in FIG. 8, data A has a defect distribution where each coordinate of address coordinates 0 to 1000 is a defect bit, and from 1000 to 50000, there is a defect bit every tenth coordinate (e.g., 1010, 1020, . . . 49990, 50000). The total number of defect bits is 5900, and out of that, 4900 are regular patterned defects with period  10 , therefore, the regular patterned defect mix rate is  
         4900/5900×83.05%.  
         [0082]    Data B has a defect distribution where each coordinate of address coordinates 0 to 10000 is a defect bit, and from 10000 to 50000, there is a defect bit every tenth coordinate (e.g. 10010, 10020, . . . , 49990, 50000). The total number of defect bits is 14000, and out of that, 4000 are regular patterned defects with period  10 , therefore, the regular patterned defect mix rate is  
         4000/14000×28.57%.  
         [0083]    As for Data C, it has a defect distribution where each coordinate of address coordinates 0 to 20000 is a defect bit, and from 20000 to 50000, there is a defect bit every tenth coordinate (e.g., 20010, 20020, . . . , 49990, 50000). The total number of defect bits is 23000, and 3000 out of 23000 are regular patterned defects with period  10 . Therefore, the regular patterned defect mix rate is  
         3000/23000×13.04%.  
         [0084]    Data D has a defect distribution where each coordinate of address coordinates 0 to 30000 is a defect bit, and from 30000 to 50000, there is a defect bit every tenth coordinate (e.g., 30010, 30020,. . . , 49990, 50000). The total number of defect bits is 32000, and out of that, 2000 are regular patterned defects with period  10 , therefore, the regular patterned defect mix rate is  
         2000/32000×6.25%.  
         [0085]    Data E has a defect distribution where each coordinate of address coordinates 0 to 40000 is a defect bit, and there is a defect bit every tenth coordinate (e.g., 40010, 40020, . . . , 49990, 50000) from 40000 to 50000. The total number of defect bits is 41000, and out of that, 1000 are regular patterned defects with period  10 . The regular patterned defect mix rate is therefore  
         1000/41000×2.44%.  
         [0086]    At the same time, as shown in the graphs of the respective expected value function T(f) in FIGS. 9 through 13, the height of the respective peaks gradually decreases as follows:  
         [0087]    in FIG. 9, which is the analyzed result of data A, T(10)=7.207,  
         [0088]    in FIG. 10, which is the analyzed result of data B, T(10)=1.734,  
         [0089]    in FIG. 11, which is the analyzed result of data C, T(10)=1.153,  
         [0090]    in FIG. 12, which is the analyzed result of data D, T(10)=1.035, and  
         [0091]    in FIG. 13, which is the analyzed result of data E, T(10)=1.005.  
         [0092]    Now, the expected value function T(f) has the following properties.  
         [0093]    (1) When the defect distribution has only regular patterned defects with period λ, T(f) has its peak at f=kλ (k=1, 2, 3, . . . ), and its maximum value is T(λ)=λ, when f=λ.  
         [0094]    (2) When the defect distribution is completely irregular patterned defects, T(f)=1 for all f.  
         [0095]    It can be regarded that in (1) above, the regular patterned defect mix rate is 100%, and in (2) above, the regular patterned defect mix rate is 0%. Accordingly, the mix rate of regular patterned defects with period f max  can be expressed as the value of  
         (T max −1)/(f max −1),  
         [0096]    where T max  is the maximum value of expected value function T(f), and f max  is the value of factor f when expected value function T(f) has the maximum value T max .  
         [0097]    In the examples of data A through data E, since f max =10, the relationship of the value of (T(10)−1)/(10−1) of each data and the value of the ratio of the number of regular patterned defects over the total number of defects is as shown in the graph in FIG. 14. The graph in FIG. 14 shows the relationship between (T max −1)/(f max −1) and the ratio of the number of regular patterned defects over the total number of defects.  
         [0098]    Note that the data that has been plotted also utilizes 35 pieces of data differing from the respective mix rates in addition to data A through E.  
         [0099]    As it is clear from the graph in FIG. 14, it can be understood that there exists the relationship of a second order equation between the value of (T(10)−1)/(10−1) in each data and the value of the number of regular patterned defects over the total number of defects.  
         [0100]    [0100]FIG. 15 is the graph where the relationship between the value of (T(10)−1)/(10−1) for each data and the squared value of the ratio of the number of regular patterned defects over the total number of defects is plotted. According to the graph in FIG. 15, it can be understood that  
         ( T (10)−1)/(10−1)={(the number of regular patterned defects)/(the total number of defects)} 2   (5)  
         [0101]    Accordingly, using maximum value T max  of expected value function T(f), and value f max  of f at that time, the following formula is obtained:  
               regular patterned defect mix rate     =           T                 max     -   1         f                 max     -   1                 (   6   )                               
 
         [0102]    From the above, the process of regular patterned defect mix rate calculation means  16  can be explained by the flowchart in FIG. 5.  
         [0103]    Specifically, the maximum value T max  of expected value function T(f) is first found in Step  51 , then f max , or the value of factor f when the expected value function T(f) takes its maximum value T max  is obtained in Step  52 .  
         [0104]    Next, whether or not T max  is greater than one is decided in Step  53 ; if it is greater than 1, then ‘regular pattern distribution’ (Step  54 ); if it is not, then ‘irregular pattern distribution’ (Step  55 ).  
         [0105]    When it is regular pattern distribution, the mix rate is found in Step  56  using the following formula:  
               regular patterned defect mix rate     =             T                 max     -   1         f                 max     -   1         .             (   7   )                               
 
         [0106]    Second Embodiment:  
         [0107]    The second embodiment of the memory LSI defect analysis device according to the present invention will be described. FIG. 2 is a diagram showing the structure of the second embodiment of the memory LSI defect analysis device according to the present invention. FIG. 7 is a flowchart showing the process flow of the second embodiment of the present invention.  
         [0108]    As shown in FIG. 2, the memory LSI defect analysis device of the second embodiment comprises storage medium  11 , data readout means  12 , address difference calculation means  13 , histogram production means  14 , and expected value function calculation means  15  as with the first embodiment described above as shown in FIG. 1, and in addition to these, regular patterned defect mix rate function calculation means  21  is provided.  
         [0109]    The second embodiment of the present invention calculates regular patterned defect mix rate function MR(f) from expected value function T(f) in the regular patterned defect mix rate function calculation means  21 .  
         [0110]    In the first embodiment described above, the formula that calculates the regular patterned defect mix rate can be applied not only for f when T(f) is maximum, but for all values of f. Therefore, regular patterned defect mix rate function MR(f) can be thought of as a function of f as with T(f), and it can be defined by the following equation:  
               MR        (   f   )       =           T        (   f   )       -   1       f   -   1                 (   8   )                               
 
         [0111]    As a result, the process of regular patterned defect mix rate function calculation means  21  can be explained using the flowchart shown in FIG. 6.  
         [0112]    In the regular patterned defect mix rate function calculation means  21 , factor f is first selected in Step  61 , the value of expected value function T(f) for that f is then found, and it is decided in Step  62  whether or not that value is greater than 1. If it is greater than 1, then the defect distribution is determined to be regularly patterned at period f (Step  63 ), and regular patterned defect mix rate function MR(f) is found in Step  64  using the following formula:  
         MR        (   f   )       =           T        (   f   )       -   1       f   -   1                               
 
         [0113]    On the other hand, if the value of T(f) is equal to or less than 1, then the defect distribution is determined, in Step  65 , not to be regularly patterned at period f, and regular patterned defect mix rate function MR(f) becomes equal to 0 (Step  66 ).  
         [0114]    Finally, it is confirmed in Step  67  that regular patterned defect mix rate MR(f) has been found for every f, and if finished then one string of steps ends, if not then it returns to Step  61 .  
         [0115]    Third Embodiment:  
         [0116]    The third embodiment of the memory LSI defect analysis device according to the present invention will be described. FIG. 3 is a diagram showing the structure of the third embodiment of the memory LSI defect analysis device according to the present invention. FIG. 7 is a flowchart showing the process flow of the third embodiment of the present invention.  
         [0117]    As shown in FIG. 3, the memory LSI defect analysis device of the third embodiment comprises storage medium  11 , data readout means  12 , address difference calculation means  13 , histogram production means  14 , and expected value function calculation means  15  as with the first embodiment described above as shown in FIG. 1, and in addition to these, baseline compensation means  31  and regular patterned defect mix rate function calculation means  32  are provided.  
         [0118]    In the third embodiment of the present invention, the baseline compensation means  31  performs baseline compensation during calculation of regular patterned defect mix rate function MR(f). The regular patterned defect mix rate function calculation means  32  then calculates, based on that compensation, regular patterned defect mix rate function MR(f).  
         [0119]    First, the baseline compensation means  31  will be described.  
         [0120]    With the equation for the regular patterned defect mix rate in the first and second embodiments described above, the expected value function T(f) of irregularly distributed defects becomes equal to 1 when there is an unlimited number of defects; or in other words, it is assumed that T(f) becomes equal to 1 for all factors f.  
         [0121]    However, since there is inevitably a limit to the number of defect data in the actual object that is analyzed, in order to calculate a more accurate mix rate, baseline compensation in the form of T(f)=1 must be performed.  
         [0122]    In order to perform this, the value of the expected value function for irregularly distributed defects in the case where there was a limit to the number of defects, is first found, and then the mix rate formula may be defined in accordance with it.  
         [0123]    In the following, the value of the expected value function for irregular patterned defects having a limited number of defects is found.  
         [0124]    To begin with, it is first assumed to have an irregular distribution of n-number of defects where a defect bit exists at each coordinate within coordinates 1 to n.  
         [0125]    Accordingly,  
         [0126]    there are n−2 defect pairs having a distance of 2: (1, 3), (2, 4), . . . , (n−2, n), and  
         [0127]    there are n−3 defect pairs having a distance of 3: (1, 4), (2, 5), . . . , (n−3, n),  
         [0128]    and similarly,  
         [0129]    there are n−d defect pairs having a distance of d: (1, d+1), (2, d+2), . . . , (n−d, n), and  
         [0130]    there are n−fi defect pairs having a distance of fi:(1, fi+1), (2, fi+2), . . . , (n−fi, n).  
         [0131]    At this point, if the number of defect pairs having multiples of distance f are counted, the counted value becomes the sum of the number of pairs in terms of distances: f, 2f, 3f, . . . . If it is defined as k=[n/f] (the largest integer not exceeding n/f), then the following equation can be given:  
         the number of pairs each having a distance that is a multiple of  f=Σ   i=1   k ( n−fi )                    f   =       ∑     ι   =   1     k          (     n   -   fi     )                   =         k        (       2      n     -   fk   -   f     )       2     .                   (     9        -        1     )                                 
         [0132]    On the other hand, the total number of defect pairs is equal to the number of combinations  n C 2 =n(n−1)/2. Accordingly, regular patterned defect ratio P(f) is given by,  
                     P        (   f   )       =       2     n        (     n   -   1     )         ×       k        (       2      n     -   fk   -   f     )       2                   =       k   (       2      n     -   fk   -   f         n        (     n   -   1     )                       (     9        -        2     )                               
 
         [0133]    Therefore, expected value function T(f) becomes  
               T        (   f   )       =       fP        (   f   )       =         fk        (       2      n     -   fk   -   f     )         n        (     n   -   1     )         .               (   10   )                               
 
         [0134]    At this point, since, if n is sufficiently large, it is regarded as  
         n=fk  
         [0135]    obtaining  
               T        (   f   )       =         n        (       2      n     -   n   -   f     )         n        (     n   -   1     )         =         n   -   f       n   -   1       .               (   11   )                               
 
         [0136]    As a result, it is possible to compensate the baseline from what is shown in FIG. 16A (before compensation, T(f)=1) into what is shown in FIG. 16B. FIG. 16B shows expected value function T(f) once it has been compensated using equation (11) above, based on the relationship to the number of defects n for factor f.  
         [0137]    In the regular patterned defect mix rate function calculation means  32 , regular patterned defect mix rate function MR(f) is calculated in the manner described below from expected value function T(f) based on the above baseline.  
                 MR        (   f   )       2     =           T        (   f   )       -       n   -   f       n   -   1           f   -       n   -   f       n   -   1           =           (     n   -   1     )          T        (   f   )         -   n   +   f       n        (     f   -   1     )                   (   12   )                               
 
         [0138]    Therefore,  
               MR        (   f   )       =             (     n   -   1     )          T        (   f   )         -   n   +   f       n        (     f   -   1     )                   (   13   )                               
 
         [0139]    However, when  
         T        (   f   )       &lt;       n   -   f       n   -   1                             
 
         [0140]    and n&lt;f, then  
           MR ( f )=0  (14).  
         [0141]    Accordingly, the processes of regular patterned defect mix rate function calculation means  32  can be described using the flowchart shown in FIG. 7.  
         [0142]    Namely, the total number of defects is found in Step  71 , the factor f is then selected in Step  61 . The value of expected value function T(f) for that f is found, and it is decided whether or not that value is equal to or greater than (n−f)/(n−1) as well as whether or not the number of defects n is equal to or greater than factor f in Step  72 .  
         [0143]    If the conditions of Step  72  have been met, then it is decided that the defect distribution is regularly patterned at period f in Step  63 , and then the value of regular patterned defect mix rate function MR(f) is found in Step  73  using the following equation:  
         MR        (   f   )       =               (     n   -   1     )          T        (   f   )         -   n   +   f       n        (     f   -   1     )         .                             
 
         [0144]    On the other hand, if the conditions of Step  72  have not been met, then it is decided in Step  64  that the defect distribution is not regularly patterned at period f, and thereby making regular patterned defect mix rate function MR(f)=0 in Step  66 .  
         [0145]    Next, and finally, it is confirmed in Step  67  whether or not regular patterned defect mix rate function MR(f) has been found for every f, and if finished then one string of steps ends. If not, then it returns to Step  61 .  
         [0146]    Expected value function T(f) has a tendency to increase with factor f, even if the mix rate of regular patterned defects does not change. However, with regular patterned defect mix rate function MR(f), since it does not have such properties, it can also be put into use as a new index for estimating the cause of defects.  
         [0147]    For example, in FIG. 17, expected value factor T(f) is found up to f=128 for certain defect data, and that data is shown as a graph. Looking at FIG. 17, it is noticeable that, since expected value function T(f) is at its maximum value 3.7 at the point where f=114, at first glance it seems as if there are the greatest number of regular patterned defects at period  114 .  
         [0148]    In contrast, in FIG. 18, regular patterned defect mix rate function MR(f) is found for the same data as one used above, and shown as a graph. Looking at FIG. 18, it is clear that there are regular patterned defects at periods (f=14, 19, 31, 38) having higher mix rates than f=114.  
         [0149]    As described above, in the conventional LSI defect analysis device, assuming it is used during the LSI development stage or trial manufacturing stage, it only outputs expected value function T(f) of regular patterned defects. The memory LSI defect analysis device according to the present invention is able to calculate regular patterned defect mix rate function MR(f) from the expected value function T(f) of regular patterned defects. Therefore, it is possible to operate the LSI defect analysis device in cases, for example where it is used in a semiconductor production line as a monitoring device which gives off an alarm when the ratio of regular patterned defects exceeds a certain preset threshold value, making it advantageous.  
         [0150]    Furthermore, the value of expected value function T(f) based on factor f tends to increase as factor f becomes larger even if the mix rate of regular patterned defects does not change. However, since expected value defect mix rate function MR(f) does not have such property, with the present invention, it is easier to spot periods having a high mix rate of regular patterned defects, it can also be put into use as a new index for estimating the cause of defects.