特開 2024-18233 測定装置および測定方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024018233

(43)【公開日】2024-02-08

(54)【発明の名称】測定装置および測定方法

(51)【国際特許分類】

   G01N 24/00 20060101AFI20240201BHJP

   G01R 33/26 20060101ALI20240201BHJP

   G01R 33/20 20060101ALI20240201BHJP

【ＦＩ】

G01N24/00 P

G01R33/26

G01R33/20

【審査請求】未請求

【請求項の数】5

【出願形態】ＯＬ

(21)【出願番号】P 2022121430

(22)【出願日】2022-07-29

【国等の委託研究の成果に係る記載事項】（出願人による申告）平成３０年度、国立研究開発法人科学技術振興機構、産学共創プラットフォーム共同研究推進プログラム委託事業、「超スマート社会実現のカギを握る革新的半導体技術を基盤としたエネルギーイノベーションの創出に関する国立大学法人京都大学による研究開発」、産業技術力強化法第１７条の適用を受ける特許出願

(71)【出願人】

【識別番号】000107804

【氏名又は名称】スミダコーポレーション株式会社

(71)【出願人】

【識別番号】504132272

【氏名又は名称】国立大学法人京都大学

(74)【代理人】

【識別番号】100114971

【弁理士】

【氏名又は名称】青木 修

(72)【発明者】

【氏名】芳井 義治

(72)【発明者】

【氏名】佐光 暁史

(72)【発明者】

【氏名】竹村 祐輝

(72)【発明者】

【氏名】水落 憲和

(57)【要約】

【課題】 被測定場の測定値におけるノイズ成分を抑制して測定精度を高くする。
【解決手段】 受光装置１３は、励起光に対応して磁気共鳴部材１から発せられる蛍光を受光し蛍光強度に対応する蛍光センサー信号を生成し、演算処理装置３１は、蛍光センサー信号または前記蛍光センサー信号から得られる検出信号に基づいて測定値を導出する。演算処理装置３１は、被測定場が磁気共鳴部材１に印加されているときの測定値を本測定値とし、被測定場が磁気共鳴部材１に印加されていないときであって本測定値の測定前に測定された測定値を前段測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定後に測定された測定値を後段測定値としたとき、前段測定値および後段測定値をそれぞれ所定の比率で本測定値から減算して、被測定場の測定値を導出する。
【選択図】 図１

【特許請求の範囲】

【請求項1】

被測定場に対応して電子スピン量子状態が変化するとともにマイクロ波で電子スピン量子操作の可能な磁気共鳴部材と、
前記マイクロ波で前記磁気共鳴部材の電子スピン量子操作を行う高周波磁場発生器と、
前記磁気共鳴部材に照射すべき励起光を出射する発光装置と、
前記磁気共鳴部材により前記励起光に対応して発せられる蛍光を受光し前記蛍光の強度に対応する蛍光センサー信号を生成する蛍光受光装置と、
前記蛍光センサー信号または前記蛍光センサー信号から得られる検出信号に基づいて測定値を導出する演算処理装置とを備え、
前記演算処理装置は、前記被測定場が前記磁気共鳴部材に印加されているときの前記測定値を本測定値とし、前記被測定場が前記磁気共鳴部材に印加されていないときであって前記本測定値の測定前に測定された前記測定値を前段測定値とし、前記被測定場が前記磁気共鳴部材に印加されていないときであって前記本測定値の測定後に測定された前記測定値を後段測定値としたとき、前記前段測定値および前記後段測定値をそれぞれ所定の比率で前記本測定値から減算して、前記被測定場の測定値を導出すること、
を特徴とする測定装置。

【請求項2】

被測定場に対応して電子スピン量子状態が変化するとともにマイクロ波で電子スピン量子操作の可能な磁気共鳴部材に対して、所定の測定シーケンスに従って、前記マイクロ波で前記磁気共鳴部材の電子スピン量子操作を行うとともに、前記磁気共鳴部材に照射すべき励起光を出射し、
前記磁気共鳴部材により前記励起光に対応して発せられる蛍光を受光し前記蛍光の強度に対応する蛍光センサー信号を生成し、
前記蛍光センサー信号または前記蛍光センサー信号から得られる検出信号に基づいて測定値を導出し、
前記被測定場が前記磁気共鳴部材に印加されているときの前記測定値を本測定値とし、前記被測定場が前記磁気共鳴部材に印加されていないときであって前記本測定値の測定前に測定された前記測定値を前段測定値とし、前記被測定場が前記磁気共鳴部材に印加されていないときであって前記本測定値の測定後に測定された前記測定値を後段測定値としたとき、前記前段測定値および前記後段測定値をそれぞれ所定の比率で前記本測定値から減算して、前記被測定場の測定値を導出すること、
を特徴とする測定方法。

【請求項3】

前記蛍光センサー信号または前記検出信号に対して、所定の窓関数を適用し、前記窓関数を適用した後の前記蛍光センサー信号または前記検出信号に基づいて前記測定値を導出し、
前記窓関数は、前記励起光の照射期間の前半における第１期間と後半における第２期間との間での、前記蛍光センサー信号または前記検出信号の積算値の差分を前記測定値として導出するためのものであり、
前記前段測定値および前記後段測定値の測定の際にそれぞれ、前記比率に対応する重み係数を乗じた前記窓関数を適用して前記前段測定値および前記後段測定値を導出し、実質的に前記前段測定値および前記後段測定値をそれぞれ前記比率で前記本測定値から減算するようにして、前記本測定値並びに前記前段測定値および前記後段測定値から前記被測定場の測定値を導出すること、
を特徴とする請求項２記載の測定方法。

【請求項4】

前記蛍光センサー信号に対して、前記励起光を分岐して得られる参照光を受光して生成される参照光センサー信号に基づくコモンモードリジェクションを行い、前記コモンモードリジェクションに基づいて前記検出信号を生成することを特徴とする請求項２または請求項３記載の測定方法。

【請求項5】

（ａ）前記被測定場が前記磁気共鳴部材に印加されているときの前記測定値と、前記被測定場が前記磁気共鳴部材に印加されていないときの前記測定値とを交互に取得し、（ｂ）前記後段測定値を、次の前記本測定値についての前記前段測定値として使用して、連続的に前記被測定場の測定値を導出することを特徴とする請求項２記載の測定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、測定装置および測定方法に関するものである。

【背景技術】

【0002】

