For manufacturing an integrated circuit, the production of a design for the circuit that comprises a plurality of mos transistors is controlled by employment of a circuit simulator. ##EQU1## are calculated in the circuit simulator for the terminal nodes of the mos transistors upon prescription of the voltages between gate and source Vgs, between drain and source Vds, and between the substrate and source Vbs in a consistent transistor model wherein drift, diffusion and short-channel effects are taken into consideration.

Patent
   5761082
Priority
Oct 01 1993
Filed
Feb 12 1997
Issued
Jun 02 1998
Expiry
Aug 30 2014
Assg.orig
Entity
Large
22
5
all paid
3. A method for manufacturing an integrated circuit having a plurality of mos transistors, comprising the steps of:
employing a circuit simulator to produce a design for the circuit to be integrated containing said plurality of mos transistors;
calculating with the circuit simulator parameters for the mos transistors including drain current, charges in the transistor, and voltages across terminals of the transistors given prescription of voltages between gate and source, between drain and source, and between substrate and source, and taking into consideration motion of charge carriers due to drift and due to diffusion as well as due to short-channel effects when making said calculations so that a transistor model is provided which is consistent in describing physical conditions in the transistor without requiring fit parameters, except for mobility of the charge carriers to designate different operating regions of each transistor;
using said design and said calculated parameters for each of said plurality of mos transistors, manufacturing said integrated circuit with said plurality of mos transistors; and
for determining the mobility μ, a current-voltage characteristic Ids (Vgs, Vds, Vbs) of a mos transistor is measured and assuming that μ in the channel of the mos transistor is location independent, calculating μ from ##EQU18## where: φSO is a surface potential at the source side,
φSL is a surface potential at the drain side,
Cox is oxide capacitance,
W is channel width,
L is channel length,
β is thermic voltage,
q is electronic charge,
LD is Debye length,
φf (y) is quasi-fermi potential,
Ids is drain-source current,
Vgs is gate-source voltage,
Vds is drain-source voltage,
Vbs is substrate-source voltage,
Nsub is dopant concentration in the substrate,
V'G =Vgs -Vfb -ΔV'G ##EQU19## Vfb is flat-band voltage, εsi is dielectric constant of silicon,
Eyy is a gradient of a lateral electrical field in the channel, and
φf (y) is quasi-fermi potential.
4. A method for manufacturing an integrated circuit having a plurality of mos transistors, comprising the steps of:
employing a circuit simulator to produce a design for the circuit to be integrated containing said plurality of mos transistors by employing a transistor model;
calculating with the circuit simulator parameters of the mos transistors based on said transistor model with said circuit simulator taking into consideration motion of charge carriers due to drift and due to diffusion as well as due to short-channel effects when making said calculations so that a transistor model is provided which is the same in describing physical conditions in the transistor for substantially all operating regions of the transistor without requiring fit parameters, except for mobility of the charge carrier, for different operating regions, or different models for different operating regions, of each transistor;
for determining the mobility μ, a current-voltage characteristic Ids (Vgs, Vds, Vbs) of a mos transistor is measured and assuming that μ in the channel of the mos transistor is location independent, calculating μ from ##EQU20## where: φSO is a surface potential at the source side,
φSL is a surface potential at the drain side,
Cox is oxide capacitance,
W is channel width,
L is channel length,
β is thermic voltage,
q is electronic charge,
LD is Debye length,
φf (y) is quasi-fermi potential,
Ids is drain-source current,
Vgs is gate-source voltage,
Vds is drain-source voltage,
Vbs is substrate-source voltage,
Nsub is dopant concentration in the substrate,
V'G =Vgs -Vfb -ΔV'G ##EQU21## Vfb is flat-band voltage, εsi is dielectric constant of silicon,
Eyy is a gradient of a lateral electrical field in the channel, and
φf (y) is quasi-fermi potential; and
using said design and said calculated parameters for each of said plurality of mos transistors, manufacturing said integrated circuit with said plurality of mos transistors.
2. A method for manufacturing an integrated circuit having a plurality of mos transistors, comprising the steps of:
employing a circuit simulator to produce a design for the circuit to be integrated containing said plurality of mos transistors;
calculating with the circuit simulator parameters for said plurality of mos transistors, said parameters taking into account motion of charge carriers due to drift and due to diffusion as well as due to short-channel effects so that a same transistor model is provided which describes physical conditions in the transistor for substantially on entire operating region of the transistors without requiring fit parameters, except for mobility, for different operating regions, or different models for different operating regions, of each transistor, said parameters including drain current Ids, overall charges Q for terminal nodes source, drain, gate, and substrate as well as for an inversion layer, and deviations ##EQU15## for the terminal nodes of said plurality of mos transistors upon prescription of voltages between gate and source Vgs, between drain and source Vds and between the substrate and source Vbs, said parameters being calculated by solving the following equations: ##EQU16## where: Cox is oxide capacitance,
φs (y) is surface potential,
q is electronic charge,
Nsub is dopant concentration in the substrate,
LD is Debye length,
Vfb is flat-band voltage,
Qi (y), Qb (y) are location-dependent charges in the inversion layer or the substrate per unit area,
Vgs is the voltage between gate and source,
εsi is dielectric constant of silicon,
Eyy is a gradient of a lateral electrical field in the channel,
β is thermic voltage,
QI, QB, QG, QS, QD are respectively overall charge in the inversion layer , in the substrate , in the gate, in the source, and in the drain,
W is channel width,
L is channel length,
Ids is drain current,
φf (y) is quasi-fermi potential,
N(y) is a location-dependent dopant concentration in the channel, and
μ is mobility;
using said design and said calculated parameters for each of said plurality of mos transistors, manufacturing said integrated circuit with said plurality of mos transistors; and
for determining the mobility μ, a current-voltage characteristic Ids (Vgs, Vds, Vbs) of a mos transistor is measured and assuming that μ in the channel of the mos transistor is location independent, calculating μ from ##EQU17## whereby φSO is a surface potential at the source side and φSL is a surface potential at the drain side.
1. A method for manufacturing an integrated circuit having a plurality of mos transistors, comprising the steps of:
employing a circuit simulator to produce a design for the circuit to be integrated containing said plurality of mos transistors;
calculating with the circuit simulator parameters for said plurality of mos transistors, said parameters taking into account motion of charge carriers due to drift and due to diffusion as well as due to short-channel effects so that a same transistor model is provided which describes physical conditions in the transistor for substantially an entire operating region of the transistors without requiring fit parameters, except for mobility, for different operating regions, or different models for different operating regions, of each transistor, said parameters including drain current Ids, overall charges Q for terminal nodes source, drain, gate, and substrate as well as for an inversion layer, and deviations ##EQU11## for the terminal nodes of said plurality of mos transistors upon prescription of voltages between gate and source Vgs, between drain and source Vds and between the substrate and source Vbs, said parameters being calculated by solving the following equations: ##EQU12## where: Cox is oxide capacitance,
φS (y) is surface potential,
q is electronic charge,
Nsub is dopant concentration in the substrate,
LD is Debye length,
Vfb is flat-band voltage,
Qi (y), Qb (y) are location-dependent charges in the inversion layer or the substrate per unit area,
Vgs is the voltage between gate and source,
εsi is dielectric constant of silicon,
Eyy is a gradient of a lateral electrical field in the channel,
β is thermic voltage,
QI, QB, QG, QS, QD are respectively overall charge in the inversion layer, in the substrate, in the gate, in the source, and in the drain,
W is channel width,
L is channel length,
Ids is drain current,
φf (y) is quasi-fermi potential,
N(y) is a location-dependent dopant concentration in the channel, and
μ is mobility;
for determining the mobility μ, a current-voltage characteristic Ids (Vgs, Vds, Vbs) of a mos transistor is measured and assuming that μ in the channel of the mos transistor is location independent, calculating μ from ##EQU13## where: φSO is a surface potential at the source side,
φSL is a surface potential at the drain side,
Cox is oxide capacitance,
W is channel width,
L is channel length,
β is thermic voltage,
q is electronic charge,
LD is Debye length,
φf (y) is quasi-fermi potential,
Ids is drain-source current,
Vgs is gate-source voltage,
Vds is drain-source voltage,
Vbs is substrate-source voltage,
Nsub is dopant concentration in the substrate,
V'G =Vgs -Vfb -ΔV'G ##EQU14## Vfb is flat-band voltage, εsi is dielectric constant of silicon,
Eyy is a gradient of a lateral electrical field in the channel, and
φf (y) is quasi-fermi potential; and
using said design and said calculated parameters for each of said plurality of mos transistors, manufacturing said integrated circuit with said plurality of mos transistors.

