Method for characterizing the upset response of CMOS circuits using alpha-particle sensitive test circuits

Abstract

A method for predicting the SEU susceptibility of a standard-cell D-latch using an alpha-particle sensitive SRAM, SPICE critical charge simulation results, and alpha-particle interaction physics. A technique utilizing test structures to quickly and inexpensively characterize the SEU sensitivity of standard cell latches intended for use in a space environment. This bench-level approach utilizes alpha particles to induce upsets in a low LET sensitive 4-k bit test SRAM. This SRAM consists of cells that employ an offset voltage to adjust their upset sensitivity and an enlarged sensitive drain junction to enhance the cell's upset rate.

Description

TECHNICAL FIELD

The present invention relates to bench-level assessment of digital IC sensitivity to SEU upsets and more specifically to a method for predicting the SEU susceptibility of a standard-cell D-latch using an alpha-particle sensitive SRAM, SPICE critical charge simulation results, and alpha-particle interaction physics.

BACKGROUND ART

The trend in electronic components has lead to the development of integrated circuits (ICs) that have higher operating frequencies and require less power per gate. This has been achieved in part by a reduction in the feature sizes that make up the transistors and wires that comprise the ICs. As a consequence, the amount of charge that is required to store a bit of information is about one-tenth of a pico-coulomb (pC), or about 600,000 electrons. Because the amount of charge is so small, cosmic rays can deposit enough charge in an IC to cause a single-event upset (SEU) and alter the state of the memory bit.

The importance of cosmic rays on the performance of integrated circuits in a space environment is evident in the upset rates of various satellites and spacecraft. For example, the Tracking and Data Relay satellite experiences a single event upset per day, which must be corrected from the ground. Such adverse experiences have caused a re-design of spacecraft, such as the Galileo spacecraft. The characterization of digital ICs to heavy ion induced upsets is essential for qualifying their use in critical systems to be used on spacecraft. This characterization usually requires evaluations at an ion source, such as a cyclotron. Such evaluations are time consuming, expensive and error prone. In recent years, laser pulses have been proposed as substitutes for the heavy ion sources. However, the laser simulations are limited by metal layers that frequently block the laser pulses from SEU-sensitive nodes and by the complexity of the ion-photon calibrations. The solution to the single event upset problem continues to be important, as the complexity of spacecraft grows, the size of integrated circuits decrease and as space systems are designed with circuits fabricated at non-radiation hardened foundries.

Thus it would be highly desireable to have a method for bench-level characterization of an IC such as a CMOS standard-cell D-latch, without requiring the use of large and expensive heavy ion sources.

SUMMARY OF THE INVENTION

The purpose of this invention is to provide a technique utilizing test structures to quickly and inexpensively characterize the SEU sensitivity of standard cell latches intended for use in a space environment. This bench-level approach utilizes alpha particles to induce upsets in a companion device, a low LET sensitive 4-k bit test SRAM. This SRAM consists of cells that employ an offset voltage to adjust their upset sensitivity and an enlarged sensitive drain junction to enhance the cell's upset rate.

The invention comprises a method for predicting the SEU susceptibility of a standard-cell D-latch using an alpha-particle sensitive SRAM, SPICE critical charge simulation results, and alpha-particle interaction physics. Measurements were made on a 1.6-.mu.m n-well CMOS 4k-bit test SRAM irradiated with an Am-241 alpha-particle source. A collection depth of 6.09 .mu.m was determined using these results and TRIM computer code. Using this collection depth and SPICE derived critical charge results on the latch design, an LET threshold of 34 Mev cm.sup.2 /mg was predicted. Heavy ion tests were then performed on the latch and an LET threshold of 41 MeV cm.sup.2 /mg was determined.

OBJECTS OF THE INVENTION

It is therefore a principal object of the present invention to provide a method for quickly and inexpensively characterizing the single-event upset (SEU) sensitivity of standard cell latches intended for use in a space environment.

It is another object of the invention to provide a method for predicting the SEU susceptability of a standard-cell D-latch using an alpha-particle sensitive SRAM, critical charge simulation results and alpha-particle interaction physics.

It is still another object of the invention to provide a bench-level technique for characterizing CMOS standard-cell D-latches using alpha-particle sensitive test circuits to obviate heavy ion sources which are expensive and time consuming to use.

BRIEF DESCRIPTION OF THE DRAWINGS

The aforementioned objects and advantages of the present invention, as well as additional objects and advantages thereof will be more fully understood hereinafter as a result of a detailed description of a preferred embodiment when taken in conjunction with the following drawings in which:

FIG. 1 is a schematic drawing of a test SRAM cell;

FIG. 2 is a schematic drawing of a D-latch in the zero state and illustrating reverse biased SEU-sensitive drain diodes;

FIG. 3 is a layout drawing of a standard-cell latch;

FIG. 4 is a graphical illustration of critical charge versus offset voltage for both an SRAM and a latch;

FIG. 5 is a graphical illustration of upset rate versus offset voltage for an SRAM in response to a 4.6 .mu.Ci Am-241 source;

FIG. 6 is a graphical illustration of the relationship between measured upset rate, the spontaneous density function, the source spectra and the SRAM critical charge response;

FIG. 7 is a graphical illustration of a test SRAM differential upset rate response versus offset voltage;

FIG. 8 is a schematic illustration of the path of alpha-particles from an alpha-particle source through the SRAM;

FIG. 9 is a graphical illustration of the D-latch alpha-particle response;

FIG. 10 is a graphical illustration of a heavy ion test performed on a latch for purposes of comparison with the inventive method herein;

FIG. 11 is a graphical illustration of the SRAM critical charge response for Node V2 of FIG. 1 calculated using SPICE-2;

FIG. 12 is a graphical illustration of the number of flipped cells for alpha particles and protons shown relative to the spontaneous bit flip response for a store time of one second;

FIG. 13 is a schematic illustration of the paths of protons from a proton source to the SRAM;

FIG. 14 is a graphical illustration of the SRAM upset probability distribution for alpha particles and protons shown relative to the spontaneous bit flip response;

FIG. 15 is a graphical illustration of charge deposited by 0.55 MeV protons in silicon which is used to calculated the over-layer thickness, .delta.X.sub.3 ; and

FIG. 16 is a graphical illustration of charge deposited by 1.0 MeV protons in silicon which is used to calculated the collection depth, .delta.X.sub.4.

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT

Test SRAM and D-Latch Designs

The test SRAM and latch designs disclosed herein were submitted to the MOS Implementation System (MOSIS) and fabricated at a 1.6-.mu.m n-well double, metal CMOS/bulk foundry.

A schematic diagram of the test SRAM cell is shown in FIG. 1. This figure does not include the read/write transistors. The pulsed current source is used to model an alpha particle strike on drain Dn2 when calculating the critical charge of the cell with SPICE. This cell differs from that of a standard six-transistor SRAM cell in three ways: (1) the source of the p-MOSFET, Mp2, is connected to an adjustable offset voltage, V.sub.o, instead of V.sub.DD to provide a control of the cells critical charge; (2) the drain area of n-MOSFET Mn2, Dn2, has been enlarged by a factor of four over minimum to enhance upset rates, thus reducing measurement time; and (3) the cell is imbalanced by widening Mn2 over minimum to enhance its SEU sensitivity versus V.sub.o response. The dimensions of the MOSFETs and their drains is given in Table 1. The bloated drain Dn2, is fabricated in a p-type substrate doped to 5.times.10.sup.14 cm.sup.3. The Dn2 junction has enhanced charge collection due to funneling compared to charge collected by the p-MOSFET drain formed in the n-well, which truncates the ion-induced plasma track.

TABLE 1______________________________________DEVICE L (.mu.m) W (.mu.m) AD (.mu.m.sup.2)______________________________________Mn1 1.6 2.4 17.92Mn2 1.6 3.2 74.88Mp1 3.2 2.4 12.16Mp2 3.2 2.4 12.16______________________________________

In operation all the memory cells are written into a "sensitive" state where Mn2 is turned ON, connecting V.sub.o to the bloated drain, Dn2. V.sub.o is then lowered from V.sub.DD =5 V for a period called the stare time. Thereafter V.sub.o is returned to V.sub.DD and the cells are read to determine the number of upset cells. This cycle is repeated at different values of V.sub.o and can be repeated a number of times at a given V.sub.o to improve the resolution of the measurement.

The latch results presented herein were obtained from a 64-cell modified transparent D-latch array fabricated through the same foundry as the test SRAM. In FIG. 2 a schematic diagram of this latch is shown with the latch in the zero state where reverse biased SEU-sensitive drain diodes are represented by said squares. The offset voltage, V.sub.o was added to sensitize the latch to alpha particles. The source of Mp2 is tied to V.sub.DD in the unmodified latch. The diodes designated with solid squares collect the most charge when struck by an ionized particle because they sit in the p-substrate. The layout of the latch cell is shown in FIG. 3. This layout includes an inverter that generates enable-bar, and inverters that buffer the V1-node and the V2-node. The sensitive diode areas at the V2-node and V3-node that are designated with solid squares in FIG. 2, are outlined with a bold line in FIG. 3. These diodes each have an area of 69 .mu.m.sup.2.