ある磁場測定装置は、窒素と格子欠陥（ＮＶセンター：Nitrogen Vacancy Center）を有するダイヤモンド構造などといったセンシング部材の電子スピン共鳴を利用した光検出磁気共鳴（ＯＤＭＲ：Optically Detected Magnetic Resonance）で磁気計測を行っている。ＯＤＭＲでは、このようなＮＶセンターを有するダイヤモンドといった磁気共鳴部材に対して、被測定磁場とは別に静磁場が印加されるとともに、所定のシーケンスでレーザー光（初期化および測定のための励起光）並びにマイクロ波が印加され、その磁気共鳴部材から出射する蛍光の光量が検出されその光量に基づいて被測定磁場の磁束密度が導出される。

【0003】

例えば、ラムゼイパルスシーケンスでは、（ａ）励起光をＮＶセンターに照射し、（ｂ）マイクロ波の第１のπ／２パルスをＮＶセンターに印加し、（ｃ）第１のπ／２パルスから所定の時間間隔ｔｔでマイクロ波の第２のπ／２パルスをＮＶセンターに印加し、（ｄ）励起光をＮＶセンターに照射してＮＶセンターの発光量を測定し、（ｅ）測定した発光量に基づいて磁束密度を導出する。また、スピンエコーパルスシーケンスでは、（ａ）励起光をＮＶセンターに照射し、（ｂ）マイクロ波の第１のπ／２パルスを被測定磁場の位相０度でＮＶセンターに印加し、（ｃ）マイクロ波のπパルスを被測定磁場の位相１８０度でＮＶセンターに印加し、（ｄ）マイクロ波の第２のπ／２パルスを被測定磁場の位相３６０度でＮＶセンターに印加し、（ｅ）励起光をＮＶセンターに照射してＮＶセンターの発光量を測定し、（ｆ）測定した発光量に基づいて磁束密度を導出する。

【0004】

あるセンサー装置は、上述のようなＮＶセンターを含むダイヤモンドセンサーを使用した核磁気共鳴で磁場測定を行っている（例えば特許文献１参照）。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開２０１９－１３８７７２号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

センシング部材の電子スピン共鳴を利用した光検出磁気共鳴においては、蛍光から得られる検出信号が微弱であるため、ノイズの影響を受けやすく測定精度が低くなってしまう。

【0007】

本発明は、上記の問題に鑑みてなされたものであり、ノイズ成分を抑制して測定精度を高くする測定装置および測定方法を得ることを目的とする。

【課題を解決するための手段】

【0008】

本発明に係る測定装置は、被測定場に対応して電子スピン量子状態が変化するとともにマイクロ波で電子スピン量子操作の可能な磁気共鳴部材と、マイクロ波で磁気共鳴部材の電子スピン量子操作を行う高周波磁場発生器と、磁気共鳴部材に照射すべき励起光を出射する発光装置と、磁気共鳴部材により励起光に対応して発せられる蛍光を受光し蛍光の強度に対応する蛍光センサー信号を生成する蛍光受光装置と、蛍光センサー信号または蛍光センサー信号から得られる検出信号に基づいて測定値を導出する演算処理装置とを備える。そして、演算処理装置は、被測定場が磁気共鳴部材に印加されているときの測定値を本測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定前に測定された測定値を前段測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定後に測定された測定値を後段測定値としたとき、前段測定値および後段測定値をそれぞれ所定の比率で本測定値から減算して、被測定場の測定値を導出する。

【0009】

本発明に係る測定方法は、（ａ）被測定場に対応して電子スピン量子状態が変化するとともにマイクロ波で電子スピン量子操作の可能な磁気共鳴部材に対して、所定の測定シーケンスに従って、マイクロ波で磁気共鳴部材の電子スピン量子操作を行うとともに、磁気共鳴部材に照射すべき励起光を出射し、（ｂ）磁気共鳴部材により励起光に対応して発せられる蛍光を受光し蛍光の強度に対応する蛍光センサー信号を生成し、蛍光センサー信号または蛍光センサー信号から得られる検出信号に基づいて測定値を導出する。そして、被測定場が磁気共鳴部材に印加されているときの測定値を本測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定前に測定された測定値を前段測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定後に測定された測定値を後段測定値としたとき、前段測定値および後段測定値をそれぞれ所定の比率で本測定値から減算して、被測定場の測定値を導出する。

【発明の効果】

【0010】

本発明によれば、ノイズ成分を抑制して測定精度を高くする測定装置および測定方法が得られる。

【図面の簡単な説明】

【0011】

【図1】図１は、本発明の実施の形態１に係る測定装置の構成を示すブロック図である。

【図2】図２は、励起光の光量に対する蛍光センサー信号のレベルの非線形性について説明する図である。

【図3】図３は、測定シーケンスの一例を示す図である。

【図4】図４は、参照光センサー信号ｒｅｆ１（ｔ），ｒｅｆ２（ｔ）、蛍光センサー信号ＰＬ（ｔ）、およびＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）について説明する図である。

【図5】図５は、窓関数の一例を示す図である。

【図6】図６は、図５に示す窓関数の周波数特性を示す図である。

【図7】図７は、実施の形態１に係る測定装置の動作について説明するフローチャートである。

【図8】図８は、測定値の取得について説明するフローチャートである。

【図9】図９は、本発明の実施の形態２に係る測定装置の構成を示すブロック図である。

【図10】図１０は、蛍光センサー信号および窓関数の一例を示す図である。

【図11】図１１は、実施の形態３における連続的な被測定場の測定について説明する図である。

【発明を実施するための形態】

【0012】

以下、図に基づいて本発明の実施の形態を説明する。

【0013】

実施の形態１．

【0014】

図１は、本発明の実施の形態１に係る測定装置の構成を示すブロック図である。図１に示す測定装置は、センサー部１０と、高周波電源１１と、発光装置１２と、受光装置１３とを備える。

【0015】

センサー部１０は、所定の位置（例えば、検査対象物体の表面上または表面上方）において、被測定場（例えば磁場の強度、向きなどといった磁場）を検出する。なお、被測定場は、単一周波数の交流場でもよいし、複数の周波数成分を有する所定周期の交流場でもよいし、直流場でもよい。

【0016】

この実施の形態では、センサー部１０は、磁気共鳴部材１、高周波磁場発生器２、および磁石３を備え、ＯＤＭＲで被測定場を検出する。

【0017】

