WO2022124053A1 - Information processing device, computer program, information processing method, and simulation information provision method - Google Patents

Information processing device, computer program, information processing method, and simulation information provision method Download PDF

Info

Publication number
WO2022124053A1
WO2022124053A1 PCT/JP2021/042597 JP2021042597W WO2022124053A1 WO 2022124053 A1 WO2022124053 A1 WO 2022124053A1 JP 2021042597 W JP2021042597 W JP 2021042597W WO 2022124053 A1 WO2022124053 A1 WO 2022124053A1
Authority
WO
WIPO (PCT)
Prior art keywords
parameter
information
model
distribution
probability density
Prior art date
Application number
PCT/JP2021/042597
Other languages
French (fr)
Japanese (ja)
Inventor
高史 佐藤
裕貴 塚本
松 辺
Original Assignee
国立大学法人京都大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 国立大学法人京都大学 filed Critical 国立大学法人京都大学
Publication of WO2022124053A1 publication Critical patent/WO2022124053A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing

Definitions

  • the present invention relates to an information processing apparatus, a computer program, an information processing method, and a simulation information providing method.
  • Non-Patent Document 1 points out that the conventional method for extracting model parameters is limited to model parameters that follow a normal distribution, and discloses a statistical extraction method that can also extract model parameters that follow a lognormal distribution. ing.
  • Non-Patent Document 1 cannot handle model parameters that follow an arbitrary distribution or model parameters whose distribution is unknown in the first place. , It is a strong practical constraint when performing a variation consideration simulation.
  • the present invention has been made in view of such circumstances, and an object of the present invention is to provide an information processing apparatus, a computer program, an information processing method, and a simulation information providing method capable of estimating the distribution of model parameters.
  • the information processing apparatus of the present embodiment is a distribution information acquisition unit that acquires distribution information representing the distribution of a first parameter relating to a product. And, based on the model information expressing the first parameter using the second parameter related to the product, the probability density information generation unit that generates the probability density information of the second parameter so that the distribution information is reproduced. And prepare.
  • the distribution of model parameters can be estimated.
  • FIG. 1 is an explanatory diagram showing an example of the configuration of the information processing apparatus 100 of the first embodiment.
  • the information processing device 100 includes a control unit 10, an input unit 11, a distribution information generation unit 12, a probability density information generation unit 13, an output unit 14, a simulation unit 15, a storage unit 16, a display panel 17, and an operation that control the entire device.
  • a unit 18 is provided.
  • the probability density information generation unit 13 includes a sampling unit 131, a conversion unit 132, and a model parameter generation unit 133.
  • the information processing apparatus 100 is difficult to measure (for example, when the measurement takes a long time or only inside the product) based on the distribution of the first parameter (distribution of the performance and characteristics of the product) which is relatively easy to measure. It is possible to estimate the distribution of the second parameter (such as when it is an observed quantity).
  • the first parameter is also called a characteristic parameter, and is relatively dependent on a so-called measuring instrument such as current, voltage, charge, frequency, impedance, length, area, volume, mass, pressure, velocity, acceleration, concentration, density, and luminosity. It can be a physical quantity that can be easily measured.
  • the second parameter is also referred to as a model parameter and is used in a model formula expressing the performance and characteristics of the product, and is a parameter that affects the performance and characteristics of the product due to an error during processing of the product. be able to.
  • Manufactured products are artificially manufactured products and include semiconductor devices, electronic parts, electrical parts, mechanical parts, optical parts, etc. used in products in various fields.
  • semiconductor devices include various devices such as MOSFETs, bipolar transistors, diodes, light emitting devices, light receiving devices, memories, logic ICs, analog ICs, sensors, etc., but in the present specification, MOSFETs (as an example of products) are used as examples of products.
  • MOSFETs as an example of products
  • a description will be given by taking (also referred to as a “device”) as an example.
  • FIG. 2 is an explanatory diagram showing an example of model parameters in MOSFET.
  • a device model also referred to as a “model” representing the voltage / current characteristics of a MOSFET, which is one of the typical circuit elements, can be represented by, for example, (Equation 1).
  • model parameter K is a current gain parameter, which is represented by the product of carrier mobility and insulating film capacity.
  • the model parameter VTH represents the threshold voltage.
  • W represents the channel width and L represents the channel length. The dimensions of the MOSFET are reflected by the channel length L and the channel width W.
  • the device model can be expressed by the equation (2) using the model parameter DELTA. Furthermore, the model parameter CLM (channel length modulation effect coefficient), model parameter MD (mobility deterioration coefficient), and model parameter MDV (mobility deterioration start voltage) are used to determine the bias voltage dependence of the channel length modulation effect and carrier mobility. Considering this, the device model can be expressed by the equation (3). As shown in equation (3), the drain current Id can be expressed as a function of the gate-source voltage Vgs and the drain-source voltage Vds. Equations (1) to (3) are model equations expressing a device model (model) as model information. The model formula is an example and has different expressions depending on the type of the product.
  • model parameters K current gain parameter
  • VTH threshold voltage
  • the input unit 11 can acquire input data. Specifically, the input unit 11 has a function as a measured value acquisition unit, and can acquire the measured value of the characteristic parameter I. Further, the input unit 11 has a function as a distribution information acquisition unit, and can acquire a probability density function Q (I) representing the distribution of characteristic parameters (also referred to as “first parameters”) related to the product. Here, I represents a characteristic parameter.
  • the control unit 10 can determine in advance whether the input unit 11 acquires the measured value of the characteristic parameter I or the probability density function Q (I).
  • the input unit 11 can acquire the measured values of the characteristic parameters I (Id1, Id2) as schematically shown in FIG.
  • the number of bias points is not limited to two, and may be three or more.
  • the distribution information generation unit 12 generates a probability density function Q (I) (also referred to as “distribution information”) representing the distribution of the characteristic parameters I (Id1, Id2) using the measured values acquired via the input unit 11. can do.
  • I probability density function
  • FIG. 4 is a schematic diagram showing an example of the probability density distribution of characteristic parameters.
  • the distribution information generation unit 12 can estimate the density distribution by using, for example, Kernel Density Estimation (KDE).
  • KDE Kernel Density Estimation
  • the kernel function and the optimum bandwidth (smoothing parameter) can be set as appropriate.
  • a standard Gaussian function can be used as the kernel function. This allows you to extrapolate the population data using kernel density estimation given the sample data for a population. That is, the probability density can be estimated even when there is no measured value of the characteristic parameter.
  • FIG. 4 shows the probability density distribution for the characteristic parameter I (Id1, Id2), and for convenience, the difference in density is expressed by the difference in pattern. For example, the central part of the elliptical contour represents the highest probability density.
  • the use of kernel density estimation is not mandatory. For example, it can be replaced by another method of approximating the density of measured values of characteristic parameters.
  • the distribution information generation unit 12 does not have to generate the probability density function Q (I).
  • the probability density information generation unit 13 is a probability of expressing the distribution of the characteristic parameter I based on a model formula (also referred to as “device model” or “model information”) expressing the characteristic parameter I using the model parameter p related to the product.
  • the probability density function P (p) also referred to as “probability density information”
  • the details of the probability density information generation unit 13 will be described later.
  • the simulation unit 15 can evaluate the performance of the simulation target including the product by using the probability density function P (p) of the model parameter p generated by the probability density information generation unit 13 and the model formula.
  • FIG. 1 shows a configuration in which the simulation unit 15 is incorporated in the information processing device 100, the simulation unit 15 may be configured as an external simulator and separated from the information processing device 100. The details of the simulation unit 15 will be described later.
  • the output unit 14 can output the probability density function P (p) of the model parameter p and the simulation result performed by the simulation unit 15 as output data.
  • the storage unit 16 can be composed of a hard disk, a semiconductor memory, or the like, and stores required information.
  • the storage unit 16 can store the input data acquired via the input unit 11. Further, the storage unit 16 can store the processing result of each processing in the information processing apparatus 100 and the information in the middle of processing. Further, the storage unit 16 can store one or a plurality of models (model formulas) in association with the model ID that identifies the device model.
  • the model ID can identify different device models of the same manufacturer, or can identify device models generated by different manufacturers.
  • the control unit 10 has a function as an update unit, and when a new model is provided for the model stored in the storage unit 16 and a new model is acquired via the input unit 11, the acquired new model is used.
  • the model stored in the storage unit 16 can be updated. This allows the latest model to be used.
  • the display panel 17 can be composed of a liquid crystal panel, an organic EL (Electro Luminescence) display, or the like.
  • the operation unit 18 is composed of, for example, a hardware keyboard, a mouse, or the like, and can operate an icon or the like displayed on the display panel 17 or input characters or the like.
  • the operation unit 18 may be configured with a touch panel.
  • the probability density information generation unit 13 includes a sampling unit 131, a conversion unit 132, and a model parameter generation unit 133.
  • the sampling unit 131 samples candidates for the model parameter p using a Markov chain.
  • a Markov chain is a stochastic process in which the future state is determined only by the present state, is irrelevant to the past state, and the possible states are discrete values.
  • the conversion unit 132 converts the sampled candidate into the candidate of the characteristic parameter I by using the model formula.
  • the model parameter generation unit 133 applies the Monte Carlo method to the probability density function Q (I) representing the distribution of the converted characteristic parameter I, and generates the model parameter p from the sampled candidates.
  • the Monte Carlo method is one of the numerical calculation methods, and is a method of finding an approximate solution by repeating trials using random numbers. Hereinafter, a method of generating model parameters will be described.
  • FIG. 5 is a schematic diagram showing the concept of an estimation method of the probability density function P (p) of the model parameter p. Since it is difficult to directly measure the model parameters of the product (for example, MOSFET, etc.) (for example, current gain parameter K, threshold voltage parameter VTH, etc.), it is not possible to know the probability density function that each model parameter follows. .. On the other hand, the characteristic parameters of the product (for example, current I) can be easily measured, and the probability density function Q (I) of the current distribution can be calculated.
  • the probability density information generation unit 13 reproduces the distribution of model parameters using the probability density function Q (I) of the current distribution. That is, the probability density information generation unit 13 approximates the distribution including the correlation of the model parameters by randomly generating a set of model parameters through a process of reproducing the distribution of the model parameters based on the probability density of the current distribution. be able to.
  • the model parameter space and the characteristic parameter space are associated with each other by using the model f (model formula).
  • the model parameter space is composed of two model parameters p1 and p2
  • the characteristic parameter space is a two-dimensional space composed of characteristic parameters Id1 and Id2 (here, currents Id1 and Id2).
  • the model parameter space and the characteristic parameter space are not limited to two dimensions, and the same procedure can be used even if the model parameter space and the characteristic parameter space are higher in dimension.
  • P be the density distribution of the model parameters
  • Q the density distribution of the measurement results of the characteristic parameters.
  • the point p (i) in the model parameter space represents the model parameter of one device and corresponds to the point I (i) in the characteristic parameter space.
  • the points p (i) and the points I (i) in the two spaces are associated with each other by the model f.
  • the points p (i) and p'in the model parameter space represent two different sets of model parameters. That is, these two points represent two different devices. Due to the difference in model parameters, the characteristic parameters (current) I (i) and I'of these two devices are also different.
  • the probability density function Q of the characteristic parameter in the characteristic parameter space is known by measuring the characteristic parameter (current) of many devices.
  • the probability density function P in the model parameter space can be defined so as to reproduce the probability density function Q through the transformation of the model f. This can be achieved by sampling the model parameters via the model f so that the density function of the sample of characteristic parameters is reproduced.
  • specific processing will be described.
  • the control unit 10 determines whether the input data is the probability density function Q (I) of the characteristic parameter I or the measured value of the characteristic parameter I (S11). When the input data is a measured value (measured value in S11), the control unit 10 acquires the measured value of the characteristic parameter I (here, the measured value of the drain electricity of the MOSFET) (S12).
  • the control unit 10 generates the probability density function Q (I) (S13) and performs the process of step S15 described later.
  • the probability density function Q (KDE) is applied to the measured value of the drain current in the J-dimensional characteristic parameter space (also referred to as "current space") stretched by the currents at the J bias points. I) can be calculated. If the number of devices N for current measurement is sufficiently large and the probability density at a certain point in the characteristic parameter space can be directly obtained from the measurement result, the provision of KDE may be omitted. Further, the probability density may be obtained by using an appropriate method other than KDE.
  • the control unit 10 calculates the current characteristic parameter (drain current) I (i) corresponding to p (i) with the current model parameter set using the model formula (S18).
  • the control unit 10 moves the model parameter set p (i) , which is the current sample in the model parameter space, using the proposed distribution q (p'/ p (i) ), and moves the next sample candidate to the next.
  • a candidate model parameter set p' is generated (S19).
  • the proposed distribution for example, a normal distribution can be used in which the current sample p (i) is centered and the variance is ⁇ p, and the sample x from this normal distribution can be obtained from p (i). Let it be the distance traveled to p'.
  • the proposed distribution is not limited to the normal distribution, and various distributions satisfying the detailed equilibrium condition can be used.
  • the detailed equilibrium condition is a condition that the existence probabilities of each are equal even if the states transition according to the transition probabilities.
  • Markov chains are used to sample model parameter candidates.
  • the control unit 10 calculates the characteristic parameter (drain current) I'corresponding to the candidate model parameter set p'using the model formula (S20).
  • P (p') / P (p (i) ) is approximated by Q (I') / Q (I (i) ).
  • the probability density function Q is known by density estimation by KDE or the like, it can be calculated immediately, whereas the probability density function P is P (p') / P (p) because it is difficult to measure. This is because it is difficult to obtain (i) ).
  • the validity of approximating the probability ratio in the model parameter space by the probability ratio in the characteristic parameter space is that the model equation f is injective (or injective) (this is usually satisfied), and the parameter variation is f.
  • the control unit 10 generates a uniform random number R (0 ⁇ R ⁇ 1) (S22).
  • the control unit 10 determines whether or not r> R (S23). When r> R (YES in S23), the control unit 10 sets the next candidate model parameter set p'as the next model parameter set p (i + 1) (p (i + 1) ⁇ p'). (S24), the process of step S27 described later is performed. If r> R (NO in S23), the control unit 10 discards the next candidate model parameter set p'(S25) and sets the current model parameter set p (i) to the next model parameter set p (i ). +1) (p (i + 1) ⁇ p (i) ) (S26). The processing of steps S23 to S26 applies the Monte Carlo method to the probability density function of the characteristic parameter in the characteristic parameter space to generate a model parameter.
  • the control unit 10 determines whether or not the counter is i ⁇ (N-1), and if i ⁇ (N-1) (YES in S27), that is, the processing is completed for all the devices. If not, 1 is added to the counter i (S28), the next model parameter set in step S24 or S26 is set as the current model parameter set in step S17, and the processing after step S17 is repeated. Each time it is repeated, a set of model parameter sets is generated.
  • control unit 10 If i ⁇ (N-1) (NO in S27), the control unit 10 generates a random set of model parameters ⁇ p (1) , p (2) , ..., P (N) ⁇ (S29). ..
  • the set of generated model parameters follows the probability density function P (p).
  • the control unit 10 evaluates the performance of the simulation target using the model formula and the generated model parameters (S30), outputs the probability density function P (p), outputs the simulation result (S31), and ends the process.
  • burn-in which eliminates the correlation caused by the algorithm by rejecting the initial sample and thinning out the sample.
  • burn-in a procedure called burn-in, which eliminates the correlation caused by the algorithm by rejecting the initial sample and thinning out the sample.
  • the model parameters were generated to reproduce the current distribution, but are not limited to the current distribution, for example, to reproduce variations in other measurable electrical properties such as capacitance. Model parameters can be defined and a highly versatile method can be provided.
  • FIG. 8 is an explanatory diagram showing an example of the configuration of the information processing apparatus 100 of the second embodiment.
  • the difference from the first embodiment illustrated in FIG. 1 is that the input data is the measured value of the characteristic parameter I, and the distribution information generation unit 12 for generating the probability density function Q (I) of the characteristic parameter It is unnecessary and includes a fitting unit 19.
  • the model parameters p1 and p2 can be, for example, the current gain parameter K and the threshold voltage parameter VTH, but are not limited thereto. Fitting is effective as a method for extracting the majority of parameters. For some parameters, it is also possible to measure directly without fitting if it is permissible to use destructive measurements. However, in this case, the device itself is destroyed, so it is desirable to guess the parameters by fitting.
  • FIG. 10 is a schematic diagram showing an example of the probability density distribution of model parameters.
  • the probability density information generation unit 13 can estimate the density distribution by using, for example, kernel density estimation (KDE: Kernel Density Estimation).
  • KDE Kernel Density Estimation
  • the kernel function and the optimum bandwidth (smoothing parameter) can be set as appropriate.
  • a standard Gaussian function can be used as the kernel function. This allows you to extrapolate the population data using kernel density estimation given the sample data for a population. That is, the probability density can be estimated even when there is no calculated value of the model parameter.
  • FIG. 10 shows the probability density distribution for the model parameters p1 and p2, and for convenience, the difference in density is expressed by the difference in pattern.
  • the central part of the elliptical contour represents the highest probability density.
  • the probability density function P (p) according to the model parameters can be obtained.
  • the probability density function P (p) can be approximated by another method, for example, a value proportional to the probability density function in pj can be obtained by using the number of parameters around a certain model parameter pj.
  • the control unit 10 acquires the measured value of the characteristic parameter I (here, the measured value of the drain current of the MOSFET) (S41).
  • the control unit 10 fits the model parameters using the model formula (S42).
  • the control unit 10 determines the presence / absence of another device (S43), and if there is another device (YES in S43), repeats the processes after step S41.
  • the control unit 10 When there is no other device (NO in S43), the control unit 10 generates a probability density function P (p) of the model parameter (S44).
  • the control unit 10 sets the initial value of the model parameter set p (1) (S45).
  • the control unit 10 moves the model parameter set p (i) , which is the current sample in the model parameter space, using the proposed distribution q (p'/ p (i) ), and moves the next sample candidate to the next.
  • the proposed distribution for example, a normal distribution can be used in which the current sample p (i) is centered and the variance is ⁇ p, and the sample x from this normal distribution can be obtained from p (i). Let it be the distance traveled to p'.
  • the proposed distribution is not limited to the normal distribution, and various distributions satisfying the detailed equilibrium condition can be used.
  • Markov chains are used to sample model parameter candidates.
  • the control unit 10 generates a uniform random number R (0 ⁇ R ⁇ 1) (S50).
  • the control unit 10 determines whether or not r> R (S51). When r> R (YES in S51), the control unit 10 sets the next candidate model parameter set p'as the next model parameter set p (i + 1) (p (i + 1) ⁇ p'). (S52), the process of step S55 described later is performed.
  • the control unit 10 discards the next candidate model parameter set p'(S53), and sets the current model parameter set p (i) to the next model parameter set p (i ). +1) (p (i + 1) ⁇ p (i) ) (S54).
  • the control unit 10 determines whether or not the counter is i ⁇ (N-1), and if i ⁇ (N-1) (YES in S55), that is, the processing is completed for all the devices. If not, 1 is added to the counter i (S56), the next model parameter set in step S52 or S54 is set as the current model parameter set in step S47, and the processing after step S47 is repeated. Each time it is repeated, a set of model parameter sets is generated.
  • the control unit 10 If i ⁇ (N-1) (NO in S55), the control unit 10 generates a random set of model parameters ⁇ p (1) , p (2) , ..., P (N) ⁇ (S57). .. The set of generated model parameters follows the probability density function P (p). The control unit 10 performs performance evaluation of the simulation target using the model formula and the generated model parameters (S58), outputs the probability density function P (p), outputs the simulation result (S59), and ends the process.
  • the probability ratio r P (p') / P (p (i) at any two points (p (i) and p') in the model parameter space. ) Is obtained, the probability ratio r can be evaluated while moving p using the proposed distribution, and model parameter candidates can be sampled using a Markov chain that repeats its adoption and rejection.
  • the information processing apparatus 100 of the second embodiment acquires the measured value of the first parameter (characteristic parameter indicating the characteristic of the product) relating to the product, and applies the measured value of the acquired first parameter to the measured value.
  • the distribution information of the second parameter can be obtained by applying the model information expressing the first parameter using the second parameter (model parameter of the product) relating to the product.
  • the information processing apparatus 100 can generate the probability density information (probability density function) of the second parameter based on the distribution information of the second parameter. More specifically, the information processing apparatus 100 samples the candidates of the second parameter using the Markov chain, applies the Monte Carlo method to the distribution information representing the distribution of the second parameter, and second from the sampled candidates. Parameters can be generated. The generated second parameter follows a probability density function.
  • FIG. 13 is a circuit diagram showing an example of an SRAM cell circuit as a simulation target.
  • the SRAM cell circuit shown in FIG. 13 uses the above-mentioned MOSFET and its device model.
  • L1 and L2 are p MOSFETs
  • D1, D2, A1 and A2 are n MOSFETs.
  • Device models polyclonal1 and EtOAc1 representing the characteristics of each of the p MOSFET and n MOSFET used in manufacturing the circuit are prepared, and photoresist1 is used for L1 and L2, and IGMP1 is used for D1, D2, A1 and A2.
  • L1 and L2 are designed using the same MOSFET with the same dimensions.
  • FIG. 14A and 14B are schematic views showing an example of the simulation result of the SRAM cell circuit.
  • FIG. 14 shows the evaluation result of the static noise margin (SNM), which is one of the circuit characteristics of the SRAM cell circuit.
  • SNM represents the amount of voltage fluctuation required for reversing the value stored in the memory cell when the voltage at the nodes N 0 and N 1 fluctuates, and is an index for evaluating the stability of the value stored in the memory cell. It is one.
  • FIG. 14A shows a state of normal operation. That is, the state in which two gaps (indicated to include the largest square whose one side length is SNM 0 and SNM 1 ) exist in two curves called butterfly curves is 0/1. Represents a normal operating state that can be stably stored. In the memory cell design, min (SNM 0 , SNM 1 ) is maximized.
  • FIG. 14B shows a state of malfunction. Due to variations in model parameters, for example, if the threshold voltage is different, the circuit characteristics will also be different. FIG. 14B shows an example in which the circuit does not operate normally due to a large variation in model parameters. In this example, SNM 0 , that is, the gap on the left side of the butterfly curve does not exist. By repeating such a simulation many times in response to variations in model parameters, it is possible to obtain the yield of the SRAM cell circuit using the static noise margin as an index.
  • the evaluation results for the generation of the probability density function of the model parameters by the information processing apparatus 100 of the first embodiment will be described.
  • Model parameters such as thresholds are often difficult to measure directly. Therefore, two methods are used as the evaluation method.
  • the first evaluation method assumes a virtual device and shows how accurately it can be estimated with respect to variations in model parameters of the virtual device.
  • the second evaluation method shows how accurately the current characteristic parameters of the actually measured device can be estimated. First, the first evaluation method will be described.
  • FIG. 15 is a schematic diagram showing variations in model parameters of 500 virtual devices.
  • the vertical axis indicates the threshold voltage VTH
  • the horizontal axis indicates the current gain K.
  • the virtual device assumes a situation where two lots are mixed and the model parameters are divided into two clusters.
  • the above equations (1) to (3) are used as the device model, and the variation of the two model parameters of the threshold voltage VTH and the current gain K is considered among the model parameters.
  • the parameter is generated using a pseudo-random number assuming that it follows a normal distribution with two peaks. Model parameters other than these two are fixed to constant values.
  • the average value, variance, and correlation for each model parameter K and VTH may be appropriately set.
  • FIG. 16A and 16B are schematic views showing an example of current variation, which is a characteristic parameter of 500 virtual devices.
  • the drain current of the device is calculated using the model formula.
  • the voltage-current characteristic data of 500 virtual devices were calculated.
  • FIG. 16A shows the distribution of drain current under two bias conditions (bias voltage).
  • bias voltage bias voltage
  • two points (Vgs, Vds) (10V, 1V) and (18V, 5V) are selected as the bias conditions of the current data constituting the characteristic parameter space. These correspond to Id1 and Id2 in FIG. 16A, respectively.
  • FIG. 16B shows the probability density distribution of the drain currents of 500 virtual devices estimated using KDE with respect to the variation of the drain currents.
  • FIG. 17 is a schematic diagram showing an example of the generation result of the model parameter by the comparative example.
  • FIG. 17 shows a scatter plot of model parameters K and VTH generated by the comparative example, and a histogram for each parameter K and VTH.
  • black circles indicate the true values of the model parameters
  • indicates the model parameters generated by the comparative example. Since the parameters generated by the comparative example assume only one normal distribution, it can be seen that the variation of the given parameters cannot be reproduced and the error is large. This situation is also clear in the histogram for each model parameter. The true distribution represented by the solid line has two peaks, whereas in the comparative example, this is approximated by a normal distribution, so that the most frequent parts of the model parameters are completely different.
  • a method based on a comparative example will be described below.
  • the model parameters are generated by the following procedure. (1) Generate a model formula that expresses the performance of the product by a function of parameters (model parameters) at the time of processing, (2) Fit the measured value (for example, current, etc.) of the product to the model formula. The variation of the parameters in the model formula is expressed by a normal distribution or the like. (3) The realization value of the parameter expressed by the normal distribution or the like is generated as a random number according to the normal distribution. Furthermore, when evaluating the performance of the product, a simulation is performed using the generated realization value and the model formula.
  • model parameters are represented as a normal distribution, but in reality, they often do not follow the normal distribution. In addition, it can be difficult to directly measure variability in model parameters. Since it is assumed that parameters that are difficult to measure and whose distribution is unknown follow a normal distribution, the generated model parameters may not be sufficiently accurate.
  • FIG. 18 is a schematic diagram showing an example of the generation result of the model parameter according to the first embodiment.
  • FIG. 18 shows a scatter plot of the model parameters K and VTH generated by the above-described embodiment, and a histogram for each of the parameters K and VTH.
  • black circles indicate the true values of the model parameters
  • indicates the model parameters generated by the present embodiment.
  • a random number created from a two-dimensional normal distribution having an average of 0 and a variance of (0.2) 2 is used as a random number (proposal distribution) when creating a next sample candidate from the current sample. board. That is, if ⁇ K and ⁇ VTH to N (0, 0.2 2 ), the next sample candidate can be generated from the current sample as follows.
  • FIG. 19A and 19B are schematic views showing an example of variation in current, which is a characteristic parameter of the fitting target.
  • FIG. 19A shows actual measurement data (current variation data) of the current under this bias condition.
  • FIG. 19B shows the probability density of the current obtained using KDE.
  • FIG. 20 is a schematic diagram showing an example of the current variation result according to the comparative example.
  • FIG. 20 is an evaluation in which variations in current characteristics are reproduced by simulation and compared with measured values (actual measurement data). It can be seen that the current simulation using the model parameters generated by the comparative example described in FIG. 17 cannot reproduce the distortion of the current distribution appearing in the measured values because it is assumed that the model parameters follow a normal distribution. ..
  • FIG. 21 is a schematic diagram showing an example of the current variation result according to the first embodiment.
  • the measured value actual measurement data
  • the current characteristic variation simulation using the generated model parameters it can be seen that in the present embodiment, more accurate current simulation is possible without assuming the distribution of model parameters.
  • the same evaluation result as in the first embodiment can be obtained also in the second embodiment.
  • FIG. 22 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the first embodiment.
  • the variation consideration simulation can be performed as follows, as shown in FIG. 22.
  • a manufacturer a company that manufactures circuits (FAB)
  • FAB circuits
  • characteristic parameters for example, voltage / current characteristics
  • the manufacturer determines the device model (model formula) of the device to be manufactured, and applies, for example, kernel density estimation to the measured characteristic parameters to generate a probability density function of the characteristic parameters.
  • the manufacturer provides the user with a device model and a probability density function of the generated characteristic parameter.
  • the measured value of the characteristic parameter may be provided.
  • the user uses the device model provided by the manufacturer and the probability density function of the characteristic parameter to use the method of the present embodiment (that is, the Markov chain in the model parameter space and the characteristic parameter space. Model parameters are generated using the model parameter generation method based on the Monte Carlo method).
  • the user acquires the measured value of the characteristic parameter from the manufacturer, the user can generate the probability density function of the characteristic parameter using the acquired measured value.
  • the user performs a circuit simulation considering variation using the provided device model and the generated model parameters.
  • the manufacturer can perform the following simulation information provision method. That is, the server that realizes the simulation information providing method acquires information about the simulation target (for example, a memory cell circuit) from the terminal on the user side, and determines the characteristics of the product (for example, MOSFET) included in the simulation target. Distribution information (probability density function) representing the distribution of the characteristic parameters to be shown is output to the terminal on the user side, and model information (for example, device model, model formula, etc.) expressing the characteristic parameters using the model parameters of the product is output to the user. It can be output to the terminal on the side.
  • the simulation information is information necessary for executing the simulation, and includes the above-mentioned distribution information (probability density function) and model information.
  • FIG. 23 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the second embodiment.
  • the variation consideration simulation can be performed as follows, as shown in FIG. 23.
  • a manufacturer for example, a company that manufactures a circuit (FAB)
  • FAB circuit
  • the manufacturer determines the model formula of the device to be used, and extracts the model parameters so as to approximate the measured voltage / current characteristics (fitting). This makes it possible to obtain the distribution of model parameters.
  • the manufacturer provides the user with a model formula and a distribution of model parameters extracted using the model formula.
  • the user applies, for example, kernel density estimation to the model parameters to generate a probability density function of the model parameters.
  • the user generates model parameters by using the probability density function of the model parameters and the model parameter generation method based on the Markov chain Monte Carlo method (MCMC) in the model parameter space.
  • MCMC Markov chain Monte Carlo method
  • the generated model parameters follow the same distribution as the model parameter distribution provided by the manufacturer, but the sample points are different.
  • the user performs a circuit simulation considering variation using the provided device model and the generated model parameters.
  • FIG. 24 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the comparative example.
  • the comparative example shows an example in the case where the information processing apparatus 100 of this embodiment is not used.
  • the variation consideration simulation can be performed as follows, as shown in FIG. 24.
  • a manufacturer for example, a company that manufactures a circuit (FAB)
  • FAB circuit
  • MOSFET MOSFET
  • the manufacturer determines the model formula of the device to be used, and extracts the model parameters so as to approximate the measured voltage / current characteristics (fitting).
  • the manufacturer approximates the distribution of the extracted model parameters with a normal distribution, and calculates the mean value and standard deviation. (3) The manufacturer provides the user with a model formula and a distribution of model parameters (mean value and standard deviation when approximated to a normal distribution) extracted using the model formula as a statistical device model. (4) The user performs a circuit simulation using the provided model formula and distribution of model parameters.
  • the variation in the characteristics can be summarized as the variation in the model parameters, and the circuit simulation can be performed by the same procedure as in the case where there is no variation. ing.
  • the second embodiment it is possible to generate a large number of new model parameters that are larger than the number of measurement data regarding the model parameters provided by the manufacturer, and the model formula and the model parameters are substantially separated. It can be performed.
  • FIG. 25 is an explanatory diagram showing another example of the configuration of the information processing apparatus of the first embodiment and the second embodiment.
  • reference numeral 200 is a computer.
  • the computer 200 includes a control unit 30, an input unit 40, an output unit 50, an external I / F (interface) unit 60, a display panel 70, an operation unit 80, and the like.
  • the control unit 30 includes a CPU 31, a ROM 32, a RAM 33, an I / F (interface) 34, and the like.
  • the input unit 40 acquires input data.
  • the output unit 50 outputs output data.
  • the I / F 34 has an interface function between the control unit 30, the input unit 40, the output unit 50, the external I / F unit 60, the display panel 70, and the operation unit 80, respectively.
  • the external I / F unit 60 can read the computer program from the recording medium M (for example, a medium such as a DVD) on which the computer program is recorded.
  • the recording medium M for example, a medium such as a DVD
  • the computer program recorded on the recording medium M is not limited to the one recorded on the freely portable medium, and the computer program transmitted via the Internet or another communication line is also included. Can be included.
  • the computer also includes a single computer equipped with a plurality of processors or a computer system composed of a plurality of computers connected via a communication network.
  • the distribution of the model parameters used in the model is accurately estimated from the distribution of characteristic parameters such as the performance of the product, which is relatively easy to measure. be able to. Since the distribution of model parameters can be estimated, more accurate variation-considered simulation can be performed using model parameters and model formulas. Further, based on the performance distribution of the product, it is possible to evaluate the performance of the product without explicitly giving the distribution of model parameters to the simulator.
  • Degree of freedom of distribution representing variation Normally, a normal distribution is assumed because of the ease of random number generation and statistical handling. However, the distribution that each model parameter follows may not always be clear and often does not follow a normal distribution. According to this embodiment, an arbitrary distribution can be expressed as the distribution of model parameters.
  • Handling of correlations between model parameters Statistical correlations are often seen between model parameters, but only linear correlations represented by Pearson's product moment correlation coefficient can be handled. According to this embodiment, any correlation between model parameters can be expressed.
  • Statistical device models may contain (often large) errors if the distribution of model parameters extracted by fitting is distorted or if the current model cannot accurately represent the measured current characteristics. .. In this embodiment, it is not necessary to use the model parameters extracted by the fitting.
  • the variation is expressed as the distribution of model parameters, it is necessary to obtain the distribution of model parameters for each device model. If the circuit simulator used by the user does not support the device model provided by the manufacturer, etc., or if it becomes necessary to perform a simulation using a different device model, the model parameters are reset according to the model formula. It is necessary to perform extraction and obtain its distribution. According to this embodiment, the expression of variation can be made independent of the device model.
  • the model parameter generation method of the present embodiment can be, for example, as follows.
  • p represents a set of second parameters that are relatively difficult to measure.
  • the second parameter is also referred to as a product model parameter.
  • I represents a set of first parameters that are easy to measure.
  • the first parameter is also referred to as a characteristic parameter indicating the characteristics of the product.
  • the other one is a characteristic parameter space, which is a space in which characteristic parameters I such as electric current are distributed.
  • the model parameter generation method of the present embodiment generates a set of model parameters in the model parameter space by repeating a process of alternately performing a Markov chain and a Monte Carlo method in two spaces via a model formula.
  • MCMC Markov chain Monte Carlo method
  • the density distribution of the random variable x is f (x)
  • a set of random variables is generated based on the probability ratio of f (x) of the density distribution, for example, ⁇ f (x t + 1 ) / f (x t ) ⁇ . It's completely different.
  • MOSFET which is one element of a semiconductor device
  • I f (p).
  • the information processing apparatus of the present embodiment has a distribution information acquisition unit that acquires distribution information that represents the distribution of the first parameter related to the product, and model information that expresses the first parameter using the second parameter related to the product. Based on the above, the probability density information generation unit that generates the probability density information of the second parameter is provided so that the distribution information is reproduced.
  • the computer program of the present embodiment acquires distribution information representing the distribution of the first parameter with respect to the product from the computer, and uses the second parameter with respect to the product to represent the first parameter based on the model information. , The process of generating the probability density information of the second parameter so that the distribution information is reproduced is executed.
  • distribution information representing the distribution of the first parameter related to the product is acquired, and the first parameter is expressed by using the second parameter related to the product, based on the model information. Probability density information of the second parameter is generated so that the distribution information is reproduced.
  • the simulation information providing method of the present embodiment acquires information about a simulation target, outputs distribution information indicating a distribution of characteristic parameters indicating the characteristics of the product included in the simulation target, and outputs model parameters of the product.
  • the model information expressing the characteristic parameter is output by using.
  • the information processing apparatus of the present embodiment includes an evaluation unit that evaluates the performance of a simulation target including the product by using the probability density information of the second parameter generated by the probability density information generation unit and the model information. ..
  • the information processing apparatus of the present embodiment has a measurement value acquisition unit that acquires the measurement value of the first parameter and a distribution information generation unit that generates distribution information representing the distribution of the first parameter using the acquired measurement value. And prepare.
  • the information processing device of the present embodiment includes a model information acquisition unit for acquiring the model information and a storage unit for storing the acquired model information.
  • the information processing apparatus of the present embodiment includes an update unit that updates the model information stored in the storage unit with the acquired model information.
  • the probability density information generation unit uses a Markov chain to sample the second parameter candidate, and the model information to sample the sampled candidate first. It is provided with a conversion unit for converting into parameter candidates and a parameter generation unit for generating a second parameter from sampled candidates by applying the Monte Carlo method to distribution information representing the converted distribution of the first parameter.
  • the first parameter includes a characteristic parameter indicating the characteristic of the product
  • the second parameter includes a model parameter of the product

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Primary Health Care (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

Provided are an information processing device, a computer program, an information processing method, and a simulation information provision method with which it is possible to estimate the distribution of model parameters. This information processing device comprises: a distribution information acquisition unit that acquires distribution information representing the distribution of a first parameter related to manufactured goods; and a probability density information generation unit that, on the basis of model information that represents the first parameter using a second parameter related to the manufactured goods, generates probability density information of the second parameter so that the distribution information can be reproduced.

Description

情報処理装置、コンピュータプログラム、情報処理方法及びシミュレーション情報提供方法Information processing equipment, computer programs, information processing methods and simulation information provision methods
 本発明は、情報処理装置、コンピュータプログラム、情報処理方法及びシミュレーション情報提供方法に関する。 The present invention relates to an information processing apparatus, a computer program, an information processing method, and a simulation information providing method.
 人工物の製造には加工限界があり、加工限界に起因するばらつき(誤差)を回避することができない。このようなばらつきは、製造物や製造物を組み込んだ製品の特性に影響を及ぼす。このため、モデルパラメータを含むモデルを生成し、生成したモデルを用いて、ばらつきの影響を事前に評価するばらつき考慮シミュレーションが広く実施されている。 There is a processing limit in the production of artificial objects, and it is not possible to avoid variations (errors) due to the processing limit. Such variability affects the characteristics of the product and the products in which it is incorporated. For this reason, a variation consideration simulation that generates a model including model parameters and evaluates the influence of variation in advance using the generated model is widely performed.
 非特許文献1には、従来のモデルパラメータの抽出方法が、正規分布に従うモデルパラメータに限定されていることを指摘し、対数正規分布に従うモデルパラメータも抽出することができる統計的抽出方法が開示されている。 Non-Patent Document 1 points out that the conventional method for extracting model parameters is limited to model parameters that follow a normal distribution, and discloses a statistical extraction method that can also extract model parameters that follow a lognormal distribution. ing.
 しかし、非特許文献1の方法は、モデルパラメータの分布が正規分布や対数正規分布に限定されているため、任意の分布に従うモデルパラメータ、あるいは、そもそも分布が分からないモデルパラメータを扱うことができないため、ばらつき考慮シミュレーションを行う場合の実用上の強い制約となっている。 However, since the distribution of model parameters is limited to a normal distribution or a lognormal distribution, the method of Non-Patent Document 1 cannot handle model parameters that follow an arbitrary distribution or model parameters whose distribution is unknown in the first place. , It is a strong practical constraint when performing a variation consideration simulation.
 本発明は斯かる事情に鑑みてなされたものであり、モデルパラメータの分布を推定することができる情報処理装置、コンピュータプログラム、情報処理方法及びシミュレーション情報提供方法を提供することを目的とする。 The present invention has been made in view of such circumstances, and an object of the present invention is to provide an information processing apparatus, a computer program, an information processing method, and a simulation information providing method capable of estimating the distribution of model parameters.
 本願は上記課題を解決する手段を複数含んでいるが、その一例を挙げるならば、本実施形態の情報処理装置は、製造物に関する第1パラメータの分布を表す分布情報を取得する分布情報取得部と、前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する確率密度情報生成部とを備える。 The present application includes a plurality of means for solving the above problems, and to give an example thereof, the information processing apparatus of the present embodiment is a distribution information acquisition unit that acquires distribution information representing the distribution of a first parameter relating to a product. And, based on the model information expressing the first parameter using the second parameter related to the product, the probability density information generation unit that generates the probability density information of the second parameter so that the distribution information is reproduced. And prepare.
 本発明によれば、モデルパラメータの分布を推定することができる。 According to the present invention, the distribution of model parameters can be estimated.
第1実施形態の情報処理装置の構成の一例を示す説明図である。It is explanatory drawing which shows an example of the structure of the information processing apparatus of 1st Embodiment. MOSFETにおけるモデルパラメータの一例を示す説明図である。It is explanatory drawing which shows an example of the model parameter in MOSFET. 特性パラメータの計測値の一例を示す模式図である。It is a schematic diagram which shows an example of the measured value of a characteristic parameter. 特性パラメータの確率密度分布の一例を示す模式図である。It is a schematic diagram which shows an example of the probability density distribution of a characteristic parameter. モデルパラメータpの確率密度関数P(p)の推定方法の概念を示す模式図である。It is a schematic diagram which shows the concept of the estimation method of the probability density function P (p) of a model parameter p. 第1実施形態の情報処理装置の処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the processing procedure of the information processing apparatus of 1st Embodiment. 第1実施形態の情報処理装置の処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the processing procedure of the information processing apparatus of 1st Embodiment. 第2実施形態の情報処理装置の構成の一例を示す説明図である。It is explanatory drawing which shows an example of the structure of the information processing apparatus of 2nd Embodiment. 計測された特性パラメータのモデル式を用いたフィッティングの一例を示す模式図である。It is a schematic diagram which shows an example of fitting using the model formula of the measured characteristic parameter. モデルパラメータの確率密度分布の一例を示す模式図である。It is a schematic diagram which shows an example of the probability density distribution of a model parameter. 第2実施形態の情報処理装置の処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the processing procedure of the information processing apparatus of 2nd Embodiment. 第2実施形態の情報処理装置の処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the processing procedure of the information processing apparatus of 2nd Embodiment. シミュレーション対象としてのSRAMセル回路の一例を示す回路図である。It is a circuit diagram which shows an example of the SRAM cell circuit as a simulation target. SRAMセル回路のシミュレーション結果の一例を示す模式図である。It is a schematic diagram which shows an example of the simulation result of the SRAM cell circuit. SRAMセル回路のシミュレーション結果の一例を示す模式図である。It is a schematic diagram which shows an example of the simulation result of the SRAM cell circuit. 500個の仮想デバイスのモデルパラメータのばらつきを示す模式図である。It is a schematic diagram which shows the variation of the model parameter of 500 virtual devices. 500個の仮想デバイスの特性パラメータである電流のばらつきの一例を示す模式図である。It is a schematic diagram which shows an example of the variation of the current which is a characteristic parameter of 500 virtual devices. 500個の仮想デバイスの特性パラメータである電流のばらつきの一例を示す模式図である。It is a schematic diagram which shows an example of the variation of the current which is a characteristic parameter of 500 virtual devices. 比較例によるモデルパラメータの生成結果の一例を示す模式図である。It is a schematic diagram which shows an example of the generation result of a model parameter by a comparative example. 第1実施形態によるモデルパラメータの生成結果の一例を示す模式図である。It is a schematic diagram which shows an example of the generation result of the model parameter by 1st Embodiment. フィッティング対象の特性パラメータである電流のばらつきの一例を示す模式図である。It is a schematic diagram which shows an example of the variation of the current which is a characteristic parameter of a fitting target. フィッティング対象の特性パラメータである電流のばらつきの一例を示す模式図である。It is a schematic diagram which shows an example of the variation of the current which is a characteristic parameter of a fitting target. 比較例による電流のばらつき結果の一例を示す模式図である。It is a schematic diagram which shows an example of the variation result of the current by the comparative example. 第1実施形態による電流のばらつき結果の一例を示す模式図である。It is a schematic diagram which shows an example of the variation result of the current by 1st Embodiment. 第1実施形態の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。It is a schematic diagram which shows an example of the use case of the variation consideration simulation in the case of 1st Embodiment. 第2実施形態の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。It is a schematic diagram which shows an example of the use case of the variation consideration simulation in the case of 2nd Embodiment. 比較例の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。It is a schematic diagram which shows an example of the use case of the variation consideration simulation in the case of a comparative example. 第1実施形態及び第2実施形態の情報処理装置の構成の他の例を示す説明図である。It is explanatory drawing which shows the other example of the structure of the information processing apparatus of 1st Embodiment and 2nd Embodiment.
(第1実施形態)
 以下、本発明をその実施の形態を示す図面に基づいて説明する。図1は第1実施形態の情報処理装置100の構成の一例を示す説明図である。情報処理装置100は、装置全体を制御する制御部10、入力部11、分布情報生成部12、確率密度情報生成部13、出力部14、シミュレーション部15、記憶部16、表示パネル17、及び操作部18を備える。確率密度情報生成部13は、サンプリング部131、変換部132、及びモデルパラメータ生成部133を備える。
(First Embodiment)
Hereinafter, the present invention will be described with reference to the drawings showing the embodiments thereof. FIG. 1 is an explanatory diagram showing an example of the configuration of the information processing apparatus 100 of the first embodiment. The information processing device 100 includes a control unit 10, an input unit 11, a distribution information generation unit 12, a probability density information generation unit 13, an output unit 14, a simulation unit 15, a storage unit 16, a display panel 17, and an operation that control the entire device. A unit 18 is provided. The probability density information generation unit 13 includes a sampling unit 131, a conversion unit 132, and a model parameter generation unit 133.
 出願人は、製造物の加工時のばらつきをモデル化する際に必要となるモデルパラメータは、どのような分布に従うか不明となる場合が多く、モデルパラメータの分布の測定と比較すると、製造物の性能や特性の分布の測定は容易であるという点に着眼した。情報処理装置100は、測定が比較的容易な第1パラメータの分布(製造物の性能や特性の分布)に基づき、測定が困難(例えば、測定に時間がかかる場合、または製造物の内部でのみ観測される量である場合など)な第2パラメータの分布を推定することができる。第1パラメータは、特性パラメータとも称し、例えば、電流、電圧、電荷、周波数、インピーダンス、長さ、面積、体積、質量、圧力、速度、加速度、濃度、密度、光度など、いわゆる計測器によって比較的容易に計測することができる物理量とすることができる。一方、第2パラメータは、モデルパラメータとも称し、製造物の性能や特性を表現するモデル式の中において使用され、製造物の加工時の誤差によって製造物の性能や特性に影響を与えるパラメータとすることができる。 Applicants are often uncertain as to what distribution the model parameters required to model the variability of the product during processing will follow, and when compared to measuring the distribution of the model parameters, the product We focused on the fact that it is easy to measure the distribution of performance and characteristics. The information processing apparatus 100 is difficult to measure (for example, when the measurement takes a long time or only inside the product) based on the distribution of the first parameter (distribution of the performance and characteristics of the product) which is relatively easy to measure. It is possible to estimate the distribution of the second parameter (such as when it is an observed quantity). The first parameter is also called a characteristic parameter, and is relatively dependent on a so-called measuring instrument such as current, voltage, charge, frequency, impedance, length, area, volume, mass, pressure, velocity, acceleration, concentration, density, and luminosity. It can be a physical quantity that can be easily measured. On the other hand, the second parameter is also referred to as a model parameter and is used in a model formula expressing the performance and characteristics of the product, and is a parameter that affects the performance and characteristics of the product due to an error during processing of the product. be able to.
 製造物は、人工的に製造される物であり、様々な分野の製品に使用される半導体デバイス、電子部品、電気部品、機械部品、光学部品などを含む。例えば、半導体デバイスには、MOSFET、バイポーラトランジスタ、ダイオード、発光デバイス、受光デバイス、メモリ、ロジックIC、アナログIC、センサなど様々なデバイスが含まれるが、本明細書では、製造物の一例としてMOSFET(以下、「デバイス」とも称する)を例に挙げて説明する。 Manufactured products are artificially manufactured products and include semiconductor devices, electronic parts, electrical parts, mechanical parts, optical parts, etc. used in products in various fields. For example, semiconductor devices include various devices such as MOSFETs, bipolar transistors, diodes, light emitting devices, light receiving devices, memories, logic ICs, analog ICs, sensors, etc., but in the present specification, MOSFETs (as an example of products) are used as examples of products. Hereinafter, a description will be given by taking (also referred to as a “device”) as an example.
 図2はMOSFETにおけるモデルパラメータの一例を示す説明図である。代表的な回路素子の一つであるMOSFETの電圧・電流特性を表すデバイスモデル(「モデル」とも称する)は、例えば、(式1)で表すことができる。 FIG. 2 is an explanatory diagram showing an example of model parameters in MOSFET. A device model (also referred to as a “model”) representing the voltage / current characteristics of a MOSFET, which is one of the typical circuit elements, can be represented by, for example, (Equation 1).
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000001
 (式1)において、太字はモデルパラメータを表し、モデルパラメータKは、電流ゲインパラメータであり、キャリア移動度と絶縁膜容量との積で表される。モデルパラメータVTHは、しきい値電圧を表す。Wはチャネル幅を表し、Lはチャネル長を表す。MOSFETの寸法はチャネル長Lとチャネル幅Wによって反映される。 In (Equation 1), bold letters represent model parameters, and model parameter K is a current gain parameter, which is represented by the product of carrier mobility and insulating film capacity. The model parameter VTH represents the threshold voltage. W represents the channel width and L represents the channel length. The dimensions of the MOSFET are reflected by the channel length L and the channel width W.
 また、線形領域と飽和領域とを緩やかに接続するため、モデルパラメータDELTAを用いて、デバイスモデルを式(2)で表すことができる。さらに、モデルパラメータCLM(チャネル長変調効果係数)、モデルパラメータMD(移動度劣化係数)、モデルパラメータMDV(移動度劣化開始電圧)を用い、チャネル長変調効果とキャリア移動度のバイアス電圧依存性を考慮すると、デバイスモデルを式(3)で表すことができる。式(3)のとおり、ドレイン電流Idは、ゲート・ソース電圧Vgsとドレイン・ソース電圧Vdsの関数として表すことができる。式(1)~(3)は、モデル情報としてのデバイスモデル(モデル)を表現するモデル式である。なお、モデル式は、一例であって、製造物の種類に応じて異なる表現となる。以下、本明細書では、便宜上、モデルパラメータK(電流ゲインパラメータ)とVTH(しきい値電圧)に着目し、他のモデルパラメータCLM、MD、MDV、DELTAについては一定値と仮定して説明する。なお、モデルパラメータは、図2に例示するパラメータに限定されるものではない。 Further, since the linear region and the saturated region are loosely connected, the device model can be expressed by the equation (2) using the model parameter DELTA. Furthermore, the model parameter CLM (channel length modulation effect coefficient), model parameter MD (mobility deterioration coefficient), and model parameter MDV (mobility deterioration start voltage) are used to determine the bias voltage dependence of the channel length modulation effect and carrier mobility. Considering this, the device model can be expressed by the equation (3). As shown in equation (3), the drain current Id can be expressed as a function of the gate-source voltage Vgs and the drain-source voltage Vds. Equations (1) to (3) are model equations expressing a device model (model) as model information. The model formula is an example and has different expressions depending on the type of the product. Hereinafter, in the present specification, for convenience, the model parameters K (current gain parameter) and VTH (threshold voltage) will be focused on, and other model parameters CLM, MD, MDV, and DELTA will be described assuming constant values. .. The model parameters are not limited to the parameters illustrated in FIG. 2.
 入力部11は、入力データを取得することができる。具体的には、入力部11は、計測値取得部としての機能を有し、特性パラメータIの計測値を取得することができる。また、入力部11は、分布情報取得部としての機能を有し、製造物に関する特性パラメータ(「第1パラメータ」とも称する)の分布を表す確率密度関数Q(I)を取得することができる。ここで、Iは特性パラメータを表す。入力部11が、特性パラメータIの計測値を取得するか、確率密度関数Q(I)を取得するかは、予め制御部10が決定しておくことができる。 The input unit 11 can acquire input data. Specifically, the input unit 11 has a function as a measured value acquisition unit, and can acquire the measured value of the characteristic parameter I. Further, the input unit 11 has a function as a distribution information acquisition unit, and can acquire a probability density function Q (I) representing the distribution of characteristic parameters (also referred to as “first parameters”) related to the product. Here, I represents a characteristic parameter. The control unit 10 can determine in advance whether the input unit 11 acquires the measured value of the characteristic parameter I or the probability density function Q (I).
 図3は特性パラメータの計測値の一例を示す模式図である。図3に示すように、半導体ウェハから、M個のデバイス(MOSFET)#1、#2、…、#Mを製造したとする。図3には、デバイス#1~#Mそれぞれのドレイン電流Idとゲート・ソース電圧Vgsとの関係を示す電圧・電流特性が図示されている。電圧・電流特性は、実際の計測値であり、ドレイン・ソース電圧Vdsが、Vds=Vds1の場合とVds=Vds2の場合について図示されている。また、デバイス#1のモデルパラメータK、VTHは、それぞれK=K1、VTH=VTH1であるとし、デバイス#2のモデルパラメータK、VTHは、それぞれK=K2、VTH=VTH2であるとし、デバイス#MのモデルパラメータK、VTHは、それぞれK=Km、VTH=VTHmであるとする。すなわち、デバイス#1~#Mの電圧・電流特性のばらつきは、各デバイスのモデルパラメータK、VTHの相違によって表現される。 FIG. 3 is a schematic diagram showing an example of measured values of characteristic parameters. As shown in FIG. 3, it is assumed that M devices (MOSFETs) # 1, # 2, ..., #M are manufactured from a semiconductor wafer. FIG. 3 shows voltage / current characteristics showing the relationship between the drain current Id and the gate / source voltage Vgs of each of the devices # 1 to #M. The voltage / current characteristics are actual measured values, and are shown for cases where the drain / source voltage Vds is Vds = Vds1 and Vds = Vds2. Further, it is assumed that the model parameters K and VTH of the device # 1 are K = K1 and VTH = VTH 1, respectively, and the model parameters K and VTH of the device # 2 are K = K2 and VTH = VTH 2, respectively. It is assumed that the model parameters K and VTH of M are K = Km and VTH = VTH, respectively. That is, the variation in the voltage / current characteristics of the devices # 1 to #M is expressed by the difference in the model parameters K and VTH of each device.
 デバイス#1の電圧・電流特性において、Vgs=Vgs1及びVds=Vds1のとき、Id=Id1となり、Vgs=Vgs2及びVds=Vds2のとき、Id=Id2となる。他のデバイス#2、…#Mも同様である。各デバイスの電圧・電流特性は、Id1=(Vgs1、Vds1)、Id2=(Vgs2、Vds2)のように関連付けることができる。1つ目のバイアス条件(Vgs1、Vds1)での電流値をId1とし、2つ目のバイアス条件(Vgs2、Vds2)での電流値をId2とすると、各デバイスのドレイン電流の分布を、電流Id1とId2で構成される電流空間(「特性パラメータ空間」とも称する)における点の分布で表すことができる。図3の例では、特性パラメータの計測値の分布(すなわち、2個のバイアス点におけるドレイン電流の分布)を図示している。入力部11は、図3に模式的に示すような、特性パラメータI(Id1、Id2)の計測値を取得することができる。なお、バイアス点の数は2個に限定されるものではなく、3個以上でもよい。 In the voltage / current characteristics of device # 1, when Vgs = Vgs1 and Vds = Vds1, Id = Id1, and when Vgs = Vgs2 and Vds = Vds2, Id = Id2. The same applies to the other devices # 2, ... # M. The voltage / current characteristics of each device can be related as Id1 = (Vgs1, Vds1), Id2 = (Vgs2, Vds2). Assuming that the current value under the first bias condition (Vgs1, Vds1) is Id1 and the current value under the second bias condition (Vgs2, Vds2) is Id2, the distribution of the drain current of each device is the current Id1. It can be represented by the distribution of points in the current space (also referred to as "characteristic parameter space") composed of and Id2. In the example of FIG. 3, the distribution of the measured values of the characteristic parameters (that is, the distribution of the drain currents at the two bias points) is illustrated. The input unit 11 can acquire the measured values of the characteristic parameters I (Id1, Id2) as schematically shown in FIG. The number of bias points is not limited to two, and may be three or more.
 分布情報生成部12は、入力部11を介して取得した計測値を用いて、特性パラメータI(Id1、Id2)の分布を表す確率密度関数Q(I)(「分布情報」とも称する)を生成することができる。 The distribution information generation unit 12 generates a probability density function Q (I) (also referred to as “distribution information”) representing the distribution of the characteristic parameters I (Id1, Id2) using the measured values acquired via the input unit 11. can do.
 図4は特性パラメータの確率密度分布の一例を示す模式図である。分布情報生成部12は、例えば、カーネル密度推定(KDE:Kernel Density Estimation)を用いて、密度分布推定を行うことができる。カーネル関数、最適なバンド幅(平滑化パラメータ)は適宜設定することができる。例えば、カーネル関数として標準的なガウス関数を用いることができる。これにより、ある母集団の標本データが与えられたとき、カーネル密度推定を用いて母集団のデータを外挿することができる。すなわち、特性パラメータの計測値がない点でも確率密度を推定できる。図4は、特性パラメータI(Id1、Id2)についての確率密度分布を示し、便宜的に模様の違いによって密度の違いを表現している。例えば、楕円状の輪郭の中央部分が最も確率密度が大きいことを表している。なお、カーネル密度推定の使用は必須ではない。例えば、特性パラメータの計測値の密度を近似する別の方法により代替できる。 FIG. 4 is a schematic diagram showing an example of the probability density distribution of characteristic parameters. The distribution information generation unit 12 can estimate the density distribution by using, for example, Kernel Density Estimation (KDE). The kernel function and the optimum bandwidth (smoothing parameter) can be set as appropriate. For example, a standard Gaussian function can be used as the kernel function. This allows you to extrapolate the population data using kernel density estimation given the sample data for a population. That is, the probability density can be estimated even when there is no measured value of the characteristic parameter. FIG. 4 shows the probability density distribution for the characteristic parameter I (Id1, Id2), and for convenience, the difference in density is expressed by the difference in pattern. For example, the central part of the elliptical contour represents the highest probability density. The use of kernel density estimation is not mandatory. For example, it can be replaced by another method of approximating the density of measured values of characteristic parameters.
 なお、入力部11が、製造物に関する特性パラメータの分布を表す確率密度関数Q(I)を取得する場合、分布情報生成部12で確率密度関数Q(I)を生成しなくてもよい。 When the input unit 11 acquires the probability density function Q (I) representing the distribution of the characteristic parameters related to the product, the distribution information generation unit 12 does not have to generate the probability density function Q (I).
 確率密度情報生成部13は、製造物に関するモデルパラメータpを用いて特性パラメータIを表現するモデル式(「デバイスモデル」、「モデル情報」ともいう)に基づいて、特性パラメータIの分布を表す確率密度関数Q(I)が再現されるように、モデルパラメータpの確率密度関数P(p)(「確率密度情報」ともいう)を生成することができる。確率密度情報生成部13の詳細は後述する。 The probability density information generation unit 13 is a probability of expressing the distribution of the characteristic parameter I based on a model formula (also referred to as “device model” or “model information”) expressing the characteristic parameter I using the model parameter p related to the product. The probability density function P (p) (also referred to as “probability density information”) of the model parameter p can be generated so that the density function Q (I) is reproduced. The details of the probability density information generation unit 13 will be described later.
 シミュレーション部15は、確率密度情報生成部13で生成したモデルパラメータpの確率密度関数P(p)及びモデル式を用いて、製造物を含むシミュレーション対象の性能を評価することができる。なお、図1では、シミュレーション部15は情報処理装置100に組み込まれた構成を示すが、シミュレーション部15を外部のシミュレータとして構成し、情報処理装置100から分離してもよい。シミュレーション部15の詳細は後述する。 The simulation unit 15 can evaluate the performance of the simulation target including the product by using the probability density function P (p) of the model parameter p generated by the probability density information generation unit 13 and the model formula. Although FIG. 1 shows a configuration in which the simulation unit 15 is incorporated in the information processing device 100, the simulation unit 15 may be configured as an external simulator and separated from the information processing device 100. The details of the simulation unit 15 will be described later.
 出力部14は、出力データとして、モデルパラメータpの確率密度関数P(p)、シミュレーション部15が行ったシミュレーション結果を出力することができる。 The output unit 14 can output the probability density function P (p) of the model parameter p and the simulation result performed by the simulation unit 15 as output data.
 記憶部16は、ハードディスク又は半導体メモリ等で構成することができ、所要の情報を記憶する。記憶部16は、入力部11を介して取得した入力データを記憶することができる。また、記憶部16は、情報処理装置100での各処理の処理結果や処理途中の情報を記憶することができる。また、記憶部16は、デバイスモデルを識別するモデルIDに対応付けて、1又は複数のモデル(モデル式)を記憶することができる。ここで、モデルIDは、同じ製造事業者の異なるデバイスモデルを識別することができ、あるいは異なる製造事業者が生成したデバイスモデルを識別することができる。 The storage unit 16 can be composed of a hard disk, a semiconductor memory, or the like, and stores required information. The storage unit 16 can store the input data acquired via the input unit 11. Further, the storage unit 16 can store the processing result of each processing in the information processing apparatus 100 and the information in the middle of processing. Further, the storage unit 16 can store one or a plurality of models (model formulas) in association with the model ID that identifies the device model. Here, the model ID can identify different device models of the same manufacturer, or can identify device models generated by different manufacturers.
 制御部10は、更新部としての機能を有し、記憶部16に記憶したモデルに対して、新しいモデルが提供され、入力部11を介して、新しいモデルを取得した場合、取得した新しいモデルで記憶部16に記憶したモデルを更新することができる。これにより、最新のモデルを用いることができる。 The control unit 10 has a function as an update unit, and when a new model is provided for the model stored in the storage unit 16 and a new model is acquired via the input unit 11, the acquired new model is used. The model stored in the storage unit 16 can be updated. This allows the latest model to be used.
 表示パネル17は、液晶パネル又は有機EL(Electro Luminescence)ディスプレイ等で構成することができる。 The display panel 17 can be composed of a liquid crystal panel, an organic EL (Electro Luminescence) display, or the like.
 操作部18は、例えば、ハードウェアキーボード、マウスなどで構成され、表示パネル17に表示されたアイコンなどの操作、文字等の入力などを行うことができる。なお、操作部18は、タッチパネルで構成してもよい。 The operation unit 18 is composed of, for example, a hardware keyboard, a mouse, or the like, and can operate an icon or the like displayed on the display panel 17 or input characters or the like. The operation unit 18 may be configured with a touch panel.
 次に、モデルパラメータの推定方法について説明する。 Next, the model parameter estimation method will be described.
 確率密度情報生成部13は、サンプリング部131、変換部132、及びモデルパラメータ生成部133を備える。サンプリング部131は、マルコフ連鎖を用いてモデルパラメータpの候補をサンプリングする。マルコフ連鎖とは、未来の状態が現在の状態だけで決定され、過去の状態と無関係であり、とりうる状態が離散的な値であるという性質を有する確率過程である。変換部132は、モデル式を用いて、サンプリングされた候補を特性パラメータIの候補に変換する。モデルパラメータ生成部133は、変換された特性パラメータIの分布を表す確率密度関数Q(I)にモンテカルロ法を適用して、サンプリングされた候補からモデルパラメータpを生成する。モンテカルロ法は、数値計算手法の一つで、乱数を用いた試行を繰り返すことにより近似解を求める手法である。以下、モデルパラメータの生成方法について説明する。 The probability density information generation unit 13 includes a sampling unit 131, a conversion unit 132, and a model parameter generation unit 133. The sampling unit 131 samples candidates for the model parameter p using a Markov chain. A Markov chain is a stochastic process in which the future state is determined only by the present state, is irrelevant to the past state, and the possible states are discrete values. The conversion unit 132 converts the sampled candidate into the candidate of the characteristic parameter I by using the model formula. The model parameter generation unit 133 applies the Monte Carlo method to the probability density function Q (I) representing the distribution of the converted characteristic parameter I, and generates the model parameter p from the sampled candidates. The Monte Carlo method is one of the numerical calculation methods, and is a method of finding an approximate solution by repeating trials using random numbers. Hereinafter, a method of generating model parameters will be described.
 図5はモデルパラメータpの確率密度関数P(p)の推定方法の概念を示す模式図である。製造物(例えば、MOSFETなど)のモデルパラメータ(例えば、電流ゲインパラメータK、しきい値電圧パラメータVTHなど)は直接の測定が困難であるため、各モデルパラメータが従う確率密度関数を知ることができない。一方、製造物の特性パラメータ(例えば、電流Iなど)は容易に測定することができ、電流分布の確率密度関数Q(I)を算出することができる。確率密度情報生成部13は、電流分布の確率密度関数Q(I)を用いてモデルパラメータの分布を再現する。すなわち、確率密度情報生成部13は、電流分布の確率密度が元となるモデルパラメータの分布を再現する処理を通じてランダムにモデルパラメータの集合を生成することにより、モデルパラメータの相関を含む分布を近似することができる。 FIG. 5 is a schematic diagram showing the concept of an estimation method of the probability density function P (p) of the model parameter p. Since it is difficult to directly measure the model parameters of the product (for example, MOSFET, etc.) (for example, current gain parameter K, threshold voltage parameter VTH, etc.), it is not possible to know the probability density function that each model parameter follows. .. On the other hand, the characteristic parameters of the product (for example, current I) can be easily measured, and the probability density function Q (I) of the current distribution can be calculated. The probability density information generation unit 13 reproduces the distribution of model parameters using the probability density function Q (I) of the current distribution. That is, the probability density information generation unit 13 approximates the distribution including the correlation of the model parameters by randomly generating a set of model parameters through a process of reproducing the distribution of the model parameters based on the probability density of the current distribution. be able to.
 図5に示すように、本実施の形態の手法は、モデルパラメータ空間と特性パラメータ空間とをモデルf(モデル式)を用いて対応付ける。本明細書では、便宜上、モデルパラメータ空間は、二つのモデルパラメータp1、p2からなり、特性パラメータ空間は、特性パラメータId1、Id2(ここでは、電流Id1、Id2)からなる、二次元空間としている。なお、モデルパラメータ空間及び特性パラメータ空間は、二次元に限定されるものではなく、より高次元であっても同様の手続きを用いることができる。また、モデルパラメータの密度分布をPとし、特性パラメータの測定結果の密度分布をQとする。 As shown in FIG. 5, in the method of this embodiment, the model parameter space and the characteristic parameter space are associated with each other by using the model f (model formula). In the present specification, for convenience, the model parameter space is composed of two model parameters p1 and p2, and the characteristic parameter space is a two-dimensional space composed of characteristic parameters Id1 and Id2 (here, currents Id1 and Id2). The model parameter space and the characteristic parameter space are not limited to two dimensions, and the same procedure can be used even if the model parameter space and the characteristic parameter space are higher in dimension. Further, let P be the density distribution of the model parameters, and let Q be the density distribution of the measurement results of the characteristic parameters.
 モデルパラメータ空間内の点p(i)は、一つのデバイスのモデルパラメータを表しており、特性パラメータ空間内の点I(i)に対応する。二つの空間の点p(i)と点I(i)とは、モデルfにより対応付けられる。モデルパラメータ空間内の点p(i)と点p′は二つの異なるモデルパラメータの組を表している。すなわち、これら2点は二つの異なるデバイスを表している。モデルパラメータの相違から、これら二つのデバイスの特性パラメータ(電流)I(i)、I′も異なっている。 The point p (i) in the model parameter space represents the model parameter of one device and corresponds to the point I (i) in the characteristic parameter space. The points p (i) and the points I (i) in the two spaces are associated with each other by the model f. The points p (i) and p'in the model parameter space represent two different sets of model parameters. That is, these two points represent two different devices. Due to the difference in model parameters, the characteristic parameters (current) I (i) and I'of these two devices are also different.
 多数のデバイスの特性パラメータ(電流)の測定により、特性パラメータ空間における特性パラメータの確率密度関数Qは既知である。本実施の形態は、モデルfの変換を介して、確率密度関数Qを再現するように、モデルパラメータ空間での確率密度関数Pを定めることができる。これは、特性パラメータのサンプルの密度関数が再現されるように、モデルf経由でモデルパラメータをサンプリングすることで達成することができる。以下、具体的な処理について説明する。 The probability density function Q of the characteristic parameter in the characteristic parameter space is known by measuring the characteristic parameter (current) of many devices. In this embodiment, the probability density function P in the model parameter space can be defined so as to reproduce the probability density function Q through the transformation of the model f. This can be achieved by sampling the model parameters via the model f so that the density function of the sample of characteristic parameters is reproduced. Hereinafter, specific processing will be described.
 図6及び図7は第1実施形態の情報処理装置100の処理手順の一例を示すフローチャートである。以下、便宜上、処理の主体を制御部10として説明する。制御部10は、入力データが、特性パラメータIの確率密度関数Q(I)又は特性パラメータIの計測値のいずれであるかを判定する(S11)。入力データが計測値である場合(S11で計測値)、制御部10は、特性パラメータIの計測値(ここでは、MOSFETのドレイン電の計測値)を取得する(S12)。 6 and 7 are flowcharts showing an example of the processing procedure of the information processing apparatus 100 of the first embodiment. Hereinafter, for convenience, the main body of the process will be described as the control unit 10. The control unit 10 determines whether the input data is the probability density function Q (I) of the characteristic parameter I or the measured value of the characteristic parameter I (S11). When the input data is a measured value (measured value in S11), the control unit 10 acquires the measured value of the characteristic parameter I (here, the measured value of the drain electricity of the MOSFET) (S12).
 制御部10は、確率密度関数Q(I)を生成し(S13)、後述のステップS15の処理を行う。ここでは、J個のバイアス点での電流によって張られるJ次元の特性パラメータ空間(「電流空間」とも称する)におけるドレイン電流の計測値にカーネル密度推定(KDE)を適用して確率密度関数Q(I)を計算することができる。なお、電流測定を行うデバイス数Nが十分多く、特性パラメータ空間のある点における確率密度が測定結果から直接求められる場合には、KDEの提供を省略してもよい。また、KDE以外の適宜の方法を用いて確率密度を求めてもよい。 The control unit 10 generates the probability density function Q (I) (S13) and performs the process of step S15 described later. Here, the probability density function Q (KDE) is applied to the measured value of the drain current in the J-dimensional characteristic parameter space (also referred to as "current space") stretched by the currents at the J bias points. I) can be calculated. If the number of devices N for current measurement is sufficiently large and the probability density at a certain point in the characteristic parameter space can be directly obtained from the measurement result, the provision of KDE may be omitted. Further, the probability density may be obtained by using an appropriate method other than KDE.
 入力データが確率密度関数である場合(S11で確率密度関数)、制御部10は、確率密度関数Q(I)を取得し(S14)、モデルパラメータセットp(1)の初期値をセットする(S15)。制御部10は、カウンタi=1とし(S16)、モデルパラメータ空間における現在のモデルパラメータセットをp(i)とする(S17)。 When the input data is a probability density function (probability density function in S11), the control unit 10 acquires the probability density function Q (I) (S14) and sets the initial value of the model parameter set p (1) (S). S15). The control unit 10 sets the counter i = 1 (S16) and sets the current model parameter set in the model parameter space to p (i) (S17).
 制御部10は、モデル式を用いて、現在のモデルパラメータセットをp(i)に対応する現在の特性パラメータ(ドレイン電流)I(i)を算出する(S18)。制御部10は、モデルパラメータ空間における現在のサンプルであるモデルパラメータセットp(i)を、提案分布q(p′/p(i))を用いて移動して、次のサンプル候補である次の候補モデルパラメータセットp′を生成する(S19)。ここで、提案分布としては、例えば、現在のサンプルp(i) を中心とし、その分散がσpであるような正規分布を用いることができ、この正規分布からのサンプルxをp(i)からp′への移動距離とする。提案分布は正規分布に限られず、詳細つり合い条件を満たす様々な分布を用いることができる。詳細つり合い条件は、遷移確率に従って状態が遷移してもそれぞれの存在確率が等しいという条件である。ここでは、モデルパラメータ空間において、マルコフ連鎖を用いてモデルパラメータの候補をサンプリングしている。制御部10は、モデル式を用いて、候補モデルパラメータセットp′に対応する特性パラメータ(ドレイン電流)I′を算出する(S20)。 The control unit 10 calculates the current characteristic parameter (drain current) I (i) corresponding to p (i) with the current model parameter set using the model formula (S18). The control unit 10 moves the model parameter set p (i) , which is the current sample in the model parameter space, using the proposed distribution q (p'/ p (i) ), and moves the next sample candidate to the next. A candidate model parameter set p'is generated (S19). Here, as the proposed distribution, for example, a normal distribution can be used in which the current sample p (i) is centered and the variance is σp, and the sample x from this normal distribution can be obtained from p (i). Let it be the distance traveled to p'. The proposed distribution is not limited to the normal distribution, and various distributions satisfying the detailed equilibrium condition can be used. The detailed equilibrium condition is a condition that the existence probabilities of each are equal even if the states transition according to the transition probabilities. Here, in the model parameter space, Markov chains are used to sample model parameter candidates. The control unit 10 calculates the characteristic parameter (drain current) I'corresponding to the candidate model parameter set p'using the model formula (S20).
 制御部10は、確率比rを、r=Q(I′)/Q(I(i))という式により算出する(S21)。ここで、Q(I′)/Q(I(i))によりP(p′)/P(p(i))を近似する。確率密度関数Qは、KDEなどによる密度推定により既知であることから即座に計算が可能であるのに対し、確率密度関数Pは、測定が困難であることによりP(p′)/P(p(i))を求めることが困難であるためである。特性パラメータ空間における確率比によりモデルパラメータ空間での確率比を近似することの妥当性は、モデル式fが単射(あるいは単写)であり(これは通常満たされる)、パラメータのばらつきが、fが線形とみなし得る程度に小さい場合には明らかである。ばらつきを考慮する範囲においてfが強い非線形性を持つ場合には、その補正を行う必要がある。制御部10は、一様乱数R(0≦R<1)を生成する(S22)。 The control unit 10 calculates the probability ratio r by the formula r = Q (I ′) / Q (I (i) ) (S21). Here, P (p') / P (p (i) ) is approximated by Q (I') / Q (I (i) ). Since the probability density function Q is known by density estimation by KDE or the like, it can be calculated immediately, whereas the probability density function P is P (p') / P (p) because it is difficult to measure. This is because it is difficult to obtain (i) ). The validity of approximating the probability ratio in the model parameter space by the probability ratio in the characteristic parameter space is that the model equation f is injective (or injective) (this is usually satisfied), and the parameter variation is f. Obviously if is small enough to be considered linear. If f has a strong non-linearity in the range where variation is taken into consideration, it is necessary to correct it. The control unit 10 generates a uniform random number R (0 ≦ R <1) (S22).
 制御部10は、r>Rであるか否かを判定する(S23)。r>Rである場合(S23でYES)、制御部10は、次の候補モデルパラメータセットp′を、次のモデルパラメータセットp(i+1)とし(p(i+1)←p′)(S24)、後述のステップS27の処理を行う。r>Rでない場合(S23でNO)、制御部10は、次の候補モデルパラメータセットp′を破棄し(S25)、現在のモデルパラメータセットp(i)を、次のモデルパラメータセットp(i+1)とする(p(i+1)←p(i))(S26)。ステップS23~S26の処理は、特性パラメータ空間において、特性パラメータの確率密度関数にモンテカルロ法を適用してモデルパラメータを生成するものである。 The control unit 10 determines whether or not r> R (S23). When r> R (YES in S23), the control unit 10 sets the next candidate model parameter set p'as the next model parameter set p (i + 1) (p (i + 1) ← p'). (S24), the process of step S27 described later is performed. If r> R (NO in S23), the control unit 10 discards the next candidate model parameter set p'(S25) and sets the current model parameter set p (i) to the next model parameter set p (i ). +1) (p (i + 1) ← p (i) ) (S26). The processing of steps S23 to S26 applies the Monte Carlo method to the probability density function of the characteristic parameter in the characteristic parameter space to generate a model parameter.
 制御部10は、カウンタi≦(N-1)であるか否かを判定し(S27)、i≦(N-1)である場合(S27でYES)、すなわち、全てのデバイスについて処理が完了していない場合、カウンタiに1を加算し(S28)、ステップS24又はS26における次のモデルパラメータセットをステップS17における現在のモデルパラメータセットとし、ステップS17以降の処理を繰り返す。繰り返しの都度、モデルパラメータセットの組が一組生成される。 The control unit 10 determines whether or not the counter is i ≦ (N-1), and if i ≦ (N-1) (YES in S27), that is, the processing is completed for all the devices. If not, 1 is added to the counter i (S28), the next model parameter set in step S24 or S26 is set as the current model parameter set in step S17, and the processing after step S17 is repeated. Each time it is repeated, a set of model parameter sets is generated.
 i≦(N-1)でない場合(S27でNO)、制御部10は、ランダムなモデルパラメータの集合{p(1)、p(2)、…、p(N)}を生成する(S29)。生成されたモデルパラメータの集合は、確率密度関数P(p)に従う。 If i≤ (N-1) (NO in S27), the control unit 10 generates a random set of model parameters {p (1) , p (2) , ..., P (N) } (S29). .. The set of generated model parameters follows the probability density function P (p).
 制御部10は、モデル式及び生成したモデルパラメータを用いて、シミュレーション対象の性能評価を実施し(S30)、確率密度関数P(p)、シミュレーション結果を出力し(S31)、処理を終了する。 The control unit 10 evaluates the performance of the simulation target using the model formula and the generated model parameters (S30), outputs the probability density function P (p), outputs the simulation result (S31), and ends the process.
 実用上は、バーンインと呼ばれる初期サンプルの棄却や、サンプルの間引きによる、アルゴリズムに起因する相関をなくす手続きを併用することがより望ましい。上述のように、本実施の形態によれば、モデルパラメータが従う分布を一切仮定することなく任意の分布のモデルパラメータセットの乱数を得ることができる。すなわち、モデルパラメータの分布を推定することができる。また、モデルパラメータ間の非線形の相関も考慮できる。上述の例では、電流分布を再現するようにモデルパラメータを生成したが、電流分布に限定されるものではなく、例えば、容量等の他の測定容易な電気的な特性のばらつきを再現するようにモデルパラメータを定めることができ、汎用性の高い手法を提供できる。 For practical purposes, it is more desirable to use a procedure called burn-in, which eliminates the correlation caused by the algorithm by rejecting the initial sample and thinning out the sample. As described above, according to the present embodiment, it is possible to obtain a random number of a model parameter set of an arbitrary distribution without assuming any distribution according to the model parameters. That is, the distribution of model parameters can be estimated. Also, non-linear correlations between model parameters can be considered. In the above example, the model parameters were generated to reproduce the current distribution, but are not limited to the current distribution, for example, to reproduce variations in other measurable electrical properties such as capacitance. Model parameters can be defined and a highly versatile method can be provided.
 これまで述べた方法に加え、より直接的に、モデルパラメータの分布に従う乱数を生成することも可能である。以下、この点について説明する。 In addition to the methods described so far, it is also possible to generate random numbers that follow the distribution of model parameters more directly. This point will be described below.
(第2実施形態)
 図8は第2実施形態の情報処理装置100の構成の一例を示す説明図である。図1に例示した第1実施形態との相違点は、入力データが特性パラメータIの計測値であり、また、特性パラメータの確率密度関数Q(I)を生成するための分布情報生成部12が不要であり、フィッティング部19を備える点である。
(Second Embodiment)
FIG. 8 is an explanatory diagram showing an example of the configuration of the information processing apparatus 100 of the second embodiment. The difference from the first embodiment illustrated in FIG. 1 is that the input data is the measured value of the characteristic parameter I, and the distribution information generation unit 12 for generating the probability density function Q (I) of the characteristic parameter It is unnecessary and includes a fitting unit 19.
 図9は計測された特性パラメータのモデル式を用いたフィッティングの一例を示す模式図である。図9に示すように、半導体ウェハから、M個のデバイス(MOSFET)#1、#2、…、#Mを製造したとする。図9には、デバイス#1~#Mそれぞれのドレイン電流Idとゲート・ソース電圧Vgsとの関係を示す電圧・電流特性が図示されている。電圧・電流特性は、実際の計測値であり、ドレイン・ソース電圧Vdsが、Vds=Vds1の場合とVds=Vds2の場合について図示されている。 FIG. 9 is a schematic diagram showing an example of fitting using a model formula of measured characteristic parameters. As shown in FIG. 9, it is assumed that M devices (MOSFETs) # 1, # 2, ..., #M are manufactured from a semiconductor wafer. FIG. 9 shows voltage / current characteristics showing the relationship between the drain current Id and the gate / source voltage Vgs of each of the devices # 1 to #M. The voltage / current characteristics are actual measured values, and are shown for cases where the drain / source voltage Vds is Vds = Vds1 and Vds = Vds2.
 フィッティング部19は、各デバイス#1~#Mについて、モデル式を用いてモデルパラメータのフィッティングを行い、モデル式を構成するモデルパラメータp1、p2(図では、モデルパラメータpk1、pk2(ここで、k=1、2、…、m))を求める。フィッティングとは、計測された特性パラメータのデータ(曲線)に最もよく当てはまるような曲線を表すモデルパラメータを求めることである。複数のデバイスに対して、モデルパラメータp1、p2のフィッティングを行うことにより、モデルパラメータ空間内におけるフィッティング済みのモデルパラメータの分布を直接得ることができる。モデルパラメータp1、p2は、例えば、電流ゲインパラメータK、しきい値電圧パラメータVTHとすることができるが、これらに限定されない。なお、フィッティングは、大多数のパラメータを抽出する手法として有効である。一部のパラメータについては、破壊的測定を用いることが許容されるならば、フィッティングを行うことなく直接測定することも可能である。ただし、この場合、デバイス自体が破壊されてしまうので、フィッティングでパラメータを推測することが望ましい。 The fitting unit 19 fits the model parameters of each device # 1 to #M using the model formula, and the model parameters p1 and p2 (in the figure, the model parameters pk1 and pk2 (here, k)) constituting the model formula. = 1, 2, ..., m)) is obtained. Fitting is to find a model parameter that represents a curve that best fits the measured characteristic parameter data (curve). By fitting the model parameters p1 and p2 to a plurality of devices, the distribution of the fitted model parameters in the model parameter space can be directly obtained. The model parameters p1 and p2 can be, for example, the current gain parameter K and the threshold voltage parameter VTH, but are not limited thereto. Fitting is effective as a method for extracting the majority of parameters. For some parameters, it is also possible to measure directly without fitting if it is permissible to use destructive measurements. However, in this case, the device itself is destroyed, so it is desirable to guess the parameters by fitting.
 図10はモデルパラメータの確率密度分布の一例を示す模式図である。確率密度情報生成部13は、例えば、カーネル密度推定(KDE:Kernel Density Estimation)を用いて、密度分布推定を行うことができる。カーネル関数、最適なバンド幅(平滑化パラメータ)は適宜設定することができる。例えば、カーネル関数として標準的なガウス関数を用いることができる。これにより、ある母集団の標本データが与えられたとき、カーネル密度推定を用いて母集団のデータを外挿することができる。すなわち、モデルパラメータの算出値がない点でも確率密度を推定できる。図10は、モデルパラメータp1、p2についての確率密度分布を示し、便宜的に模様の違いによって密度の違いを表現している。例えば、楕円状の輪郭の中央部分が最も確率密度が大きいことを表している。上述のように、複数デバイスのモデルパラメータに対してKDEを実行すれば、モデルパラメータの従う確率密度関数P(p)を得ることができる。または、あるモデルパラメータpjの周囲のパラメータ数などをもって、pjにおける確率密度関数に比例する値を得ることができるなど、確率密度関数P(p)を別の方法で近似することも可能である。 FIG. 10 is a schematic diagram showing an example of the probability density distribution of model parameters. The probability density information generation unit 13 can estimate the density distribution by using, for example, kernel density estimation (KDE: Kernel Density Estimation). The kernel function and the optimum bandwidth (smoothing parameter) can be set as appropriate. For example, a standard Gaussian function can be used as the kernel function. This allows you to extrapolate the population data using kernel density estimation given the sample data for a population. That is, the probability density can be estimated even when there is no calculated value of the model parameter. FIG. 10 shows the probability density distribution for the model parameters p1 and p2, and for convenience, the difference in density is expressed by the difference in pattern. For example, the central part of the elliptical contour represents the highest probability density. As described above, by executing KDE on the model parameters of a plurality of devices, the probability density function P (p) according to the model parameters can be obtained. Alternatively, the probability density function P (p) can be approximated by another method, for example, a value proportional to the probability density function in pj can be obtained by using the number of parameters around a certain model parameter pj.
 図11及び図12は第2実施形態の情報処理装置100の処理手順の一例を示すフローチャートである。制御部10は、特性パラメータIの計測値(ここでは、MOSFETのドレイン電流の計測値)を取得する(S41)。制御部10は、モデル式を用いてモデルパラメータのフィッティングを実施する(S42)。制御部10は、他のデバイスの有無を判定し(S43)、他のデバイスがある場合(S43でYES)、ステップS41以降の処理を繰り返す。 11 and 12 are flowcharts showing an example of the processing procedure of the information processing apparatus 100 of the second embodiment. The control unit 10 acquires the measured value of the characteristic parameter I (here, the measured value of the drain current of the MOSFET) (S41). The control unit 10 fits the model parameters using the model formula (S42). The control unit 10 determines the presence / absence of another device (S43), and if there is another device (YES in S43), repeats the processes after step S41.
 他のデバイスがない場合(S43でNO)、制御部10は、モデルパラメータの確率密度関数P(p)を生成する(S44)。制御部10は、モデルパラメータセットp(1)の初期値をセットする(S45)。制御部10は、カウンタi=1とし(S46)、モデルパラメータ空間における現在のモデルパラメータセットをp(i)とする(S47)。 When there is no other device (NO in S43), the control unit 10 generates a probability density function P (p) of the model parameter (S44). The control unit 10 sets the initial value of the model parameter set p (1) (S45). The control unit 10 sets the counter i = 1 (S46) and sets the current model parameter set in the model parameter space to p (i) (S47).
 制御部10は、モデルパラメータ空間における現在のサンプルであるモデルパラメータセットp(i)を、提案分布q(p′/p(i))を用いて移動して、次のサンプル候補である次の候補モデルパラメータセットp′を生成する(S48)。ここで、提案分布としては、例えば、現在のサンプルp(i) を中心とし、その分散がσpであるような正規分布を用いることができ、この正規分布からのサンプルxをp(i)からp′への移動距離とする。提案分布は正規分布に限られず、詳細つり合い条件を満たす様々な分布を用いることができる。ここでは、モデルパラメータ空間において、マルコフ連鎖を用いてモデルパラメータの候補をサンプリングしている。 The control unit 10 moves the model parameter set p (i) , which is the current sample in the model parameter space, using the proposed distribution q (p'/ p (i) ), and moves the next sample candidate to the next. Generate a candidate model parameter set p'(S48). Here, as the proposed distribution, for example, a normal distribution can be used in which the current sample p (i) is centered and the variance is σp, and the sample x from this normal distribution can be obtained from p (i). Let it be the distance traveled to p'. The proposed distribution is not limited to the normal distribution, and various distributions satisfying the detailed equilibrium condition can be used. Here, in the model parameter space, Markov chains are used to sample model parameter candidates.
 制御部10は、確率比rを、r=P(p′)/P(p(i))という式により算出する(S49)。制御部10は、一様乱数R(0≦R<1)を生成する(S50)。制御部10は、r>Rであるか否かを判定する(S51)。r>Rである場合(S51でYES)、制御部10は、次の候補モデルパラメータセットp′を、次のモデルパラメータセットp(i+1)とし(p(i+1)←p′)(S52)、後述のステップS55の処理を行う。 The control unit 10 calculates the probability ratio r by the formula r = P (p') / P (p (i) ) (S49). The control unit 10 generates a uniform random number R (0 ≦ R <1) (S50). The control unit 10 determines whether or not r> R (S51). When r> R (YES in S51), the control unit 10 sets the next candidate model parameter set p'as the next model parameter set p (i + 1) (p (i + 1) ← p'). (S52), the process of step S55 described later is performed.
 r>Rでない場合(S51でNO)、制御部10は、次の候補モデルパラメータセットp′を破棄し(S53)、現在のモデルパラメータセットp(i)を、次のモデルパラメータセットp(i+1)とする(p(i+1)←p(i))(S54)。制御部10は、カウンタi≦(N-1)であるか否かを判定し(S55)、i≦(N-1)である場合(S55でYES)、すなわち、全てのデバイスについて処理が完了していない場合、カウンタiに1を加算し(S56)、ステップS52又はS54における次のモデルパラメータセットをステップS47における現在のモデルパラメータセットとし、ステップS47以降の処理を繰り返す。繰り返しの都度、モデルパラメータセットの組が一組生成される。 If r> R (NO in S51), the control unit 10 discards the next candidate model parameter set p'(S53), and sets the current model parameter set p (i) to the next model parameter set p (i ). +1) (p (i + 1) ← p (i) ) (S54). The control unit 10 determines whether or not the counter is i ≦ (N-1), and if i ≦ (N-1) (YES in S55), that is, the processing is completed for all the devices. If not, 1 is added to the counter i (S56), the next model parameter set in step S52 or S54 is set as the current model parameter set in step S47, and the processing after step S47 is repeated. Each time it is repeated, a set of model parameter sets is generated.
 i≦(N-1)でない場合(S55でNO)、制御部10は、ランダムなモデルパラメータの集合{p(1)、p(2)、…、p(N)}を生成する(S57)。生成されたモデルパラメータの集合は、確率密度関数P(p)に従う。制御部10は、モデル式及び生成したモデルパラメータを用いて、シミュレーション対象の性能評価を実施し(S58)、確率密度関数P(p)、シミュレーション結果を出力し(S59)、処理を終了する。 If i≤ (N-1) (NO in S55), the control unit 10 generates a random set of model parameters {p (1) , p (2) , ..., P (N) } (S57). .. The set of generated model parameters follows the probability density function P (p). The control unit 10 performs performance evaluation of the simulation target using the model formula and the generated model parameters (S58), outputs the probability density function P (p), outputs the simulation result (S59), and ends the process.
 確率密度関数P(p)やその近似値を用いれば、モデルパラメータ空間内の任意の2点(p(i)およびp′)における確率比r=P(p′)/P(p(i))が得られるから、提案分布を用いてpを移動しながら確率比rを評価し、その採択と棄却を繰り返すマルコフ連鎖をもちいてモデルパラメータの候補をサンプリングできる。 Using the probability density function P (p) or its approximation, the probability ratio r = P (p') / P (p (i) at any two points (p (i) and p') in the model parameter space. ) Is obtained, the probability ratio r can be evaluated while moving p using the proposed distribution, and model parameter candidates can be sampled using a Markov chain that repeats its adoption and rejection.
 上述のように、第2実施形態の情報処理装置100は、製造物に関する第1パラメータ(製造物の特性を示す特性パラメータ)の計測値を取得し、取得した第1パラメータの計測値に、当該製造物に関する第2パラメータ(製造物のモデルパラメータ)を用いて第1パラメータを表現するモデル情報を当てはめて、第2パラメータの分布情報を求めることができる。情報処理装置100は、第2パラメータの分布情報に基づいて第2パラメータの確率密度情報(確率密度関数)を生成することができる。より具体的には、情報処理装置100は、マルコフ連鎖を用いて第2パラメータの候補をサンプリングし、第2パラメータの分布を表す分布情報にモンテカルロ法を適用して、サンプリングされた候補から第2パラメータを生成することができる。生成した第2パラメータは、確率密度関数に従う。 As described above, the information processing apparatus 100 of the second embodiment acquires the measured value of the first parameter (characteristic parameter indicating the characteristic of the product) relating to the product, and applies the measured value of the acquired first parameter to the measured value. The distribution information of the second parameter can be obtained by applying the model information expressing the first parameter using the second parameter (model parameter of the product) relating to the product. The information processing apparatus 100 can generate the probability density information (probability density function) of the second parameter based on the distribution information of the second parameter. More specifically, the information processing apparatus 100 samples the candidates of the second parameter using the Markov chain, applies the Monte Carlo method to the distribution information representing the distribution of the second parameter, and second from the sampled candidates. Parameters can be generated. The generated second parameter follows a probability density function.
 次に、シミュレーション部15が行う回路シミュレーションについて説明する。 Next, the circuit simulation performed by the simulation unit 15 will be described.
 図13はシミュレーション対象としてのSRAMセル回路の一例を示す回路図である。図13に示すSRAMセル回路は、前述のMOSFET及びそのデバイスモデルを用いる。SRAMセル回路において、L1、L2はpMOSFETであり、D1、D2、A1、A2はnMOSFETである。回路を製造する際に用いるpMOSFETとnMOSFETそれぞれの特性を表すデバイスモデルPMOS1、NMOS1を用意し、L1、L2にはPMOS1を、D1、D2、A1、A2にはNMOS1を用いる。ここで、SRAMセル回路の対称性から、L1、L2は同寸法で、同じMOSFETを用いて設計する。すなわち、回路図上では、モデルパラメータを含めて全く同じモデルが用いられることとなる。同様に、D1、D2は互いに寸法が同じであり、A1、A2も互いに寸法が同じである。なお、回路シミュレーションを行う際には、モデル定義文、素子定義文を含むネットリストを生成しておく。 FIG. 13 is a circuit diagram showing an example of an SRAM cell circuit as a simulation target. The SRAM cell circuit shown in FIG. 13 uses the above-mentioned MOSFET and its device model. In the SRAM cell circuit, L1 and L2 are p MOSFETs, and D1, D2, A1 and A2 are n MOSFETs. Device models polyclonal1 and EtOAc1 representing the characteristics of each of the p MOSFET and n MOSFET used in manufacturing the circuit are prepared, and photoresist1 is used for L1 and L2, and IGMP1 is used for D1, D2, A1 and A2. Here, due to the symmetry of the SRAM cell circuit, L1 and L2 are designed using the same MOSFET with the same dimensions. That is, on the circuit diagram, the exact same model including the model parameters is used. Similarly, D1 and D2 have the same dimensions as each other, and A1 and A2 also have the same dimensions as each other. When performing a circuit simulation, a netlist including a model definition statement and an element definition statement is generated.
 図14A及び図14BはSRAMセル回路のシミュレーション結果の一例を示す模式図である。図14は、SRAMセル回路の回路特性の一つである静的ノイズマージン(SNM)についての評価結果を示す。SNMは、節点N、Nの電圧が変動するときに、メモリセルが記憶している値の反転に必要な電圧変動量を表し、メモリセルに記憶した値の安定性を評価する指標の一つである。図14Aは正常動作の状態を示す。すなわち、バタフライカーブと呼ばれる2つの曲線に二つの隙間(一辺の長さがSNM、SNMである最大の正方形を内包するように示している)が存在する状態が、0/1のいずれをも安定して記憶できる正常動作状態を表している。メモリセルの設計においては、min(SNM、SNM1)の最大化が図られる。 14A and 14B are schematic views showing an example of the simulation result of the SRAM cell circuit. FIG. 14 shows the evaluation result of the static noise margin (SNM), which is one of the circuit characteristics of the SRAM cell circuit. SNM represents the amount of voltage fluctuation required for reversing the value stored in the memory cell when the voltage at the nodes N 0 and N 1 fluctuates, and is an index for evaluating the stability of the value stored in the memory cell. It is one. FIG. 14A shows a state of normal operation. That is, the state in which two gaps (indicated to include the largest square whose one side length is SNM 0 and SNM 1 ) exist in two curves called butterfly curves is 0/1. Represents a normal operating state that can be stably stored. In the memory cell design, min (SNM 0 , SNM 1 ) is maximized.
 これに対して、図14Bは動作不良の状態を示す。モデルパラメータのばらつきによって、例えば、しきい値電圧が異なれば回路特性も異なることになる。図14Bは、モデルパラメータのばらつきが大きいために回路が正常動作しない例を示している。この例では、SNM、すなわち、バタフライカーブの左側の隙間が存在しない状態となっている。このようなシミュレーションを、モデルパラメータのばらつきに対応させて多数回繰り返してシミュレーションを行うことで、静的ノイズマージンを指標とするSRAMセル回路の歩留まりを求めることができる。 On the other hand, FIG. 14B shows a state of malfunction. Due to variations in model parameters, for example, if the threshold voltage is different, the circuit characteristics will also be different. FIG. 14B shows an example in which the circuit does not operate normally due to a large variation in model parameters. In this example, SNM 0 , that is, the gap on the left side of the butterfly curve does not exist. By repeating such a simulation many times in response to variations in model parameters, it is possible to obtain the yield of the SRAM cell circuit using the static noise margin as an index.
 次に、第1実施形態の情報処理装置100によるモデルパラメータの確率密度関数の生成に対する評価結果について説明する。しきい値に代表されるモデルパラメータは直接測定することが困難な場合が多い。そこで、評価方法として2つの方法を用いる。第1の評価方法は、仮想デバイスを想定し、仮想デバイスのモデルパラメータのばらつきに対して、どの程度精度良く推定できるかを示す。第2の評価方法は、実際に測定したデバイスの電流特性パラメータに対して、どの程度精度良く推定できるかを示す。まず、第1の評価方法について説明する。 Next, the evaluation results for the generation of the probability density function of the model parameters by the information processing apparatus 100 of the first embodiment will be described. Model parameters such as thresholds are often difficult to measure directly. Therefore, two methods are used as the evaluation method. The first evaluation method assumes a virtual device and shows how accurately it can be estimated with respect to variations in model parameters of the virtual device. The second evaluation method shows how accurately the current characteristic parameters of the actually measured device can be estimated. First, the first evaluation method will be described.
 図15は500個の仮想デバイスのモデルパラメータのばらつきを示す模式図である。図において、縦軸はしきい値電圧VTHを示し、横軸は電流ゲインKを示す。仮想デバイスは、2つのロットが混在し、モデルパラメータが2つのクラスタに分かれている状況を想定した。以下では、デバイスモデルとして、前述の式(1)~(3)を用い、モデルパラメータのうち、しきい値電圧VTH、電流ゲインKの2つのモデルパラメータのばらつきを考慮する。 FIG. 15 is a schematic diagram showing variations in model parameters of 500 virtual devices. In the figure, the vertical axis indicates the threshold voltage VTH, and the horizontal axis indicates the current gain K. The virtual device assumes a situation where two lots are mixed and the model parameters are divided into two clusters. In the following, the above equations (1) to (3) are used as the device model, and the variation of the two model parameters of the threshold voltage VTH and the current gain K is considered among the model parameters.
 いずれのモデルパラメータについても、2つのピークを持つ正規分布に従うものとして、疑似乱数を用いてパラメータを生成する。これら2つ以外のモデルパラメータは、一定値に固定する。乱数を用いてモデルパラメータを生成する際、それぞれのモデルパラメータK、VTHについての平均値、分散、相関は適宜設定すればよい。500個の仮想デバイスのモデルパラメータとする2変量正規乱数を500個生成した。500個のうち、70%の350個が一方のクラスタ、30%の150個が他方のクラスタである。 For each model parameter, the parameter is generated using a pseudo-random number assuming that it follows a normal distribution with two peaks. Model parameters other than these two are fixed to constant values. When generating model parameters using random numbers, the average value, variance, and correlation for each model parameter K and VTH may be appropriately set. We generated 500 bivariate normal random numbers as model parameters for 500 virtual devices. Of the 500, 70% of 350 are one cluster and 30% of 150 are the other cluster.
 図16A及び図16Bは500個の仮想デバイスの特性パラメータである電流のばらつきの一例を示す模式図である。図15に例示した、生成した500個のモデルパラメータについて、モデル式を用いてデバイスのドレイン電流を計算する。ここでは、500個の仮想デバイスの電圧電流特性データを計算した。図16Aは、2つのバイアス条件(バイアス電圧)におけるドレイン電流の分布を示す。特性パラメータ空間を構成する電流データのバイアス条件として、ここでは(Vgs、Vds)=(10V、1V)、(18V、5V)の2点を選択している。これらはそれぞれ図16AのId1、Id2に相当する。 16A and 16B are schematic views showing an example of current variation, which is a characteristic parameter of 500 virtual devices. For the generated 500 model parameters exemplified in FIG. 15, the drain current of the device is calculated using the model formula. Here, the voltage-current characteristic data of 500 virtual devices were calculated. FIG. 16A shows the distribution of drain current under two bias conditions (bias voltage). Here, two points (Vgs, Vds) = (10V, 1V) and (18V, 5V) are selected as the bias conditions of the current data constituting the characteristic parameter space. These correspond to Id1 and Id2 in FIG. 16A, respectively.
 図16Bは、ドレイン電流のばらつきに対して、KDEを用いて推定された仮想デバイス500個のドレイン電流の確率密度分布を示す。 FIG. 16B shows the probability density distribution of the drain currents of 500 virtual devices estimated using KDE with respect to the variation of the drain currents.
 図17は比較例によるモデルパラメータの生成結果の一例を示す模式図である。図17は、比較例によって生成されたモデルパラメータK、VTHの散布図、パラメータK、VTHごとのヒストグラムを示す。図中、黒丸はモデルパラメータの真値を示し、◇は比較例によって生成されたモデルパラメータを示す。比較例によって生成されたパラメータは1つの正規分布しか仮定していないため、与えたパラメータのばらつきを再現できず誤差が大きくなっている事がわかる。この状況は、モデルパラメータごとのヒストグラムでも明らかである。実線で表されている真の分布が2つのピークを持つのに対し、比較例ではこれを正規分布で近似しているため、モデルパラメータのうち頻度の高い部分が完全に異なっている。ここで、以下に比較例による手法について説明する。 FIG. 17 is a schematic diagram showing an example of the generation result of the model parameter by the comparative example. FIG. 17 shows a scatter plot of model parameters K and VTH generated by the comparative example, and a histogram for each parameter K and VTH. In the figure, black circles indicate the true values of the model parameters, and ◇ indicates the model parameters generated by the comparative example. Since the parameters generated by the comparative example assume only one normal distribution, it can be seen that the variation of the given parameters cannot be reproduced and the error is large. This situation is also clear in the histogram for each model parameter. The true distribution represented by the solid line has two peaks, whereas in the comparative example, this is approximated by a normal distribution, so that the most frequent parts of the model parameters are completely different. Here, a method based on a comparative example will be described below.
 比較例では、概ね以下の手順でモデルパラメータを生成している。(1)製造物の性能を、加工時のパラメータ(モデルパラメータ)の関数等により表現するモデル式を生成する、(2)製造物の測定値(例えば、電流など)がモデル式にフィッティングするようにモデル式中のパラメータのばらつきを正規分布等で表す、(3)正規分布等で表されたパラメータの実現値を、その正規分布に従う乱数として生成する。さらに、製造物の性能を評価する場合には、生成した実現値とモデル式を用いてシミュレーションを行う。 In the comparative example, the model parameters are generated by the following procedure. (1) Generate a model formula that expresses the performance of the product by a function of parameters (model parameters) at the time of processing, (2) Fit the measured value (for example, current, etc.) of the product to the model formula. The variation of the parameters in the model formula is expressed by a normal distribution or the like. (3) The realization value of the parameter expressed by the normal distribution or the like is generated as a random number according to the normal distribution. Furthermore, when evaluating the performance of the product, a simulation is performed using the generated realization value and the model formula.
 比較例では、モデルパラメータを正規分布として表しているが、実際には、正規分布に従わない場合も多い。さらに、モデルパラメータのばらつきを直接測定することが困難な場合がある。測定困難で分布が不明なパラメータを、正規分布に従うと仮定するので、生成したモデルパラメータは、十分な精度が得られない場合がある。 In the comparative example, the model parameters are represented as a normal distribution, but in reality, they often do not follow the normal distribution. In addition, it can be difficult to directly measure variability in model parameters. Since it is assumed that parameters that are difficult to measure and whose distribution is unknown follow a normal distribution, the generated model parameters may not be sufficiently accurate.
 図18は第1実施形態によるモデルパラメータの生成結果の一例を示す模式図である。図18は、前述の本実施の形態によって生成されたモデルパラメータK、VTHの散布図、パラメータK、VTHごとのヒストグラムを示す。図中、黒丸はモデルパラメータの真値を示し、△は本実施の形態によって生成されたモデルパラメータを示す。なお、本実施の形態の手法において、現サンプルから次サンプル候補を作る時の乱数(提案分布)として、平均が0、分散が(0.2)の2次元正規分布から作成した乱数を用いた。すなわち、ΔK、ΔVTH ~N(0、0.22)とすると、次サンプル候補は、現サンプルから次のように生成することができる。
 (K′、VTH′)=(K(i) 、VTH(i) +(ΔK、ΔVTHT
サンプル数は55,000個である。ただし、最初の5000サンプルはバーインとして捨て、50,000個を残した。さらに、100サンプル間隔でサンプルを間引いて、最終的に500サンプルを得た。
FIG. 18 is a schematic diagram showing an example of the generation result of the model parameter according to the first embodiment. FIG. 18 shows a scatter plot of the model parameters K and VTH generated by the above-described embodiment, and a histogram for each of the parameters K and VTH. In the figure, black circles indicate the true values of the model parameters, and Δ indicates the model parameters generated by the present embodiment. In the method of this embodiment, a random number created from a two-dimensional normal distribution having an average of 0 and a variance of (0.2) 2 is used as a random number (proposal distribution) when creating a next sample candidate from the current sample. board. That is, if Δ K and Δ VTH to N (0, 0.2 2 ), the next sample candidate can be generated from the current sample as follows.
(K', VTH') T = (K (i) , VTH (i) ) T + ( ΔK, ΔVTH ) T
The number of samples is 55,000. However, the first 5000 samples were discarded as burn-ins, leaving 50,000. Further, the samples were thinned out at intervals of 100 samples to finally obtain 500 samples.
 図18に示すように、本実施の形態の手法で生成されたモデルパラメータのヒストグラムでは、モデルパラメータの分布を仮定していないが、結果として2つのピークを良く再現していることが分かる。これにより、本実施の形態は、誤差の小さいばらつき考慮シミュレーションに有用であることがわかる。 As shown in FIG. 18, in the histogram of the model parameters generated by the method of this embodiment, the distribution of the model parameters is not assumed, but as a result, it can be seen that the two peaks are well reproduced. From this, it can be seen that this embodiment is useful for a variation consideration simulation with a small error.
 次に、第2の評価方法について説明する。 Next, the second evaluation method will be described.
 図19A及び図19Bはフィッティング対象の特性パラメータである電流のばらつきの一例を示す模式図である。第1の評価方法の場合と同様、簡単のために、モデルパラメータのうち、しきい値電圧VTH、電流ゲインKの2つのモデルパラメータのばらつきを考慮する。フィッティング対象とするバイアス条件は(Vgs、Vds)=(18V、1V)、(12V、3V)の2点とし、それぞれId1、Id2とする。図19Aは、このバイアス条件の電流の実測データ(電流ばらつきデータ)を示す。図19Bは、KDEを用いて得られた電流の確率密度を示す。 19A and 19B are schematic views showing an example of variation in current, which is a characteristic parameter of the fitting target. As in the case of the first evaluation method, for the sake of simplicity, the variation of the two model parameters of the threshold voltage VTH and the current gain K is considered among the model parameters. The bias conditions to be fitted are (Vgs, Vds) = (18V, 1V) and (12V, 3V) at two points, which are Id1 and Id2, respectively. FIG. 19A shows actual measurement data (current variation data) of the current under this bias condition. FIG. 19B shows the probability density of the current obtained using KDE.
 図20は比較例による電流のばらつき結果の一例を示す模式図である。図20は、電流特性ばらつきをシミュレーションにより再現して計測値(実測データ)と比較して評価したものである。図17において説明した比較例によって生成されたモデルパラメータを用いた電流シミュレーションは、モデルパラメータが正規分布に従うという仮定がされているので、計測値に現れる電流分布の歪みを再現できていないことが分かる。 FIG. 20 is a schematic diagram showing an example of the current variation result according to the comparative example. FIG. 20 is an evaluation in which variations in current characteristics are reproduced by simulation and compared with measured values (actual measurement data). It can be seen that the current simulation using the model parameters generated by the comparative example described in FIG. 17 cannot reproduce the distortion of the current distribution appearing in the measured values because it is assumed that the model parameters follow a normal distribution. ..
 図21は第1実施形態による電流のばらつき結果の一例を示す模式図である。本実施の形態の場合、生成されたモデルパラメータによる電流特性ばらつきシミュレーションは、計測値(実測データ)を十分に再現できていることが分かる。すなわち、本実施の形態では、モデルパラメータの分布を仮定せずに、より正確な電流シミュレーションが可能であることが分かる。なお、第2実施形態においても第1実施形態と同様の評価結果を得ることができる。 FIG. 21 is a schematic diagram showing an example of the current variation result according to the first embodiment. In the case of this embodiment, it can be seen that the measured value (actual measurement data) can be sufficiently reproduced in the current characteristic variation simulation using the generated model parameters. That is, it can be seen that in the present embodiment, more accurate current simulation is possible without assuming the distribution of model parameters. In addition, the same evaluation result as in the first embodiment can be obtained also in the second embodiment.
 次に、第1実施形態及び第2実施形態の情報処理装置100を用いた場合のばらつき考慮シミュレーションのユースケースについて説明する。 Next, a use case of the variation consideration simulation when the information processing apparatus 100 of the first embodiment and the second embodiment is used will be described.
 図22は第1実施形態の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。第1実施形態の情報処理装置100を用いる場合、ばらつき考慮シミュレーションは、図22に示すように、次のように行うことができる。
(1)製造事業者(回路を製造する企業(FAB))が多数のデバイス(例えば、MOSFETなど)を製造し、製造したデバイスの特性パラメータ(例えば、電圧・電流特性)を測定する。
(2)製造事業者は、製造するデバイスのデバイスモデル(モデル式)を決定するとともに、測定した特性パラメータに対して、例えば、カーネル密度推定を適用して、特性パラメータの確率密度関数を生成する。
(3)製造事業者は、デバイスモデル及び生成した特性パラメータの確率密度関数をユーザに提供する。なお、特性パラメータの確率密度関数に代えて、特性パラメータの計測値を提供してもよい。
(4)ユーザは、製造事業者から提供されたデバイスモデル、及び特性パラメータの確率密度関数を用いて、本実施の形態の手法(すなわち、モデルパラメータ空間でのマルコフ連鎖、及び特性パラメータ空間でのモンテカルロ法に基づくモデルパラメータの生成手法)を用いて、モデルパラメータを生成する。なお、ユーザは、製造事業者から特性パラメータの計測値を取得する場合、取得した計測値を用いて特性パラメータの確率密度関数を生成することができる。
(5)ユーザは、提供されたデバイスモデル、及び生成したモデルパラメータを用いて、ばらつき考慮の回路シミュレーションを行う。
FIG. 22 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the first embodiment. When the information processing apparatus 100 of the first embodiment is used, the variation consideration simulation can be performed as follows, as shown in FIG. 22.
(1) A manufacturer (a company that manufactures circuits (FAB)) manufactures a large number of devices (for example, MOSFETs) and measures characteristic parameters (for example, voltage / current characteristics) of the manufactured devices.
(2) The manufacturer determines the device model (model formula) of the device to be manufactured, and applies, for example, kernel density estimation to the measured characteristic parameters to generate a probability density function of the characteristic parameters. ..
(3) The manufacturer provides the user with a device model and a probability density function of the generated characteristic parameter. In addition, instead of the probability density function of the characteristic parameter, the measured value of the characteristic parameter may be provided.
(4) The user uses the device model provided by the manufacturer and the probability density function of the characteristic parameter to use the method of the present embodiment (that is, the Markov chain in the model parameter space and the characteristic parameter space. Model parameters are generated using the model parameter generation method based on the Monte Carlo method). When the user acquires the measured value of the characteristic parameter from the manufacturer, the user can generate the probability density function of the characteristic parameter using the acquired measured value.
(5) The user performs a circuit simulation considering variation using the provided device model and the generated model parameters.
 ここでは、製造事業者は、以下のようなシミュレーション情報提供方法を行うことができる。すなわち、シミュレーション情報提供方法を実現するサーバは、ユーザ側の端末から、シミュレーション対象(例えば、メモリセル回路など)に関する情報を取得し、シミュレーション対象に含まれる製造物(例えば、MOSFETなど)の特性を示す特性パラメータの分布を表す分布情報(確率密度関数)をユーザ側の端末に出力し、製造物のモデルパラメータを用いて特性パラメータを表現するモデル情報(例えば、デバイスモデル、モデル式など)をユーザ側の端末に出力することができる。ここで、シミュレーション情報は、シミュレーションを実行するのに必要となる情報であり、前述の分布情報(確率密度関数)、モデル情報を含む。 Here, the manufacturer can perform the following simulation information provision method. That is, the server that realizes the simulation information providing method acquires information about the simulation target (for example, a memory cell circuit) from the terminal on the user side, and determines the characteristics of the product (for example, MOSFET) included in the simulation target. Distribution information (probability density function) representing the distribution of the characteristic parameters to be shown is output to the terminal on the user side, and model information (for example, device model, model formula, etc.) expressing the characteristic parameters using the model parameters of the product is output to the user. It can be output to the terminal on the side. Here, the simulation information is information necessary for executing the simulation, and includes the above-mentioned distribution information (probability density function) and model information.
 図23は第2実施形態の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。第2実施形態の情報処理装置100を用いる場合、ばらつき考慮シミュレーションは、図23に示すように、次のように行うことができる。
(1)製造事業者(例えば、回路を製造する企業(FAB))は、デバイス(例えば、MOSFET)を作製し、デバイスの電圧・電流特性を計測し、計測データを収集する。
(2)製造事業者は、使用するデバイスのモデル式を決定し、計測した電圧・電流特性を近似するよう(フィッティング)、モデルパラメータを抽出する。これにより、モデルパラメータの分布を得ることができる。
(3)製造事業者は、モデル式、及び当該モデル式を用いて抽出したモデルパラメータの分布をユーザに提供する。
(4)ユーザは、モデルパラメータに対して、例えば、カーネル密度推定を適用して、モデルパラメータの確率密度関数を生成する。
(5)ユーザは、モデルパラメータの確率密度関数、及びモデルパラメータ空間でのマルコフ連鎖モンテカルロ法(MCMC:Markov chain Monte Carlo)に基づくモデルパラメータの生成手法を用いて、モデルパラメータを生成する。ここで、生成されるモデルパラメータは、製造事業者から提供されたモデルパラメータの分布と同じ分布に従うが、サンプル点は異なる。
(6)ユーザは、提供されたデバイスモデル、及び生成したモデルパラメータを用いて、ばらつき考慮の回路シミュレーションを行う。
FIG. 23 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the second embodiment. When the information processing apparatus 100 of the second embodiment is used, the variation consideration simulation can be performed as follows, as shown in FIG. 23.
(1) A manufacturer (for example, a company that manufactures a circuit (FAB)) manufactures a device (for example, MOSFET), measures the voltage / current characteristics of the device, and collects measurement data.
(2) The manufacturer determines the model formula of the device to be used, and extracts the model parameters so as to approximate the measured voltage / current characteristics (fitting). This makes it possible to obtain the distribution of model parameters.
(3) The manufacturer provides the user with a model formula and a distribution of model parameters extracted using the model formula.
(4) The user applies, for example, kernel density estimation to the model parameters to generate a probability density function of the model parameters.
(5) The user generates model parameters by using the probability density function of the model parameters and the model parameter generation method based on the Markov chain Monte Carlo method (MCMC) in the model parameter space. Here, the generated model parameters follow the same distribution as the model parameter distribution provided by the manufacturer, but the sample points are different.
(6) The user performs a circuit simulation considering variation using the provided device model and the generated model parameters.
 図24は比較例の場合のばらつき考慮シミュレーションのユースケースの一例を示す模式図である。比較例は、例えば、本実施の形態の情報処理装置100を用いない場合の例を示す。比較例の場合、ばらつき考慮シミュレーションは、図24に示すように、次のように行うことができる。
(1)製造事業者(例えば、回路を製造する企業(FAB))は、デバイス(例えば、MOSFET)を作製し、デバイスの電圧・電流特性を計測し、計測データを収集する。
(2)製造事業者は、使用するデバイスのモデル式を決定し、計測した電圧・電流特性を近似するよう(フィッティング)、モデルパラメータを抽出する。製造事業者は、抽出したモデルパラメータの分布を正規分布で近似して、その平均値と標準偏差を算出する。
(3)製造事業者は、モデル式、及び当該モデル式を用いて抽出したモデルパラメータの分布(正規分布と近似したときの平均値と標準偏差)を、統計的デバイスモデルとしてユーザに提供する。
(4)ユーザは、提供されたモデル式及びモデルパラメータの分布を用いて回路シミュレーションを行う。
FIG. 24 is a schematic diagram showing an example of a use case of the variation consideration simulation in the case of the comparative example. The comparative example shows an example in the case where the information processing apparatus 100 of this embodiment is not used. In the case of the comparative example, the variation consideration simulation can be performed as follows, as shown in FIG. 24.
(1) A manufacturer (for example, a company that manufactures a circuit (FAB)) manufactures a device (for example, MOSFET), measures the voltage / current characteristics of the device, and collects measurement data.
(2) The manufacturer determines the model formula of the device to be used, and extracts the model parameters so as to approximate the measured voltage / current characteristics (fitting). The manufacturer approximates the distribution of the extracted model parameters with a normal distribution, and calculates the mean value and standard deviation.
(3) The manufacturer provides the user with a model formula and a distribution of model parameters (mean value and standard deviation when approximated to a normal distribution) extracted using the model formula as a statistical device model.
(4) The user performs a circuit simulation using the provided model formula and distribution of model parameters.
 上述の比較例の場合、ユーザが用いる回路シミュレータでは、製造事業者から提供されたモデル式を読み込む必要があり、モデル式とモデルパラメータはセットとして用いる必要がある。例えば、モデル式が使用できない場合は、モデルパラメータを提供されても、ユーザは、回路シミュレーションを実行することができない。また、特性のばらつきがモデルパラメータの分布として表現されることから、モデル式毎にモデルパラメータの分布を求める必要がある。このため、異なるモデル式を用いて回路シミュレーションを行う必要が生じた場合には、モデル式に合わせてモデルパラメータの再抽出を行う必要があり、ばらつきの表現がモデル式から独立していること、すなわち、モデル式とモデルパラメータとが分離していることが望ましい。第1実施形態によれば、図22に例示するように、特性のばらつきを、そのモデルパラメータのばらつきとして要約することができ、ばらつきがない場合と同様の手続きで回路シミュレーションを行うことを可能にしている。また、第2実施形態によれば、製造事業者から提供されるモデルパラメータに関する計測データ数よりも多くの新たなモデルパラメータを多数生成することができ、実質的にモデル式とモデルパラメータとの分離を行うことができる。 In the case of the above comparative example, in the circuit simulator used by the user, it is necessary to read the model formula provided by the manufacturer, and the model formula and the model parameter need to be used as a set. For example, if a model expression is not available, the user will not be able to perform a circuit simulation even if the model parameters are provided. Moreover, since the variation in characteristics is expressed as the distribution of model parameters, it is necessary to obtain the distribution of model parameters for each model formula. Therefore, when it becomes necessary to perform a circuit simulation using a different model formula, it is necessary to re-extract the model parameters according to the model formula, and the expression of variation is independent of the model formula. That is, it is desirable that the model formula and the model parameters are separated. According to the first embodiment, as illustrated in FIG. 22, the variation in the characteristics can be summarized as the variation in the model parameters, and the circuit simulation can be performed by the same procedure as in the case where there is no variation. ing. Further, according to the second embodiment, it is possible to generate a large number of new model parameters that are larger than the number of measurement data regarding the model parameters provided by the manufacturer, and the model formula and the model parameters are substantially separated. It can be performed.
 図25は第1実施形態及び第2実施形態の情報処理装置の構成の他の例を示す説明図である。図25において、符号200は、コンピュータである。コンピュータ200は、制御部30、入力部40、出力部50、外部I/F(インタフェース)部60、表示パネル70、及び操作部80などを備える。制御部30は、CPU31、ROM32、RAM33、I/F(インタフェース)34などを備える。 FIG. 25 is an explanatory diagram showing another example of the configuration of the information processing apparatus of the first embodiment and the second embodiment. In FIG. 25, reference numeral 200 is a computer. The computer 200 includes a control unit 30, an input unit 40, an output unit 50, an external I / F (interface) unit 60, a display panel 70, an operation unit 80, and the like. The control unit 30 includes a CPU 31, a ROM 32, a RAM 33, an I / F (interface) 34, and the like.
 入力部40は、入力データを取得する。出力部50は、出力データを出力する。I/F34は、制御部30と、入力部40、出力部50、外部I/F部60、表示パネル70、及び操作部80それぞれとの間のインタフェース機能を有する。 The input unit 40 acquires input data. The output unit 50 outputs output data. The I / F 34 has an interface function between the control unit 30, the input unit 40, the output unit 50, the external I / F unit 60, the display panel 70, and the operation unit 80, respectively.
 外部I/F部60は、コンピュータプログラムを記録した記録媒体M(例えば、DVDなどのメディア)からコンピュータプログラムを読み取ることが可能である。 The external I / F unit 60 can read the computer program from the recording medium M (for example, a medium such as a DVD) on which the computer program is recorded.
 なお、図示していないが、記録媒体Mに記録されたコンピュータプログラムは、持ち運びが自由なメディアに記録されたものに限定されるものではなく、インターネット又は他の通信回線を通じて伝送されるコンピュータプログラムも含めることができる。また、コンピュータには、複数のプロセッサを搭載した1台のコンピュータ、あるいは、通信ネットワークを介して接続された複数台のコンピュータで構成されるコンピュータシステムも含まれる。 Although not shown, the computer program recorded on the recording medium M is not limited to the one recorded on the freely portable medium, and the computer program transmitted via the Internet or another communication line is also included. Can be included. The computer also includes a single computer equipped with a plurality of processors or a computer system composed of a plurality of computers connected via a communication network.
 上述のとおり、本実施の形態によれば、比較的測定が容易な、製造物の性能などの特性パラメータの分布から、モデル(モデル式)に使用されているモデルパラメータの分布を精度よく推定することができる。モデルパラメータの分布を推定できるので、モデルパラメータ及びモデル式を用いて、より正確なばらつき考慮シミュレーションを行うことができる。また、製造物の性能分布をもとに、モデルパラメータの分布をシミュレータに陽に与えることなく、製造物の性能評価を可能とする。 As described above, according to the present embodiment, the distribution of the model parameters used in the model (model formula) is accurately estimated from the distribution of characteristic parameters such as the performance of the product, which is relatively easy to measure. be able to. Since the distribution of model parameters can be estimated, more accurate variation-considered simulation can be performed using model parameters and model formulas. Further, based on the performance distribution of the product, it is possible to evaluate the performance of the product without explicitly giving the distribution of model parameters to the simulator.
 また、本実施の形態によれば、以下のような、ばらつき考慮シミュレーションの精度劣化の原因を解決することができる。すなわち、
(1)ばらつきを表す分布の自由度:通常、乱数生成の容易さや統計的な扱いの容易さから、正規分布が仮定される。しかし、各モデルパラメータが従う分布は必ずしも明らかでない場合があり、正規分布に従わない場合もしばしばある。本実施の形態によれば、モデルパラメータの分布として任意の分布が表現できる。
(2)モデルパラメータ間の相関の扱い:モデルパラメータ間には統計的な相関が見られる場合が多いが、ピアソンの積率相関係数に代表される線形相関のみしか扱うことができない。本実施の形態によれば、モデルパラメータ間の任意の相関を表現できる。
(3)フィッティングによって抽出されたモデルパラメータの分布に歪みがある場合、あるいは電流モデルが測定された電流特性を正しく表現できない場合には、統計的デバイスモデルが(しばしば大きな)誤差を含む場合がある。本実施の形態ではフィッティングによって抽出されたモデルパラメータを使用する必要はない。
(4)ばらつきがモデルパラメータの分布として表現されているので、デバイスモデルごとにモデルパラメータの分布を求める必要がある。ユーザが使用する回路シミュレータが、製造事業者等から提供されるデバイスモデルをサポートしない場合や、異なるデバイスモデルを用いてシミュレーションを行う必要が生じた場合には、モデル式に合わせてモデルパラメータの再抽出を行って、その分布を求める必要がある。本実施の形態によれば、ばらつきの表現をデバイスモデルから独立させることができる。
Further, according to the present embodiment, it is possible to solve the following causes of deterioration in accuracy of the variation consideration simulation. That is,
(1) Degree of freedom of distribution representing variation: Normally, a normal distribution is assumed because of the ease of random number generation and statistical handling. However, the distribution that each model parameter follows may not always be clear and often does not follow a normal distribution. According to this embodiment, an arbitrary distribution can be expressed as the distribution of model parameters.
(2) Handling of correlations between model parameters: Statistical correlations are often seen between model parameters, but only linear correlations represented by Pearson's product moment correlation coefficient can be handled. According to this embodiment, any correlation between model parameters can be expressed.
(3) Statistical device models may contain (often large) errors if the distribution of model parameters extracted by fitting is distorted or if the current model cannot accurately represent the measured current characteristics. .. In this embodiment, it is not necessary to use the model parameters extracted by the fitting.
(4) Since the variation is expressed as the distribution of model parameters, it is necessary to obtain the distribution of model parameters for each device model. If the circuit simulator used by the user does not support the device model provided by the manufacturer, etc., or if it becomes necessary to perform a simulation using a different device model, the model parameters are reset according to the model formula. It is necessary to perform extraction and obtain its distribution. According to this embodiment, the expression of variation can be made independent of the device model.
 前述のとおり、本実施の形態のモデルパラメータ生成方法は、例えば、以下のとおりとすることができる。製造物のモデル式を、I=f(p)とする。ここで、pは、比較的測定が困難な第2パラメータの集合を表す。第2パラメータは、製造物のモデルパラメータとも称する。Iは、測定が容易な第1パラメータの集合を表す。第1パラメータは、製造物の特性を示す特性パラメータとも称する。ここで、2つの空間を観念することができる。1つはモデルパラメータ空間であり、モデルパラメータpが分布する空間である。他の1つは特性パラメータ空間であり、電流などの特性パラメータIが分布する空間である。本実施の形態のモデルパラメータ生成方法は、モデル式を経由して、2つの空間でマルコフ連鎖とモンテカルロ法とを交互に行う処理を繰り返すことにより、モデルパラメータ空間におけるモデルパラメータの集合を生成するものである。一般的なマルコフ連鎖モンテカルロ法(MCMC)は、確率変数xの密度分布をf(x)としたときに、マルコフ連鎖を用いて確率変数xをサンプリング(xt+1 =x+Δ)し、密度分布のf(x)の確率比、例えば、{f(xt+1 )/f(x)}に基づいて、確率変数の集合を生成するものであるが、本実施の形態とは全く相違する。 As described above, the model parameter generation method of the present embodiment can be, for example, as follows. The model formula of the product is I = f (p). Here, p represents a set of second parameters that are relatively difficult to measure. The second parameter is also referred to as a product model parameter. I represents a set of first parameters that are easy to measure. The first parameter is also referred to as a characteristic parameter indicating the characteristics of the product. Here, we can think of two spaces. One is the model parameter space, which is the space in which the model parameter p is distributed. The other one is a characteristic parameter space, which is a space in which characteristic parameters I such as electric current are distributed. The model parameter generation method of the present embodiment generates a set of model parameters in the model parameter space by repeating a process of alternately performing a Markov chain and a Monte Carlo method in two spaces via a model formula. Is. In the general Markov chain Monte Carlo method (MCMC), when the density distribution of the random variable x is f (x), the random variable x is sampled (x t + 1 = x t + Δ) using the Markov chain. A set of random variables is generated based on the probability ratio of f (x) of the density distribution, for example, {f (x t + 1 ) / f (x t )}. It's completely different.
 上述の実施の形態では、製造物として半導体デバイスの一素子であるMOSFETを例に挙げて説明したが、製造物はMOSFETに限定されるものではなく、上述のように、I=f(p)というモデル式で表現できるものであれば、様々な分野の製造物であってもよい。 In the above-described embodiment, MOSFET, which is one element of a semiconductor device, has been described as an example as a product, but the product is not limited to MOSFET, and as described above, I = f (p). As long as it can be expressed by the model formula, it may be a product in various fields.
 本実施の形態の情報処理装置は、製造物に関する第1パラメータの分布を表す分布情報を取得する分布情報取得部と、前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する確率密度情報生成部とを備える。 The information processing apparatus of the present embodiment has a distribution information acquisition unit that acquires distribution information that represents the distribution of the first parameter related to the product, and model information that expresses the first parameter using the second parameter related to the product. Based on the above, the probability density information generation unit that generates the probability density information of the second parameter is provided so that the distribution information is reproduced.
 本実施の形態のコンピュータプログラムは、コンピュータに、製造物に関する第1パラメータの分布を表す分布情報を取得し、前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する、処理を実行させる。 The computer program of the present embodiment acquires distribution information representing the distribution of the first parameter with respect to the product from the computer, and uses the second parameter with respect to the product to represent the first parameter based on the model information. , The process of generating the probability density information of the second parameter so that the distribution information is reproduced is executed.
 本実施の形態の情報処理方法は、製造物に関する第1パラメータの分布を表す分布情報を取得し、前記製造物に関する第2パラメータを用いて前記第1パラメータが表現されるモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する。 In the information processing method of the present embodiment, distribution information representing the distribution of the first parameter related to the product is acquired, and the first parameter is expressed by using the second parameter related to the product, based on the model information. Probability density information of the second parameter is generated so that the distribution information is reproduced.
 本実施の形態のシミュレーション情報提供方法は、シミュレーション対象に関する情報を取得し、前記シミュレーション対象に含まれる製造物の特性を示す特性パラメータの分布を表す分布情報を出力し、前記製造物のモデルパラメータを用いて前記特性パラメータを表現するモデル情報を出力する。 The simulation information providing method of the present embodiment acquires information about a simulation target, outputs distribution information indicating a distribution of characteristic parameters indicating the characteristics of the product included in the simulation target, and outputs model parameters of the product. The model information expressing the characteristic parameter is output by using.
 本実施の形態の情報処理装置は、前記確率密度情報生成部で生成した第2パラメータの確率密度情報及び前記モデル情報を用いて、前記製造物を含むシミュレーション対象の性能を評価する評価部を備える。 The information processing apparatus of the present embodiment includes an evaluation unit that evaluates the performance of a simulation target including the product by using the probability density information of the second parameter generated by the probability density information generation unit and the model information. ..
 本実施の形態の情報処理装置は、前記第1パラメータの計測値を取得する計測値取得部と、取得した計測値を用いて前記第1パラメータの分布を表す分布情報を生成する分布情報生成部とを備える。 The information processing apparatus of the present embodiment has a measurement value acquisition unit that acquires the measurement value of the first parameter and a distribution information generation unit that generates distribution information representing the distribution of the first parameter using the acquired measurement value. And prepare.
 本実施の形態の情報処理装置は、前記モデル情報を取得するモデル情報取得部と、取得したモデル情報を記憶する記憶部とを備える。 The information processing device of the present embodiment includes a model information acquisition unit for acquiring the model information and a storage unit for storing the acquired model information.
 本実施の形態の情報処理装置は、取得したモデル情報で前記記憶部に記憶したモデル情報を更新する更新部を備える。 The information processing apparatus of the present embodiment includes an update unit that updates the model information stored in the storage unit with the acquired model information.
 本実施の形態の情報処理装置において、前記確率密度情報生成部は、マルコフ連鎖を用いて前記第2パラメータの候補をサンプリングするサンプリング部と、前記モデル情報を用いて、サンプリングされた候補を第1パラメータの候補に変換する変換部と、変換された第1パラメータの分布を表す分布情報にモンテカルロ法を適用して、サンプリングされた候補から第2パラメータを生成するパラメータ生成部とを備える。 In the information processing apparatus of the present embodiment, the probability density information generation unit uses a Markov chain to sample the second parameter candidate, and the model information to sample the sampled candidate first. It is provided with a conversion unit for converting into parameter candidates and a parameter generation unit for generating a second parameter from sampled candidates by applying the Monte Carlo method to distribution information representing the converted distribution of the first parameter.
 本実施の形態の情報処理装置において、前記第1パラメータは、前記製造物の特性を示す特性パラメータを含み、前記第2パラメータは、前記製造物のモデルパラメータを含む。 In the information processing apparatus of the present embodiment, the first parameter includes a characteristic parameter indicating the characteristic of the product, and the second parameter includes a model parameter of the product.
 10、30 制御部
 11、40 入力部
 12 分布情報生成部
 13 確率密度情報生成部
 131 サンプリング部
 132 変換部
 133 モデルパラメータ生成部
 14、50 出力部
 15 シミュレーション部
 16 記憶部
 17、70 表示パネル
 18、80 操作部
 19 フィッティング部
 31 CPU
 32 ROM
 33 RAM
 34 I/F
 60 外部I/F部
 100 情報処理装置
 200 コンピュータ
 
10, 30 Control unit 11, 40 Input unit 12 Distribution information generation unit 13 Probability density information generation unit 131 Sampling unit 132 Conversion unit 133 Model parameter generation unit 14, 50 Output unit 15 Simulation unit 16 Storage unit 17, 70 Display panel 18, 80 Operation unit 19 Fitting unit 31 CPU
32 ROM
33 RAM
34 I / F
60 External I / F section 100 Information processing device 200 Computer

Claims (12)

  1.  製造物に関する第1パラメータの分布を表す分布情報を取得する分布情報取得部と、
     前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する確率密度情報生成部と
     を備える情報処理装置。
    A distribution information acquisition unit that acquires distribution information that represents the distribution of the first parameter related to the product,
    A probability density information generation unit that generates probability density information of the second parameter so that the distribution information can be reproduced based on model information expressing the first parameter using the second parameter related to the product. Information processing device to be equipped.
  2.  前記第1パラメータの計測値を取得する計測値取得部と、
     取得した計測値を用いて前記第1パラメータの分布を表す分布情報を生成する分布情報生成部と
     を備える請求項1に記載の情報処理装置。
    A measurement value acquisition unit that acquires the measurement value of the first parameter,
    The information processing apparatus according to claim 1, further comprising a distribution information generation unit that generates distribution information representing the distribution of the first parameter using the acquired measured values.
  3.  前記確率密度情報生成部は、
     マルコフ連鎖を用いて前記第2パラメータの候補をサンプリングするサンプリング部と、
     前記モデル情報を用いて、サンプリングされた候補を第1パラメータの候補に変換する変換部と、
     変換された第1パラメータの分布を表す分布情報にモンテカルロ法を適用して、サンプリングされた候補から第2パラメータを生成するパラメータ生成部と
     を備える、
     請求項1又は請求項2に記載の情報処理装置。
    The probability density information generation unit is
    A sampling unit that samples the candidates for the second parameter using a Markov chain,
    A conversion unit that converts sampled candidates into candidates for the first parameter using the model information, and
    It is provided with a parameter generator that applies the Monte Carlo method to the distribution information representing the converted distribution of the first parameter and generates the second parameter from the sampled candidates.
    The information processing apparatus according to claim 1 or 2.
  4.  製造物に関する第1パラメータの計測値を取得する計測値取得部と、
     取得した第1パラメータの計測値に、前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報を当てはめて前記第2パラメータの確率密度情報を生成する確率密度情報生成部と
     を備える情報処理装置。
    A measurement value acquisition unit that acquires the measurement value of the first parameter related to the product,
    A probability density information generation unit that generates probability density information of the second parameter by applying model information expressing the first parameter to the acquired measured value of the first parameter using the second parameter related to the product. Information processing device to be equipped.
  5.  前記確率密度情報生成部は、
     マルコフ連鎖を用いて前記第2パラメータの候補をサンプリングするサンプリング部と、
     前記第2パラメータの分布を表す分布情報にモンテカルロ法を適用して、サンプリングされた候補から第2パラメータを生成するパラメータ生成部と
     を備える、
     請求項4に記載の情報処理装置。
    The probability density information generation unit is
    A sampling unit that samples the candidates for the second parameter using a Markov chain,
    It is provided with a parameter generation unit that generates a second parameter from sampled candidates by applying the Monte Carlo method to the distribution information representing the distribution of the second parameter.
    The information processing apparatus according to claim 4.
  6.  前記確率密度情報生成部で生成した第2パラメータの確率密度情報及び前記モデル情報を用いて、前記製造物を含むシミュレーション対象の性能を評価する評価部を備える、
     請求項1から請求項5のいずれか一項に記載の情報処理装置。
    It is provided with an evaluation unit for evaluating the performance of a simulation target including the product by using the probability density information of the second parameter generated by the probability density information generation unit and the model information.
    The information processing apparatus according to any one of claims 1 to 5.
  7.  前記モデル情報を取得するモデル情報取得部と、
     取得したモデル情報を記憶する記憶部と
     を備える請求項1から請求項6のいずれか一項に記載の情報処理装置。
    The model information acquisition unit that acquires the model information and
    The information processing apparatus according to any one of claims 1 to 6, further comprising a storage unit for storing acquired model information.
  8.  取得したモデル情報で前記記憶部に記憶したモデル情報を更新する更新部を備える、
     請求項7に記載の情報処理装置。
    It is provided with an update unit that updates the model information stored in the storage unit with the acquired model information.
    The information processing apparatus according to claim 7.
  9.  前記第1パラメータは、前記製造物の特性を示す特性パラメータを含み、
     前記第2パラメータは、前記製造物のモデルパラメータを含む、
     請求項1から請求項8のいずれか一項に記載の情報処理装置。
    The first parameter includes a characteristic parameter indicating the characteristic of the product.
    The second parameter includes the model parameter of the product.
    The information processing apparatus according to any one of claims 1 to 8.
  10.  コンピュータに、
     製造物に関する第1パラメータの分布を表す分布情報を取得し、
     前記製造物に関する第2パラメータを用いて前記第1パラメータを表現するモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する、
     処理を実行させるコンピュータプログラム。
    On the computer
    Obtain distribution information representing the distribution of the first parameter for the product,
    Based on the model information expressing the first parameter using the second parameter related to the product, the probability density information of the second parameter is generated so that the distribution information is reproduced.
    A computer program that executes processing.
  11.  製造物に関する第1パラメータの分布を表す分布情報を取得し、
     前記製造物に関する第2パラメータを用いて前記第1パラメータが表現されるモデル情報に基づいて、前記分布情報が再現されるように前記第2パラメータの確率密度情報を生成する、
     情報処理方法。
    Obtain distribution information representing the distribution of the first parameter for the product,
    Based on the model information in which the first parameter is expressed using the second parameter related to the product, the probability density information of the second parameter is generated so that the distribution information is reproduced.
    Information processing method.
  12.  シミュレーション対象に関する情報を取得し、
     前記シミュレーション対象に含まれる製造物の特性を示す特性パラメータの分布を表す分布情報を出力し、
     前記製造物のモデルパラメータを用いて前記特性パラメータを表現するモデル情報を出力する、
     シミュレーション情報提供方法。
     
    Get information about the simulation target and
    Distribution information representing the distribution of characteristic parameters indicating the characteristics of the product included in the simulation target is output.
    Outputs model information expressing the characteristic parameter using the model parameter of the product.
    Simulation information provision method.
PCT/JP2021/042597 2020-12-10 2021-11-19 Information processing device, computer program, information processing method, and simulation information provision method WO2022124053A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020-205323 2020-12-10
JP2020205323A JP2022092491A (en) 2020-12-10 2020-12-10 Information processing device, computer program, information processing method, and simulation information providing method

Publications (1)

Publication Number Publication Date
WO2022124053A1 true WO2022124053A1 (en) 2022-06-16

Family

ID=81972846

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/042597 WO2022124053A1 (en) 2020-12-10 2021-11-19 Information processing device, computer program, information processing method, and simulation information provision method

Country Status (2)

Country Link
JP (1) JP2022092491A (en)
WO (1) WO2022124053A1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007108843A (en) * 2005-10-11 2007-04-26 Toshiba Corp Semiconductor device design support method, semiconductor device design support system and semiconductor device design support program
WO2007049555A1 (en) * 2005-10-24 2007-05-03 Kyoto University Cmos model making apparatus, its method, program of that method and recording medium
JP2011070430A (en) * 2009-09-25 2011-04-07 Toshiba Corp Variation distribution simulation device and method
JP2012074574A (en) * 2010-09-29 2012-04-12 Hitachi Ltd Control system for processing apparatus and method for controlling processing apparatus

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007108843A (en) * 2005-10-11 2007-04-26 Toshiba Corp Semiconductor device design support method, semiconductor device design support system and semiconductor device design support program
WO2007049555A1 (en) * 2005-10-24 2007-05-03 Kyoto University Cmos model making apparatus, its method, program of that method and recording medium
JP2011070430A (en) * 2009-09-25 2011-04-07 Toshiba Corp Variation distribution simulation device and method
JP2012074574A (en) * 2010-09-29 2012-04-12 Hitachi Ltd Control system for processing apparatus and method for controlling processing apparatus

Also Published As

Publication number Publication date
JP2022092491A (en) 2022-06-22

Similar Documents

Publication Publication Date Title
US11551153B2 (en) Localized learning from a global model
Blume-Kohout et al. Robust, self-consistent, closed-form tomography of quantum logic gates on a trapped ion qubit
Chorin et al. Discrete approach to stochastic parametrization and dimension reduction in nonlinear dynamics
Póczos et al. Distribution-free distribution regression
JP6574527B2 (en) Time-series data feature extraction device, time-series data feature extraction method, and time-series data feature extraction program
Tal Calibration: Modelling the measurement process
Hacine-Gharbi et al. Low bias histogram-based estimation of mutual information for feature selection
JP2012058972A (en) Evaluation prediction device, evaluation prediction method, and program
US20180285729A1 (en) Reservoir computing system
US11216534B2 (en) Apparatus, system, and method of covariance estimation based on data missing rate for information processing
Harrell et al. Resampling, validating, describing, and simplifying the model
CN114462613B (en) Quantum computer performance characterization method and device, electronic equipment and medium
Liu et al. Misspecification analysis of two‐phase gamma‐Wiener degradation models
Demos et al. Comparing nested data sets and objectively determining financial bubbles’ inceptions
WO2022124053A1 (en) Information processing device, computer program, information processing method, and simulation information provision method
Bergot et al. A study on the optimization of the deployment of targeted observations using adjoint‐based methods
JP7231829B2 (en) Machine learning program, machine learning method and machine learning apparatus
Ghosh et al. On generalized conditional cumulative past inaccuracy measure
Liu et al. Field reliability prediction based on degradation data and environmental data
KR102311412B1 (en) Parameter estimation method based on error term independence correction, and its program implementation and performance system
CN111027448B (en) Video behavior category identification method based on time domain inference graph
Vanck et al. Using hyperbolic cross approximation to measure and compensate covariate shift
JP6985687B2 (en) PU classification device, PU classification method, and PU classification program
Legler et al. Machine learning for estimating parameters of a convective-scale model: A comparison of neural networks and random forests
Sharma et al. Estimating Hamiltonian fluctuations from quantum time averages

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21903157

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21903157

Country of ref document: EP

Kind code of ref document: A1