This is a continuation of application Ser. No. 08/298,169, filed Aug. 30, 1994, abandoned now.

Modern integrated circuits frequently comprise 106 or more transistors.

A circuit diagram is first produced on the basis of a functional description of the circuit for manufacturing such complex integrated circuits. The properties of the circuit are calculated with the assistance of a circuit simulator. Given deviations of the switching behavior from the required functional description, modifications of the circuit diagram are undertaken and the modified circuit diagram is again calculated in the circuit simulator. This cycle is run through until the properties of the circuit diagram agree with the required specifications. The production of a layout then occurs. The masks for the process sequence, with which the integrated circuit is ultimately realized, are produced on the basis of the layout.

The relevant electrical circuit quantities for MOS transistors are calculated in the circuit simulator with the assistance of a transistor model for the terminal nodes of the MOS transistor, i.e. source, drain, gate, and substrate. How well the switching behavior of the integrated circuit realized at the end coincides with the production with the circuit simulator is critically dependent on the quality of the transistor model.

For checking transistor models, simulation calculations of ring oscillators are compared to measurements of the ring oscillators. It turns out that the ring oscillator frequencies identified in the simulation calculations do not agree with the measured frequencies. The discrepancy becomes greater with increasing reduction of the structural sizes of the MOS transistors.

