JP7077481B2 - Mass spectrometer and mass spectrometry method - Google Patents

Mass spectrometer and mass spectrometry method Download PDF

Info

Publication number
JP7077481B2
JP7077481B2 JP2021514815A JP2021514815A JP7077481B2 JP 7077481 B2 JP7077481 B2 JP 7077481B2 JP 2021514815 A JP2021514815 A JP 2021514815A JP 2021514815 A JP2021514815 A JP 2021514815A JP 7077481 B2 JP7077481 B2 JP 7077481B2
Authority
JP
Japan
Prior art keywords
value
ionic strength
mass
unit
mass spectrometer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2021514815A
Other languages
Japanese (ja)
Other versions
JPWO2020213283A1 (en
Inventor
博教 山下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi High Tech Corp
Original Assignee
Hitachi High Tech Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi High Tech Corp filed Critical Hitachi High Tech Corp
Publication of JPWO2020213283A1 publication Critical patent/JPWO2020213283A1/en
Application granted granted Critical
Publication of JP7077481B2 publication Critical patent/JP7077481B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/62Investigating 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
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01JELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J49/00Particle spectrometers or separator tubes
    • H01J49/26Mass spectrometers or separator tubes
    • H01J49/34Dynamic spectrometers
    • H01J49/42Stability-of-path spectrometers, e.g. monopole, quadrupole, multipole, farvitrons

Description

本発明は、質量分析装置および質量分析方法に関する。 The present invention relates to a mass spectrometer and a mass spectrometry method.

質量分析装置で試料成分を定量する場合、目的の成分に由来するイオンのm/z値に着目し、そのイオン量の経時変化を測定する選択イオンモニタリング(SIM:Selected Ion Monitoring)や選択反応モニタリング(SRM:Selected Reaction Monitoring)と呼ばれる測定方法が一般に利用される。 When quantifying sample components with a mass spectrometer, focus on the m / z value of ions derived from the target component, and select ion monitoring (SIM: Selected Ion Monitoring) or selective reaction monitoring to measure the change in the amount of ions over time. A measurement method called (SRM: Selected Reaction Monitoring) is generally used.

これらの方法は、目的のm/z値に対し装置の制御条件を固定しているため、目的のイオンに対する感度は良いがm/z値のずれを直接測定することができない。 In these methods, since the control conditions of the apparatus are fixed to the target m / z value, the sensitivity to the target ion is good, but the deviation of the m / z value cannot be directly measured.

あらかじめm/z値が既知の試料を測定し装置を校正しておく必要があるが、温度や湿度などの変化が電極の部材や電気回路などに影響を与え、校正したm/z値がずれる場合がある。 It is necessary to measure a sample with a known m / z value and calibrate the device in advance, but changes in temperature and humidity affect the electrode members and electric circuits, and the calibrated m / z value shifts. In some cases.

従って、SIMやSRMによる測定の前後に試料そのものまたはm/z値が既知の成分をスキャン測定し、m/z値がずれていないことを確認している。もしくは、定量測定中に、SIMやSRMによる測定と、既知成分のm/z値の近くのスキャン測定を時間的に分割することで、m/z値に関する情報を得ている。いずれの方法もSIMやSRMで測定しているm/z値がずれていないであろうことを間接的に証明しようとしている。 Therefore, before and after the measurement by SIM or SRM, the sample itself or the component whose m / z value is known is scanned and measured, and it is confirmed that the m / z value does not deviate. Alternatively, during the quantitative measurement, information on the m / z value is obtained by temporally dividing the measurement by SIM or SRM and the scan measurement near the m / z value of the known component. Both methods indirectly try to prove that the m / z value measured by SIM or SRM will not deviate.

例えば、特許文献1には、ロックマスと呼ばれる手法により、「外部検量物質又は参照化合物を用いて四重極マスフィルターの質量精度のずれを補正する手法」が記述されている(ここで外部検量物質や参照化合物はm/z値が既知の成分であり、質量精度のずれとは(ばらつきではなく)m/z値の真値からのずれを意味している)。 For example, Patent Document 1 describes "a method for correcting a deviation in mass accuracy of a quadrupole mass filter using an external calibration substance or a reference compound" by a method called rock mass (here, external calibration). A substance or a reference compound is a component having a known m / z value, and a deviation in mass accuracy means a deviation from the true value of the m / z value (not variation)).

特表2014-508937号公報Japanese Patent Publication No. 2014-508937

SIMやSRMによって測定する場合、試料成分に由来するイオンのm/z値に対し装置の制御条件を固定している。そのため、試料成分を正しく指定のm/z値で測定しているとの証憑を得られない。 When measuring by SIM or SRM, the control conditions of the device are fixed to the m / z value of the ions derived from the sample component. Therefore, it is not possible to obtain a voucher that the sample component is correctly measured at the specified m / z value.

特許文献1に記載の技術にあっては、質量精度のずれを補正する手法が開示されているが、SIMやSRMの測定結果から、正しく指定のm/z値で測定できているかを直接判定することはできない。つまり、温度や湿度の変化により目的のm/z値に対する装置制御条件が変化した場合、それを直接認識することが困難である。 In the technique described in Patent Document 1, a method for correcting a deviation in mass accuracy is disclosed, but it is directly determined from the measurement results of SIM or SRM whether or not the measurement can be performed correctly at a specified m / z value. Can't be done. That is, when the device control conditions for the target m / z value change due to changes in temperature and humidity, it is difficult to directly recognize them.

m/z値対する装置制御条件が変化した場合、本来測定できるはずのイオンを取り逃がし、感度の低下を招く。さらには、定量測定そのものを無駄にする状況もあり得る。 When the device control conditions for the m / z value change, ions that should be able to be measured are missed, resulting in a decrease in sensitivity. Furthermore, there may be situations where the quantitative measurement itself is wasted.

本発明の目的は、制御条件の変化によるm/z値のずれを検知し、定量分析に必要な感度を維持しつつ、m/z値のずれに対する証憑を与え、m/z値のずれにより減少したイオン量を補うことができる質量分析装置および質量分析方法を実現することである。 An object of the present invention is to detect a deviation of the m / z value due to a change in control conditions, to provide a voucher for the deviation of the m / z value while maintaining the sensitivity required for quantitative analysis, and to obtain a voucher for the deviation of the m / z value. It is to realize a mass spectrometer and a mass spectrometry method that can compensate for the reduced amount of ions.

上記目的を達成するため、本発明は次のように構成される。 In order to achieve the above object, the present invention is configured as follows.

試料成分をイオン化するイオン化部と、前記イオン化部により供給されたイオンを選別する質量分析部と、前記質量分析部により選別されたイオンの強度を測定する検出部と、前記イオン化部、前記質量分析部及び前記検出部の動作を制御する制御部とを備える質量分析装置において、前記制御部は、目的成分に由来するm/z値の近傍の複数個所のm/z値を選択するため、前記質量分析部の制御条件を切り替える質量分析制御部と、前記質量分析制御部が選択したm/z値のイオン強度を検出するイオン強度検出部と、前記イオン強度検出部が検出したイオン強度の変化を計算する成分量の変化計算部と、成分量一定の状態での前記複数個所のm/z値のイオン強度の比率を求め、近傍パターンを計算する近傍パターン計算部と、前記近傍パターンに基づき測定したイオン強度の不足を補填するイオン強度補正部と、を備え、目的成分に由来するm/z値のイオン強度の経時変化を測定する。 An ionization unit that ionizes sample components, a mass analysis unit that selects the ions supplied by the ionization unit, a detection unit that measures the intensity of the ions selected by the mass analysis unit, the ionization unit, and the mass analysis. In a mass analyzer including a unit and a control unit that controls the operation of the detection unit, the control unit selects m / z values at a plurality of locations near the m / z value derived from the target component. A mass analysis control unit that switches the control conditions of the mass analysis unit, an ion intensity detection unit that detects the ion intensity of the m / z value selected by the mass analysis control unit, and a change in the ion intensity detected by the ion intensity detection unit. Based on the neighborhood pattern calculation unit that calculates the ratio of the ion intensity of the m / z values at the plurality of locations in a state where the component amount is constant, and the neighborhood pattern calculation unit that calculates the neighborhood pattern. It is provided with an ion intensity correction unit for compensating for the insufficient ion intensity measured, and measures the change over time in the ion intensity of the m / z value derived from the target component.

また、試料成分をイオン化し、イオン化されたイオンを選別し、選別したイオンの強度を測定する質量分析方法において、目的成分に由来するm/z値の近傍の複数個所のm/z値のイオン強度を連続して検出する近傍探索ステップと、検出したイオン強度の変化を計算する成分量の変化計算ステップと、成分量一定の状態での前記目的成分に由来するm/z値の近傍の複数個所のm/z値のイオン強度の比率を求める近傍パターン計算ステップと、前記近傍パターンに基づき測定したイオン強度の不足を補填するイオン強度補正ステップと、を備え、目的成分に由来するm/z値のイオン強度の経時変化を測定する。 Further, in a mass analysis method in which a sample component is ionized, ionized ions are selected, and the intensity of the selected ions is measured, ions having m / z values at a plurality of locations near the m / z value derived from the target component are used. A neighborhood search step that continuously detects the intensity, a component amount change calculation step that calculates the change in the detected ion intensity, and a plurality of vicinitys of the m / z value derived from the target component in a state where the component amount is constant. It is provided with a neighborhood pattern calculation step for obtaining the ratio of the ion intensity of the m / z value at each location and an ion intensity correction step for compensating for the shortage of the ion intensity measured based on the neighborhood pattern, and the m / z derived from the target component. Measure the change in ionic intensity of the value over time.

制御条件の変化によるm/z値のずれを検知し、定量分析に必要な感度を維持しつつ、m/z値のずれに対する証憑を与え、m/z値のずれにより減少したイオン量を補うことができる質量分析装置および質量分析方法を実現することができる。 Detects m / z value deviations due to changes in control conditions, provides vouchers for m / z value deviations while maintaining the sensitivity required for quantitative analysis, and compensates for the amount of ions reduced due to m / z value deviations. It is possible to realize a mass spectrometer and a mass spectrometry method that can be used.

質量分析装置の概略構成図である。It is a schematic block diagram of a mass spectrometer. 質量分析部の校正の説明図である。It is explanatory drawing of the calibration of the mass spectrometry part. 四重極マスフィルターの制御の説明図である。It is explanatory drawing of the control of a quadrupole mass filter. 制御部の動作フローチャートである。It is an operation flowchart of a control unit. 制御部の概略ブロック図である。It is a schematic block diagram of a control part. 四重極マスフィルターの校正時のマススペクトルを例示する図である。It is a figure which illustrates the mass spectrum at the time of calibration of a quadrupole mass filter. 成分量の変化計算の説明図である。It is explanatory drawing of the change calculation of the component amount. m/z値と近傍パターンの変化を説明する図である。It is a figure explaining the change of the m / z value and the neighborhood pattern. m/z値のずれとスコアの計算例を示す図である。It is a figure which shows the deviation of m / z value and the calculation example of a score. スコアとm/z値のずれの関係を示す図である。It is a figure which shows the relationship between the score and the deviation of m / z value. 条件設定画面の例を示す図である。It is a figure which shows the example of the condition setting screen. 補正の結果を示すログ出力の例を示す図である。It is a figure which shows the example of the log output which shows the result of correction. スキャン測定による近傍探索の説明図である。It is explanatory drawing of the neighborhood search by scan measurement. スキャン測定による成分量の変化計算の説明図である。It is explanatory drawing of the change calculation of the component amount by a scan measurement. スキャン測定によるm/z値と近傍パターンの変化を説明する図である。It is a figure explaining the change of the m / z value and the neighborhood pattern by a scan measurement. 接線による成分量の変化計算の説明図である。It is explanatory drawing of the change calculation of the component amount by a tangent line.

本発明を実施するための形態について、図面を参照して説明する。なお、本発明の実施の形態は、ここに示す実施例に限定されるものではない。 A mode for carrying out the present invention will be described with reference to the drawings. The embodiments of the present invention are not limited to the examples shown here.

(実施例1)
実施例1として、目的成分に由来するイオンのm/z値の近傍探索において、その両端に着目する質量分析装置及び質量分析方法を開示する。
(Example 1)
As Example 1, a mass spectrometer and a mass spectrometry method that focus on both ends of the m / z value of an ion derived from a target component in the neighborhood search will be disclosed.

図1は実施例1の質量分析装置の概略構成図である。 FIG. 1 is a schematic configuration diagram of the mass spectrometer of the first embodiment.

図1において、試料供給装置110は質量分析装置120に試料成分を供給する。例えば、液体クロマトグラフやガスクロマトグラフ、試料を加熱することで成分を気化させる装置などが試料供給装置である。また、イメージング質量分析などにおいて試料を保持する機構も、試料供給装置の一種である。 In FIG. 1, the sample supply device 110 supplies sample components to the mass spectrometer 120. For example, a liquid chromatograph, a gas chromatograph, a device for vaporizing components by heating a sample, and the like are sample supply devices. A mechanism for holding a sample in imaging mass spectrometry is also a kind of sample supply device.

質量分析装置120は、イオン化部121と、質量分析部122と、検出部123と、制御部124とを備える。試料供給装置110から供給される試料成分を、イオン化部121によりイオン化し、質量分析部122で、イオン化部121から供給されたイオンのうち所定のm/z値のイオンを選別し、検出部123においてイオン強度として測定する。 The mass spectrometer 120 includes an ionization unit 121, a mass spectrometry unit 122, a detection unit 123, and a control unit 124. The sample component supplied from the sample supply device 110 is ionized by the ionization unit 121, and the mass spectrometry unit 122 selects ions having a predetermined m / z value from the ions supplied from the ionization unit 121, and the detection unit 123. It is measured as an ion intensity in.