Spice Analysis

The critical charge, Q.sub.c of the test SRAM cell and the D-latch cell were determined as a function of V.sub.o using MOSIS supplied parameters and SPICE. The parasitic nodal capacitances were modeled by fixed metal and polysilicon interconnect capacitances and by the drain depletion capacitances given by their areas and peripheries and by the SPICE LEVEL 2 junction capacitance parameters. A triangle current pulse with a 1:19 rise:fall shape was used to upset these cells. For a given pulse height, the transient simulation was evaluated out to 100 ns where the response was compared to V.sub.DD /2 to determine if the cell had flipped. The pulse height was adjusted using the binary search algorithm until the difference in charge (area under current pulse) between successive simulation runs differed by less than 1 fC. The resulting critical charge versus V.sub.o response for the V2-node of the test SRAM cell and for the average of the V2-node and V3-node for the D-latch is shown in FIG. 4. These results were found to be invariant with current pulse widths up to 500 ps. Since the alpha particle pulse width is about 200 ps, the response of these circuits exceeds that of the alpha particle pulse width.

As V.sub.o decreases the critical charge of the SRAM cell and D-latch decreases approximately linearly and goes to zero at V.sub.o =V.sub.os. For the SRAM cell V.sub.osS =1.8 V and for the latch V.sub.osL =2.4 V. DC SPICE simulations of these bistable circuits show that the sensitive state and metastable state approach each other as V.sub.o decreases until they are equal at V.sub.o =V.sub.os. For V.sub.o .ltoreq.V.sub.os the circuit spontaneously goes to the flipped state. The circuit at this point is no longer bistable. V.sub.os is a function of the threshold voltages and geometries of the n- and p-MOSFETs. Process induced dispersion in the MOSFET parameters gives rise to a dispersion in V.sub.os for an array of cells. The slopes of the critical charge curves in FIG. 4 are defined by their upset capacitances, C.sub.U =.delta.Q.sub.c /.delta.V.sub.o. The upset capacitance of the SRAM cell is C.sub.US =56 fC/V and of the latch cell is C.sub.UL =850 fC/V. SPICE is a well-known computer program for simulating and analyzing circuits which, for example, is described in detail in the text entitled "SPICE: A Guide To Circuit Simulation And Analysis Using PSPICE", second edition by Paul W. Twimenga, copyright 1992, Prentice Hall.

Test SRAM Alpha-Particle Response

The response of the SRAM to alpha particles from a 4.6 .mu.Ci Am-241 source placed 0.508 cm above the SRAM in air is shown in FIG. 5. Also shown, is the spontaneous response which was measured without the alpha-particle source present. As V.sub.o is lowered from V.sub.o, the SRAM begins to upset at V.sub.o <3.2 V in a region termed the tail region. As V.sub.o is further lowered, the upset rate increases more slowly; this region is termed the periphery hit region. Finally, as V.sub.o reaches V.sub.os =1.8 V, the upset rate increases very rapidly in what is termed the spontaneous flip region.

In the periphery hit region, cells upset when the alpha particle hit the Dn2 drain or come close to, but miss, the Dn2 drain. This effect can be explained by a delayed field-funneling effect and to diffused charge that is collected by Dn2. Since the cell critical charge decreases as V.sub.o decreases, the alpha particles can strike further from the Dn2 depletion edge and still upset the cell.

The upset rate R, displayed in FIG. 5, was obtained by dividing the average measured number of upsets N.sub.a that occurred in a stare time t by t, or R=N.sub.a /t, where ##EQU1## and the number of samples, k=200. The stare time in this case was fixed at 2 seconds. The spontaneous upset data was also divided by t so that both responses saturate at N.sub.t /t=2048, where N.sub.t is the number of cells in the memory cell (4096 in this case). The response of the SRAM cells is given by the detector equation: ##EQU2## whose homogeneous solution with initial condition at t=0 of no tripped cells (N=0) is:

N=N.sub.t (1-.epsilon..sup.-.sigma..PHI.) (3) where .sigma. is the device cross section (Dn2 drain area), and .PHI. is the total integrated flux whose particles deposit charge exceeding the critical charge Q.sub.c required to trip a cell.