It has been proposed to modify the transistor models by introducing fit parameters which realistically reflect the physical conditions below strong inversion region, even given reduced structural size of the transistors; see, for example, H. J. Park, "Charge Sheet and Non-Quasistatic MOSFET Models for SPICE", Ph.D. dissertation Memorandum No. UCB/ERL M89/20, Electronics Research Lab., Dep. EECS, University of California, Berkeley, Feb. 24, 1989, incorporated herein.

The fit parameters are identified by comparing simulation calculations to measurements. Transient measurements can be required for this purpose.

Since the required fit parameters are non-physical, no reliable predictions for the fit parameters can be made for a modified transistor technology or reduced structural size. Renewed measurements are required for every transistor design.

It is an object of the invention to specify a method for manufacturing an integrated circuit wherein the production of a design for the circuit is controlled upon employment of a circuit simulator for which a matching of fit parameters to measurements is not required.

According to the method of the invention for manufacturing an integrated circuit having a plurality of MOS transistors, a circuit simulator is employed to produce a design for the circuit to be integrated. The circuit simulator is utilized to calculate drain current Ids and overall charges Q for terminal nodes source, drain, gate, and substrate as well as for an inversion layer and derivations ##EQU2## for the terminal nodes of the MOS transistors for upon prescribed voltages between gate and source Vgs, between drain and source Vds and between the substrate and source Vbs. These parameters are calculated by solving the following equations: ##EQU3## where: Cox is oxide capacitance,

φs (y) is surface potential,

q is electronic charge,

Nsub is dopant concentration in the substrate,

LD is Debye length,

Vfb is flat-band voltage,

Qi (y), Qb (y) are location-dependent charges in the inversion layer or the substrate per unit area,

Vgs is the voltage between gate and source,

εsi is the dielectric constant of silicon,

Eyy is a gradient of a lateral electrical field in the channel,

β is thermic voltage,

QI, QB, QG, QS, QD are respectively overall charge in the inversion layer, in the substrate, in the gate, in the source, and in the drain,

W is channel width,

L is channel length,

Ids is drain current,

φf is quasi-fermi potential,

N(y) is a location-dependent dopant concentration in the channel, and

μ is mobility.

In the method of the invention, the electrical quantities of the transistors are calculated according to a transistor model that takes the motion of the charge carriers due to drift and due to diffusion as well as due to the short-channel effect into consideration. This transistor model is consistent, i.e. it describes the physical conditions in the transistor without requiring fit parameters except for the mobility. The mobility μ of the charge carriers can be determined from a standard mobility mode, or by fitting to the measured current Ids between source and drain.