磁気共鳴部材１は、結晶構造を有し、被測定場（ここでは、磁場）に対応して電子スピン量子状態が変化するとともに、結晶格子における欠陥および不純物の配列方向に応じた周波数のマイクロ波で（ラビ振動に基づく）電子スピン量子操作の可能な部材である。つまり、磁場の測定位置に、磁気共鳴部材１が配置される。

【0018】

この実施の形態では、磁気共鳴部材１は、複数（つまり、アンサンブル）の特定カラーセンターを有する光検出磁気共鳴部材である。この特定カラーセンターは、ゼーマン分裂可能なエネルギー準位を有し、かつ、ゼーマン分裂時のエネルギー準位のシフト幅が互いに異なる複数の向きを取り得る。

【0019】

ここでは、磁気共鳴部材１は、単一種別の特定カラーセンターとして複数のＮＶ（Nitrogen Vacancy）センターを含むダイヤモンドなどの部材である。ＮＶセンターの場合、基底状態がｍｓ＝０，＋１，－１の三重項状態であり、ｍｓ＝＋１の準位およびｍｓ＝－１の準位がゼーマン分裂する。ＮＶセンターが、ｍｓ＝＋１およびｍｓ＝－１の準位の励起状態から基底状態へ遷移する際に、所定の割合で蛍光を伴い、残りの割合のＮＶセンターは、励起状態（ｍｓ＝＋１またはｍｓ＝－１）から基底状態（ｍｓ＝０）へ無輻射で遷移する。

【0020】

なお、磁気共鳴部材１に含まれるカラーセンターは、ＮＶセンター以外のカラーセンターでもよい。

【0021】

高周波磁場発生器２は、マイクロ波を磁気共鳴部材１に印加して、磁気共鳴部材１の電子スピン量子操作を行う。例えば、高周波磁場発生器２は、板状コイルであって、マイクロ波を放出する略円形状のコイル部と、そのコイル部の両端から延び基板に固定される端子部とを備える。高周波電源１１は、そのマイクロ波の電流を生成して高周波磁場発生器２に導通させる。そのコイル部は、その両端面部分において、磁気共鳴部材１を挟むように所定の間隔で互いに平行な２つの電流を導通させ、上述のマイクロ波を放出する。ここでは、コイル部は板状コイルであるが、表皮効果により、コイル部の端面部分をマイクロ波の電流が流れるため、２つの電流が形成される。これにより、空間的に均一な強度のマイクロ波が磁気共鳴部材１に印加される。

【0022】

ＮＶセンターの場合、ダイヤモンド結晶において、欠陥（空孔）（Ｖ）および不純物としての窒素（Ｎ）によってカラーセンターが形成されており、ダイヤモンド結晶内の欠陥（空孔）（Ｖ）に対して、隣接する窒素（Ｎ）の取り得る位置（つまり空孔と窒素との対の配列方向）は４種類あり、それらの配列方向のそれぞれに対応するゼーマン分裂後のサブ準位（つまり、基底からのエネルギー準位）が互いに異なる。したがって、マイクロ波の周波数に対する静磁場によるゼーマン分裂後の蛍光強度の特性において、それぞれの向きｉ（ｉ＝１，２，３，４）に対応して、互いに異なる４つのディップ周波数対（ｆｉ＋，ｆｉ－）が現れる。ここでは、この４つのディップ周波数対のうちのいずれかのディップ周波数に対応して、上述のマイクロ波の周波数（波長）が設定される。

【0023】

また、磁石３は、磁気共鳴部材１に静磁場（直流磁場）を印加し、磁気共鳴部材１内の複数の特定カラーセンター（ここでは、複数のＮＶセンター）のエネルギー準位をゼーマン分裂させる。ここでは、磁石３は、リング型の永久磁石であり、例えば、フェライト磁石、アルニコ磁石、サマコバ磁石などである。

【0024】

この実施の形態では、上述の静磁場の印加方向は、上述の被測定磁場の印加方向と同一となり、上述の静磁場の印加によって、上述のディップ周波数での蛍光強度変化が増強され、感度が高くなる。

【0025】

さらに、この実施の形態では、磁気共鳴部材１は、上述のマイクロ波で電子スピン量子操作の可能な複数のカラーセンター（ここでは、ＮＶセンター）を備え、磁石３は、磁気共鳴部材１の所定領域（励起光の照射領域）に対して略均一な静磁場を印加する。例えば、その所定領域における静磁場の強度についての最大値と最低値との差分や比率が所定値以下となるように静磁場が印加される。

【0026】

また、磁気共鳴部材１において、上述の欠陥および不純物の配列方向が、上述の静磁場の向き（および印加磁場の向き）に略一致するように、磁気共鳴部材１の結晶が形成され、磁気共鳴部材１の向きが設定される。

【0027】

さらに、この実施の形態では、励起光を磁気共鳴部材１に照射するために、発光装置１２から磁気共鳴部材１までの光学系が設けられており、また、磁気共鳴部材１からの蛍光を検出するために、磁気共鳴部材１から受光装置１３までの光学系が設けられている。

【0028】

発光装置１２は、光源としてのレーザーダイオードなどを備え、その光源で、磁気共鳴部材１に照射すべき励起光として、所定波長のレーザー光を出射する。また、受光装置１３は、受光素子としてのフォトダイオードやフォトトランジスターなどを備え、磁気共鳴部材１により励起光に対応して発せられる蛍光を受光し、その蛍光の強度に対応する蛍光センサー信号ＰＬを生成する。この蛍光は、例えば複合放物面型集光器（ＣＰＣ）などの光学系によって受光装置１３へ向けて集光される。

【0029】

ここで、測定原理について説明する。

【0030】

上述の励起光の強度Ｉは、次式のように、本来の強度Ｉ_{ｌａｓｅｒ}とノイズ成分の強度Ｉ_{ｎｏｉｓｅ}との和となっている。なお、このノイズ成分は、発光装置１２の電源電圧のふらつきや光源の発光量のふらつきなどに起因して発生し、例えばｋＨｚオーダー程度から１００ｋＨｚオーダー程度の範囲の周波数を有する。

【0031】

Ｉ＝Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}

【0032】

