The "median" method for the reduction of noise and trigger jitter on waveform data.

The "median" method for the reduction of the effect of noise and trigger jitter on waveform data is described. The effectiveness of this method was examined using simulations and experiments and for typical jitter and noise observed in electrical pulse metrology, is shown to provide reconstructed waveforms with transition durations that accurately match those of the input signal. Also, for aberrations, an upper bound on the error in the amplitude of the aberration is provided.

Key words: aberration; median; median method; noise; overshoot; transition duration; trigger jitter; waveform.

1. Introduction

The effect of trigger jitter on the measurement of signals is to add signal-dependent noise and to reduce the measurement bandwidth or, equivalently, to slow the step response of the measurement instrumentation. This reduction or slowing results in an increase in the transition duration, [t.sub.sys], of the step response of the measurement system. [t.sub.sys] can be approximated using:

[t.sub.sys] = [square root of ([t.sub.instr.sup.2] + [t.sub.jit.sup.2])], (1)

where [t.sub.instr] is the transition duration of the step response of the measurement instrument and [t.sub.jit] is the transition duration of the cumulative jitter distribution. The relationship between [t.sub.jit] and the measured value of the rms jitter, [[sigma].sub.jit], can be found either by calculation and tables [1] or by waveform simulation and extraction of waveform parameters, both of which yield [t.sub.jit] = 2.56 [[sigma].sub.jit].

If [t.sub.jit] < 0.1 [t.sub.instr], then in most cases [t.sub.jit] can be ignored. However, in high-speed pulse metrology, [t.sub.jit] is often a large fraction of [t.sub.instr]. For example, in some samplers with 3 dB attenuation bandwidths of 80 GHz, [t.sub.jit] [approximately equal to] 1 ps and [t.sub.instr] [approximately equal to] 5 ps. Jitter in this case must be removed to get an accurate estimate of pulse parameters and an accurate reconstruction of the input signal. Since jitter effectively acts as a lowpass filter [2], it can be removed via deconvolution. Although uncertainties associated with the deconvolution process can be assigned to jitter deconvolution [3], it would be advantageous to pulse metrology if an alternative method with lower uncertainties could be developed. Furthermore, the presence of noise increases the uncertainties associated with the deconvolution or waveform reconstruction process [4,5]. Also, typical waveform reconstruction methods employ a regularization operator that comprises a second order difference operator. This difference operator typically causes ringing near waveform transitions [5] and, thereby, may increase the amplitudes of overshoot, undershoot, and other waveform aberrations. Noise in the measurement is reduced by signal averaging; this reduction is proportional to [M.sup.-1/2], where M is the number of waveforms averaged.

Each of the m (m = 1, 2,..., M) waveforms, [f.sub.m][[t.sub.n]] (n = 1, 2,..., N, where N is the number of samples in a waveform) is described to first order by:

[f.sub.m][[t.sub.n]] = [f.sub.0][[t.sub.n]] + [[[partial derivative][f.sub.0]]/[[partial derivative]t]])[.sub.t.sub.n] x [j.sub.m][[t.sub.n]] + [s.sub.m][[t.sub.n]], (2)

where [f.sub.0] is the jitter-free, noise-free waveform, [j.sub.m] is the trigger jitter corresponding to the mth waveform, [s.sub.m] is the noise corresponding to the mth waveform, [t.sub.n] is discrete time with index n, and [[partial derivative][f.sub.0]]/[[partial derivative]t] is the derivative of the signal in the vicinity of [t.sub.n]. The average of Eq. (2), <f[[t.sub.n]][.sub.M]>, is the value typically provided by oscilloscopes or other waveform recorders:

<f[[t.sub.n]]>[.sub.M] = [f.sub.0][[t.sub.n]] + [g.sub.M][[t.sub.n]] + <s[[t.sub.n]]>[.sub.M], (3)