制御部124は、入出力装置130からの指示に従い、所定のm/z値を選別するよう質量分析部122を制御し、m/z値に対応するイオン強度のデータを入出力装置130に返す。 The control unit 124 controls the mass spectrometry unit 122 to select a predetermined m / z value according to the instruction from the input / output device 130, and returns the ionic strength data corresponding to the m / z value to the input / output device 130. ..

ここで、m/z値を所定の範囲でスキャンしつつ、m/z値とイオン強度の並び、すなわちマススペクトルを得る手法をスキャン測定と呼ぶ。また、質量分析部122で選別するm/z値を固定し、イオン強度の経時変化、すなわちクロマトグラムを得るのがSIMやSRMと呼ばれる定量測定である。 Here, a method of obtaining an arrangement of m / z value and ionic strength, that is, a mass spectrum while scanning the m / z value in a predetermined range is called scan measurement. Further, it is a quantitative measurement called SIM or SRM to fix the m / z value selected by the mass spectrometer 122 and obtain a change with time of ionic strength, that is, a chromatogram.

図2は質量分析部122の校正を説明する図である。 FIG. 2 is a diagram illustrating calibration of the mass spectrometer 122.

一般に、質量分析部122の校正では、m/z値が既知の試料成分を測定し、測定値と既知のm/z値が一致するよう質量分析部122の制御条件を調整する。図2において、m/z値がMの試料成分についてMを与えるピークの範囲(例えばM±2Da)をスキャンして得たピークが校正前のピーク210であり、そこから得られるm/z値をM’とした。 Generally, in the calibration of the mass spectrometer 122, a sample component having a known m / z value is measured, and the control conditions of the mass spectrometer 122 are adjusted so that the measured value and the known m / z value match. In FIG. 2, the peak obtained by scanning the range of peaks giving M to the sample component having an m / z value of M (for example, M ± 2Da) is the peak 210 before calibration, and the m / z value obtained from the peak 210. Was set to M'.

ここで、m/z値がMとなるように、質量分析部122の制御値を調整して得られたピークが校正後のピーク220である。また、両端に矢印を記載した直線230が示す幅は、校正後のピーク220の半値幅、すなわちピーク高さの半値におけるピークの横幅である。 Here, the peak obtained by adjusting the control value of the mass spectrometer 122 so that the m / z value becomes M is the peak 220 after calibration. Further, the width indicated by the straight line 230 with arrows at both ends is the half-value width of the peak 220 after calibration, that is, the width of the peak at the half-value of the peak height.

なお、校正に際し、m/z値が既知の試料成分の供給量およびイオン化部121における効率は一定である必要がある。イオン化される試料成分の量が変化すると、検出部123で得られるイオン強度が変化し、図2におけるM’やMの値を正しく測定することができない。 At the time of calibration, it is necessary that the supply amount of the sample component having a known m / z value and the efficiency in the ionization unit 121 are constant. When the amount of the sample component to be ionized changes, the ionic strength obtained by the detection unit 123 changes, and the values of M'and M in FIG. 2 cannot be measured correctly.

次に、質量分析部122として四重極マスフィルターを採用した場合を想定し、制御の概要を説明する。 Next, assuming a case where a quadrupole mass filter is adopted as the mass spectrometer 122, an outline of control will be described.

四重極マスフィルターの場合、それは一対の対向電極で構成される。これらに印加する電圧Φ(ファイ)は次式(1)で与えられる。 In the case of a quadrupole mass filter, it consists of a pair of counter electrodes. The voltage Φ (phi) applied to these is given by the following equation (1).

Figure 0007077481000001
Figure 0007077481000001

ただし、式(1)において、Φは四重極への印加電圧、Uは直流電圧(DC)、Vは交流電圧の最大値(RF)、cosは三角関数、ωは角周波数、tは時間である。 However, in equation (1), Φ is the voltage applied to the quadrupole, U is the DC voltage (DC), V is the maximum value of the AC voltage (RF), cos is a trigonometric function, ω is the angular frequency, and t is the time. Is.

式(1)において、校正で調整する制御値は、直流電圧Uと交流電圧の最大値Vである。 In the formula (1), the control values adjusted by calibration are the DC voltage U and the maximum value V of the AC voltage.

次に、図3を用いて四重極マスフィルターの制御について説明する。図3は校正点のm/z値であるM、M、M(M<M<M)と、校正で求めたRFとDCの関係を示すグラフである。Next, the control of the quadrupole mass filter will be described with reference to FIG. FIG. 3 is a graph showing the relationship between M 1 , M 2 , and M 3 (M 1 <M 2 <M 3 ), which are the m / z values of the calibration points, and RF and DC obtained by calibration.

図3において、質量分析装置120の制御部124は、各校正点を結んだ線分(校正曲線)310を参照し、校正点M、M、Mから内挿または外挿により、制御すべきm/z値に対応する点を求め、RFとDCの値を決定する。一般に定量の目的成分のm/z値は既知であるため、定量成分そのものを利用して校正してもよい。In FIG. 3, the control unit 124 of the mass spectrometer 120 refers to the line segment (calibration curve) 310 connecting the calibration points, and controls the calibration points M 1 , M 2 , and M 3 by interpolation or extrapolation. Find the point corresponding to the m / z value to be, and determine the RF and DC values. Since the m / z value of the target component for quantification is generally known, the quantification component itself may be used for calibration.

なお、装置の校正とは、RFやDCなどの値を調整することで、目的とするピークのm/z値や、半値幅などで代表されるピーク形状を、目的とする値に合わせる作業となる。 In addition, calibration of the device is the work of adjusting the m / z value of the target peak or the peak shape represented by the full width at half maximum to the target value by adjusting the values of RF and DC. Become.

図4は制御部124の動作フローチャートであり、図5は制御部124の概略ブロック図である。 FIG. 4 is an operation flowchart of the control unit 124, and FIG. 5 is a schematic block diagram of the control unit 124.

以下、図4の動作フローチャートに従い、図5などを参照しつつ制御動作の詳細を説明する。 Hereinafter, the details of the control operation will be described with reference to FIG. 5 and the like according to the operation flowchart of FIG.

(1)初期設定ステップ410
初期設定として、定量分析の条件や、近傍探索の条件などを入出力装置130において設定し、制御部124に転送する。図5の制御部124における各処理ブロックは、初期設定された値を用いて処理を進める。
(1) Initial setting step 410
As initial settings, the conditions for quantitative analysis, the conditions for neighborhood search, and the like are set in the input / output device 130 and transferred to the control unit 124. Each processing block in the control unit 124 of FIG. 5 proceeds with processing using the initially set value.

初期設定の詳細は、近傍探索ステップ420以降の説明に含める。 The details of the initial setting are included in the description after the neighborhood search step 420.

また、入力装置130の表示部における条件設定画面の例を図11に示す。詳細は動作フローチャートの説明の後で説明する。 Further, FIG. 11 shows an example of a condition setting screen in the display unit of the input device 130. Details will be described after the explanation of the operation flowchart.

(2)近傍探索ステップ420
近傍探索は従来のSIMやSRMなどの定量測定に替わるものである。その処理は、図5のマスフィルター制御部(質量分析制御部)510とイオン強度検出部520で実施する。
(2) Neighborhood search step 420
The neighborhood search replaces the conventional quantitative measurement such as SIM and SRM. The processing is carried out by the mass filter control unit (mass spectrometry control unit) 510 and the ion intensity detection unit 520 of FIG.

マスフィルター制御部510は、目的成分に由来するm/z値の近傍の探索点に対応する2または複数個所のm/z値を選択するため、マスフィルターである質量分析装部122の制御条件を順次切り替える。また、マスフィルター制御部(質量分析制御部)510は、近傍探索の探索点を、目的成分に由来するm/z値-0.5から+0.5の範囲においてスキャン測定することも可能である。 Since the mass filter control unit 510 selects two or a plurality of m / z values corresponding to the search points in the vicinity of the m / z value derived from the target component, the control conditions of the mass spectrometric unit 122, which is a mass filter. Are switched sequentially. Further, the mass filter control unit (mass spectrometry control unit) 510 can also scan and measure the search point of the neighborhood search in the range of the m / z value of −0.5 to +0.5 derived from the target component. ..

イオン強度検出部520は、マスフィルター制御部510で選択したm/z値のイオンを検出部123から検出する。 The ion intensity detection unit 520 detects ions having an m / z value selected by the mass filter control unit 510 from the detection unit 123.

まず、近傍探索における探索点と基本的な考え方について図6を用いて説明する。 First, the search points and the basic idea in the neighborhood search will be described with reference to FIG.

近傍探索は、従来のSIMやSRMが測定対象とする目的成分のm/z値に替わり、その近傍のm/z値に対するイオン強度の経時変化を求める処理である。ここでは、目的成分のm/z値も含む3点を探索点とする例を示す。なお、実施例1は、複数の探索点の両端に着目する方法であり、必ずしもm/z値に対応する探索点は必要ではない。 The neighborhood search is a process for obtaining a change over time in ionic strength with respect to the m / z value in the vicinity thereof, instead of the m / z value of the target component to be measured by the conventional SIM or SRM. Here, an example is shown in which three points including the m / z value of the target component are set as search points. The first embodiment is a method of paying attention to both ends of a plurality of search points, and the search points corresponding to the m / z values are not always required.

図6は、四重極マスフィルターの校正時のマススペクトル610を例示する図である。これは、目的成分のm/z値Mの前後をスキャン測定して得られるものであり、図6ではそのピークトップの近傍に着目して示している。 FIG. 6 is a diagram illustrating a mass spectrum 610 at the time of calibration of a quadrupole mass filter. This is obtained by scanning and measuring before and after the m / z value M of the target component, and is shown in FIG. 6 focusing on the vicinity of the peak top.

マススペクトル610から、ピークトップ630の左側の点620に対応するM-Δ、ピークトップ630のM、およびピークトップ630の右側の点640に対応するM+Δの3か所のm/z値を探索点として設定する。 From the mass spectrum 610, search for three m / z values of M-Δ corresponding to the point 620 on the left side of the peak top 630, M of the peak top 630, and M + Δ corresponding to the point 640 on the right side of the peak top 630. Set as a point.

すなわち、近傍探索として、M-Δ、M、M+Δの3点に対応するイオン強度の経時変化を求める。ここで、Δ(例えば0.1Da)と探索点は、マススペクトルのピークトップの近傍において、ピークの幅などを参考にして初期設定する。校正時のピークの幅が広ければΔを大きく、狭ければ小さい値を設定するなど、最適なΔの値はピーク形状に依存する。なお、目的成分に由来するm/z値(探索点)の近傍のm/z値のうち、目的成分に由来するm/z値の両端に位置するm/z値は、理想的にはピークトップから数%低下した位置であり、最大でもピークの半値以上となるm/z値が望ましい。 That is, as a neighborhood search, the change with time of the ionic strength corresponding to the three points of M−Δ, M, and M + Δ is obtained. Here, Δ (for example, 0.1 Da) and the search point are initially set in the vicinity of the peak top of the mass spectrum with reference to the peak width and the like. The optimum value of Δ depends on the peak shape, such as setting a large value for Δ if the width of the peak at the time of calibration is wide and a small value if it is narrow. Of the m / z values in the vicinity of the m / z value (search point) derived from the target component, the m / z values located at both ends of the m / z value derived from the target component are ideally peaks. It is desirable that the m / z value is a few percent lower than the top and is at least half the peak value.

例えば、前者の例は、半値幅0.7のピークにおいて、Mのイオン強度を100%とし、m/z値(M-Δ)のイオン強度が98.6%となるΔ値0.05Daを採用する。また後者の例は、半値幅1.0のピークにおいて、m/z値M-Δのイオン強度が50%となるΔ値0.5Daを採用する。 For example, in the former example, at the peak with a half width of 0.7, the ionic strength of M is 100%, and the ionic strength of m / z value (M−Δ) is 98.6%. adopt. Further, in the latter example, a Δ value of 0.5 Da at which the ionic strength of the m / z value M−Δ is 50% at the peak of the half width of 1.0 is adopted.

なお、近傍探索では、初期設定であらかじめ指定された探索点の一式(この例ではM-Δ、M、M+Δ)のイオン強度を順に求める工程を1サイクルとし、それを繰り返す。 In the neighborhood search, the process of sequentially obtaining the ionic strength of a set of search points (M−Δ, M, M + Δ in this example) specified in advance in the initial setting is set as one cycle, and the process is repeated.

校正時は基本的に成分量が一定であるが、実試料の測定では成分量が変化する状況で測定する。 At the time of calibration, the amount of components is basically constant, but when measuring an actual sample, the amount of components changes.

そこで、1サイクルの各点の測定の間に成分量が増加している例を、図6の点660(P)、点670(P)、および点680(P)で示した。それぞれの点は対応するサイクルを明確にするために線分690で結んでいる。ここで、Pは測定点を示し、jはサイクル内での探索点の番号かつサイクル内の測定順である。Therefore, an example in which the amount of the component increases during the measurement of each point in one cycle is shown at points 660 (P 1 ), 670 (P 2 ), and 680 (P 3 ) in FIG. Each point is connected by a line segment 690 to clarify the corresponding cycle. Here, P j indicates a measurement point, and j is the number of the search point in the cycle and the measurement order in the cycle.

