(57)【要約】【課題】遅延部の機械的移動による影響を抑制することができるテラヘルツ波の時間波形を取得するための装置及び方法を提供する。【解決手段】テラヘルツ波の発生部（１０１）と、テラヘルツ波の検出部（１０２）と、発生部と検出部とに照射される励起光の光路長差を変えるための第１の遅延部（１０３）と、励起光を発生部と検出部とに到達させるまでの時間差を変えるための第２の遅延部（１０４）を有する装置で、第２の遅延部における第１の光路長差で、第１の遅延部を用いて第１の時間波形を取得する。第２の遅延部を用いて、第１の光路長差を第２の光路長差に変える。第２の光路長差で、第１の遅延部を用いて第２の時間波形を取得する。取得された第１及び第２の時間波形を、第１及び第２の光路長差それぞれから予め定められた光路長差に合わせる。予め定められた光路長差に合わせられた第１及び第２の時間波形を加算平均する。【選択図】図１
2010281700.tif 000002
【特許請求の範囲】
【請求項１】
テラヘルツ波の時間波形を取得するための方法であって、
テラヘルツ波を発生させるための発生部と、該テラヘルツ波を検出するための検出部と、該発生部と該検出部とにそれぞれ照射される励起光の光路長差を変えるための第１の遅延部と、励起光を該発生部と該検出部とにそれぞれ到達させるまでの時間差を変えるための第２の遅延部と、を有する装置において、該第２の遅延部における第１の光路長差で、前記第１の遅延部を用いて第１の時間波形を取得するステップと、
前記第２の遅延部を用いて、前記第１の光路長差を該第１の光路長差とは異なる第２の光路長差に変えるステップと、
前記第２の光路長差で、前記第１の遅延部を用いて第２の時間波形を取得するステップと、
前記取得された第１及び第２の時間波形を、前記第１及び第２の光路長差それぞれから予め定められた光路長差に合わせるステップと、
前記予め定められた光路長差に合わせられた第１及び第２の時間波形を加算平均するステップと、を含むことを特徴とする方法。
【請求項２】
前記発生部から発生されたテラヘルツ波を検体に照射し、
前記検体を透過あるいは該検体を反射したテラヘルツ波を前記検出部で検出し、
前記加算平均により、前記第１の遅延部の移動に由来する周波数成分の抑制された時間波形を用いて、前記検体の情報を取得することを特徴とする請求項１に記載の方法。
【請求項３】
テラヘルツ波の時間波形を取得するための装置であって、
テラヘルツ波を発生させるための発生部と、
前記発生部から発生されたテラヘルツ波を検出するための検出部と、
前記発生部と前記検出部とにそれぞれ照射される励起光の光路長差を変えるための第１の遅延部と、
励起光を前記発生部と前記検出部とにそれぞれ到達させるまでの時間差を変えるための第２の遅延部と、
前記第２の遅延部における第１及び第２の光路長差それぞれで、前記第１の遅延部を用いて第１及び第２の時間波形を取得するための処理部と、を有し、
前記処理部は、前記取得された第１及び第２の時間波形を、前記第１及び第２の光路長差それぞれから予め定められた光路長差に合わせ、前記予め定められた光路長差に合わせられた第１及び第２の時間波形を加算平均することを特徴とする装置。
【請求項４】
前記発生部から発生されたテラヘルツ波を検体に照射し、
前記検体を透過あるいは該検体を反射したテラヘルツ波を前記検出部で検出し、
前記加算平均により、前記第１の遅延部の移動に由来する周波数成分の抑制された時間波形を用いて、前記検体の情報を取得することを特徴とする請求項３に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、テラヘルツ波の時間波形を取得するための装置及び方法に関する。特には、時間領域でのテラヘルツ波を測定するテラヘルツ波測定装置（THz-TDS装置、THz-Time Domain Spectroscopy装置）及び方法に関する。
【背景技術】
【０００２】
テラヘルツ波は、０.０３THz以上３０THz以下の範囲のうちの任意の周波数帯域の成分を有する電磁波である。この様な周波数帯域には、生体分子をはじめとして、様々な物質の構造や状態に由来する特徴的な吸収が多く存在する。この様な特徴を活かして、非破壊にて物質の分析や同定などを行う検査技術が開発されている。また、Ｘ線に替わる安全なイメージング技術や高速な通信技術などへの応用が期待されている。
【０００３】
テラヘルツ波の時間波形は、多くの場合、サブピコ秒のパルス形状である。この様なパルスを実時間で取得することは困難であるため、THz-TDS装置は、フェムト秒オーダのパルス幅を有する励起光によってサンプリング計測を行っている。このサンプリングは、テラヘルツ波を発生する発生部とこの波を検出する検出部に到達する励起光の時間差を調整することで実現される。例えば、この時間差は、折り返し光学系を有するステージ（遅延光学部とも呼ぶ）を励起光の伝播経路に挿入し、折り返し量を調整することで取得する。この発生部又は検出部としては、半導体膜に微小間隙を有するアンテナパターンを設けた光伝導素子を用いる場合がある。よって、諸所で、発生部又は検出部がこうした光伝導素子であるとして説明を進める。
【０００４】
近年、前記THz-TDS装置の測定感度が向上するに従い、遅延光学部に用いたステージの振動の影響が指摘されている（特許文献１参照）。遅延光学部に用いるステージが振動すると、励起光の光軸が揺れる。この結果、前記光伝導素子の微小間隙に照射される単位面積あたりの光量が変化し、信号のＳＮ比の悪化を招くことが予想される。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８-０１４８１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記の如き遅延光学部のステージとしては、数１０mmの移動を行う直動ステージを用いることがある。そして、直動ステージとして、ステッピングモータと精密ボールネジを組み合わせ、テーブルを移動するものが使用される場合がある。この様なステージでは、テーブルを機械的に移動するため、例えば、ボールネジのネジピッチ等の形状に応じて、テーブルの面がステージの設置基準面に対して変動することが起こり得る。また、雰囲気などの設置環境により、遅延光学部を支持する筐体や部材が歪むことや、場所固有の外部振動を拾うことにより、テーブルの面がステージの設置基準面に対して変動することも起こり得る。また、遅延光学部を支持する筐体や部材が光学的に歪む場合も、テーブルの面がステージの設置基準面に対して変動することも起こり得る。テーブルの移動によって、この変動が周期的な振動として現れる場合、励起光の光軸も周期的に揺れ、テラヘルツ波の時間波形に周期的な信号成分が重畳する。この様な時間波形をフーリエ変換すると、図２（ｄ）に示す様に、検出部が検出したテラヘルツ波のスペクトル（９１９）と共に、遅延光学部の振動による偽スペクトル（９２０）が現れる。例えば、テラヘルツ波の時間波形に数１００Hzの振動成分が重畳すると、４THzから６THz付近に偽スペクトルが現れ、測定装置の測定帯域を制限し、分析性能を低下させる要因となる。
【課題を解決するための手段】
【０００７】
前記課題に鑑み、本発明のテラヘルツ波の時間波形を取得するための方法は、
テラヘルツ波を発生させるための発生部と、該テラヘルツ波を検出するための検出部と、該発生部と該検出部とにそれぞれ照射される励起光の光路長差を変えるための第１の遅延部と、励起光を該発生部と該検出部とにそれぞれ到達させるまでの時間差を変えるための第２の遅延部と、を有する装置において、該第２の遅延部における第１の光路長差で、前記第１の遅延部を用いて第１の時間波形を取得するステップと、
前記第２の遅延部を用いて、前記第１の光路長差を該第１の光路長差とは異なる第２の光路長差に変えるステップと、
前記第２の光路長差で、前記第１の遅延部を用いて第２の時間波形を取得するステップと、
前記取得された第１及び第２の時間波形を、前記第１及び第２の光路長差それぞれから予め定められた光路長差に合わせるステップと、
前記予め定められた光路長差に合わせられた第１及び第２の時間波形を加算平均するステップと、を含むことを特徴とする。
【０００８】
また、前記課題に鑑み、本発明のテラヘルツ波の時間波形を取得するための装置は、
テラヘルツ波を発生させるための発生部と、
前記発生部から発生されたテラヘルツ波を検出するための検出部と、
前記発生部と前記検出部とにそれぞれ照射される励起光の光路長差を変えるための第１の遅延部と、
励起光を前記発生部と前記検出部とにそれぞれ到達させるまでの時間差を変えるための第２の遅延部と、
前記第２の遅延部における第１及び第２の光路長差それぞれで、前記第１の遅延部を用いて第１及び第２の時間波形を取得するための処理部と、を有し、
前記処理部は、前記取得された第１及び第２の時間波形を、前記第１及び第２の光路長差それぞれから予め定められた光路長差に合わせ、前記予め定められた光路長差に合わせられた第１及び第２の時間波形を加算平均することを特徴とする。
【発明の効果】
【０００９】
本発明によれば、測定した複数のテラヘルツ波時間波形を、第２の遅延部における複数の光路長差それぞれから予め定められた光路長差に合わせ、予め定められた光路長差に合わせられた複数の時間波形を加算平均する。こうして、テラヘルツ波の時間波形に重畳する装置由来の振動について、この振動による周波数成分の位相を変化させて測定を行っている。そして、これらの時間波形を加算平均することで、こうした振動が時間波形に与える影響を抑制している。この結果、例えば、第１の遅延部の周期的な振動由来の偽スペクトルを抑制できるので、測定帯域が広がり、分析性能を向上することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】本発明によるテラヘルツ波測定装置及び方法の一実施形態の概略構成、及び第１の遅延部と第２の遅延部の動作を説明する図。
【図２】テラヘルツ波の時間波形を処理する本発明における工程例、及び従来技術における遅延光学部の振動による偽スペクトルを説明する図。
【図３】本発明の実施例１におけるテラヘルツ波測定装置及び方法の動作フローを説明する図。
【図４】本発明の実施例２におけるテラヘルツ波測定装置及び方法の動作フローを説明する図。
【図５】本発明の実施例３におけるテラヘルツ波測定装置及び方法の動作フローを説明する図。
【図６】第１の遅延部と第２の遅延部の配置の他の例を示す図。
【発明を実施するための形態】
【００１１】
以下、本発明の実施形態を説明する。本発明のテラヘルツ波の時間波形を取得するための装置及び方法において重要なことは、次の点である。測定した複数のテラヘルツ波時間波形を、第２の遅延部における複数の光路長差それぞれから予め定められた光路長差に合わせ、予め定められた光路長差に合わせられた複数の時間波形を加算平均することである。
【００１２】
前記考え方に基づき、テラヘルツ波の時間波形を取得するための本発明の技術の基本的な実施形態は、上述した様な構成を有する。より具体的には、本発明の方法は、第１の遅延部を備える装置において、テラヘルツ波の時間波形を取得するための方法であって、次の第１から第５のステップを含む。ここで、第１の遅延部は、テラヘルツ波を発生する発生部と検出する検出部とに夫々到達する励起光の時間差を、励起光が伝播する光路長差により調整する。第１のステップでは、前記光路長差を第１の光路長差（第１の値）に設定した第１の状態で、第１の遅延部の移動を用いて第１の時間波形を測定する。第２のステップでは、前記第１の状態から、前記光路長差を第１の値とは異なる第２の光路長差（第２の値）に設定した第２の状態に移す。第３のステップでは、前記第２の状態で、第１の遅延部の移動を用いて第２の時間波形を測定する。第４のステップでは、前記第１及び第２のテラヘルツ波時間波形を構成するデータ列について、前記光路長差が基準の値である基準の状態における第１の遅延部の予め定められた位置からのデータ列に変換する。すなわち、前記光路長差を変えるための第２の遅延部における複数の光路長差それぞれから予め定められた光路長差に合わせる。第５のステップでは、前記第４のステップを経た時間波形について加算平均を行う。また、テラヘルツ波の時間波形を取得するための本発明の装置は、第１の遅延部と第２の遅延部と処理部を有する。第１の遅延部は、所定の移動を行うことで、テラヘルツ波を発生する発生部とテラヘルツ波を検出する検出部とに夫々到達する励起光の時間差を、励起光が伝播する光路長差の変化により変化させる。第２の遅延部は、第１の遅延部が前記所定の移動を行うところの状態における前記光路長差を調整する。処理部は、第２の遅延部により前記光路長差が夫々異なる値に調整された複数の状態において第１の遅延部が前記所定の移動を行うことで測定される複数の時間波形を加算平均する。また、処理部は、前記複数の時間波形のデータ列について、前記光路長差が基準の値である基準の状態における第１の遅延部の予め定められた位置からのデータ列に変換し、変換したデータ列について加算平均を行う。
【００１３】
前記基本的な実施形態に基づいて、次に述べる様なより具体的な実施形態が可能である。
前記第１の遅延部が所定の移動を行う状態において、前記光路長差は第２の遅延部によって設定することができる。また、テラヘルツ波取得技術を用いて、発生部からのテラヘルツ波を検体に照射し検体を透過或いは検体で反射したテラヘルツ波を検出部で検出して、前記第１の遅延部の移動に由来する実時間における予め定めた周波数成分の信号を抑制した状態で検体の情報を取得することができる。
【００１４】
前記第１の遅延部は、時間波形測定毎に、時間波形測定のためにホームポジションである所定の出発位置から所定の移動を行うものである。これは、機械的な移動をする可動部を含むものである。後述する例の如く、可動部が折り返しミラーの様な光学系である場合は第１の遅延光学部と称する。前記第２の遅延部は、時間波形測定毎に、第１の遅延部の前記所定の移動の前記所定の出発位置における励起光の光路長差を調整するものである。故に、前記所定の移動の各位置において、第１の遅延部による光路長差はこの調整前のものからこの調整量だけずれる。第２の遅延部は、機械的な移動をする可動部を含むものであっても、励起光が伝搬する光路中の時定数（屈折率）を電気光学結晶の電気的制御などで変化させるものなどであってもよい。要するに、第２の遅延部は、第１の遅延部の出発位置における励起光の光路長差を調整できるものであればら、どの様なものでもよい。第２の遅延部についても、後述する例の如く、可動部を備え、この可動部が折り返しミラーの様な光学系である場合は第２の遅延光学部と称する。光路長差を第１の値に設定した第１の状態や第２の値に設定した第２の状態などの「状態」について、可動部を含む第２の遅延部の当該可動部の位置の移動で光路長差の調整を行う場合は、後述する例の如く、第１の位置や第２の位置などと「位置」と称する。第２の遅延部の調整量は、後述する例における測定位置出力部などにより定めることができる。この調整量は、第１の遅延部の移動に由来する実時間における予め定めた周波数成分の信号を効果的に抑制すると言う観点から、前記予め定められた位置における光路長差の値を含むＶ・Ｔの光路長差の調整範囲にあるのが好ましい。ここで、Ｖは、第１の遅延部の移動の動作速度、Ｔは、第１の遅延部の移動に由来する実時間における予め定めた周波数成分の周期である。
【００１５】
次に、本発明の思想を実施し得る形態について、図面を参照して説明する。図１に示すテラヘルツ波測定装置の実施形態は、THz-TDS装置の基本構成を踏襲しており、時間領域でテラヘルツ波の時間波形を測定するものである。テラヘルツ波の時間波形を測定する部分は、発生部（１０１）、検出部（１０２）、第１の遅延光学部（１０３）、処理部（１０６）、バイアス印加部（１０７）、電流検出部（１０８）、レーザ源（１０９）を含む。本実施形態では、第１の遅延光学部（１０３）に由来する周期的な振動成分を抑制するために、この基本構成に加え、第２の遅延光学部（１０４）と測定位置出力部（１０５）を備える。
【００１６】
発生部（１０１）は、テラヘルツ波を発生する部分である。発生部（１０１）では、瞬時電流を利用する手法やキャリヤのバンド間遷移を利用する手法を用いることができる。前者の手法としては、半導体や有機結晶の表面にレーザ光を照射してテラヘルツ波を発生する手法や、半導体薄膜上に金属電極でアンテナパターンを形成した光伝導素子に電界を印加してレーザ光を照射する方法などが適用できる。また、PINダイオードが適用できる。利得構造を利用する手法としては、半導体量子井戸構造を用いる手法が適用できる。本実施形態では、光伝導素子を用いる。検出部（１０２）は、テラヘルツ波の電界強度を検出する部分である。検出部（１０２）では、光伝導による電界強度に対応した電流を検出する手法、電気光学効果を用いて電場を検出する手法、磁気光学効果を用いて磁場を検出する手法などを用いることができる。光伝導による電流を検出する手法としては、光伝導素子が適用できる。電気光学効果を用いて電場を検出する手法としては、直交偏光子と電気光学結晶を使った手法が適用できる。磁気光学効果を用いて磁場を検出する手法としては、直交偏光子と磁気光学結晶を使った手法が適用できる。検出部（１０２）に入射するテラヘルツ波は、検出部（１０２）に焦点を結ぶことで、単位面積あたりの強度を増やし、検出感度を高められる。本実施形態では、光伝導素子を用いる。
【００１７】
レーザ源（１０９）は、超短パルスレーザを出力する部分である。発生部（１０１）と検出部（１０２）は、この超短パルスレーザによってキャリヤが励起されることで動作する。図１（ａ）に示す様に、超短パルスレーザは、L１とL２の二つの光路に分岐される。光路L１を通る超短パルスレーザは、発生部（１０１）に入力する。光路L２を通る超短パルスレーザは、第１の遅延光学部（１０３）を介して検出部（１０２）に入力する。
【００１８】
第１の遅延光学部（１０３）は、テラヘルツ波の時間波形のうちの、検出部（１０２）でサンプリングする位置を調整する部分である。具体的には、例えば、発生部（１０１）に入力される励起光L１に対し、検出部（１０２）に入力される励起光L２のタイミングを遅延させる。この励起光L１とL２の時間差の調整法としては、励起光が伝播する光路長を直接調整する手法と実効的な光路長を調整する手法がある。前者の手法としては、折り返し光学系と可動部を用いる方法がある。後者の手法としては、励起光が伝搬する光路中の時定数を変化させる方法がある。例えば、厚みの変化する光学結晶を機械的移動機構で光路中に出し入れする方法などがある。何れの手法を用いるものであれ、第１の遅延光学部（１０３）は、機械的な可動部を備える。図１（ａ）は、１段の折り返し光学系と可動部としての直動ステージとを用いる例を示している。これに対して、後述の第２の遅延光学部（１０４）は、励起光L１とL２の時間差を調整できるものであれば、機械的な可動部を備えるものに限らず、どの様なものでもよい。例えば、電気光学結晶に印加する電圧を変化させて、電気光学結晶を挿入した光路の時定数を変化させる様な手段を用いることもできる。
【００１９】
バイアス印加部（１０７）は、発生部（１０１）を駆動するためのバイアスを供給する部分である。発生部（１０１）として光伝導素子を用いる場合、アンテナパターンを含む金属電極にバイアスを印加する。特に、電流検出部（１０８）がロックイン検出系を含む場合、バイアス印加部（１０７）が供給するバイアスをロックイン検出系の参照信号と同等の周波数で変調する方法を用いることができる。電流検出部（１０８）は、電流信号を測定可能な電圧信号に変換する部分である。検出部（１０２）として光伝導素子を用いる場合、電流検出部（１０８）は、検出部（１０２）から出力される電流信号を電圧信号に変換する。前述した様に、検出部（１０２）から出力される信号が微小な場合、ロックイン検出系を含むことがある。ただし、本発明においてこの電流検出部（１０８）は必須でなく、検出部（１０２）の信号に応じて、この信号を、処理部（１０６）が処理し易い信号に調整し得る手段に置き換えることができる。
【００２０】
処理部（１０６）は、テラヘルツ波の時間波形を構築する部分である。第１の遅延光学部（１０３）の調整量（本実施形態では直動ステージの移動量）と電流検出部（１０８）の出力を参照し、時間波形を構築する。そして、処理部（１０６）は、複数の時間波形を加算平均して信号のSN比の向上を図る。ここで、周波数スペクトルを求める場合、この時間波形をフーリエ変換する。THz-TDS装置を分析装置として用いる場合は、図１（ａ）に示す様に配置されるサンプルにテラヘルツ波を照射したときの時間波形の変化を求める。また、処理部（１０６）は、照射されるテラヘルツ波とサンプルとの相対的な位置をモニタすることで、サンプルを可視化することができる。以上の構成は、THz-TDS装置として知られる一般的な構成である。
【００２１】
図１（ａ）に示す如く、本発明の実施形態では、第１の遅延光学部（１０３）由来の振動成分を抑制するための構成として、第２の遅延光学部（１０４）と測定位置出力部（１０５）を備える。本明細書では、この振動成分をテラヘルツ波時間波形における所定の周波数成分とも表現する。
【００２２】
第２の遅延光学部（１０４）は、テラヘルツ波の時間波形を測定する位置を定める部分である。第１の遅延光学部（１０３）と同じく、第２の遅延光学部（１０４）は、励起光L１とL２が発生部（１０１）と検出部（１０２）に到達する時間差を調整する。詳細には、第１の遅延光学部（１０３）がテラヘルツ波時間波形のサンプリングの点を調整するのに対し、第２の遅延光学部（１０４）は、時間波形の測定毎に、第１の遅延光学部（１０３）で決められる各サンプリングの位置を僅かに変化させる。詳細な動作については後述する。測定位置出力部（１０５）は、第２の遅延光学部（１０４）の移動量を調整してこの移動量を処理部（１０６）に出力する部分である。本装置は、テラヘルツ波を繰り返し測定し、これらのテラヘルツ波の時間波形を加算平均するものであるが、測定位置出力部（１０５）は、テラヘルツ波の測定毎に、第２の遅延光学部（１０４）の移動量を調整する。
【００２３】
次に、本実施形態の装置の動作を説明する。本装置は、テラヘルツ波の時間波形の測定を複数回行い、これらの時間波形を加算平均する装置である。図１（ａ）において、励起光L１の照射によって発生部（１０１）はテラヘルツ波（T１）を発生する。テラヘルツ波の伝播経路中にサンプルがある場合、テラヘルツ波（T１）は、サンプルによる吸収や散乱の影響を含んだテラヘルツ波（T２）になり、検出部（１０２）に入射する。検出部（１０２）に入射されるテラヘルツ波（T２）は、第１の遅延光学部（１０３）を経た励起光（L２）によってサンプリング計測される。サンプリングされる点は、第１の遅延光学部（１０３）で励起光L１とL２が伝播する光路長差を変化させることによって調整される。サンプリングされる点の信号は、電流検出部（１０８）を経て処理部（１０６）に送られ、第１の遅延光学部（１０３）の調整量と共にプロットされ、時間波形を構築する。ここまでの動作は、一般的なTHz-TDS装置の動作である。本実施形態の装置は、この時間波形を複数回測定する間に、第２の遅延光学部（１０４）を調整して、所定の周波数成分の信号（第１の遅延光学部（１０３）由来の励起光の振動）を抑制している。この部分の動作を以下に説明する。
【００２４】
図１（ｂ）に示す様に、第２の遅延光学部（１０４）は第１の遅延光学部（１０３）を構成するステージ上に設置され、第２の遅延光学部（１０４）上に、第１の遅延光学部（１０３）を構成する折り返しミラーが設置されている。図１（ａ）の処理部（１０６）は、第２の遅延光学部（１０４）と第１の遅延光学部（１０３）との位置を決めて図１（ｂ）の基準の位置Ａ（２１０）を定める。位置Ａ（２１０）は、例えば、時間波形のピークが生じる位置であり、最終的に、処理部（１０６）はこの位置Ａ（２１０）におけるテラヘルツ波の時間波形を提供する。測定位置出力部（１０５）は、測定毎に第２の遅延光学部（１０４）を調整することで、位置Ａ（２１０）を含む範囲において、測定する位置を変化させる。例えば、図１（ｂ）によると、測定回数がｎ-１の時、第２の遅延光学部（１０４）は第１の遅延光学部（１０３）を構成する折り返しミラーを第１の位置（２１１）に調整する。そして、この状態で、第１の遅延光学部（１０３）を用いて第１のテラヘルツ波の時間波形（２１４）を取得する。また、測定回数がｎの時、第２の遅延光学部（１０４）は折り返しミラーを第１の位置（２１１）とは異なる第２の位置（２１２）に調整する。そして、この状態で、第１の遅延光学部（１０３）を用いて第２のテラヘルツ波の時間波形（２１５）を取得する。
【００２５】
この時、第１のテラヘルツ波時間波形（２１４）と第２のテラヘルツ波時間波形（２１５）を比較すると、テラヘルツ波パルスのピークが生じる位置は異なるが、第１の遅延光学系（１０３）に由来する振動成分の信号の位相は殆ど変化がない。これは、振動成分の信号が第１の遅延光学部（１０３）に由来するものであるため、第１の遅延光学部（１０３）のステージの駆動開始位置や駆動速度（即ち、測定毎の第１の遅延光学部（１０３）の振動態様）を変化させない限り固定されるためである。図１（ｂ）では、説明の都合上、この所定の周波数成分の信号をテラヘルツ波パルスと独立して表現しているが、実際には、テラヘルツ波の時間波形全体に対して重畳している。
【００２６】
テラヘルツ波の時間波形の測定毎に、例えば、第１の遅延光学部（１０３）を１０mm/s程度の速度Ｖで動かす場合、ボールネジの形状等を勘案すると、第１の遅延光学部（１０３）は数kHzの周波数で振動することが多い。この時、テラヘルツ波の時間波形について、第１の遅延光学部（１０３）が数μmから数１０μm移動する毎にサンプリングを行うとすると、サンプリング周波数が振動成分の周波数よりも小さくなる。具体的には、エイリアシング（Aliasing）の影響で振動成分は数１００Hzの信号に変換される。この数１００Hzの振動成分の信号は、テラヘルツ波の時間波形に換算すると５THz程度の周波数成分となる。典型的には、図２（ｄ）に示す様に、遅延光学部の振動による偽スペクトル（９２０）となる。本明細書では、テラヘルツ波の時間波形に換算される前の信号成分（時間軸が第１の遅延光学部（１０３）の移動速度で定義される信号成分）を実時間における信号成分とする。この実時間における周波数成分の周期をＴとする時、第２の遅延光学部（１０４）による光路長差の調整範囲（２１３）はＶ・Ｔとするのが好ましい。図１（ｂ）では、１段の折り返し光学系を用いているため、第２の遅延光学部（１０４）は、（１/２）・Ｖ・Ｔの範囲で動かす。
【００２７】
上述したテラヘルツ波の時間波形を処理する工程を説明する。測定毎に第２の遅延光学部（１０４）で測定する位置を変化させると、処理部（１０６）では、図２（a）の様に、パルスのピーク位置が異なる時間波形が構築される。処理部（１０６）では、これらの時間波形を加算平均するために、パルスのピーク位置を位置Ａ（２１０）に調整する。具体的には、テラヘルツ波時間波形を構成するデータ列について、パルスのピーク位置が位置Ａ（２１０）となる様にデータ列を構成する時間軸を移動する。図２（ｂ）は、このデータ列の変換ステップを経た後の各時間波形を示す。図２（ｂ）によると、この変換ステップを経た時間波形は、パルスのピーク位置が位置Ａとなるが、第１の遅延光学部（１０３）に由来する振動成分の位相は異なっている。この様に、振動成分の位相は異なっているため、これらの波形を加算平均すると、図２（ｃ）の様に、振動成分が重畳した従来の時間波形（３１７）とは異なり、所定の周波数成分を抑制した時間波形（３１８）を取得できる。
【００２８】
第２の遅延光学部（１０４）によって調整される光路長差の調整範囲（２１３）をＶ・Ｔにすると、振動成分の位相変化を１波長以内にすることができる。よって、振動の繰り返しの影響（例えば、１つ目と２つ目の振動波形について相対的に同じ位置となること）を防止することができる。そのため、加算平均する際の振動成分の抑制をより効率的に行うことができる。すなわち、単位測定時間あたりの周期的な振動成分の抑制効果を向上することが可能となる。
【００２９】
また、位置Ａ（２１０）に対し、第１の位置（２１１）と第２の位置（２１２）を擬似ランダムに設定してもよい。この時、各テラヘルツ波の時間波形を加算平均する際、変換ステップを経た振動成分の位相の相関性が小さい（例えば、ほぼ０）ため、所定の周波数成分の抑制効果が高まる。本実施形態では、位置Ａ（２１０）は、光路長差の調整範囲（２１３）の中心付近に設定しているが、これに限るものではなく、例えば、光路長差の調整範囲（２１３）の境界にあってもよい。
【００３０】
また、本実施形態では、第２の遅延光学部（１０４）は第１の遅延光学部（１０３）の上に設置されているが、これに限らない。例えば、図６（ａ）の様に、励起光L２の伝播光路に第１の遅延光学部（１０３）と第２の遅延光学部（１０４）が個々別個に配置されていてもよい。この場合、励起光L１の伝播光路に個々に配置されてもよい。また、図６（ｂ）の様に、励起光L１とL２の伝播光路に夫々分離して配置される構成も可能である。この様な構成を取ることで、装置形態や設置環境に応じて柔軟に対応することができる。
【００３１】
以上の様な構成によって、本実施形態のテラヘルツ波測定装置は、第１の遅延光学部（１０３）の形状に由来する周期的な振動成分を抑制することができる。そのため、周期的な振動由来の偽スペクトルを抑制できるので、測定装置の測定帯域が広がり、分析性能を向上することが可能となる。
【００３２】
以下、図を用いて、動作フローなどに特徴を持つ実施例を説明する。
（実施例１）
図３は、実施例１におけるテラヘルツ波測定装置の動作フローを示す。本実施例の基本的構成は前記実施形態と共通するので、その部分の説明は省略する。図３において、テラヘルツ波測定装置が測定を開始すると、処理部（１０６）は、位置Ａ（２１０）を決定する（S４０１）。第１の遅延光学部（１０３）は位置Ａ（２１０）に移動され、第１の遅延光学部（１０３）は、この時の位置データを記憶部に格納する。より具体的には、第２の遅延光学部（１０４）の位置は固定され、第１の遅延光学部の位置が位置Ａに移動される。この後、第２の遅延光学部（１０４）は、時間波形の測定毎に、位置Ａを基準として移動される。記憶部には、位置Ａの時の各遅延光学部の位置が記憶される。記憶部は処理部（１０６）の内部に備えられていてもよいし、外部に備えられていてもよい。
【００３３】
次に、測定位置出力部（１０５）が、第２の遅延光学部（１０４）を第１の位置（２１１）に移動する（S４０２）。測定位置出力部（１０５）は、この時の測定位置出力部（１０５）の位置データを記憶部に格納する。第２の遅延光学部（１０４）の可動範囲は光路長差の調整範囲（２１３）内である。
【００３４】
テラヘルツ波測定装置は、この状態で、第１の遅延光学部（１０３）によってサンプリング点を移動し、第１のテラヘルツ波時間波形（２１４）を取得する（S４０３）。これは、従来のテラヘルツ波測定装置と同じ動作である。取得された波形データは、第１の遅延光学部（１０３）の位置と電流検出部（１０８）の強度信号の配列である。この位置データは、テラヘルツ波領域のサンプリングの時間間隔に相当する。そのため、位置のデータ列は、テラヘルツ波領域の時間のデータ列に変換してもよい。この第１のテラヘルツ波時間波形（２１４）の波形データは記憶部に格納される。
【００３５】
処理部（１０６）は記憶部を参照し、波形データを構成する位置データ列を、位置Ａ（２１０）における位置データ列に変換する（S４０８）。具体的には、位置Ａ（２１０）と第１の位置（２１１）のデータを記憶部より参照し、位置Ａ（２１０）と第１の位置（２１１）の差分を求める。そして、波形データの位置データ列について、この差分の値で補正する。位置データ列が時間データ列の場合、位置データ列は、差分の値を時間に換算して補正される。変換された波形データは、再び記憶部に格納される。
【００３６】
次に、測定位置出力部（１０５）は、第２の遅延光学部（１０４）を第２の位置（２１２）に移動する（S４０４）。測定位置出力部（１０５）は、この時の第２の遅延光学部（１０４）の位置データを記憶部に格納する。ここでも、前記第１の位置（２１１）におけるステップと同様なステップが実行される（S４０５、S４０８）。こうして、第２のテラヘルツ波時間波形（２１５）について、変換された波形データが記憶部に格納される。
【００３７】
続いて、テラヘルツ波測定装置は時間波形の取得回数をモニタし、測定回数が所定回数に達したかを判断する（S４０６）。測定回数が所定回数に満たない場合、第２の位置（２１２）を第１の位置（２１１）に再定義し、（S４０４）から（S４０５）の工程を実行する（S４０７）。測定回数が所定回数に達した場合、処理部（１０６）は記憶部に格納された複数の波形データを参照し、測定回数によって加算平均を行う（S４０９）。この工程は、（S４０６）の工程の判断を待たず、波形データが記憶部に格納される度に逐次行う形態にしてもよい。加算平均した波形データは、作業者の望むデータ形式に処理され、提示される（S４１０）。テラヘルツ波測定装置は、これら一連の工程が終了したら、測定を終了する。
【００３８】
本実施例のこの様な工程によっても、前記実施形態の説明のところで述べたものと同様な効果が奏される。
【００３９】
（実施例２）
実施例２は、実施例１に示した動作フローの変形例に係る。具体的には、実施例１とは、テラヘルツ波測定装置の測定終了を判断する基準が異なる。ここでも、これまでの説明と共通する部分の説明は省略する。図４は、本実施例におけるテラヘルツ波測定装置の動作フローを示す。テラヘルツ波測定装置が測定を開始すると、処理部（１０６）は、位置Ａ（２１０）を決定する（S４０１）。ここでの第１の遅延光学部（１０３）の位置と第２の遅延光学部（１０４）の位置との情報は記憶しておく。次に、測定位置出力部（１０５）は、第２の遅延光学部（１０４）を第１の位置（２１１）に移動する（S４０２）。この状態において、テラヘルツ波測定装置は、第１の遅延光学部（１０３）の所定の移動によってサンプリング点を移動し、第１のテラヘルツ波時間波形（２１４）を取得する（S４０３）。次に、処理部（１０６）は記憶部に記憶されている前記情報を参照し、波形データを構成する位置データ列を、位置Ａ（２１０）における位置データ列に変換する（S４０８）。続いて、測定位置出力部（１０５）は、第２の遅延光学部（１０４）を第２の位置（２１２）に移動する（S４０４）。この状態で、テラヘルツ波測定装置は、第１の遅延光学部（１０３）の所定の移動によってサンプリング点を移動し、第２のテラヘルツ波時間波形（２１５）を取得する（S４０５）。ここでも、処理部（１０６）は記憶部の情報を参照し、波形データを構成する位置データ列を、位置Ａ（２１０）における位置データ列に変換する（S４０８）。
【００４０】
本実施例では、第２のテラヘルツ波時間波形（２１５）を測定した段階で、逐次加算平均を行った（S４０９）後、テラヘルツ波測定装置は、加算平均したテラヘルツ波時間波形のSN比を求め、所望のSN比が得られたかを判断する（S５０６）。時間波形が所望のSN比に満たない場合、テラヘルツ波測定装置は、第２の位置（２１２）を第１の位置（２１１）に再定義し、（S４０４）から（S４０９）の工程を行う（S４０７）。時間波形が所望のSN比に達した場合、テラヘルツ波測定装置は、繰り返し測定を停止する。そして、加算平均した波形データは、作業者の望むデータ形式に処理され、提示される（S４１０）。テラヘルツ波測定装置は、これら一連の工程が終了したら、測定を終了する。
【００４１】
本実施例のこの様な工程によって、テラヘルツ波測定装置は、所望のSN比に達したら測定を終了することができる。そのため、装置は、常に同じSN比の時間波形を提供できるので、複数の時間波形の比較を行う際、精度が向上する。また、所望のSN比に達したら測定を終了するため、余分な測定を行うことがなく、作業効率が向上する。
【００４２】
（実施例３）
実施例３は、実施例２に示した動作フローの変形例に係る。具体的には、実施例２とは、テラヘルツ波測定装置の測定終了を判断する基準が異なる。ここでも、これまでの説明と共通する部分の説明は省略する。図５は、本実施例におけるテラヘルツ波測定装置の動作フローを示す。テラヘルツ波測定装置が測定を開始すると、処理部（１０６）は、位置Ａ（２１０）を決定する（S４０１）。測定位置出力部（１０５）は、第２の遅延光学部（１０４）を第１の位置（２１１）に移動する（S４０２）。テラヘルツ波測定装置は、第１の遅延光学部（１０３）によってサンプリング点を移動し、第１のテラヘルツ波時間波形（２１４）を取得する（S４０３）。処理部（１０６）は記憶部の情報を参照し、波形データを構成する位置データ列を、位置Ａ（２１０）における位置データ列に変換する（S４０８）。測定位置出力部（１０５）は、第２の遅延光学部（１０４）を第２の位置（２１２）に移動する（S４０４）。装置は、第１の遅延光学部（１０３）によってサンプリング点を移動し、第２のテラヘルツ波時間波形（２１５）を取得する（S４０５）。処理部（１０６）は記憶部の情報を参照し、波形データを構成する位置データ列を、位置Ａ（２１０）における位置データ列に変換する（S４０８）。
【００４３】
本実施例では、第２のテラヘルツ波時間波形（２１５）を測定した段階で、逐次加算平均を行う（S４０９）。装置は、加算平均したテラヘルツ波時間波形のSN比を求め、所望のSN比が得られたかを判断する（S５０６）。更に、装置は、時間波形の取得回数をモニタし、測定回数が所定回数に達したかを判断する（S４０６）。時間波形のSN比が所望のSN比に満たなく、測定回数が所定回数に満たない場合、装置は、第２の位置（２１２）を第１の位置（２１１）に再定義し、（S４０４）から（S４０９）の工程を行う（S４０７）。
【００４４】
時間波形が所望のSN比に達するか、測定回数が所定回数に達した場合、装置は繰り返し測定を停止する。そして、加算平均した波形データは、作業者の望むデータ形式に処理され、提示される（S４１０）。テラヘルツ波測定装置は、これら一連の工程が終了したら、測定を終了する。
【００４５】
本実施例のこの様な工程によって、テラヘルツ波測定装置は、何らかの理由で所望のSN比に達することが困難な場合であっても、一定時間で測定を終了することができる。本実施例によれば、終了条件を複数用意することで、汎用性が高まる。より具体的には、例えば、過剰な測定時間がかかることを防止し、測定者へのストレスを軽減することができる。また、不測の事態にも対応が容易になる。
【符号の説明】
【００４６】
１０１  発生部：１０２  検出部：１０３、７０３、８０３  第１の遅延光学部（第１の遅延部）：１０４、７０４、８０４  第２の遅延光学部（第２の遅延部）：１０５  測定位置出力部：１０６  処理部：２１０  位置Ａ（第１の遅延部の予め定められた位置、予め定められた光路長差の状態）：２１１  第１の位置（第１の状態）：２１２  第２の位置（第２の状態）：２１３  光路長差の調整範囲：２１４  第１の時間波形：２１５  第２の時間波形：３１８  本発明の時間波形：L１、L２  励起光（光路）
【図１】
2010281700.tif 000003
【図２】
2010281700.tif 000004
【図３】
2010281700.tif 000005
【図４】
2010281700.tif 000006
【図５】
2010281700.tif 000007
【図６】
2010281700.tif 000008