JP2021047130A - Spectrometry device and method - Google Patents
Spectrometry device and method Download PDFInfo
- Publication number
- JP2021047130A JP2021047130A JP2019170884A JP2019170884A JP2021047130A JP 2021047130 A JP2021047130 A JP 2021047130A JP 2019170884 A JP2019170884 A JP 2019170884A JP 2019170884 A JP2019170884 A JP 2019170884A JP 2021047130 A JP2021047130 A JP 2021047130A
- Authority
- JP
- Japan
- Prior art keywords
- light
- observation data
- spectrum
- time
- probe
- 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.)
- Pending
Links
Images
Landscapes
- Spectrometry And Color Measurement (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
Description
本発明は、分光測定装置及び方法に関する。 The present invention relates to a spectroscopic measuring device and a method.
測定対象の分子振動スペクトルを測定する分光測定装置として、フーリエ変換型コヒーレントアンチストークスラマン散乱(FT-CARS)分光方法を用いたものが知られている。FT-CARS分光方法では、ポンプ光及びストークス光と、ポンプ光及びストークス光に対して遅延させたプローブ光を照射する。このプローブ光の遅延時間を変化させながら、プローブ光の照射によって測定対象から放出されたアンチストークス光の時間領域インターフェログラムを取得し、この時間領域インターフェログラムに対してデータ解析手法としてのフーリエ変換を施し測定対象の分子振動スペクトルを得ている。 As a spectroscopic measuring device for measuring the molecular vibration spectrum to be measured, a spectroscopic measuring device using a Fourier transform coherent anti-Stoke Raman scattering (FT-CARS) spectroscopic method is known. In the FT-CARS spectroscopy method, pump light and Stokes light and probe light delayed with respect to the pump light and Stokes light are irradiated. While changing the delay time of this probe light, the time domain interferogram of the anti-Stokes light emitted from the measurement target by the irradiation of the probe light is acquired, and Fourier as a data analysis method for this time domain interferogram is obtained. The molecular vibration spectrum to be measured is obtained by conversion.
また、ポンプ光及びストークス光となる参照光に対するプローブ光の遅延時間を高速に変化させるFT-CARS分光方法を用いた分光測定装置が、特許文献1によって知られている。特許文献1の分光測定装置では、光源から放出されるパルス状のレーザ光をビームスプリッタで参照光とプローブ光とに分割して、参照光を光路長が固定された第1アームに伝搬させ、プローブ光を光路長が回転駆動されるスキャニングミラーによって増減される第2アームに伝搬させる構成である。 Further, Patent Document 1 discloses a spectroscopic measurement apparatus using an FT-CARS spectroscopic method in which the delay time of a probe light with respect to a reference light serving as a pump light and a Stokes light is changed at high speed. In the spectroscopic measurement device of Patent Document 1, a pulsed laser beam emitted from a light source is split into a reference light and a probe light by a beam splitter, and the reference light is propagated to a first arm having a fixed optical path length. The probe light is propagated to the second arm whose optical path length is increased or decreased by the scanning mirror driven by rotation.
ところで、上記のような分光測定装置においては、より高い測定速度が望まれている。しかしながら、光源からのパルス状のレーザ光の繰り返し周波数を一定にしたまま、測定速度を高くするためにスキャニングミラーの走査速度(回転速度)を高くすると、サンプリング間隔が増大する。このサンプリング間隔が、標本化定理で定まる間隔よりも大きくなると、エイリアシングによるアーチフェクトが発生し、時間領域インターフェログラムのサンプリング結果から適切な分子振動スペクトルが得られないという問題があった。 By the way, in the spectroscopic measuring apparatus as described above, a higher measurement speed is desired. However, if the scanning speed (rotation speed) of the scanning mirror is increased in order to increase the measurement speed while keeping the repetition frequency of the pulsed laser beam from the light source constant, the sampling interval increases. If this sampling interval is larger than the interval determined by the sampling theorem, an artifact due to aliasing occurs, and there is a problem that an appropriate molecular vibration spectrum cannot be obtained from the sampling result of the time domain interferogram.
本発明は、上記事情を鑑みてなされたものであり、スペクトルの測定速度の向上に資する分光測定装置及び方法を提供することを目的とする。 The present invention has been made in view of the above circumstances, and an object of the present invention is to provide a spectroscopic measuring device and a method that contribute to improving the measurement speed of a spectrum.
本発明の分光測定装置は、ポンプ光とこのポンプ光に対して遅延するプローブ光とからなるパルス対を生成し、ポンプ光の光路長に対してプローブ光の光路長を時間的に非線形に変化させるパルス対生成部と、前記パルス対生成部からのパルス対が繰り返し照射される測定対象からプローブ光が照射されるごとに得られる観測光の強度を時間領域インターフェログラムの観測データとして取得する観測データ取得部と、前記観測データ取得部によって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成部とを備えるものである。 The spectroscopic measuring apparatus of the present invention generates a pulse pair consisting of a pump light and a probe light delayed with respect to the pump light, and changes the optical path length of the probe light in a non-linear manner with respect to the optical path length of the pump light. The intensity of the observation light obtained each time the probe light is irradiated from the measurement target to which the pulse pair generation unit to be generated and the pulse pair from the pulse pair generation unit are repeatedly irradiated is acquired as the observation data of the time region interferogram. Spectral spectrum generation that matches the observation data acquisition unit and the observation data acquired by the observation data acquisition unit and obtains a spectrum that minimizes the L1 norm or L0 norm as a spectrum obtained by converting the time region interferogram into the frequency region. It is equipped with a part.
本発明の分光測定方法は、ポンプ光に対する遅延時間を時間的に非線形に変化させて、ポンプ光とプローブ光とからなるパルス対を繰り返し生成するパルス対生成ステップと、前記パルス対生成ステップで生成されるパルス対が繰り返し照射される測定対象からプローブ光が照射されるごとに得られる観測光の強度を時間領域インターフェログラムの観測データとして取得する観測データ取得ステップと、前記観測データ取得ステップで取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成ステップとを有するものである。 The spectral measurement method of the present invention is generated by a pulse pair generation step of repeatedly generating a pulse pair consisting of a pump light and a probe light by changing the delay time with respect to the pump light in a non-linear manner, and the pulse pair generation step. In the observation data acquisition step of acquiring the intensity of the observation light obtained each time the probe light is irradiated from the measurement target to which the pulse pair is repeatedly irradiated as the observation data of the time domain interferogram, and the observation data acquisition step. It has a spectral spectrum generation step of obtaining a spectrum that fits the acquired observation data and minimizes the L1 norm or the L0 norm as a spectrum obtained by converting a time domain interferogram into a frequency domain.
本発明の分光測定装置は、参照光に対する信号光の位相を時間的に非線形に変化させて参照光と信号光とを合波した干渉光を生成する干渉計と、前記干渉計からの干渉光を測定対象に照射することにより得られる観測光の強度を一定の時間間隔で、時間領域インターフェログラムの観測データとして順次に取得する観測データ取得部と、前記観測データ取得部によって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成部とを備えるものである。 The spectral measuring apparatus of the present invention includes an interferometer that generates interference light by combining the reference light and the signal light by changing the phase of the signal light with respect to the reference light in a non-linear manner in time, and the interference light from the interferometer. The observation data acquisition unit that sequentially acquires the intensity of the observation light obtained by irradiating the measurement target with the observation data of the time domain interferogram at regular time intervals, and the observation acquired by the observation data acquisition unit. It is provided with a spectral spectrum generator that obtains a spectrum that fits the data and minimizes the L1 norm or L0 norm as a spectrum obtained by converting a time domain interferogram into a frequency domain.
本発明の分光測定方法は、参照光に対する信号光の位相を時間的に非線形に変化させて参照光と信号光とを合波した干渉光を生成する干渉光生成ステップと、前記干渉光生成ステップで生成される干渉光を測定対象に照射することにより得られる観測光の強度を一定の時間間隔で時間領域インターフェログラムの観測データとして順次に取得する観測データ取得ステップと、前記観測データ取得ステップによって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成ステップとを有するものである。 The spectral measurement method of the present invention includes an interference light generation step of generating interference light in which the reference light and the signal light are combined by changing the phase of the signal light with respect to the reference light in a non-linear manner, and the interference light generation step. An observation data acquisition step in which the intensity of the observation light obtained by irradiating the measurement target with the interference light generated in the above is sequentially acquired as observation data of the time region interferogram at regular time intervals, and the observation data acquisition step. It has a spectral spectrum generation step of obtaining a spectrum that fits the observation data acquired by the above and minimizes the L1 norm or the L0 norm as a spectrum obtained by converting a time region interferogram into a frequency region.
本発明によれば、繰り返し生成されるパルス対におけるポンプ光に対するプローブ光の遅延時間を時間的に非線形に変化させ、プローブ光が照射されるごとに測定対象から得られる観測光の強度を取得した時間領域インターフェログラムの観測データから、時間領域インターフェログラムを周波数領域に変換したスペクトルとしてL1ノルムまたはL0ノルムを最小化するものを求めるようにしたので、標本化定理で定まる間隔よりも大きなサンプリング間隔で少ないサンプリング数でサンプリングしても適切なスペクトルが得られるため、測定速度を向上することができる。 According to the present invention, the delay time of the probe light with respect to the pump light in the repeatedly generated pulse pair is changed in a non-linear manner in time, and the intensity of the observed light obtained from the measurement target is obtained each time the probe light is irradiated. From the observation data of the time domain interferogram, we tried to find the spectrum that minimizes the L1 norm or L0 norm as the spectrum converted from the time domain interferogram to the frequency domain, so the sampling is larger than the interval determined by the sampling theorem. Since an appropriate spectrum can be obtained even if sampling is performed with a small number of samplings at intervals, the measurement speed can be improved.
本発明によれば、参照光に対する信号光の位相を時間的に非線形に変化させて参照光と信号光とを合波した干渉光を照射した測定対象から得られる観測光の強度を一定の時間間隔で取得した時間領域インターフェログラムの観測データから時間領域インターフェログラムを周波数領域に変換したスペクトルとしてL1ノルムまたはL0ノルムを最小化するものを求めるようにしたので、標本化定理で定まる間隔よりも大きなサンプリング間隔で少ないサンプリング数でサンプリングしても適切なスペクトルが得られるため、測定速度を向上することができる。 According to the present invention, the intensity of the observed light obtained from the measurement target irradiated with the interference light obtained by combining the reference light and the signal light by changing the phase of the signal light with respect to the reference light in a non-linear manner for a certain period of time. Since the spectrum that minimizes the L1 norm or L0 norm is obtained as the spectrum obtained by converting the time domain interferogram into the frequency domain from the observation data of the time domain interferogram acquired at intervals, the interval determined by the sampling theorem is used. However, since an appropriate spectrum can be obtained even when sampling with a small number of samplings at a large sampling interval, the measurement speed can be improved.
図1において、分光測定装置10は、コヒーレント反ストークスラマン散乱分光法により分子振動スペクトルを得るのである。すなわち、この分光測定装置10では、ポンプ光とプローブ光とからなるパルス対を測定対象TGに照射することによって測定対象TGから放出される観測光としてのアンチストークス光の時間領域インターフェログラムをサンプリングし、その結果から分子振動スペクトルを得る。なお、この例では、ポンプ光は、測定対象TGに対してストークス光としても作用する。
In FIG. 1, the
分光測定装置10は、パルス対生成部12、このパルス対生成部12に含まれる光源13、補償器14、ロングパスフィルタ15、第1対物レンズ16、第2対物レンズ17、ショートパスフィルタ18、光検出器21、ローパスフィルタ22及び分光スペクトル生成部としてのコンピュータ23を有する。また、分光測定装置10は、ポンプ光に対するプローブ光の遅延時間を特定するための遅延特定部25を備えている。
The
光源13は、コヒーレント性が高く広帯域のスペクトルを有する非常に短いパルス幅の光パルスを一定の繰り返し周波数で出力する、例えば超短パルスレーザが用いられている。この例では、光源13として、フェムト秒レーザを用いており、中心周波数が792nm、帯域幅が47nm、パルス幅が17fsの光パルスを75MHzの繰り返し周波数で放出する。なお、光源13としては、繰り返し周波数の高い種々のパルスレーザ(固体レーザ、ファイバーレーザ、色素レーザ等)を用いることができる。
As the
パルス対生成部12は、光源13の他に、ビームスプリッタ27、第1アームA1及び第2アームA2を有する。ビームスプリッタ27は、光源13からの光パルスをポンプ光とプローブ光に分割し、ポンプ光を第1アームA1に、プローブ光を第2アームA2にそれぞれ射出する。また、ビームスプリッタ27は、第1アームA1及び第2アームA2でそれぞれ反射して戻ってくるポンプ光とプローブ光を同一の光路上に射出する。
The pulse
ビームスプリッタ27として、例えば、偏光無依存型のものが用いられており、光源13からの光パルスの一部をポンプ光として第1アームA1に向けて反射し、残りの光パルスをプローブ光として第2アームA2に向けて透過する。また、ビームスプリッタ27は、第1アームA1からのポンプ光を透過し、第2アームA2からのプローブ光を反射することで、ポンプ光及びプローブ光を同一の経路に射出する。ポンプ光とプローブ光とは、いずれも光源13からの光パルスと同様の光パルスである。ビームスプリッタ27としては、偏光無依存型ビームスプリッタに代えて、偏光ビームスプリッタを用いてもよい。
As the
第1アームA1は、第1ミラー31を有している。第1ミラー31は、固定された平面鏡であり、その鏡面がビームスプリッタ27で反射されるポンプ光の光路と垂直になるように配されている。これにより、第1アームA1では、ポンプ光が垂直に第1ミラー31に入射し、第1ミラー31で反射されたポンプ光が、同じ経路を逆方向に通ってビームスプリッタ27に入射する。
The first arm A1 has a
光路長変化部としての第2アームA2は、プローブ光をポンプ光に対して遅延させるとともに、繰り返し生成されるプローブ光の遅延時間の時間的な変化が非線形になるように、プローブ光の光路長を時間の変化に対して非線形に変化させる。第2アームA2は、このように、プローブ光の光路長を増減することで、ポンプ光に対するプローブ光の遅延時間を擬似ランダムに変化させる。 The second arm A2 as the optical path length changing portion delays the probe light with respect to the pump light, and the optical path length of the probe light is so that the time change of the delay time of the repeatedly generated probe light becomes non-linear. Is changed non-linearly with respect to the change of time. In this way, the second arm A2 changes the delay time of the probe light with respect to the pump light in a pseudo-random manner by increasing or decreasing the optical path length of the probe light.
この例における第2アームA2は、レゾナントスキャナ32、曲面ミラー33及び第2ミラー34を有している。第2アームA2では、ビームスプリッタ27からのプローブ光がレゾナントスキャナ32、曲面ミラー33及び第2ミラー34で順次に反射された後、これと同じ経路を逆方向に通ってビームスプリッタ27に入射するように構成されている。
The second arm A2 in this example has a
レゾナントスキャナ32は、スキャニングミラー32a、回転軸32b、駆動部32cを有する。スキャニングミラー32aは、回転軸32bを介して駆動部32cで駆動され、所定の周波数で回転軸32bを中心に共振振動する。これにより、スキャニングミラー32aを所定の振れ角で周期運動させ、第2アームA2におけるプローブ光の光路長を変化させる。レゾナントスキャナ32は、例えばスキャニングミラー32aが12kHzで振動するが、振動数は特に限定されない。なお、レゾナントスキャナ32としては、他の形式のレゾナントスキャナを用いてもよい。
The
上記スキャニングミラー32aは、共振振動によって、図1中に二点鎖線で示すように、プローブ光の入射部分がビームスプリッタ27に最も近づいた第1揺動位置と、図1中に実線で示すように、この第1揺動位置から所定の振れ角だけ回転して第1揺動位置よりも入射部分がビームスプリッタ27から離れた第2揺動位置との間で周期的に揺動する。このスキャニングミラー32aの揺動角度(回転角度)に応じて、プローブ光が反射される方向が変化する。揺動角度は、回転軸32bを中心としてスキャニングミラー32aが回転したときの、例えば第1揺動位置からの回転角度である。なお、この例では、スキャニングミラー32aの共振振動による変位が大きい部分すなわちスキャニングミラー32aの回転軸32bから距離が大きい外縁近傍部分にプローブ光が入射するように、レゾナントスキャナ32が位置決めされている。
The
この例では、曲面ミラー33と第2ミラー34とによりミラー部が構成される。曲面ミラー33は、焦点距離fを有する凹面鏡であり、第2ミラー34は、平面鏡である。これらの曲面ミラー33及び第2ミラー34は、固定されている。スキャニングミラー32aと曲面ミラー33の間隔と、曲面ミラー33と第2ミラー34の間隔とは、それぞれ曲面ミラー33の焦点距離fと同じにしてある。曲面ミラー33には、スキャニングミラー32aからのプローブ光が入射する。曲面ミラー33と第2ミラー34とは、スキャニングミラー32aの揺動角度によらず、スキャニングミラー32aで反射されるプローブ光が曲面ミラー33で反射されて第2ミラー34に垂直に入射するように調整されている。これにより、第2アームA2では、第2ミラー34で反射されたプローブ光が、入射時と同じ経路を逆方向に通って、ビームスプリッタ27に入射する。
In this example, the
第2アームA2におけるプローブ光の光路長は、スキャニングミラー32aが第1揺動位置にあるときに、第1アームA1におけるポンプ光の光路長と等しく、第1揺動位置から第2揺動位置に向かって揺動角度が増大するのにしたがって漸増する。第1アームA1におけるポンプ光の光路長は、ビームスプリッタ27で反射されてからビームスプリッタ27に入射するまでのポンプ光の経路の長さである。また、第2アームA2におけるプローブ光の光路長は、ビームスプリッタ27を透過してからビームスプリッタ27に入射するまでのプローブ光の経路の長さである。
The optical path length of the probe light in the second arm A2 is equal to the optical path length of the pump light in the first arm A1 when the
この例では、スキャニングミラー32aが第1揺動位置と第2揺動位置との間で変化したときのプローブ光の光路長の変化幅は1mm程度である。プローブ光の光路長の変化幅は、スキャニングミラー32aに対するプローブ光の入射位置を変えることで調整できる。
In this example, the change width of the optical path length of the probe light when the
第2アームA2では、上記のような構成により、スキャニングミラー32aが第1揺動位置から第2揺動位置に向かって揺動角度が増大する過程では、第2アームA2におけるプローブ光の光路長が増大してポンプ光の光路長との差が漸増する。スキャニングミラー32aが第2揺動位置から第1揺動位置に向かって揺動角度が減少する過程では、第2アームA2におけるプローブ光の光路長が増大してポンプ光の光路長との差が漸減する。このように、プローブ光の光路長は、スキャニングミラー32aの周期的に揺動によって、周期的にプローブ光の光路長が増減する。
In the second arm A2, the optical path length of the probe light in the second arm A2 is in the process in which the swing angle of the
第1アームA1からビームスプリッタ27に入射するポンプ光及び第2アームA2からビームスプリッタ27に入射するプローブ光は、ビームスプリッタ27を透過してまたはビームスプリッタ27で反射されて、同一の経路に射出される。このようにして、パルス対生成部12は、同軸上にポンプ光とプローブ光が並び、プローブ光がポンプ光に遅延したコリニアなパルス対を生成する。
The pump light incident on the
なお、この例では、第2アームA2に曲面ミラー33を用いているが、この曲面ミラー33に代えて集光レンズを用いてもよい。この場合、スキャニングミラー32aで反射されたプローブ光が、スキャニングミラー32aの揺動角度によらず集光レンズを介して第2ミラー34に垂直に入射するように配置する。また、この場合、集光レンズと同等のレンズを第1アームA1の光路上に配置するのがよい。
In this example, the
分光測定装置10は、後述するように圧縮センシングの手法を利用することによって、アンチストークス光の時間領域インターフェログラムを、標本化定理で必要とされるサンプリング周波数よりも低いサンプリング周波数でサンプリングし、標本化定理によって必要とされるサンプリング間隔よりも疎に測定された時間領域インターフェログラムから分子振動スペクトルを再構成する。
The
光パルス(パルス対)の出力ごとに、アンチストークス光の時間領域インターフェログラムのサンプリングを行うので、時間領域インターフェログラムをサンプリングするサンプリング周波数と、光源13による光パルスの繰り返し周波数とは一致する。したがって、光源13による光パルスの繰り返し周波数を、標本化定理で必要とされるサンプリング周波数よりも低くできる。換言すれば、分光測定装置10は、サンプリング周波数(=光パルスの繰り返し周波数)を同じにしながら、そのサンプリング周波数に基づいて標本化定理から決まる分子振動スペクトルの波数よりもの大きな波数の分子振動スペクトルまで測定することができ、分子振動スペクトルの波数の上限を拡張できる。
Since the time domain interferogram of the anti-Stokes light is sampled for each output of the optical pulse (pulse pair), the sampling frequency for sampling the time domain interferogram matches the repetition frequency of the optical pulse by the
また、別の観点からは、光源13による光パルスの繰り返し周波数を維持しながら、レゾナントスキャナ32のスキャニングミラー32aの速度を高くすることで、サンプリング間隔が標本化定理で必要とされる間隔より大きくなっても分子振動スペクトルを再構成することができることを意味する。したがって、レゾナントスキャナ32のスキャニングミラー32aの速度を高くして測定速度を向上することができる。この例でも、そのようにしている。
From another point of view, by increasing the speed of the
パルス対生成部12からのパルス対は、補償器14へ入射する。補償器14は、パルス対のポンプ光及びプローブ光の各々の群速度分散を補償する。補償器14は、ポンプ光及びプローブ光の各々のパルス幅を、測定対象TGの位置で、すなわち測定対象TGに照射されたときに最小にする。このように、補償器14によってポンプ光及びプローブ光の群速度分散を補償することで、測定対象TGにおいてラマン散乱を生じさせ易くし、より確実に分子振動スペクトルを得る。
The pulse pair from the pulse
この例では、補償器14として、鏡面が互いに平行に対面して配された第1チャープミラー14aと第2チャープミラー14bとからなる、いわゆるチャープミラー対を用いている。補償器14の一端から入射するポンプ光及びプローブ光は、それぞれ第1チャープミラー14aと第2チャープミラー14bとの間で繰り返し反射されながら伝播し、補償器14の他端から射出される。
In this example, as the
なお、補償器14は、ポンプ光及びプローブ光の群速度分散を補償してパルス幅を小さくすることができれば特に限定されず、他の形式の群速度分散を補償できる素子等を用いることができる。例えば、回折格子対、フォトニック結晶ファイバー(Photonic Crystal Fiber: PCF)、チャープFBG(Chirped Fiber Bragg Grating)等を用いてもよい。チャープミラー対を含め、これらは群速度分散すなわち二次の分散までを補償するが、例えばプリズム対やグリズム対等のように高次(三次以上)の分散までを補償する補償器14を用いることも好ましい。
The
補償器14からのパルス対は、ロングパスフィルタ15を経て、第1対物レンズ16に入射する。ロングパスフィルタ15は、ポンプ光及びプローブ光のスペクトルにおける短波長側の裾の光を透過しないように、カットオフ波長が設定されている。このロングパスフィルタ15により、測定対象TGで生じるアンチストークス光と同程度の波長の光をポンプ光及びプローブ光からカットし、測定対象TGからのアンチストークス光を検出し易くしてより確実に分子振動スペクトルを得られるようにしている。なお、この例では、ロングパスフィルタ15のカットオフ波長は、750nmである。
The pulse pair from the compensator 14 passes through the
第1対物レンズ16の焦点位置に測定対象TGが配置される。第1対物レンズ16は、入射するポンプ光及びプローブ光を集光して測定対象TGに照射する。したがって、測定対象TGには、ポンプ光とプローブ光とが順番に照射される。この照射では、プローブ光は、第2アームA2で与えられた遅延時間だけポンプ光から遅れて測定対象TGに照射される。
The measurement target TG is arranged at the focal position of the first
測定対象TGにポンプ光が照射されると、測定対象TGの分子に対してポンプ光に含まれるある周波数の光が物質を励起する実質的なポンプ光となり、このポンプ光と異なるある周波数の光がストークス光となり、これらが測定対象TGに照射されたことと同じになる。この結果、測定対象TGにおいてポンプ光とストークス光の周波数差と同じ振動周波数を有する分子振動が誘起される。 When the TG to be measured is irradiated with pump light, the light of a certain frequency contained in the pump light with respect to the molecule of the TG to be measured becomes a substantial pump light that excites a substance, and the light of a certain frequency different from this pump light. Becomes Stokes light, which is the same as irradiating the TG to be measured. As a result, molecular vibration having the same vibration frequency as the frequency difference between the pump light and the Stokes light is induced in the measurement target TG.
ポンプ光に遅れてプローブ光が測定対象TGに照射されると、観測用のプローブ光と誘起された分子振動とが作用し、プローブ光に含まれる光の周波数がシフトして散乱光が放出される。この散乱光には、プローブ光の周波数に対して、実質的なポンプ光とストークス光とによって誘起された分子振動の周波数だけ周波数が増加したアンチストークス光と、プローブ光の周波数に対して、実質的なポンプ光とストークス光によって誘起された分子振動の周波数だけ周波数が減少したストークス光とが含まれている。 When the probe light is applied to the TG to be measured after the pump light, the probe light for observation and the induced molecular vibration act to shift the frequency of the light contained in the probe light and emit scattered light. To. The scattered light includes anti-Stokes light whose frequency is increased by the frequency of the molecular vibration induced by the substantial pump light and Stokes light with respect to the frequency of the probe light, and substantially the frequency of the probe light. It includes the typical pump light and the Stokes light whose frequency is reduced by the frequency of the molecular vibration induced by the Stokes light.
測定対象TGから放出されるアンチストークス光及びストークス光の強度は、ポンプ光に対するプローブ光の遅延時間に応じて変化する。このため、パルス対生成部12によってプローブ光の遅延時間を変化させたパルス対を繰り返し測定対象TGに照射することによって、アンチストークス光又はストークス光の時間領域インターフェログラムが生成される。ポンプ光及びプローブ光は、広帯域な光パルスであり、広い範囲の波長成分を有している。このため、様々な周波数のポンプ光により、種々の分子に分子振動を誘起し、それら分子振動に対応したアンチストークス光及びストークス光が放出される。このため、観測される時間領域インターフェログラムには、広い波数領域の分子振動スペクトルが含まれる。
The intensities of the anti-Stokes light and the Stokes light emitted from the TG to be measured vary depending on the delay time of the probe light with respect to the pump light. Therefore, the time domain interferogram of anti-Stokes light or Stokes light is generated by repeatedly irradiating the measurement target TG with a pulse pair in which the delay time of the probe light is changed by the pulse
第2対物レンズ17は、第1対物レンズ16の焦点位置を挟んで第1対物レンズ16と対向し、その焦点位置が第1対物レンズ16の焦点位置と一致するように配置されている。測定対象TGからのアンチストークス光及びストークス光は、測定対象TGを透過したパルス対とともに、第2対物レンズ17によってコリメートされ、ショートパスフィルタ18に入射する。
The second
ショートパスフィルタ18は、ロングパスフィルタ15のカットオフ波長以下の短い波長がカットオフ波長に設定されている。パルス対のポンプ光とプローブ光は、ロングパスフィルタ15によって、そのカットオフ波長以下の波長成分の大部分がカットされている。このため、測定対象TGを透過したポンプ光とプローブ光は、そのほとんどがショートパスフィルタ18によってカットされる。
In the
また、ストークス光は、プローブ光の特定の周波数成分が低周波側に周波数シフトしたものであるので、そのほとんどがショートパスフィルタ18でカットされる。これに対して、プローブ光に対して周波数が増加したアンチストークス光は、そのほとんどがショートパスフィルタ18を透過する。このようにして、測定対象TGからの光のうちアンチストークス光のみがショートパスフィルタ18を透過する。この例では、ショートパスフィルタ18のカットオフ波長は、750nmである。
Further, since the Stokes light is a probe light in which a specific frequency component is frequency-shifted to the low frequency side, most of the Stokes light is cut by the
光検出器21は、測定対象TGからのアンチストークス光の強度を検出するものであり、ショートパスフィルタ18を透過したアンチストークス光が入射する。この光検出器21は、アンチストークス光を受光し、そのアンチストークス光の強度を電気信号(以下、光電信号と称する)に変換する。この例では、光検出器21は、特に限定されず、例えばアバランシェフォトダイオードや光電子増倍管PINフォトダイオードを用いたもの等を用いることができる。光電信号は、ローパスフィルタ22を介してコンピュータ23に送られる。ローパスフィルタ22は、光電信号高周波ノイズを除去する。
The
コンピュータ23は、サンプリング部としてのA/D変換器23aを備えている。このA/D変換器23aは、ローパスフィルタ22からの光電信号が入力され、光電信号をサンプリングする。このサンプリングでは、A/D変換器23aは、光電信号をデジタル変換し、得られる観測データを検出時間と共にコンピュータ23のメモリ23bに保存する。この例では、光検出器21、ローパスフィルタ22、A/D変換器23aにより観測データ取得部が構成される。
The
A/D変換器23aは、光源13からの光パルスの出力に同期してサンプリングを行う。すなわち、光パルスの繰り返し周期でサンプリングし、得られる観測データを順次に保存する。これにより、アンチストークス光の時間領域インターフェログラムの強度をサンプリングした複数の観測データが保存される。観測データとともに保存される検出時間は、例えば光源13が光パルスを射出する時刻(ポンプ光が照射される時刻)である。
The A /
ところで、上述のように構成される第2アームA2では、プローブ光の単位時間当たりの光路長の変化量がスキャニングミラー32aの速度に応じて変わる。また、スキャニングミラー32aの速度が変化する。例えば、第1揺動位置及び第2揺動位置の近傍では、スキャニングミラー32aの速度変化が大きく、第1揺動位置と第2揺動位置の中間位置近傍では速度変化が小さい。また、相対的に、第1揺動位置及び第2揺動位置の近傍では、スキャニングミラー32aの速度が小さく、中間位置近傍では速度が大きい。このため、第2アームA2によって、プローブ光の光路長が時間に対して非線形(例えば正弦波状)に変化する。
By the way, in the second arm A2 configured as described above, the amount of change in the optical path length per unit time of the probe light changes according to the speed of the
プローブ光の光路長の時間的な変化が非線形であるため、図2に示すように、ポンプ光に対するプローブ光の遅延時間(Pump-probe delay)は、光路長の時間的な変化と同様に、時間的に非線形に変化し、この例では正弦波状に変化する。現実の時間(Lab time)において、光源13から光パルスを一定の間隔で出力することによって、ポンプ光を基準とするプローブ光の測定対象TGへの照射タイミング及びそれによる測定対象TGからアンチストークス光が射出されるタイミングは時間的に非線形に変化する。このため、例えばスキャニングミラー32aが第1揺動位置から第2揺動位置に回転している期間あるいは第2揺動位置から第1揺動位置に回転している期間における、ポンプ光に対するプローブ光の遅延時間は、ランダムなものとみなすことができる擬似ランダムとなる。なお、スキャニングミラー32aが第1揺動位置から第2揺動位置に回転している期間では、遅延時間は漸増し、第2揺動位置から第1揺動位置に回転している期間では、遅延時間は漸減する。
Since the temporal change of the optical path length of the probe light is non-linear, as shown in FIG. 2, the delay time (Pump-probe delay) of the probe light with respect to the pump light is similar to the temporal change of the optical path length. It changes non-linearly with time, and in this example it changes like a sinusoidal shape. By outputting light pulses from the
一方、A/D変換器23aは、光パルスの出力に同期して一定のサンプリング間隔でサンプリングを行うが、例えばA/D変換器23aのサンプリングタイミング、サンプリング時間(光電信号を読み取っている時間)、ローパスフィルタ22の特性の調整等によって、光検出器21が受光したアンチストークス光の強度に対応した光電信号についてサンプリングする。このため、A/D変換器23aによるサンプリングは、アンチストークス光が放出されたタイミング、すなわちプローブ光が測定対象TGに照射されるタイミングでサンプリングすることと同じである。したがって、各パルス対のポンプ光の照射タイミングを基準にすると、A/D変換器23aによる実質的なサンプリングタイミングは、時間的に非線形に変化し擬似ランダムなものとなり、遅延時間の変化が等間隔ではない擬似ランダムな観測データが取得される。
On the other hand, the A /
上記のように、分光測定装置10では、第2アームA2におけるプローブ光の光路長の時間的な変化を非線形にすることにより、A/D変換器23aによって、アンチストークス光の時間領域インターフェログラムがランダムにサンプリングされ、遅延時間を非線形に変化させた観測データを得る。後述する再構成部23cによる圧縮センシングの手法を用いた分子振動スペクトルの再構成では、アンチストークス光の時間領域インターフェログラムに対するサンプリング結果である観測データの遅延時間がランダムないし疑似ランダムである場合に良好な再構成結果が得られる。
As described above, in the
光スペクトル生成部としてのコンピュータ23の再構成部23cは、時間領域インターフェログラムをサンプリングした観測データから分子振動スペクトルを再構成する。この再構成部23cは、コンピュータ23にインストールした所定のプログラムによって実現されている。
The reconstruction unit 23c of the
パルス対生成部12によって生成されるポンプ光に対するプローブ光の遅延は、スキャニングミラー32aの第1揺動位置からの揺動角度に応じて周期的に変化し、半周期で最大となる。このため、再構成部23cは、例えば、第1揺動位置から第2揺動位置に回転している期間あるいは第2揺動位置から第1揺動位置に回転している期間の1組の観測データを観測ベクトルとしてメモリ23bから取り出し、詳細を後述するように、その観測データを用いて測定対象TGの分子振動スペクトルを求める。求めた分子振動スペクトルは、例えばモニタ等に表示される。
The delay of the probe light with respect to the pump light generated by the pulse
再構成部23cは、分子振動スペクトルを求める際には、観測データに対応する遅延時間を使用する。観測データに対応した遅延時間を算出可能にするために遅延特定部25が設けられている。
The reconstruction unit 23c uses the delay time corresponding to the observed data when obtaining the molecular vibration spectrum. A
遅延特定部25は、光源41、ビームスプリッタ42、43及び光検出器44を有している。光源41は、光源13からの光パルスとは波長域の異なるCW(continuous wave)レーザ光を出力する。例えば、光源41が出力するCWレーザ光の波長は、1064nmである。ビームスプリッタ42は、光源13とパルス対生成部12との間の光パルスの光路上に配置されている。このビームスプリッタ42は、光パルスを透過し、CWレーザ光を反射する。これにより、光源41からのCWレーザ光は、光源13からの光パルスの光路を通ってパルス対生成部12に入射する。
The
パルス対生成部12に入射するCWレーザ光は、光源13からの光パルスと同様に、ビームスプリッタ27によって、第1アームA1を伝播するCWレーザ光と第2アームA2を伝播するCWレーザ光とに分割される。第1アームA1を伝播するCWレーザ光と第2アームA2を伝播するCWレーザ光は、ビームスプリッタ27で合波することにより、これらが互いに干渉した干渉波(以下、CW干渉波と称する)となってパルス対生成部12から出力される。
The CW laser light incident on the pulse
上記のように生成されるCW干渉波は、第1アームA1と第2アームA2との光路長の差に応じて強度が変化する。したがって、CW干渉波のインターフェログラムは、強度が時間に対して周期的に変化する正弦波のような形状となり、その強度変化の周波数が、スキャニングミラー32aの速度の変化に対応して変化する。
The intensity of the CW interference wave generated as described above changes according to the difference in the optical path length between the first arm A1 and the second arm A2. Therefore, the interferogram of the CW interference wave has a shape like a sine wave whose intensity changes periodically with time, and the frequency of the intensity change changes in response to the change in the speed of the
ビームスプリッタ43は、パルス対生成部12と補償器14との間のパルス対の光路上に配置されている。このビームスプリッタ43は、パルス対を透過し、CW干渉波を反射することで、パルス対とCW干渉波とを分離する。ビームスプリッタ43で分離されたCW干渉波は、光検出器44に入射する。光検出器44は、CW干渉波を受光し、CWレーザ干渉波の強度を電気的なCW干渉信号に変換する。
The beam splitter 43 is arranged on the optical path of the pulse pair between the
光検出器44からのCW干渉信号は、コンピュータ23内のA/D変換器23dに入力される。A/D変換器23dは、入力されるCW干渉信号を光源13の光パルスの出力に同期してサンプリングし、CW干渉信号をデジタル変換したCW干渉データを検出時間と共にコンピュータ23のメモリ23bに保存する。これにより、CW干渉波のインターフェログラムの強度をサンプリングした複数のCW干渉データが検出時間によって観測データと対応づけられて保存される。
The CW interference signal from the
なお、遅延特定部25は、ポンプ光に対するプローブ光の遅延時間を特定できるものであればどのような構成でもよい。例えば、ポンプ光とプローブ光とのそれぞれ一部を分割して受光し、それらの受光タイミングから遅延時間を求めてもよい。
The
遅延時間の算出は、測定されたCW干渉信号を解析することによって行う。まず、CW干渉信号を時間軸に関して逆フーリエ変換し、CW干渉信号のスペクトルを得る。得られるスペクトルには周波数原点、及び正と負の周波数領域にそれぞれ1つずつ合計3つのピークを有する。CW干渉信号のスペクトルに対してフィルター関数を施すことによって、正の周波数を有するピーク成分のみを抽出し、抽出したピーク成分にフーリエ変換を行う。これにより、CW干渉信号を、複素数成分を有するインターフェログラムとして復元することができる。得られる複素インターフェログラムの虚部と実部の比の逆正接を計算することで、遅延時間をCWレーザの波長に対する位相として求められる。このように求めた位相に波長を乗じた値を値「2π×c(cは光速)」で除すことで遅延時間を算出する。 The delay time is calculated by analyzing the measured CW interference signal. First, the CW interference signal is inverse-Fourier-transformed with respect to the time axis to obtain the spectrum of the CW interference signal. The resulting spectrum has a frequency origin and a total of three peaks, one in each of the positive and negative frequency domains. By applying a filter function to the spectrum of the CW interference signal, only the peak component having a positive frequency is extracted, and the extracted peak component is subjected to Fourier transform. Thereby, the CW interference signal can be restored as an interferogram having a complex number component. By calculating the inverse tangent of the ratio of the imaginary part to the real part of the obtained complex interferogram, the delay time can be obtained as the phase with respect to the wavelength of the CW laser. The delay time is calculated by dividing the value obtained by multiplying the phase obtained in this way by the wavelength by the value "2π × c (c is the speed of light)".
なお、上記遅延時間の算出の手法は一例であり、CW干渉信号から遅延時間を求めることができるのであればどのような計算手法を用いてもよい。例えば、CW干渉信号が0点を交差する地点をCWレーザの波長分の遅延に相当する部分として遅延を算出してもよい。 The method for calculating the delay time is an example, and any calculation method may be used as long as the delay time can be obtained from the CW interference signal. For example, the delay may be calculated by setting the point where the CW interference signal intersects the 0 points as a portion corresponding to the delay corresponding to the wavelength of the CW laser.
次に再構成部23cによる分子振動スペクトルの再構成について説明する。再構成部23cは、M個の観測データを1組の観測データとして、この1組の観測データから構成されM次元ベクトルとなる観測ベクトル(y)からN次元ベクトルである分子振動スペクトル(x)を再構成する。M、Nは、いずれも1以上の整数であり、M<Nである。この再構成では、分子振動スペクトルが、それに作用する観測行列Φ(M×N行列)を用い、y=Φxで表される線形変換によるものとして、圧縮センシングの手法により、より低次元なM次元ベクトルの観測データから高次元の分子振動スペクトルを復元する。すなわち、分子振動スペクトルはスパース性を有することを事前情報として用い、観測データに適合し、分子振動スペクトルのL1ノルムが最小化される分子振動スペクトルを求める。なお、この例では、分子振動スペクトルのL1ノルムを用いているが、L0ノルムを用いてもよい。 Next, the reconstruction of the molecular vibration spectrum by the reconstruction unit 23c will be described. The reconstruction unit 23c uses M observation data as a set of observation data, and the molecular vibration spectrum (x) which is an N-dimensional vector from the observation vector (y) which is composed of the set of observation data and becomes an M-dimensional vector. To reconstruct. Both M and N are integers of 1 or more, and M <N. In this reconstruction, the molecular vibrational spectrum uses the observation matrix Φ (M × N matrix) that acts on it, and is based on the linear transformation represented by y = Φx. Restore high-dimensional molecular vibrational spectra from vector observations. That is, using the fact that the molecular vibration spectrum has sparseness as prior information, the molecular vibration spectrum that matches the observation data and minimizes the L1 norm of the molecular vibration spectrum is obtained. In this example, the L1 norm of the molecular vibration spectrum is used, but the L0 norm may be used.
上記のように、圧縮センシングの手法によって時間領域インターフェログラムを周波数領域の分子振動スペクトルに再構成する場合、アンチストークス光の時間領域インターフェログラムをサンプリングした観測データに対応する遅延時間がランダムになっていることで良好な結果が得られる。この例では、上述のように、プローブ光を遅延させる第2アームA2の光路長を非線形に変化させて、そのプローブ光の照射によるアンチストークス光をサンプリングしているため、アンチストークス光の時間領域インターフェログラムをサンプリングした観測データに対応する遅延時間が疑似ランダムでありランダムとみなせる。このように、観測データの遅延時間をランダムにするための構成としては、レゾナントスキャナ32の特性を利用しているので、それ以外に特別な構成を設ける必要はない。レゾナントスキャナ32を用いた分光測定装置10の構成は、構成部材を少なくする上で有利である。
As described above, when the time domain interferogram is reconstructed into the molecular vibration spectrum of the frequency domain by the compression sensing method, the delay time corresponding to the observed data obtained by sampling the time domain interferogram of anti-Stokes light is randomly generated. Good results can be obtained. In this example, as described above, the optical path length of the second arm A2 that delays the probe light is non-linearly changed, and the anti-Stokes light generated by the irradiation of the probe light is sampled. The delay time corresponding to the observed data obtained by sampling the interferogram is pseudo-random and can be regarded as random. As described above, since the characteristic of the
以下に、分子振動スペクトルを求める手順の一例を説明する。以下の手順では、式(q)を満たすxRとして分子振動スペクトルを算出するための具体的な手続きを示す。なお、分子振動スペクトルを求める手順は、これに限定されるものではない。 An example of the procedure for obtaining the molecular vibration spectrum will be described below. In the following procedure, a specific procedure for calculating the molecular vibration spectrum as x R satisfying the equation (q) is shown. The procedure for obtaining the molecular vibration spectrum is not limited to this.
(1) 1組の観測データy(j)と、これに対応する、すなわち検出時間が同じ1組のCW干渉データw(j)をメモリ23bから取り出す。
(2) 各観測データy(j)について、対応するCW干渉データw(j)を用いて、ポンプ光に対するプローブ光の遅延時間τ(j)を上述のように算出する。
(3) 式(a)により、観測ベクトル(1組の観測データ)yに対応する観測行列Φを構成する。式(a)中の行列Aは、フーリエ変換行列Wのうち、特定の行成分のみを抽出する行列、すなわち観測ベクトル中のまばらに観測した信号成分のみを用いて分子振動スペクトルを計算できるようにする行列であり、0と1のみを要素として有するM行N列のバイナリ行列である。フーリエ変換行列Wは、不等間隔離散フーリエ変換行列であり、手順(2)で求めた不等間隔(疑似ランダム)な1組の遅延時間で構成されるデータに対してフーリエ変換を施す行列であり、各列ベクトルが異なる周期で振動する複素指数関数となっているN行N列の行列である。
(1) A set of observation data y (j) and a set of CW interference data w (j) corresponding to the observation data y (j) having the same detection time are taken out from the memory 23b.
(2) For each observation data y (j), the delay time τ (j) of the probe light with respect to the pump light is calculated as described above using the corresponding CW interference data w (j).
(3) The observation matrix Φ corresponding to the observation vector (a set of observation data) y is constructed by the equation (a). The matrix A in the equation (a) is a matrix that extracts only specific row components from the Fourier transform matrix W, that is, the molecular vibration spectrum can be calculated using only the sparsely observed signal components in the observation vector. It is an M-row and N-column binary matrix having only 0 and 1 as elements. The Fourier transform matrix W is an unequal-interval discrete Fourier transform matrix, and is a matrix that performs a Fourier transform on data composed of a set of unequal-interval (pseudo-random) delay times obtained in step (2). It is an N-by-N matrix in which each column vector is a complex exponential function that vibrates in different periods.
(4) N次元ベクトルである分子振動スペクトル(x)の初期値(x(0))として0ベクトルを代入する(式(b))。なお、x(k)(k=0,1,2,3・・・)は、分子振動スペクトルを推定する後述の繰り返し計算におけるk番目の計算ステップにおける分子振動スペクトル(推定値)である。 (4) Substitute the 0 vector as the initial value (x (0) ) of the molecular vibration spectrum (x) which is an N-dimensional vector (Equation (b)). Note that x (k) (k = 0, 1, 2, 3 ...) Is the molecular vibration spectrum (estimated value) in the k-th calculation step in the iterative calculation described later for estimating the molecular vibration spectrum.
(5) ステップサイズγを式(c)により算出する。ステップサイズγは、後述の繰り返し計算において、次の(k+1)番目の計算ステップにおける分子振動スペクトルx(k+1)を推定するにあたり、現在のk番目の計算ステップにおいて推定されている分子振動スペクトルx(k)と観測行列Φから計算される勾配方向へとベクトルを加算するが、その加算するベクトルの大きさを特徴づける量である。なお、式(c)の右辺の分母は、無限大ノルムである。 (5) The step size γ is calculated by the formula (c). The step size gamma, in iterative calculation described later, the next (k + 1) -th computation Upon estimating the molecular vibration spectra x (k + 1) in step is estimated in the current k-th calculation step molecular vibration spectrum x ( The vector is added in the gradient direction calculated from k) and the observation matrix Φ, and it is an amount that characterizes the magnitude of the added vector. The denominator on the right side of equation (c) is an infinite norm.
(6)上記の式(q)において、推定される分子振動スペクトルのスパース性を特徴づけるハイパーパラメータであるλを再構成部23cに入力する。この例では、予め決めておいたハイパーパラメータλを再構成部23cに入力するが、データの再構成の精度などを指標として自動的に算出されたハイパーパラメータλを自動入力することもできる。例えば異なるハイパーパラメータλに対して、観測ベクトルを訓練データとテストデータに分割し、交差検定を用いて推定された分子振動スペクトルの妥当性を評価し、その妥当性が最大となるようなハイパーパラメータλの値を用いてもよい。 (6) In the above equation (q), λ, which is a hyperparameter that characterizes the sparsity of the estimated molecular vibration spectrum, is input to the reconstruction unit 23c. In this example, the hyperparameter λ determined in advance is input to the reconstruction unit 23c, but the hyperparameter λ automatically calculated by using the accuracy of data reconstruction as an index can also be automatically input. For example, for different hyperparameters λ, the observation vector is divided into training data and test data, the validity of the molecular vibration spectrum estimated using cross-validation is evaluated, and the hyperparameters that maximize the validity. You may use the value of λ.
(7) 式(d)により、値r(k)を算出する。r(k)は、x(k)を勾配方向にステップサイズγで特徴づけられる大きさだけ変化させたベクトルである。式(d)中のΦHは、観測行列Φのエルミート共役行列である。 (7) The value r (k) is calculated by the equation (d). r (k) is a vector obtained by changing x (k ) in the gradient direction by a magnitude characterized by the step size γ. Φ H in the equation (d) is the Hermitian conjugate matrix of the observation matrix Φ.
(8) 式(e)により、xl (k+1)を求める。xl (k+1)は、(k+1)番目の計算ステップにおける分子振動スペクトルx(k+1)の第l成分である。また、式(e)中の右辺式(e)の右辺における関数Sは、信号(分子振動スペクトル)の疎性を促進し、繰り返し計算の収束先を式(q)の最適解にするための軟判定しきい値関数である。また、rl (k)は、r(k)の第l成分である。 (8) Obtain x l (k + 1) by the equation (e). x l (k + 1) is the first component of the molecular vibration spectrum x (k + 1) in the (k + 1) th calculation step. Further, the function S on the right side of the right-hand side equation (e) in the equation (e) promotes the sparseness of the signal (molecular vibration spectrum) and makes the convergence destination of the iterative calculation the optimum solution of the equation (q). It is a soft judgment threshold function. Further, r l (k) is the l-th component of r (k).
(9) 値aを十分に小さい値として、「|x(k+1)−x(k)|2<a」の条件を判定する。すなわち、k番目の計算ステップで推定した分子振動スペクトル(x(k))とk+1番目の計算ステップで推定した分子振動スペクトル(x(k+1))の差のL2ノルムが所定の値aより小さいか否かを判定する。この条件を満たしていない場合には、値kをインクリメント(K←K+1)して、上記手順(7)、手順(8)を行い、再度この手順(9)により条件を判定する。条件を満たしている場合には、|x(k+1)|を求めるべき分子振動スペクトルとする。 (9) With the value a being a sufficiently small value, the condition of "| x (k + 1) -x (k) | 2 <a" is determined. That is, is the L2 norm of the difference between the molecular vibration spectrum (x (k) ) estimated in the kth calculation step and the molecular vibration spectrum (x (k + 1)) estimated in the k + 1th calculation step smaller than the predetermined value a? Determine if not. If this condition is not satisfied, the value k is incremented (K ← K + 1), the above procedure (7) and procedure (8) are performed, and the condition is determined again by this procedure (9). When the condition is satisfied, | x (k + 1) | is used as the molecular vibrational spectrum to be obtained.
以上のようにして、分子振動スペクトルが再構成される。なお、再構成では、繰り返し生成されるパルス対のポンプ光に対するプローブ光の遅延時間が時間に対して線形的に変化しているようにする時間軸の補正等は不要である。 As described above, the molecular vibration spectrum is reconstructed. In the reconstruction, it is not necessary to correct the time axis so that the delay time of the probe light with respect to the pump light of the repeatedly generated pulse pair changes linearly with time.
図3は、トルエンを測定対象TGとした場合におけるアンチストークス光の時間領域インターフェログラムを示している。図3において、黒丸は、分光測定装置10において、スキャニングミラー32aを共振振動させた場合に取得された観測データを示している。一方、実線で示す時間領域インターフェログラムは、分光測定装置10と同じ構成の分光測定装置において、レゾナントスキャナのスキャニングミラーを分光測定装置10の1/15の速度で共振振動させた場合に取得された観測データ(以下、参考観測データと称する)を便宜的に連続した線で接続して描いたものである。観測データは、標本化定理を満足するサンプリング間隔よりも大きなサンプリング間隔となっており、参考観測データは、アンチストークス光の時間領域インターフェログラムに対して、標本化定理を十分に満足するサンプリング間隔で取得されたものになっている。なお、図3に示されるグラフの横軸は、ポンプ光の光路長に対するプローブ光の光路長の差である。
FIG. 3 shows a time domain interferogram of anti-Stokes light when toluene is used as the measurement target TG. In FIG. 3, black circles indicate observation data acquired when the
観測データを用いた上述の手順により、図4において符号G1に示す分子振動スペクトルが再構成された。一方、参考観測データを用いて、データ解析手法としてのフーリエ変換を行うことにより、符号G3に示す分子振動スペクトルが得られた。フーリエ変換する場合には、ポンプ光に対するプローブ光の遅延時間の変化する間隔が等間隔となるように補正を行った。これらの結果より、圧縮センシングの手法を用いた分光測定装置10により、適切な分子振動スペクトルが再構成できることがわかる。なお、比較として、観測データを用いて、参考観測データと同様にフーリエ変換を行った場合には、符号G2に示す分子振動スペクトルとなり、適切な分子振動スペクトルが得られなかった。
The molecular vibration spectrum shown by reference numeral G1 in FIG. 4 was reconstructed by the above-mentioned procedure using the observation data. On the other hand, the molecular vibration spectrum shown by reference numeral G3 was obtained by performing the Fourier transform as a data analysis method using the reference observation data. In the case of Fourier transform, correction was performed so that the interval at which the delay time of the probe light changes with respect to the pump light changes at equal intervals. From these results, it can be seen that an appropriate molecular vibration spectrum can be reconstructed by the
上記の例では、光路長増減部としての第2アームは、レゾナントスキャナを用いてプローブ光の光路長を増減しているが、光路長増減部は、プローブ光の光路長の時間的な変化が非線形すなわちポンプ光に対するプローブ光の遅延時間を非線形に変化させ、遅延時間が擬似ランダムとみなせるような構成であればよい。また、遅延時間の擬似ランダム化については、厳密でなくてよい。 In the above example, the second arm as the optical path length increasing / decreasing part increases / decreases the optical path length of the probe light by using a resonant scanner, but the optical path length increasing / decreasing part causes a temporal change in the optical path length of the probe light. The configuration may be non-linear, that is, the delay time of the probe light with respect to the pump light may be changed non-linearly so that the delay time can be regarded as pseudo-random. Also, the pseudo-randomization of the delay time does not have to be strict.
レゾナントスキャナに代えて直線的に移動するミラーによりプローブ光の遅延時間を非線形に変化させる例を図5に示す。この例では、第2アームA2には、平面鏡であるミラー51が設けられている。ミラー51には、ビームスプリッタ27で分割されたプローブ光が垂直に入射し、ミラー51で反射したプローブ光がビームスプリッタ27に戻るように構成されている。ミラー51は、プローブ光の進行方向に沿って移動自在に設けられており、駆動部51aによって、ビームスプリッタ27に近づいた第1位置とこの第1位置よりもビームスプリッタ27から離れた第2位置との間で往復動される。駆動部51aは、第1位置と第2位置との間での移動速度の時間的な変化が非線形となるようにミラー51を移動する。これにより、ポンプ光に対するプローブ光の遅延時間が非線形に変化させ、遅延時間が擬似ランダムとする。
FIG. 5 shows an example in which the delay time of the probe light is changed non-linearly by a mirror that moves linearly instead of the resonant scanner. In this example, the second arm A2 is provided with a
上記では、コヒーレントアンチストークスラマン散乱(CARS)で生じたアンチストークス光から分子振動スペクトルを得る分光測定装置について説明しているが、測定対象から放出された散乱光としてのストークス光を用いて時間領域インターフェログラムを作成し、当該時間領域インターフェログラムから分子振動スペクトルを取得するようにしてもよい。 In the above, the spectroscopic measurement device that obtains the molecular vibration spectrum from the anti-Stokes light generated by coherent anti-Stoke Raman scattering (CARS) is described, but the time region using the Stokes light as the scattered light emitted from the measurement target. An interferogram may be created and the molecular vibrational spectrum may be obtained from the time region interferogram.
本発明は、物質を励起させるポンプ光と観測するプローブ光とからなるパルス対を、ポンプ光に対するプローブ光の遅延時間を変化させながら繰り返し生成して測定対象に照射し、プローブ光の照射で得られる時間領域インターフェログラムからスペクトルを得るいわゆるポンプ−プローブ分光法を用いた各種の分光測定装置に適用することができる。例えば、測定対象の吸光特性に基づく分子振動スペクトルの取得にも利用できる。観測光は、散乱光に限らず、透過光、反射光でもよい。また、本発明は、テラヘルツ時間領域分光法等にも利用することができる。 In the present invention, a pulse pair consisting of a pump light that excites a substance and a probe light to be observed is repeatedly generated while changing the delay time of the probe light with respect to the pump light and irradiated to the measurement target, and obtained by irradiation with the probe light. It can be applied to various spectroscopic measuring devices using so-called pump-probe spectroscopy for obtaining a spectrum from a time region interferogram. For example, it can also be used to acquire a molecular vibration spectrum based on the absorption characteristics of the measurement target. The observed light is not limited to scattered light, but may be transmitted light or reflected light. The present invention can also be used for terahertz time region spectroscopy and the like.
また、パルス対生成部の光源を赤外領域の光源へと変更することにより、フーリエ変換赤外分光光度計(FT−IR)と同様に、赤外線吸収スペクトルを測定することが可能である。図6は、測定対象に赤外線を照射し、赤外線吸収スペクトルを得る分光測定装置の例を示す。なお、以下に詳細を説明する他は、上記の実施形態と同様であるので、実質的に同じ構成部材には、同一の符号を付してその詳細な説明を省略する。 Further, by changing the light source of the pulse pair generation unit to a light source in the infrared region, it is possible to measure the infrared absorption spectrum in the same manner as the Fourier transform infrared spectrophotometer (FT-IR). FIG. 6 shows an example of a spectroscopic measuring device that irradiates a measurement target with infrared rays to obtain an infrared absorption spectrum. In addition, since it is the same as the above-described embodiment except that the details will be described below, substantially the same constituent members are designated by the same reference numerals and the detailed description thereof will be omitted.
分光測定装置50は、干渉計52、この干渉計52に含まれる光源53、光検出器21、ローパスフィルタ22、コンピュータ23及び遅延特定部25から構成される。光源53は、広帯域のスペクトルを有する赤外光を連続的に出力するものが用いられている。光源53としては、量子カスケードレーザーや光パラメトリック増幅器、光パラメトリック発振器などのコヒーレント光源や、高輝度セラミック光源、タングステン・ヨウ素ランプ、グローバ光源、発光ダイオードのようなインコヒーレントな光源を用いることができる。干渉計52は、光源53が異なる他は、上記実施形態のパルス対生成部12(図1参照)と同じ構成である。
The
光源53からのレーザ光は、ビームスプリッタ27によって、参照光と信号光とに分割される。この例では、ビームスプリッタ27で反射された光が参照光となって、この参照光が第1アームA1を伝播し、第1ミラー31で反射されてビームスプリッタ27に再入射する。また、ビームスプリッタ27を透過した光が信号光となって、この信号光が第2アームA2を伝播し、レゾナントスキャナ32のスキャニングミラー32a、曲面ミラー33及び第2ミラー34で順次に反射された後、これと同じ経路を逆方向に通ってビームスプリッタ27に再入射する。ビームスプリッタ27に再入射する参照光と信号光とは合波されて干渉光として射出される。
The laser light from the
干渉計52では、スキャニングミラー32aの揺動によって、第2アームA2における信号光の光路長が第1アームA1における参照光の光路長に対して増減される。これにより、参照光に対する信号光の位相が変化する干渉光が生成される。この例では、スキャニングミラー32aの揺動によって、信号光の光路長を参照光の光路長に対して、時間的に非線形に変化させて、参照光に対する信号光の位相を時間的に非線形に変化させる。参照光に対する信号光の位相を時間的に非線形に変化させて干渉光を生成する構成は、上記実施形態と同様に特に限定されない。干渉計52で生成される干渉光は、測定対象TGに照射される。
In the
光検出器21は、測定対象TGからの観測光としての透過光の強度を検出するものであり、干渉光が測定対象TGを透過した透過光が入射する。光検出器21は、透過光を受光し、その透過光の強度を光電信号に変換する。
The
コンピュータ23のA/D変換器23aは、光電信号を順次にサンプリングし、得られる観測データを検出時間と共にコンピュータ23のメモリ23bに保存する。これにより、透過光の時間領域インターフェログラムの強度をサンプリングした複数の観測データを得る。このA/D変換器23aは、一定のサンプリング間隔で光電信号をサンプリングするようにコンピュータ23によって制御される。上記実施形態と同様に、スキャニングミラー32aによる参照光に対する信号光の遅延時間を時間的に非線形に変化させているから、遅延時間を時間的に非線形に変化させた複数の観測データがA/D変換器23aによって順次にサンプリングされて保存される。
The A /
分光測定装置50においても、圧縮センシングの手法を利用して、透過光の時間領域インターフェログラムを、周波数領域の赤外線吸収スペクトルに再構成する。このため、A/D変換器23aによるサンプリング間隔は、標本化定理によって必要とされるサンプリング間隔よりも疎であってよい。したがって、レゾナントスキャナ32のスキャニングミラー32aの速度を高くして測定速度を向上することができる。
The
また、A/D変換器23dは、A/D変換器23aによるサンプリングと同じタイミングで、遅延特定部25から入力されるCW干渉信号をサンプリングし、得られるCW干渉データを検出時間と共にメモリ23bに保存する。なお、遅延特定部25の光源41は、光源53の赤外光の波長域の異なるレーザ光を出力するものが用いられる。
Further, the A /
再構成部23cは、メモリ23bに保存されている観測データ及びCW干渉データを用いて、上記の実施形態と同じ手順で赤外線吸収スペクトルを復元する。すなわち、赤外線吸収スペクトルがスパース性を有することを事前情報として用い、観測データに適合し、赤外線吸収スペクトルのL1ノルムが最小化される赤外線吸収スペクトルを求めることにより、赤外線吸収スペクトルを復元する。なお、このときに、CW干渉データから求められる遅延時間は、参照光に対する信号光の位相の遅れ量に相当する。 The reconstruction unit 23c restores the infrared absorption spectrum by the same procedure as in the above embodiment using the observation data and the CW interference data stored in the memory 23b. That is, the infrared absorption spectrum is restored by using the fact that the infrared absorption spectrum has sparseness as prior information and obtaining an infrared absorption spectrum that matches the observation data and minimizes the L1 norm of the infrared absorption spectrum. At this time, the delay time obtained from the CW interference data corresponds to the amount of phase delay of the signal light with respect to the reference light.
上記では、測定対象からの透過光を観測光とした場合について説明したが、金属上に配された測定対象から反射されてくる反射光を観測光としたり、測定対象に密着させた赤外光に対して透明な高屈折率媒質に対して干渉光を入射させ、測定対象と高屈折率媒質との界面で全反射する全反射光を観測光としたりしてもよい。 In the above, the case where the transmitted light from the measurement target is used as the observation light has been described, but the reflected light reflected from the measurement target arranged on the metal can be used as the observation light or the infrared light which is brought into close contact with the measurement target. Interfering light may be incident on a transparent high refractive index medium, and the total reflected light totally reflected at the interface between the measurement target and the high refractive index medium may be used as the observation light.
10、50 分光測定装置
12 パルス対生成部
13、53 光源
21 光検出器
23a A/D変換器
23c 再構成部
27 ビームスプリッタ
31 第1ミラー
32 レゾナントスキャナ
32a スキャニングミラー
34 第2ミラー
52 干渉計
A1 第1アーム
A2 第2アーム
10, 50
Claims (6)
前記パルス対生成部からのパルス対が繰り返し照射される測定対象からプローブ光が照射されるごとに得られる観測光の強度を時間領域インターフェログラムの観測データとして取得する観測データ取得部と、
前記観測データ取得部によって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成部と
を備えることを特徴とする分光測定装置。 A pulse pair generator that generates a pulse pair consisting of a pump light and a probe light that is delayed with respect to the pump light, and changes the optical path length of the probe light in a non-linear manner with respect to the optical path length of the pump light.
An observation data acquisition unit that acquires the intensity of the observation light obtained each time the probe light is irradiated from the measurement target to which the pulse pair from the pulse pair generation unit is repeatedly irradiated as observation data of the time domain interferogram.
It is provided with a spectral spectrum generation unit that obtains a spectrum that matches the observation data acquired by the observation data acquisition unit and minimizes the L1 norm or the L0 norm as a spectrum obtained by converting the time domain interferogram into the frequency domain. A featured spectroscopic measuring device.
一定の繰り返し周期で光パルスを出力する光源と、
前記光源からの光パルスをポンプ光とプローブ光とに分割するビームスプリッタと、
ポンプ光を反射して前記ビームスプリッタに再入射させる第1ミラーを有し、ポンプ光の光路長が固定された第1アームと、
第1揺動位置と第2揺動位置との間で揺動し、揺動角度に応じて前記ビームスプリッタから入射するプローブ光の反射方向を変化させるとともに、前記ビームスプリッタとの間のプローブ光の光路長を変化させるスキャニングミラーと、前記スキャニングミラーからのプローブ光を反射して前記スキャニングミラーを介して前記ビームスプリッタに再入射させるミラー部とを有し、時間的に非線形に変化する前記スキャニングミラーの速度によりプローブ光の光路長を時間的に非線形に変化させて、ポンプ光に対して遅延するプローブ光を生成する第2アームと
を有することを特徴とする請求項1に記載の分光測定装置。 The pulse pair generator
A light source that outputs an optical pulse at a constant repetition cycle,
A beam splitter that splits the light pulse from the light source into pump light and probe light,
A first arm having a first mirror that reflects the pump light and re-incidents it into the beam splitter, and a fixed optical path length of the pump light.
It swings between the first swing position and the second swing position, changes the reflection direction of the probe light incident from the beam splitter according to the swing angle, and the probe light between the beam splitter and the beam splitter. A scanning mirror that changes the optical path length of the above, and a mirror portion that reflects probe light from the scanning mirror and re-incidents it into the beam splitter via the scanning mirror, and the scanning changes non-linearly with time. The spectroscopic measurement according to claim 1, further comprising a second arm that generates a probe light that is delayed with respect to the pump light by changing the optical path length of the probe light non-linearly with respect to the speed of the mirror. apparatus.
前記パルス対生成ステップで生成されるパルス対が繰り返し照射される測定対象からプローブ光が照射されるごとに得られる観測光の強度を時間領域インターフェログラムの観測データとして取得する観測データ取得ステップと、
前記観測データ取得ステップで取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成ステップと
を有することを特徴とする分光測定方法。 A pulse pair production step that repeatedly generates a pulse pair consisting of a pump light and a probe light by changing the delay time with respect to the pump light in a non-linear manner.
The observation data acquisition step of acquiring the intensity of the observation light obtained each time the probe light is irradiated from the measurement target to which the pulse pair generated in the pulse pair production step is repeatedly irradiated as the observation data of the time region interferogram. ,
Having a spectrum spectrum generation step of matching the observation data acquired in the observation data acquisition step and obtaining a spectrum that minimizes the L1 norm or the L0 norm as a spectrum obtained by converting the time domain interferogram into the frequency domain. A characteristic spectral measurement method.
前記干渉計からの干渉光を測定対象に照射することにより得られる観測光の強度を一定の時間間隔で時間領域インターフェログラムの観測データとして順次に取得する観測データ取得部と、
前記観測データ取得部によって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成部と
を備えることを特徴とする分光測定装置。 An interferometer that generates interference light by combining the reference light and signal light by changing the phase of the signal light with respect to the reference light in a non-linear manner in time.
An observation data acquisition unit that sequentially acquires the intensity of the observation light obtained by irradiating the measurement target with the interference light from the interferometer as observation data of the time region interferogram at regular time intervals.
It is provided with a spectral spectrum generation unit that obtains a spectrum that matches the observation data acquired by the observation data acquisition unit and minimizes the L1 norm or the L0 norm as a spectrum obtained by converting the time domain interferogram into the frequency domain. A featured spectroscopic measuring device.
前記干渉光生成ステップで生成される干渉光を測定対象に照射することにより得られる観測光の強度を一定の時間間隔で時間領域インターフェログラムの観測データとして順次に取得する観測データ取得ステップと、
前記観測データ取得ステップによって取得された観測データに適合し、L1ノルムまたはL0ノルムを最小化するスペクトルを、時間領域インターフェログラムを周波数領域に変換したスペクトルとして求める分光スペクトル生成ステップと
を有することを特徴とする分光測定方法。
An interference light generation step in which the phase of the signal light with respect to the reference light is changed in a non-linear manner in time to generate interference light in which the reference light and the signal light are combined.
An observation data acquisition step in which the intensity of the observation light obtained by irradiating the measurement target with the interference light generated in the interference light generation step is sequentially acquired as observation data of the time region interferogram at regular time intervals.
Having a spectral spectrum generation step of matching the observation data acquired by the observation data acquisition step and obtaining a spectrum that minimizes the L1 norm or the L0 norm as a spectrum obtained by converting the time domain interferogram into the frequency domain. A characteristic spectral measurement method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019170884A JP2021047130A (en) | 2019-09-19 | 2019-09-19 | Spectrometry device and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019170884A JP2021047130A (en) | 2019-09-19 | 2019-09-19 | Spectrometry device and method |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2021047130A true JP2021047130A (en) | 2021-03-25 |
Family
ID=74878299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019170884A Pending JP2021047130A (en) | 2019-09-19 | 2019-09-19 | Spectrometry device and method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2021047130A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113074813A (en) * | 2021-03-30 | 2021-07-06 | 青岛青源峰达太赫兹科技有限公司 | Terahertz time-domain spectroscopy system and step length self-adaptive adjustment method thereof |
-
2019
- 2019-09-19 JP JP2019170884A patent/JP2021047130A/en active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113074813A (en) * | 2021-03-30 | 2021-07-06 | 青岛青源峰达太赫兹科技有限公司 | Terahertz time-domain spectroscopy system and step length self-adaptive adjustment method thereof |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4804820B2 (en) | Optical tomographic image display system | |
US8693004B2 (en) | Dual-etalon cavity ring-down frequency-comb spectroscopy with broad band light source | |
US20150205090A1 (en) | Methods for optical amplified imaging using a two-dimensional spectral brush | |
JP7070567B2 (en) | Spectroscopic measurement method, spectroscopic measurement device and wideband pulse light source unit | |
US20100027019A1 (en) | Optical Coherence Tomographic Imaging Method and Optical Coherence Tomographic Imaging Apparatus | |
US9163988B2 (en) | Detection systems and methods using coherent anti-stokes Raman spectroscopy | |
EP3488505B1 (en) | Acoustic resonance spectroscopy method and system | |
JP5984693B2 (en) | Optical coherence tomography apparatus and optical coherence tomography method | |
WO2017119389A1 (en) | Fourier transform-type spectroscopic device | |
CN110274880A (en) | A kind of optical spectrum detecting method and system of high-precision spatial resolution | |
WO2020075441A1 (en) | Light source for spectroscopic analysis, spectroscopic analysis device, and spectroscopic analysis method | |
JP2007278768A (en) | Microscope device | |
US20140009743A1 (en) | Apparatus and method for optical coherence tomography | |
US9128059B2 (en) | Coherent anti-stokes raman spectroscopy | |
CN111122535A (en) | Hyperspectral rapid imaging measurement system for molecular vibration mode | |
JP2021047130A (en) | Spectrometry device and method | |
JP2020159972A (en) | Broadband pulse light source device, spectral measurement device and spectral measurement method | |
KR100996638B1 (en) | Terahertz Frequency Comb Fourier Transform Spectrometer and Spectroscopic Method | |
Jin et al. | A concept of multi-mode high spectral resolution lidar using Mach-Zehnder interferometer | |
JP7272652B2 (en) | Two-dimensional spectroscopy and two-dimensional spectroscopic device | |
KR102340037B1 (en) | Apparatus and method for identifying for traces gases | |
WO2023112909A1 (en) | Time-stretch optical measurement instrument and time-stretch spectroscopy | |
Hansen | Mid-Infrared Supercontinuum based Spectroscopic OCT | |
JP2006300664A (en) | Fourier spectral device and measuring timing detection method | |
JP2019196949A (en) | Optical interference device and OCT device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20191010 |