また、蛍光センサー信号ＰＬのレベルは、基本的には励起光の強度Ｉが大きいほど高くなる。ただし、被測定場による電子スピン量子状態の変化に起因して、測定時の励起光照射開始時点の蛍光強度が小さくなっており、その後、被測定場による電子スピン量子状態の変化の影響がなくなるまで、蛍光強度が徐々に大きくなっていく。そのため、その被測定場による電子スピン量子状態の変化に起因する蛍光センサー信号ＰＬのレベル変動分α_ｃｏｎｔ（ｔ）に比例する検出信号が、被測定場を示す信号として導出される。蛍光センサー信号ＰＬは、次式で表される。ここで、当該蛍光強度は比較的低いため、蛍光センサー信号ＰＬのレベルは、非線形性を有する。図２は、励起光の光量に対する蛍光センサー信号のレベルの非線形性について説明する図である。例えば図２に示すように、蛍光センサー信号ＰＬのレベルの傾きは、蛍光強度が小さいほど、大きくなっている。

【0033】

ＰＬ（ｔ）＝α（ｔ）×ｆ（Ｉ）＝（α_ｉｎｉｔ＋α_ｃｏｎｔ（ｔ））×Ｉ＝（α_ｉｎｉｔ＋α_ｃｏｎｔ（ｔ））×ｆ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）

【0034】

ここで、α_ｉｎｉｔは、励起光の強度Ｉに対して比例する部分（被測定場による電子スピン量子状態の変化に影響を受けない部分）を示す係数であり、ｆは、励起光の強度Ｉに対する非線形性を示す関数である。例えば、ｆは、Ｎ次式（Ｎ≧２）であり、この実施の形態では、ｆは２次式（ｆ（ｘ）＝ｘ＋γ・ｘ^２，γは定数）である。この場合の蛍光センサー信号ＰＬは、次式で表される。

【0035】

なお、この関数ｆや定数γは、受光装置１３のセンサーについての実験などにより予め導出される。また、この関数ｆは、指数関数などといった他の関数形を有していてもよい。

【0036】

ＰＬ（ｔ）＝α（ｔ）×ｆ（Ｉ）＝（α_ｉｎｉｔ＋α_ｃｏｎｔ（ｔ））×Ｉ＝（α_ｉｎｉｔ＋α_ｃｏｎｔ（ｔ））×｛（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＋γ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）^２｝＝α_ｉｎｉｔ ｛（γＩ_{ｌａｓｅｒ} ^２＋Ｉ_{ｌａｓｅｒ}）＋（２γＩ_{ｌａｓｅｒ}＋１）Ｉ_{ｎｏｉｓｅ}＋γＩ_{ｎｏｉｓｅ} ^２｝＋α_ｃｏｎｔ（ｔ）｛（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＋γ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）^２｝

【0037】

他方、励起光から分岐した参照光の参照光センサー信号のレベルｒｅｆは、励起光の強度Ｉに比例するため、次式のように表される。ここで、当該参照光の強度は比較的高いため、参照光センサー信号ＰＬのレベルは、参照光強度に対して線形（比例）となる。

【0038】

ｒｅｆ＝ｒｅｆ１＝ｒｅｆ２＝β×Ｉ＝β×（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）

【0039】

ここで、βは、定数である。なお、ｒｅｆ１，ｒｅｆ２については後述する。

【0040】

具体的には、ＰＬは、テーラー展開し、Ｉ_{ｎｏｉｓｅ}＜＜Ｉ_{ｌａｓｅｒ}であるため３次項以降を無視することで、上述の式で示される。また、以下のように、上述のＰＬおよびｒｅｆをＩ_{ｎｏｉｓｅ}でテーラー展開（マクローリン展開）し、３次項以降を無視して近似することで、そのＰＬおよびｒｅｆからＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）が導出される。なお、以下の式において、ＰＬ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）は、ＰＬが（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）の関数であることを示し、ＰＬ（Ｉ_{ｌａｓｅｒ}）は、Ｉ_{ｎｏｉｓｅ}を０としたときのＰＬである。同様に、ｒｅｆ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）は、ｒｅｆが（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）の関数であることを示し、ｒｅｆ（Ｉ_{ｌａｓｅｒ}）は、Ｉ_{ｎｏｉｓｅ}を０としたときのｒｅｆである。なお、ＰＬ’は、Ｉ_{ｎｏｉｓｅ}＝０での、Ｉ_{ｎｏｉｓｅ}についてのＰＬの１次微分係数であり、ｒｅｆ’は、Ｉ_{ｎｏｉｓｅ}＝０での、Ｉ_{ｎｏｉｓｅ}についてのｒｅｆの１次微分係数である。

【0041】

ＰＬ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＝ＰＬ（Ｉ_{ｌａｓｅｒ}）＋ＰＬ’（Ｉ_{ｌａｓｅｒ}）Ｉ_{ｎｏｉｓｅ}＋・・・

【0042】

ｒｅｆ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＝ｒｅｆ（Ｉ_{ｌａｓｅｒ}）＋ｒｅｆ’（Ｉ_{ｌａｓｅｒ}）Ｉ_{ｎｏｉｓｅ}＋・・・

【0043】

ここで、Ｉ_{ｎｏｉｓｅ}＜＜Ｉ_{ｌａｓｅｒ}であるため、Ｉ_{ｎｏｉｓｅ}の３次項以降は無視できる。そのため、上述のＰＬおよびｒｅｆは、以下の式に近似される。

【0044】

ＰＬ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＝ＰＬ（Ｉ_{ｌａｓｅｒ}）＋ＰＬ’（Ｉ_{ｌａｓｅｒ}）Ｉ_{ｎｏｉｓｅ}

【0045】

ｒｅｆ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＝ｒｅｆ（Ｉ_{ｌａｓｅｒ}）＋ｒｅｆ’（Ｉ_{ｌａｓｅｒ}）Ｉ_{ｎｏｉｓｅ}

【0046】

ここで、ＰＬ（ｔ）に対してコモンモードリジェクションが実行され、ＣＭＲ信号が生成される。具体的には、ＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）は、α_ｉｎｉｔについてのＩ_{ｎｏｉｓｅ}の影響が除去されるように、次式で導出される。

【0047】

ＣＭＲ＿ＳＩＧ（ｔ）＝ＰＬ－（ＰＬ’（Ｉ_{ｌａｓｅｒ}）／ｒｅｆ’（Ｉ_{ｌａｓｅｒ}））ｒｅｆ

【0048】