参考として、成分量が増加している状況でスキャン測定した場合に得られるマススペクトルを点線650で示した(実際にはスキャン測定は行わない)。成分量が右肩上がりであることを受け、マススペクトル650はマススペクトル610に比べて右側が高くなっており、Pのイオン強度がPに対して大きくなっていることがわかる。すなわち、成分量が一定でない場合、装置のm/z値のずれを正しく検出することができない。As a reference, the mass spectrum obtained by scanning measurement in a situation where the amount of components is increasing is shown by a dotted line 650 (actually, scanning measurement is not performed). It can be seen that the mass spectrum 650 is higher on the right side than the mass spectrum 610 due to the rising component amount, and the ionic strength of P 3 is higher than that of P 1 . That is, if the amount of components is not constant, the deviation of the m / z value of the device cannot be detected correctly.

ここで、マススペクトル610を得た際の成分量をQ0、ピークトップ630のイオン強度をIBとする。また、マススペクトル610における各探索点jについて、ピークトップとのイオン強度比をRとする。Here, the amount of components when the mass spectrum 610 is obtained is Q0, and the ionic strength of the peak top 630 is IB. Further, for each search point j in the mass spectrum 610, the ionic strength ratio with the peak top is R j .

例えば、探索点1のイオン強度比(点620のイオン強度/IB)がR、探索点2のイオン強度比(点630のイオン強度/IB)がR(=1)、そして探索点3のイオン強度比(点640のイオン強度/IB)がRとする。各探索点Pにおいて、成分量をQ、イオン強度比をIとすると、イオン強度は成分量に比例する前提で次式(2)が成り立つ。For example, the ionic strength ratio of search point 1 (ionic strength of point 620 / IB) is R 1 , the ionic strength ratio of search point 2 (ionic strength of point 630 / IB) is R 2 (= 1), and the search point 3 Ionic strength ratio ( ion strength at point 640 / IB) is R3. Assuming that the component amount is Q j and the ionic strength ratio is I j at each search point P j , the following equation (2) holds on the premise that the ionic strength is proportional to the component amount.

Figure 0007077481000002
Figure 0007077481000002

ただし、式(2)において、Q0は校正時の成分量、IBは校正時のピークトップのイオン強度、下付添字jは探索点の番号、Rは探索点jのイオン強度比(成分量一定でのイオン強度のピークトップ比)、Iは探索点jのイオン強度、Qは探索点jの成分量である。However, in the formula (2), Q0 is the component amount at the time of calibration, IB is the ionic strength of the peak top at the time of calibration, the subscript j is the number of the search point, and R j is the ion intensity ratio (component amount) of the search point j. The peak top ratio of the ionic strength at a constant level), I j is the ionic strength of the search point j, and Q j is the component amount of the search point j.

ここで、Rは校正時に求めることができ、初期設定において後の処理で参照できるように登録する。装置のm/z値がずれた場合、そのずれに応じてRの値も変化する。その場合、(I/Q)とRが比例することから、探索点数が十分かつピークトップを含む必要があるが、(I/Q)をその最大値で規格化することでRを再計算することができる。Here, R j can be obtained at the time of calibration, and is registered so that it can be referred to in a later process in the initial setting. When the m / z value of the device deviates, the value of Rj also changes according to the deviation. In that case, since (I j / Q j ) and R j are proportional, it is necessary that the number of search points is sufficient and the peak top is included . R j can be recalculated.

(3)成分量の変化計算ステップ430
近傍探索のサイクルで得られる各探索点のイオン強度から、成分量の変化を計算する処理であるステップ430について説明する。これは図5の成分量の変化計算部530において実施する。成分量の変化計算部530は、イオン強度検出部520が検出したイオンの強度変化を計算する。その詳細を、図7を参照しつつ説明する。成分量の変化計算部530は、イオン強度の変化を複数のサイクルにより計算し、前後のサイクルから2次式等で成分量の変化を計算する。
(3) Change calculation step 430 of component amount
Step 430, which is a process of calculating the change in the amount of components from the ionic strength of each search point obtained in the neighborhood search cycle, will be described. This is carried out in the change calculation unit 530 of the component amount in FIG. The component amount change calculation unit 530 calculates the ion intensity change detected by the ion intensity detection unit 520. The details will be described with reference to FIG. 7. The component amount change calculation unit 530 calculates the change in ionic strength by a plurality of cycles, and calculates the change in the component amount by a quadratic expression or the like from the cycles before and after.

図7は、成分量の変化計算の説明図であり、1サイクル3点の近傍探索を連続して実施した様子を示す図である。図7の横軸は時間、縦軸はイオン強度であり、1サイクル内の測定点を直線で結んだ。また、点線710、720、730は、それぞれ探索点1、2、3の経時変化を示している。また、各サイクルに対し測定順に1から番号を与え、C~Cと表記した。FIG. 7 is an explanatory diagram of the change calculation of the component amount, and is a diagram showing a state in which the neighborhood search of three points in one cycle is continuously performed. The horizontal axis of FIG. 7 is time and the vertical axis is ionic strength, and the measurement points in one cycle are connected by a straight line. Further, the dotted lines 710, 720, and 730 indicate changes with time of the search points 1, 2, and 3, respectively. In addition, numbers from 1 were given to each cycle in the order of measurement, and they were described as C 1 to C 5 .

なお、本発明の実施例1においては、装置のm/z値のずれを検出するものであり、図7に示した例では説明のためピークトップをPから、P側にずらしている。In the first embodiment of the present invention, the deviation of the m / z value of the apparatus is detected, and in the example shown in FIG. 7, the peak top is shifted from P 2 to the P 1 side for explanation. ..

前述のイオン強度比Rは、成分量が一定の条件で求めているため、各探索点すなわちm/z値におけるイオンの測定効率を示している。装置のm/z値がずれている場合、ずれ量に従いRの値も変化するが、ずれているピークトップのm/z値に近い探索点の方がRの値も大きくなる。一連のサイクルにおいて、Rの値が維持されると考えると、ある時間領域において最も多くのイオンを計測した探索点が、最もピークトップに近いといえる。Since the above-mentioned ion intensity ratio Rj is obtained under a constant condition for the amount of components, it indicates the measurement efficiency of ions at each search point, that is, the m / z value. When the m / z value of the device is deviated, the value of R j also changes according to the amount of deviation, but the value of R j is larger at the search point closer to the m / z value of the deviated peak top. Considering that the value of Rj is maintained in a series of cycles, it can be said that the search point that measured the largest number of ions in a certain time domain is the closest to the peak top.

例えば、図7では点線710、720、730が与える面積を考えれば、P、P、Pの順となるのは自明であり、ここではPが最もピークトップに近い。For example, in FIG. 7, considering the area given by the dotted lines 710, 720, and 730, it is obvious that the order is P 2 , P 1 , and P 3 , where P 2 is the closest to the peak top.

ここで、Rはそのm/z値におけるイオンの測定効率であれば、各探索点において同じ時間範囲における面積Aの比とRの比は一致する。そこに着目し、面積AからRを求めるのが実施例1の方法である。Here, if R j is the measurement efficiency of ions at the m / z value, the ratio of the area A j and the ratio of R j in the same time range at each search point are the same. Focusing on this, the method of Example 1 is to obtain R j from the area A j .

面積を求めるには、前述したように、個々の探索点のイオン強度では判断できない。ここでは、サイクルごとの時間範囲において、イオン強度の変化を示す関数を求め、その積分値から面積を得る。 As mentioned above, the area cannot be determined by the ionic strength of each search point. Here, the function showing the change in ionic strength is obtained in the time range for each cycle, and the area is obtained from the integrated value.

例えば、k番目のサイクル(C)に着目すると、探索点jのイオン量の変化は前後のサイクルの各点にも着目し、(Tj,k-1,j,k-1)、(Tj,k,j,k)、(Tj,k+1,j,k+1)の3点から2次式を求めることができる。ここで、Tj,kは探索した時間、Ij,kはイオン強度である。For example, focusing on the k-th cycle (C k ), the change in the amount of ions at the search point j also focuses on each point in the previous and next cycles (T j, k-1, I j, k-1 ). A quadratic equation can be obtained from the three points (T j, k, I j, k ) and (T j, k + 1, I j, k + 1 ). Here, T j and k are the searched time, and I j and k are the ionic strength.

このようにして探索点jごとに求めた2次式について、Cの時間区間での積分値を得る。For the quadratic equation obtained for each search point j in this way, the integral value in the time interval of C k is obtained.

次に、k番目のサイクルにおける探索点jの変化を示す2次式(3)の例を示す。 Next, an example of the quadratic equation (3) showing the change of the search point j in the kth cycle is shown.

Figure 0007077481000003
Figure 0007077481000003

ただし、式(3)において、下付添字jは探索点の番号、下付添字kはサイクルの番号、Tj,kは探索点jサイクルkにおける時間、Ij,kは探索点jサイクルkにおけるイオン強度、tは時間、f(j,k,t)は探索点j、サイクルk、時間tにおけるイオン強度を示す関数である。However, in the equation (3), the subscript j is the search point number, the subscript k is the cycle number, T j and k are the time at the search point j cycle k, and I j and k are the search point j cycle k. Ionic strength, t is time, and f (j, k, t) is a function indicating the ionic strength at the search point j, cycle k, and time t.

ここではラグランジュ補間を採用し2次式を例に説明したが、1次式の利用も考えられる。すなわち、前後の点を直線(1次式)で結び、イオン強度の変化を示す関数f(j,k,t)とする。また、その他の補間法を採用してもよい。 Here, Lagrange interpolation is adopted and a quadratic equation is used as an example, but the use of a linear equation is also conceivable. That is, the points before and after are connected by a straight line (linear equation) to obtain a function f (j, k, t) indicating the change in ionic strength. Further, other interpolation methods may be adopted.

関数f(j,k,t)について、サイクルkの時間範囲における面積をAとする。このAを例えば式(4)で求める。For the function f (j, k, t), let A j be the area in the time range of cycle k. This Aj is obtained by, for example, the equation (4).

Figure 0007077481000004
Figure 0007077481000004

ただし、式(4)において、tCkmaxはサイクルkの時間範囲の最大値、tCkminはサイクルkの時間範囲の最小値、Aは(サイクルkの)探索点jに対応する面積である。However, in the equation (4), t Ckmax is the maximum value in the time range of the cycle k, t Ckmin is the minimum value in the time range of the cycle k, and A j is the area corresponding to the search point j (of the cycle k).

このようにして求めた面積Aの比は、Rの比に一致する。例えば、探索点が3点の場合、A:A:A=R:R:Rとなる。また、装置のm/z値がずれている場合、Rはそれぞれの探索点におけるイオン強度比を示す。The ratio of the area A j thus obtained corresponds to the ratio of R j . For example, when there are three search points, A 1 : A 2 : A 3 = R 1 : R 2 : R 3 . Further, when the m / z value of the apparatus deviates, R j indicates the ionic strength ratio at each search point.

なお、イオン強度の値が小さい場合は、面積の比は不安定であり、イオン強度が大きくなると真の値に近づく。そこで、複数のサイクルにおけるAの値を順次Sに加算し、後の近傍パターン計算450において利用する。Sは面積Aの積算値である。When the ionic strength value is small, the area ratio is unstable, and when the ionic strength increases, it approaches the true value. Therefore, the values of Aj in a plurality of cycles are sequentially added to Sj and used in the subsequent neighborhood pattern calculation 450. S j is an integrated value of the area A j .

(4)近傍パターン確定判定ステップ440
ステップ440では、成分量の変化計算ステップ430で求めた面積Aの積算値Sが安定したかどうかの判定を行う。
(4) Neighborhood pattern confirmation determination step 440
In step 440, it is determined whether or not the integrated value Sj of the area Aj obtained in the change calculation step 430 of the component amount is stable.

例えば、すべての積算値Sの合計値が近傍パターン確定判定のしきい値を超えた場合、近傍パターンが安定したとみなしてYと判定し、イオン強度補正ステップ460に進む。近傍パターン確定判定ステップ440にて、近傍パターン確定判定のしきい値を越えていない場合は、近傍パターンは安定していないとみなしてNと判定し、近傍パターン計算ステップ450に進む。For example, when the total value of all the integrated values S j exceeds the threshold value for determining the neighborhood pattern, it is considered that the neighborhood pattern is stable and it is determined as Y, and the process proceeds to the ion intensity correction step 460. If the threshold value for the neighborhood pattern determination determination is not exceeded in the neighborhood pattern determination determination step 440, the neighborhood pattern is regarded as unstable and determined to be N, and the process proceeds to the neighborhood pattern calculation step 450.

別法として、積算値Sに替わって、イオン強度Iを積算した値を利用することもできる。さらに、当該サイクルのSと一つ前のサイクルで得たSについて、探索点ごとに差の絶対値を求め、その合計がしきい値以下になることで判定する方法なども考えられる。Alternatively, instead of the integrated value S j , the integrated value of the ionic strength I j can be used. Further, a method of obtaining the absolute value of the difference between S j of the cycle and S j obtained in the previous cycle for each search point and determining the total value is equal to or less than the threshold value can be considered.

(5)近傍パターン計算ステップ450
近傍パターン計算ステップ450は、図5における近傍パターン計算部540において実施する。
(5) Neighborhood pattern calculation step 450
The neighborhood pattern calculation step 450 is performed by the neighborhood pattern calculation unit 540 in FIG.