The upset capacitance C.sub.U, the spontaneous upset voltage V.sub.os, and the hole-electron pair charge-energy factor for silicon K=44.2 fC/MeV enable one to calibrate the offset voltage V.sub.o axis in terms of the critical energy E.sub.c, by

E.sub.c =(C.sub.U /K)(V.sub.o -V.sub.os), (4) or for this SRAM, E.sub.c =1.267(V.sub.o -1.8) MeV.

FIG. 6 illustrates the relationship between upset rate R, the SRAM spontaneous density n, the source spectra .PHI., and the SPICE derived critical charge curve. The spontaneous density, n=dR/dV.sub.o, is the derivative of the spontaneous response data shown in FIG. 5, and was found to be Gaussian with a mean and standard deviation of V.sub.os =1.8 .+-.0.011 V. The area under the n(V.sub.o) curve is equal to N.sub.t. The exponential periphery hit response is not shown in FIG. 6 because it was subtracted from the measured R. In the analysis to follow, the spontaneous density, n, is approximated as an impulse function. This is justified for the SRAM data shown in FIG. 5 because the slope of n/t=4096 cells/34 mV sec is much steeper than the SRAM alpha-particle response.

The tail region response shown in FIG. 5 was smoothed and differentiated using a twenty-five point least squares difference algorithm. This differential response, shown in FIG. 7, is composed of a number of peaks plus the onset of the peripheral hit response. These peaks are due to the selective degradation of the alpha-particle energies as the particles pass through the various regions found in the bloated drain Dn2. The major peak, shown in FIG. 7, is due to the alpha particles that pass through the thickest regions; those covered with metal 2 (M2). This M2 region is 38.2 .mu.m.sup.2. In the region covered with M2, the alpha particles lose more energy than in other regions so that they have the highest linear energy transfer (LET).

A schematic view of the path of the alpha particle through the silicon is shown in FIG. 8. The thickness X.sub.4 -X.sub.2 is obtained by first converting the full-width-half-maximum of the peak in the response curve shown in FIG. 7, V.sub.of, to the voltage straggle, V.sub.o.sigma., of the alpha particles at the end of the collection region in the silicon using V.sub.o.sigma. =V.sub.of /2.35. Then the voltage straggle is converted to energy straggle using

E.sub..sigma. =(C.sub.u /K)V.sub.o.sigma. (5) Finally, X.sub.4 -X.sub.2 =13.58 .mu.m was calculated using the range-straggle tables for alpha particles in silicon.

The silicon collection layer thickness .delta.X.sub.4 was determined from the offset voltage shift .delta.V.sub.o =V.sub.op -V.sub.os =0.98 V. This was converted to an energy shift using

.delta.E=(C.sub.u /K)V.sub.o (6) .delta.X.sub.4 =6.09 .mu.m was calculated using the energy range tables for alpha particles in silicon. Finally the over-layer thickness, .delta.X.sub.3 was calculated by .delta.X.sub.3 =(X.sub.4 -X.sub.2)-.delta.X.sub.4 =7.59 .mu.m.

D-Latch Alpha Particle Response

The alpha-particle response of the latch array was obtained in the same manner as for the SRAM except that the Am-241 source was placed 0.558 cm above the latch chip in air. The results, shown in FIG. 9, include the spontaneous data obtained before and after exposing the part to alpha particles. The shift observed in these responses is due to total dose effects such as the charging of gate oxides. The average of the pre- and post-radiation spontaneous data was used in the analysis of this part.

The dispersion in the latch alpha-particle response, shown in FIG. 9, is primarily that of the spontaneous response because the offset voltage dispersion due to the alpha-particle energy spectrum is much smaller than the dispersion in the spontaneous response. This is a consequence of the reduced energy detection sensitivity of the latch (K/C.sub.UL =0.052 V/MeV for the latch versus K/C.sub.US =0.789 V/MeV for the SRAM). The dispersion in the latch alpha-particle response, V.sub.os =34 mV, was found by fitting the data to a Gaussian function. As shown in FIG. 9 the shift in the offset voltage, .delta.V.sub.o, was obtained between the points on the average spontaneous curve and the alpha-particle curve where, on the average, half of the latch cells are flipping (N.sub.t /2=32). Using the SPICE derived upset capacitance for the latch, C.sub.UL =850 fC/V, the collected charge required to upset a latch cell is .delta.Q.sub.cL =50.2 .+-.28.9 fC which agrees with that obtained with the SRAM, Q.sub.cS. The dispersion in Q.sub.cL, however, is large compared to the dispersion in Q.sub.cS. Since Q.sub.cL =Q.sub.cS, it is valid to apply the collection depth dX.sub.4 and the over layer thickness obtained with the SRAM to the analysis of the latch.