ここで、ＰＬ’（Ｉ_{ｌａｓｅｒ}）／ｒｅｆ’（Ｉ_{ｌａｓｅｒ}）＝α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／βであるので、ＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）は、次式のようになる。なお、次式では、時間ｔに拘わらず定数である項（後述）を除去するために、定数α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２が加算されている。

【0049】

ＣＭＲ＿ＳＩＧ（ｔ）＝ＰＬ（ｔ）－α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／β×ｒｅｆ（ｔ）＋α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２＝α_ｃｏｎｔ（ｔ）｛（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＋γ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）^２｝

【0050】

ここで、ＰＬ（ｔ）における、α_ｉｎｉｔ｛（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＋γ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）^２｝については、十分小さい項γＩ_{ｎｏｉｓｅ} ^２を無視し、α_ｉｎｉｔ｛（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）＋γ（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）^２｝＝α_ｉｎｉｔ（γＩ_{ｌａｓｅｒ}＋１）Ｉ_{ｌａｓｅｒ}＋α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）Ｉ_{ｎｏｉｓｅ}＝α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）×（Ｉ_{ｌａｓｅｒ}＋Ｉ_{ｎｏｉｓｅ}）－α_ｉｎｉｔγＩ_{ｌａｓｅｒ}としている。

【0051】

また、レベル変動成分信号ＳＤ（ｔ）は、次式のように導出される。これにより、α_ｃｏｎｔ（ｔ）についてのＩ_{ｎｏｉｓｅ}の影響が除去される。

【0052】

ＳＤ（ｔ）＝α_ｃｏｎｔ（ｔ）＝ＣＭＲ＿ＳＩＧ（ｔ）／（ｒｅｆ（ｔ）＋γ・ｒｅｆ（ｔ）^２）

【0053】

上述のＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）またはレベル変動成分信号ＳＤ（ｔ）が検出信号（被測定場を示す信号）として導出され、その検出信号のピーク値（ｔ＝０の値）、時間積分値、後述の積算値の差分などが被測定場の強度と相関があるため、検出信号のピーク値、時間積分値、後述の積算値の差分などと被測定場の強度との対応関係を予め実験などで特定しておき、その対応関係を示す計算式やテーブルを使用して、検出信号のピーク値、時間積分値、後述の積算値の差分などから被測定場の強度が導出される。

【0054】

このような測定原理に基づき、以下に述べる構成が設けられている。

【0055】

図１に示す測定装置は、さらに、発光装置１２から磁気共鳴部材１までの励起光の光路上に、光学素子としての光分離部２１，２２を備える。光分離部２１，２２は、それぞれ、励起光の一部を、励起光から分岐させ、参照光として別の方向に出射する。例えば、光分離部２１，２２は、偏向無依存型のビームスプリッターである。

【0056】

また、図１に示す測定装置は、参照光を受光し、その参照光の強度に対応する参照光センサー信号ｒｅｆ１，ｒｅｆ２（上述のｒｅｆ）を生成する受光装置２３，２４を備える。

【0057】

この実施の形態では、励起光から２つの参照光が別々に生成され、２つの参照光センサー信号ｒｅｆ１，ｒｅｆ２が生成される。参照光センサー信号ｒｅｆ１は、後述のコモンモードリジェクションに使用され、参照光センサー信号ｒｅｆ２は、デジタイズされ、検出信号ＳＤの生成に使用される。

【0058】

さらに、図１に示す測定装置は、アナログ演算回路としてのＣＭＲ演算部２５を備える。ＣＭＲ演算部２５は、蛍光センサー信号ＰＬに対して、参照光センサー信号ｒｅｆ１に基づくコモンモードリジェクションを行い、そのコモンモードリジェクションに基づくＣＭＲ信号ＣＭＲ＿ＳＩＧを生成する。具体的には、ＣＭＲ演算部２５は、係数部２５ａと、オフセット除去部２５ｂと、差動アンプ２５ｃとを備える。係数部２５ａは、所定係数α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／βを参照光センサー信号ｒｅｆ１に乗算する。オフセット除去部２５ｂは、上述のオフセット分を示す定数α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２を係数部２５ａの出力信号ｒｅｆ１×α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／βから減算する。差動アンプ２５ｃは、蛍光センサー信号ＰＬと、オフセット除去部２５ｂの出力信号ｒｅｆ１×α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／β－α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２との差分を演算し、その演算結果をＣＭＲ信号ＣＭＲ＿ＳＩＧとして出力する。なお、図１では、オフセット除去部２５ｂは、上述のオフセット分を示す定数α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２を係数部２５ａの出力信号ｒｅｆ１×α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／βから減算しているが、その代わりに、上述のオフセット分を示す定数α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２を蛍光センサー信号ＰＬに加算してもよいし、あるいは、上述のオフセット分を示す定数α_ｉｎｉｔγＩ_{ｌａｓｅｒ} ^２を差動アンプ２５ｃの出力信号に加算してもよい。

【0059】

なお、アナログ演算回路としての係数部２５ａを設けてもよいし、係数部２５ａを設けずに、受光装置２３のゲインを調整して、所定係数α_ｉｎｉｔ（２γＩ_{ｌａｓｅｒ}＋１）／βを乗算した参照光センサー信号を出力するようにしてもよい。

【0060】

さらに、図１に示す測定装置は、ＣＭＲ信号ＣＭＲ＿ＳＩＧおよび参照センサー信号ｒｅｆ２をそれぞれデジタイズするアナログデジタル変換器２６，２７、および測定装置の制御および信号処理を行う演算処理装置３１を備える。

【0061】

アナログデジタル変換器２６，２７は、所定のビット数かつ所定のサンプリング周期（速度）で、ＣＭＲ信号ＣＭＲ＿ＳＩＧおよび参照センサー信号ｒｅｆ２をそれぞれデジタイズし、デジタイズしたＣＭＲ信号ＣＭＲ＿ＳＩＧおよび参照センサー信号ｒｅｆ２を演算処理装置３１に出力する。

【0062】

演算処理装置３１は、例えばコンピューターを備え、信号処理プログラムをコンピューターで実行して、各種処理部として動作する。この実施の形態では、演算処理装置３１は、そのコンピューターを測定制御部４１および演算部４２として動作させ、また、不揮発性の記憶装置４３を備える。

【0063】

記憶装置４３には、信号処理プログラムが記憶されており、そのコンピューターは、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）などを備え、信号処理プログラをＲＡＭにロードしてＣＰＵで実行することで、測定制御部４１および演算部４２として動作する。

【0064】