In transistor models previously employed in circuit simulators, the diffusion amount in the charge carrier transport is neglected. This, however, dominates the sub-threshold region. The dependency of the charges on the applied voltage is overestimated when the diffusion amount is left out of consideration. The physical conditions are realistically reproduced, since the method of the invention employs a transistor model that takes both the drift as well as the diffusion contribution to the charge carrier transport into consideration and which, beyond, takes short-channel effects into consideration.

FIG. 1 is a graph illustrating calculated capacitances dependent on applied voltage between source and drain;

FIG. 2 is a graph illustrating a response function Vout of a ring oscillator dependent on time t; and

FIG. 3 is a block diagram of the method of the invention for manufacturing an integrated circuit with a plurality of MOS transistors.

FIG. 1 shows calculated intrinsic capacitances for a 0.65 μm channel length dependent on the applied voltage between source and drain Vds. The solid-line curves 11, 12, 13, and 14 show the result of the calculations of the capacitance Cgs between gate and source according to the transistor model that is employed in the method of the invention, whereby the voltage Vgs between gate and source amounts to 1, 2, 3 or, respectively, 4 volts. For comparison, the broken-line curves 21, 22, 23, and 24 show the corresponding capacitances Cgs according to a transistor model as employed in known methods and which leaves the diffusion amount and the short-channel effect out of consideration. The solid-line curves 31, 32, 33, and 34 show the result of the calculation of capacitances Cgd between gate and drain according to the transistor model that is employed in the method of the invention. The voltage Vds between gate and source again amounts to 1, 2, 3, or 4 volts. For comparison, the broken-line curves 41 are shown, these again reproducing the calculated capacitance Cgd between gate and drain according to a conventional transistor model.

FIG. 2 shows the response function Vout of a ring oscillator, dependent on the time t. The curve referenced 1 is thus the result of a calculation upon employment of the consistent transistor model that is employed in the method of the invention. The curve referenced 2 shows the result of the calculation according to a conventional transistor model.

The consistent transistor model employed in the method of the invention is based on charges that are stored in the gate, in the inversion layer, in the drain, in the source, and in the substrate. Qg, Qi, Qd, Qs, and Qb reference these charges per unit area. Valid for Qi and Qb as a function of the location y along the channel are the following equations: ##EQU4##

Cox is the oxide capacitance, q is the electronic charge, Nsub is the dopant concentration in the substrate, Ld is the Debye length, and Vfb is the flat-band voltage. φs (y) is the location-dependent surface potential that is identified by solving the Poisson equation and the drift-diffusion equation. β is the thermal voltage. εsi is the dielectric constant of silicon. Eyy is the gradient of the lateral electrical field. Vgs is the voltage between gate and source.

The term ΔV'G takes the short-channel effect into consideration. Given reduction of the channel length, the lateral electrical field increases, so that the lateral electrical field can even dominate the vertical electrical field for short-channel transistors. Eyy can no longer be left out of consideration in this case. The value for Eyy is calculated from the measured threshold voltage of the MOS transistor.

The overall charge Qi in the inversion layer and QB in the substrate are calculated by integration over the channel direction y from source to drain: ##EQU5##

W is the channel width and L is the channel length. The following is valid for the overall charge on the gate QG, source QS, and drain QD :

QG =-(QI +QB) (8)

QI =QS +QD (9)

The partitioning of QI, into QS and QD is undertaken by the following equations in accordance with S. Y. Oh et al, IEEE J. Solid-State Circuits, SC 15, 1980, pp. 636-643, incorporated herein: ##EQU6##

Given the assumption that the current in the channel flows along the surface (charge-sheet approximation), the following is valid for the current Ids between drain and source: ##EQU7##

The mobility of the charge carriers is referenced μ. φf (y) is the quasi-Fermi potential that is described taking both the diffusion amount as well as the drift amount into consideration. The term ##EQU8## thus describes the amount due to drift, and the term ##EQU9## describes the diffusion amount.

N(y) is the location-dependent dopant concentration in the inversion layer.

