WO2025004999A1 - Data processing device, mass spectrometry system, and mass spectrometry data processing method - Google Patents
Data processing device, mass spectrometry system, and mass spectrometry data processing method Download PDFInfo
- Publication number
- WO2025004999A1 WO2025004999A1 PCT/JP2024/022568 JP2024022568W WO2025004999A1 WO 2025004999 A1 WO2025004999 A1 WO 2025004999A1 JP 2024022568 W JP2024022568 W JP 2024022568W WO 2025004999 A1 WO2025004999 A1 WO 2025004999A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- calibration
- voltage
- quadrupole
- ions
- processing unit
- Prior art date
Links
- 238000012545 processing Methods 0.000 title claims abstract description 118
- 238000004949 mass spectrometry Methods 0.000 title claims abstract description 31
- 238000003672 processing method Methods 0.000 title claims description 8
- 150000002500 ions Chemical class 0.000 claims abstract description 172
- 238000012937 correction Methods 0.000 claims abstract description 75
- 238000005259 measurement Methods 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims description 47
- 238000004458 analytical method Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 description 48
- 239000012535 impurity Substances 0.000 description 36
- 238000005086 pumping Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 14
- 239000000523 sample Substances 0.000 description 12
- 230000006870 function Effects 0.000 description 7
- 239000000538 analytical sample Substances 0.000 description 6
- 239000000126 substance Substances 0.000 description 6
- 239000012634 fragment Substances 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010494 dissociation reaction Methods 0.000 description 1
- 230000005593 dissociations Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000005040 ion trap Methods 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/26—Mass spectrometers or separator tubes
- H01J49/34—Dynamic spectrometers
- H01J49/40—Time-of-flight spectrometers
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/26—Mass spectrometers or separator tubes
- H01J49/34—Dynamic spectrometers
- H01J49/42—Stability-of-path spectrometers, e.g. monopole, quadrupole, multipole, farvitrons
Definitions
- the present invention relates to a data processing device, a mass spectrometry system, and a mass spectrometry data processing method.
- ions with a known m/z are measured and the relationship between a specific parameter (the voltage applied to the electrodes in a quadrupole mass filter) and m/z is determined.
- a specific parameter the voltage applied to the electrodes in a quadrupole mass filter
- m/z is determined by using such techniques.
- Patent Document 1 discloses a mass analysis control device, a mass analysis device, a mass analysis control method, and a mass analysis method, in which "an ionization control unit 1 controls an ionization unit 21 so as to ionize the analytical sample and the mass calibration sample into analytical sample ions and mass calibration sample ions.
- An ion dissociation control unit 2 controls an ion trap 23 so as to trap the analytical sample ions and mass calibration sample ions and selectively dissociate the mass calibration sample ions into a plurality of calibration fragment ions without selecting the analytical sample ions as precursor ions.
- An analysis control unit 3 controls a mass analyzer 24 so as to perform mass analysis of the analytical sample ions and the plurality of calibration fragment ions.
- the mass-to-charge ratio of the analytical sample ions is calibrated by a calibration unit 4 based on the mass-to-charge ratios of the plurality of calibration fragment ions" (see abstract).
- Patent document 2 discloses a method for measuring mass spectra under multiple conditions with different ionization polarities, etc., and performing mass axis calibration if there is a difference from the reference data.
- the problem that this invention aims to solve is to determine an accurate mass axis when performing mass axis calibration, even when there are impurity ions whose m/z values are adjacent to the calibration point.
- the present invention was made in light of this background, and its objective is to enable highly accurate mass analysis.
- the present invention comprises a plot processing unit that plots a voltage obtained as a result of measuring a plurality of calibration ions, whose m/z is known, using a mass spectrometer, and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes, and a correction processing unit that corrects the calibration points by correcting the outliers if there are any outliers for each of the calibration points.
- a plot processing unit that plots a voltage obtained as a result of measuring a plurality of calibration ions, whose m/z is known, using a mass spectrometer, and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes
- a correction processing unit that corrects the calibration points by correcting the outliers if there are any outliers for each of the calibration points.
- the present invention makes it possible to perform mass analysis with high accuracy.
- FIG. 1 is a diagram illustrating an example of the configuration of a mass spectrometry system according to a first embodiment.
- FIG. 2 illustrates an example of a hardware configuration of a data processing device.
- FIG. 1 is a diagram showing the configuration of a quadrupole mass filter (part 1).
- FIG. 2 is a diagram showing the configuration of a quadrupole mass filter (part 2).
- FIG. 13 illustrates the operation of a quadrupole mass filter and the relationship between the quadrupole RF voltage and the quadrupole DC voltage.
- FIG. 1 is a graph showing the relationship between the quadrupole RF voltage and the ion signal intensity.
- FIG. 1 is a graph showing the relationship between m/z and quadrupole RF voltage (part 1).
- FIG. 1 illustrates an example of a hardware configuration of a data processing device.
- FIG. 1 is a diagram showing the configuration of a quadrupole mass filter (part 1).
- FIG. 2 is a diagram showing
- FIG. 2 is a graph showing the relationship between m/z and quadrupole RF voltage (part 2).
- FIG. 3 is a graph showing the relationship between m/z and quadrupole RF voltage (part 3).
- 3 is a flowchart showing the procedure of a mass spectrometry data processing method according to the first embodiment.
- FIG. 13 is a table showing specific values of calibration points.
- FIG. 13 is a graph showing deviations.
- FIG. 13 is a diagram illustrating a voltage calibration technique.
- 10 is a flowchart showing the procedure of a mass spectrometry data processing method according to a second embodiment.
- FIG. 4 shows a graph of the relationship between m/z and quadrupole RF voltage.
- FIG. 13 is a diagram illustrating an example of the configuration of a mass spectrometry system according to a third embodiment.
- FIG. 1 is a diagram showing an example of the configuration of a mass spectrometry system 1 according to the first embodiment.
- the mass spectrometry system 1 includes a mass spectrometry device 100, a voltage control device 300, DC power supplies 301 and 303, and an RF power supply 302.
- the mass spectrometry system 1 further includes a data processing device 200 to which an output device 201 is connected.
- the dashed lines connecting the voltage control device 300 and the DC power supplies 301 and 303, the RF power supply 302, and the data processing device 200, and the dashed line connecting the detector 152 and the data processing device 200 indicate signal lines.
- the solid line connecting the DC power supply 303 and the ion guide 130, and the solid lines connecting the RF power supply 302 and the quadrupole mass filter 140 and the DC power supply 301 indicate electric wires to which voltages are applied.
- the mass spectrometer 100 includes an ion source 151 , a first differential pumping section 101 , a second differential pumping section 102 , and an analyzing section 103 .
- Ions generated in the ion source 151 are introduced into the first differential pumping section 101 through the aperture 121.
- the ion source 151 is composed of an electrospray ion source, an atmospheric pressure chemical ion source, an atmospheric pressure photoion source, an atmospheric pressure matrix-assisted laser desorption ion source, or the like, and operates at atmospheric pressure or at a vacuum lower than atmospheric pressure.
- the first differential pumping section 101 is evacuated by a pump 111. As a result, the first differential pumping section 101 is maintained at a vacuum level of 10 Pa to 500 Pa. Ions that pass through the first differential pumping section 101 are introduced into the second differential pumping section 102 through an aperture 122.
- the second differential pumping section 102 is evacuated by a pump 112. As a result, the second differential pumping section 102 is maintained at a vacuum level of 0.1 Pa to 10 Pa.
- An ion guide 130 that focuses ions is installed in the second differential pumping section 102. Ions focused by the ion guide 130 are introduced through an aperture 123 into the analysis section 103 in which a quadrupole mass filter 140 is installed.
- the ion guide 130 is supplied with a static voltage (hereinafter referred to as DC voltage) from a DC power supply 303.
- the RF power supply 302 applies a composite wave of the radio frequency voltage (hereinafter referred to as RF voltage) generated by the RF power supply 302 and the DC voltage sent from the DC power supply 301 to the quadrupole mass filter 140.
- the ion guide 130 and the quadrupole mass filter 140 are connected via a dielectric 153 such as a capacitor.
- RF voltage is supplied to the quadrupole rod electrodes 141 (see FIGS. 3A and 3B) of the quadrupole mass filter 140 from an RF power supply 302 controlled by a voltage control device 300.
- the high-frequency component of the RF voltage is then supplied from the quadrupole rod electrodes 141 to the ion guide rod electrodes of the ion guide 130 through the dielectric 153.
- the number of power supplies can be reduced compared to a configuration in which RF voltage is supplied separately to the ion guide 130 and the quadrupole mass filter 140, making the mass spectrometry system 1 less expensive.
- the analytical section 103 is evacuated by the pump 113. This keeps the analytical section 103 at a pressure of 1E-3 Pa or less.
- ions are separated according to m/z.
- Ions that pass through the quadrupole mass filter 140 are detected by the detector 152.
- the detector 152 an electron multiplier tube, a type that combines a scintillator and a photomultiplier tube, or a multichannel plate can generally be used.
- the detection intensity of the ions detected by the detector 152 is converted into an electrical signal (output signal) and sent to the data processing device 200.
- the output signal output from the detector 152 is converted into digital data with a fixed sampling period (typically 1 ⁇ s to 1000 ⁇ s) and then sent to the data processing device 200.
- the conversion to digital data is performed by an analog-to-digital converter (ADC) not shown or a pulse counting unit not shown.
- ADC analog-to-digital converter
- the data processing device 200 stores the sent digital data in the memory unit 220.
- a plot processing unit 211 In addition to the storage unit 220, a plot processing unit 211, a correction processing unit 212, and a measurement processing unit 213 are executed.
- the plot processing unit 211 plots the calibration points 521 (FIGS. 6 to 7B) on the m/z-quadrupole RF voltage relationship graphs shown in FIGS. 6 to 7B.
- the correction processing unit 212 corrects the calibration points 521 based on the calibration points 521 plotted by the plot processing unit 211 .
- the measurement processing unit 213 controls the measurement by the mass spectrometer 100 .
- the storage unit 220 is made up of a memory, HD, etc., and can store and hold information such as the data of the spectrum 501 (see FIG. 5), as well as numerical values and relational expressions required for correction.
- the plot processing unit 211, correction processing unit 212, and measurement processing unit 213 each have a calculation function and temporary memory for temporarily storing numerical values required for calculation.
- the data processing device 200 also has a function to control the voltage control device 300 that controls each electrode, etc., and a function to output information to the output device 201.
- the output device 201 is made up of a display, printer, etc.
- the output device 201 outputs information such as the spectrum 501 itself, the m/z corresponding to the peak of the spectrum 501, the signal intensity, and the presence or absence of the substance to be measured.
- the storage unit 220 is mounted on the data processing device 200, but this is not limiting.
- the storage unit 220 may be installed as a device separate from the data processing device 200, such as a database.
- FIG. 2 is a diagram illustrating an example of the hardware configuration of the data processing device 200.
- the data processing device 200 includes a memory 231 including a RAM or the like, a calculation device 232 including a CPU, a GPU, or the like, and a storage device 233 including a HD, an SSD, or the like.
- the data processing device 200 also includes an input device 234 including a keyboard, a mouse, or the like, an output device 201, and a communication device 235.
- the storage device 233 corresponds to the memory unit 220.
- the output device 201 is the output device 201 shown in FIG. 1.
- a program stored in the storage device 233 is loaded into the memory 231.
- the loaded program is then executed by the calculation device 232.
- This embodies the plot processing unit 211, the correction processing unit 212, and the measurement processing unit 213 shown in FIG. 1.
- FIGS. 3A and 3B are diagrams showing the configuration of the quadrupole mass filter 140.
- Fig. 3A shows a perspective view of the quadrupole mass filter 140
- Fig. 3B shows a cross section of the quadrupole mass filter 140 and a diagram showing voltage control for the quadrupole mass filter 140.
- the quadrupole mass filter 140 is composed of four quadrupole rod electrodes 141 (141a to 141d).
- an RF voltage and a DC voltage are applied to each quadrupole rod electrode 141.
- the RF voltage is an AC voltage generated by an RF power supply 302 controlled by a voltage control device 300.
- the DC voltage is a DC voltage generated by a DC power supply 301 controlled by the voltage control device 300.
- a composite wave in which a DC voltage is applied to an RF voltage is applied to each of the quadrupole rod electrodes 141.
- the composite wave of the RF voltage and DC voltage is applied so that it is in phase with each other between adjacent quadrupole rod electrodes 141 and in phase with each other between opposing quadrupole rod electrodes 141.
- RF voltages of opposite phases are applied between the pair of quadrupole rod electrodes 141a, 141c and the pair of quadrupole rod electrodes 141b, 141d.
- the DC voltage is a direct current voltage generated by the DC power supply 301 controlled by the voltage control device 300.
- the DC voltage applied to the quadrupole rod electrodes 141a and 141c is VDC1
- the DC voltage applied to the quadrupole rod electrodes 141b and 141d is -VDC1.
- the RF voltage and DC voltage applied to each quadrupole rod electrode 141 are appropriately referred to as the quadrupole RF voltage and the quadrupole DC voltage.
- the typical voltage amplitude of the quadrupole RF voltage is several hundreds of V to several kV, and the frequency is about 500 kHz to 2 MHz.
- the quadrupole DC voltage is about several tens of V to several hundreds of V.
- FIG. 4 is a diagram illustrating the operation of the quadrupole mass filter 140 and the relationship between the quadrupole RF voltage and the quadrupole DC voltage.
- the m/z range of ions capable of stable orbital motion in the quadrupole mass filter 140 depends on the amplitude of the quadrupole RF voltage and the value of the quadrupole DC voltage. Only ions present inside the stability regions R1 to R3 shown in FIG. 4 can pass through the quadrupole mass filter 140.
- the stability region R1 is the region within the line R1a
- the stability region R2 is the region within the line R2a
- the stability region R3 is the region within the line R3a.
- the stability regions R1 to R3 are different for each m/z of the ions.
- the stability regions R1 to R3 are arranged in the relationship shown in FIG. 4 from ions with small m/z to ions with large m/z. That is, the stability region R1 is the stability region of ions having a certain m/z.
- the stability region R2 is the stability region of ions with a different m/z from the ions in the stability region R1
- the stability region R3 is the stability region of ions with a different m/z from the ions in the stability regions R1 and R2.
- the quadrupole RF voltage and quadrupole DC voltage are set near the apex of the stability region R1-R3 of a certain m/z, it is possible to transmit only ions having that m/z.
- the quadrupole RF voltage is scanned so as to pass near the apex of the stability region R1-R3 of each m/z ion, as shown by the scan line L1 in Figure 4.
- the spectrum 501 shown in Figure 5 can be obtained. In other words, it becomes possible to detect ions having each m/z.
- Fig. 5 is a graph showing the relationship between the quadrupole RF voltage and the ion signal intensity
- Fig. 6 is a graph showing the relationship between the m/z and the quadrupole RF voltage.
- the horizontal axis indicates the quadrupole RF voltage
- the vertical axis indicates the ion signal intensity.
- Fig. 5 also shows a spectrum 501 of calibration ions.
- the calibration ions are ions generated from a calibration sample, and their m/z is known.
- Fig. 5 also shows a spectrum 502 of impurity ions.
- the impurity ions are ions derived from impurities due to contamination or the like.
- the quadrupole RF voltage refers to the amplitude value of the quadrupole RF voltage.
- the quadrupole RF voltage is the voltage obtained as a result of measuring multiple calibration ions with the mass spectrometer 100.
- spectrum 502 of impurity ions is detected in close proximity to spectrum 501a of calibration ions.
- a scan range 511 which is a voltage range.
- the scan range 511 is the scan width of the quadrupole RF voltage performed for ion detection.
- the scan range 511 will be described later.
- the m/z-quadrupole RF voltage relationship graph shown in FIG. 6 is a coordinate system having an m/z axis and a quadrupole RF voltage axis as coordinate axes.
- the m/z-quadrupole RF voltage relationship graph is a coordinate system having voltage and m/z as coordinate axes.
- Calibration point 521 shows the relationship between the quadrupole RF voltage corresponding to the peak value of the ion signal intensity in spectrum 501 shown in FIG. 5 and m/z. As described above, the m/z of the calibration ions is known, so it is possible to plot an m/z-quadrupole RF voltage relationship graph as shown in FIG. 6. Calibration point 521 is a coordinate plot of the voltage obtained as a result of measuring multiple calibration ions with mass spectrometer 100 and the m/z of each calibration ion.
- spectrum 502 due to impurity ions is detected near spectrum 501a due to calibration ions. Furthermore, the peak value of spectrum 501a is lower than the peak value of spectrum 502. Furthermore, in FIG. 5, spectrum 501a and spectrum 502 fall within the same scan range 511. This indicates that spectrum 501a and spectrum 502 are detected as one spectrum 501. This may cause the data processing device 200 to mistake the impurity ions for calibration ions and make a judgment.
- calibration point 521b derived from impurity ions is plotted close to calibration point 521a derived from spectrum 501a.
- the m/z of the impurity ions is linked to the quadrupole RF voltage to which the m/z of the calibration ions that generate spectrum 501a should be linked.
- the m/z of the measured substance is output based on the m/z-quadrupole RF voltage relationship graph as shown in FIG. 6 and the obtained quadrupole RF voltage.
- the m/z of the impurity ions is linked to the quadrupole RF voltage to which the m/z of the calibration ions that produce spectrum 501a should be linked, an incorrect m/z will be output.
- the quadrupole DC voltage is determined so that the range of the spectrum 501 produced by the calibration ions is a specified (constant) value.
- the quadrupole mass filter 140 operates to selectively transmit ions having a specific m/z while keeping the quadrupole RF voltage fixed for a certain period of time.
- the measurement processing unit 213 rescans the quadrupole RF voltage in a specified range and obtains the peak value of the spectrum 501 shown in FIG. 5. Determining the quadrupole RF voltage corresponding to the m/z of the calibration ions based on this peak value is called voltage calibration.
- the first scan is a scan in the presence of impurity ions as described below, and the rescan is a scan after correction of the calibration point 521 (see FIG. 6) described below has been performed.
- a low scan rate is, for example, about 0.1 Da/s to 100 Da/s (1.66 ⁇ 10 ⁇ 25 g/s to 1.66 ⁇ 10 ⁇ 22 g/s).
- a low scan rate it takes a long time to acquire a spectrum 501 in a predetermined range.
- the measurement processing unit 213 scans the entire quadrupole RF voltage range at a low scan rate, it will take an enormous amount of time. Therefore, there is a disadvantage in that the time required for voltage calibration is also long.
- a spectrum 501 is acquired by scanning only the vicinity of the calibration ions at a low scan speed, as shown in the scan range 511 shown by the dashed line in Figure 5. In this way, it is often possible to measure both the accuracy of the voltage calibration and the time required for the voltage calibration.
- FIGS. 7A and 7B are graphs showing the relationship between m/z and quadrupole RF voltage, in which Fig. 7A shows an example in which no impurity ions are detected, and Fig. 7B shows an example in which peaks are confused due to impurity ions.
- the m/z of the calibration ions stably passing through the quadrupole mass filter 140 is proportional to the (amplitude of) the quadrupole RF voltage. Therefore, when no impurity ions are detected, as shown in Fig. 7A, the slope (ai) of the line connecting the calibration points 521 (mi, xi) shown in formula (1) is approximately the same value for all the calibration points 521.
- the slope (ai) is the slope between each of the calibration points 521.
- xi (i1, ..., N) is the value of the quadrupole RF voltage
- i is the number of the calibration point 521.
- the calibration ion (calibration point 521a) is confused with the adjacent impurity ion (calibration point 521b).
- the calibration point 521b of the impurity ion is plotted as an outlier.
- (mk, xk) indicates the coordinates of the calibration point 521 in the m/z-quadrupole RF voltage relationship graph. Due to the confusion between the calibration ion and the adjacent impurity ion, the correct calibration point 521a is recognized as the calibration point 521b of the impurity ion, and the corresponding quadrupole RF voltage is shifted.
- the deviation (Di) is defined by the following formula (2). In this way, the deviation (Di) is the deviation with respect to each slope (ai).
- a_mean is the average of ai defined by the following formula (3).
- the deviation does not have to be the deviation from the average of ai.
- the correction processing unit 212 determines whether a peak has been confused based on whether the standard deviation (SD) defined by the following formula (4) is equal to or greater than a threshold value (Sth). In other words, whether there is an outlier for each calibration point 521 is determined based on whether the standard deviation (SD) defined by the formula (4) is equal to or greater than a threshold value (Sth).
- (flowchart) 8 is a flowchart showing the procedure of the mass spectrometry data processing method according to the first embodiment, with reference to FIGS.
- the plot processing unit 211 performs a plotting process of plotting the measurement results of the calibration ions by the mass spectrometer 100 on an m/z-quadrupole RF voltage relationship graph (S101: first step).
- the plot processing unit 211 performs the plotting based on the quadrupole RF voltage obtained as a result of the scan in the scan range 511 shown in FIG. 5 and the m/z of the calibration ions that is known in advance.
- the correction processing unit 212 calculates (S102) all of the slopes (ai) between the respective obtained calibration points 521.
- the slopes (ai) are calculated according to the formula (1).
- the correction processing unit 212 calculates the deviation (Di) for each of the inclinations (ai) (S103).
- the deviation (Di) is calculated according to the formula (2).
- the correction processing unit 212 calculates the standard deviation (SD) of the slope calculated in step S102 (S104).
- the standard deviation (SD) is calculated according to the formula (4).
- the correction processing unit 212 determines whether or not the standard deviation is greater than a preset standard deviation threshold value (Sth) (SD>Sth) (S105). In step S105, it is determined whether or not there is an outlier for each calibration point 521. If the standard deviation is equal to or smaller than the threshold value of the standard deviation (S105 ⁇ No), the data processing device 200 ends the process.
- the correction processing unit 212 calculates the deviation (Di) for all between the calibration points 521 (S121). The correction processing unit 212 calculates the deviation (Di) according to formula (2).
- the correction processing unit 212 searches for the deviation (Dmax) having the largest value from among the deviations (Di) (S122). Then, it is determined whether the maximum deviation (Dmax) is greater than a deviation threshold value (Dth) (S123).
- step S141 If the deviation (Dmax) having the largest value is equal to or smaller than the deviation threshold value (Dth) (S123 ⁇ No), the correction processing unit 212 advances the process to step S141. If the deviation (Dmax) having the largest value is greater than the deviation threshold value (Dth) (S123 ⁇ Yes), the calibration point is corrected (S124: second step). At this time, the correction processing unit 212 performs the correction by replacing the calibration point 521 corresponding to the deviation (Dmax) having the largest value. The detailed procedure of the correction will be described later. In step S124, the outlier of the calibration point 521 is corrected, and the calibration point 521 is corrected.
- the correction processing unit 212 uses the corrected calibration points 521 to recalculate all the slopes (ai) between the obtained calibration points 521 (S131).
- the correction processing unit 212 re-calculates the standard deviation (SD) of the inclination calculated in step S131 (S132). Note that, prior to step S132, the calculation of the deviation (Di) is performed in the same manner as in step S103A.
- the correction processing unit 212 judges whether or not the standard deviation (SD) calculated in step S132 is greater than a threshold value (Sth) of the standard deviation (SD>Sth) (S133).
- the threshold value (Sth) of the standard deviation used in step S133 is the same as that used in step S105.
- the correction processing unit 212 adds 1 (h++) to the loop count (h) (S134) and returns the process to step S112.
- the measurement processing unit 213 performs voltage calibration using the corrected calibration point 521 (S141). After that, the correction processing unit 212 returns the process to step S102. However, after step S141 is performed, the process does not necessarily have to return to step S102.
- FIG. 9A is a diagram showing a table of specific values of the calibration points 521.
- the table shown in FIG. 9A has items for the number of the calibration point 521, the reference m/z, and the quadrupole RF voltage.
- the reference m/z is the m/z of the calibration ion used.
- the quadrupole RF voltage item has items for the correct value, the measured value, and the corrected value.
- the correct value is the quadrupole RF voltage in the case where the calibration ion is not confused with the impurity ion.
- the last row of the table in Figure 9A shows the standard deviation.
- the standard deviation is calculated for each of the correct value, the measured value, and the corrected value. Note that the standard deviation shown in Figure 9A is the standard deviation of the slope between the calibration points (ai shown in equation (1)).
- the correct value differs from the actual measured value because the calibration point 251 number "6" has been confused with an impurity ion.
- FIG. 9B is a diagram showing a graph of the deviation (Di) based on FIG. 9A.
- the horizontal axis represents the deviation number
- the vertical axis represents the magnitude of the deviation.
- the correct quadrupole RF voltage value (x1) is calculated based on the first deviation (D1)
- the correct quadrupole RF voltage value (x1) is calculated using the following formula (11).
- the correct quadrupole RF voltage value (xN) is calculated based on the Nth (last) deviation (DN)
- the correct quadrupole RF voltage value (xN) is calculated using the following formula (12).
- the correct quadrupole RF voltage value (xi) is calculated using the following equations (13) and (14).
- the correction processing unit 212 corrects the calibration point 512 based on the slope (ai).
- step numbers refer to the step numbers of the processing shown in FIG. 8.
- the correction processing unit 212 searches for the deviation (Di) having the largest value from among all the deviations.
- i is the deviation number shown in Fig. 9B.
- the correction processing unit 212 determines whether or not the deviation value found in step S122 is equal to or greater than the deviation threshold value (Dth). It is assumed that D5 (reference numeral 601) shown in FIG. 9B is equal to or greater than the deviation threshold value.
- the correction processing unit 212 corrects the calibration point 521 by replacing the calibration point 521 corresponding to the deviation (Dmax) having the largest value. Specifically, the quadrupole RF voltage (xi or xi+1) of the calibration point 521 (mi or mi+1) having the slope (ai) corresponding to the maximum deviation (Di) is replaced with a value corrected from the previous and next calibration points 521. The replacement is performed by the quadrupole RF voltage calculated according to equations (11) to (14).
- the correction processing unit 212 recalculates the tilt (ai) using the quadrupole RF voltage corrected in step S124.
- the correction processing unit 212 recalculates the standard deviation (SD) of the slope recalculated in step S131. In the case of the example shown in Fig. 9A, the standard deviation (SD) after correction is "0.099".
- the correction processing unit 212 compares the corrected standard deviation with the standard deviation threshold value (Sth).
- step S141 the measurement processing unit 213 performs the process of step S141.
- step S141 the measurement processing unit 213 performs voltage calibration for the corrected calibration point 521, using the corrected quadrupole RF voltage as the initial value (third step). Specifically, the measurement processing unit 213 performs the above-mentioned scan with the corrected quadrupole RF voltage as the initial value to reacquire the spectrum 501. Then, based on the reacquired spectrum 501, the measurement processing unit 213 estimates the quadrupole RF voltage corresponding to the peak value of the calibration ion shown in the spectrum 501a by correction. In this way, the correct quadrupole RF voltage can be acquired.
- the process performed in step S141 will be described later.
- step S133 if the corrected standard deviation is greater than the threshold value (S133 ⁇ Yes), the correction processing unit 212 repeats the processes from step S112 onwards until the standard deviation becomes equal to or less than the threshold value. However, if the number of loops is greater than a predetermined number (hmax), the correction processing unit 212 outputs an error (S151).
- equation (1) which shows the slope between calibration points
- the smaller the m/z difference (mi+1 - mi) the more sensitively it reflects the deviation in the amplitude of the quadrupole RF voltage.
- a calibration sample is selected so that the calibration points 521 have uniform m/z spacing.
- the m/z spacing is not uniform. In this way, when the m/z spacing of the calibration points 521 is not uniform, the m/z difference (mi+1 - mi) between the calibration points may be weighted when calculating the standard deviation.
- voltage calibration (S141 in FIG. 8) is performed, but this is not limited to the above.
- the correction result i.e., the quadrupole RF voltage after correction
- the voltage calibration may be performed as a process separate from the correction of the calibration point 521.
- (Voltage calibration) 10 is a diagram showing a voltage calibration method, with reference to FIG. Fig. 10 is an enlarged view of spectrum 501a and the vicinity of spectrum 502 in the quadrupole RF voltage-ion signal intensity relationship graph shown in Fig. 5.
- the process shown in Fig. 10 is performed by the measurement processing unit 213 shown in Fig. 1 in step S141 in Fig. 8.
- the measurement is started in a narrow scan range 511a centered on the quadrupole RF voltage (reference symbol 531) calculated by the corrected calibration point 521.
- the scan range 511 is gradually widened from the narrow scan range 511 like the scan range 511a centered on the quadrupole RF voltage (reference symbol 531) estimated by the correction of the calibration point 521.
- the number "6" in FIG. 9A shows that the difference between the quadrupole RF voltage by the correction of the calibration point 521 and the correct quadrupole RF voltage "111977” (V) is only 19 (V).
- the quadrupole RF voltage estimated by the correction of the calibration point 521 is not necessarily the true center of the spectrum 501a to be subjected to the voltage calibration. Therefore, in FIG. 10, the quadrupole RF voltage (reference symbol 531) calculated by the corrected calibration point 521 is shown at a position shifted from the center of the spectrum 501a.
- the scan range 551 is changed each time a measurement is performed.
- the scan range 511 is gradually expanded from scan range 511a to scan range 511c.
- the scan range 511 is then gradually expanded until the calibration ion spectrum 501a enters scan range 511c.
- the measurement processing unit 213 measures the calibration ions in the mass spectrometer 100 (see figure) in a scan range 511 centered on the quadrupole RF voltage calculated using the corrected calibration point 521. At this time, the measurement processing unit 213 changes the scan range 511 so that it becomes gradually wider each time a measurement is performed, as shown by scan ranges 511a to 511c in Figure 10.
- the first embodiment when performing voltage calibration, even if there is an impurity ion having an m/z adjacent to the calibration point 521, it is possible to correct the peak confusion between the impurity ion and the calibration ion. This makes it possible to determine an accurate quadrupole RF voltage derived from the calibration ion. As a result, when performing mass analysis of an unknown substance, highly accurate mass analysis is possible. Furthermore, according to the first embodiment, correction of the calibration point 521 and voltage calibration can be performed even if there is only one peak confusion.
- FIG. 11 is a flow chart showing the procedure of the mass spectrometry data processing method according to the second embodiment.
- Fig. 12 is a graph showing the relationship between m/z and quadrupole RF voltage.
- the step numbers are those shown in Fig. 11.
- the same step numbers are used for steps that are the same as those in Fig. 8, and descriptions thereof will be omitted as appropriate.
- the correction processing unit 212 calculates the slope and intercept 711 by applying the least squares method to each of the calibration points 521. Then, the correction processing unit 212 compares the standard deviation and the threshold value based on the calculated slope.
- the slope (ak) between the calibration points 521 is calculated one by one.
- the correction processing unit 212 obtains a straight line 701 close to all the calibration points 521 as shown in FIG. 12 by the least squares method according to the following formula (21). That is, the correction processing unit 212 calculates the slope and intercept 711 of the straight line 701 by applying the least squares method to each of the calibration points 521 (S102A and S131A). That is, the straight line 701 has a slope and intercept 711 calculated by applying the least squares method to each of the calibration points 521. In this way, in the second embodiment, the slope between the respective calibration points 521 is unified to a slope obtained by the least squares method.
- the correction processing unit 212 calculates the deviation (Di) (S103A).
- the deviation (Di) is defined as the length of the perpendicular line 721 from each calibration point 521 to the straight line 701 obtained by the least squares method.
- the length of the perpendicular line 721 is the distance from the straight line 701 having the slope and intercept 711 for each calibration point 521.
- the correction processing unit 212 calculates the distance (length of the perpendicular line 721) from the straight line 701 having the slope and intercept 711 for each calibration point 512 as the deviation (Di).
- the correction processing unit 212 calculates the standard deviation (SD) of the deviation (Di) calculated in step S103A (S104A).
- the correction processing unit 212 corrects the calibration point 521 based on the slope.
- SD standard deviation
- Sth preset threshold
- xi_cor indicates the value of the corrected quadrupole RF voltage
- b is the intercept 711 of the quadrupole RF voltage axis calculated by the least squares method.
- the corrected calibration point 521 is corrected to lie on the straight line 701 calculated by the least squares method.
- step S121A in FIG. 11 the deviation (Di) is calculated in the same manner as in step S103A.
- the slope and intercept 711 are calculated in the same manner as in step S102A.
- step S131A is the same process as step S102A
- step S132A is the same process as step S104A.
- step S131A and step S132A the same process as in step S103A is performed.
- the misidentification of the peaks of impurity ions and calibration ions is corrected based on the slope calculated by using the least squares method for the calibration point 521.
- the slope between the calibration points 521 is calculated at once by the least squares method, compared to the first embodiment in which the slope between the calibration points 521 is calculated one by one. Therefore, the method shown in the second embodiment has simpler calculations and can be performed quickly compared to the method shown in the first embodiment.
- FIG. 13 is a diagram showing an example of the configuration of a mass spectrometry system 1a according to the third embodiment.
- a mass spectrometer 100a in the mass spectrometer system 1a shown in FIG. 13 is similar to the mass spectrometer system 1 shown in FIG. 1 in the following respects.
- the analysis section 103 is provided with a time-of-flight mass spectrometer 400 .
- a power supply 410 that applies a pulse voltage is connected to a pushing electrode 401 that constitutes a time-of-flight mass spectrometer 400.
- This power supply 410 is controlled by a voltage control device 300.
- FIG. 13 the same components as those in FIG. 1 are denoted by the same reference numerals.
- ions generated in the ion source 151 are introduced into the first differential pumping section 101 through the aperture 121, as in FIG. 1.
- the first differential pumping section 101 is evacuated by the pump 111, as in FIG. 1. This allows the first differential pumping section 101 to be maintained at a vacuum level of 10 Pa to 500 Pa.
- the ions that have passed through the first differential pumping section 101 are introduced into the second differential pumping section 102 through the aperture 122.
- the second differential pumping section 102 is evacuated by the pump 112, as in FIG. 1. This allows the second differential pumping section 102 to be maintained at a vacuum level of 0.1 Pa to 10 Pa. Also, as in FIG.
- the second differential pumping section 102 is provided with an ion guide 130 that focuses the ions.
- the ions focused by the ion guide 130 pass through the aperture 123 and are introduced into the analysis section 103 in which the time-of-flight mass spectrometer 400 is provided.
- the analysis section 103 is evacuated by a pump 113. This keeps the analysis section 103 at a pressure of 1E-5 Pa or less.
- the time-of-flight mass spectrometer 400 is composed of a pusher electrode 401, an acceleration electrode 402, a reflector electrode 403, and a detector 152.
- a pulse voltage is applied to the pusher electrode 401 by the power supply 410.
- the pusher electrode 401 bends the traveling direction of some of the ions incident from the ion guide 130 to a direction perpendicular to the incident direction of the ions.
- the bent ions head toward the acceleration electrode 402.
- the ions whose traveling direction is bent are accelerated by the acceleration electrode 402 and reflected by the reflector electrode 403 (white arrow in Figure 13).
- the ions reflected by the reflector electrode 403 are detected by the detector 152.
- a multi-channel plate type detector 152 is generally used as the detector 152.
- the detection intensity of the ions is converted into an electrical signal (output signal) by the detector 152 and sent to the data processing device 200.
- the output signal output from the detector 152 is converted into digital data at a fixed sampling period and then sent to the data processing device 200.
- the conversion to digital data is performed by an analog-to-digital converter (ADC) (not shown) or a pulse counting unit (not shown).
- ADC analog-to-digital converter
- the data processing device 200 stores the sent digital data in the memory unit 220.
- the amplitude of the pulse voltage applied to the pusher electrode 401 is used instead of the quadrupole RF voltage of the first embodiment.
- the pulse voltage is a voltage obtained as a result of measuring a plurality of calibration ions with the mass spectrometry device 100.
- the quadrupole RF voltage axis of the m/z-quadrupole RF voltage relationship graphs shown in FIGS. 6 to 7B becomes the pulse voltage axis.
- the peak of the impurity ion may be confused with the peak of the calibration ion.
- a pulse voltage instead of the quadrupole RF voltage of the first embodiment, it is possible to perform processing similar to that of the first and second embodiments. This makes it possible to correct the confusion of the peaks. This makes it possible to obtain an accurate m/z when performing mass analysis on an unknown substance.
- the present invention is not limited to the above-described embodiments, and includes various modified examples.
- the above-described embodiments have been described in detail to clearly explain the present invention, and are not necessarily limited to having all of the configurations described. It is also possible to replace part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment. It is also possible to add, delete, or replace part of the configuration of each embodiment with other configurations.
- the above-mentioned configurations, functions, plot processing unit 211, correction processing unit 212, measurement processing unit 213, storage unit 220, etc. may be realized in hardware by designing some or all of them as an integrated circuit, for example.
- the above-mentioned configurations, functions, etc. may be realized in software by a processor such as a CPU interpreting and executing a program that realizes each function. Information such as the program, table, file, etc.
- a recording device such as an SSD (Solid State Drive), or a recording medium such as an IC (Integrated Circuit) card, an SD (Secure Digital) card, or a DVD (Digital Versatile Disc), in addition to being stored in an HD (Hard Disk).
- IC Integrated Circuit
- SD Secure Digital
- DVD Digital Versatile Disc
- control lines and information lines shown are those that are considered necessary for the explanation, and not all control lines and information lines in the product are necessarily shown. In reality, it can be considered that almost all components are interconnected.
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Electron Tubes For Measurement (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
Description
本発明は、データ処理装置、質量分析システム及び質量分析データ処理方法の技術に関する。 The present invention relates to a data processing device, a mass spectrometry system, and a mass spectrometry data processing method.
質量分析装置では、m/zが既知のイオンを測定し、特定のパラメータ(四重極質量フィルタでは電極に印可する電圧)とm/zの関係を定めることが行われる。m/zが未知のサンプルの測定では、このような技術を使用することによってm/zが決定される。 In a mass spectrometer, ions with a known m/z are measured and the relationship between a specific parameter (the voltage applied to the electrodes in a quadrupole mass filter) and m/z is determined. When measuring a sample with an unknown m/z, the m/z is determined by using such techniques.
特許文献1には、「分析試料および質量校正試料を分析試料イオンおよび質量校正試料イオンとしてイオン化させるようにイオン化制御部1によりイオン化部21が制御される。分析試料イオンおよび質量校正試料イオンを捕捉させ、分析試料イオンをプリカーサイオンとして選択せずに質量校正試料イオンを選択的に複数の校正用フラグメントイオンに解離させるようにイオントラップ23がイオン解離制御部2により制御される。分析試料イオンおよび複数の校正用フラグメントイオンの質量分析が行われるように質量分析器24が分析制御部3により制御される。質量分析により得られる質量分析スペクトルにおいて複数の校正用フラグメントイオンの質量電荷比に基づいて分析試料イオンの質量電荷比が校正部4により校正される」質量分析制御装置、質量分析装置、質量分析制御方法および質量分析方法が開示されている(要約参照)。
特許文献2には、イオン化の極性などが異なる複数の条件で質量スペクトルを測定し、参照データとの差異がある場合に質量軸校正を行う方法が開示されている。
本発明で解決しようとする課題は、質量軸校正を行う際に、校正点に対しm/zが隣接する夾雑物イオンが存在している場合でも正確な質量軸を決定することである。 The problem that this invention aims to solve is to determine an accurate mass axis when performing mass axis calibration, even when there are impurity ions whose m/z values are adjacent to the calibration point.
特許文献1の記載の技術では、試料イオンや夾雑物イオンが共存しているため、校正用フラグメントイオンと夾雑物イオンのm/zが隣接している場合、ピークの取り違えにより正確な質量軸が決定できない場合がある。
In the technology described in
特許文献2に記載の技術では、校正用試料イオンと夾雑物イオンのm/zが隣接している場合、隣接するスペクトルのピークのいずれが校正用試料イオンであるかを判断することができない。そのため、特許文献2に記載の技術では、参照データとの差異から異常を検知できても、校正用試料イオンであるという確証を得ることができない。
In the technology described in
このような背景に鑑みて本発明がなされたのであり、本発明は、精度の高い質量分析を可能とすることを課題とする。 The present invention was made in light of this background, and its objective is to enable highly accurate mass analysis.
前記した課題を解決するため、本発明は、m/zが既知である、複数の校正用イオンを質量分析装置で測定した結果として取得される電圧と、それぞれの校正用イオンの前記m/zと、を、前記電圧及び前記m/zのそれぞれを座標軸として有する座標に校正点としてプロットするプロット処理部と、それぞれの前記校正点に関して外れ値があれば、前記外れ値を補正することで、前記校正点の補正を行う補正処理部と、を有する。
その他の解決手段は実施形態中において適宜記載する。
In order to solve the above-mentioned problems, the present invention comprises a plot processing unit that plots a voltage obtained as a result of measuring a plurality of calibration ions, whose m/z is known, using a mass spectrometer, and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes, and a correction processing unit that corrects the calibration points by correcting the outliers if there are any outliers for each of the calibration points.
Other solutions will be described in the embodiments as appropriate.
本発明によれば、精度の高い質量分析が可能となる。 The present invention makes it possible to perform mass analysis with high accuracy.
次に、本発明を実施するための形態(「実施形態」という)について、適宜図面を参照しながら詳細に説明する。 Next, a form for implementing the present invention (referred to as an "embodiment") will be described in detail with reference to the drawings as appropriate.
<第1実施形態>
[質量分析システム1]
図1は、第1実施形態に係る質量分析システム1の構成例を示す図である。
質量分析システム1は、質量分析装置100、電圧制御装置300、DC電源301,303、RF電源302を有している。さらに、質量分析システム1は、出力装置201が接続されているデータ処理装置200を有している。ちなみに、図1において、電圧制御装置300とDC電源301,303、RF電源302、データ処理装置200とを結ぶ破線、及び、検出器152とデータ処理装置200を結ぶ破線は信号線を示している。また、DC電源303と、イオンガイド130とを結ぶ実線、RF電源302と四重極質量フィルタ140及びDC電源301とを結ぶ実線は、電圧が印加される電線を示している。
First Embodiment
[Mass spectrometry system 1]
FIG. 1 is a diagram showing an example of the configuration of a
The
質量分析装置100は、イオン源151、第1差動排気部101、第2差動排気部102、分析部103を有している。
イオン源151で生成されたイオンは細孔121を介して第1差動排気部101に導入される。イオン源151は、エレクトロスプレーイオン源、大気圧化学イオン源、大気圧光イオン源、大気圧マトリックス支援レーザー脱離イオン源等によって構成されており、大気圧、もしくは、大気圧より低真空で動作する。
The
Ions generated in the
第1差動排気部101は、ポンプ111によって排気されている。これにより、第1差動排気部101は、10Paから500Paの真空度に維持されている。第1差動排気部101を通過したイオンは、細孔122を介して第2差動排気部102に導入される。第2差動排気部102は、ポンプ112によって排気されている。これにより、第2差動排気部102は、0.1Pa~10Paの真空度に維持されている。また、第2差動排気部102には、イオンを収束させるイオンガイド130が設置されている。イオンガイド130で収束されたイオンは細孔123を介して四重極質量フィルタ140が設置されている分析部103へ導入される。
The first
なお、イオンガイド130には、DC電源303から静電圧(以降、DC電圧と称する)が印加される。また、RF電源302は、RF電源302が生成する高周波電圧(以降、RF電圧と称する)と、DC電源301から送られたDC電圧の合成波を四重極質量フィルタ140に印加する。
The
また、図1に示すように、イオンガイド130と、四重極質量フィルタ140とはコンデンサ等の誘電体153を介して接続している。前記したように、RF電圧は、電圧制御装置300で制御されるRF電源302から四重極質量フィルタ140の四重極ロッド電極141(図3A及び図3B参照)に供給される。そして、RF電圧の高周波成分は、四重極ロッド電極141から誘電体153を通してイオンガイド130のイオンガイドロッド電極に供給される。このような構成を有することにより、イオンガイド130と四重極質量フィルタ140とに対し、個別にRF電圧を供給する構成に比べて電源の数を少なくすることができ、質量分析システム1を安価にすることができる。
Also, as shown in FIG. 1, the
分析部103は、ポンプ113によって、排気されている。これにより、分析部103は、1E-3Pa以下の圧力に維持されている。四重極質量フィルタ140では、イオンがm/zに応じて分離される。
The
四重極質量フィルタ140を通過したイオンは検出器152で検出される。検出器152として、電子増倍管や、シンチレータと光電子増倍管とを組み合わせたタイプや、マルチチャンネルプレートが一般に使用可能である。検出器152により検出されたイオンの検出強度は、電気的な信号(出力信号)に変換され、データ処理装置200へ送られる。この際、検出器152から出力された出力信号は、一定のサンプリング周期(典型的には1μs~1000μs)のデジタルデータに変換された後、データ処理装置200へ送られる。デジタルデータへの変換は、図示しないアナログ-デジタルコンバータ(ADC)や、図示しないパルスカウンティングユニットで行われる。データ処理装置200は、送られたデジタルデータを記憶部220に蓄積する。
Ions that pass through the
データ処理装置200では、記憶部220の外に、プロット処理部211、補正処理部212、測定処理部213が実行されている。
プロット処理部211は、図6~図7Bに示すm/z-四重極RF電圧関係グラフに校正点521(図6~図7B)をプロットする。
補正処理部212は、プロット処理部211によってプロットされた校正点521を基に校正点521の補正を行う。
測定処理部213は、質量分析装置100による測定を制御する。
In the
The
The
The
記憶部220はメモリやHD等からなり、スペクトル501(図5参照)のデータの他、補正を行なうために必要な数値や関係式等の情報を記憶して保持することができる。また、プロット処理部211、補正処理部212、測定処理部213は、演算機能や演算に必要な数値を一時保存するための一時メモリを有する。さらに、データ処理装置200は、これらの情報を蓄積、変換する以外にも、各電極等を制御する電圧制御装置300等をコントロールする機能や、出力装置201に情報を出力する機能を有する。出力装置201は、ディスプレイや、プリンタ等で構成されている。出力装置201は、スペクトル501そのものや、スペクトル501のピークに対応するm/zや、信号強度、測定対象の物質の有無等の情報を出力する。
The
なお、図1に示す例では、記憶部220がデータ処理装置200に搭載されている構成であるが、これに限らない。例えば、記憶部220は、データベース等、データ処理装置200とは別の装置として設置されてもよい。
In the example shown in FIG. 1, the
[ハードウェア構成]
図2は、データ処理装置200のハードウェア構成例を示す図である。
データ処理装置200は、RAM等で構成されるメモリ231、CPUや、GPU等で構成される演算装置232、HDや、SSD等で構成される記憶装置233を備える。また、データ処理装置200は、キーボードや、マウス等で構成される入力装置234、出力装置201、通信装置235を備える。
[Hardware configuration]
FIG. 2 is a diagram illustrating an example of the hardware configuration of the
The
データ処理装置200に記憶部220(図1参照)が搭載されている場合、記憶装置233は記憶部220に該当する。出力装置201は、図1に示す出力装置201である。
If the
また、記憶装置233に格納されているプログラムがメモリ231にロードされる。そして、ロードされたプログラムが演算装置232によって実行される。これによって、図1に示すプロット処理部211、補正処理部212、測定処理部213が具現化する。
In addition, a program stored in the
[四重極質量フィルタ140]
図3A及び図3Bは四重極質量フィルタ140の構成を示す図である。図3Aには、四重極質量フィルタ140の斜視図が示されており、図3Bには、四重極質量フィルタ140の断面と、四重極質量フィルタ140に対する電圧制御の示す図が示されている。
図3A及び図3Bに示すように、四重極質量フィルタ140は4本の四重極ロッド電極141(141a~141d)により構成される。そして、図3Bに示すように、それぞれの四重極ロッド電極141には、RF電圧とDC電圧が印加される。RF電圧は、電圧制御装置300で制御されるRF電源302で生成される交流電圧である。また、DC電圧は、電圧制御装置300で制御されるDC電源301で生成される直流電圧である。実際には、RF電圧にDC電圧が印加された合成波が四重極ロッド電極141のそれぞれに印加される。
[Quadrupole mass filter 140]
3A and 3B are diagrams showing the configuration of the quadrupole
As shown in Figures 3A and 3B, the quadrupole
そして、RF電圧及びDC電圧の合成波は、隣接する四重極ロッド電極141の間で逆相となり、対向する四重極ロッド電極141の間で同相となるよう印加される。つまり、図3Bに示す例では、四重極ロッド電極141a,141cのペアと、四重極ロッド電極141b,141dのペアとの間に逆相のRF電圧が印加される。
The composite wave of the RF voltage and DC voltage is applied so that it is in phase with each other between adjacent
また、前記したように、DC電圧は、電圧制御装置300で制御されるDC電源301で生成される直流電圧である。ここで、四重極ロッド電極141a,141cに印加されるDC電圧をVDC1とすると、四重極ロッド電極141b,141dに印加されるDC電圧は-VDC1である。それぞれの四重極ロッド電極141に印可されるRF電圧とDC電圧を、四重極RF電圧、四重極DC電圧と適宜称する。四重極RF電圧の典型的な電圧振幅は数100V~数kVであり、周波数は500kHz~2MHz程度である。四重極DC電圧は、数10V~数100V程度である。
As described above, the DC voltage is a direct current voltage generated by the
(四重極質量フィルタ140の動作)
次に、図4を参照して四重極質量フィルタ140の動作について説明する。
図4は、四重極質量フィルタ140の動作と、四重極RF電圧及び四重極DC電圧との関係を示す図である。
四重極質量フィルタ140内で安定な軌道運動が可能なイオンのm/z範囲は、四重極RF電圧の振幅と、四重極DC電圧の値に依存している。図4に示す安定領域R1~R3の内側に存在するイオンのみが、四重極質量フィルタ140を透過することができる。ここで、安定領域R1は線R1aの線内の領域であり、安定領域R2は線R2aの線内の領域であり、安定領域R3は線R3aの線内の領域である。イオンのm/zごとに安定領域R1~R3は異なる。そして、安定領域R1~R3は、m/zの小さいイオンから大きいイオンまで図4に示した関係で並んでいる。つまり、安定領域R1は、あるm/zを有するイオンの安定領域である。同様に、安定領域R2は、安定領域R1のイオンとは異なるm/zを有するイオンの安定領域であり、安定領域R3は、安定領域R1,R2のイオンとは異なるm/zを有するイオンの安定領域である。
(Operation of the quadrupole mass filter 140)
Next, the operation of the quadrupole
FIG. 4 is a diagram illustrating the operation of the quadrupole
The m/z range of ions capable of stable orbital motion in the quadrupole
あるm/zの安定領域R1~R3の頂点近傍に四重極RF電圧と四重極DC電圧とを設定すれば、そのm/zを有するイオンのみを透過させることが可能である。また、図4中に示すスキャンラインL1のように各m/zのイオンの安定領域R1~R3の頂点近傍を通るように、四重極RF電圧をスキャンすることが行われる。この際、四重極RF電圧と四重極DC電圧の関係を維持しつつ、四重極RF電圧をスキャンすることにより、図5に示すスペクトル501を得ることができる。つまり、各m/zを有するイオンの検出が可能となる。
If the quadrupole RF voltage and quadrupole DC voltage are set near the apex of the stability region R1-R3 of a certain m/z, it is possible to transmit only ions having that m/z. In addition, the quadrupole RF voltage is scanned so as to pass near the apex of the stability region R1-R3 of each m/z ion, as shown by the scan line L1 in Figure 4. In this case, by scanning the quadrupole RF voltage while maintaining the relationship between the quadrupole RF voltage and the quadrupole DC voltage, the
(校正点補正処理)
次に、図5~図7Bを参照して、校正点補正処理を示す。
図5は、四重極RF電圧と、イオン信号強度の関係を示す四重極RF電圧-イオン信号強度関係グラフを示す図である。また、図6は、m/zと、四重極RF電圧との関係を示すm/z-四重極RF電圧関係グラフを示す図である。
図5では、横軸に四重極RF電圧が示され、縦軸にイオン信号強度が示されている。そして、図5には、校正用イオンによるスペクトル501が示されている。校正用イオンとは、校正用試料から生成されるイオンであり、そのm/zは既知である。また、図5には、夾雑物イオンによるスペクトル502が検出されている。夾雑物イオンは、コンタミネーション等による夾雑物に由来するイオンである。
(Calibration point correction process)
Next, the calibration point correction process will be described with reference to FIGS. 5 to 7B.
Fig. 5 is a graph showing the relationship between the quadrupole RF voltage and the ion signal intensity, and Fig. 6 is a graph showing the relationship between the m/z and the quadrupole RF voltage.
In Fig. 5, the horizontal axis indicates the quadrupole RF voltage, and the vertical axis indicates the ion signal intensity. Fig. 5 also shows a
なお、以降の説明において、四重極RF電圧とは、四重極RF電圧の振幅値を示すものとする。また、四重極RF電圧は、複数の校正用イオンを質量分析装置100で測定した結果として取得される電圧である。
In the following description, the quadrupole RF voltage refers to the amplitude value of the quadrupole RF voltage. The quadrupole RF voltage is the voltage obtained as a result of measuring multiple calibration ions with the
図5に示すように、夾雑物イオンによるスペクトル502が、校正用イオンによるスペクトル501aに対し、近接した状態で検出されている。
As shown in FIG. 5,
また、図5において、電圧範囲であるスキャン範囲511が示されている。スキャン範囲511は、イオン検出のために行われる四重極RF電圧のスキャン幅である。スキャン範囲511については後記する。
Also shown in FIG. 5 is a
また、図6に示すm/z-四重極RF電圧関係グラフでは、校正点521がプロットされている。図6に示すようにm/z-四重極RF電圧関係グラフは、座標軸として、m/z軸と、四重極RF電圧軸を有している座標である。つまり、m/z-四重極RF電圧関係グラフは、電圧及びm/zのそれぞれを座標軸として有する座標である。
Also, in the m/z-quadrupole RF voltage relationship graph shown in FIG. 6,
校正点521は、図5に示すスペクトル501におけるイオン信号強度のピーク値に対応する四重極RF電圧と、m/zとの関係を示している。前記したように、校正用イオンのm/zは既知であるため、図6に示すようなm/z-四重極RF電圧関係グラフの描画が可能である。校正点521は、複数の校正用イオンを質量分析装置100で測定した結果として取得される電圧と、それぞれの校正用イオンのm/zと、が座標にプロットされたものである。
図5では、前記したように夾雑物イオンによるスペクトル502が、校正用イオンによるスペクトル501aの近傍で検出されている。また、スペクトル501aのピーク値が、スペクトル502のピーク値より低い。また、図5では、スペクトル501aと、スペクトル502とが同一のスキャン範囲511に収まっている。これは、スペクトル501aと、スペクトル502とが1つのスペクトル501として検出されていることを示している。これにより、データ処理装置200が、夾雑物イオンを校正用イオンと取り違えて判定するおそれがある。
In FIG. 5, as described above,
そのため、図6に示すように、スペクトル501aに由来する校正点521aの近傍に、夾雑物イオンに由来する校正点521bが近接してプロットされる。
As a result, as shown in FIG. 6,
このように、校正用イオンに対して、隣接するm/zを有する夾雑物イオンが存在する場合がある。このような場合、図5のスペクトル501a,502に示すように、校正用イオンに由来するスペクトル501aのピークと、夾雑物イオンに由来するスペクトル502のピークが近接する。このように、データ処理装置200が、校正用イオンと夾雑物イオンにおけるピークを取り違える可能性がある。このような、校正用イオンと夾雑物イオンのピークを取り違えが生じると、図6に示すように四重極RF電圧の値が、実際の校正用イオンの四重極RF電圧と異なる値となる。すなわち、本来、スペクトル501aを生じる校正用イオンのm/zが結び付けられるべき四重極RF電圧に、夾雑物イオンのm/zが結び付けられてしまう。m/zが未知の物質を測定する際に、図6に示すようなm/z-四重極RF電圧関係グラフと、得られる四重極RF電圧を基に、測定された物質のm/zが出力される。しかし、スペクトル501aを生じる校正用イオンのm/zが結び付けられるべき四重極RF電圧に、夾雑物イオンのm/zが結び付けられてしまっていると、誤ったm/zが出力されてしまう。
In this way, there may be impurity ions having adjacent m/z to the calibration ions. In such a case, as shown in
なお、校正用イオンによるスペクトル501の範囲が所定の(一定の)値になるように四重極DC電圧は決定されている。
The quadrupole DC voltage is determined so that the range of the
また、スペクトル501を取得する際に、四重極質量フィルタ140をイオンが抜けてから検出器152で検出されるまでのタイムラグがある。このため、検出器152でイオンが検出されたタイミングにおける四重極RF電圧の値と、イオンが四重極質量フィルタ140を通過するタイミングでの四重極RF電圧の値が異なる現象が生じる。そのため、スキャン速度により検出されるイオンのm/zが変化する。
In addition, when acquiring
質量分析時において、四重極質量フィルタ140は、一定時間、四重極RF電圧を固定した状態で、特定のm/zを有するイオンを選択的に透過させる動作を行う。特定のm/zを有するイオンを透過させる四重極RF電圧を決定するため、測定処理部213(図1参照)は、所定の範囲で四重極RF電圧を再スキャンし、図5に示すスペクトル501のピーク値を取得する。このピーク値を基に、校正用イオンのm/zに対応する四重極RF電圧を決定することを電圧校正と称する。なお、1回目のスキャンは、後記するように夾雑物イオンが存在する状態でのスキャンであり、再スキャンは、後記する校正点521(図6参照)の補正が行われた後のスキャンである。
During mass analysis, the quadrupole
高精度で電圧校正を行うためには低いスキャン速度で質量分析が行われることが望ましい。低いスキャン速度とは、例えば、0.1Da/s~100Da/s(1.66×10-25g/s~1.66×10-22g/s)程度である。一方、低いスキャン速度では、所定の範囲のスペクトル501を取得する時間が長くなる。つまり、測定処理部213が、すべての四重極RF電圧の範囲を低いスキャン速度でスキャンすると、膨大な時間が要されることになる。そのため、電圧校正に要する時間も長くなるというデメリットがある。
In order to perform voltage calibration with high accuracy, it is desirable to perform mass analysis at a low scan rate. A low scan rate is, for example, about 0.1 Da/s to 100 Da/s (1.66×10 −25 g/s to 1.66×10 −22 g/s). On the other hand, at a low scan rate, it takes a long time to acquire a
そこで、図5の一点鎖線で示すスキャン範囲511のように校正用イオンの近傍のみを低いスキャン速度でスキャンしてスペクトル501を取得することが行われる。このようにして、電圧校正の精度と電圧校正にかかる時間の両立を測ることが多い。
Therefore, a
続いて、図7A及び図7Bを参照して、本実施形態における校正点補正処理を説明する。
図7A及び図7Bは、m/z-四重極RF電圧関係グラフを示す図である。図7Aは、夾雑物イオンが検出されていない例を示し、図7Bは、夾雑物イオンによるピークの取り違えが生じている例を示す。
図7Aに示すように、四重極質量フィルタ140を安定的に透過する校正用イオンのm/zは四重極RF電圧(の振幅)に比例する。したがって、夾雑物イオンが検出されていない場合、図7Aに示すように、式(1)で示される校正点521の間(mi,xi)を結ぶ直線の傾き(ai)はすべての校正点521の間でほぼ同じ値となる。ここで、傾き(ai)は校正点521それぞれの間の傾きである。
Next, the calibration point correction process in this embodiment will be described with reference to FIGS. 7A and 7B.
7A and 7B are graphs showing the relationship between m/z and quadrupole RF voltage, in which Fig. 7A shows an example in which no impurity ions are detected, and Fig. 7B shows an example in which peaks are confused due to impurity ions.
As shown in Fig. 7A, the m/z of the calibration ions stably passing through the quadrupole
ai = (xi+1 - xi)/(mi+1 - mi) ・・・(1) ai = (xi+1 - xi)/(mi+1 - mi)...(1)
ここで、mi(i=1,・・・,N)は、m/zの値であり、xi(i1,・・・,N)は、四重極RF電圧の値である。そして、iは校正点521の番号である。
Here, mi (i = 1, ..., N) is the m/z value, xi (i1, ..., N) is the value of the quadrupole RF voltage, and i is the number of the
これに対して、図7Bでは、(mk,xk)において、校正用イオン(校正点521a)と、隣接する夾雑物イオンの取り違えが生じている(校正点521b)。図6及び図7Bに示すように、夾雑物イオンの校正点521bが外れ値としてプロットされる。ちなみに、(mk,xk)とは、m/z-四重極RF電圧関係グラフにおける校正点521の座標を示す。校正用イオンと、隣接する夾雑物イオンの取り違えにより、正しい校正点521aが、夾雑物イオンの校正点521bとして認識され、対応する四重極RF電圧がずれる。そのため、mkの両側において傾き(ak-1,ak)の値が変化し、その結果、aiの偏差(Di)が大きくなる。偏差(Di)は以下の式(2)で定義される。このように、偏差(Di)は、それぞれの傾き(ai)に関する偏差である。
In contrast, in FIG. 7B, at (mk, xk), the calibration ion (
Di = (ai - a_mean)2 ・・・(2) Di = (ai - a_mean) 2 ...(2)
ここで、a_meanは、以下の式(3)で定義されるaiの平均である。ただし、偏差は、aiの平均に対する偏差でなくてもよい。例えば、偏差(Di)は、Di=ai - ai-1等で定義されてもよい。 Here, a_mean is the average of ai defined by the following formula (3). However, the deviation does not have to be the deviation from the average of ai. For example, the deviation (Di) may be defined as Di = ai - ai-1, etc.
図7Aに示すように、夾雑物イオンと、校正用イオンとの取り違えが生じていない場合、隣接する校正点間の関係、つまり隣接する校正点521の間の傾き(ai)(i=1,…,N-1)の値は、ほぼ同じ値となる。これを利用して、本実施形態では、夾雑物イオンとピークを取り違えた校正点521の補正が行われる。
As shown in FIG. 7A, if there is no confusion between impurity ions and calibration ions, the relationship between adjacent calibration points, that is, the value of the slope (ai) (i = 1, ..., N-1) between adjacent calibration points 521, will be approximately the same. Using this, in this embodiment, correction is performed for
具体的には、本実施形態において、補正処理部212(図1参照)は、以下の式(4)で定義される標準偏差(SD)がしきい値(Sth)以上となるか否かでピークの取り違えが生じているか否かを判断する。換言すれば、式(4)で定義される標準偏差(SD)がしきい値(Sth)以上となるか否かで、それぞれの校正点521に関して外れ値があるか否かが判定される。
Specifically, in this embodiment, the correction processing unit 212 (see FIG. 1) determines whether a peak has been confused based on whether the standard deviation (SD) defined by the following formula (4) is equal to or greater than a threshold value (Sth). In other words, whether there is an outlier for each
(フローチャート)
図8は、第1実施形態に係る質量分析データ処理方法の手順を示すフローチャートである。適宜、図1、図6等を参照する。
まず、プロット処理部211は、質量分析装置100による校正用イオンの測定結果をm/z-四重極RF電圧関係グラフにプロットするプロット処理を行う(S101:第1のステップ)。この際、プロット処理部211は、図5に示すスキャン範囲511によるスキャンの結果、得られる四重極RF電圧と、予めわかっている校正用イオンのm/zとを基にプロットを行う。
(flowchart)
8 is a flowchart showing the procedure of the mass spectrometry data processing method according to the first embodiment, with reference to FIGS.
First, the
その後、補正処理部212は、得られている校正点521それぞれの間の傾き(ai)をすべて算出する(S102)。傾き(ai)は式(1)に従って算出される。
そして、補正処理部212は、それぞれの傾き(ai)に関する偏差(Di)を算出する(S103)。偏差(Di)は、式(2)に従って算出される。
続いて、補正処理部212は、ステップS102で算出した傾きの標準偏差(SD)を算出する(S104)。標準偏差(SD)は、式(4)に従って算出される。
そして、補正処理部212は、標準偏差が、予め設定されている標準偏差のしきい値(Sth)より大きい(SD>Sth)か否かを判定する(S105)。ステップS105によって、それぞれの校正点521に関して外れ値があるか否かが判定される。
標準偏差が、標準偏差のしきい値以下の場合(S105→No)、データ処理装置200は処理を終了する。
Thereafter, the
Then, the
Next, the
Then, the
If the standard deviation is equal to or smaller than the threshold value of the standard deviation (S105→No), the
標準偏差が標準偏差のしきい値より大きい場合(S105→Yes)、補正処理部212はループ回数(h)を「0」に設定する(h=0)(S111)。ステップS105で「Yes」が判定されることは、補正処理部212が、それぞれの校正点521に関して外れ値がある、と判定したことになる。つまり、補正処理部212は、偏差(Di)の標準偏差(SD)が、予め設定されている閾値より大きい場合、外れ値があると判定する。
続いて、補正処理部212は、ループ回数(h)が予め設定されている最大ループ回数(hmax)より大きい(h>hmax)か否かを判定する(S112)。
ループ回数(h)が最大ループ回数(hmax)より大きい場合(S112→Yes)、補正処理部212は、エラー出力を行う(S151)。
If the standard deviation is greater than the standard deviation threshold value (S105→Yes), the
Next, the
If the loop count (h) is greater than the maximum loop count (hmax) (S112→Yes), the
ループ回数(h)が最大ループ回数(hmax)以下の場合(S112→No)、補正処理部212は、校正点521の間のすべてについて偏差(Di)を算出する(S121)。補正処理部212は、式(2)に従って偏差(Di)を算出する。
If the number of loops (h) is equal to or less than the maximum number of loops (hmax) (S112→No), the
そして、補正処理部212は、偏差(Di)の中から、最も大きな値を有する偏差(Dmax)を探索する(S122)。
そして、最も大きな値を有する偏差(Dmax)が偏差のしきい値(Dth)より大きいか否かを判定する(S123)。
Then, the
Then, it is determined whether the maximum deviation (Dmax) is greater than a deviation threshold value (Dth) (S123).
最も大きな値を有する偏差(Dmax)が偏差のしきい値(Dth)以下である場合(S123→No)、補正処理部212はステップS141へ処理を進める。
最も大きな値を有する偏差(Dmax)が偏差のしきい値(Dth)より大きい場合(S123→Yes)、校正点の補正を行う(S124:第2のステップ)。この際、補正処理部212は、最も大きな値を有する偏差(Dmax)に該当する校正点521を置換することによって補正する。補正の詳細な手順については後記する。ステップS124によって、校正点521の外れ値が補正され、校正点521の補正が行われる。
If the deviation (Dmax) having the largest value is equal to or smaller than the deviation threshold value (Dth) (S123→No), the
If the deviation (Dmax) having the largest value is greater than the deviation threshold value (Dth) (S123→Yes), the calibration point is corrected (S124: second step). At this time, the
そして、補正処理部212は、補正後の校正点521を用いて、得られている校正点521の間の傾き(ai)を、再度、すべて算出する(S131)。
続いて、補正処理部212は、ステップS131で算出した傾きの標準偏差(SD)を、再度、算出する(S132)。なお、ステップS132の前段において、ステップS103Aと同様の偏差(Di)の算出が行われる。
そして、補正処理部212は、ステップS132で算出した標準偏差(SD)が標準偏差のしきい値(Sth)より大きい(SD>Sth)か否かを判定する(S133)。ステップS133で使用される標準偏差のしきい値(Sth)は、ステップS105で使用されたものと同じものである。
Then, the
Next, the
Then, the
標準偏差(SD)が標準偏差のしきい値(Sth)より大きい場合(S133→Yes)、補正処理部212は、ループ回数(h)を1加算(h++)し(S134)、ステップS112へ処理を戻す。
If the standard deviation (SD) is greater than the standard deviation threshold value (Sth) (S133 → Yes), the
標準偏差(SD)が標準偏差のしきい値(Sth)以下の場合(S133→Yes)、測定処理部213は、補正した校正点521を用いて電圧校正を行う(S141)。その後、補正処理部212は、ステップS102へ処理を戻す。ただし、ステップS141が行われた後、必ずしもステップS102へ処理を戻さなくてもよい。
If the standard deviation (SD) is equal to or less than the standard deviation threshold value (Sth) (S133→Yes), the
(具体例)
続いて、図6、図8~図9Bを参照しつつ、図8で示す質量分析データ処理方法の具体例を示す。
図9Aは、校正点521の具体的な値の表を示す図である。
図9Aに示す表では、校正点521の番号、参照用のm/z、四重極RF電圧の項目を有する。参照用としてのm/zは、使用される校正用イオンのm/zである。また、四重極RF電圧の項目は、正しい値、実測値、補正後の値の各項目を有する。正しい値は、校正用イオンと、夾雑物イオンとの取り違えが生じていない場合における四重極RF電圧である。また、補正後の値は、前記した校正点補正が行われた後における四重極RF電圧である。取得された校正点521が正しければ(ピークの取り違えが生じていなければ)、正しい値=実測値=補正後の値となる。
(Specific example)
Next, a specific example of the mass spectrometry data processing method shown in FIG. 8 will be described with reference to FIGS. 6 and 8 to 9B.
FIG. 9A is a diagram showing a table of specific values of the calibration points 521.
The table shown in FIG. 9A has items for the number of the
また、図9Aに示す表の最後の行には、標準偏差が示されている。標準偏差は、正しい値、実測値、補正後の値のそれぞれについて算出されている。なお、図9Aに示す標準偏差は、校正点間の傾き(式(1)で示されるai)の標準偏差である。 The last row of the table in Figure 9A shows the standard deviation. The standard deviation is calculated for each of the correct value, the measured value, and the corrected value. Note that the standard deviation shown in Figure 9A is the standard deviation of the slope between the calibration points (ai shown in equation (1)).
図9Aに示す表の例では、校正点251の番号「6」において夾雑物イオンとの取り違えが生じているため、正しい値と、実測値が異なっている。 In the example table shown in Figure 9A, the correct value differs from the actual measured value because the calibration point 251 number "6" has been confused with an impurity ion.
図9Bは、図9Aに基づく偏差(Di)のグラフを示す図である。
図9Bに示すグラフでは、横軸が偏差の番号、縦軸が偏差の大きさである。
図9Aの例で示すように、校正点521(図6参照)の番号「6」で取り違えが発生しているため、偏差の番号「5」と「6」(即ち、番号「6」の校正点521の両側の傾き)の偏差が大きくなっている。
FIG. 9B is a diagram showing a graph of the deviation (Di) based on FIG. 9A.
In the graph shown in FIG. 9B, the horizontal axis represents the deviation number, and the vertical axis represents the magnitude of the deviation.
As shown in the example of Figure 9A, a mix-up has occurred at calibration point 521 (see Figure 6) number "6", resulting in a large deviation between deviation numbers "5" and "6" (i.e., the slope on both sides of
ここで、図9Aのような校正点521のデータが得られた場合における校正点補正の具体的な手順を示す。
本実施形態における補正は偏差のずれが大きい校正点521を以下の式を用いて線形補正することで行われる(適宜、図7Bを参照)。
Here, a specific procedure for correcting the calibration point when data of the
In this embodiment, the correction is performed by linearly correcting the
まず、1番目(最初)の偏差(D1)に基づいて正しい四重極RF電圧の値(x1)が算出される場合、以下の式(11)で正しい四重極RF電圧の値(x1)が算出される。 First, when the correct quadrupole RF voltage value (x1) is calculated based on the first deviation (D1), the correct quadrupole RF voltage value (x1) is calculated using the following formula (11).
D1: x1=a2・m1 + b2・・・(11) D1: x1=a2・m1+b2...(11)
また、N番目(最後)の偏差(DN)に基づいて正しい四重極RF電圧の値(xN)が算出される場合、以下の式(12)で正しい四重極RF電圧の値(xN)が算出される。 Furthermore, when the correct quadrupole RF voltage value (xN) is calculated based on the Nth (last) deviation (DN), the correct quadrupole RF voltage value (xN) is calculated using the following formula (12).
DN: xN=aN-1・mN + bN-1・・・(12) DN: xN=aN-1・mN + bN-1...(12)
i(i=2~N-1)番目の偏差Diが偏差の閾値より大きい場合、以下の式(13)、(14)により正しい四重極RF電圧の値(xi)が算出される。 If the i-th (i = 2 to N-1) deviation Di is greater than the deviation threshold, the correct quadrupole RF voltage value (xi) is calculated using the following equations (13) and (14).
Di+1>Di-1の場合:
xi+1=ai+2・mi+1 + bi+2・・・(13)
Di+1≦Di-1の場合:
xi=ai-1・mi + bi-1・・・(14)
When Di+1>Di-1:
xi+1=ai+2・mi+1+bi+2...(13)
If Di+1≦Di−1:
xi=ai-1・mi+bi-1...(14)
式(11)~(14)において、xi,mi(i=1,・・・,N)は、図7A及び図7Bに示されるものであり、ai(i=1,・・・,N)は、式(1)によって算出されるものである。また、bi(i=1,・・・,N)は、図7A、図7Bに示すm/z-四重極RF電圧関係グラフ上の座標(mi、xi)を通り、傾き(ai)を有する直線における四重極RF電圧軸の切片である。このように、補正処理部212は、傾き(ai)を基に、校正点512の補正を行う。
In equations (11) to (14), xi and mi (i = 1, ..., N) are shown in Figures 7A and 7B, and ai (i = 1, ..., N) is calculated by equation (1). Furthermore, bi (i = 1, ..., N) is the intercept of the quadrupole RF voltage axis of a straight line that passes through the coordinates (mi, xi) on the m/z-quadrupole RF voltage relationship graph shown in Figures 7A and 7B and has a slope (ai). In this way, the
以下では、図8のステップS122~S141で行われる処理について図7A、図7B、図9A、図9Bを参照して具体的に説明する。なお、以下の記載において、ステップ番号は図8に示す処理のステップ番号である。 Below, the processing performed in steps S122 to S141 in FIG. 8 will be described in detail with reference to FIGS. 7A, 7B, 9A, and 9B. Note that in the following description, the step numbers refer to the step numbers of the processing shown in FIG. 8.
前記したように、例として、図9Aの番号「6」の校正点521(m/z:1172.145)において、ピークの取り違えが生じた場合の動作を以下に説明する。図9Aの番号「6」では、統一原子質量単位で約1.0Da(1.66×10-24g)高い夾雑物イオンのピークとの取り違えが生じている。 As described above, the operation when a peak is confused at calibration point 521 (m/z: 1172.145) of number "6" in Fig. 9A will be described below as an example. At number "6" in Fig. 9A, the peak is confused with the peak of an impurity ion that is about 1.0 Da (1.66 x 10-24 g) higher in unified atomic mass units.
(S122)補正処理部212は、すべての偏差の中から最も大きな値を有する偏差(Di)を検索する。ここで、iは図9Bに示す偏差の番号である。図9Bに示す例では、符号601で示されるD5(i=5)=0.186が最大である。
(S123)補正処理部212は、ステップS122で検索された偏差の値が偏差のしきい値(Dth)以上であるか否かを判定する。図9Bに示すD5(符号601)は偏差のしきい値以上であるものとする。
(S122) The
(S123) The
(S124)補正処理部212は、最も大きな値を有する偏差(Dmax)に該当する校正点521を置換することによって校正点521の補正を行う。具体的には、最大の偏差(Di)に対応する傾き(ai)の校正点521(mi又はmi+1)の四重極RF電圧(xi又はxi+1)を、前後の校正点521から補正した値に置換する。置換は、式(11)~(14)に従って算出される四重極RF電圧によって行われる。
(S124) The
図9Bに示す例では、D5=0.186が偏差の最大であり、D6>D4であるので、補正処理部212は、番号「6」の校正点521の四重極RF電圧(xi+1)をxi+1=ai+2・mi+1 + bi+2(式(13))で求めた値に置換する。これによって、図9Aに示すように、番号「6」を有する校正点521の四重極RF電圧は「111958」に補正される。
In the example shown in FIG. 9B, D5 = 0.186 is the maximum deviation, and D6 > D4, so the
(S131)補正処理部212は、ステップS124で補正後の四重極RF電圧を用いて傾き(ai)を再計算する。
(S132)補正処理部212は、ステップS131で再計算した傾きの標準偏差(SD)を再計算する。図9Aに示す例の場合、補正後の標準偏差(SD)は「0.099」となる。
(S131) The
(S132) The
(S133)補正処理部212は、補正後の標準偏差を標準偏差のしきい値(Sth)と比較する。図9Aに示す例の場合、補正後の標準偏差(SD)は、前記したように「0.099」で標準偏差のしきい値(Sth=0.1714)以下となる。
(S133) The
ステップS133で補正後の標準偏差がしきい値以下となった場合(S133→No)、測定処理部213はステップS141の処理を行う。ステップS141において、測定処理部213は、補正した校正点521について、補正後の四重極RF電圧を初期値として電圧校正を実施する(第3のステップ)。具体的には、測定処理部213は、補正後の四重極RF電圧を初期値として、前記したスキャンを行うことでスペクトル501を再取得する。そして、測定処理部213は、再取得したスペクトル501を基に、補正によって、スペクトル501aで示される校正用イオンのピーク値に対応する四重極RF電圧を推定する。このようにして、正しい四重極RF電圧を取得することができる。なお、ステップS141で行われる処理については後記する。
If the corrected standard deviation is equal to or less than the threshold value in step S133 (S133→No), the
また、ステップS133で、補正後の標準偏差がしきい値より大きくなった場合(S133→Yes)、補正処理部212は、標準偏差がしきい値以下になるまでステップS112以降の処理を繰り返す。ただし、ループ回数が所定回数(hmax)より大きくなった場合、補正処理部212はエラー出力を行う(S151)。
Also, in step S133, if the corrected standard deviation is greater than the threshold value (S133 → Yes), the
図9Aを参照すると、番号「6」の校正点521について、補正後の四重極RF電圧は、前記したように「111958」(V)である。この四重極RF電圧と、正しい四重極RF電圧「111977」(V)との差は、わずか19(V)(=約0.2Da(0.332×10-24g))である。これにより、補正後の四重極RF電圧を初期値とした電圧校正では、前記したように正しい校正用イオンに由来するピーク値を取得することができる。また、図9Bの破線及び白丸に示すように、偏差も小さくなる。
9A, for
校正点間の傾きを示す式(1)では、m/z差(mi+1 - mi)が小さいほど四重極RF電圧の振幅のずれを鋭敏に反映する。校正点521が均一のm/zの間隔となるよう校正用試料が選ばれる場合が多い。しかし、m/zの間隔が均一ではない校正用試料が選ばれる場合がある。このように、校正点521m/zの間隔が均一でない場合、標準偏差を算出する際に校正点間のm/z差(mi+1 - mi)に重みづけが行われてもよい。 In equation (1), which shows the slope between calibration points, the smaller the m/z difference (mi+1 - mi), the more sensitively it reflects the deviation in the amplitude of the quadrupole RF voltage. In many cases, a calibration sample is selected so that the calibration points 521 have uniform m/z spacing. However, there are cases where a calibration sample is selected in which the m/z spacing is not uniform. In this way, when the m/z spacing of the calibration points 521 is not uniform, the m/z difference (mi+1 - mi) between the calibration points may be weighted when calculating the standard deviation.
また、本実施形態では、校正点521の補正が行われた後、電圧校正(図8のS141)が行われているが、これに限らない。校正点521の補正が行われた後、補正の結果(つまり、補正後の四重極RF電圧)のみが記憶部220に保存され、電圧校正は校正点521の補正とは別の処理として行われてもよい。
In addition, in this embodiment, after the
(電圧校正)
図10は、電圧校正の手法を示す図である。適宜、図1を参照する。
図10では、図5に示す四重極RF電圧-イオン信号強度関係グラフのうち、スペクトル501aと、スペクトル502の近傍を拡大した図である。また、図10に示す処理は、図1に示す測定処理部213が、図8のステップS141で行うものである。
補正により校正点521が補正された後の電圧校正では、補正した校正点521によって算出される四重極RF電圧(符号531)を中心とした狭いスキャン範囲511aで測定が開始される。つまり、校正点521の補正によって推測される四重極RF電圧(符号531)を中心として、スキャン範囲511aのように狭いスキャン範囲511から徐々にスキャン範囲511が広げられる。なお、図9Aの番号「6」では、校正点521の補正による四重極RF電圧と、正しい四重極RF電圧「111977」(V)との差は、わずか19(V)であることが示されている。換言すれば、校正点521の補正によって推定される四重極RF電圧は、電圧校正の対象となるスペクトル501aにおける真の中心とは限らない。そのため、図10では、補正した校正点521によって算出される四重極RF電圧(符号531)が、スペクトル501aの中心からずれた位置に示されている。
(Voltage calibration)
10 is a diagram showing a voltage calibration method, with reference to FIG.
Fig. 10 is an enlarged view of
In the voltage calibration after the
スキャン範囲551の変更は、測定が行われるごとに行われる。つまり、スキャン範囲511a~スキャン範囲511cへと徐々にスキャン範囲511が広げられる。そして、校正用イオンのスペクトル501aがスキャン範囲511cに入るまで徐々にスキャン範囲511が拡大される。このように狭いスキャン範囲511から測定が開始されることで、スペクトル502aをスキャンすることを防止することができる。従って、スペクトル501aのピークが、隣接する夾雑物イオンによるスペクトル502のピークと取り違えられることを防止することができる。また、スペクトル501aの真の中心を得ることが可能となる。
The scan range 551 is changed each time a measurement is performed. In other words, the
このように、電圧校正において、測定処理部213は、補正した校正点521によって算出される四重極RF電圧を中心とするスキャン範囲511で、質量分析装置100(図による校正用イオンの測定を行う。この際、測定処理部213は、図10のスキャン範囲511a~511cに示すように、測定が行われるたびに、スキャン範囲511が、徐々に広い範囲となるようスキャン範囲511を変更する。
In this way, in voltage calibration, the
なお、校正点521によってスペクトル501aのピーク値に対応する四重極EF電圧が推測されているにもかかわらず、図10に示すような再スキャンが行われている。これは、前記したように、図9Aの校正点521の番号「6」に示すように、補正によって推測される四重極RF電圧値は、必ずしも正しい値を示しているとは限らないためである。
Note that even though the quadrupole EF voltage corresponding to the peak value of
第1実施形態によれば、電圧校正を行う際に、校正点521と隣接するm/zを有する夾雑物イオンが存在している場合でも、夾雑物イオンと校正用イオンのピークの取り違えを補正することができる。これによって、校正用イオンに由来する、正確な四重極RF電圧の決定が可能となる。この結果、未知の物質を質量分析する際、精度の高い質量分析が可能となる。また、第1実施形態によれば、ピークの取り違えが1つのみでも校正点521の補正、及び、電圧校正を行うことができる。
According to the first embodiment, when performing voltage calibration, even if there is an impurity ion having an m/z adjacent to the
<第2実施形態>
以下に第2実施形態について説明する。なお、第2実施形態における質量分析システム1の構成は、図1に示す構成と同様であるため、第2実施形態での図示及び説明を省略する。図12において、図7Bと同様の構成については同一の符号を付し、説明を省略する。
Second Embodiment
The second embodiment will be described below. Note that the configuration of the
以下、図11及び図12を参照して、第2実施形態で行われる処理について説明する。
図11は、第2実施形態に係る質量分析データ処理方法の手順を示すフローチャートである。また、図12は、m/z-四重極RF電圧関係グラフを示す図である。以下の説明において、ステップ番号は、図11に示すステップ番号である。また、図11において、図8と同様のステップについては、同一のステップ番号を付し、適宜説明を省略する。
第2実施形態において、補正処理部212は、校正点521のそれぞれに対して最小二乗法を適用することで傾き及び切片711を算出する。そして、補正処理部212は、算出した傾きに基づいて標準偏差としきい値を比較する。つまり、第1実施形態では、図8のステップS102及びステップS131では、校正点521の間の傾き(ak)が1つ1つ算出されていた。これに対し、第2実施形態において、補正処理部212は、以下の式(21)による最小二乗法によって、図12に示すように、すべての校正点521に近接する直線701を求める。つまり、補正処理部212は、校正点521のそれぞれに対して最小二乗法を適用することで、直線701の傾き及び切片711を算出する(S102A及びS131A)。つまり、直線701は、校正点521のそれぞれに対して最小二乗法を適用することで算出される傾き及び切片711を有する。このように、第2実施形態では、それぞれの校正点521の間の傾きが、最小二乗法で求められる傾きで統一される。
Hereinafter, the processing performed in the second embodiment will be described with reference to FIG. 11 and FIG.
Fig. 11 is a flow chart showing the procedure of the mass spectrometry data processing method according to the second embodiment. Fig. 12 is a graph showing the relationship between m/z and quadrupole RF voltage. In the following description, the step numbers are those shown in Fig. 11. In Fig. 11, the same step numbers are used for steps that are the same as those in Fig. 8, and descriptions thereof will be omitted as appropriate.
In the second embodiment, the
続いて、補正処理部212は、偏差(Di)を算出する(S103A)。第2実施形態において、偏差(Di)は、それぞれの校正点521から最小二乗法で求められる直線701への垂線721の長さとして定義される。垂線721の長さは、校正点521のそれぞれに対する傾き及び切片711を有する直線701との距離である。このように、ステップS103Aにおいて、補正処理部212は、校正点512のそれぞれに対して、傾き及び切片711を有する直線701との距離(垂線721の長さ)を偏差(Di)として算出する。
Then, the
その後、補正処理部212は、ステップS103Aで算出した偏差(Di)の標準偏差(SD)を算出する(S104A)。
Then, the
そして、補正処理部212は、偏差(Di)の標準偏差(SD)が予め設定されている閾値(SD>Sth)より大きい場合(S105→Yes)、傾きを基に、校正点521の補正を行う。傾きを基に、校正点521の補正を行う手順は、第1実施形態と同様である。
Then, if the standard deviation (SD) of the deviation (Di) is greater than a preset threshold (SD>Sth) (S105→Yes), the
式(21)において、「slope」は傾きを示し、mi,xiは、図7A及び図7Bに示すようにm/z及び四重極RF電圧を示す。また、mの上にバーが示されている変数はmiの平均を示す。また、xの上にバーが示されている変数はxi(四重極RF電圧)の平均を示す。 In equation (21), "slope" indicates the slope, and mi and xi indicate the m/z and quadrupole RF voltage as shown in Figures 7A and 7B. Furthermore, the variable with a bar above m indicates the average of mi. Furthermore, the variable with a bar above x indicates the average of xi (quadrupole RF voltage).
そして、校正点521の補正は、以下の式(22)によって行われる。
Then, the correction of the
xi_cor=slope・mi+b・・・(22) xi_cor=slope・mi+b...(22)
式(22)において、xi_corは、補正された四重極RF電圧の値を示し、bは最小二乗法によって算出される四重極RF電圧軸の切片711である。すなわち、補正された校正点521は、最小二乗法による直線701に乗るよう補正される。
In equation (22), xi_cor indicates the value of the corrected quadrupole RF voltage, and b is the
なお、図11のステップS121Aでは、ステップS103Aと同様の手順で偏差(Di)が算出される。また、図11では図示省略されているが、ステップS121Aの前段ではステップS102Aと同様の手順で傾き及び切片711が算出される。また、ステップS131AはステップS102Aと同様の処理であり、ステップS132AはステップS104Aと同様の処理である。さらに、図示省略されているが、ステップS131Aと、ステップS132Aとの間には、ステップS103Aと同様の処理が行われる。 In step S121A in FIG. 11, the deviation (Di) is calculated in the same manner as in step S103A. Although not shown in FIG. 11, before step S121A, the slope and intercept 711 are calculated in the same manner as in step S102A. Furthermore, step S131A is the same process as step S102A, and step S132A is the same process as step S104A. Furthermore, although not shown, between step S131A and step S132A, the same process as in step S103A is performed.
このように、第2実施形態では、校正点521について最小二乗法を用いることにより算出される傾きを基に夾雑物イオンと校正用イオンのピークの取り違えが補正される。
In this way, in the second embodiment, the misidentification of the peaks of impurity ions and calibration ions is corrected based on the slope calculated by using the least squares method for the
そして、第2実施形態に示す手法では、校正点521の間の傾きを1つ1つ算出する第1実施形態と比較して、校正点521の間の傾きが最小二乗法によって1度に算出される。そのため、第2実施形態に示す手法は、第1実施形態に示す手法と比較して計算が単純となり高速に実施出来る。 In the method shown in the second embodiment, the slope between the calibration points 521 is calculated at once by the least squares method, compared to the first embodiment in which the slope between the calibration points 521 is calculated one by one. Therefore, the method shown in the second embodiment has simpler calculations and can be performed quickly compared to the method shown in the first embodiment.
<第3実施形態>
次に、図13を参照して、分析部103に飛行時間型質量分析装置400が設置されている質量分析システム1aに、第1実施形態及び第2実施形態に示す校正点補正処理を適用する例を示す。
図13は、第3実施形態に係る質量分析システム1aの構成例を示す図である。
図13に示す質量分析システム1aに示す質量分析装置100aが図1に示す質量分析システム1と頃なる点は以下の通りである。
(A1)分析部103に飛行時間型質量分析装置400が設けられている。
(A2)飛行時間型質量分析装置400を構成する押出電極401にパルス電圧を印加する電源410が接続されている。この電源410は、電圧制御装置300によって制御される。
図13において、図1と同様の構成については、同一の符号を付す。
Third Embodiment
Next, referring to FIG. 13, an example will be shown in which the calibration point correction process shown in the first and second embodiments is applied to a
FIG. 13 is a diagram showing an example of the configuration of a
A
(A1) The
(A2) A
In FIG. 13, the same components as those in FIG. 1 are denoted by the same reference numerals.
質量分析装置100aにおいて、イオン源151で生成されたイオンは、図1と同様、細孔121を介して第1差動排気部101に導入される。そして、図1と同様、第1差動排気部101は、ポンプ111で排気されている。これにより、第1差動排気部101は、10Paから500Paの真空度に維持されている。第1差動排気部101を通過したイオンは、細孔122を介して第2差動排気部102に導入される。図1と同様、第2差動排気部102は、ポンプ112で排気されている。これにより、第2差動排気部102は、0.1Pa~10Paの真空度に維持されている。また、図1と同様、第2差動排気部102には、イオンを収束させるイオンガイド130が設置されている。イオンガイド130で収束されたイオンは細孔123を通過して飛行時間型質量分析装置400が設置されている分析部103に導入される。図1と同様、分析部103は、ポンプ113で排気されている。これによって、分析部103は、1E-5Pa以下の圧力に維持されている。
In the
飛行時間型質量分析装置400は、押出電極401、加速電極402、リフレクタ電極403および検出器152から構成される。押出電極401には電源410によってパルス電圧が印可される。パルス電圧が印加されることにより、押出電極401は、イオンガイド130から入射されたイオンの一部の進行方向を、イオンの入射方向に対して垂直方向に屈曲させる。屈曲されたイオンは加速電極402の方へ向かう。そして、進行方向を屈曲されたイオンは加速電極402で加速され、リフレクタ電極403で反射される(図13の白抜矢印)。リフレクタ電極403で反射されイオンは、検出器152で検出される。押出電極401によって屈曲されたイオンが検出器152に到達するまでの時間、つまり、イオンの飛行時間は、イオンのm/zに依存する。そのため、飛行時間と検出器152で検出されたイオンの信号強度をプロットすることでイオンのスペクトル501(図5参照)を得ることができる。
The time-of-
飛行時間型質量分析装置400では検出器152として、マルチチャンネルプレート型の検出器152が一般に用いられている。検出器152によりイオンの検出強度は電気的な信号(出力信号)に変換され、データ処理装置200へと送られる。この際、図1と同様、検出器152から出力された出力信号は、一定のサンプリング周期のデジタルデータに変換された後、データ処理装置200へ送られる。図1と同様、デジタルデータへの変換は、図示しないアナログ-デジタルコンバータ(ADC)や、図示しないパルスカウンティングユニットで行われる。データ処理装置200は、送られたデジタルデータを記憶部220に蓄積する。
In the time-of-
図13に示すような質量分析システム1aでも、校正用イオンを測定し、得られたスペクトル501(図5参照)について、第1実施形態や、第2実施形態に示す手法を適用することができる。この際、第1実施形態の四重極RF電圧の代わりに、押出電極401に印加されるパルス電圧の振幅が用いられる。つまり、第3実施形態では、パルス電圧が、複数の校正用イオンを質量分析装置100で測定した結果として取得される電圧である。また、第3実施形態において、図6~図7Bに示されるm/z-四重極RF電圧関係グラフの四重極RF電圧軸は、パルス電圧軸となる。
In the
質量分析システム1aにおいて、校正点521に対してm/zが隣接する夾雑物イオンが存在している場合、夾雑物イオンのピークと校正用イオンのピークの取り違えが生じる可能性がある。しかしながら、第1実施形態の四重極RF電圧の代わりにパルス電圧をもちいることで、第1実施形態及び第2実施形態と同様の処理を行うことができる。これにより、ピークの取り違えを補正することができる。これにより、未知の物質を質量分析した際に正確なm/zを得ることができる。
In the
本発明は前記した実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、前記した実施形態は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明したすべての構成を有するものに限定されるものではない。また、ある実施形態の構成の一部を他の実施形態の構成に置き換えることが可能であり、ある実施形態の構成に他の実施形態の構成を加えることも可能である。また、各実施形態の構成の一部について、他の構成の追加・削除・置換をすることが可能である。 The present invention is not limited to the above-described embodiments, and includes various modified examples. For example, the above-described embodiments have been described in detail to clearly explain the present invention, and are not necessarily limited to having all of the configurations described. It is also possible to replace part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment. It is also possible to add, delete, or replace part of the configuration of each embodiment with other configurations.
また、前記した各構成、機能、プロット処理部211、補正処理部212、測定処理部213、記憶部220等は、それらの一部又はすべてを、例えば集積回路で設計すること等によりハードウェアで実現してもよい。また、図2に示すように、前記した各構成、機能等は、CPU等のプロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、HD(Hard Disk)に格納すること以外に、メモリや、SSD(Solid State Drive)等の記録装置、又は、IC(Integrated Circuit)カードや、SD(Secure Digital)カード、DVD(Digital Versatile Disc)等の記録媒体に格納することができる。
Furthermore, the above-mentioned configurations, functions,
また、各実施形態において、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしもすべての制御線や情報線を示しているとは限らない。実際には、ほとんどすべての構成が相互に接続されていると考えてよい。 In addition, in each embodiment, the control lines and information lines shown are those that are considered necessary for the explanation, and not all control lines and information lines in the product are necessarily shown. In reality, it can be considered that almost all components are interconnected.
1,1a 質量分析システム
100,100a 質量分析装置
140 四重極質量フィルタ
141,141a~141d 四重極ロッド電極
152 検出器
200 データ処理装置
211 プロット処理部
212 補正処理部
213 測定処理部
400 飛行時間型質量分析装置
501 スペクトル
501a スペクトル
502 スペクトル
511,511a~511c スキャン範囲(電圧範囲)
521 校正点
521a 校正点(校正用イオンによる校正点)
521b 校正点(外れ値)
701 直線(傾きを含む)
711 切片
721 垂線(偏差)
S101 プロット処理(第1のステップ)
S124 補正(第2のステップ)
S141 質量校正(第3のステップ)
1, 1a
521
521b Calibration point (outlier)
701 Straight line (including inclination)
S101 Plot process (first step)
S124 Correction (Second step)
S141 Mass calibration (third step)
Claims (10)
それぞれの前記校正点に関して外れ値があれば、前記外れ値を補正することで、前記校正点の補正を行う補正処理部と、
を有するデータ処理装置。 a plot processing unit that plots a voltage obtained as a result of measuring a plurality of calibration ions, each having a known m/z, by a mass spectrometer, and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes;
a correction processing unit that corrects the calibration points by correcting the outliers if there are outliers for the calibration points;
A data processing device having:
を有することを特徴とする請求項1に記載のデータ処理装置。 2. The data processing apparatus according to claim 1, further comprising a measurement processing unit that performs measurement of the calibration ions by the mass spectrometer within a voltage range centered on the voltage calculated based on the corrected calibration point.
前記測定が行われるたびに、前記電圧範囲が、徐々に広い範囲となるよう前記電圧範囲を変更する
ことを特徴とする請求項2に記載のデータ処理装置。 The measurement processing unit includes:
3. The data processing apparatus according to claim 2, wherein the voltage range is changed so that the voltage range becomes gradually wider each time the measurement is performed.
前記校正点それぞれの間の傾きを算出し、
それぞれの前記傾きに関する偏差を算出し、
前記偏差の標準偏差が、予め設定されている閾値より大きい場合、前記外れ値があると判定し、前記傾きを基に、前記校正点の補正を行う
ことを特徴とする請求項1に記載のデータ処理装置。 The correction processing unit:
Calculating a slope between each of the calibration points;
Calculating the deviation for each of the slopes;
2 . The data processing device according to claim 1 , wherein, when the standard deviation of the deviations is greater than a preset threshold value, it is determined that the outlier exists, and the calibration point is corrected based on the slope.
前記校正点のそれぞれに対して最小二乗法を適用することで、傾き及び切片を算出し、
前記校正点のそれぞれに対して、前記傾き及び前記切片を有する直線との距離を偏差として算出し、前記偏差の標準偏差が予め設定されている閾値より大きい場合、前記傾きを基に、前記校正点の補正を行う
ことを特徴とする請求項1に記載のデータ処理装置。 The correction processing unit:
Calculating a slope and an intercept by applying a least squares method to each of the calibration points;
2. The data processing device according to claim 1, further comprising: a calculation unit that calculates a deviation for each of the calibration points from a straight line having the slope and the intercept; and, if a standard deviation of the deviation is greater than a preset threshold, corrects the calibration points based on the slope.
m/zが既知である、複数の校正用イオンを質量分析装置で測定した結果として取得される電圧と、それぞれの校正用イオンの前記m/zと、を、前記電圧及び前記m/zのそれぞれを座標軸として有する座標に校正点としてプロットするプロット処理部と、
それぞれの前記校正点に関して外れ値があれば、前記外れ値を補正することで、前記校正点の補正を行う補正処理部と、
を有する質量分析システム。 A mass spectrometer;
a plot processing unit that plots a voltage obtained as a result of measuring a plurality of calibration ions, each having a known m/z, by a mass spectrometer, and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes;
a correction processing unit that corrects the calibration points by correcting the outliers if there are outliers for the calibration points;
A mass spectrometry system comprising:
を有することを特徴とする請求項6に記載の質量分析システム。 7. The mass spectrometry system according to claim 6, further comprising a measurement processing unit that performs measurement of the calibration ions by the mass spectrometer within a voltage range centered on the voltage calculated based on the corrected calibration point.
分析部に、四重極質量フィルタ又は飛行時間型質量分析装置
を有することを特徴とする請求項6に記載の質量分析システム。 The mass spectrometer is
7. The mass spectrometry system according to claim 6, wherein the analysis section comprises a quadrupole mass filter or a time-of-flight mass analyzer.
m/zが既知である、複数の校正用イオンを質量分析装置で測定した結果として取得される電圧と、それぞれの校正用イオンの前記m/zと、を、前記電圧及び前記m/zのそれぞれを座標軸として有する座標に校正点としてプロットする第1のステップと、
それぞれの前記校正点に関して外れ値があれば、前記外れ値を補正することで、前記校正点の補正を行う第2のステップと、
を実行する質量分析データ処理方法。 A data processing device comprising:
a first step of plotting a voltage obtained as a result of measuring a plurality of calibration ions, each having a known m/z, by a mass spectrometer and the m/z of each calibration ion as calibration points on a coordinate system having the voltage and the m/z as coordinate axes;
a second step of correcting the calibration points by correcting any outliers associated with the calibration points;
A mass spectrometry data processing method for performing
を実行することを特徴とする請求項9に記載の質量分析データ処理方法。 10. The method for processing mass spectrometry data according to claim 9, further comprising: a third step of measuring the calibration ions by the mass spectrometer within a voltage range centered on the voltage calculated based on the corrected calibration point.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2023108722 | 2023-06-30 | ||
JP2023-108722 | 2023-06-30 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2025004999A1 true WO2025004999A1 (en) | 2025-01-02 |
Family
ID=93938599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2024/022568 WO2025004999A1 (en) | 2023-06-30 | 2024-06-21 | Data processing device, mass spectrometry system, and mass spectrometry data processing method |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2025004999A1 (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015503825A (en) * | 2011-12-27 | 2015-02-02 | ディーエイチ テクノロジーズ デベロップメント プライベート リミテッド | Method for extracting ions having a low M / Z ratio from an ion trap |
US20170221692A1 (en) * | 2014-08-05 | 2017-08-03 | DH Technologies Development Pte Ltd. | Band pass extraction from an ion trapping device and tof mass spectrometer sensitivity enhancement |
JP2019086467A (en) * | 2017-11-09 | 2019-06-06 | 株式会社島津製作所 | MASS SPECTROMETRY CONTROLLER, MASS SPECTROMETER, MASS SPECTROMETRY CONTROL METHOD, AND MASS SPECTROMETRY METHOD |
WO2020066009A1 (en) * | 2018-09-28 | 2020-04-02 | 株式会社島津製作所 | Mass spectrometry device evaluation method, mass spectrometry device calibration method, analysis method, mass spectrometry device, and mass spectrometry reagent |
WO2021220671A1 (en) * | 2020-04-28 | 2021-11-04 | 株式会社日立ハイテク | Mass spectrometry device control method, mass spectrometry system, and voltage control device |
-
2024
- 2024-06-21 WO PCT/JP2024/022568 patent/WO2025004999A1/en unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015503825A (en) * | 2011-12-27 | 2015-02-02 | ディーエイチ テクノロジーズ デベロップメント プライベート リミテッド | Method for extracting ions having a low M / Z ratio from an ion trap |
US20170221692A1 (en) * | 2014-08-05 | 2017-08-03 | DH Technologies Development Pte Ltd. | Band pass extraction from an ion trapping device and tof mass spectrometer sensitivity enhancement |
JP2019086467A (en) * | 2017-11-09 | 2019-06-06 | 株式会社島津製作所 | MASS SPECTROMETRY CONTROLLER, MASS SPECTROMETER, MASS SPECTROMETRY CONTROL METHOD, AND MASS SPECTROMETRY METHOD |
WO2020066009A1 (en) * | 2018-09-28 | 2020-04-02 | 株式会社島津製作所 | Mass spectrometry device evaluation method, mass spectrometry device calibration method, analysis method, mass spectrometry device, and mass spectrometry reagent |
WO2021220671A1 (en) * | 2020-04-28 | 2021-11-04 | 株式会社日立ハイテク | Mass spectrometry device control method, mass spectrometry system, and voltage control device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7202473B2 (en) | Mass spectrometer | |
US8431886B2 (en) | Estimation of ion cyclotron resonance parameters in fourier transform mass spectrometry | |
US20080201095A1 (en) | Method for Calibrating an Analytical Instrument | |
US20070136017A1 (en) | Method for calibrating mass spectrometry (ms) and other instrument systems and for processing ms and other data | |
CN113490994B (en) | Self-calibration of high-resolution mass spectrometry | |
US20250046592A1 (en) | Adaptive intrinsic lock mass correction | |
US7381948B2 (en) | Mass spectroscope and method of calibrating the same | |
US11573213B2 (en) | Method for time-alignment of chromatography-mass spectrometry data sets | |
WO2019150576A1 (en) | Mass spectroscope and mass calibration method for mass spectroscope | |
US11515139B2 (en) | Method for determining a parameter to perform a mass analysis of sample ions with an ion trapping mass analyser | |
JP2005521030A (en) | Calibration method of mass spectrometer | |
CA2464004C (en) | Mass spectrometer | |
CN115461838A (en) | Laser desorption ionization mass analysis device and laser power adjustment method | |
US6623935B2 (en) | Deconvolution method and apparatus for analyzing compounds | |
EP2483641B1 (en) | Systems and methods for maintaining the precision of mass measurement | |
WO2025004999A1 (en) | Data processing device, mass spectrometry system, and mass spectrometry data processing method | |
US12159776B2 (en) | Quadrupole mass spectrometer data to enable new hardware operating regimes | |
JP4950029B2 (en) | Mass spectrometer | |
EP4473537A1 (en) | Compound identification using a mass spectrum | |
JP6591565B2 (en) | Mass spectrometer and ion detection method thereof | |
WO2024161624A1 (en) | Mass spectrometry method and mass spectrometry device | |
CN118817868A (en) | Mass analysis method and mass analysis device | |
EP3664124A1 (en) | Spectrum calculation processing device, spectrum calculation processing method, spectrum calculation processing program, ion trap mass spectrometry system and ion trap mass spectrometry method | |
WO2022091047A1 (en) | Compound identification by mass spectrometry | |
CN116959951A (en) | Quality analysis device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 24831860 Country of ref document: EP Kind code of ref document: A1 |