For the latch, the LET threshold for V.sub.DD =5 V operation can be predicted from the SPICE derived upset capacitance, C.sub.UL, the measured mean spontaneous offset voltage, V.sub.os, and the collection depth obtained from the SRAM alpha-particle response, .delta.X.sub.4 : ##EQU3## where, C.sub.UL =850 fC/V, V.sub.os =2.49 V, K=44.2 fC/MeV, .rho.=2320 mg/cm.sup.3 for Si, and .delta.X.sub.4 =6.09.times.10.sup.-4 cm.

D-Latch Heavy Ion Response

Heavy ion tests were performed at the Brookhaven National Laboratory on the latch array with V.sub.DD =V.sub.o =5 V. These measurements were made by gating the ion beam on and then off when a fluence, F=1.13.times.10.sup.7 ions/cm.sup.2 cos .theta. was reached, where .theta. is the angle of incidence. The stare time was set so that an average of two upsets would occur each test cycle. This assured that essentially all N.sub.t =64 latch cells in the array were susceptible to upset and allowed the effective cross section per cell to be calculated by: ##EQU4## where N is the total number of upsets that occur while the ion beam is gated on.

Using the average over layer thickness, .delta.X.sub.3 =4.3 .mu.m, and the collection depth, .delta.X.sub.4 =6.1 .mu.m obtained from the SRAM alpha particle measurements, the energies of the ions used in the test were computed from the TRIM at the entry point of the Si collection layer, X.sub.3, and at the exit point of the Si collection layer, X.sub.4. The average LET of the ion in the collection layer was then calculated by: LET=[E(X.sub.3)-E(X.sub.4)].rho..delta.X.sub.4. These results are shown in Table 2 and in FIG. 9.

TABLE 2______________________________________ PLOT ANGLE, .theta. .sigma. LETION SYMBOL (degrees) (.mu.m.sup.2) (MeV/cm.sup.2 mg)______________________________________I-127 .quadrature. 0 232 58.8Br-79 45 209 57.3Br-79 .largecircle. 35 144 49.5Br-79 0 63.8 40.3Fe-56 X 45 33.0 37.5______________________________________

A test SRAM, designed to be upset-sensitive (LET soft) and having an electrically adjustable critical charge has been characterized using alpha particles using an inexpensive laboratory setup. The SRAM alpha particle results and SPICE simulation results facilitated in predicting the LET threshold of a standard-cell D-latch. The predicted LET=34 MeV cm.sup.2 /mg compares favorably with heavy ion test value of 41 MeV cm.sup.2 /mg.

Preferred Analysis Method

A preferred method for calibration of the SRAM detector will now be described in conjunction with FIGS. 11 to 16. The procedure requires determining the detector's overlayer thickness, .delta.X.sub.3, and collection depth, .delta.X.sub.4. These were determined using 0.55 and 1.0 MeV protons from the Caltech Tandem Van de Graaff. The analysis requires knowledge of the proton charge deposition versus range and this was calculated using particle range physics from TRIM. The validity of the analysis was verified with energy straggling measurements. Finally, the SPICE circuit simulation program was used to compute the relationship between the charge deposited in the memory cell and the cell offset voltage.

The critical charge, Q.sub.c, of the SRAM cell was determined as a function of V.sub.o using MOSIS supplied SPICE parameters. The parasitic nodal capacitances were modeled by fixed metal and polysilicon interconnect capacitances and by the drain depletion capacitances using their areas and peripheries. The SPICE simulation used the level-2-model. A triangle current pulse with a 1:19 rise:fall shape was used to upset these cells. The location of the pulse generator in the memory cell is shown in FIG. 1. For a given pulse height, the transient simulation was examined at 100 ns where the response was compared to V.sub.DD /2 to determine if the cell had flipped. The current pulse height was adjusted using a binary search algorithm until the difference in charge (area under current pulse) between successive simulation runs differed by less than 1 fC. These results were found to be invariant with current pulse widths up to 500 ps. Since the proton width is about 200 ps, the response of these circuits exceeds that of the proton current pulse.

SPICE simulations provide the Q.sub.c versus V.sub.o curve shown in FIG. 11 which is well approximated by the straight line relationship:

Q.sub.c =Cu(V.sub.o -V.sub.os.mu.) where V.sub.os.mu. =1.8 V is the mean offset voltage in the spontaneous flip range and the slope is the upset capacitance, C.sub.u =dQ.sub.c /dV.sub.o =56 fC/V.

In the SRAM test sequence, all memory cells are written into the sensitive or stare state where Mn2 is turned OFF and Mp2 is turned ON which connects V.sub.o to the bloated drain, Dn2; see FIG. 1. V.sub.o is then lowered from 5 V for a period called the stare time which can last a few seconds to several days. Thereafter, V.sub.o is returned to 5 V and the cells read to determine the number of upsets. This cycle is repeated for different values of V.sub.o.

Test results are shown in FIG. 12 for protons and alpha particles and compared against the spontaneous flip response. The overall behavior indicates that at high V.sub.o the response curve is determined by particle energy straggling. At intermediate V.sub.o, the response is determined by the collection of particles outside the area of drain Dn2. This is called the peripheral hit region of the curve. Finally, at low V.sub.o, the cells flip spontaneously.

Layer Thickness Analysis

The following analysis is used to determine the SRAM detector Dn2 overlayer material thickness, .delta.X.sub.3, and charge collection region thickness, .delta.X.sub.4. A schematic view of these layers is shown in FIG. 13.

The analysis requires the determination of offset voltage, V.sub.op, at the peak of the particle upset distribution. In this technique the particle upset distribution is Gaussianly distributed as seen in the cumulative probability plot shown in FIG. 14. The analysis follows from the SRAM detector equation:

dN/dt=.sigma..phi.(N.sub.t -N) where .sigma. is the area of Dn2,.phi. is the flux, N is the measured number of flipped cells, and N.sub.t is the total number of cells in the SRAM. Evaluating for N.fwdarw.0 dN/dt N.fwdarw.0=R.sub.o =.sigma..phi.N.sub.t

The particle flux at diode, Dn2, is described by the Gaussian or normal distribution expressed by the error function: ##EQU5## where .phi..sub.m is the maximum particle flux entering the SRAM, V.sub.op is the offset voltage in the vacinity of the upset peak, V.sub.op.mu. is the mean or peak value, and V.sub.op.sigma. is the standard deviation for the upset distribution. Note that .phi.=.phi..sub.m /2 at V.sub.op =V.sub.op.mu. which states that only half the particle flux entering the SRAM can upset cells.

Combining the above equations leads to the following expression for the number of flipped cells in the vacinity of the particle upset distribution: ##EQU6##

The probability distribution given in percent and plotted in FIG. 14 was determined from:

P(V.sub.oi >V.sub.op)=100.multidot.(N-0.5)/N.sub.t where N is the number of flipped cells at V.sub.op and for this SRAM N.sub.t =4096 cells.

The calculation of N can take two approaches. The first approach is used when the experimenter knows the beam flux, Q.sub.m. In this case, N.tbd.R.sub.o /.sigma..phi..sub.m where .sigma.=68.8 .mu.m for Dn2. The second approach is used when the experimenter knows the beam fluence, F.sub.m, and the number of upsets, N.sub.f, observed during the time the beam is on. That is N=N.sub.f /.sigma.F.sub.m. This approach is useful when the flux is variable and the fluence can be monitored.

The data points in FIG. 14 allow an analysis of the Gaussian nature of the energy dispersion of the particles as they lose energy in the silicon. At the peak of the upset distribution, at the 50 percent point, only half of the particles can deposit sufficient charge to flip the cells. This is defined as Q.sub.cp and is given by:

Q.sub.cp =C.sub.u (V.sub.op.mu. -V.sub.os.mu.) where Q.sub.cp values for the four curves shown in FIG. 14 are listed in Table 3.

The SRAM overlayer thickness, .delta.X.sub.3, was determined using a 0.55 MeV proton beam which stopped within the collection layer. This is crucial for it means that all the charge deposited in the collection layer, .delta.X.sub.4, is collected by the diode Dn2. The charge deposited by the 0.55 MeV proton beam is shown in FIG. 15. This curve was plotted using TRIM. The analysis for .delta.X.sub.4 begins by determining V.sub.op.mu. for the 0.55 MeV proton upset data shown in FIG. 14; the results are listed in Table 3. Then Q.sub.cp is determined from FIG. 11. Finally, the Q.sub.cp value is subtracted from the End of Range, EOR, of the charge deposition curve shown in FIG. 15. Thus .delta.X.sub.3 =4.32 .mu.m.