近傍パターン計算ステップ450では、近傍パターン計算部540が、成分量の変化計算部530において求めた探索点jにおける面積Aの積算値Sに基づき、その比率である近傍パターンを計算する。あるいは、近傍パターン計算部540は、成分量一定の状態での近傍のm/z値(目的成分に由来するm/z値の近傍の複数個所のm/z値)のイオン強度の比率を求め、近傍パターンを計算する。また、近傍パターン計算部540は、2次式で計算された成分量の変化から、単位成分量あたりのイオン強度を近傍パターンとすることもできる。さらに、近傍パターンからスコアを求め、対応するRやm/z値のずれを求める。In the neighborhood pattern calculation step 450, the neighborhood pattern calculation unit 540 calculates the neighborhood pattern, which is the ratio thereof, based on the integrated value Sj of the area Aj at the search point j obtained by the change calculation unit 530 of the component amount. Alternatively, the neighborhood pattern calculation unit 540 obtains the ratio of the ionic strength of the m / z value in the vicinity (m / z values at a plurality of locations near the m / z value derived from the target component) in a state where the component amount is constant. , Calculate the neighborhood pattern. Further, the neighborhood pattern calculation unit 540 can also set the ionic strength per unit component amount as the neighborhood pattern from the change in the component amount calculated by the quadratic equation. Further, the score is obtained from the neighborhood pattern, and the deviation of the corresponding Rj and m / z values is obtained.

簡単のため、近傍探索が2点の場合について、図8を用いて説明する。 For the sake of simplicity, the case where the neighborhood search has two points will be described with reference to FIG.

図8は、四重極マスフィルターのm/z値のずれとイオン強度の関係を示すグラフであり、m/z値と近傍パターンの変化を説明する図である。 FIG. 8 is a graph showing the relationship between the deviation of the m / z value of the quadrupole mass filter and the ionic strength, and is a diagram for explaining the change between the m / z value and the neighborhood pattern.

図8において、校正直後のm/z値Mのマススペクトル810、m/z値が0.05ずれたマススペクトル820、m/z値が0.10ずれたマススペクトル830に対し、近傍探索で得られる2点(点840と点850)をそれぞれプロットし直線で結んだ。このように、m/z値のずれが大きくなるに従い、2点間のイオン強度の差も大きくなる。 In FIG. 8, the mass spectrum 810 with the m / z value M immediately after calibration, the mass spectrum 820 with the m / z value deviated by 0.05, and the mass spectrum 830 with the m / z value deviated by 0.10 are searched in the vicinity. The two obtained points (point 840 and point 850) were plotted and connected by a straight line. In this way, as the deviation of the m / z value increases, the difference in ionic strength between the two points also increases.

すなわち、ピークトップ付近のマススペクトルの形状が安定している場合、近傍パターンからm/z値のずれを求めることができる。特に、近傍パターンの両端に着目し、その比から求めたスコアとm/z値のずれを対応づけることができる。 That is, when the shape of the mass spectrum near the peak top is stable, the deviation of the m / z value can be obtained from the neighborhood pattern. In particular, it is possible to focus on both ends of the neighborhood pattern and associate the score obtained from the ratio with the deviation of the m / z value.

また、マススペクトルの形状は四重極マスフィルターの校正時に確定する。そのため、近傍パターンの差とm/z値のずれの関係を校正時に計算することが可能となる。 The shape of the mass spectrum is determined when the quadrupole mass filter is calibrated. Therefore, it is possible to calculate the relationship between the difference in the neighborhood pattern and the deviation of the m / z value at the time of calibration.

図9は、m/z値のずれとスコアの計算例を示す図であり、マススペクトルのピークが半値幅0.7の場合について、m/z値のずれ(m/z Error)910に対するスコアの計算方法について例示する表である。 FIG. 9 is a diagram showing a deviation of the m / z value and a calculation example of the score. When the peak of the mass spectrum has a half width of 0.7, the score with respect to the deviation of the m / z value (m / z Error) 910 is shown. It is a table exemplifying the calculation method of.

図9において、イオン強度(I:Ion Intensity)920において、上段1はm/z値がM-0.05の探索点1であり、下段2はm/z値がM+0.05の探索点2であり、それぞれのm/z値のずれに対するイオン強度の値を示している。なお、説明のためイオン強度の最大値を100として規格化した値を記載している。 In FIG. 9, in the ionic strength (I: Ion Integrity) 920, the upper row 1 is the search point 1 having an m / z value of M-0.05, and the lower row 2 is the search point 2 having an m / z value of M + 0.05. It shows the value of the ionic strength with respect to the deviation of each m / z value. For the sake of explanation, the value standardized with the maximum value of ionic strength as 100 is described.

ここで、m/z値のずれが0.00、0.05、0.10の場合は、図8におけるマススペクトル810、820、830に対応する。このような数値は、校正時に実測したマススペクトルから求めることができる。もしくは半値幅から求めたガウシアン関数を適用することで、m/z値とイオン強度の関係を決定することもできる。 Here, when the deviation of the m / z value is 0.00, 0.05, 0.10, it corresponds to the mass spectra 810, 820, and 830 in FIG. Such a numerical value can be obtained from the mass spectrum actually measured at the time of calibration. Alternatively, the relationship between the m / z value and the ionic strength can be determined by applying the Gaussian function obtained from the full width at half maximum.

さらに、図8に示すように、m/z値がずれたマススペクトルを仮定することで、m/z値のずれと近傍探索で得られる各点のイオン強度の比率などを求めることができる。 Further, as shown in FIG. 8, by assuming a mass spectrum in which the m / z value is deviated, the ratio of the m / z value deviation and the ionic strength of each point obtained by the neighborhood search can be obtained.