Given the assumption that the mobility μ is independent of the position in the channel, the following is valid for the drain current Ids in the charge-sheet model: ##EQU10##

The value of the surface potential at the source side φSO and of the surface potential at the drain side φSL are calculated by solving the Poisson equation and equation 13.

By algebraic transformations, the equations for the overall charges Q can be described as functions of φSO and φSL.

The mobility μ can be calculated from a standard mobility model or from measured characteristics of the MOS transistors.

The quantities Ids required for the circuit simulation as well as the overall charges Q and the intrinsic capacitances, which are the derivations of the charge at the terminal node according to the applied voltages Vgs between gate and source, Vds between drain and source, and Vbs between substrate and source, are described by algebraic transformation as analytic expressions that are functions of φSO and φSL. φSO and φSL are directly dependent on technological parameters such as the oxide thickness TXO and the dopant concentration in the substrate Nsub.

As shown in FIG. 3, according to the method of the invention, an integrated circuit is manufactured having a plurality of MOS transistors by employing a circuit simulator 100 to produce a design for the integrated circuit. Into circuit simulator 100 circuit parameters are input, including the input parameters for the MOS transistor as previously described. The circuit simulator 100 outputs the circuit design parameters including the MOS transistor parameters, which are then used in the manufacture of the integrated circuit having the plurality of MOS transistors as shown at 200 in FIG. 3.

Although various minor changes and modifications might be suggested by those skilled in the art, it will be apparent that I wish to include within the scope of the patent warranted hereon all such changes and modifications as reasonably come within my contribution to the art.

Miura-Mattausch, Mitiko