The charge collection thickness, .delta.X.sub.4, was measured with 1.0 MeV proton beam which had a range greater than .delta.X.sub.3 +.delta.X.sub.4. The charge deposited by the 1.0 MeV proton beam is shown in FIG. 16. The analysis begins by determining V.sub.op.mu. for the 1.0 MeV proton upset data shown in FIG. 14. Then Q.sub.cp was determined from FIG. 11. Finally .delta.X.sub.4 was determined by adding Q.sub.cp to the overlayer charge as determined from .delta.X.sub.3 ; the technique is illustrated in FIG. 16. Thus .delta.X.sub.4 =6.64 .mu.m.

The charge collection depth, .delta.X.sub.4, was also determined for a 4.7 MeV alpha particle and a collection depth of 6.33 .mu.m was determined. This value is slightly smaller than the collection layer thickness determined for the 1.0 MeV proton. This is due to the fact that heavy particles have a shorter range.

The linear energy transfer, LET, in MeV.multidot.cm.sup.2 /mg for the particles can now be calculated using:

LET=C.sub.u (V.sub.op.mu. -V.sub.os.mu.)/K.sub..rho. .delta.X.sub.4 where for silicon K=44.2 fC/MeV, .rho.=2320 mg/cm.sup.3, C.sub.u is in fC/V, V.sub.o is in Volts, and .delta.X.sub.4 is in cm. LET values are listed in Table 3. The LET value for V.sub.op.mu. =5 V is 2.88 MeV cm.sup.2 /mg for an estimated .delta.X.sub.4 =6.00 .mu.m. The conversion factor K is determined from: K=1.6.times.10.sup.-19 (Coul./electron).times.10.sup.15 (fCoul./Coul.).times.10.sup.6 (eV/MeV)/3.62(eV/eh-pair).

Layer Thickness Dispersion Analysis

The result of calculating the errors in .delta.X.sub.3 and .delta.X.sub.4, that is, .delta.X.sub..sigma. and .delta.X.sub..sigma., is shown in Table 3. It shows that as the 0.55 and 1.0 MeV protons pass through the various regions that the energy dispersion increases. This leads to an uncertainty in the thickness of each layer. The thickness dispersion is calculated from the energy dispersion. For the over-layer the depth dispersion is:

.delta.X.sub.3.sigma. =.delta.E.sub.3 X.sub..sigma. .multidot.K.multidot.d.delta.X.sub.3 /dQ.sub.d where .delta.E.sub.3 X.sub..sigma. is the energy dispersion due to over-layer thickness variations and from FIG. 15, dQ.sub.d /d.delta.X.sub.3 =3.85 fC/.mu.m. The collection layer depth dispersion is: .delta.X.sub.4.sigma. =.delta.E.sub.4 X.sub..sigma. .multidot.K.multidot.d.delta.X.sub.4 /dQ.sub.d where .delta.E.sub.4 X.sub..sigma. is the energy dispersion due to collection depth variations and from FIG. 16, dQ.sub.d /d.delta.X.sub.4 =3.125 fC/.mu.m.

The evaluation of .delta.X.sub.3.sigma. requires an evaluation of .delta.E.sub.3.sigma.. Since the 0.55 MeV protons stop in the charge collection region, the energy dispersion is determined by the Au scattering foil and the Si over-layer and not the Si collection layer. The energy dispersion for the 0.55 MeV protons is evaluated at the peak of the probability distribution as defined at the 50 percent point in FIG. 14, it is calculated from E.sub.op.sigma. (0.55)=(C.sub.u /K)V.sub.op.sigma. (0.55). This energy dispersion consists of the following components:

E.sub.op.sigma. (0.55)=[(E.sub.os.sigma.).sup.2 +(E.sub.l.sigma.).sup.2 +(.delta.E.sub.3.sigma.).sup.2 ].sup.1/2 where E.sub.os.sigma. =(C.sub.u /K)V.sub.os.sigma. is the instrument function energy dispersion, E.sub.l.sigma. is the source Au scattering foil energy dispersion, and .delta.E.sub.3.sigma. is .delta.E.sub.3.sigma. =[(.delta..sub.E.sub.3 S.sub..sigma.).sup.2 +(.delta.E.sub.3 X.sub..sigma.).sup.2 ].sup.1/2 where .delta.E.sub.3 S.sub..sigma. is over-layer energy dispersion due to particle straggling and .delta.E.sub.3 X.sub..sigma. is the over-layer energy dispersion due to over-layer thickness fluctuations. The E.sub.l.sigma. =0.012 MeV was determined from Rutherford scattering theory for a 1.2-.mu.m thick Au foil. .delta.E.sub.3 S.sub..sigma. =0.009 MeV was determined from Rutherford scattering theory for a 4.32-.mu.m thick over layer. .delta.E.sub.3 X.sub..sigma. is now calculated from the above two equations and .delta.X.sub.3.sigma. =0.19 .mu.m was then determined.