測定制御部４１は、所定の測定シーケンスに従って、（ａ）高周波電源１１および発光装置１２を制御し、（ｂ）上述のようにデジタイズされたＣＭＲ信号ＣＭＲ＿ＳＩＧおよび参照光センサー信号ｒｅｆ２を取得してＲＡＭや記憶装置４３に記憶し、演算部４２に、被測定場の測定値を導出させる。

【0065】

この測定シーケンスは、被測定場の周波数などに従って設定される。例えば、被測定場が比較的高周波数の交流場である場合には、この測定シーケンスには、スピンエコーパルスシーケンス（ハーンエコーシーケンスなど）が適用される。ただし、測定シーケンスは、これに限定されるものではない。また、例えば、被測定場が比較的低周波数の交流場である場合、被測定場の１周期において、複数回、ラムゼイパルスシーケンス（つまり、直流場の測定シーケンス）で、その物理場の測定を行い、それらの測定結果に基づいて、被測定場（強度、波形など）を特定するようにしてもよい。

【0066】

図３は、測定シーケンスの一例を示す図である。図３は、スピンエコーパルスシーケンスの場合の、被測定磁場に対するマイクロ波パルスのタイミング、および励起光の照射タイミング（初期化と測定の２回）を示している。図３に示すように、各測定について、励起光の照射期間において、蛍光が検出される。

【0067】

図４は、参照光センサー信号ｒｅｆ１（ｔ），ｒｅｆ２（ｔ）、蛍光センサー信号ＰＬ（ｔ）、およびＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）について説明する図である。

【0068】

図４に示すように、参照光センサー信号ｒｅｆ１（ｔ），ｒｅｆ２（ｔ）は、照射期間の略矩形のパルス信号となり、蛍光センサー信号ＰＬ（ｔ）は、照射期間において徐々に立ち上がり一定レベルに収束するパルス信号となる。そして、コモンモードリジェクションによって、ＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）が得られる。

【0069】

ここでは、演算部４２は、上述のように、デジタイズされたＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）およびデジタイズされた参照光センサー信号ｒｅｆ２（ｔ）から、α_ｃｏｎｔ（ｔ）に比例するレベル変動成分信号ＳＤ（ｔ）を計算して生成し、その信号ＳＤ（ｔ）に基づいて被測定場の測定値（ここでは、磁束密度や磁場の波形など）を導出する。あるいは、演算部４２は、上述のように、デジタイズされたＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）を検出信号とし、その検出信号（ＣＭＲ信号）に基づいて被測定場の測定値を導出してもよい。

【0070】

具体的には、演算部４２は、被測定場が磁気共鳴部材１に印加されているときの測定値を本測定値とし、被測定場が磁気共鳴部材１に印加されていないときであって本測定値の測定前に測定された測定値を前段測定値とし、被測定場が磁気共鳴部材１に印加されていないときであって本測定値の測定後に測定された測定値を後段測定値としたとき、前段測定値および後段測定値をそれぞれ所定の比率で本測定値から減算して、被測定場の測定値を導出する。つまり、同一の測定方法で、本測定値、前段測定値、および後段測定値が取得され、前段測定値および後段測定値で本測定値が補正され、補正後の本測定値が、その時点での被測定場の測定値とされる。

【0071】

例えば、前段測定値Ｓａ、本測定値Ｓｏ、および後段測定値Ｓｂの重み係数を－０．５、＋１．０、－０．５として、前段測定値、本測定値、および後段測定値にそれぞれ重み係数を乗じたものの合計を、被測定場の測定値Ｓ（Ｓ＝－０．５×Ｓａ＋Ｓｏ－０．５×Ｓｂ）とする。ここでは、前段測定値Ｓａの重み係数と後段測定値Ｓｂの重み係数とは互いに同一である。なお、この比率は、前段測定値および後段測定値の重み係数と本測定値の重み係数との間で正負が反転しており、かつ、前段測定値、本測定値、および後段測定値の重み係数の合計がゼロであれば、他の比率でもよい。

【0072】

さらに、実施の形態１では、演算部４２は、検出信号（ＣＭＲ＿ＳＩＧ（ｔ）またはＳＤ（ｔ））に対して、所定の窓関数を適用し、その窓関数を適用した後の蛍光センサー信号ＰＬ（ｔ）または検出信号（ＣＭＲ＿ＳＩＧ（ｔ）またはＳＤ（ｔ））に基づいて各測定値（前段測定値、本測定値、および後段測定値）を導出する。

【0073】

図５は、窓関数の一例を示す図である。図５に示すように、上述の窓関数は、期間Ｐ１と期間Ｐ２とで正負が反転し、絶対値が同一である波形として表される。ここでは、窓関数は、励起光の照射期間の前半における第１期間Ｐ１と後半における第２期間Ｐ２との間での、検出信号（ＣＭＲ＿ＳＩＧ（ｔ）またはＳＤ（ｔ））の積算値の差分を導出するためのものである。

【0074】

また、図６は、図５に示す窓関数の周波数特性を示す図である。検出信号がＣＭＲ信号である場合には、演算部４２は、（ａ）励起光の照射期間の前半部分における期間Ｐ１および後半部分における期間Ｐ２において、それぞれ、互いに同一複数回得られる（例えば１０００回などといった所定サンプリング数の）、デジタイズされたＣＭＲ信号の値を積算し、（ｂ）期間Ｐ１についてのＣＭＲ信号の積算値（総和や平均）と期間Ｐ２についてのＣＭＲ信号の積算値（総和や平均）との差分を計算してＣＭＲ信号におけるノイズ成分を除去する。例えば図４に示すように、前半部分における期間Ｐ１は、照射開始時刻（ｔ＝０）から所定時間長の期間であり、後半部分における期間Ｐ２は、照射終了時刻（ｔ＝ｔｅ）まで所定時間長の期間である。ここでは、Ｐ１とＰ２の時間長は同一である。そして、前半部分の積算値から後半部分の積算値を減算して、その減算結果の値を、窓関数適用後のＣＭＲ信号の値とすることで、ＣＭＲ信号のノイズ成分（約１０ｋＨｚ以上といった高周波ノイズ成分）が抑制される。なお、検出信号がＳＤ信号である場合も同様である。

【0075】

