(57)【特許請求の範囲】
【請求項１】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)周回軌道に沿ってイオンを多重周回させることなく又は多重周回させる場合でも異種のイオンの追いつき・追い越しが起こらないことが保証される周回数で以てイオンを飛行させる第１測定モードで被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて非追い越し飛行時間スペクトルを取得する第１測定モード実行ステップと、
b)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する第２測定モードで前記被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて追い越し飛行時間スペクトルを取得する第２測定モード実行ステップと、
c)第２測定モードの実行により取得された追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と、目的のマススペクトル上の各ピークの強度を要素とする行列に乗じることで、該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と等しくなる係数行列を求め、前記目的のマススペクトルを求めるために前記係数行列の逆行列を演算する演算処理ステップと、
を有し、前記演算処理ステップは、前記逆行列の正則化一般逆行列を用い、推定される目的のマススペクトルと、前記非追い越し飛行時間スペクトルで近似した真のマススペクトルとの自乗誤差を最小にするように一般逆行列を求めることを特徴とする質量分析方法。
【請求項２】
請求項１に記載の質量分析方法であって、
ｍ行ｎ列の係数行列をＡとしたときに次式で示す正則化一般逆行列を用い、該式中の正則化パラメータλをλ＝σ（ノイズの標準偏差）、ペナルティ行列ＰをＰ＝（Ｘ0Ｘ0T）-1（ただし、Ｘ0は非追い越し飛行時間スペクトルから求まるマススペクトル）とすることにより、推定される目的のマススペクトルと非追い越し飛行時間スペクトルから求まるマススペクトルとの平均自乗誤差を最小にするようにしたことを特徴とする質量分析方法。
ｍ≦ｎである場合
Ａ+＝ＡT（ＡＡT＋λ2Ｐ）-1
ここでＰはｍ行ｍ列
ｍ＞ｎである場合
Ａ+＝（ＡTＡ＋λ2Ｐ）-1ＡT
ここでＰはｎ行ｎ列
【請求項３】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行ステップと、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、前記目的のマススペクトルを求めるために該係数行列の逆行列を正則化一般逆行列を用いて演算する演算処理ステップと、
を有することを特徴とする質量分析方法。
【請求項４】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行ステップと、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、該係数行列の各要素と追い越し飛行時間スペクトル上の各ピークの強度とに基づき、点拡がり関数に対するベイズ的逐次近似解法を用いて目的のマススペクトル上の各質量電荷比におけるピーク強度を算出する演算処理ステップと、
を有することを特徴とする質量分析方法。
【請求項５】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)周回軌道に沿ってイオンを多重周回させることなく又は多重周回させる場合でも異種のイオンの追いつき・追い越しが起こらないことが保証される周回数で以てイオンを飛行させる第１測定モードで被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて非追い越し飛行時間スペクトルを取得する第１測定モード実行手段と、
b)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する第２測定モードで前記被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて追い越し飛行時間スペクトルを取得する第２測定モード実行手段と、
c)第２測定モードの実行により取得された追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と、目的のマススペクトル上の各ピークの強度を要素とする行列に乗じることで、該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と等しくなる係数行列を求め、前記目的のマススペクトルを求めるために前記係数行列の逆行列を演算する手段であって、前記逆行列の正則化一般逆行列を用い、推定される目的のマススペクトルと前記非追い越し飛行時間スペクトルで近似した真のマススペクトルとの自乗誤差を最小にするように一般逆行列を求める演算処理手段と、
を備えることを特徴とする質量分析装置。
【請求項６】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行手段と、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、前記目的のマススペクトルを求めるために該係数行列の逆行列を正則化一般逆行列を用いて演算する演算処理手段と、
を備えることを特徴とする質量分析装置。
【請求項７】
イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行手段と、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、該係数行列の各要素と追い越し飛行時間スペクトル上の各ピークの強度とに基づき、点拡がり関数に対するベイズ的逐次近似解法を用いて目的のマススペクトル上の各質量電荷比におけるピーク強度を算出する演算処理手段と、
を備えることを特徴とする質量分析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、試料由来のイオンを閉じた周回軌道に沿って繰り返し飛行させることでイオンを質量電荷比（m/z値）に応じて分離して検出する多重周回飛行時間型の質量分析装置、及び多重周回飛行時間型質量分析装置を用いた質量分析方法に関する。
【背景技術】
【０００２】
一般に、飛行時間型質量分析装置（Time of Flight Mass Spectrometer、以下「ＴＯＦＭＳ」と略す）では、一定のエネルギーを与えることで加速したイオンはそれぞれ質量電荷比に応じた飛行速度を持つ、という原理に基づき、そうしたイオンが一定距離を飛行するのに要する飛行時間を計測し、その飛行時間を質量電荷比に換算することにより質量分析を行う。したがって、質量分解能を向上させるためには飛行距離を長くすればよいが、直線的な飛行距離を延ばすには装置を大型化する必要がある。そこで、飛行距離の延伸と装置の小形化とを両立させるため、略円形状、略楕円形状、略８の字形状など様々な態様の閉軌道に沿ってイオンを繰り返し飛行させるようにした多重周回飛行時間型質量分析装置（Multi Turn ＴＯＦＭＳ、以下「ＭＴ−ＴＯＦＭＳ」と略す）が開発されている（例えば特許文献１、２など参照）。
【０００３】
また、同様の目的で、上記のような周回軌道ではなく反射電場によりイオンを複数回反射させる往復軌道とすることで飛行距離を延ばすようにした、多重反射飛行時間型質量分析装置も開発されている。多重周回飛行時間型と多重反射飛行時間型とではイオン光学系が相違するものの、質量分解能を向上させるための基本的な原理は同じである。そこで、本明細書では、「多重周回飛行時間型」は「多重反射飛行時間型」を包含するものとする。
【０００４】
上述のようにＭＴ−ＴＯＦＭＳは高い質量分解能を達成することができるものの、イオンの飛行経路が閉軌道であることを原因とする欠点が存在する。即ち、周回数が増加するに伴い、低質量電荷比で速度の大きなイオンが高質量電荷比で速度の小さなイオンに追いつき、追い越してしまうという問題である。このように異なる質量電荷比のイオンの追い越しが生じると、取得された飛行時間スペクトル上では、観測されるピーク毎にそのピークに対応するイオンの周回数が異なる、即ち、飛行距離が異なる、という状態が起こる。こうなると、イオンの質量電荷比と飛行距離とを一意的に決定することができないため、飛行時間スペクトルを直接的にマススペクトルに換算することができなくなる。
【０００５】
上記欠点のため、従来のＭＴ−ＴＯＦＭＳでは、イオン源で生成された試料由来のイオンの中で上記のような追い越しの起こらない質量範囲に限定したイオンを選別し、選別したイオンを周回軌道に乗せて所定周回数だけ飛行させた後に検出する、という制御を行うのが一般的である。このような手法では、高質量分解能のマススペクトルを得ることはできるものの、そのマススペクトルの質量範囲はかなり限られたものとなる。これは、一度の分析で比較的広い質量範囲のマススペクトルを得られる、というＴＯＦＭＳの利点に反する。
【０００６】
そこで、飛行時間を質量電荷比に換算するために、飛行時間スペクトル上に現れるピークの周回数を推定する様々な手法が従来提案されている。例えば特許文献３には、同一試料に対し条件の相違する複数回の質量分析を実行して得られた結果を比較することで、飛行時間スペクトル上に現れるピークの周回数を推定する方法が提案されている。こうした手法は有効ではあるものの、データ処理が複雑になることは避けられない。また、特に試料に含まれる成分の数が多い場合に周回数の推定が困難である。
【０００７】
特許文献４には、周回軌道を周回しているイオンを該軌道から出射させて検出器に向かわせるタイミングを変えた（つまり同一種のイオンの周回数が相違する）複数の飛行時間スペクトルを取得し、それら飛行時間スペクトルの多重相関関数を計算することによって単一周回数の飛行時間スペクトルを再構成する方法が提案されている。この方法では、複数個の飛行時間スペクトルデータに対し、周回軌道上の飛行時間がＴであるイオンの信号強度Ｇ（Ｔ）を求めるために、次の(1)式が用いられる。
【数１】
0005533255.tif 000002
ここで、Ｆj（j＝１、２、…、ｒ）は測定データから探した周回数Ｎjのイオンの信号強度、ｙは飛行時間のずれ、ｙlはこのずれ時間の下限値、ｙuはずれ時間の上限値、Ｈは変数Ｆjの値により決まる関数である。関数Ｈの具体例として算術平均、最小値、相乗平均、調和平均などが挙げられているが、Ｆjが偶然大きな値を持った偽ピークを排除するためには、様々な大きさのＦjの中で大きな値ではなく小さな値のほうがＨに大きく反映されるような、関数Ｈの定義が望ましいとされている。
【０００８】
上記特許文献４によれば、周回軌道からのイオンの出射タイミングによっては、周回軌道上を飛行しているイオンの一部が検出されない無効期間が生じることが指摘されている。しかしながら、特にそのための対策については提案されていない。この無効期間は、周回軌道からイオンを出射させるためのゲート電極が有限の長さを有し、周回しているイオンを出射させようとする時点でゲート電極を通過中であるイオンは適切な方向（つまり検出器により検出可能な方向）に出射されないために生じるものである。上述したように、大きな値ではなく小さな値のＦjを重視して定義されたＨに従った処理を行った場合、具体的には相乗平均や調和平均が採用された場合には、ピーク強度Ｆjの中に１つでも強度０のものが存在すれば再構成される飛行時間スペクトル中にそのピークは含まれなくなる。即ち、上記のような無効期間のために周回数の異なる複数回の質量分析の中で観測されなかったイオンについては、そのピークは再構成スペクトルには現れなくなってしまう。そのため、イオンの検出見逃しが発生し、試料中に含まれる筈の成分が含まれないとの誤判断が下される原因となり得る。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２００６−２２８４３５号公報
【特許文献２】特開２００８−２７６８３号公報
【特許文献３】特開２００５−１１６３４３号公報
【特許文献４】特開２００５−７９０４９号公報
【特許文献５】国際公開第２０１０／０５２７５６号パンフレット
【発明の概要】
【発明が解決しようとする課題】
【００１０】
上述した従来技術の問題に対し本願発明者は、特許文献５にて新規の質量分析方法を提案している。この方法は、多重周回飛行時間スペクトル上で検出されたイオン強度値を、イオン追い越しのない非周回モードで測定された飛行時間スペクトルを利用して、周回数が補正された各質量電荷比に分配することにより、マススペクトルとして再構成するものである。この方法では、周回軌道からのイオン出射タイミングを変えた多重周回飛行時間スペクトルからそれぞれ再構成した複数の再構成マススペクトルが得られるから、これらを同時に表示することができる。したがって、分析者はその複数の再構成マススペクトルを見比べることで、周回軌道からの出射タイミングが原因で正常に検出されなかったイオンについてもその存在を確認することができるという大きな利点がある。
【００１１】
飛行するイオンの種類が比較的少ない場合には、上記の新規方法は優れているものの、イオンの種類が多くなると問題が生じる場合がある。詳しくは後述するが、本願発明者の検討によれば、イオンの種類が多く多重周回飛行時間スペクトル上で複数のイオン由来のピークが重なったときに、その重なったピークのイオン強度値の分配が適切に行われず、再構成マススペクトルの精度が落ちる場合があることが判明した。
【００１２】
本発明は上記課題を解決するためになされたものであり、その目的とするところは、多重周回により得られる複数の飛行時間スペクトルを再構成する際に、周回軌道からの出射タイミングが原因で正常に検出されなかったイオンについてもその存在の確認を可能とすることで、成分の見逃し・見落としを防止することができ、特に多種のイオンが存在する複雑な飛行時間スペクトルに対しても、ピーク強度を正しく求め高精度のマススペクトルを得ることができる多重周回飛行時間型の質量分析装置及び質量分析方法を提供することである。
【課題を解決するための手段】
【００１３】
上記課題を解決するために成された第１発明は、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)周回軌道に沿ってイオンを多重周回させることなく又は多重周回させる場合でも異種のイオンの追いつき・追い越しが起こらないことが保証される周回数で以てイオンを飛行させる第１測定モードで被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて非追い越し飛行時間スペクトルを取得する第１測定モード実行ステップと、
b)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する第２測定モードで前記被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて追い越し飛行時間スペクトルを取得する第２測定モード実行ステップと、
c)第２測定モードの実行により取得された追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と、目的のマススペクトル上の各ピークの強度を要素とする行列に乗じることで、該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と等しくなる係数行列を求め、前記目的のマススペクトルを求めるために前記係数行列の逆行列を演算する演算処理ステップと、
を有し、前記演算処理ステップは、前記逆行列の正則化一般逆行列を用い、推定される目的のマススペクトルと、前記非追い越し飛行時間スペクトルで近似した真のマススペクトルとの自乗誤差を最小にするように一般逆行列を求めることを特徴としている。
【００１４】
第１発明に係る質量分析方法では、具体的には、ｍ行ｎ列の係数行列をＡとしたときに次式で示す正則化一般逆行列を用い、該式中の正則化パラメータλをλ＝σ（ノイズの標準偏差）、ペナルティ行列ＰをＰ＝（Ｘ0Ｘ0T）-1（ただし、Ｘ0は非追い越し飛行時間スペクトルから求まるマススペクトル）とすることにより、推定される目的のマススペクトルと非追い越し飛行時間スペクトルから求まるマススペクトルとの平均自乗誤差を最小にすることができる。
ｍ≦ｎである場合
Ａ+＝ＡT（ＡＡT＋λ2Ｐ）-1
ここでＰはｍ行ｍ列
ｍ＞ｎである場合
Ａ+＝（ＡTＡ＋λ2Ｐ）-1ＡT
ここでＰはｎ行ｎ列
【００１５】
また上記課題を解決するために成された第２発明は、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行ステップと、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、前記目的のマススペクトルを求めるために該係数行列の逆行列を正則化一般逆行列を用いて演算する演算処理ステップと、
を有することを特徴としている。
【００１６】
多重周回飛行時間型質量分析装置では、周回軌道からのイオンの出射時間が決まると、飛行時間と質量電荷比との一対多の関係が理論的に求まる。即ち、或る飛行時間に対応してどのような質量電荷比のイオンが現れる（ただし、周回数は質量電荷比によって異なる）可能性があるのかが分かるから、飛行時間スペクトル上の或る飛行時間におけるピーク強度は、マススペクトル上の複数の質量電荷比におけるピーク強度の線形結合で表すことができる。このため、追い越し飛行時間スペクトル上の各飛行時間におけるピーク強度を要素とする行列Ｙと目的のマススペクトル上の各質量電荷比におけるピーク強度を要素とする行列Ｘとは、係数行列をＡとしてＹ＝ＡＸで表すことができ、Ａは値が１である幾つかの要素を含むスパースな行列である。マススペクトルを求めることは行列Ｘの要素を求めることであるから、係数行列Ａの逆行列Ａ-1を求めれば観測した飛行時間スペクトルからマススペクトルが得られる。
【００１７】
上記のような逆行列を求めるために正則化一般逆行列を用いるのが有効であるが、第１発明においては、特に正則化一般逆行列のペナルティ行列に非追い越し飛行時間スペクトルに基づく制約を付ける。具体的には、逆行列の正則化一般逆行列を用い、推定される目的のマススペクトルと非追い越し飛行時間マススペクトルで近似した真のマススペクトルとの自乗誤差を最小にするように一般逆行列を求める。これにより、目的のマススペクトルを正確に推定することができる。
【００１８】
一方、或る１つの質量電荷比のイオンは飛行に伴い時間的に拡がって複数の飛行時間で観測されることから、第２発明では、上記拡がりを与える点拡がり関数を係数行列に導入する。そして、周回数の異なる２以上の多重周回飛行時間スペクトルを連立行列式とし、これを解くことで一般化逆行列を求めてマススペクトルを算出する。これにより、同種の（同一質量電荷の）イオンを異なる条件の下で測定した複数の飛行時間スペクトル上のピークに対し、矛盾のないマススペクトル上のピークを得ることができ、正確なマススペクトルの算出が可能である。ただし、通常、連立行列式の係数行列は巨大な行列になるため、ハードウエアの制約などから一般化逆行列を計算するのが困難である場合が多い。そこで、これに代わる解法としてベイズ的逐次近似解法を利用することができる。
【００１９】
即ち、上記課題を解決するために成された第３発明は、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置を用いた質量分析方法であって、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行ステップと、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、該係数行列の各要素と追い越し飛行時間スペクトル上の各ピークの強度とに基づき、点拡がり関数に対するベイズ的逐次近似解法を用いて目的のマススペクトル上の各質量電荷比におけるピーク強度を算出する演算処理ステップと、
を有することを特徴としている。
【００２０】
点拡がり関数に対するベイズ的逐次近似解法は、一般にイメージプロセッシングなどに用いられている手法であるが、これを利用して目的のマススペクトル上の各質量電荷比におけるピーク強度を算出することができる。具体的には、マススペクトル上でm/z＝ｉに真のピークが生じる確率をｆ(ｉ)、m/z＝ｉのピークが飛行時間スペクトル上の時間tof＝ｋで観測される点拡がり関数をｈ(ｋ｜ｉ)としたとき、次の漸化式が得られる。
ｆr+1(ｉ）＝ｆr(ｉ)Σk｛ｈ(ｋ｜ｉ)ｆ(ｉ)・ｈ(ｋ)／｛Σjｆr(ｊ)・ｈ(ｋ｜ｊ)｝
ｒ→∞のときにｆr(ｉ)→ｆ(ｉ)に収束するとすれば、上記式でｒ→∞として求めたｆr(ｉ)が求めるｆ(ｉ)、つまり目的の質量電荷比におけるピーク強度となる。もちろん、実用的にはｒ→∞にする必要はなく、ｒを適宜な値とした段階で近似を打ち切ることができる。
【００２１】
また上記課題を解決するために成された第４発明は第１発明に係る質量分析方法を実施するための装置であって、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)周回軌道に沿ってイオンを多重周回させることなく又は多重周回させる場合でも異種のイオンの追いつき・追い越しが起こらないことが保証される周回数で以てイオンを飛行させる第１測定モードで被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて非追い越し飛行時間スペクトルを取得する第１測定モード実行手段と、
b)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する第２測定モードで前記被測定試料に対する質量分析を実行し、検出器により得られる検出信号に基づいて追い越し飛行時間スペクトルを取得する第２測定モード実行手段と、
c)第２測定モードの実行により取得された追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と、目的のマススペクトル上の各ピークの強度を要素とする行列に乗じることで、該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列と等しくなる係数行列を求め、前記目的のマススペクトルを求めるために前記係数行列の逆行列を演算する手段であって、前記逆行列の正則化一般逆行列を用い、推定される目的のマススペクトルと前記非追い越し飛行時間スペクトルで近似した真のマススペクトルとの自乗誤差を最小にするように一般逆行列を求める演算処理手段と、
を備えることを特徴としている。
【００２２】
また上記課題を解決するために成された第５発明は第２発明に係る質量分析方法を実施するための装置であって、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行手段と、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、前記目的のマススペクトルを求めるために該係数行列の逆行列を正則化一般逆行列を用いて演算する演算処理手段と、
を備えることを特徴としている。
【００２３】
また上記課題を解決するために成された第６発明は第３発明に係る質量分析方法を実施するための装置であって、イオン源から出発した各種イオンを周回軌道に沿って複数回繰り返し飛行させた後に検出器に導入し、その検出信号に基づいてマススペクトルを取得する多重周回飛行時間型の質量分析装置において、
a)前記周回軌道上でイオンの追い越しが生じるように多重周回させた後の所定の時点以降にイオンを周回軌道から離脱させて検出器に導入する測定モードによる質量分析を、周回軌道からのイオンの離脱のタイミングを変えて被測定試料に対して２回以上実行し、検出器により得られる検出信号に基づいて２以上の追い越し飛行時間スペクトルを取得する測定実行手段と、
b)取得された２以上の追い越し飛行時間スペクトルの中の１つと目的のマススペクトルとの関係を、目的のマススペクトル上の各ピークの強度を要素とする行列に、飛行時間における拡がりを与える点拡がり関数を導入した係数行列を乗じて該追い越し飛行時間スペクトル上の各ピークの強度を要素とする行列を求める関係に整理し、前記取得された２以上の追い越し飛行時間スペクトルに対する前記関係を連立させた係数行列を求め、該係数行列の各要素と追い越し飛行時間スペクトル上の各ピークの強度とに基づき、点拡がり関数に対するベイズ的逐次近似解法を用いて目的のマススペクトル上の各質量電荷比におけるピーク強度を算出する演算処理手段と、
を備えることを特徴としている。
【発明の効果】
【００２４】
第１乃至第３発明に係る質量分析方法及び第４乃至第６発明に係る質量分析装置によれば、周回軌道に沿ってイオンを多数回繰り返し飛行させることで長い飛行距離を確保して得られる飛行時間スペクトルに基づいて、高い質量分解能及び質量精度であって、しかもピーク強度も正確なマススペクトルを求めることができる。また、被測定試料に含まれる成分の数が多い場合や取得される多重周回飛行時間スペクトルの数が少ないような場合であっても、正確で高い質量分解能のマススペクトルを得ることができる。
【図面の簡単な説明】
【００２５】
【図１】本発明の一実施例によるＭＴ−ＴＯＦＭＳの概略構成図。
【図２】本実施例のＭＴ−ＴＯＦＭＳにおける分析手順を示すフローチャート。
【図３】係数行列の構造の概念図。
【図４】自乗誤差最小法により推定されたマススペクトルを示す図。
【図５】点拡がり関数を導入した場合の係数行列の構造の概念図。
【図６】ベイズ的逐次近似法により推定されたマススペクトルを示す図。
【図７】周回数１の飛行時間スペクトルをシミュレートした結果を示す図。
【図８】図７の飛行時間スペクトルの横軸をm/zに変換して作成したマススペクトルを示す図。
【図９】図７に示した飛行時間スペクトルにおける最大ピーク近辺の拡大図（ａ）及び図８に示したマススペクトルにおける最大ピーク近辺の拡大図（ｂ）。
【図１０】m/z150であるイオンが100周回した後のタイミングで周回軌道からイオンが出射される場合をシミュレートした飛行時間スペクトルを示す図。
【図１１】図１０に示した条件におけるm/zと周回数との関係を示す図。
【図１２】図１０に示した条件におけるm/zとtof値との関係を示す図。
【図１３】適切に変換が行われた場合のマススペクトルを示す図。
【図１４】図１３に示したマススペクトルにおける最大ピーク付近の拡大図。
【図１５】従来方法の一つを用いて図１０に示した飛行時間スペクトルを再構成して得られたマススペクトルを示す図。
【図１６】m/z179とm/z138の２種のイオンのピークが重なったtof値のピークの強度分配を示す図。
【発明を実施するための形態】
【００２６】
本発明に係る質量分析方法を適用する質量分析装置の一実施例を、添付図面を参照して説明する。図１は本実施例のＭＴ−ＴＯＦＭＳの概略構成図である。
【００２７】
イオン源１では試料分子がイオン化され、生成された各種イオンは所定のエネルギーを付与されて飛行を開始する。なお、イオン源１は例えば３次元四重極型イオントラップなどのように、外部で生成された各種イオンを一時的に保持し、それらイオンに所定のタイミングで一斉にエネルギーを付与して飛行を開始させるものでもよい。
【００２８】
イオン源１を出発点として飛行を開始したイオンは、ゲート電極２により形成される偏向電場を介し周回軌道５に導入される。この周回軌道５は、複数のトロイダル扇形電極４によりそれぞれ生成される電場の作用で形成されるものである。なお、周回軌道５の形状はこれに限られるものではなく、８の字形状など、様々な形状が実現可能である。また、完全に閉じた周回軌道でなくても、直線状や曲線状の往復軌道や、徐々に位置がずれてゆく螺旋軌道などでもよい。
【００２９】
イオンは周回軌道５を１周又は複数周周回した後にゲート電極２により形成される偏向電場を介して周回軌道５から離脱（出射）され、外側に設けられた検出器８に到達して検出される。この例では、周回軌道５への入射用と周回軌道５から出射用のゲート電極２とは共通になっているが、これは別であってもよいし、扇形電極４の一部がゲート電極２の代わりに使用される構成とすることもできる。
【００３０】
検出器８による検出信号はデータ処理部９に入力され、ここで、イオンがイオン源１を出射した時点から最終的に検出器８に到達するまでの飛行時間が計測され、飛行時間スペクトルが作成される。さらに後述するようなデータ解析処理が実行されることで、飛行時間スペクトルからマススペクトルが作成される。入出射電圧発生部１１はゲート電極２に所定のタイミングで、周回軌道５へイオンを入射するための偏向電圧、及び周回軌道５からイオンを出射するための偏向電圧を印加するものである。周回飛行電圧発生部１２は、複数の扇形電極４にそれぞれ所定の電圧を印加することで扇形電場を形成させるものである。これら電圧発生部１１、１２やイオン源１はいずれも制御部１０の制御の下に動作する。また制御部１０には、ユーザが分析のための各種パラメータを入力する入力部１３、及び、マススペクトル等の分析結果を表示する表示部１４が接続されている。
【００３１】
図１において、入射軌道６はイオン源１から周回軌道５への入射点までのイオンの飛行経路であり、その距離がＬinである。出射軌道７は周回軌道５からのイオン出射点から検出器８までのイオンの飛行経路であり、その距離はＬoutである。また、周回軌道５の１周の飛行距離、つまり周回長はＬturnである。
【００３２】
本実施例のＭＴ−ＴＯＦＭＳにおいて、イオン源１から出射された各種イオンは、入射軌道６を通り、入出射電圧発生部１１により制御されるゲート電極２を介して周回軌道５に導入され、周回飛行電圧発生部１２により制御される扇形電極４により形成される扇形電場に従って周回軌道５を１乃至複数回周回飛行する。イオン源１からの各種イオンの出射はほぼ同時であるが、質量電荷比が小さなイオンほど飛行速度が速いため、周回軌道５に沿ってイオンが周回を重ねるに従い、異なる質量電荷比のイオンの間隔は開いてゆく。周回軌道５上をイオンが周回している状態で、入出射電圧発生部１１からゲート電極２への印加電圧がイオン出射用電圧に切り替えられると、ゲート電極２へ到達する順にイオンは周回軌道５から離れ、出射軌道７を通って検出器８に到達する。
【００３３】
周回軌道５は閉じているから、周回軌道５に沿ってイオンが周回を重ねてゆくと、或る時点で、導入されたイオンの中で飛行速度の最も速い低質量電荷比のイオンが飛行速度の最も遅い高質量電荷比のイオンに追いつき、追い越す。その時点までは、速度最小のイオンがゲート電極２を通過してから速度最大のイオンがゲート電極２に達するまでの間にゲート電極２への印加電圧を出射用電圧に切り替えることで、質量電荷比の小さなイオンから順に周回軌道５を離脱させて検出器８に到達させることができる。このようなことが保証される測定モードが第１測定モードであり、第１測定モードの下で得られた飛行時間スペクトルを非追い越し飛行時間スペクトルという。この場合、飛行時間スペクトル上に現れる全てのピークに対応したイオンの周回数は同一である。
【００３４】
これに対し、周回軌道５上で１種類のイオンでも追い越しが発生すると、その時点以降、検出器８に到達するイオンの順序は質量電荷比の小さい順とはならない。この測定モードを第２測定モードといい、第２測定モードの下で得られた飛行時間スペクトルを追い越し飛行時間スペクトルという。この場合、飛行時間スペクトル上に現れるピークに対応したイオンの周回数は同一とはならず、周回遅れになったイオンの周回数は相対的に少ない。但し、質量分解能は飛行距離、つまりは周回数に依存するから、周回数が多いほど質量分解能は高まる。したがって、正確な質量電荷比を算出するには、非追い越し飛行時間スペクトルではなく追い越し飛行時間スペクトルからマススペクトルを求めることが必要である。そこで、本実施例のＭＴ−ＴＯＦＭＳではそのために特徴的な測定及びデータ処理を実行する。それについて、図２〜図６を参照し具体例を挙げつつ説明する。
【００３５】
本発明に係る質量分析方法を説明する前に、上述した特許文献５で提案した質量分析方法において、導出されるマススペクトルの精度が低下されるケースについて説明する。
ＭＴ−ＴＯＦＭＳにおける飛行時間tofと質量電荷比m/zとの関係は、次の(2)式又は(3)式で表すことができる。また、このとき周回数ｎturnは(4)式で計算することができる。
【数２】
0005533255.tif 000003
ここで、ｕは原子質量[kg]、ｅは電気素量[C]、Ｖ0はイオン加速電圧[Ｖ]、Lin、Lout、Lturnは図１中の各距離[ｍ]、Ｔgは周回軌道５上にあるイオンを離脱させるためにゲート電極２を開放する時間[sec]、を表す。
【００３６】
上記(2)式でｎturn＝１として、周回数１における飛行時間スペクトルをシミュレートした結果を図７に示す。また図８は図７の飛行時間スペクトルの横軸をm/zに変換して作成したマススペクトルである。図７に示した飛行時間スペクトルではイオンの追い越しはない。したがって、図８に示したマススペクトルのピーク強度は正確であるものの質量分解能は低い。図９（ａ）は図７に示した飛行時間スペクトルにおいて最大ピークであるtof＝30.6989[μsec]近辺の拡大図、図９（ｂ）は図８に示したマススペクトルにおける最大ピークであるm/z200近辺の拡大図である。図９（ｂ）から、ピークの半値幅が拡がっていて質量分解能が低いことが分かる。
【００３７】
m/z150であるイオンが100周回した後のタイミングでゲート電極２が開放されて周回軌道５からイオンが出射される場合をシミュレートした飛行時間スペクトルを図１０に示す。この飛行時間スペクトルの縦軸は図７と同じ相対強度であり、本来は100が最大であるにも拘わらず、強度が100を超えているピークが存在する。このことから、この飛行時間スペクトルでは、m/zが異なる複数のイオン由来のピークが重なって（同時に）観測されていることが推測できる。
【００３８】
このとき、上記(4)式を用いてイオンの周回数を計算すると、図１１に示すように、m/z100〜200の質量範囲のイオンは86周回〜126周回していることが分かる。また、上記(2)式、(3)式により質量電荷比m/zと飛行時間tofとの関係を求めると、図１２に示すようになる。即ち、m/z値に応じてtof値が決まる一方、一つのtof値に対して離散的な複数のm/z値が対応していることが分かる。図１１に示したようにm/zが異なるイオン由来のピークが重なった状態である飛行時間スペクトルから、ピークが重なったイオンを分離してそれぞれ周回数を適切に求めることができれば、図１３に示すようなマススペクトルが得られる。図１２及び図１３における最大ピーク付近の拡大図を図１４（ａ）及び（ｂ）に示す。図１４（ｂ）に示すように、この場合にはマススペクトルにおけるピークの半値幅が非常に狭くなっており、質量分解能が高いことが分かる。
【００３９】
特許文献５に提案の方法を用いて、図１０に示した多重周回飛行時間スペクトルを再構成して得られたマススペクトルを図１５に示す。このマススペクトルを図１３に示したマススペクトルと比較すれば分かるように、消滅してしまうピークはないもののピーク強度は正しく推定されていない。この原因は、飛行時間スペクトル上で異なるm/zのイオンが重なっているピークの強度が、マススペクトル上ではその重なっているイオンだけでなく誤って別のイオンに分配されてしまう場合があるためである。一例として、m/z179とm/z138の２種のイオンのピークが重なったtof値＝1049.0[μsec]のピークの強度分配を図１６に示す。図１６中のint.が処理の結果、分配される強度である。この図で分かるように、このピークの強度は本来m/z179とm/z138のみに割り振られるべきであるが、実際にはそれ以外に、m/z148.9712、m/z191.0519など別の質量電荷比のイオンにも割り振られてしまっている。被測定試料に含まれるイオンの種類が多くなるほど、こうした不適切な強度分配が生じる可能性が高くなり、その結果、マススペクトルの強度の精度が下がってしまう。
【００４０】
これに対し、本発明に係る質量分析方法で採用したデータ処理では、被測定試料に含まれるイオンの種類が多い場合であっても、ピークの相対強度を高い精度で推定することができる。図２は本実施例のＭＴ−ＴＯＦＭＳにおける特徴的な分析手順を示すフローチャートである。
【００４１】
未知の被測定試料に対する測定を実行する際には、制御部１０による制御の下で、第１測定モード、つまり異種のイオンの追い越しが発生しないことが保証されている条件の下で被測定試料に対する質量分析を実行し、非追い越し飛行時間スペクトルを取得する（ステップＳ１）。非追い越し飛行時間スペクトルでは、異種のイオンの追いつきや追い越しは生じていないため、飛行速度の大きい順、即ち質量電荷比の小さい順にイオンのピークが並ぶ。
【００４２】
被測定試料に含まれる成分が全く未知の場合には、イオン源１から出発した各種イオンを周回軌道５に沿って全く周回させずに検出器８に導入すればよい。図１の構成では、ゲート電極２に偏向電圧を印加しないことで、各種イオンをそのまま直進通過させ、入射軌道６→出射軌道７に直接導くことにより、周回軌道５をパスすればよい。また、被測定試料に含まれる成分が推測可能である場合で、各種イオンを周回軌道に乗せて１乃至少数回周回させてもイオンの追いつき、追い越しが発生しないことが確実である場合には、そうした周回を許可してもよい。
【００４３】
次に、同じ被測定試料に対し、第２の測定モード、つまり、各種イオンが周回軌道５を多数回周回した後にゲート電極２を経て検出器８に向かうようにゲート電極２へ印加する電圧を制御して質量分析を実行し、追い越し飛行時間スペクトルを取得する（ステップＳ２）。この追い越し飛行時間スペクトルでは、異種のイオンの追い越しが生じているために、質量電荷比の小さい順にピークは並ばず、各ピークの周回数も同一ではない。
【００４４】
いま、上記シミュレーションと同様に、追い越し飛行時間スペクトルが図１０に示すように得られた場合を考える。前述のように、この飛行時間スペクトル上でm/z100〜200であるイオンの周回数は86回〜122回であり、このとき、質量電荷比m/zと周回数ｎturnとの関係は図１１に示すようになる。したがって、図１０に示したように観測される飛行時間スペクトル上の各ピークは、異なる周回数のイオンに由来するピークが加算されたものと考えることができる。図１２に示したtof値とm/z値との関係に基づいて一例を挙げると、tof＝1049.0[μsec]のイオン強度ｙ(1049.0)は、m/zのイオン強度をｘ(m/z)として、次の(5)式で表される。
ｙ(1049.0)＝ｘ(179.0000)＋ｘ(138.0001)＋ｘ(148.9712)＋ｘ(191.0519)＋…    …(5)
同様に、全てのtof値におけるピーク強度ｙ（長さｍのベクトル：ｍはtof値におけるサンプル点数）は、m/zでのピーク強度ｘ（長さｎのベクトル：ｎはm/zにおけるサンプル点数）のいくつかの加算により表される。したがって、次の(6)式による行列式が成り立つ。ここで、Ｘはｘ(m/z)を要素とする行列、Ｙはｙ(tof)を要素とする行列である。
Ｙ＝ＡＸ    …(6)
Ａは係数行列であり、値が１である要素をいくつか含むスパース（疎）な行列である。この係数行列Ａの構造を概念的に示したのが図３である。係数行列Ａの行の数は測定したtofのサンプル数、列の数は換算後のマススペクトル上のm/zのサンプル数で、図３中の塗り潰し部分が値が１である要素を表している。データ処理部９では、追い越し飛行時間スペクトルを構成するスペクトルデータに基づいて、係数行列Ａを算出する（ステップＳ３）。
【００４５】
(6)式を変形すれば次の(7)式となる。即ち、飛行時間スペクトル上で観測されるピーク強度ｙ(tof)にｍ行ｎ列である係数行列Ａの逆行列Ａ-1を乗じれば、目的とする質量電荷比m/zにおけるピーク強度ｘ(m/z)を求めることができる。
Ｘ＝Ａ-1Ｙ    …(7)
したがって、飛行時間スペクトルからマススペクトルを求める作業は、係数行列Ａから逆行列Ａ-1を計算する作業に帰着される。一般に、或る逆行列を求める解法にはいくつかの方法がある。
具体的には、上記(7)式における逆行列Ａ-1を求めるために、次の(8)式に示すムーア・ペンローズ（Moore-Penrose）の一般逆行列（擬似逆行列）Ａ+を利用して、Ｘのノルムを最小とする制約をつけた最小ノルム解を計算する方法がある。
Ａ+＝ＡT （ＡＡT) -1    …(8)
より一般的には、次の(9)式に示す正則化一般逆行列を用い、最小ノルム制約と同時に最小自乗制約を付けることができる。ここで、λは正則化パラメータ、Ｐはペナルティ行列である。
ｍ≦ｎである場合
Ａ+＝ＡT（ＡＡT＋λ2Ｐ）-1    …(9)
ここでＰはｍ行ｍ列
ｍ＞ｎである場合
Ａ+＝（ＡTＡ＋λ2Ｐ）-1ＡT    …(9)
ここでＰはｎ行ｎ列
このときの評価関数は次の(10)式となる。
Ｓ(x)＝｜ｙ−Ａｘ｜2＋λ2ｘTＰｘ    …(10)
上記(8)式は(9)式においてλ＝０とした特異なケースである。
【００４６】
この実施例では平均自乗誤差最小化法のアルゴリズムにより正則化一般逆行列を算出する。即ち、(9)及び(10)式において、λ＝σ（ノイズの標準偏差）、Ｐ＝<ｘ0ｘ0T>−１（＝Ｒxx−１とおく）とすると、推定マススペクトルｘと真のマススペクトルｘ0との平均自乗誤差を最小にすることができる。このときの逆行列と評価関数は次式の通りである。
ｍ≦ｎである場合
Ａ+＝ＡT（ＡＡT＋σ2Ｒxx-1）-1    …(11)
ｍ＞ｎである場合
Ａ+＝（ＲxxＡTＡ＋σ2Ｉ）-1ＲxxＡT    …(11)
Ｓ(x)＝<（ｘ0−ｘ）T（ｘ0−ｘ）>    …(12)
Ｒxxつまり真のマススペクトルｘ0は本来未知であるが、図８に示したような非追い越し飛行時間スペクトルから求まるマススペクトルは質量分解能は低いものの、各ピークの相対強度はほぼ真のマススペクトルを反映していると考えることができる。そこで、ここでは、ステップＳ１で得られた非追い越し飛行時間スペクトルから求まるマススペクトルで真のマススペクトルｘ0を近似する。これにより、(11),(12)式に基づいて、ステップＳ３で求めた係数行列Ａの逆行列を得ることができる（ステップＳ４）。
【００４７】
データ処理部９では、逆行列が求まったならば、(7)式に基づいて各質量電荷比におけるピーク強度を計算し（ステップＳ５）、マススペクトルを作成して表示部１４に表示する（ステップＳ６）。
【００４８】
図１０に示した追い越し飛行時間スペクトルに対し、(11)式に基づいてマススペクトルを推定した結果を図４に示す。ただし、本例では、演算処理のためのコンピュータのメモリの制約によりＲxxを計算することができなかったため、図１０に示した飛行時間スペクトルからピークをピッキングし、それらのピークのみに対するマススペクトル上のピークを推定した。図４に示したマススペクトル上のピークは、図１３に示した正しい推定マススペクトルと比べて、全て一様に強度が小さくなっているものの、その相対強度は正しく求められていることが分かる。実際にマススペクトルを表示する際には、縦軸を相対強度で表示することが一般的であるため、上記の点は実用上問題とはならない。
以上のように、平均自乗誤差最小化法を用いたデータ処理によれば、非追い越し飛行時間スペクトルから求めたマススペクトルと同様のピーク強度を有し、追い越し飛行時間スペクトルで達成される高い時間分解能を活かして高い質量分解能のマススペクトルを作成することができる。
【００４９】
次に、本発明に係る別の実施例に基づく質量分析方法を説明する。この方法では、第１測定モードによる非追い越し飛行時間スペクトルを取得する代わりに、同種のイオンの周回数が異なるように周回軌道５からのイオンの出射タイミング（ゲート電極２を開放する時間）を変えた、２以上の追い越し飛行時間スペクトルを取得する。この場合、２以上の追い越し飛行時間スペクトルを構成するデータから係数行列を算出するが、その際に、点拡がり関数（Point Spread Function）を導入する。
即ち、或る一つのm/zのイオンは飛行中に拡がって、幾つかのtof値で観測されることになる。このイオンの拡がりは点拡がり関数に従い、ここでは、点拡がり関数は標準偏差を質量分解能で表現するガウス関数であると捉える。この場合、飛行時間スペクトル上の或るtof値におけるピークの強度ｙ(tof)は、近傍の点からの拡がりが重なったものである。そこで、点拡がり関数をｈ(ｕ)として、上記(5)式を以下の(13)式のように変形する。
ｙ(1049.0)＝Σｈ(ｕ)×（179.0000−ｕ）＋Σｈ(ｕ)×（138.0001−ｕ）＋Σｈ(ｕ)×（148.9712−ｕ）＋Σｈ(ｕ)×（191.0519−ｕ）＋…      …(13)
したがって、点拡がり関数を導入した係数行列Ａの要素は１だけではなく、ｈ(ｕ)の和が１となる確率密度の値を含んでいる。この場合の係数行列Ａの構造の概念を図５に示す。ここで、図５中の塗り潰しの濃淡は要素の値の大小を表している。
【００５０】
また、１つの追い越し飛行時間スペクトルと１つの（唯一の）マススペクトルとの関係を表す式は追い越し飛行時間スペクトルと同数だけできるから、それらを連立方程式として考える。例えば、いま測定により得られた追い越し飛行時間スペクトルが２つであってそれらをＹ1、Ｙ2とし、それぞれの係数行列をＡ1、Ａ2とすると、行列式(14)により、周回数を補正したマススペクトルＸを求めることができる。
【数３】
0005533255.tif 000004
以上より、イオンが多重周回して測定された追い越し飛行時間スペクトルから、周回数の相違により生じる飛行距離の相違を補正しマススペクトルを求める問題は、(14)式においてＹからＸ（又はｙ(tof）からｘ(m/z)）を求める逆問題と捉えることができる。これは(6)式の関係と同じである。この場合、係数行列は正方行列ではないので、(9)式に示す正則化一般逆行列を用いるのが一般的であり、これを計算することが可能であればマススペクトルを得ることができる。しかしながら、一般には、図５に示したように要素が複雑で巨大な行列から正則化一般逆行列を求めることはハードウエアの制約から困難である。例えば、図１０に示した飛行時間スペクトルをＹとする場合の係数行列Ａを計算すると、238263行194014列の巨大な行列となり、汎用のコンピュータではメモリの制約から逆行列を計算することはできない。
【００５１】
そこで、ここでは点拡がり関数を導入した係数行列に対する逆行列を計算するために、点拡がり関数に対するベイズ的逐次近似解法を利用する。即ち、m/z＝ｉに真のピークが生じる確率をｆ(ｉ)、m/z＝ｉで生じたピークがtof＝ｋで観測される確率をｈ(ｋ｜ｉ)、tof＝ｋで観測されたピークが実はm/z＝ｉで起こった条件付確率をＦ(i｜ｋ)とすると、周知のベイズの定理より(15)式が成り立つ。
Ｆ（i｜ｋ）＝ｆ(ｉ)・ｈ(ｋ｜ｉ)／｛Σjｆ(ｊ)・ｈ(ｋ｜ｊ)｝    …(15)
また、元のピークｆ(ｉ)は、
ｆ(ｉ)＝ΣkＦ(ｉ｜ｋ)・ｈ(ｋ)    …(16)
である。(15)式を(16)式に代入すると、
ｆ(ｉ）＝Σk｛ｈ(ｋ｜ｉ)ｆ(ｉ)・ｈ(ｋ)／｛Σjｆ(ｊ)・ｈ(ｋ｜ｊ)｝    …(17)
となる。(17)式の右辺のｆ(ｉ)をｆr(ｉ)、左辺のｆ(ｉ)をｆr+1(ｉ)とすると、次の漸化式(18)が得られる。
ｆr+1(ｉ）＝ｆr(ｉ)Σk｛ｈ(ｋ｜ｉ)・ｈ(ｋ)／｛Σjｆr(ｊ)・ｈ(ｋ｜ｊ)｝    …(18)
【００５２】
このｆr(ｉ)がｒ→∞のときにｆr(ｉ)→ｆ(ｉ)に収束するとすると、(18)式でｒ→∞として求めたｆr(ｉ)がｆ(ｉ)になると考えることができる。ただし、初期値としてｒ＝０ではｆ0(ｉ)＝１とする。ここでは、(17)式中のｈ(ｋ｜ｉ)をm/z＝ｉのピークがtof＝ｋで観測される点拡がり関数（ガウス関数）、ｈ(ｋ)をtof＝ｋで観測されたピーク強度として、ｒを順に大きくしながらm/z＝ｉにおけるピーク強度ｆ(ｉ)を収束させる。実際には、ｒの数、つまり繰り返し計算回数が所定値になった段階でそのm/z＝ｉに対する計算を打ち切っても十分な精度を確保することができる。
【００５３】
以上説明したベイズ的逐次近似解法で、図７と図１０に示した飛行時間スペクトルの周回数補正を行った結果であるマススペクトルを図６に示す。図６のマススペクトル上のピーク強度はピーク面積として求められるので、図１３に示したマススペクトルとは異なるが、相対強度は正しく求められているので実用上問題ない。
【００５４】
このベイズ的逐次近似解法では上記の平均自乗誤差最小化法とは異なり、質量分解能が低い非追い越し飛行時間スペクトルを必要とせず、質量分解能が相対的に高い追い越し飛行時間スペクトルのみからマススペクトルを求めることができる。このため、処理に供する飛行時間スペクトルの数が少ない場合でも、例えば上記例のように２つであっても、質量分解能が高いマススペクトルを求めることができる。
【００５５】
なお、上記実施例は本発明の一例であり、本発明の趣旨の範囲で適宜変形、修正、追加を行っても本願請求の範囲に包含されることは明らかである。
【符号の説明】
【００５６】
１…イオン源
２…ゲート電極
４…トロイダル扇形電極
５…周回軌道
６…入射軌道
７…出射軌道
８…検出器
９…データ処理部
１０…制御部
１１…入出射電圧発生部
１２…周回飛行電圧発生部
１３…入力部
１４…表示部
【図１】
0005533255.tif 000005
【図２】
0005533255.tif 000006
【図３】
0005533255.tif 000007
【図９】
0005533255.tif 000008
【図１１】
0005533255.tif 000009
【図１４】
0005533255.tif 000010
【図１５】
0005533255.tif 000011
【図１６】
0005533255.tif 000012
【図４】
0005533255.tif 000013
【図５】
0005533255.tif 000014
【図６】
0005533255.tif 000015
【図７】
0005533255.tif 000016
【図８】
0005533255.tif 000017
【図１０】
0005533255.tif 000018
【図１２】
0005533255.tif 000019
【図１３】
0005533255.tif 000020