The evaluation of .delta.X.sub.4.sigma. requires an evaluation of .delta.E.sub.4.sigma.. Since the 1.0 MeV protons pass through the charge collection region, the energy dispersion is determined by the Au scattering foil, over-layer and collection layer. The energy dispersion for the 1.0 MeV protons is evaluated at the peak of the probability distribution as defined at the 50 percent point in FIG. 14; it is calculated from E.sub.op.sigma. (1.0)=(C.sub.u /K) V.sub.op.sigma. (1.0). This dispersion consists of the following components:

E.sub.op.sigma. (1.0)=[(E.sub.os.sigma.).sup.2 +(E.sub.l.sigma.).sup.2 +(.delta.E.sub.3.sigma.).sup.2 +(.delta.E.sub.4.sigma.).sup.2 ].sup.1/2 where .delta.E.sub.4.sigma. is .delta.E.sub.4.sigma. =[(.delta.E.sub.4 S.sub..sigma.).sup.2 +(.delta.E.sub.4 X.sub..sigma.).sup.2 ].sup.1/2 where .delta.E.sub.4 S.sub..sigma. is collection layer energy dispersion due to particle straggling and .delta.E.sub.4 X.sub..sigma. is the energy dispersion due to collection layer thickness fluctuations. The .delta.E.sub.4 S.sub..sigma. =0.011 MeV was determined from Rutherford scattering theory for a 6.64-.mu.m thick collection layer. .delta.E.sub.4 X.sub..sigma. =0.025 MeV is now calculated from the above two equations and .delta.X.sub.4.sigma. =0.36 .mu.m was then determined. TABLE 3______________________________________.delta.X.sub.3, .delta.X.sub.4, and LET Values LET Q.sub.cp V.sub.op.mu. V.sub.op.sigma. .delta.X3 .delta.X4 MeV .multidot.PARTICLE fC V V .mu.m .mu.m cm.sup.2 /mg______________________________________0.55 MeV 10.47 1.987 0.021 4.32 -- --PROTON1.00 MeV 16.18 2.089 0.030 4.32 6.64 0.24PROTON4.70 MeV 50.40 2.66 0.143 4.32 6.33 0.77ALPHA-- 178.00 5.00 -- 4.32 6.00 2.88SPONTANEOUS V.sub.os.mu. V.sub.os.sigma. 1.80 0.012______________________________________

Claims

1. A method of predicting the SEU susceptibility of a SRAM; the method comprising the steps of:

a) providing an enhanced SEU sensitivity test SRAM having an electrically adjustable critical charge;

b) determining the critical charge as a function of offset voltage to ascertain the upset capacitance C.sub.US of said test SRAM;

c) exposing the test SRAM to a first known level of alpha particle energy to ascertain the overlayer thickness of said SRAM, said first known level being selected to stop the alpha particle energy in the collection layer of said SRAM;

d) calculating the overlayer thickness .delta.X.sub.3 of said test SRAM based upon the charge deposited by said first known level;

e) exposing the test SRAM to a second known level of alpha particle energy to ascertain the charge collection thickness of said SRAM, said second known level being selected so alpha particle energy stops beyond the collection layer;

f) calculating the collection layer thickness .delta.X.sub.4 of said test SRAM based upon charge deposited by said second known level and the calculated thickness of said overlayer; and

g) calculating the linear energy transfer (LET) for said test SRAM: ##EQU7## where V.sub.op.mu. =the peak value of offset voltage at maximum particle upset distribution;

V.sub.os.mu. =the mean offset voltage in the spontaneous flip range;

K=the hole-election pair charge-energy factor for silicon;

.rho.=the density of silicon;

C.sub.US =upset capacitance for SRAM.

2. The method recited in claim 1 wherein in step a) said SRAM is made more sensitive to incident-ion-induced SEU by providing at least one pull-down field effect transistor having a bloated drain surface area and at least one pull-up field effect transistor having a source connected to an offset voltage.

3. The method recited in claim 1 wherein said offset voltage is determined by measuring the offset voltage at which about one-half of the total number of latch cells are upset spontaneously without deliberate ion exposure.