1 Field of the Invention
The present invention relates to a method of analyzing mass spectra obtained by a multi-turn time-of-flight mass spectrometer.
2. Description of Related Art
The resolution R of a time-of-flight mass spectrometer is given by
where dT is the time width of each ion peak and T is the flight time. Accordingly, if the time width dT is constant, the resolution can be increased by increasing the flight time T. To achieve this object, multi-turn time-of-flight mass spectrometers in which increased flight distances are used to improve the resolution R have been proposed in the field of mass spectrometry using time-of-flight mass spectrometers.
A multi-turn time-of-flight mass spectrometer can increase the flight distance of ions while avoiding increase in size of the instrument. Therefore, the small instrument can obtain mass spectra with high mass resolution. To increase the flight time T, it is necessary to increase the flight distance. This, in turn, increases the size of the instrument. Accordingly, it has been contemplated to increase the flight distance by making ions rotate along the same trajectory many times.
Multi-turn time-of-flight mass spectrometers are disclosed in the following publications: Japanese Patent Laid-Open No. H11-135060; Japanese Patent Laid-Open No. H11-135061; and Japanese Patent Laid-Open No. H11-195398.
With this instrument, however, as ions make many revolutions, faster ions overrun slower ions on the circulating trajectory. After this phenomenon, peaks of ions having undergone different numbers of turns are mixed on the obtained flight time spectrum. This spectrum is hereinafter referred to as a heterogeneous turn number spectrum. On the other hand, an ordinary time-of-flight spectrum is referred to as a single turn number spectrum. In a heterogeneous turn number spectrum, the numbers of turns of ions giving peaks are unknown in cases where there are peaks of unknown masses. Therefore, it is difficult to obtain a mass spectrum from a time-of-flight spectrum.
Accordingly, in order that ions make many turns and a single turn number spectrum be obtained, the velocity range of ions to be measured simultaneously may be restricted to prevent overrunning of ions spatially.
In this method, however, the mass range of ions that can be measured simultaneously is restricted. Where ions in a wide range should be measured, it is necessary to splice together many single turn number spectra. Consequently, the sensitivity is substantially deteriorated. This is detrimental to the measuring time.
In view of the foregoing, it is an object of the present invention to provide a multi-turn time-of-flight mass spectrometer in which the foregoing drawbacks have been alleviated. It is a more specific object of the present invention to provide a method of reconstructing a single turn number spectrum from plural heterogeneous turn number spectra obtained by a multi-turn time-of-flight mass spectrometer.
A method of multi-turn time-of-flight mass analysis for achieving the above-described objects in accordance with the present invention is implemented using a time-of-flight mass spectrometer having a portion for producing ions in a pulsed manner, a multi-turn ion trajectory portion for causing ions to travel in the same circulating trajectory many times, an ejection portion for ejecting the ions out of the ion trajectory portion, and an ion detection portion for detecting the ejected ions. After turning many times in the trajectory portion, the ions are taken out and dispersed according to flight time. Thus, a mass spectrum is obtained. The method of multi-turn time-of-flight mass analysis starts with recording plural heterogeneous turn number spectra F1(t), F2(t), . . . , Fq(t) containing plural ion peaks that might be different in number of turns. The spectra are obtained with different ion residence times, i.e., taken from entry to departure. A definable multi-correlation function G(t) is obtained from the heterogeneous turn number spectra. In this way, a single turn number spectrum is reconstructed. The multi-correlation function G(t) is given by
where Nj(T) (i=1, 2, . . . , q) is an integer determined by the time T (herein also referred to as one turn time) taken for ions to make one turn in the circulating trajectory, yu is an upper limit value of the deviation time from one turn time T, yi is a lower limit value of the deviation time from one turn time T, and H is a function determined by the value of Fj{Nj(t)×T+y} (j=1, 2, . . . , r).
In one feature of the present invention, the range bounded by the limit values yu and yl is wider than a range in which ions of one turn time T are forecast to be detected.
In another feature of the present invention, the limit values yu and yl are functions of one turn time T.
In a further feature of the present invention, if it is assumed that ions of one turn time T are contained in a heterogeneous turn number spectrum Fj(t), the integer Nj(T) is a forecast number of turns of the ions.
In an additional feature of the present invention, the function H provides a calculation of the arithmetic mean of Fj{Nj(t)×T+y} (j=1, 2, . . . , r).
In still another feature of the present invention, the function H provides a calculation of the minimum value of Fj{Nj(t)×T+y} (j=1, 2, . . . , r).
In yet another feature of the present invention, the function H provides a calculation of the geometric mean of Fj{Nj(t)×T+y} (j=1, 2, . . . , r).
In an additional feature of the present invention, the function H provides a calculation of the harmonic mean of Fj{Nj(t)×T+y} (j=1, 2, . . . , r).
Other objects and features of the present invention will appear in the course of the description thereof, which follows.
As can be understood from the description of the present embodiment, “travel in the circulating trajectory” means “circulating travel in a closed trajectory”, such as a circular, elliptical, or figure 8-shaped trajectory as shown in the references identified above. In addition, the “travel in the circulating trajectory” means reciprocation along a non-closed trajectory as shown in
Ions emitted from the pulsed ion source 1 pass through the first ion mirror 2 and enter the circulating trajectory portion (
It is to be understood that this instrument has been described to illustrate the principle of the present invention. The present invention can be applied to any type of multi-turn TOF mass spectrometer having a circular, elliptical, or figure 8-shaped trajectory as shown in the references identified above.
The flight time of ions is now discussed. Let L1 be the effective distance between the pulsed ion source 1 and the incident surface S11 of the first ion mirror 2. Let L2 be the effective distance between the incident surface S11 and the exit surface S22 of the second ion mirror 3. Let S3 be the effective distance from the exit surface S22 to the ion detector 4. Let L4 be the effective flight distance that ions must travel to reciprocate once between the first ion mirror 2 and the second ion mirror 3. The effective distances referred to herein are distances determined in cases where times necessary for acceleration and deceleration are converted into corresponding distances. Therefore, it cannot be said that the effective distances are strictly identical for all ions. However, the effective distances are approximately given by geometrical conditions.
The velocity V of an ion and one turn time T (i.e., period) taken for the ion to make one revolution in the circulating trajectory have the relation
The total flight time S of the ion is a function of the number of turns N and the period T and given by
To detect ions normally after they make N turns in the circulating trajectory, the conditions of the times at which the ion mirrors are turned on and off are now discussed. It is assumed that the ions are emitted at time t=0. Let t1 be the time at which the first ion mirror 2 is turned on. Thus, we have
At this time, the relation V=L4/T holds and so the following relation must hold:
Let t2 be the time at which the second ion mirror 3 is turned off. The following relation must hold:
a(N, T)<t2<b(N, T) (7)
where a (N, T) is the instant of time when the ions have just started the Nth normal trip along the circulating trajectory after leaving the second ion mirror 3, and b (N, T) is the instant of time when the ions are about to enter the second ion mirror 3 after finishing the Nth trip along the circulating trajectory. More specifically, they are given by
where d1 is the effective flight distance within the second ion mirror 3 and d2 is the length of the second ion mirror 3, i.e., the flight distance achieved when the ion mirror voltage is off. Since a(N+1, T) and b(N, T) are not equal, there is an insensitive period in which even ions traveling in the circulating trajectory are not detected normally, depending on the value of the time t2.
A case is discussed in which measurement of a heterogeneous turn number spectrum is repeated q times with varying ion residence time (i.e., from entry to departure). The period from the emission from the pulsed ion source 1 to the instant when the second ion mirror 3 is turned off is set to different values s1, s2, . . . , sq corresponding to the q measurements. The obtained spectra are individually recorded. As a result, q heterogeneous turn number spectra are derived. These spectra are respectively indicated by F1(t), F2(t), . . . , Fq(t).
It is now demonstrated that an ion of velocity V providing the period T is contained in these spectra. The necessary condition that the ion of period T be contained in F1(t) is that an integer N1 capable of satisfying the following inequality exists:
a(N1, T)<s1<b(N1,T) (10)
where s1 is the time at which the second ion mirror 3 is turned off. If this integer N1 is present, it can be identified as a single value. At this time, Eq. (3) indicates that a peak of the ion of the period T is detected near time h1 indicated by
h1=(N1+D)×T (11)
from the heterogeneous turn number spectrum F1(t). The reason why the time is expressed using the expression “near” is that the obtained values of the effective flight distances of L1, L2, and so on are only rough values. Normally, their accurate values are unknown. By applying the same theory to other times sj (j=1, 2, . . . , q), the possibility that peaks of ions of period T are contained can be determined. Also, the number of turns Nj in a case where a peak is detected can be determined.
Let r be the number of heterogeneous turn number spectra which have been judged that they might contain peaks of ions of the period T. Since a subscript is used at will, if a possible spectrum is given by Fj (j=1, 2, . . . , r), the generality of the following discussion is not lost.
At this time, if there are ions of period T, peaks of the ions of interest should be present close to time hj given by
hj=(Nj+D)×T (12)
Let y be the deviation time from the time hj at which a peak should appear. This deviation time y should be common to every ion regardless of the integer Nj, for the following reason. For the same ion, the period T is the same for all trips in the circulating trajectory. The causes of deviation are setting of the effective times, delay in the pulsed circuit, deviation of the time of ion emission, and so on, which are not related to the number of turns. Accordingly, every ion peak is present in F1 (h1+y), F2 (h2+y), . . . , Fr (hr+y). It is difficult, however, to know the deviation time y in advance. However, if the deviation time is varied in appropriate increments, there is an equality y=y0 at which every Fj (hj+y0) has a large value. Conversely, if such y0 is detected, it is demonstrated that ions of period T are present. Note that y0 has a finite spread corresponding to the width of the ion peak. Accordingly, the strength G(T) of an ion of period T is defined in terms of the following multi-correlation function:
where the function H (F1, F2, . . . , Fr) is a function of a variable Fj. Various features can be given only depending on the method of determining the function. The expression is in an integral form. Obviously, in actual calculations, the expression may be replaced by a form of discrete summation. The lower limit value ya and upper limit value yb of the deviation time from the time hj at which a peak originating from an ion of period T appears may be determined according to actual situations. That is, setting the range somewhat wider than the range of possible error of hj is advantageous in terms of the calculation time. In principle, setting the range to a wider range will present no problems. Furthermore, Eq. (13) can be rewritten into the form:
A specific form of the function H (F1, F2, . . . , Fr) is now discussed. The simplest definition of the function H is to take the arithmetic mean of Fj's. That is,
In this case, each Fj has a large value at y=y0 and so G(T) that is an integrated value of the range (Yu<y<yl) containing y0 also has a large value. In this method, however, if Fj happens to have a large value, G(T) is equal to the product of Fj and 1/r, i.e., has a large value. That is, pseudo-peaks will also appear. This complicates the spectrum. Hence, it is highly likely that it is difficult to interpret the spectrum.
To eliminate the pseudo-peaks, the calculation for taking the minimum peak Fj from spectral peaks F1, F2, . . . , Fr may be defined as the function H. In this case, however, the spectral peaks other than Fj are not reflected in the results of calculation. Therefore, there is the disadvantage that other observed data all become useless.
One method of compensating the drawback with the arithmetic averaging method is to take the geometric mean. That is,
H(F1, F2, . . . , Fr)=r{square root}{square root over (F1×F2× . . . ×Fr)} (18)
If the number r increases to a considerably large value, the probability that every Fj happens to have a large value decreases considerably. Accordingly, in cases where y is other than y0, at least one value is quite small. In the geometric averaging method, if any one value becomes null, G(T) takes the form of the product of Fj's and becomes equal to zero irrespective of the values of the other Fj's. Consequently, pseudo-peaks due to accidental coincidence are eliminated. As a result, no pseudo-peaks are finally contained in a single turn number spectrum G(T) reconstructed from plural heterogeneous turn number spectra Fj(t).
If the harmonic mean is taken instead of the geometric mean, the same advantages can be expected as the geometric averaging method. That is,
Also, in this case, if any peak having an intensity of 0 is contained in peak intensities Fj's, the harmonic mean found from the peak is 0. Thus, it is judged as a pseudo-peak. Finally, such pseudo-peaks are not contained in a single turn number spectrum G(T) reconstructed from plural heterogeneous turn number spectra Fj(t).
Besides, various definitions of the function H are conceivable. A definition in which small values are reflected to a larger extent in the function H out of various values of Fj's is desirable. Generally, as the number of turns along the trajectory increases, the ion intensity decreases. Therefore, it is advantageous to normalize the Fj's before the multi-correlation function is taken.
The definition of G(T) has been described so far for the case where the period T is determined. The period T of ions traveling in the circulating trajectory differs according to different mass numbers of ions. Therefore, the value of the period T is varied as a final phase. G(T) is found within a required range. Thus, a single turn-number spectrum is obtained. In this way, a single turn-number spectrum can be reconstructed from plural heterogeneous turn number spectra Fj(t).
This single turn-number spectrum G(T) is a spectrum regarding the time taken for ions to make one reciprocation along the circulating trajectory. The other portions, e.g., the contribution of the time corresponding to D×T, are eliminated.
The principle of the present invention is described using a specific example of simulation. A figure 8-shaped multi-turn time-of-flight mass spectrometer as shown in
When the ions entering the circulating trajectory from the point b make a half revolution along the figure 8-shaped circulating trajectory (over a distance L2 of 1 m), the ions reach the point c where the gate that opens at given timing is mounted. The gate can be opened at arbitrary timing. In this way, ions can be taken out from the circulating trajectory after the ions make a half turn, after 1.5 turns, after 2.5 turns, after 3.5 turns, and after other given number of turns. The flight distance L4 of the ions which go from the point b through the point c and reach the point b again along the figure 8-shaped circulating trajectory (i.e., make one turn) is just twice the distance L2, i.e., 2 m.
The ions taken out of the circulating trajectory from the point c via the gate intact move a distance of L3 (=0.2 m) straightly and are detected by an ion detector placed at a point d.
Four kinds of ions, i.e., lithium ions (m/z=6), aluminum ions (m/z=27), manganese ions (m/z=55), and uranium ions (m/z=238), are assumed as ions traveling in the optical system described above. The timings at which the gate placed at the point c is opened are 13.86 μs, 40.69 μs, and 67.51 μs after the ions leave the point a. At these timings, the aluminum ions reach the point c after making 1.5 turns, 4.5 turns, and 7.5 turns.
When the gate at the point c was opened after a lapse of 13.86 μs since the departure of the ions, the flight times of the ion peaks detected by the ion detector at the point d were four, i.e., 14.75 μs, 15.39 μs, 17.26 μs, and 21.06 μs, respectively.
When the gate at the point c was opened after a lapse of 40.69 μs since emission of ions, the flight times of ion peaks detected by the ion detector at the point d were four, i.e., 41.58 μs, 43.81 μs, 44.89 μs, and 46.53 μs, respectively.
When the gate at the point c was opened after a lapse of 67.51 μs since emission of ions, the flight times of ion peaks detected by the ion detector at the point d were four, i.e., 68.41 μs, 70.19 μs, 70.36 μs, and 72.11 μs, respectively.
These are summarized into
The results of the calculation reveal that three peaks of 15.39 μs (a2), 44.89 μs (b3), and 70.19 μs (c2) of the mass spectra (a), (b), and (c) obtained after lapses of 13.86 μs, 40.69 μs, and 67.51 μs, respectively, since ejection of ions represent peaks produced when ions of mass 6 made 3.5 turns, 10.5 turns, and 16.5 turns, respectively, and that three peaks of 14.75 μs (a1), 41.58 μs (b1), and 68.41 μs (c1) represent peaks produced when ions of mass 27 made 1.5 turns, 4.5 turns, and 7.5 turns, respectively. Furthermore, three peaks of 21.06 μs (a4), 46.58 μs (b4), and 72.11 μs (c4) represent peaks produced when ions of mass 55 made 1.5 turns, 3.5 turns, and 5.5 turns, respectively. In addition, three peaks of 17.26 Is (a3), 43.81 its (b2), and 70.36 μs (C3) represent peaks produced when ions of mass 238 made 0.5 turn, 1.5 turns, and 2.5 turns, respectively. The results are listed in FIGS. 6(a), 6(b), and 6(c).
In
Where ions of greatly different masses travel in the same circulating trajectory in this way, lighter ions overrun heavier ions, producing a quite complex mass spectrum. Heterogeneous turn number spectra can be reconstructed into a single turn number spectrum by computationally extracting peaks containing time elements that are N times as large as the period T and normalizing the different turn numbers N.
If peaks due to noise are accidentally introduced into a heterogeneous turn number spectrum, the noise-derived peaks can be eliminated with a quite high probability by taking the geometric mean or harmonic mean of plural peaks obtained from plural heterogeneous turn number spectra. Therefore, noise-derived peaks hardly appear in the finally obtained single turn number spectrum.
The method of analysis of the present invention can be widely used in analysis of heterogeneous turn number spectra obtained by a multi-turn TOF mass spectrometer.
Having thus described my invention with the detail and particularity required by the Patent Laws, what is desired protected by Letters Patent is set forth in the following claims.
Number | Date | Country | Kind |
---|---|---|---|
2003-311317 | Sep 2003 | JP | national |