Continuous monochrometer drift compensation of a spectral monochromator

• US 4,916,645 A
• Filed: 06/02/1988
• Issued: 04/10/1990
• Est. Priority Date: 06/02/1988
• Status: Expired due to Term

• Current Assignee: Perkinelmer Incorporated
• Original Assignee: Perkinelmer Incorporated
• Inventors: Wuest, Charles E., Collins, John B.
• Primary Examiner(s): Lall, P. S.
• Assistant Examiner(s): TRANS, VINCENT N
• Application Number:
  US07/201,250
• Time in Patent Office: 677 Days
• Field of Search: 364/571.04, 364/571.01, 364/571.02, 364/525, 364/559, 364/167.01, 73/1 R, 356/319, 356/308, 356/326, 356/331, 356/332, 356/334
• US Class Current: 356/319
• CPC Class Codes: G01J 2003/063 Step motor
• CPC Class Codes: G01J 2003/066 Microprocessor control of functions, e.g. slit, scan, bandwidth during scan
• CPC Class Codes: G01J 3/18 using diffraction elements, e.g. grating gratings per se G02B
• CPC Class Codes: G01J 3/28 Investigating the spectrum using colour filters G01J3/51
• CPC Class Codes: G01J 3/32 Investigating bands of a spectrum in sequence by a single detector
• CPC Class Codes: G01N 21/73 using plasma burners or torches
• CPC Class Codes: G01N 21/74 using flameless atomising, e.g. graphite furnaces

Abstract

Iterative compensation of drift of peak positions of spectral lines is effected in a spectral monochromator including a grating, a detector of spectral fractions of a spectral band, a stepper motor for varying relative orientation of the grating and the detector, and a computer. Computer-defined spectral windows each encopasses a spectral band and has an initial spectral center. Each window is scanned to determine a peak spectral position. Calculation is made for determining a spectral position error of the peak position from the initial center for each corresponding window. A functional average of the offsets for the peaks is calculated as a linear function of window position, and a revised spectral center for each window is calculated as being equal to the initial center plus the functional average for the window position determined from the linear function. Each window is shifted correspondingly. The step of successively scanning through each window is repeated to determine a new peak position for each corresponding band, whereby each new peak position is maintained near the spectral center of each corresponding window.

Claims

1. A method of continuously compensating for drift of peak positions of spectral lines in a spectral monochromator including radiation means for generating elemental radiation for atomic elements, optical dispersion means receptive of the radiation for producing a series of spectral bands each having a peak spectral position associated with maximum radiation intensity, detector means selectively receptive of a spectral function of a spectral band for producing an intensity signal representative of the intensity of the spectral fraction being received, scanning means responsive of a scanning signal for varying relative orientation of the dispersion means and the detector means such that the series of spectral bands is scanned across the detector means in spectral positions corresponding to successive spectral fractions, and computer means receptive of the intensity signal for generating the scanning signal and calculating spectral position of each peak;

the method comprising the steps of, in sequence;

Establishing a series of spectral windows each encompassing a selected spectral band and having an initial window center, successively scanning through each of the windows so as to determine a maximum peak spectral position for a corresponding band, calculating a spectral position error of the peak position from the initial center for each corresponding window, calculating a drift offset equal to a functional average of the position errors, calculating a revised spectral center for each window equal to the initial center plus the corresponding drift offset, shifting each window to its revised center, and repeating the step of successively scanning through each of the windows so as to determine a new peak position for each corresponding band, whereby each new peak position is maintained near a spectral center of each corresponding window.

2. A method according to claim 1 wherein the dispersion means includes a dispersion member and the scanning means includes a stepper motor receptive of the scanning signal and operatively connected to the dispersion member.

3. A method according to claim 1 further comprising successively repeating in sequence the steps of calculating a spectral position error of the peak position from the initial center for each window, calculating a drift offset, calculating a revised spectral center, shifting each window, and repeating the step of successively scanning through each of the windows, such as to determine a further new peak position for each corresponding band, whereby each further new peak position is maintained near the spectral center of each corresponding window for successive spectral scans.

4. A method according to claim 1 further comprising the step of utilizing each position error for calculating the drift offset only if such offset meets a predetermined criterion.

5. A method according to claim 4 wherein the step of utilizing each position error comprises measuring a peak intensity for each peak and utilizing each position error for calculating the drift offset only if such peak intensity associated with such position error exceeds a predetermined threshold intensity.

6. A method according to claim 5 wherein the threshold intensity for each peak is determined as equal to a predetermined threshold percentage of a standard sample peak intensity.

7. A method according to claim 4 wherein the step of utilizing each position error comprises calculating a curvature parameter for each peak, and utilizing the corresponding position error for calculating the simple function only if the corresponding curvature parameter meets a predetermined curvature criterion.

8. A method according to claim 7 further comprising determining for each window a parabolic line computed from the parabolic equation Y=AX₂ +BX+C where A, B and C are parabolic parameters, X represents successive spectral positions of spectral fractions in the corresponding window, Y represents the intensity of each corresponding spectral fraction, and the computation of the parabolic line is made by fitting intensities for at least three intensity points consisting of the peak and adjacent spectral fractions adjacent to and on either side of the peak, and wherein the curvature parameter is defined as parameter A and is less than a predetermined curvature value.

9. A method according to claim 8 wherein each of the intensity points has an intensity error from the parabola, and the method further comprises calculating an error ratio defined as the sum of the squares of the intensity errors divided by the curvature, and utilizing the corresponding offset for calculating the simple function only if the corresponding error ratio exceeds a predetermined error ratio.

10. A method according to claim 4 wherein the step of utilizing each position error comprises detecting each peak that is in a spectral fraction abutting an edge of a window, measuring a peak intensity for each such abutting spectral fraction, measuring corresponding intensities of at least two spectral fractions adjacent the abutting spectral fraction, and utilizing a corresponding position error for calculating the simple function only if the abutting intensity exceeds a predetermined threshold intensity and if the adjacent spectral fractions have successively decreasing intensities away from the peak abutting fraction.

11. A method according to claim 1 wherein the functional average is a simple function of nominal window position.

12. A method according to claim 11 wherein the step of calculating a drift offset comprises computing a linear function of offsets vs. nominal window position, and the functional average for each nominal window is determined from the linear function.

13. A method according to claim 11 further comprising calculating a standard deviation for the simple function, rejecting offsets that exceed a preselected multiple of the standard deviation, and recalculating the simple function without the rejected offsets prior to shifting each window.

14. A method according to claim 11 further comprising additional steps of calculating a merit value for the simple function and shifting each window only if the merit value meets a predetermined criterion.

15. A method according to claim 14 wherein the additional steps comprise calculating a correlation coefficient for the simple function, and shifting each window only if the correlation coefficient exceeds a predetermined correlation coefficient.