Patent Priority Assignee Title
6426644, Mar 22 1999 SEMICONDUCTOR PHYSICS LABORATORY, INC Apparatus and method for determining the active dopant profile in a semiconductor wafer
6714027, Jun 15 1999 Seiko Epson Corporation Method and apparatus for calculating the electrical characteristics of materials of thin film transistors
6812047, Mar 08 2000 Applied Materials, Inc Evaluating a geometric or material property of a multilayered structure
6812717, Mar 05 2001 SEMICONDUCTOR PHYSICS LABORATORY, INC Use of a coefficient of a power curve to evaluate a semiconductor wafer
6845347, Jul 21 2000 Winbond Electronics Corporation Method for modeling an integrated circuit including a DRAM cell
6851097, Sep 14 2000 Cadence Design Systems, INC MOSFET modeling for IC design accurate for high frequencies
6885458, Mar 22 1999 SEMICONDUCTOR PHYSICS LABORATORY, INC Apparatus and method for determining the active dopant profile in a semiconductor wafer
6906801, Jun 10 1998 Applied Materials, Inc Measuring a property of a layer in multilayered structure
6911349, Feb 16 2001 Applied Materials, Inc Evaluating sidewall coverage in a semiconductor wafer
6940592, Oct 09 2001 Applied Materials, Inc Calibration as well as measurement on the same workpiece during fabrication
6958814, Mar 01 2002 Applied Materials, Inc Apparatus and method for measuring a property of a layer in a multilayered structure
6963393, Sep 23 2002 SEMICONDUCTOR PHYSICS LABORATORY, INC Measurement of lateral diffusion of diffused layers
6971791, Mar 01 2002 Applied Materials, Inc Identifying defects in a conductive structure of a wafer, based on heat transfer therethrough
7130055, Mar 05 2001 SEMICONDUCTOR PHYSICS LABORATORY, INC Use of coefficient of a power curve to evaluate a semiconductor wafer
7141440, Mar 01 2002 Applied Materials, Inc. Apparatus and method for measuring a property of a layer in a multilayered structure
7213218, Dec 18 2003 SEMICONDUCTOR TECHNOLOGY Simulation model for design of semiconductor device, thermal drain noise analysis method, simulation method, and simulation apparatus
7313770, Sep 14 2000 Cadence Design Systems, Inc. MOSFET modeling for IC design accurate for high frequencies
7465591, Mar 08 2000 Applied Materials, Inc Evaluating a geometric or material property of a multilayered structure
7675372, Aug 09 2006 Qualcomm Incorporated Circuit simulator parameter extraction using a configurable ring oscillator
7735033, Sep 14 2000 Cadence Design Systems, Inc. MOSFET modeling for IC design accurate for high frequencies
8250508, Jun 18 2008 Renesas Electronics Corporation Method and apparatus for analysis and design of a semiconductor device using impurity concentration distribution
9032348, Dec 07 2012 Arizona Board of Regents on behalf of Arizona State University Physics-based reliability model for large-scale CMOS circuit design
Patent Priority Assignee Title
4907180, May 04 1987 Hewlett-Packard Company Hardware switch level simulator for MOS circuits
5051911, Nov 25 1988 Fujitsu Limited Apparatus for effecting simulation of a logic circuit and method for producing a semiconductor device using the simulation approach
5245543, Dec 21 1990 Texas Instruments Incorporated Method and apparatus for integrated circuit design
5345401, Apr 27 1992 Renesas Electronics Corporation Logic circuit simulator for verifying circuit operations having MOS transistors
5440160, Jan 28 1992 THUNDERBIRD TECHNOLOGIES, INC High saturation current, low leakage current fermi threshold field effect transistor
//////////
Executed onAssignorAssigneeConveyanceFrameReelDoc
Aug 22 1994MIURA-MATTAUSCH, MITIKOSiemens AktiengesellschaftASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0275750988 pdf
Feb 12 1997Siemens Aktiengesellschaft(assignment on the face of the patent)
May 22 2006Siemens AktiengesellschaftInfineon Technologies AGASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0177180394 pdf
Sep 05 2006Infineon Technologies AGMosaid Technologies IncorporatedASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS 0273980190 pdf
Dec 23 2011658276 N B LTD ROYAL BANK OF CANADAU S INTELLECTUAL PROPERTY SECURITY AGREEMENT FOR NON-U S GRANTORS - SHORT FORM0275120196 pdf
Dec 23 2011Mosaid Technologies IncorporatedROYAL BANK OF CANADAU S INTELLECTUAL PROPERTY SECURITY AGREEMENT FOR NON-U S GRANTORS - SHORT FORM0275120196 pdf
Dec 23 2011658868 N B INC ROYAL BANK OF CANADAU S INTELLECTUAL PROPERTY SECURITY AGREEMENT FOR NON-U S GRANTORS - SHORT FORM0275120196 pdf
Jun 11 2014ROYAL BANK OF CANADACONVERSANT IP N B 868 INC RELEASE OF SECURITY INTEREST0334840344 pdf
Jun 11 2014ROYAL BANK OF CANADACONVERSANT IP N B 276 INC RELEASE OF SECURITY INTEREST0334840344 pdf
Jun 11 2014ROYAL BANK OF CANADACONVERSANT INTELLECTUAL PROPERTY MANAGEMENT INC RELEASE OF SECURITY INTEREST0334840344 pdf
Date Maintenance Fee Events
Jan 31 2001ASPN: Payor Number Assigned.
Nov 14 2001M183: Payment of Maintenance Fee, 4th Year, Large Entity.
Dec 26 2001REM: Maintenance Fee Reminder Mailed.
Dec 21 2005REM: Maintenance Fee Reminder Mailed.
May 08 2006M1552: Payment of Maintenance Fee, 8th Year, Large Entity.
May 08 2006M1555: 7.5 yr surcharge - late pmt w/in 6 mo, Large Entity.
May 23 2006ASPN: Payor Number Assigned.
May 23 2006RMPN: Payer Number De-assigned.
Nov 04 2009M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
Jun 02 20014 years fee payment window open
Dec 02 20016 months grace period start (w surcharge)
Jun 02 2002patent expiry (for year 4)
Jun 02 20042 years to revive unintentionally abandoned end. (for year 4)
Jun 02 20058 years fee payment window open
Dec 02 20056 months grace period start (w surcharge)
Jun 02 2006patent expiry (for year 8)
Jun 02 20082 years to revive unintentionally abandoned end. (for year 8)
Jun 02 200912 years fee payment window open
Dec 02 20096 months grace period start (w surcharge)
Jun 02 2010patent expiry (for year 12)
Jun 02 20122 years to revive unintentionally abandoned end. (for year 12)