(57)【要約】      （修正有）【課題】スペクトルデータから任意の２つの分子のピーク間の組成を推定することができる質量分析データ処理装置、質量分析システム及び質量分析データ処理方法を提供する。【解決手段】質量分析データ処理装置は、演算部と、探索部と、を備えている。演算部は、マススペクトルデータから選択された２つの分子のピーク間の質量差Δｍ／ｚを算出、又は質量差Δｍ／ｚを取得する。探索部は、質量差Δｍ／ｚの範囲で、２つの分子のピーク間の原子の組み合わせを推定する。さらに、探索部は、２つの分子のうち一方の分子から減少する原子群と増加する原子群の質量差が質量差Δｍ／ｚと一致する原子の組み合わせを探索する。【選択図】図２
2019002696.tif 000002
【特許請求の範囲】
【請求項１】
マススペクトルデータから選択された２つの分子のピーク間の質量差を算出、又は前記質量差を取得する演算部と、
前記質量差の範囲で、２つの分子のピーク間の原子の組み合わせを推定する探索部と、を備え、
前記探索部は、２つの分子のうち一方の分子から減少する原子群と増加する原子群の質量差が前記質量差と一致する原子の組み合わせを探索する
質量分析データ処理装置。
【請求項２】
前記探索部は、予め設定された原子の種類と、原子の数の最小値と最大値の範囲で原子の組み合わせを探索し、
前記最小値には、０以下の負数を設定可能である
請求項１に記載の質量分析データ処理装置。
【請求項３】
前記探索部が探索した原子の組み合わせに対して絞り込みを行う絞り込み部を備えた
請求項１又は２に記載の質量分析データ処理装置。
【請求項４】
前記絞り込み部は、減少する原子の組み合わせの不飽和度と、増加する原子の組み合わせの不飽和度に基づいて絞り込み処理を行う
請求項３に記載の質量分析データ処理装置。
【請求項５】
試料の質量分析を行い、マススペクトルデータを生成する質量分析計と、
前記質量分析計から前記マススペクトルデータを取得する質量分析データ処理装置と、を備え、
前記質量分析データ処理装置は、
前記マススペクトルデータから選択された２つの分子のピーク間の質量差を算出、又は前記質量差を取得する演算部と、
前記質量差の範囲で、２つの分子のピーク間の原子の組み合わせを推定する探索部と、を備え、
前記探索部は、２つの分子のうち一方の試料から減少する原子群と増加する原子群の質量差が前記質量差と一致する原子の組み合わせを探索する
質量分析システム。
【請求項６】
マススペクトルデータから選択された２つの分子のピーク間の質量差を算出、又は前記質量差を取得する演算工程と、
前記質量差の範囲で、２つの分子のピーク間の原子の組み合わせを推定する探索工程と、を含み、
前記探索工程では、２つの分子のうち一方の分子から減少する原子群と増加する原子群の質量差が前記質量差と一致する原子の組み合わせを探索する
質量分析データ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、導入された試料を解析するために用いられる質量分析データ処理装置、質量分析システム及び質量分析データ処理方法に関する。
【背景技術】
【０００２】
従来から、質量分析データ処理装置では、質量分析計が分析したマススペクトルを用いて種々のデータ処理が行われている。例えば、質量分析データ処理装置では、マススペクトルデータのピーク値から組成が未知の分子の組成の推定が行われている。この組成の推定では、対象となる分子のピーク値（精密質量）から原子の組み合わせを探索することで行われている。そのため、対象となる分子の精密質量の値が大きくなるほど、原子の組み合わせの数が膨大になり、候補の絞り込みが困難となっていた。
【０００３】
また、特許文献１には、マススペクトルで検出された任意のイオンと他のイオンとの質量差を算出し、質量差が所望分子の組み合わせと一致するイオン群から選択する技術が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−９２４９５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に記載された技術は、選択されたイオンとの質量差が所望分子の組み合わせと一致するイオンを探索する技術であり、マススペクトルデータから任意の２つの分子のピーク間の組成の推定は行われていなかった。
【０００６】
本発明の目的は、上記の問題点を考慮し、マススペクトルデータから任意の２つの分子のピーク間の組成を推定することができる質量分析データ処理装置、質量分析システム及び質量分析データ処理方法を提供することにある。
【課題を解決するための手段】
【０００７】
上記課題を解決し、本発明の目的を達成するため、本発明の質量分析データ処理装置は、演算部と、探索部と、を備えている。演算部は、マススペクトルデータから選択された２つの分子のピーク間の質量差を算出、又は質量差を取得する。探索部は、質量差の範囲で、２つの分子のピーク間の原子の組み合わせを推定する。さらに、探索部は、２つの分子のうち一方の分子から減少する原子群と増加する原子群の質量差が質量差と一致する原子の組み合わせを探索する。
【０００８】
また、本発明の質量分析システムは、試料の質量分析を行い、マススペクトルデータを生成する質量分析計と、質量分析計からマススペクトルデータを取得する質量分析データ処理装置と、を備えている。質量分析データ処理装置としては、上述した質量分析データ処理装置が用いられる。
【０００９】
さらに、本発明の質量分析データ処理方法は、下記（１）から（２）に示す工程を含んでいる。
（１）マススペクトルデータから選択された２つの分子のピーク間の質量差を算出、又は質量差を取得する演算工程。
（２）質量差の範囲で、２つの分子のピーク間の原子の組み合わせを推定する探索工程。
また、探索工程では、２つの分子のうち一方の分子から減少する原子群と増加する原子群の質量差が質量差と一致する原子の組み合わせを探索する。
【発明の効果】
【００１０】
本発明の質量分析データ処理装置、質量分析システム及び質量分析データ処理方法によれば、任意の２つの分子のピーク間の組成を推定することができる。
【図面の簡単な説明】
【００１１】
【図１】実施の形態例にかかる質量分析システムを示すブロック図である。
【図２】実施の形態例にかかる質量分析データ処理方法を示すフローチャートである。
【図３】実施の形態例にかかる質量分析データ処理方法と従来例の条件値を示すテーブルである。
【図４】実施の形態例にかかる質量分析データ処理方法において用いられるマススペクトルデータの一例を示すものである。
【図５】実施の形態例にかかる質量分析データ処理方法の概要を示す説明図である。
【図６】実施の形態例にかかる質量分析データ処理方法における絞り込み処理を示すフローチャートである。
【図７】実施の形態例にかかる質量分析データ処理システムの表示装置に表示される表示例である。
【発明を実施するための形態】
【００１２】
以下、本発明の質量分析データ処理装置、質量分析システム及び質量分析データ処理方法の実施の形態例について、図１〜図７を参照して説明する。なお、各図において共通の部材には、同一の符号を付している。また、説明は以下の順序で行うが、本発明は、必ずしも以下の形態に限定されるものではない。
【００１３】
１．質量分析システムの構成
まず、本発明の実施の形態例（以下、「本例」という。）にかかる質量分析システムについて図１を参照して説明する。
図１は、本例の質量分析システムの構成を示すブロック図である。
【００１４】
図１に示す質量分析システム１００は、導入された試料を解析するために用いられるシステムである。図１に示すように、質量分析システム１００は、質量分析計（ＭＳ）１と、質量分析データ処理装置１０と、を備えている。質量分析計１と質量分析データ処理装置１０は、無線、又は有線のネットワーク（ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）、インターネット、専用線等）を介して接続され、互いにデータの送受信が可能である。
【００１５】
質量分析計１は、導入された試料をイオン化し、イオンの質量電荷比（ｍ／ｚ）毎の検出強度を検出し、マススペクトルデータを生成する装置である。この質量分析計１は、導入された試料をイオン化するイオン源と、イオン源で生成されたイオンを質量に応じて分離させる質量分離部と、質量分離部により分離されたイオンを検出する検出部とを備えている。
【００１６】
イオン源によるイオン化法としては、電子イオン化（electron ionization：ＥＩ）法、化学イオン化（chemical ionization：ＣＩ）法、高速原子衝撃（fast atom bombardment：ＦＡＢ）法、エレクトロスプレーイオン化（electrospray ionization：ＥＳＩ）法、大気圧化学イオン化（atmospheric pressure chemical ionization：ＡＰＣＩ）法や、マトリックス支援レーザ脱離イオン化（matrix-assisted laser desorption/ ionization：ＭＡＬＤＩ）法等その他各種のイオン化法が適用されるものである。なお、本例のイオン源のイオン化法としては、ＭＡＬＤＩ法が用いられている。
【００１７】
また、質量分離部としては、磁場型、四重極型、イオントラップ型、フーリエ変換イオンサイクトロン共鳴型、飛行時間型等や、これらを複数組み合わせたもの等その他各種の型の質量分離部が適用されるものである。なお、本例の質量分離部としては、飛行時間型が用いられている。
【００１８】
検出部は、検出したイオンの検出強度をアナログ信号に変換して、質量分析データ処理装置１０の制御部２に送信する。
【００１９】
質量分析データ処理装置１０は、制御部２と、記憶部３と、入力部４と、表示装置５とを備えている。制御部２は、質量分析計１の制御を行うコントロール部２ａと、質量分析計１から質量分析データを取得する取込部２ｂと、データ処理部２ｃと、表示装置５を制御する表示制御部２ｄとを有している。
【００２０】
コントロール部２ａには、入力部４が接続されている。入力部４としては、例えばキーボードやスイッチ等の各種入力手段が適用される。取込部２ｂは、質量分析計１からマススペクトルデータを所得する。そして、取込部２ｂは、取得した質量分析データをデータ処理部２ｃに送信する。
【００２１】
データ処理部２ｃには、取得した質量分析データに対して演算処理を行う。データ処理部２ｃは、演算部１１と、探索部１２と、絞り込み部１３とを有している。演算部１１は、取込部２ｂが取得した質量分析データに対して演算処理を行う。探索部１２は、演算部１１が演算処理を行った情報と、入力部４に入力された情報や、記憶部３に格納された情報に基づいて、組成を推定する。
【００２２】
絞り込み部１３は、探索部１２が探索した組成を、予め設定された条件に基づいて、絞り込み処理を行う。絞り込み部１３は、絞り込まれた組成の候補を表示制御部２ｄや記憶部３に送信する。
【００２３】
また、表示制御部２ｄは、データ処理部２ｃで演算処理されたデータや、取込部２ｂが取得した質量分析データ等を表示装置５に表示させるための処理を行う。
【００２４】
記憶部３には、制御部２から送信された各種データや、組成推定の際に用いられる原子の精密質量が質量電荷比（ｍ／ｚ値）として格納されている。
【００２５】
質量分析データ処理装置１０としては、質量分析計１と一体に設けられた制御装置を適用してもよく、あるいは、外部の携帯情報処理端末や、ＰＣ（パーソナルコンピューター）等を適用してもよい。
【００２６】
２．質量分析システムを用いたデータ処理方法
次に、上述した構成を有する質量分析システム１００を用いたデータ処理方法（組成推定処理）の一例について図２〜図６を参照して説明する。
図２は、データ処理方法を示すフローチャートである。図３は、データ処理方法で用いられる条件値を示すテーブルである。図４は、データ処理方法で用いられるマススペクトルデータである。図５は、本例のデータ処理方法の概要を示す説明図である。
【００２７】
図２に示すように、まず、使用者は、入力部４を介して条件値を入力する（ステップＳ１１）。ステップＳ１１の処理で入力される条件値としては、例えば、図３に示すテーブルの、「質量（ｍ／ｚ）の誤差」、絞り込み処理で用いられる「ＤＢＥ値の範囲」、組み合わせの探索（組成推定）で用いられる原子の種類と、その原子の数の「最小値」及び「最大値」である。また、従来のデータ処理方法では、絞り込み処理で用いられる「電子の数」、「電荷」も入力される。入力部４に入力された条件値は、記憶部３に格納される。
【００２８】
本例では、条件値を使用者が入力部４を介して入力する例を説明したが、これに限定されるものではない。例えば、予め図３に示す条件値を記憶部３に格納し、データ処理を行う際に、制御部２のデータ処理部２ｃが、記憶部３から取得してもよい。
【００２９】
次に、制御部２の取込部２ｂは、質量分析計１からマススペクトルデータを取得する（ステップＳ１２）。なお、ステップＳ１２におけるマススペクトルデータを取得する工程は、ステップＳ１１の条件値を入力する工程の前に行ってもよい。
【００３０】
次に、マススペクトルデータから、２つの分子のピーク（モノアイソートピーク）を選択する（ステップＳ１３）。ステップＳ１３の処理は、例えば、表示装置５に図４に示すマススペクトルを表示し、使用者に差の組成推定を行う２つの分子のピークＡ、Ｂを選択させることで行われる。または、差の組成推定を行う２つの分子のピークの質量電荷比（ｍ／ｚ値）又は精密質量を入力部４に入力してもよい。
【００３１】
次に、データ処理部２ｃは、選択された２つのピークＡ、Ｂの質量電荷比（ｍ／ｚ値）、すなわち任意の２つの分子の精密質量を抽出する。そして、データ処理部２ｃの演算部１１は、２つの分子の精密質量の差（質量差）を演算する（ステップＳ１４）。具体的には、演算部１１は、選択された２つのピークＡ、Ｂにおけるｍ／ｚ値の差であるΔｍ／ｚを算出する（ステップＳ１４）。なお、演算部１１は、取込部２ｂや入力部４を介して、質量差Δｍ／ｚを取得してもよい。そして、演算部１１で算出された質量差Δｍ／ｚは、記憶部３に格納されると共に探索部１２に出力される。
【００３２】
次に、探索部１２は、ステップＳ１１で入力された条件値と、記憶部３に格納された組み合わせ探索に用いられる原子の精密質量に基づいて、質量差Δｍ／ｚの範囲内で、原子の組み合わせ（以下、単に「組」という。）を探索する（ステップＳ１５）。
【００３３】
図５に示すように、ピークＡとピークＢの間では、増加する原子だけではなく、減少する原子も存在する。そして、選択したピークＢから減少した原子群（負数の原子群）の精密質量（質量電荷比）を減少質量Ｄ１とし、ピークＢから増加した原子群（正数の原子群）の精密質量（質量電荷比）を増加質量Ａ１とした場合、Ａ１−Ｄ１＝Δｍ／ｚとなる。したがって、探索部１２は、減少質量Ｄ１と増加質量Ａ１の差が、質量差Δｍ／ｚと一致する原子の組を探索する。
【００３４】
ステップＳ１５の探索処理で用いられる原子の種類と、原子の数の最大値と最小値は、ステップＳ１１で入力した図３に示す条件値である。
【００３５】
なお、従来のデータ処理方法では、０から組成推定を行う分子のピークＡまでのｍ／ｚ値の範囲Ｈ１内で原子の組の探索処理が行われていた。そのため、図３に示すように、従来例のデータ処理方法で用いられる条件値では、組み合わせの探索で用いられる原子の数は、０以下の負数を含まない０以上の整数が入力され、その「最小値」は、０となっている。さらに、従来のデータ処理方法では、図５に示すように、探索を行う範囲Ｈ１が広いため、組成推定の誤差も大きくなっていた。
【００３６】
これに対して、本例のデータ処理では、２つのピークＡ、Ｂにおける質量差Δｍ／ｚの範囲内で原子の組の探索処理が行われている。上述したように、本例のデータ処理の探索部１２は、減少する原子も考慮して、原子の組の探索処理を行っている。これにより、任意の２つの分子のピーク間の差の組成を推定することができる。
【００３７】
また、減少する原子も考慮するため、図３に示すように、ステップＳ１１の処理で入力される条件値には、原子の数の「最小値」に０以下の数として、負数も入力可能になっている。本例では、Ｃの数の範囲が−６〜＋１００（最小値−６、最大値＋１００）、Ｈの数の範囲が−１５〜＋２００（最小値−１５、最大値＋２００）、Ｏの数の範囲が−２〜＋１５（最小値−２、最大値＋１５）、Ｎの数の範囲が−４〜＋１０（最小値−４、最大値＋１０）、Ｐの数の範囲が０〜＋２（最小値０、最大値＋２）に設定されている。これらの原子の精密質量は、予め記憶部３に格納されている。なお、原子の種類及び原子の数の範囲は、組成推定する対象に応じて任意に変更可能である。
【００３８】
さらに、本例のデータ処理方法では、原子の組の探索処理の範囲が、２つのピークＡ、Ｂ間であるため、探索を行う範囲を従来のデータ処理方法の範囲Ｈ１よりも狭めることができ、組成推定の誤差を小さくすることができる。さらに、組成推定を行う範囲が狭まるため、原子の組み合わせの数を減らすことができ、候補の絞り込みを容易に行うことができる。
【００３９】
例えば、チロシンの組成推定を行う場合、チロシンの精密質量、すなわち実測したｍ／ｚ値は、１８２．０８１１７である。そのため、従来例の組成推定を行う範囲のｍ／ｚ値は、１８２．０８１１７となる。
【００４０】
これに対して、チロシンとアルギニンの差の組成推定を行う本例のデータ処理の場合、ピークＡとしてチロシンのピーク値が選択され、ピークＢとしてアルギニンのピーク値が選択される。ここで、アルギニンのチロシンの精密質量、すなわち実測したｍ／ｚ値は、１７５．１１８８９である。そして、チロシンとアルギニンの質量差Δｍ／ｚは、６．９６２２８となる。したがって、本例でのデータ処理方法で組成推定を行う範囲のｍ／ｚ値は、６．９６２２８となり、上述した従来の組成推定で行う範囲よりも狭まったことが分かる。
【００４１】
そして、探索部１２は、探索した原子の組を、表示制御部２ｄ、絞り込み部１３や、記憶部３等に出力する。なお、探索部１２が出力する原子の組は、負数の原子群と正数の原子群が混在した原子の組である。例えば、探索部１２からＣ３Ｈ−３Ｏ１Ｎ−３という原子の組が出力された場合、Ｃが３つ、Ｏが１つ増加し、Ｈが３つ、Ｎが３つ減少したことを示す。
【００４２】
次に、データ処理部２ｃは、探索した原子の組を、減少した原子群である負数の原子群と、増加した原子群である正数の原子群に分ける（ステップＳ１６）。ステップＳ１５の処理において、例えば、Ｃ３Ｈ−３Ｏ１Ｎ−３という原子の組が探索された場合、負数の原子群としてＨ３とＮ３、正数の原子群としてＣ３とＯ１に分けられる。
【００４３】
そして、ステップＳ１６で分けられた、負数の原子群と、正数の原子群に基づいて、絞り込み部１３は、絞り込み処理を行う（ステップＳ１７）。なお、本例では、負数の原子群と、正数の原子群に分けた例を説明したが、これに限定されるものではなく、負数の原子群と正数の原子群が混在した原子の組を用いて絞り込み処理を行ってもよい。
【００４４】
［絞り込み処理］
ここで、図６を参照してステップＳ１７の絞り込み処理の一例について説明する。
図６は、絞り込み処理の一例を示すフローチャートである。
【００４５】
図６に示すように、絞り込み部１３は、ステップＳ１６の処理で分けられた負数の原子群のＤＢＥ（Ｄｏｕｂｌｅ  Ｂｏｎｄ  Ｅｑｕｉｖａｌｅｎｔ）値と、正数の原子群のＤＢＥ値を算出する（ステップＳ２１）。すなわち、ステップＳ２１の処理では、負数の原子群と正数の原子群のそれぞれの不飽和度を算出する。
【００４６】
ＤＢＥ値は、原子の組において、炭素原子の数をｗ、水素原子の数をｘ、ハロゲン原子の数をｙ、窒素原子の数をｚとすると、下記式１から算出される。なお、酸素や硫黄のような第１６族元素の数は、カウントしない。ただし、ヒドロキシ基のような酸素原子に結合した水素元素は、カウントする。
［式１］
ＤＢＥ＝（２＋２ｗ−ｘ−ｙ＋ｚ）／２
【００４７】
例えば、Ｃ３Ｈ−３Ｏ１Ｎ−３という原子の組が探索された場合、負数の原子群としてＨ３とＮ３、正数の原子群としてＣ３とＯ１に分けられるため、負数の原子群のＤＢＥ値は、（２−３＋３）／２から１となり、正数の原子群のＤＢＥ値は、（２＋２×３）／２から４となる。
【００４８】
次に、絞り込み部１３は、ステップＳ２１の処理で算出された負数の原子群のＤＢＥ値と、正数の原子群のＤＢＥ値の差を算出する（ステップＳ２２）。そして、負数の原子群のＤＢＥ値が１、正数の原子群のＤＢＥ値が４の場合、ＤＢＥ値の差は、４−１から３となる。
【００４９】
次に、絞り込み部１３は、ステップＳ２２の処理で算出されたＤＢＥ値の差が、ステップＳ１１の処理で設定した範囲内であるか否かを判断する（ステップＳ２３）。ここで用いられる範囲とは、図３に示す「ＤＢＥの範囲」であり、本例では、−１．５〜２５．０に設定されている。なお、ＤＢＥの範囲は、組成推定する対象によって任意に変更可能である。
【００５０】
そして、ステップＳ２３の処理において、ＤＢＥ値の差が範囲外であると判断した場合（ステップＳ２３のＮＯ判定）、絞り込み部１３は、絞り込みを行った原子の組を組成の候補から排除する（ステップＳ２６）。そして、絞り込み処理が完了し、図２に示すステップＳ１８の処理に戻る。
【００５１】
これに対して、ステップＳ２３の処理において、ＤＢＥ値の差が範囲内であると判断した場合（ステップＳ２３のＹＥＳ判定）、絞り込み部１３は、負数の原子群のＤＢＥ値と正数の原子群のＤＢＥ値のｏｄｄ／ｅｖｅｎが同一であるか否かを判断する（ステップＳ２４）。
【００５２】
ここで、ｅｖｅｎ（偶数）では、全ての原子の手（原子価）が結合されており、ｏｄｄ（奇数）では、原子に結合されていない手を有することを示している。なお、ＤＢＥ値は、式１に示すように２で除算されている。そのため、ｅｖｅｎは整数となり、ｏｄｄ（奇数）は整数に０．５を加えた値となる。
【００５３】
上述したように、Ｃ３Ｈ−３Ｏ１Ｎ−３という原子の組では、負数の原子群のＤＢＥ値は１で、正数の原子群のＤＢＥ値は４であるため、負数の原子群のＤＢＥ値及び正数の原子群のＤＢＥ値は、それぞれ整数である。そのため、減少した原子群の原子価の状態と、増加した原子群の原子価の状態が一致していると判断できる。
【００５４】
また、本例のステップＳ２４の処理では、負数の原子群のＤＢＥ値と正数の原子群のＤＢＥ値のｏｄｄ／ｅｖｅｎがそれぞれ同一であるか否かを判断した例を説明したが、これに限定されるものではない。ここで、負数の原子群のＤＢＥ値と正数の原子群のＤＥＢ値のｏｄｄ／ｅｖｅｎが同一の場合、正数の原子群のＤＢＥ値から負数の原子群のＤＢＥ値を減算すると必ず整数となる。逆に、ｏｄｄ／ｅｖｅｎが異なる場合には、必ず半整数、すなわち整数に０．５を加えた値となる。そのため、正数の原子群のＤＢＥ値から負数の原子群のＤＢＥ値を減算することで、ステップＳ２４の処理と同一の判定を行うことができる。
【００５５】
ステップＳ２４の処理において、負数の原子群のＤＢＥ値と正数の原子群のＤＢＥ値のｏｄｄ／ｅｖｅｎが同一でないと判断した場合（ステップＳ２４のＮＯ判定）、絞り込み部１３は、絞り込みを行った原子の組を組成の候補から排除する（ステップＳ２６）。そして、絞り込み処理が完了し、図２に示すステップＳ１８の処理に戻る。
【００５６】
これに対して、ステップＳ２４の処理において、負数の原子群のＤＢＥ値と正数の原子群のＤＢＥ値のｏｄｄ／ｅｖｅｎが同一であると判断した場合（ステップＳ２４のＹＥＳ判定）、絞り込み部１３は、絞り込みを行った原子の組を組成の候補に登録する（ステップＳ２５）。すなわち、絞り込み部１３は、絞り込みを行った原子の組を記憶部３に格納する。これにより、絞り込み部１３による絞り込み処理が完了し、図２に示すステップＳ１８の処理に戻る。
【００５７】
次に、図２に示すように、データ処理部２ｃは、探索部１２によって質量差Δｍ／ｚの範囲における全ての原子の組の探索が終了したか否かを判断する（ステップＳ１８）。ステップＳ１８の処理において、データ処理部２ｃがまだ探索が終了していないと判断した場合（ステップＳ１８のＮＯ判定）、ステップＳ１５の処理に戻り、原子の組の探索処理を続ける。
【００５８】
また、ステップＳ１８の処理において、データ処理部２ｃが全ての原子の組の探索が終了したと判断した場合（ステップＳ１８のＹＥＳ判定）、表示制御部２ｄは、ステップＳ１７の処理で絞り込まれた原子の組の候補を、表示装置５に表示させる。これにより、質量分析システム１００を用いたデータ処理が終了する。
【００５９】
上述したように、本例のデータ処理方法によれば、組成が未知の２つの分子のピーク間の組成を容易に推定することができる。これに対して、従来のデータ処理方法では、組成が未知の２つの分子のピーク間の組成を推定するためには、２つの分子の組成をそれぞれ推定する必要があった。従来のデータ処理方法では、さらに推定された２つの分子の組成の組み合わせから、２つの分子のピーク間の組成を推定する必要があり、組成の候補が膨大となるだけではなく、候補の絞り込みが困難なものとなっていた。
【００６０】
さらに、本例のデータ処理方法では、２つの分子の内、一方の組成が既知の場合、既知の分子の組成に、上述したデータ処理方法で取得した２つの分子のピーク間の組成の候補を組み合わせることで、残りの未知の分子の組成を推定することができる。また、組成推定を行う範囲が狭まるため、原子の組み合わせの数を減らすことができ、候補の絞り込みを容易に行うことができる。
【００６１】
図７は、表示装置５に表示される組成推定処理の結果画面の表示例を示す説明図である。
図７に示すように、表示装置５の表示画面５０には、精密質量差表示領域５１と、条件値表示領域５２と、組成候補リスト表示領域５３が表示される。精密質量差表示領域５１には、上述したデータ処理方法のステップＳ１４の処理において演算部１１の演算工程で演算された２つの分子の質量差Δｍ／ｚが表示される。
【００６２】
条件値表示領域５２には、上述したデータ処理方法のステップＳ１１の処理で設定された条件値が表示される。条件値表示領域５２には、「質量（ｍ／ｚ）の誤差」が表示される分析誤差表示領域５２ａと、絞り込み処理で用いられる「ＤＢＥ値の範囲」が表示されるＤＢＥ値表示領域５２ｂが表示されている。さらに、条件値表示領域５２には、原子の組み合わせの探索に用いられる原子の種類と、その原子の数の「最小値（Ｍｉｎ）」と、「最大値（Ｍａｘ）」を表示する絞り込み原子リスト５２ｃが表示される。
【００６３】
組成候補リスト表示領域５３には、上述したデータ処理方法の組成正推定処理の結果が表示される。組成候補リスト表示領域５３には、組成候補原子リスト５４と、組成候補精密質量リスト５５と、組成候補ＤＢＥ値リスト５６と、推定質量差表示領域５７と、誤差表示領域５８が表示される。
【００６４】
組成候補原子リスト５４には、上述したデータ処理方法で組成推定された原子の組み合わせが表示される。組成候補原子リスト５４は、組成推定された原子の組み合わせのうち減少した原子群が表示される減少原子群リスト５４ａと、組成推定された原子の組み合わせのうち増加した原子群が表示される増加原子群リスト５４ｂとを有している。
【００６５】
組成候補精密質量リスト５５は、減少原子群精密質量表示領域５５ａと、増加原子群精密質量表示領域５５ｂとを有している。減少原子群精密質量表示領域５５ａには、減少した原子群の精密質量が表示される。増加原子群精密質量表示領域５５ｂには、増加した原子群の精密質量が表示される。
【００６６】
組成候補ＤＢＥ値リスト５６には、上述した絞り込み処理におけるステップＳ２２の処理で算出されたＤＢＥ値が表示される。組成候補ＤＢＥ値リスト５６は、減少原子群ＤＢＥ値表示領域５６ａと、増加原子群ＤＢＥ値表示領域５６ｂとを有している。減少原子群ＤＢＥ値表示領域５６ａには、減少した原子群のＤＢＥ値、すなわち負数の原子群のＤＢＥ値が表示される。増加原子群ＤＢＥ値表示領域５６ｂには、増加した原子群のＤＢＥ値、すなわち正数の原子群のＤＢＥ値が表示される。
【００６７】
推定質量差表示領域５７は、組成推定された原子の組み合わせにおける増加した原子群の精密質量と、減少した原子群の精密質量の差（以下、「推定質量差」という）が表示される。この推定質量差は、演算部１１又は探索部１２で演算される。また、誤差表示領域５８には、推定質量差表示領域５７に表示された推定質量差における、演算部１１で演算された質量差Δｍ／ｚに対する誤差が表示される。誤差としては、質量差Δｍ／ｚと推定質量差との差分や、質量差Δｍ／ｚに対する推定質量差の比率等が用いられる。そして、誤差は、演算部１１又は探索部１２で演算される。
【００６８】
なお、表示装置５の表示画面５０に表示される表示例としては、上述した図７に示すものに限定されるものではなく、例えば、組成候補リスト表示領域５３のみが表示されていてもよい。さらに、組成候補リスト表示領域５３には、組成候補原子リスト５４と、誤差表示領域５８のみが表示されていてもよい。表示装置５の表示画面５０には、その他各種の表示例が表示されるものである。さらに、絞り込み処理を行わず、探索部１２で探索した全ての原子の組み合わせを処理結果として表示装置５に表示させてもよい。
【００６９】
なお、本発明は上述しかつ図面に示した実施例に限定されるものではなく、特許請求の範囲に記載した発明の要旨を逸脱しない範囲内で種々の変形実施が可能である。
【００７０】
上述した実施の形態例では、ＤＢＥ値を用いた絞り込み処理を適用した例を説明したが、絞り込み処理はこれに限定されるものではない。絞り込み処理としては、例えば、水素と炭素の比率から絞り込みを行ってもよく、あるいは減少した原子の数と、増加した原子の数の比率から絞り込みを行ってもよくその他各種の絞り込み方法を適用できるものである。
【符号の説明】
【００７１】
１…質量分析計、  ２…制御部、  ２ａ…コントロール部、  ２ｂ…取込部、  ２ｃ…データ処理部、  ２ｄ…表示制御部、  ３…記憶部、  ４…入力部、  ５…表示装置、  １０…質量分析データ処理装置、  １１…演算部、  １２…探索部、  １３…絞り込み部、５０…表示画面、  １００…質量分析システム、  Δｍ／ｚ…質量差
【図１】
2019002696.tif 000003
【図２】
2019002696.tif 000004
【図３】
2019002696.tif 000005
【図４】
2019002696.tif 000006
【図５】
2019002696.tif 000007
【図６】
2019002696.tif 000008
【図７】
2019002696.tif 000009