図9において、イオン強度比(R:Ion Intensity Ratio)930は、イオン強度の値を最大強度(この場合は100)で規格化した数値である。また、イオン強度の測定予想比(R’)940は両端の探索点のイオン強度に着目し、大きい方を1として規格化し、R’max-R’の値を差(Dif)とした。この例は2点であるため、R’のmaxは2となる。さらに差(Dif)を1000倍しスコア(Score)950としている。In FIG. 9, the ionic strength ratio (R: Ion Integrity Ratio) 930 is a numerical value obtained by normalizing the ionic strength value at the maximum intensity (100 in this case). Further, the expected measurement ratio (R') 940 of the ionic strength was standardized with the larger one as 1 and the value of R'max -R'1 as the difference (Dif), paying attention to the ionic strength of the search points at both ends. Since this example has two points, the max of R'is 2. Further, the difference (Dif) is multiplied by 1000 to obtain a score (Score) of 950.

図9においてR’は実測値から算出できる比率であり、そこから求めたスコアから、m/z値のずれや、Rを決定できることを示している。In FIG. 9, R'is a ratio that can be calculated from the measured value, and it is shown that the deviation of the m / z value and R j can be determined from the score obtained from the ratio.

このようにして求めた、スコアとm/z値のずれの関係を図10に示す。図10の横軸は図9におけるスコア(Score)であり、縦軸はm/z値のずれ(m/z Error)である。また、線分1010はそれぞれのプロットを内挿、または外挿したものである。このようにして求めて作成したグラフを示す図10は、半値幅0.7におけるスコアとm/z値のずれの関係を示している。 FIG. 10 shows the relationship between the score and the deviation of the m / z value obtained in this way. The horizontal axis of FIG. 10 is the score (Score) in FIG. 9, and the vertical axis is the deviation of the m / z value (m / z Error). In addition, the line segment 1010 is an interpolated or extrapolated plot of each. FIG. 10 showing the graph thus obtained and created shows the relationship between the score and the deviation of the m / z value in the half width of 0.7.

ここで整理すると、初期設定において、図9に示すスコアとm/z値のずれ、およびRの関係を登録しておく。さらに、近傍パターン計算で求めた両端のイオン強度比(R’とR’max)からスコアを計算することで、初期設定で登録したRやm/z値のずれを決定し、後述するイオン強度補正ステップ460やm/z値補正ステップ470で参照する。また、Sは積算前の初期値に戻す。To summarize here, in the initial setting, the deviation between the score and the m / z value shown in FIG. 9, and the relationship between Rj are registered. Furthermore, by calculating the score from the ionic strength ratios ( R'1 and R'max ) at both ends obtained by the neighborhood pattern calculation, the deviation of the Rj and m / z values registered in the initial settings is determined, which will be described later. It is referred to in the ion intensity correction step 460 and the m / z value correction step 470. Further, S j is returned to the initial value before integration.

この例では、高m/z値側にずれている場合を挙げたが、当然低m/z値側につても計算が可能である。また、近傍探索するピークトップ部分が左右対称であれば、片側のみの計算で済む。スコアの符号が逆になることで、どちら側にずれているか判断できる。 In this example, the case where the value is shifted to the high m / z value side is mentioned, but of course, the calculation can be performed even on the low m / z value side. Further, if the peak top portion for neighborhood search is symmetrical, only one side can be calculated. By reversing the sign of the score, it is possible to determine which side is deviated.

(6)イオン強度補正ステップ460
続いて、イオン強度補正ステップ460について説明する。これは図5におけるイオン強度補正部550での処理に相当する。イオン強度補正部550は、近傍パターン計算部540で計算された近傍パターンに基づき測定したイオン強度の不足を補填する。
(6) Ionic strength correction step 460
Subsequently, the ionic strength correction step 460 will be described. This corresponds to the processing in the ion intensity correction unit 550 in FIG. The ionic strength correction unit 550 compensates for the lack of ionic strength measured based on the neighborhood pattern calculated by the neighborhood pattern calculation unit 540.

ここでは、イオン強度比Rの値に基づいて測定したイオン強度を補正する方法を説明する。Here, a method of correcting the ionic strength measured based on the value of the ionic strength ratio Rj will be described.

探索点jに対する本来のイオン強度は、I/Rで求めることができる。すなわち、サイクル内のイオン強度の合計は全ての探索点の合計、Σ(I/R)として表すことができる。この値は、一般的な定量測定におけるSIMやSRMのクロマトグラム1点のイオン強度に相当する。The original ionic strength with respect to the search point j can be obtained by I j / R j . That is, the total ionic strength in the cycle can be expressed as the total of all search points, Σ (I j / R j ). This value corresponds to the ionic strength of one point on the chromatogram of SIM or SRM in a general quantitative measurement.

なお、各Rの値は、初期設定ステップ410で初期化値が設定され、近傍パターン計算ステップ450で更新される。さらに、m/z値補正が適用された場合は、初期値に戻される。The value of each R j is set as an initialization value in the initial setting step 410, and is updated in the neighborhood pattern calculation step 450. Further, when the m / z value correction is applied, it is returned to the initial value.

(7)m/z値補正ステップ470
m/z値補正ステップ470は、図5におけるm/z値補正部560において実施する。
(7) m / z value correction step 470
The m / z value correction step 470 is carried out by the m / z value correction unit 560 in FIG.

m/z値補正ステップ470では、近傍パターン計算ステップ450で求めたm/z値のずれに基づき、m/z値を補正する。つまり、近傍パターン計算ステップ450で求めた近傍パターンに基づき質量分析制御部510の制御条件を修正する。近傍パターンが計算されていない場合は、m/z値補正を実施しない。 In the m / z value correction step 470, the m / z value is corrected based on the deviation of the m / z value obtained in the neighborhood pattern calculation step 450. That is, the control conditions of the mass spectrometry control unit 510 are modified based on the neighborhood pattern obtained in the neighborhood pattern calculation step 450. If the neighborhood pattern has not been calculated, the m / z value correction is not performed.

具体的には、次回の測定から、図3に示したように校正で得たm/z値と図3に示すRFおよびDCの関係に従い、ピークトップに対応するであろうずれたm/z値に相当するRFとDCの値を入出力装置130から出力する。また、m/z値の補正は四重極マスフィルターの質量精度などの仕様に応じ、しきい値を設けることも考えられる。 Specifically, from the next measurement, according to the relationship between the m / z value obtained by calibration as shown in FIG. 3 and RF and DC shown in FIG. 3, the shifted m / z that will correspond to the peak top. The RF and DC values corresponding to the values are output from the input / output device 130. Further, for the correction of the m / z value, it is conceivable to set a threshold value according to the specifications such as the mass accuracy of the quadrupole mass filter.

m/z値のずれを補正した場合、イオン強度補正ステップ460(イオン強度補正部550による処理)で参照するイオン強度比Rの値をm/z値のずれがない校正直後の数値に戻す。When the deviation of the m / z value is corrected, the value of the ion intensity ratio Rj referred to in the ion intensity correction step 460 (processing by the ion intensity correction unit 550) is returned to the value immediately after calibration without the deviation of the m / z value. ..

(8)終了判定ステップ480
図4の動作フローチャートには終了判定ステップ480とあるが、本発明の実施例1は目的成分の定量測定においてm/z値のずれに対応するものであり、その効果はピーク形状が大きく変わらなければ、一連の定量測定が続く限り継続できる。新たに四重極マスフィルターの校正がなされ、初期設定ステップ410が実施された場合、実質的な終了となる。
(8) End determination step 480
Although the operation flowchart of FIG. 4 has the end determination step 480, the first embodiment of the present invention corresponds to the deviation of the m / z value in the quantitative measurement of the target component, and the effect is that the peak shape does not change significantly. For example, it can be continued as long as a series of quantitative measurements continues. When the quadrupole mass filter is newly calibrated and the initial setting step 410 is performed, the process is substantially completed.

以上、本発明における実施例1を、動作フローチャートに従い説明した。なお、この実施例1において、初期設定ステップ410の際に用いる条件設定画面(入出力装置130の表示画面)について図11を用いて説明する。図11は条件設定画面の例を示す図である。図11に示した例においては、説明のため定量測定の目的成分を一つに限定している。実際の条件設定画面では、複数の目的成分ごとに条件を設定するなどの応用が考えられる。 The first embodiment of the present invention has been described above according to the operation flowchart. In the first embodiment, the condition setting screen (display screen of the input / output device 130) used in the initial setting step 410 will be described with reference to FIG. FIG. 11 is a diagram showing an example of a condition setting screen. In the example shown in FIG. 11, the target component for quantitative measurement is limited to one for the sake of explanation. On the actual condition setting screen, applications such as setting conditions for each of a plurality of target components can be considered.

図11において、条件設定画面は目的成分のマススペクトルを規定したTarget1110、近傍探索の条件を規定したSearch1120、m/z値補正の条件を規定したCorrection1130に大別される。それぞれにおけるパラメータの意味を詳述する。 In FIG. 11, the condition setting screen is roughly classified into Target 1110 that defines the mass spectrum of the target component, Search 1120 that defines the conditions for neighborhood search, and Direction 1130 that defines the conditions for m / z value correction. The meaning of each parameter will be described in detail.

Target1110は、定量の目的成分の定義を行う表示領域を示し、Top1111は目的成分のピークトップのm/z値である。また、Width1112は目的成分のマススペクトルのピークにおける半値幅であり、Dwell Time1113は目的成分のイオン強度の取り込み時間であり、定量処理を行う際のクロマトグラムの1点の値を与える。 Target 1110 indicates a display area for defining the target component for quantification, and Top 1111 is the m / z value of the peak top of the target component. Further, Wide 1112 is the half width at the peak of the mass spectrum of the target component, and Dwell Time 1113 is the uptake time of the ionic strength of the target component, which gives a value of one point in the chromatogram when performing quantitative processing.

Search1120は、近傍探索条件を設定する領域を示し、Point Number1121は、近傍探索の探索点の点数であり、Dwell Time/Point Numberは近傍探索の1点の取り込み時間である。また、Delta1122は探索点のm/z値の間隔である。また、Threshold1123は近傍パターン確定判定のしきい値であり、実施例1ではすべてのSの合計値がこの値を超えた場合、安定したとみなす。The Search 1120 indicates a region for setting the neighborhood search condition, the Point Number 1121 is the score of the search point of the neighborhood search, and the Dwell Time / Point Number is the capture time of one point of the neighborhood search. Further, Delta1122 is an interval of m / z values of search points. Further, Threshold 1123 is a threshold value for determining the neighborhood pattern, and in Example 1, when the total value of all S js exceeds this value, it is considered to be stable.

Correction1130は、m/z値の補正条件を設定する領域を示す。また、Step1131はm/z値補正の刻みを示し、m/z値のずれがこの値以上になった場合に、補正を実施する。また、Maximum1132はm/z値補正の最大値であり、校正直後のm/z値に対し、これ以上の補正はエラー扱いとするしきい値である。 Direction 1130 indicates an area for setting a correction condition for the m / z value. Further, Step 1131 indicates a step of m / z value correction, and when the deviation of the m / z value becomes greater than or equal to this value, the correction is performed. Further, Maximum1132 is the maximum value of m / z value correction, and is a threshold value for which any further correction is treated as an error with respect to the m / z value immediately after calibration.

図11に示したような条件設定された、WidthやDeltaなどの値などから、図9に示すScore950とm/z Error910、およびイオン強度比Rの関係を求める。ここで用いるマススペクトルは、四重極マスフィルターの校正で得たものを利用できる。マススペクトルのピーク形状をガウス関数(Gaussian function)で代用する場合、半値幅Wのピークに対応する標準偏差σ、およびその際のガウス関数の式は次のような式(5)、(6)に表すことができる。The relationship between the Score 950, the m / z Error 910, and the ionic strength ratio Rj shown in FIG. 9 is obtained from the values such as Wide and Delta set under the conditions shown in FIG. As the mass spectrum used here, the one obtained by calibration of the quadrupole mass filter can be used. When the peak shape of the mass spectrum is substituted by the Gaussian function, the standard deviation σ corresponding to the peak of the half width W and the Gaussian function at that time are the following equations (5) and (6). Can be expressed in.

Figure 0007077481000005
Figure 0007077481000005

Figure 0007077481000006
Figure 0007077481000006

ただし、式(5)、式(6)において、σは標準偏差、Wは半値幅、lnは自然対数、xは探索点のm/z値、yはイオン強度、Mは目的成分のm/z値、Eはm/z値のずれ、πは円周率、eはネイピア数である。 However, in equations (5) and (6), σ is the standard deviation, W is the half-value width, ln is the natural logarithm, x is the m / z value of the search point, y is the ion intensity, and M is the target component m /. The z value and E are the deviations of the m / z value, π is the pi, and e is the Napier number.

上記式(5)、(6)により、イオン強度yを求め、最大値を100として規格化すれば、図9に示したイオン強度(I:Ion Intensity)920に対応する値等を得ることができる。 If the ionic strength y is obtained from the above equations (5) and (6) and standardized with the maximum value set to 100, a value corresponding to the ionic strength (I: Ion Integrity) 920 shown in FIG. 9 can be obtained. can.

図4に示す初期設定ステップ410では、図11に示した画面により設定された入力を受け、例えば、図6のマススペクトル610と探索点の表示や、図8のようなシミュレーション結果を表示することにより、より効果的にパラメータを設定することができる。 In the initial setting step 410 shown in FIG. 4, the input set by the screen shown in FIG. 11 is received, and for example, the mass spectrum 610 and the search point in FIG. 6 are displayed, and the simulation result as shown in FIG. 8 is displayed. Therefore, the parameters can be set more effectively.

また、図12は、補正の状況を示すログ出力の例を示す図である。 Further, FIG. 12 is a diagram showing an example of log output showing the status of correction.

図12において、1行目1210はログのタイトル、2行目1220は目的成分のm/z値を示す。3行目1230はログの項目を示し、左から順にログ日時、スコア、m/z値、m/z値のずれを意味する。4行目1240以下は、ログの実体を示す。 In FIG. 12, the first line 1210 shows the title of the log, and the second line 1220 shows the m / z value of the target component. The third line 1230 indicates the log item, and means the deviation of the log date / time, the score, the m / z value, and the m / z value in order from the left. The fourth line below 1240 indicates the substance of the log.

図12に示した例では、目的成分のm/z値693.10に対し、2018/09/01のずれていない状況から、2018/09/13に図11においてStepで規定したm/z値0.05のずれに相当するスコア55を超えたため、693.15から目的とするm/z値を693.10に補正している。2018/09/18のほぼ誤差なしの状況から2018/09/25には逆方向の低m/z値側にずれたため補正している。補正を実施した場合は、スコアの横にアスタリスクを付け識別している。 In the example shown in FIG. 12, the m / z value specified by Step in FIG. 11 is set on 2018/09/13 from the situation where the m / z value 693.10 of the target component is not deviated from 2018/09/01. Since the score 55 corresponding to the deviation of 0.05 was exceeded, the target m / z value was corrected to 693.10 from 693.15. From the situation of 2018/09/18 with almost no error, it is corrected because it deviates to the low m / z value side in the opposite direction on 2018/09/25. When corrections are made, an asterisk is added next to the score to identify them.

本発明の実施例1は上記のように構成されているので、制御条件の変化によるm/z値のずれを検知し、定量分析に必要な感度を維持しつつ、m/z値のずれに対する証憑を与え、m/z値のずれにより減少したイオン量を補うことができる質量分析装置および質量分析方法を実現することができる。 Since the first embodiment of the present invention is configured as described above, it detects the deviation of the m / z value due to the change of the control condition, maintains the sensitivity required for quantitative analysis, and with respect to the deviation of the m / z value. It is possible to realize a mass spectrometer and a mass spectrometry method that can provide a voucher and compensate for the amount of ions reduced due to the deviation of the m / z value.

さらに、実施例1によれば、m/z値のずれを自動的に補正することができる質量分析装置および質量分析方法を実現することができる。 Further, according to the first embodiment, it is possible to realize a mass spectrometer and a mass spectrometry method that can automatically correct the deviation of the m / z value.

(実施例2)
次に、本発明の実施例2について説明する。
(Example 2)
Next, Example 2 of the present invention will be described.

実施例2は、近傍探索においてピークトップ近傍をスキャン測定する例である。なお、スキャン測定の両端に着目すれば実施例1が適用できる。ここでは、スキャン測定において、ピークトップに対応する探索点が含まれているものとして、スキャン測定で得たマススペクトルから、Rやm/z値のずれを求める方法を説明する。The second embodiment is an example of scanning and measuring the vicinity of the peak top in the neighborhood search. The first embodiment can be applied by paying attention to both ends of the scan measurement. Here, it is assumed that the search point corresponding to the peak top is included in the scan measurement, and a method of obtaining the deviation of the Rj and m / z values from the mass spectrum obtained by the scan measurement will be described.

実施例2における質量分析装置の概略構成や動作フローチャートは実施例1と同様である。ここでは、実施例1と実施例2との違いを中心に説明する。 The schematic configuration and operation flowchart of the mass spectrometer in the second embodiment are the same as those in the first embodiment. Here, the difference between the first embodiment and the second embodiment will be mainly described.

(1)初期設定ステップ410
実施例1と同様の初期設定を実施する。ただし、実施例2では、図11に示す条件設定画面においてm/z値の間隔Delta(Δ)を、探索に用いるスキャン測定の始点と終点のm/z値の差と読み替える。
(1) Initial setting step 410
The same initial settings as in Example 1 are carried out. However, in the second embodiment, in the condition setting screen shown in FIG. 11, the m / z value interval Delta (Δ) is read as the difference between the m / z values of the start point and the end point of the scan measurement used for the search.

(2)近傍探索ステップ420
近傍探索の探索点についての考え方を図13に示す。図13はスキャン測定による近傍探索の説明図である。
(2) Neighborhood search step 420
FIG. 13 shows the concept of the search point of the neighborhood search. FIG. 13 is an explanatory diagram of neighborhood search by scan measurement.

図13において、校正時の成分量が一定の条件で取得できるマススペクトル1310にて、目的成分のm/z値M-Δ/2からM+Δ/2までスキャン測定すると、理想的には点1320から点1330に至るマススペクトル1340を得ることができる。このスキャン範囲M-Δ/2からM+Δ/2が探索の範囲となる。 In FIG. 13, when the mass spectrum 1310 in which the amount of components at the time of calibration can be obtained by scanning is measured from the m / z value M−Δ / 2 to M + Δ / 2 of the target component, ideally from the point 1320. A mass spectrum of 1340 up to point 1330 can be obtained. The search range is from this scan range M−Δ / 2 to M + Δ / 2.

また、探索点の点数は初期設定で指定されたものであり、その取得時間は同じく指定された取り込み時間(Dwell Time)を点数で割ったものとなる。探索点の点数(Point Number)が2または3のように少ない場合は実施例1と同等とみなすことができるが、ここではさらに大きい値を設定することを想定する。例えば、Δを0.2、点数200、Dwell Time100msとした場合、1点あたり500μs、各探索点のm/z値の刻みは0.001 (0.2/200)Daとなる。 Further, the score of the search point is the one specified in the initial setting, and the acquisition time is the same as the designated capture time (Dwell Time) divided by the score. When the number of search points (Point Number) is as small as 2 or 3, it can be regarded as equivalent to that of the first embodiment, but here, it is assumed that a larger value is set. For example, when Δ is 0.2, the score is 200, and the Dwell Time is 100 ms, 500 μs per point, and the m / z value step of each search point is 0.001 (0.2 / 200) Da.

なお、1回のスキャン測定を1サイクルとし、それを繰り返す。また、実施例1と同様に、成分量が増加する場合に得られるマススペクトル1350において、近傍探索では点1360から点1370に至るマススペクトル1380を得る。ここで、Pが最初の探索点であり、Pmaxが最後の探索点となる。上の例ではmaxは200となる。ここで探索点jのイオン強度をIとする。In addition, one scan measurement is regarded as one cycle, and it is repeated. Further, as in the first embodiment, in the mass spectrum 1350 obtained when the amount of components increases, the mass spectrum 1380 from the point 1360 to the point 1370 is obtained in the neighborhood search. Here, P 1 is the first search point and P max is the last search point. In the above example, max is 200. Here, let I j be the ionic strength of the search point j.

なお、実施例2の場合は1点あたりの測定時間が短いため、マススペクトル1340や1380のイオン強度はばらつきの大きいものとなる。そこで、マススペクトル1340や1380を、回帰分析等で求めた2次関数などの近似式に置き換えてもよい。 In the case of Example 2, since the measurement time per point is short, the ionic strengths of the mass spectra 1340 and 1380 have large variations. Therefore, the mass spectra 1340 and 1380 may be replaced with an approximate expression such as a quadratic function obtained by regression analysis or the like.

実施例2では実施例1と同様、イオン強度比Rを算出することができる。マススペクトル1340のピークトップのイオン強度をIB、それを得た成分量をC0とし、マススペクトル1340の各探索点jについて、ピークトップとのイオン強度比をRとすると、実施例1の式(2)を求めることができる。すなわち、実施例1の基本的な考え方は実施例2においても適用することができる。In Example 2, the ionic strength ratio Rj can be calculated as in Example 1. Assuming that the ionic strength of the peak top of the mass spectrum 1340 is IB, the amount of the obtained component is C0, and the ionic strength ratio with the peak top is Rj for each search point j of the mass spectrum 1340, the formula of Example 1 is used. (2) can be obtained. That is, the basic idea of the first embodiment can be applied to the second embodiment as well.

(3)成分量の変化計算ステップ430
実施例2における成分量の変化計算ステップ430について、図14を参照しつつ説明する。
(3) Change calculation step 430 of component amount
The change calculation step 430 of the component amount in Example 2 will be described with reference to FIG.

図14は、スキャン測定による成分量の変化計算の説明図である。 FIG. 14 is an explanatory diagram of calculation of change in the amount of components by scan measurement.

図14において、各サイクルのPからPmaxまでが近傍探索の結果となる。また、点線1410は探索点1(P)の経時変化を示し、点線1420は最後の探索点(Pmax)の経時変化を示している。In FIG. 14, P 1 to P max of each cycle are the results of the neighborhood search. Further, the dotted line 1410 shows the change with time of the search point 1 (P 1 ), and the dotted line 1420 shows the change with time of the last search point (P max ).

実施例1では、探索点の両端であるPとPmaxの経時変化を2次式等で補間し、それらの面積AとAmaxを、順次SとSmaxに加算し、後の近傍パターン計算に利用している。In the first embodiment, the temporal changes of P 1 and P max at both ends of the search point are interpolated by a quadratic equation or the like, and the areas A 1 and A max are sequentially added to S 1 and S max , and later. It is used for neighborhood pattern calculation.

実施例2では、面積値が最大となる点線1430に着目した方法を提供する。ここで、面積値が最大となるとは、イオン強度比Rの値が最大(=1)になる点、すなわち成分量が一定の場合のビークトップに相当するm/z値に対応する点を結んだものに他ならない。すなわち、点線1430は成分量の変化を示している。また、一連のサイクルにおけるスキャン測定の結果と点線1430の接点は、測定した時点での目的成分のm/z値に対応するため、m/z値のずれを検出することができる。また、成分量Qの変化とマススペクトルのイオン強度IからRを求めることで、イオン強度補正も可能となる。In the second embodiment, a method focusing on the dotted line 1430 having the maximum area value is provided. Here, the maximum area value means that the value of the ionic strength ratio Rj becomes the maximum ( = 1), that is, the point corresponding to the m / z value corresponding to the beak top when the component amount is constant. It is nothing but a tied thing. That is, the dotted line 1430 shows the change in the amount of components. Further, since the contact point between the scan measurement result and the dotted line 1430 in a series of cycles corresponds to the m / z value of the target component at the time of measurement, the deviation of the m / z value can be detected. Further, by obtaining R j from the change in the component amount Q j and the ionic strength I j of the mass spectrum, the ionic strength can be corrected.

図15は、参考としてスキャン測定によるm/z値と近傍パターンの変化を例示した図である。図8に示した例と図15に示した例とを比較すると、図15に示した例の方が、より高密度にm/z値とイオン強度の関係をとらえていることがわかる。なお、スキャン測定の場合、1点あたりの測定時間が短くなるため、イオン強度のばらつきが大きくなる特徴がある。 FIG. 15 is a diagram illustrating changes in the m / z value and the neighborhood pattern by scan measurement for reference. Comparing the example shown in FIG. 8 with the example shown in FIG. 15, it can be seen that the example shown in FIG. 15 captures the relationship between the m / z value and the ionic strength at a higher density. In the case of scan measurement, since the measurement time per point is shortened, there is a feature that the variation in ionic strength becomes large.

次に、接線によって点線1430に相当する成分量の変化を求める方法について図16を参照しつつ説明する。 Next, a method of obtaining a change in the amount of the component corresponding to the dotted line 1430 by the tangent line will be described with reference to FIG.

図16は、接線による成分量の変化計算の説明図である。成分量の変化計算部530がスキャン測定によるマススペクトルに外接する接線により成分量の変化を求める。 FIG. 16 is an explanatory diagram of the calculation of the change in the amount of components due to the tangent line. The change calculation unit 530 of the component amount obtains the change of the component amount by the tangent line circumscribing the mass spectrum by the scan measurement.

図16はサイクルkにおける成分量の変化を、二本の接線1610と1620で表している。すなわち、サイクルk-1(Ck-1)とサイクルk(C)のマススペクトル(スキャン測定によるマススペクトル)に外接する接線1610と、サイクルkとサイクルk+1(Ck+1)のマススペクトルに外接する接線1620を引き、それらの交点1630を得る。FIG. 16 shows the change in the amount of components in cycle k by two tangents 1610 and 1620. That is, the tangent line 1610 tangent to the mass spectrum of cycle k-1 (C k-1 ) and cycle k (C k ) (mass spectrum measured by scan measurement), and the tangent to the mass spectrum of cycle k and cycle k + 1 (C k + 1 ). Draw tangents 1620 to obtain their intersections 1630.

このようにして、サイクルが進むに連れ、スペクトルの接戦の交点と交点を結ぶ線分を決定することができる。それを、図14の点線1430に相当するもの、すなわち成分量Qの変化とする。In this way, as the cycle progresses, the line segment connecting the intersections of the close battles of the spectrum can be determined. It corresponds to the dotted line 1430 in FIG. 14, that is, the change in the component amount Qj .

ここで、マススペクトルのイオン強度I/Qが、成分量一定の条件におけるイオン強度となり、ピークトップを1としたイオン強度比Rを求めることができる。Here, the ionic strength I j / Q j of the mass spectrum becomes the ionic strength under the condition that the amount of components is constant, and the ionic strength ratio R j with the peak top as 1 can be obtained.

なお、この処理において求めたRの値は、不安定であることが考えられる。そこで、実施例1のAに替わり、各サイクルで求めたI/Qの値にそのサイクルでの交点のイオン強度を掛けた値を重みとして掛けた数値をAとする。複数のサイクルにおけるAの値を順次Sに加算し、後の近傍パターン計算で利用する。It is considered that the value of Rj obtained in this process is unstable. Therefore, instead of A j in Example 1, the value obtained by multiplying the value of I j / Q j obtained in each cycle by the ionic strength at the intersection in that cycle is taken as the weight. The value of A j in a plurality of cycles is sequentially added to S j and used in the later neighborhood pattern calculation.

ここでは接線を求める方法を示したが、例えばサイクルkにおいて、サイクルk-1、サイクルk、サイクルk+1のマススペクトル(スキャン測定によるマススペクトル)に外接する2次関数などを採用してもよい。サイクルk-1、サイクルk、サイクルk+1に対応するマススペクトルの最大値を与える点、3か所を最初の候補として2次関数を仮定する。それが同じサイクル内のマススペクトルのピーク内に入り、別の場所から出るようなら、それらの点の中間にある点を新たな接点候補とし、最終的に3サイクルのマススペクトルと接する2次式を求める。 Although the method of obtaining the tangent line is shown here, for example, in cycle k, a quadratic function circumscribing the mass spectrum (mass spectrum measured by scan) of cycle k-1, cycle k, and cycle k + 1 may be adopted. A quadratic function is assumed with three points as the first candidates that give the maximum value of the mass spectrum corresponding to the cycle k-1, the cycle k, and the cycle k + 1. If it enters within the peak of the mass spectrum in the same cycle and exits from another location, a point in the middle of those points is used as a new contact candidate, and finally a quadratic equation that contacts the mass spectrum of 3 cycles. Ask for.

(4)近傍パターン確定判定ステップ440
近傍パターン確定判定ステップ440では、実施例2において成分量の変化計算ステップ430で求めた、積算値Sが安定したかどうかを判定する。
(4) Neighborhood pattern confirmation determination step 440
In the neighborhood pattern determination determination step 440, it is determined whether or not the integrated value Sj obtained in the component amount change calculation step 430 in Example 2 is stable.

例えば、すべての積算値Sの合計値が近傍パターン確定判定ステップ440のしきい値を超えた場合、近傍パターンが安定しているとみなしてY(Yes)、そうでない場合はN(No)と判定する。For example, if the total value of all the integrated values S j exceeds the threshold value of the neighborhood pattern determination determination step 440, it is considered that the neighborhood pattern is stable, Y (Yes), and if not, N (No). Is determined.

別法として、積算値Sに替わって、イオン強度Iを積算した値を利用することもできる。さらに、当該サイクルのSと一つ前のサイクルで得たSについて、探索点ごとに差の絶対値を求め、その合計がしきい値以下になることで判定する方法なども考えられる。Alternatively, instead of the integrated value S j , the integrated value of the ionic strength I j can be used. Further, a method of obtaining the absolute value of the difference between S j of the cycle and S j obtained in the previous cycle for each search point and determining the total value is equal to or less than the threshold value can be considered.

(5)近傍パターン計算ステップ450
近傍パターン計算ステップ450では、成分量の変化計算で求めたSに対し、それらの最大値を1として格化した値を求めてRとする。さらに、最大値に相当するjの値に対応するm/z値が、現在のm/z値であり、目的成分のm/z値とのずれを求めることができる。また、Rは、イオン強度補正ステップ460で利用でき、m/z値のずれはm/z値補正ステップ470で利用する。
(5) Neighborhood pattern calculation step 450
In the neighborhood pattern calculation step 450, for S j obtained by the calculation of the change in the amount of components, the value rated with the maximum value of 1 is obtained and used as R j . Further, the m / z value corresponding to the value of j corresponding to the maximum value is the current m / z value, and the deviation from the m / z value of the target component can be obtained. Further, R j can be used in the ionic strength correction step 460, and the deviation of the m / z value is used in the m / z value correction step 470.

ここで、Sの値はクリアする。Here, the value of S j is cleared.

(6)イオン強度補正ステップ460
イオン強度補正ステップ460は、実施例1と同様にして行う。なお、補正後に用いる各点のイオン強度Iは2次式などで近似する前の値を採用することもできる。
(6) Ionic strength correction step 460
The ionic strength correction step 460 is performed in the same manner as in Example 1. The ionic strength Ij of each point used after the correction may be a value before being approximated by a quadratic equation or the like.

(7)m/z値補正ステップ470
m/z値補正ステップ470は、実施例1と同様である。m/z値を補正した場合、イオン強度補正ステップ460で参照するRを、m/z値がずれていない校正直後の数値に戻す。
(7) m / z value correction step 470
The m / z value correction step 470 is the same as that of the first embodiment. When the m / z value is corrected, Rj referred to in the ion intensity correction step 460 is returned to the value immediately after calibration in which the m / z value is not deviated.

(8)終了判定ステップ480
終了判定ステップ480は実施例1と同様である。
(8) End determination step 480
The end determination step 480 is the same as that of the first embodiment.

以上、目的成分のm/z値Mの近傍探索にスキャン測定を採用した例を実施例2として示した。 As described above, the example in which the scan measurement is adopted for the neighborhood search of the m / z value M of the target component is shown as Example 2.

実施例2においても、実施例1と同様な効果を得ることができる。 In Example 2, the same effect as in Example 1 can be obtained.

(実施例3)
次に、本発明の実施例3について説明する。
(Example 3)
Next, Example 3 of the present invention will be described.

実施例3は、近傍探索においてピークトップ近傍をスキャン測定する実施例2の別法である。 Example 3 is another method of Example 2 in which the vicinity of the peak top is scanned and measured in the neighborhood search.

(1)初期設定ステップ410
実施例2と同様である。
(1) Initial setting step 410
It is the same as Example 2.

(2)近傍探索ステップ420
実施例2と同様である。
(2) Neighborhood search step 420
It is the same as Example 2.

(3)成分量の変化計算ステップ430
実施例1で示したように、各サイクルの面積比はRの比に一致する。そこで、実施例3においても、実施例1と同様にして、各探索点jについて面積Aと、複数のサイクルにおけるAの積算値Sを求める。ここで、Sを最大にするjに対応するm/z値が、その時点でのビークトップに対応する。また、SはRに対応する。
(3) Change calculation step 430 of component amount
As shown in Example 1, the area ratio of each cycle corresponds to the ratio of Rj . Therefore, also in the third embodiment, the area A j and the integrated value S j of A j in a plurality of cycles are obtained for each search point j in the same manner as in the first embodiment. Here, the m / z value corresponding to j that maximizes S j corresponds to the beak top at that time. Further, S j corresponds to R j .

(4)近傍パターン確定判定ステップ440
近傍パターン確定判定ステップ440では、成分量の変化計算ステップ430で求めた面積Aの積算値Sが安定したかどうかの判定を行う。
(4) Neighborhood pattern confirmation determination step 440
In the neighborhood pattern determination determination step 440, it is determined whether or not the integrated value Sj of the area Aj obtained in the change calculation step 430 of the component amount is stable.

例えば、すべての積算値Sの合計値が近傍パターン確定判定ステップ440のしきい値を超えた場合、近傍パターンが安定したとみなしてY(Yes)、そうでない場合はN(No)と判定する。For example, if the total value of all the integrated values S j exceeds the threshold value in the neighborhood pattern determination determination step 440, it is determined that the neighborhood pattern is stable and Y (Yes), otherwise it is determined as N (No). do.

別法として、積算値Sに替わって、イオン強度Iを積算した値を利用することもできる。さらに、当該サイクルのSと一つ前のサイクルで得たSについて、探索点ごとに差の絶対値を求め、その合計がしきい値以下になることで判定する方法なども考えられる。Alternatively, instead of the integrated value S j , the integrated value of the ionic strength I j can be used. Further, a method of obtaining the absolute value of the difference between S j of the cycle and S j obtained in the previous cycle for each search point and determining the total value is equal to or less than the threshold value can be considered.

(5)近傍パターン計算ステップ450
近傍パターン計算ステップ450では、成分量の変化計算で求めたSに対し、それらの最大値を1として格化した値を求めてRとする。さらに、最大値に相当するjの値に対応するm/z値が、現在のm/z値であり、目的成分のm/z値とのずれを求めることができる。また、Rは、イオン強度補正ステップ460で利用でき、m/z値のずれはm/z値補正ステップ470で利用する。
(5) Neighborhood pattern calculation step 450
In the neighborhood pattern calculation step 450, for S j obtained by the calculation of the change in the amount of components, the value rated with the maximum value of 1 is obtained and used as R j . Further, the m / z value corresponding to the value of j corresponding to the maximum value is the current m / z value, and the deviation from the m / z value of the target component can be obtained. Further, R j can be used in the ionic strength correction step 460, and the deviation of the m / z value is used in the m / z value correction step 470.

ここで、Sの値はクリアする。Here, the value of S j is cleared.

(6)イオン強度補正ステップ460
イオン強度補正ステップ460は、実施例2と同様にして行う。
(6) Ionic strength correction step 460
The ionic strength correction step 460 is performed in the same manner as in Example 2.

(7)m/z値補正ステップ470
m/z値補正ステップ470は、実施例2と同様にして行う。
(7) m / z value correction step 470
The m / z value correction step 470 is performed in the same manner as in the second embodiment.

(8)終了判定ステップ480
終了判定ステップ480は、実施例2と同様にして行う。
(8) End determination step 480
The end determination step 480 is performed in the same manner as in the second embodiment.

実施例3においても、実施例1と同様な効果を得ることができる。 In Example 3, the same effect as in Example 1 can be obtained.

上述した実施例1、2、3は、それぞれ目的成分に由来するイオンの検出に関するものである。もし目的成分が量的に少ない場合、m/z値のずれが目的成分と同様である他の成分に対して本発明を適用し、目的成分のm/z値も併せて補正する応用が考えられる。 The above-mentioned Examples 1, 2 and 3 relate to the detection of ions derived from the target component, respectively. If the amount of the target component is small, the present invention may be applied to other components having the same deviation of the m / z value as the target component, and the m / z value of the target component may also be corrected. Will be.

例えば、目的成分を構成する幾つかの炭素を安定同位体である炭素13(6個の陽子と7個の中性子で構成)に置き換えたものを内部標準物質として、定量測定において目的成分と同時に測定する。その際、目的成分に由来するイオンについては従来どおりSIMやSRMにより測定し、内部標準物質のイオンに対し本発明を適用し、イオン強度補正やm/z値の補正を実施する。その際、目的成分に由来するイオンのm/z値のずれも内部標準物質と同じとみなし、目的成分に対応するイオンのm/z値のずれに対応するイオン強度補正やm/z値補正を適用する。 For example, carbon 13 (composed of 6 protons and 7 neutrons), which is a stable isotope, is used as an internal standard to replace some of the carbons that make up the target component, and it is measured at the same time as the target component in quantitative measurement. do. At that time, the ions derived from the target component are measured by SIM or SRM as before, and the present invention is applied to the ions of the internal standard substance to correct the ionic strength and the m / z value. At that time, the deviation of the m / z value of the ion derived from the target component is regarded as the same as that of the internal standard substance, and the ionic strength correction and the m / z value correction corresponding to the deviation of the m / z value of the ion corresponding to the target component are considered. To apply.

上述した例は、内部標準物質を用いる例であるが、必ずしも内部標準物質を用いなくともよい。つまり、内部標準物質を用いることなく、イオン強度検出部520は、目的成分のm/z値と同じずれを示すイオンの近傍を探索し、近傍パターン計算部により、探索したイオンのイオン強度の近傍パターンを求め、その結果に基づいて、目的成分のイオン強度をイオン強度補正部550により補正し、m/z値補正部560によりm/z値を補正するように構成することができる。 The above-mentioned example is an example in which an internal standard substance is used, but it is not always necessary to use an internal standard substance. That is, without using the internal standard substance, the ionic strength detection unit 520 searches for the vicinity of the ion showing the same deviation as the m / z value of the target component, and the neighborhood pattern calculation unit searches for the vicinity of the ionic strength of the searched ion. A pattern can be obtained, and based on the result, the ionic strength of the target component can be corrected by the ionic strength correction unit 550, and the m / z value can be corrected by the m / z value correction unit 560.

上述した本発明の実施例1~3により、SIMやSRMによる定量測定において次の効果を得ることができる。 According to Examples 1 to 3 of the present invention described above, the following effects can be obtained in quantitative measurement by SIM or SRM.

(1)m/z値のずれにより減少したイオン強度の補填
本来の定量測定では、目的成分のピークトップに対応するm/z値のイオン強度の経時変化を求めることが理想である。本発明によれば、温度・湿度等の影響により装置のm/z値がずれた場合でも、測定しているm/z値近傍の探索点のパターンから、ずれている状況を検出できる。さらに、測定したイオン強度を、本来測定できるはずのイオン強度に変換することが可能となる。m/z値のずれに起因する試料成分の単位量あたりのイオン強度の変動を防ぐことは、定量できる濃度範囲や精度の維持において重要である。
(1) Compensation for the ionic strength reduced due to the deviation of the m / z value In the original quantitative measurement, it is ideal to obtain the time course of the ionic strength of the m / z value corresponding to the peak top of the target component. According to the present invention, even if the m / z value of the device is deviated due to the influence of temperature, humidity, etc., the deviated situation can be detected from the pattern of the search points near the m / z value being measured. Further, the measured ionic strength can be converted into the ionic strength that should be able to be measured originally. It is important to prevent the fluctuation of the ionic strength per unit amount of the sample component due to the deviation of the m / z value in maintaining the quantifiable concentration range and accuracy.

(2)m/z値がずれていない証憑の取得
SIMやSRMによる定量測定において、定量値とともにm/z値がずれていない証憑を直接得ることができる。
(2) Acquisition of voucher that does not deviate from m / z value In quantitative measurement by SIM or SRM, a voucher that does not deviate from m / z value can be directly obtained together with the quantitative value.

(3)m/z値のずれの自動補正
目的成分の近傍パターンから、m/z値のずれを検出することができる。すなわち、リアルタイムにm/z値のずれ補正を可能とする。装置のm/z値が徐々にずれるような場合でも、より最適な定量性能を維持することが可能となる。
(3) Automatic correction of m / z value deviation The m / z value deviation can be detected from the neighborhood pattern of the target component. That is, it is possible to correct the deviation of the m / z value in real time. Even when the m / z value of the device gradually deviates, it is possible to maintain more optimum quantitative performance.

副次的な効果として、例えば定量測定中にスキャン測定を挿入するような処理が不要となるとこから、より目的成分に由来するイオンの検出に装置の時間を割くことが可能となり、定量の感度向上等にも寄与する。 As a side effect, for example, since the process of inserting a scan measurement during quantitative measurement becomes unnecessary, it becomes possible to devote more time to the device to detect ions derived from the target component, and the sensitivity of quantitative measurement. It also contributes to improvement.

なお、上述した例においては、m/z値のずれを自動補正するように構成したが、m/z値のずれを自動補正することを省略しても本発明は成立する。 In the above-mentioned example, the deviation of the m / z value is automatically corrected, but the present invention is established even if the deviation of the m / z value is omitted.

つまり、本発明によれば、制御条件の変化によるm/z値のずれを検知し、定量分析に必要な感度を維持しつつ、m/z値のずれに対する証憑を与え、m/z値のずれにより減少したイオン量を補うことができる質量分析装置および質量分析方法を実現することができる。 That is, according to the present invention, the deviation of the m / z value due to the change of the control condition is detected, the sensitivity required for the quantitative analysis is maintained, and the voucher for the deviation of the m / z value is given, and the deviation of the m / z value is given. It is possible to realize a mass spectrometer and a mass spectrometry method that can compensate for the amount of ions reduced due to the displacement.

110:試料供給装置、120:質量分析装置、121:イオン化部、122:質量分析部、123:検出部、124:制御部、130:入出力装置、210:校正前のピーク、220:校正後のピーク、230:半値幅、310:校正曲線、410:初期設定ステップ、420:近傍探索ステップ、430:成分量の変化計算ステップ、440:近傍パターン確定判定ステップ、450:近傍パターン計算ステップ、460:イオン強度補正ステップ、470:m/z値補正ステップ、480:終了判定ステップ、510:マスフィルター制御部、520:イオン強度検出部、530:成分量の変化計算部、540:近傍パターン計算部、550:イオン強度補正部、560:m/z値補正部、610:校正時のマススペクトル、620:m/z値M-Δの探索点(探索点1)、630:m/z値Mの探索点(ピークトップ)(探索点2)、640:m/z値M+Δの探索点(探索点3)、650:成分量増加の場合のマススペクトルの例、660:探索点1の測定値、670:探索点2の測定値(ピークトップの測定値)、680:探索点3の測定値、690:1サイクル内の探索点を結んだ線分、710:探索点1の変化、720:探索点2の変化、730:探索点3の変化、810:校正直後のマススペクトル、820:高m/z値側に0.05ずれた状態、830:高m/z値側に0.10ずれた状態、840:探索点1の測定値、850:探索点2の測定値、910:m/z Error(m/z値のずれ)、920:I(イオン強度)、930:R(イオン強度比)、940:R’(イオン強度の測定予想比)、950:Score(スコア)、1010:Scoreとm/z値のずれの関係を示す線分、1110:Target(目的成分の定義領域)、1111:Top(ピートップのm/z値)、1112:Width(ピークの半値幅)、1113:Dwell Time(取り込み時間)、1120:Search(近傍探索条件の定義領域)、1121:Point Number(探索点の点数)、1122:Delta(探索点の間隔)、1123:Threshold(近傍パターン確定判定のしきい値)、1130:Correction(m/z値の補正条件の定義領域)、1131:Step(補正の刻み)、1132:Maximum(補正の最大値)、1210:ログのタイトル、1220:目的成分のm/z値、1230:ログの項目、1240:ログの実体、1310:成分量が一定の場合のマススペクトル、1320:m/z値M-Δ/2の探索点(探索点1)、1330:m/z値M+Δ/2の探索点(探索点max)、1340:成分量が一定の場合の測定されるマススペクトル、1350:成分量増加の場合のマススペクトル、1360:探索点1の測定値、1370:探索点maxの測定値、1380:成分量増加の場合の測定されるマススペクトル、1410:探索点1の変化、1420:探索点maxの変化、1430:ピークトップに対応する探索点の変化、1510:校正直後のマススペクトル、1520:高m/z値側に0.05ずれた状態、1530:高m/z値側に0.10ずれた状態、1540:探索点1の測定値、1550:探索点2の測定値、1560:ピークトップ、1610:接線1、1620:接線2、1630:接線1と接線2の交点 110: Sample supply device, 120: Mass analyzer, 121: Ionization section, 122: Mass analyzer, 123: Detection section, 124: Control section, 130: Input / output device, 210: Peak before calibration, 220: After calibration Peak, 230: Half price width, 310: Calibration curve, 410: Initial setting step, 420: Neighborhood search step, 430: Component amount change calculation step, 440: Neighborhood pattern determination determination step, 450: Neighborhood pattern calculation step, 460 : Ion intensity correction step, 470: m / z value correction step, 480: End determination step, 510: Mass filter control unit, 520: Ion intensity detection unit, 530: Component amount change calculation unit, 540: Neighborhood pattern calculation unit 550: Ion intensity correction unit, 560: m / z value correction unit, 610: Mass spectrum at the time of calibration, 620: M / z value M-Δ search point (search point 1), 630: m / z value M Search point (peak top) (search point 2), 640: search point with m / z value M + Δ (search point 3), 650: example of mass spectrum when the amount of components increases, 660: measured value at search point 1. , 670: measured value of search point 2 (measured value of peak top), 680: measured value of search point 3, line segment connecting search points in 690: 1 cycle, 710: change of search point 1, 720: Change of search point 2, 730: Change of search point 3, 810: Mass spectrum immediately after calibration, 820: 0.05 deviation to the high m / z value side, 830: 0.10 to the high m / z value side Displaced state, 840: measured value at search point 1, 850: measured value at search point 2, 910: m / z Error (deviation of m / z value), 920: I (ion intensity), 930: R (ion). Intensity ratio), 940: R'(expected measurement ratio of ion intensity), 950: Score (score), 1010: Line line showing the relationship between Score and m / z value deviation 1110: Target (definition region of target component) ), 1111: Top (m / z value of pea top), 1112: Wide (half price width of peak), 1113: Dwell Time (uptake time), 1120: Search (definition area of neighborhood search condition), 1121: Point Number. (Score of search points), 1122: Delta (interval of search points), 1123: Threat (threshold value for determination of neighborhood pattern determination), 1130: Direction (definition area of m / z value correction condition), 1131: Step. (Correction step) 1132: Maximum (maximum correction value), 1210: Log Title, 1220: m / z value of target component, 1230: log item, 1240: log substance, 1310: mass spectrum when component amount is constant, 1320: search point of m / z value M-Δ / 2. (Search point 1), 1330: Search point with m / z value M + Δ / 2 (search point max), 1340: Mass spectrum measured when the component amount is constant, 1350: Mass spectrum when the component amount increases, 1360: Measured value of search point 1, 1370: Measured value of search point max, 1380: Mass spectrum measured when the amount of components increases, 1410: Change of search point 1, 1420: Change of search point max, 1430: Change of search point corresponding to peak top, 1510: Mass spectrum immediately after calibration, 1520: 0.05 deviation to high m / z value side, 1530: 0.10 deviation to high m / z value side, 1540: measured value of search point 1, 1550: measured value of search point 2, 1560: peak top, 1610: tangent line 1, 1620: tangent line 2, 1630: intersection of tangent line 1 and tangent line 2.