この実施の形態では、演算部４２は、例えば図５に示すように、前段測定値および後段測定値の測定の際にそれぞれ、その比率に対応する重み係数を乗じた窓関数を適用して前段測定値および後段測定値を導出し、実質的に前段測定値および後段測定値をそれぞれ上述の比率で本測定値から減算するようにして、本測定値並びに前段測定値および後段測定値から被測定場の測定値を導出する。

【0076】

さらに、図５に示すように、前段測定値、本測定値、および後段測定値の重み係数を－０．５、＋１．０、－０．５とした場合の、一連の窓関数の周波数特性は、例えば図６に示すような周波数特性となり、数ｋＨｚ以下のノイズ成分（１／ｆノイズなど）が抑制される。なお、図６における破線（比較例）は、前段測定値および後段測定値を適用しない場合の周波数特性を示しており、前段測定値および後段測定値を適用した場合、前段測定値および後段測定値を適用しない場合より数ｋＨｚ以下のノイズ成分（１／ｆノイズなど）の減衰率が大きくなっている。

【0077】

また、この実施の形態では、アナログデジタル変換器２６は、アナログデジタル変換器２７に比べ高速に動作し、アナログデジタル変換器２７は、アナログデジタル変換器２６に比べ高精度のデジタイズを行う。

【0078】

例えば、アナログデジタル変換器２６は、例えば２００Ｍサンプル／秒で、入力アナログ信号を２０ビットのデジタル信号へ変換し、アナログデジタル変換器２７は、例えば１００ｋサンプル／秒で、入力アナログ信号を２４ビットのデジタル信号へ変換する。

【0079】

なお、アナログデジタル変換器２６によりアナログデジタル変換されたＣＭＲ信号は、比較的早く変化するため、高速なアナログデジタル変換器２６でサンプリングされ、上述のように、複数サンプリングに基づくノイズ除去が行われる。他方、上述のように、アナログデジタル変換器２７によりアナログデジタル変換された参照光センサー信号ｒｅｆ２は、ＣＭＲ信号における、励起光のノイズ成分Ｉ_{ｎｏｉｓｅ}の影響を理論上消すための演算に使用され、検出される蛍光センサー信号ＰＬにおけるレベル変動分α_ｃｏｎｔ（ｔ）とノイズ成分強度Ｉ_{ｎｏｉｓｅ}との積の項の、蛍光センサー信号ＰＬの電圧レベルに対する相対的な電圧レベルに応じた精度でその演算をする必要があるため、比較的高精度のアナログデジタル変換器２７でサンプリングされる。

【0080】

次に、当該実施の形態に係る測定装置の動作について説明する。図７は、実施の形態１に係る測定装置の動作（つまり、測定方法）について説明するフローチャートである。

【0081】

被測定場（ここでは、交流磁場Ｂａｃ）の測定位置にセンサー部１０が配置される。なお、センサー部１０を走査しつつ複数の測定位置でそれぞれ測定を行うようにしてもよい。

【0082】

その後、まず、測定制御部４１は、被測定場がオフされた状態（つまり、被測定場が磁気共鳴部材１に印加されていない状態）で測定を行い、前段測定値を取得する（ステップＳ１）。

【0083】

ここで、測定値の取得について説明する。図８は、測定値の取得について説明するフローチャートである。

【0084】

測定制御部４１は、所定の測定シーケンスに従って、発光装置１２に励起光を発光させたり、高周波磁場発生器２にマイクロ波を送出させたりする（ステップＳ１１）。

【0085】

これにより、測定時の励起光の照射期間において、受光装置１３から蛍光センサー信号ＰＬ（アナログ信号）が出力され、受光装置２３，２４から参照光センサー信号ｒｅｆ１，ｒｅｆ２（アナログ信号）がそれぞれ出力される（ステップＳ１２）。そして、ＣＭＲ演算部２５により、蛍光センサー信号ＰＬ（ｔ）および参照光センサー信号ｒｅｆ１からＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）が出力されアナログデジタル変換器２６によりデジタイズされる（ステップＳ１３）。

【0086】

測定制御部４１が、このＣＭＲ信号ＣＭＲ＿ＳＩＧ（ｔ）（デジタル信号）を取得すると、演算部４２は、ＣＭＲ信号に対して上述の窓関数を適用し、窓関数適用後のＣＭＲ信号の積算値（例えば、上述の期間Ｐ１での積算値と期間Ｐ２での積算値との差分）を上述の測定値として導出し、図示せぬメモリーなどに記憶する（ステップＳ１４）。

【0087】

このようにして、測定値（ここでは、前段測定値、ただし、本測定値および後段測定値についても同様）が取得される。なお、検出信号として、ＳＤ（ｔ）が使用される場合には、上述のように、参照光センサー信号ｒｅｆ２およびＣＭＲ信号に基づいて、ＳＤ（ｔ）が導出され、ＳＤ（ｔ）に対して上述の窓関数が適用され、窓関数適用後のＳＤ（ｔ）の積算値が上述の測定値として導出される。

【0088】

次に、測定制御部４１は、被測定場がオンされた状態（つまり、被測定場が磁気共鳴部材１に印加されている状態）で上述の測定を行い、本測定値を取得する（ステップＳ２）。

【0089】

その後、測定制御部４１は、被測定場がオフされた状態で上述の測定を行い、後段測定値を取得する（ステップＳ３）。

【0090】

そして、演算部４２は、後段測定値の取得後、本測定値、並びに、本測定値の直前および直後の前段測定値および後段測定値に基づいて、被測定場の測定値を導出する（ステップＳ４）。具体的には、被測定場が磁気共鳴部材１に印加されていないため、前段測定値および後段測定値は、測定系に由来するノイズ成分（１／ｆノイズなど）を示す。そのため、本測定値から、特定の比率で前段測定値および後段測定値を減算することで、本測定値における、それらのノイズ成分が抑制される。

【0091】

以上のように、上記実施の形態１によれば、高周波磁場発生器２は、被測定場に対応して電子スピン量子状態が変化するとともにマイクロ波で電子スピン量子操作の可能な磁気共鳴部材１の電子スピン量子操作をマイクロ波で行う。発光装置１２は、磁気共鳴部材１に照射すべき励起光を出射する。受光装置１３は、磁気共鳴部材１によりその励起光に対応して発せられる蛍光を受光しその蛍光の強度に対応する蛍光センサー信号を生成する。演算処理装置３１は、蛍光センサー信号から得られる検出信号に基づいて測定値を導出する。演算処理装置３１は、被測定場が磁気共鳴部材１に印加されているときの測定値を本測定値とし、被測定場が磁気共鳴部材１に印加されていないときであって本測定値の測定前に測定された測定値を前段測定値とし、被測定場が磁気共鳴部材に印加されていないときであって本測定値の測定後に測定された測定値を後段測定値としたとき、前段測定値および後段測定値をそれぞれ所定の比率で本測定値から減算して、被測定場の測定値を導出する。