where <s[[t.sub.n]][.sub.M]> is the mean value of the noise at [t.sub.n] and [g.sub.M][[t.sub.n]] is the amplitude contribution from jitter at [t.sub.n]. Based on our observations and the physical cause of the distributions of s and j, both distributions are normal with zero mean and are independent and identically distributed for each [t.sub.n] and each waveform. The two distributions, one caused by jitter and the other by noise, can be dealt with independently. For noise, therefore, we can write:

<s[[t.sub.n]]>[.sub.M] = <s[[t.sub.n]]> = <s>. (4)

For noise and for large M, the median, mode, and mean of <s[[t.sub.n]][.sub.M]> will yield the same value for a symmetric unimodal distribution, which is <s>, where <s> = 0 for the situation described herein. This is not the case for the effect of jitter because the coefficient of j is signal (data) dependent [(see Eq. (2)]. In the case where the mean estimator is used to approximate the signal, jitter will act as a low pass filter [2] thereby reducing the bandwidth of the signal and increasing its transition duration. The median and mode do not involve averaging, and they do not bandlimit the signal as strongly as does the mean. However, as will be shown later, the ability to accurately reconstruct waveform aberrations using the median is limited by the duration of the aberrations relative to [[sigma].sub.jit]. The use of the word aberration and how to compute its values are defined in Ref. [6].

The method we have examined to reduce the effects of noise and trigger jitter on the reconstructed waveform is based on taking the median value for each instant, [t.sub.n], of the set of M waveforms. Therefore, instead of letting the oscilloscope average a set of M waveforms to yield one waveform that is typically displayed on the oscilloscope, a set of M unaveraged waveforms are acquired by a computer to provide a two-dimensional array of M unaveraged waveforms each with N samples. At each [t.sub.n] the median is computed, and the results recorded to yield the reconstructed waveform. We also examined using the mode to remove jitter, but the results for the pulse parameter values showed greater variation than that from using the median. Also the waveform root-sum-of-squares error was greater for the mode than for the median. Although we have implemented the median method on a computer after acquisition of all M waveforms, it could also be implemented in the firmware of the oscilloscope prior to acquisition of the resultant (median reconstructed) waveform, just as signal averaging is presently implemented.

The median method has also been examined by other researchers [7,8]. As they [7,8] both point out, the median method works well for monotonic signals but introduces errors when the signal is not monotonic (that is, exhibits a local extremum). The duration of this "monotonicity," however, need only be long compared to [[sigma].sub.jit] (this is discussed in Sec. 5). The median method causes the magnitudes of the local extrema to be underestimated. In Sec. 5, we provide an upper bound for this error and show that it is small compared to the amplitude of the aberration when the duration of the aberration is small compared to [[sigma].sub.jit]. In [8], the author opts for a jitter deconvolution method rather than the median to improve the accuracy of the reconstructed spectrum (the higher frequency components are increased). However, the nature of the regularization operator in waveform reconstructions using deconvolution will cause a false increase in the magnitude of the high-frequency components simply because of the spurious ringing introduced near transitions (see Sec. 5). The magnitude of this ringing will depend on, for example, the accuracy of the estimate of the jitter function and signal noise. An algorithm for removing noise and jitter from waveforms using spline approximations was also examined [9]. No pulse parameter values were provided in [9], consequently, a quantitative comparison was not possible. However, by comparing figures of waveforms generated with similar conditions of noise and jitter relative to pulse amplitude and transition duration (or pulse duration), it appears the waveforms extracted using the median method exhibit less departure from the target values in the vicinities of aberrations than does the spline method.

2. Simulation Results of Ideal Step Waveforms

We examined the median method for jitter and noise reduction by using simulations and experiments. In the first set of simulations, we used a jittered ideal step with additive noise and examined the error in the reconstructed waveforms as a function of noise and jitter. The additive noise and trigger jitter were both modeled as having a zero-mean normal distribution. The amplitude, [A.sub.p], of the step was 1 amplitude unit (au); the rms additive noise, [[sigma].sub.n], as a percentage of [A.sub.p], ranged between 0% and 50% (or, 0 au to 0.5 au); and the rms jitter, in units of a sampling interval (si), ranged from 0 si to 50 si. Figure 1 shows an example of one set of 128 jittered noisy waveforms (for the set shown, [[sigma].sub.jit] = 50 si, [[sigma].sub.n] = 0.01 [A.sub.p]). Figure 2 shows the waveform reconstructed using the median method, the input (ideal step) waveform, and the waveform that is the result of the average of the waveforms shown in Fig. 1. As can be seen from Fig. 2, the reconstructed step closely resembles the input step even for large [[sigma].sub.jit]. Figure 3 shows the rms error value of the median reconstructed steps as a function of noise and trigger jitter. The rms error is given by:

[FIGURE 1 OMITTED]

e = {[1/K][K.summation over (k=1)][square root of ([1/N][N.summation over (n=1)]([s.sub.n,k] - [s.sub.n,ideal])[.sup.2])]}[.sup.1/2], (5)

where K is the number of iterations (100 in this case), [s.sub.n] is the value of the step at index n, [s.sub.n,ideal] is the value of ideal step at index n, and N is the number of samples (1000 in this case). Since the ideal step had a 50% duty factor and a low and high state value of 0 au and 1 au, its rms value is 0.707 au. Comparing this value to the rms error shown in Fig. 3, it can be seen that the median reconstruction method significantly reduces noise and jitter with respect to the input waveforms.

3. Simulation Results of Step-like Waveforms

In the second set of simulations, we used the same ideal step but filtered it with a Butterworth filter (order = 3, low-frequency-pass cut off = 0.125/si, and sampling frequency = 2/si) to give the signal features (such as aberrations) typically observed in pulse waveforms. The pulse parameters of transition duration (between the instants corresponding to 10% and 90% of pulse amplitude), [t.sub.10-90], and overshoot, OS, were compared. The process of determining the pulse parameters starts with obtaining the bimodal histogram of the waveform amplitude values, per methods described in [6]. The methods to compute transition duration and overshoot are given in [6].

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

The noise-free waveform had a [t.sub.10-90] = 5.86 si and an OS of 0.08 au. Overshoot was examined because it is the aberration with the largest amplitude and would be the most obviously sensitive to jitter. The duration of the overshoot (between its 50% reference level instants, as defined in Ref. [6]) was approximately 1 si. Figures 4 and 5 show the effect of trigger jitter and noise on the transition duration of the averaged and median reconstructed waveforms. As can be seen, the transition duration of the averaged waveform is very sensitive to jitter and noise. The median-reconstructed waveform provides waveforms that accurately reproduce the transition duration of the input signal for [[sigma].sub.jit] [less than or equal to] 10 si and [[sigma].sub.n] [less than or equal to] 0.1 au. Figures 6 and 7 show the effect of trigger jitter and noise on the overshoot of the averaged and reconstructed waveforms. As can be seen, the overshoot of the averaged waveform is very sensitive to jitter and noise. The overshoot in the median reconstructed waveform is also sensitive to jitter and noise. The large error in overshoot for large jitter ([[sigma].sub.jit] > 1 si) is caused by the fact that both the positive and negative transitions of the overshoot contribute to the set of data values in the vicinity of the overshoot. This effect is more clearly shown in Fig. 8 where different median-reconstructed waveforms are shown along with the averaged waveforms. In this figure, the duration, [T.sub.OS], (between 50% reference levels) of the overshoot is approximately 5 si and the transition duration of the cumulative jitter distribution is 2.56 [[sigma].sub.jit]. The results of Fig. 8 indicate that integrity of the reconstructed waveform is better for [[sigma].sub.jit] < [T.sub.OS]. OS in the median reconstructed waveform is not sensitive to noise for [[sigma].sub.n] [less than or equal to] 0.1 au. The primary purpose for showing Figs. 4 and 6 is as a reference for the effects of averaging on overshoot and transition duration in the presence of noise and trigger jitter.

[FIGURE 4 OMITTED]

Figures 9 and 10 show the effect of the number of waveforms, M, used in the simulation to compute the median for [[sigma].sub.jit] = 2 si and 4 si. The generated waveforms were noise free. Each datum in Figs. 9 and 10 is the result of repeating the simulation 100 times for each set of M waveforms. The transition duration values show a monotonic decrease in value, for both levels of jitter, to the target value of 5.86 si. The overshoot values, on the other hand, do not exhibit a common monotonic change approaching the target value of 0.08 au. For [[sigma].sub.jit] [greater than or equal to] 4 si, increasing M does not result in the correct value of 0.08 au for overshoot. As discussed earlier, this is because the duration of the overshoot is less than the transition duration of the cumulative jitter distribution. For [[sigma].sub.jit] = 2 si, on the other hand, the target value of 0.08 au for overshoot is obtained. Figures 9 and 10 also show that M [greater than or equal to] 512 is required to achieve convergence to the correct value. Furthermore, the variation in the values of overshoot and transition duration are significant for M < 128.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

4. Experimental Results

Figures 11 to 18 show the results of experiments using two different sampling oscilloscopes (samplers), sampler A and sampler B, one each from a different manufacturer with different levels of jitter. The jitter was varied by adjusting the trigger level on the oscilloscope. All trigger jitter values were obtained using the oscilloscope firmware. The 3 dB attenuation bandwidths for both samplers are nominally 50 GHz. The pulses produced by the pulse generator have a nominal transition duration of 12 ps. The sampling intervals were 2 ps. There were M = 128 waveforms in each set of waveforms.

Figure 11 shows the results of 128 waveforms acquired using sampler A with a jitter level of 9.5 ps rms. The transition duration of the cumulative jitter distribution is about 24.3 ps. Figure 12 shows the oscilloscope-averaged waveform, the mean computed from each of the waveforms shown in Fig. 11, and the waveform reconstructed using the median method. Table 1 shows the transition duration values computed for these three cases. Figures 13 and 14 show the same results as Figs. 11 and 12 except for an rms jitter value of 1.5 ps. The transition duration of the oscilloscope-averaged waveform was not equal to that obtained from the computer-averaged waveform. This may in part be explained by the fact that the waveforms represented by the oscilloscope-averaged waveform were not the same as those used for the computer-averaged waveform. This discrepancy may not have been as large if the number (128) of waveforms to generate the averaged waveforms shown in Figs. 11-14 was larger. Figures 9 and 10 show that around 500 waveforms should be averaged to converge to a consistent value. However, for the low jitter case, the differences between the transition durations of the computer- and scope-averaged waveforms are within the uncertainties of the measurement [3]. The differences in the transition durations of the median-reconstructed waveforms for both levels of jitter are within measurement uncertainty.

Figures 15 through 18 and Table 2 show the same results for Sampler B as Figs. 11 through 14 and Table 1 did for Sampler A, except with different levels of trigger jitter. The same observation made for Sampler A and the low jitter case is the same for Sampler B (see Tables 1 and 2). For the large jitter case for Sampler B, however, the transition duration of oscilloscope-averaged waveforms was greater than that of the computer-averaged waveforms. Also, the transition durations of the median-reconstructed waveforms were faster for Sampler A than for Sampler B, indicating that the bandwidth of Sampler A is greater than that of Sampler B. This is corroborated by the low jitter results (see Tables 1 and 2 first row) where it can be approximated that the contribution of jitter to [t.sub.sys] is negligible and, given the same pulse source, the difference in transition duration is caused by the difference in [t.sub.instr] of the samplers.

Table 3 shows the results of limiting the number, M, of waveforms used to compute the median on the value of transition duration. This is analogous to the simulation results shown in Fig. 10 except that only one set of M waveforms is used in Table 3 whereas 100 sets of M waveforms were used to produce the results shown in Fig. 10. Table 3, however, still shows the importance of selecting an appropriate value of M, namely one that is adequate to obtain reproducible results.

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

[FIGURE 13 OMITTED]

[FIGURE 14 OMITTED]

[FIGURE 15 OMITTED]

[FIGURE 16 OMITTED]

[FIGURE 17 OMITTED]

[FIGURE 18 OMITTED]

5. Aberration Amplitudes

As mentioned earlier, a waveform reconstructed using the median method will exhibit an underestimate of the magnitude of aberrations. Typically, however, it is difficult to estimate the uncertainty in overshoot (OS) and undershoot (US) (both terms are defined in Ref. [6]) of the reconstructed waveform because OS and US in this waveform are affected by a variety of input signal and measurement instrument characteristics [3]. Consequently, we developed an empirical method to estimate the uncertainty in OS and US of the median reconstructed waveforms using the OS and US values from the acquired waveform. This uncertainty is estimated by the difference, [d.sub.max], between the jitter-caused attenuation of a sinewave and that of a constant signal. The sinewave has a half-period equal to the duration, [T.sub.abb], of the OS or US from the acquired waveform. [d.sub.max] is given by:

[d.sub.max] = [[V.sub.abb][[integral].sub.t=0.sup.k[T.sub.PD]](1-cos[2[pi][f.sub.abb]t])dt]/[[[integral].sub.t=0.sup.k[T.sub.PD]]dt] = [V.sub.abb][[k[T.sub.PD] - [[sin(2[pi][f.sub.abb]k[T.sub.PD]]/[2[pi][f.sub.abb]]]]/[k[T.sub.PD]]] = [V.sub.abb][1-[[sin(2[pi][f.sub.abb]k[T.sub.PD])]/[2[pi][f.sub.abb]k[T.sub.PD]]]], (6)

where [T.sub.PD] is the width of the probability density function of the jitter, [V.sub.abb] is the amplitude of the aberration, k is an empirically-determined factor adjusted to ensure [d.sub.max] is greater than the measured error (in situations where it is possible to measure the error), [f.sub.abb] = 0.5/[T.sub.abb], and [T.sub.abb] is the duration of the aberration. As a check for Eq. (6), if [T.sub.PD] = 0 (no jitter), then [d.sub.max] = 0 (no uncertainty contribution from the jitter), as would be expected. On the other hand, if [T.sub.PD] = N/[T.sub.abb], where N is an integer, then [d.sub.max] would approach [V.sub.abb] as N increased. This also makes sense because as N increases, jitter would tend to cause the aberration to disappear in the observed waveform. The expanded uncertainties for OS and US should include [d.sub.max] as well as contributions from other sources [3]. Simulations were performed that yielded 25 sets of data for four different levels of noise and different values of [T.sub.abb]. Each set of data was the result of the average of 256 waveforms. As shown in Table 4, the results from the analysis of this simulated data show that Eq. (6) provides an overestimate of the errors in OS values. With respect to fidelity of reconstruction of the aberrations in the waveform, "monotonicity" is satisfied if [[sigma].sub.jit] [less than or equal to] 0.2 [T.sub.abb]. The value of k = 4 was used in Table 4 to ensure that uncertainty in OS for the noisier signals would be less than [d.sub.max]. This method for finding a limit in the OS or US uncertainty when [[sigma].sub.jit]/[T.sub.abb] = 2 did not work because this amount of jitter effectively removed the aberration from the observed waveform.

A comparison of waveforms reconstructed using the median method and jitter deconvolution are shown in Figs. 19-21. Figures 19 and 20 show waveforms reconstructed using the median and jitter deconvolution methods respectively. These figures show that for large jitter, the median method provides a much more realistic reconstruction of the input signal than does the jitter deconvolution method. This is corroborated by the distorted spectra caused by the jitter deconvolution method (see Fig. 21). For small levels of jitter, [t.sub.jit] [less than or equal to] [t.sub.instr], the jitter deconvolution and median methods appear to work equally well for removal of jitter from the waveforms.

[FIGURE 19 OMITTED]

[FIGURE 20 OMITTED]

[FIGURE 21 OMITTED]

6. Discussion

The simulation results of Secs. 2 and 3 and experimental results of Sec. 4 can be compared by multiplying the simulation horizontal axis values and corresponding waveform parameters by two and then changing the simulation horizontal axis units from si to ps. Accordingly, the jitter of 2 si becomes 4 ps and the simulated transition duration of 5.86 si becomes 11.72 ps, which corresponds to the transition duration of typical high-speed pulse generators. The simulated trigger jitter of 4 ps is at least three times larger than what is typically observed. Consequently, the difficulty in reconstructing the aberrations of a waveform that was experienced in the simulation would not be experienced in reconstructing actual signals used in waveform metrology. In general, the simulation and measurement results show that the median method for reducing the effects of trigger jitter and noise on the acquired waveform is effective. However, to accurately reproduce waveform aberrations, the duration of the aberrations should not be greater than the transition duration of the cumulative jitter distribution.

7. References

[1] Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, M. Abramowitz and I. A. Stegun, eds., National Bureau of Standards, U.S. Department of Commerce, U.S. Government Printing Office, Washington, DC, USA, June 1964.

[2] W. L. Gans, The measurement and deconvolution of time jitter in equivalent-time waveform samplers, IEEE Trans. Instrum. Meas. IM-32, 126-133 (1983).

[3] N. G. Paulter and D. R. Larson, Pulse parameter uncertainty analysis, Metrologia 93, 143-155 (2002).

[4] N. S. Nahman and M. E. Guillaume, Deconvolution of time domain waveforms in the presence of noise, National Bureau of Standards Technical Note 1047, U.S. Department of Commerce, Washington, DC (1981).

[5] N. G. Paulter, A causal regularizing deconvolution filter for optimal waveform reconstruction, IEEE Trans. Instrum. Meas. 43, 740-747 (1994).

[6] IEEE Standards on Transitions, Pulses, and Related Waveforms, Std 181-2003, The Institute of Electrical and Electronics Engineers, New York, NY, USA, July 2003.

[7] T. M. Souders, D. R. Flach, C. Hagwood, and G. L. Yang, The effect of timing jitter in sampling systems, IEEE Trans. Instrum. Meas. 39, 80-85 (1990).

[8] J. Verspecht, Compensation of timing jitter-induced distortion of sampled waveforms, IEEE Trans. Instrum. Meas. 43, 726-732 (1994).

[9] M. G. Cox, P. M. Harris, and D. A. Humphreys, An algorithm for removal of noise and jitter in signals and its application to picosecond electrical measurement, Numerical Algorithms 5, 491-508 (1993).

About the authors: N. G. Paulter is an electronics engineer in the Quantum Electrical Metrology Division of the NIST Electronics and Electrical Engineering Laboratory and leader of the Pulse Metrology and Time Domain Measurements project. D. R. Larson is an electronics engineer in the Quantum Electrical Metrology Division and Calibration Leader for the 65XXX series of high-speed-pulse-related special tests. The National Institute of Standards and Technology is an agency of the Technology Administration, U.S. Department of Commerce.

N. G. Paulter and D. R. Larson

National Institute of Standards and Technology, Gaithersburg, MD 20899-0001

[e-mail]

[e-mail]

Accepted: July 18, 2005

Available online: http://www.nist.gov/jres

Table 1. Transition duration values computed for waveforms taken from
Sampler A

             Averaged,   Averaged,
Jitter (ps)  scope (ps)  computer (ps)  Median (ps)

1.5          12.9        12.7           12.2
9.5          23.1        26.3           12.6

Table 2. Transition duration values computed for waveforms taken from
Sampler B

             Averaged,   Averaged,
Jitter (ps)  scope (ps)  computer (ps)  Median (ps)

0.9          13.3        13.4           13.3
4.7          16.8        15.6           13.2

Table 3. Transition duration values for waveforms acquired using
Samplers A and B. [t.sub.meas] is the measured transition duration

     Sampler ([[sigma].sub.jit])
     A (1.5 ps)         A (9.5 ps)
M    [t.sub.meas] (ps)  [t.sub.meas] (ps)

  2  11.4               10.7
  4  12.2               19.1
  8  12.5                6.4
 16  11.8                9.5
 32  12.4                9.3
 64  12.5                9.4
128  12.2               12.6

     Sampler ([[sigma].sub.jit])
     B (0.9 ps)         B (4.7 ps)
M    [t.sub.meas] (ps)  [t.sub.meas] (ps)

  2  11.5               10.3
  4  13.0                9.7
  8  13.1                8.4
 16  13.1               11.2
 32  13.4               11.3
 64  13.3               13.0
128  13.3               13.2

Table 4. The effect of [[sigma].sub.jit] on OS for a given [T.sub.abb]
with and without noise. The four rows in each cell of the table are the
OS value obtained from the simulation, [d.sub.max]; the estimated value,
O[S.sub.sm], plus [d.sub.max]; and the standard deviation of the
observed values, [[sigma].sub.OS]. For the data in this table, k = 4 and
[T.sub.abb] = 5.293 si. The bottom row shows the value of OS for the
noise-free, jitter-free input signal

                          [[sigma].sub.jit]/[T.sub.abb]
[[sigma].sub.n](au)       0.1                  0.2

0       O[S.sub.sm]       0.05747              0.05747
        [d.sub.max]       0.00530              0.02057
        O[S.sub.sm]+      0.06277              0.07804
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2x          [+ or -] 2x[10.sup.-12]
                            [10.sup.-13]
0.001   O[S.sub.sm]       0.05748              0.05744
        [d.sub.max]       0.00524              0.01992
        O[S.sub.sm]+      0.06272              0.07736
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2.7x        [+ or -] [1.610.sup.-4]
                            [10.sup.-5]
0.01    O[S.sub.sm]       0.05755              0.04688
        [d.sub.max]       0.00505              0.01697
        O[S.sub.sm]+      0.06260              0.06385
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2x          [+ or -] [7.410.sup.-4]
                            [10.sup.-4]
0.1     O[S.sub.sm]       0.05576              0.05055
        [d.sub.max]       0.00516              0.02262
        O[S.sub.sm]+      0.06082              0.07317
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 9x          [+ or -] 2x[10.sup.-3]
                            [10.sup.-3]
Input                     0.05474              0.05474
  OS
  (au)

                          [[sigma].sub.jit]/[T.sub.abb]
                          OS (an)
[[sigma].sub.n] (au)      0.5                  1.0

0       O[S.sub.sm]       0.05747              0.01472
        [d.sub.max]       0.02057              0.06558
        O[S.sub.sm]+      0.07804              0.08030
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2x          [+ or -] 1.3x[10.sup.-2]
                            [10.sup.-12]
0.001   O[S.sub.sm]       0.05744              0.01472
        [d.sub.max]       0.01992              0.06558
        O[S.sub.sm]+      0.07736              0.08030
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 1.7x        [+ or -] 2.2x[10.sup.-2]
                            [10.sup.-4]
0.01    O[S.sub.sm]       0.04688              0.01634
        [d.sub.max]       0.01647              0.05896
        O[S.sub.sm]+      0.06335              0.07530
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 7.4x        [+ or -] 1.6x[10.sup.-2]
                            [10.sup.-4]
0.1     O[S.sub.sm]       0.05055              0.0367
        [d.sub.max]       0.02262              0.1351
        O[S.sub.sm]+      0.07317              0.1718
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2x          [+ or -] 1.2x[10.sup.-2]
                            [10.sup.-2]
Input                     0.05474              0.05474
  OS
  (au)

                          [[sigma].sub.jit]/[T.sub.abb]
[[sigma].sub.n] (au)      2.0

0       O[S.sub.sm]       0.00495
        [d.sub.max]       0.00001
        O[S.sub.sm]+      0.00496
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 1.7x[10.sup.-7]
0.001   O[S.sub.sm]       0.00495
        [d.sub.max]       0.00001
        O[S.sub.sm]+      0.00496
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 2.4x[10.sup.-7]
0.01    O[S.sub.sm]       0.00498
        [d.sub.max]       0.00930
        O[S.sub.sm]+      0.01428
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 1x[10.sup.-2]
0.1     O[S.sub.sm]       0.02595
        [d.sub.max]       0.00280
        O[S.sub.sm]+      0.02575
          [d.sub.max]
        [[sigma].sub.OS]  [+ or -] 7x[10.sup.-3]
Input                     0.05474
  OS
  (au)