Claims (15)

試料成分をイオン化するイオン化部と、前記イオン化部により供給されたイオンを選別する質量分析部と、前記質量分析部により選別されたイオンの強度を測定する検出部と、前記イオン化部、前記質量分析部及び前記検出部の動作を制御する制御部とを備える質量分析装置において、
前記制御部は、
目的成分に由来するm/z値の近傍の複数個所のm/z値を選択するため、前記質量分析部の制御条件を切り替える質量分析制御部と、
前記質量分析制御部が選択したm/z値のイオン強度を検出するイオン強度検出部と、
前記イオン強度検出部が検出したイオン強度の変化を計算する成分量の変化計算部と、
成分量一定の状態での前記複数個所のm/z値のイオン強度の比率を求め、近傍パターンを計算する近傍パターン計算部と、
前記近傍パターンに基づき測定したイオン強度の不足を補填するイオン強度補正部と、
を備え、目的成分に由来するm/z値のイオン強度の経時変化を測定する質量分析装置。
An ionization unit that ionizes sample components, a mass spectrometry unit that selects ions supplied by the ionization unit, a detection unit that measures the intensity of ions selected by the mass spectrometry unit, an ionization unit, and the mass spectrometry. In a mass spectrometer including a unit and a control unit that controls the operation of the detection unit.
The control unit
A mass spectrometric control unit that switches the control conditions of the mass spectrometric unit in order to select multiple m / z values in the vicinity of the m / z value derived from the target component.
An ion intensity detecting unit that detects the ion intensity of the m / z value selected by the mass spectrometry control unit, and an ion intensity detecting unit.
A change calculation unit for the amount of components that calculates the change in ionic strength detected by the ion intensity detection unit, and a change calculation unit.
A neighborhood pattern calculation unit that calculates the ratio of the ionic strength of the m / z values at the plurality of locations in a state where the amount of components is constant and calculates the neighborhood pattern.
An ionic strength correction unit that compensates for the lack of ionic strength measured based on the neighborhood pattern, and
A mass spectrometer that measures the change over time in the ionic strength of the m / z value derived from the target component.
請求項1に記載の質量分析装置において、
前記制御部は、前記近傍パターンに基づき前記質量分析制御部の制御条件を修正するm/z値補正部を備えることを特徴とする質量分析装置。
In the mass spectrometer according to claim 1,
The control unit is a mass spectrometer including an m / z value correction unit that corrects the control conditions of the mass spectrometry control unit based on the neighborhood pattern.
請求項1に記載の質量分析装置において、
前記目的成分に由来するm/z値の近傍のm/z値のうち、前記目的成分に由来するm/z値の両端のm/z値は、ピークトップから数%低下した位置であり、最大でもピークの半値以上であることを特徴とする質量分析装置。
In the mass spectrometer according to claim 1,
Of the m / z values in the vicinity of the m / z value derived from the target component, the m / z values at both ends of the m / z value derived from the target component are located at positions lowered by several percent from the peak top. A mass spectrometer characterized by having a maximum of more than half of the peak.
請求項1に記載の質量分析装置において、
前記質量分析制御部は、前記目的成分に由来するm/z値の近傍の複数個所のm/z値を選択するため、前記目的成分に由来するm/z値-0.5から+0.5の範囲にてスキャン測定可能であることを特徴とする質量分析装置。
In the mass spectrometer according to claim 1,
Since the mass spectrometric control unit selects m / z values at a plurality of locations near the m / z value derived from the target component, the m / z value derived from the target component is −0.5 to +0.5. A mass spectrometer characterized by being capable of scan measurement within the range of.
請求項1に記載の質量分析装置において、
前記成分量の変化計算部は、イオン強度の変化を複数のサイクルにより計算し、前後のサイクルから補間法で成分量の変化を計算することを特徴とする質量分析装置。
In the mass spectrometer according to claim 1,
The component amount change calculation unit is a mass spectrometer characterized in that the change in ionic strength is calculated by a plurality of cycles and the change in the component amount is calculated from the previous and next cycles by an interpolation method.
請求項5に記載の質量分析装置において、
前記近傍パターン計算部は、前記補間法で計算された成分量の変化から、単位成分量あたりのイオン強度を近傍パターンとすることを特徴とする質量分析装置。
In the mass spectrometer according to claim 5,
The neighborhood pattern calculation unit is a mass spectrometer characterized in that the ionic strength per unit component amount is set as a neighborhood pattern from the change in the component amount calculated by the interpolation method.
請求項2に記載の質量分析装置において、
前記近傍パターン計算部は、登録された近傍パターンとスコア、m/z値のずれの関係を参照しスコアを決定し、前記イオン強度補正部は、前記スコアに対応するイオン強度比からイオン強度を補正し、前記m/z値補正部は、前記スコアに対応するm/z値のずれを補正することを特徴とする質量分析装置。
In the mass spectrometer according to claim 2,
The neighborhood pattern calculation unit determines the score by referring to the relationship between the registered neighborhood pattern, the score, and the deviation of the m / z value, and the ion intensity correction unit determines the ion intensity from the ion intensity ratio corresponding to the score. The m / z value correction unit corrects the mass spectrometer, and corrects the deviation of the m / z value corresponding to the score.
請求項4に記載の質量分析装置において、
成分量の変化計算部は、前記スキャン測定によるマススペクトルに外接する接線により成分量の変化を求めることを特徴とする質量分析装置。
In the mass spectrometer according to claim 4,
The component amount change calculation unit is a mass spectrometer characterized in that the change in the component amount is obtained by a tangent line tangent to the mass spectrum obtained by the scan measurement.
請求項4に記載の質量分析装置において、
成分量の変化計算部は、前記スキャン測定によるマススペクトルに外接する2次関数により成分量の変化を求めることを特徴とする質量分析装置。
In the mass spectrometer according to claim 4,
The component amount change calculation unit is a mass spectrometer characterized in that the change in the component amount is obtained by a quadratic function circumscribing the mass spectrum obtained by the scan measurement.
請求項2に記載の質量分析装置において、
前記イオン強度検出部は、前記目的成分に由来するm/z値と同じずれを示すイオンの近傍を探索し、前記近傍パターン計算部により、前記探索したイオンのイオン強度の近傍パターンを計算し、目的成分のイオン強度を前記イオン強度補正部により補正し、m/z値補正部によりm/z値を補正することを特徴とする質量分析装置。
In the mass spectrometer according to claim 2,
The ion intensity detection unit searches for the vicinity of ions showing the same deviation as the m / z value derived from the target component, and the neighborhood pattern calculation unit calculates the vicinity pattern of the ion intensity of the searched ions. A mass spectrometer characterized in that the ion intensity of a target component is corrected by the ion intensity correction unit and the m / z value is corrected by the m / z value correction unit.
試料成分をイオン化し、イオン化されたイオンを選別し、選別したイオンの強度を測定する質量分析方法において、
目的成分に由来するm/z値の近傍の複数個所のm/z値のイオン強度を連続して検出する近傍探索ステップと、
検出したイオン強度の変化を計算する成分量の変化計算ステップと、
成分量一定の状態での前記目的成分に由来するm/z値の近傍の複数個所のm/z値のイオン強度の比率を求める近傍パターン計算ステップと、
前記近傍パターンに基づき測定したイオン強度の不足を補填するイオン強度補正ステップと、
を備え、目的成分に由来するm/z値のイオン強度の経時変化を測定することを特徴とする質量分析方法。
In a mass spectrometry method in which a sample component is ionized, ionized ions are selected, and the intensity of the selected ions is measured.
A neighborhood search step that continuously detects the ionic strength of m / z values at multiple locations near the m / z values derived from the target component, and
The change calculation step of the component amount to calculate the change of the detected ionic strength, and
A neighborhood pattern calculation step for obtaining the ratio of the ionic strength of the m / z value at a plurality of locations near the m / z value derived from the target component in a state where the component amount is constant, and
An ionic strength correction step that compensates for the lack of ionic strength measured based on the neighborhood pattern, and
A mass spectrometric method comprising the above, and measuring the change with time of the ionic strength of the m / z value derived from the target component.
請求項11に記載の質量分析方法において、
近傍の複数個所のm/z値のずれに基づき、検出したm/z値を補正するm/z値補正ステップをさらに備えることを特徴とする質量分析方法。
In the mass spectrometry method according to claim 11,
A mass spectrometric method further comprising an m / z value correction step for correcting a detected m / z value based on a deviation of m / z values at a plurality of nearby locations.
請求項11に記載の質量分析方法において、
前記目的成分に由来するm/z値の近傍のm/z値のうち、前記目的成分に由来するm/z値の両端のm/z値は、ピークトップから数%低下した位置であり、最大でもピークの半値以上であることを特徴とする質量分析方法。
In the mass spectrometry method according to claim 11,
Of the m / z values in the vicinity of the m / z value derived from the target component, the m / z values at both ends of the m / z value derived from the target component are located at positions lowered by several percent from the peak top. A mass spectrometric method characterized by having a maximum of more than half of the peak.
請求項11に記載の質量分析方法において、
前記目的成分に由来するm/z値の近傍の複数個所のm/z値を選択するため、前記目的成分に由来するm/z値-0.5から+0.5の範囲にてスキャン測定可能であることを特徴とする質量分析方法。
In the mass spectrometry method according to claim 11,
Since the m / z values in the vicinity of the m / z value derived from the target component are selected, the m / z value derived from the target component can be scanned and measured in the range of -0.5 to +0.5. A mass spectrometric method characterized by being.
請求項11に記載の質量分析方法において、
イオン強度の変化を複数のサイクルにより計算し、前後のサイクルから2次式で成分量の変化を計算することを特徴とする質量分析方法。
In the mass spectrometry method according to claim 11,
A mass spectrometric method characterized in that changes in ionic strength are calculated by a plurality of cycles, and changes in the amount of components are calculated by a quadratic equation from the cycles before and after.
JP2021514815A 2019-04-16 2020-03-05 Mass spectrometer and mass spectrometry method Active JP7077481B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2019077796 2019-04-16
JP2019077796 2019-04-16
PCT/JP2020/009456 WO2020213283A1 (en) 2019-04-16 2020-03-05 Mass spectrometer and mass spectrometry method