【0092】

これにより、本測定値の直前および直後の両方においてノイズ成分の測定値として前段測定値および後段測定値を導出し、前段測定値および後段測定値によって本測定値を補正しているので、ノイズ成分が抑制され測定精度が高くなる。

【0093】

実施の形態２．

【0094】

図９は、本発明の実施の形態２に係る測定装置の構成を示すブロック図である。図１０は、蛍光センサー信号および窓関数の一例を示す図である。

【0095】

実施の形態２では、コモンモードリジェクションを行わずに、アナログデジタル変換器２６によって蛍光センサー信号ＰＬ（ｔ）がデジタイズされ、演算部４２は、蛍光センサー信号ＰＬ（ｔ）に対して上述の窓関数を適用し、その窓関数を適用した後の蛍光センサー信号ＰＬ（ｔ）に基づいて各測定値（前段測定値、本測定値、および後段測定値）を導出する。ここでは、窓関数は、励起光の照射期間の前半における第１期間Ｐ１と後半における第２期間Ｐ２との間での、蛍光センサー信号ＰＬ（ｔ）の積算値の差分を導出するためのものである。なお、上述の検出信号と同様に、窓関数適用後の蛍光センサー信号ＰＬ（ｔ）の積算値の差分は被測定場の強度と相関があるため、両者の対応関係を予め実験などで特定しておき、その対応関係を示す計算式やテーブルを使用して、窓関数適用後の蛍光センサー信号ＰＬ（ｔ）の積算値の差分から被測定場の強度が導出される。

【0096】

そして、演算部４２は、後段測定値の取得後、本測定値、並びに、本測定値の直前および直後の前段測定値および後段測定値に基づいて、被測定場の測定値を導出する。

【0097】

なお、実施の形態２に係る測定装置についてのその他の構成および動作については、実施の形態１と同様であるので、その説明を省略する。

【0098】

実施の形態３．

【0099】

図１１は、実施の形態３における連続的な被測定場の測定について説明する図である。実施の形態３では、例えば図１１に示すように、（ａ）被測定場が磁気共鳴部材１に印加されているときの測定値と、被測定場が磁気共鳴部材１に印加されていないときの測定値とを（ここでは、所定の時間間隔で繰り返し）交互に取得し、（ｂ）（第ｉ回の測定）の後段測定値Ｓｂ（ｉ）を、次の（つまり、第（ｉ＋１）回の測定の）本測定値Ｓｏ（ｉ＋１）についての前段測定値Ｓａ（ｉ＋１）として使用して、連続的に被測定場の測定値Ｓ（ｉ），Ｓ（ｉ＋１），・・・を導出する。なお、各測定値Ｓａ，Ｓｂ，Ｓｏの取得（つまり、測定）については、実施の形態１または２と同様である。

【0100】

したがって、（第ｉ回の測定についての）被測定場が磁気共鳴部材１に印加されていないときの測定値は、後段測定値Ｓｂ（ｉ）および前段測定値Ｓａ（ｉ＋１）として使用されるため、少なくとも、第（ｉ＋１）回の測定が完了するまで図示せぬメモリーで保持される。

【0101】

なお、実施の形態３に係る測定装置についてのその他の構成および動作については、実施の形態１または２と同様であるので、その説明を省略する。

【0102】

以上のように、上記実施の形態３によれば、後段測定値Ｓｂ（ｉ）を、次の本測定値Ｓｏ（ｉ＋１）についての前段測定値Ｓａ（ｉ＋１）として使用するため、連続的に複数回の測定を行う際の所要時間が短くなる。すなわち、所定時間あたりの測定回数が多くなり、測定値を平均化した場合、測定値のランダムノイズを（測定回数の平方根に比例して）低減させることができ、測定感度が向上する。

【0103】

なお、上述の実施の形態に対する様々な変更および修正については、当業者には明らかである。そのような変更および修正は、その主題の趣旨および範囲から離れることなく、かつ、意図された利点を弱めることなく行われてもよい。つまり、そのような変更および修正が請求の範囲に含まれることを意図している。

【0104】

例えば、上記実施の形態では、光検出磁気共鳴に基づき磁場測定が行われるが、温度測定なども同様に実行可能である。また、磁場測定により得られる磁場に基づく電流測定を行うことも可能である。

【0105】

また、上記実施の形態では、コモンモードリジェクション用の参照光とは別に、上述の参照光センサー信号を使用した検出信号の演算のための参照光を励起光から分岐させているが、コモンモードリジェクション用に分岐した参照光を、上述の参照光センサー信号を使用した検出信号の演算のための参照光としても使用してもよい。また、上述の参照光センサー信号ｒｅｆ１，ｒｅｆ２を別々に生成しているが、１つの参照光センサー信号ｒｅｆを生成し、その参照光センサー信号ｒｅｆを上述の参照光センサー信号ｒｅｆ１，ｒｅｆ２として使用するようにしてもよい。

【0106】

また、上記実施の形態では、検出信号としてＣＭＲ信号またはレベル変動成分信号が使用されるが、それらとは異なる信号（被測定場を示す信号）を検出信号として使用してもよい。

【0107】

また、上記実施の形態において、差動アンプ２５ｂは、アナログデジタル変換器２６とは別に設けられていてもよいし、アナログデジタル変換器２６に内蔵されていてもよい。

【0108】

また、上記実施の形態では、窓関数は、矩形状の波形（つまり、期間Ｐ１，Ｐ２のそれぞれにおいて一定な値）を有しているが、他の形状の波形でもよい。つまり、窓関数は、期間Ｐ１，Ｐ２のそれぞれにおいて変化する値を取ってもよい。

【産業上の利用可能性】

【0109】

本発明は、例えば、光検出磁気共鳴を利用した測定装置に適用可能である。

【符号の説明】

【0110】

１ 磁気共鳴部材
２ 高周波磁場発生器
３ 磁石
１２ 発光装置
１３ 受光装置（蛍光受光装置の一例）
２３ 受光装置
２４ 受光装置
３１ 演算処理装置

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】