Publications (2)

Publication Number Publication Date
JPWO2020213283A1 JPWO2020213283A1 (en) 2021-11-25
JP7077481B2 true JP7077481B2 (en) 2022-05-30

Family

ID=72838279

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021514815A Active JP7077481B2 (en) 2019-04-16 2020-03-05 Mass spectrometer and mass spectrometry method

Country Status (2)

Country Link
JP (1) JP7077481B2 (en)
WO (1) WO2020213283A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230402273A1 (en) 2021-01-22 2023-12-14 Hitachi High-Tech Corporation Mass Spectrometer

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070200060A1 (en) 2006-02-28 2007-08-30 Russ Charles W Iv Pulsed internal lock mass for axis calibration
JP2014508937A (en) 2011-03-07 2014-04-10 マイクロマス・ユーケイ・リミテッド Dynamic resolution correction of a quadrupole mass spectrometer.
US20160211125A1 (en) 2012-12-05 2016-07-21 Tofwerk Ag Method of calibrating mass-to-charge ratio measurements obtained from a mass spectrometer connected in fluid communication with an analysis system delivering a temporally changing sample

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070200060A1 (en) 2006-02-28 2007-08-30 Russ Charles W Iv Pulsed internal lock mass for axis calibration
JP2014508937A (en) 2011-03-07 2014-04-10 マイクロマス・ユーケイ・リミテッド Dynamic resolution correction of a quadrupole mass spectrometer.
US20160211125A1 (en) 2012-12-05 2016-07-21 Tofwerk Ag Method of calibrating mass-to-charge ratio measurements obtained from a mass spectrometer connected in fluid communication with an analysis system delivering a temporally changing sample

Also Published As

Publication number Publication date
JPWO2020213283A1 (en) 2021-11-25
WO2020213283A1 (en) 2020-10-22

Similar Documents

Publication Publication Date Title
CN108603867B (en) Peak detection method and data processing apparatus
Liland et al. Optimal choice of baseline correction for multivariate calibration of spectra
US8368010B2 (en) Quadrupole mass spectrometer
CN102262076B (en) Method for laser-induced breakdown spectroscopy element concentration determination based on spectral line combination
JP7077481B2 (en) Mass spectrometer and mass spectrometry method
US10600625B2 (en) Method of calibrating a mass spectrometer
CN105021546A (en) Method for measuring chemical elements by whole-spectrum direct-reading spectrometer
JP4730482B2 (en) Quadrupole mass spectrometer
JP2004093241A (en) Gas sensor characteristic compensation device and gas concentration measurement device
KR102007241B1 (en) System and method for gas analysis value correction of the standard insulating oil
CN107314999B (en) Liquid cathode glow discharge spectral analysis method based on multiple linear regression method
CN110849827B (en) Spectrum calibration method, device and equipment for echelle spectrometer and storage medium
JP2001324476A (en) Inductively-coupled plasma mass spectrometeric analysis method
TWI589871B (en) Method for determining the maximum of the mass peak in mass spectrometry
US20220373522A1 (en) Waveform Analytical Method and Waveform Analytical Device
CN113567602A (en) Detection method and detection device
JPH11162400A (en) Process gas monitor device
JP2019090757A (en) Processing device, inspection device and processing method
US20220163391A1 (en) Diagnostic testing method for a spectrometer
EP4062141B1 (en) A method of detecting a deviation in operation of a pressure sensor for internal combustion piston engine and computer readable memory device for detecting a deviation in operation of a pressure sensor
CN116613613B (en) Laser calibration method, device, equipment and storage medium of laser
US20230154735A1 (en) Mass spectrometer isolation profile analyser
US8754361B1 (en) Systems and methods for adjusting a mass spectrometer output
Abbott et al. Tune Computation via Model Fitting to Swept Machine Response Measurement
JP6505268B1 (en) Mass spectrometer and mass spectrometry method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210720

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220419

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220518

R150 Certificate of patent or registration of utility model

Ref document number: 7077481

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150