JP6229923B2 - テラヘルツ波検出装置と方法 - Google Patents
テラヘルツ波検出装置と方法 Download PDFInfo
- Publication number
- JP6229923B2 JP6229923B2 JP2013081055A JP2013081055A JP6229923B2 JP 6229923 B2 JP6229923 B2 JP 6229923B2 JP 2013081055 A JP2013081055 A JP 2013081055A JP 2013081055 A JP2013081055 A JP 2013081055A JP 6229923 B2 JP6229923 B2 JP 6229923B2
- Authority
- JP
- Japan
- Prior art keywords
- terahertz wave
- light
- excitation light
- detection signal
- terahertz
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims description 140
- 238000000034 method Methods 0.000 title description 7
- 230000005284 excitation Effects 0.000 claims description 85
- 230000010287 polarization Effects 0.000 claims description 78
- 230000000737 periodic effect Effects 0.000 claims description 52
- 230000003287 optical effect Effects 0.000 claims description 43
- 239000013307 optical fiber Substances 0.000 claims description 28
- 239000013078 crystal Substances 0.000 claims description 16
- 239000000835 fiber Substances 0.000 claims description 15
- 238000000926 separation method Methods 0.000 claims description 13
- GQYHUHYESMUTHG-UHFFFAOYSA-N lithium niobate Chemical compound [Li+].[O-][Nb](=O)=O GQYHUHYESMUTHG-UHFFFAOYSA-N 0.000 claims description 5
- 230000003111 delayed effect Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 239000000758 substrate Substances 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 238000005286 illumination Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 15
- 230000008878 coupling Effects 0.000 description 8
- 238000010168 coupling process Methods 0.000 description 8
- 238000005859 coupling reaction Methods 0.000 description 8
- 230000003993 interaction Effects 0.000 description 8
- 230000003321 amplification Effects 0.000 description 7
- 238000003199 nucleic acid amplification method Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 229910013641 LiNbO 3 Inorganic materials 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 229910052734 helium Inorganic materials 0.000 description 1
- 239000001307 helium Substances 0.000 description 1
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000004038 photonic crystal Substances 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/02—Details
- G01J1/04—Optical or mechanical part supplementary adjustable parts
- G01J1/0407—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings
- G01J1/0425—Optical elements not provided otherwise, e.g. manifolds, windows, holograms, gratings using optical fibers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/12—Generating the spectrum; Monochromators
- G01J3/18—Generating the spectrum; Monochromators using diffraction elements, e.g. grating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/42—Absorption spectrometry; Double beam spectrometry; Flicker spectrometry; Reflection spectrometry
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/35—Non-linear optics
- G02F1/355—Non-linear optics characterised by the materials used
- G02F1/3558—Poled materials, e.g. with periodic poling; Fabrication of domain inverted structures, e.g. for quasi-phase-matching [QPM]
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/35—Non-linear optics
- G02F1/365—Non-linear optics in an optical waveguide structure
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
- G02F1/35—Non-linear optics
- G02F1/353—Frequency conversion, i.e. wherein a light beam is generated with frequency components different from those of the incident light beams
- G02F1/3544—Particular phase matching techniques
- G02F1/3548—Quasi phase matching [QPM], e.g. using a periodic domain inverted structure
-
- G—PHYSICS
- G02—OPTICS
- G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
- G02F2203/00—Function characteristic
- G02F2203/13—Function characteristic involving THZ radiation
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- Optics & Photonics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は、非線形光学効果を用いて微弱なテラヘルツ波を検出するテラヘルツ波検出装置と方法に関する。
テラヘルツ波（ｔｅｒａｈｅｒｔｚ ｗａｖｅ）とは周波数が０．１〜１０ＴＨｚ（１ＴＨｚ＝１０１２Ｈｚ）の領域、すなわち波長が０．０３〜３ｍｍのミリ波から遠赤外線領域の電磁波（ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ ｗａｖｅ）を意味する。
テラヘルツ波は、電波天文学、材料科学、生体分子分光学などの基礎学術分野から、セキュリティ、情報通信、環境、医療などの実用分野に至る幅広い分野での応用が期待されている。
しかし、テラヘルツ波は、中近赤外線、可視光、紫外線などの光（周波数１×１０１３〜１０１５Ｈｚ）と電波（周波数１０３〜１０１２Ｈｚ）の間に挟まれた周波数帯域の電磁波であり、光学と電子工学という既存の技術がそのままでは適用できない問題点がある。
非線形光学効果を用いてテラヘルツ波を検出する検出手段は、例えば特許文献１〜３に開示されている。
また本願発明に関連する論文は、例えば非特許文献１〜５に開示されている。
また本願発明に関連する論文は、例えば非特許文献１〜５に開示されている。
非特許文献１は、周期分極反転素子を用いたテラヘルツ波発生に関する論文であり、本発明とは異なる。
また、非特許文献２〜５は、バルク結晶を用いたテラヘルツ波検出に関するものであり、周期分極反転素子を用いたテラヘルツ波検出とは異なる。
特許文献１の「単色波長可変型テラヘルツ波発生／検出システム及び方法」は、非線形光学効果を有するバルク結晶を用いたテラヘルツ波検出手段である。
特許文献２の「テラヘルツ波発生素子、テラヘルツ波検出素子、及びテラヘルツ時間領域分光装置」は、ピコ秒（１０−１２秒）より短いパルス幅の超短パルスレーザー光をバルク結晶に入射してテラヘルツ波検出を行うものである。
特許文献３の「テラヘルツ光検出素子および光学設備」は、屈折率の異なる層を交互に組み合わせた結晶（フォトニック結晶）を用いるものである。
また、非線形光学効果以外のテラヘルツ波検出器として、熱検出型テラヘルツ波検出器、テラヘルツ二乗検出器、量子検出器が知られている。
熱検出型テラヘルツ波検出器は、ボロメータ、焦電効果検出器、ゴーレイセルなどが該当し、テラヘルツ波を熱エネルギーとして検出するものである。非線形光学効果を用いる本発明とはまったく別である。代表的な熱型検出器については以下のとおりである。
シリコンボロメータなど４Ｋの極低温下で動作させる検出器は、比較的検出感度が高いが、液体ヘリウムを使用する必要があって、実応用で汎用的に利用できない。一方、焦電効果検出器やゴーレイセルは、常温で動作するものの、ボロメータに比べ検出感度が２ケタ以上低く、テラヘルツ波光源の高出力化が困難であるため、利用上問題となる場合が多い。
また、これらの検出器は、基本的に応答速度がマイクロ秒〜ミリ秒と遅く、時間分解分光などの高度な計測には使用できない。
シリコンボロメータなど４Ｋの極低温下で動作させる検出器は、比較的検出感度が高いが、液体ヘリウムを使用する必要があって、実応用で汎用的に利用できない。一方、焦電効果検出器やゴーレイセルは、常温で動作するものの、ボロメータに比べ検出感度が２ケタ以上低く、テラヘルツ波光源の高出力化が困難であるため、利用上問題となる場合が多い。
また、これらの検出器は、基本的に応答速度がマイクロ秒〜ミリ秒と遅く、時間分解分光などの高度な計測には使用できない。
テラヘルツ二乗検出器は、例えばショットキ・ダイオードなどが該当する。
高速動作が可能なＧａＡｓなどの半導体を用いたショットキ・ダイオードなどは、室温で動作し、ナノ秒より短い時間のパルス計測なども可能なテラヘルツ波検出器であるが、アンテナを介してテラヘルツ波を検出する構造であることから、アンテナの性能に大きく依存する。特に、アンテナ設計は、特定の周波数領域で最適となるように行うため、１〜３ＴＨｚなど広帯域にわたり高効率にテラヘルツ波検出を実現することはできない。また、テラヘルツ波は、波長がマイクロ波などと比べ短く、数百μｍ以下であるため、波長に応じて縮小され、作製時のエラーが性能に大きく寄与してしまう。
さらには、高周波数まで応答するショットキ・ダイオードは、ウィスカーアンテナが用いられ、針状のアンテナを検出器に接触させて使用する。この仕組みでは、振動などによる機械的なショックで、接触が外れる場合があり、安定したテラヘルツ波計測には問題となる。
高速動作が可能なＧａＡｓなどの半導体を用いたショットキ・ダイオードなどは、室温で動作し、ナノ秒より短い時間のパルス計測なども可能なテラヘルツ波検出器であるが、アンテナを介してテラヘルツ波を検出する構造であることから、アンテナの性能に大きく依存する。特に、アンテナ設計は、特定の周波数領域で最適となるように行うため、１〜３ＴＨｚなど広帯域にわたり高効率にテラヘルツ波検出を実現することはできない。また、テラヘルツ波は、波長がマイクロ波などと比べ短く、数百μｍ以下であるため、波長に応じて縮小され、作製時のエラーが性能に大きく寄与してしまう。
さらには、高周波数まで応答するショットキ・ダイオードは、ウィスカーアンテナが用いられ、針状のアンテナを検出器に接触させて使用する。この仕組みでは、振動などによる機械的なショックで、接触が外れる場合があり、安定したテラヘルツ波計測には問題となる。
量子検出器は、量子ドット検出器、半導体光伝導検出器などが該当する。
量子検出器は、感度が高く応答速度も速いが、極低温下での動作温度となるため、実応用上では、汎用的に用いられる検出器ではなく、天文など究極の性能が求められる限られた応用での利用にとどまる。
量子検出器は、感度が高く応答速度も速いが、極低温下での動作温度となるため、実応用上では、汎用的に用いられる検出器ではなく、天文など究極の性能が求められる限られた応用での利用にとどまる。
Ｙ． Ｓａｓａｋｉ， Ａ． Ｙｕｒｉ， Ｋ． Ｋａｗａｓｅ， ａｎｄ Ｈ． Ｉｔｏ， "Ｔｅｒａｈｅｒｔｚ‐ｗａｖｅ ｓｕｒｆａｃｅ‐ｅｍｉｔｔｅｄ ｄｉｆｆｅｒｅｎｃｅ ｆｒｅｑｕｅｎｃｙ ｇｅｎｅｒａｔｉｏｎ ｉｎ ｓｌａｎｔ‐ｓｔｒｉｐｅ‐ｔｙｐｅ ｐｅｒｉｏｄｉｃａｌｌｙ ｐｏｌｅｄ ＬｉＮｂＯ３ ｃｒｙｓｔａｌ，" Ａｐｐｌ． Ｐｈｙｓ． Ｌｅｔｔ． ８１， ３３２３（２００２）．
Ａ． Ａ． Ｂａｂｉｎ， Ｖ． Ｎ． Ｐｅｔｒｙａｋｏｖ， ａｎｄ Ｇ． Ｉ． Ｆｒｅｉｄｍａｎ， "Ｕｓｅ ｏｆ ｓｔｉｍｕｌａｔｅｄ ｓｃａｔｔｅｒｉｎｇ ｂｙ ｐｏｌａｒｉｔｏｎｓ ｉｎ ｄｅｔｅｃｔｉｏｎ ｏｆ ｓｕｂｍｉｌｌｉｍｅｔｅｒ ｒａｄｉａｔｉｏｎ，" Ｓｏｖｉｅｔ Ｊｏｕｒｎａｌ ｏｆ Ｑｕａｎｔｕｍ Ｅｌｅｃｔｒｏｎｉｃｓ １３， ９５８‐９６０ （１９８３）．
Ｙ． Ｊ． Ｄｉｎｇ， ａｎｄ Ｗ． Ｓｈｉ， "Ｅｆｆｉｃｉｅｎｔ ＴＨｚ ｇｅｎｅｒａｔｉｏｎ ａｎｄ ｆｒｅｑｕｅｎｃｙ ｕｐｃｏｎｖｅｒｓｉｏｎ ｉｎ ＧａＰ ｃｒｙｓｔａｌｓ，" Ｓｏｌｉｄ‐Ｓｔａｔｅ Ｅｌｅｃｔｒｏｎｉｃｓ ５０， １１２８‐１１３６ （２００６）．
Ｒ． Ｇｕｏ， Ｓ． Ｏｈｎｏ， Ｈ． Ｍｉｎａｍｉｄｅ，Ｔ． Ｉｋａｒｉ， ａｎｄ Ｈ． Ｉｔｏ， "Ｈｉｇｈｌｙ ｓｅｎｓｉｔｉｖｅ ｃｏｈｅｒｅｎｔ ｄｅｔｅｃｔｉｏｎ ｏｆ ｔｅｒａｈｅｒｔｚ ｗａｖｅｓ ａｔ ｒｏｏｍ ｔｅｍｐｅｒａｔｕｒｅ ｕｓｉｎｇ ａ ｐａｒａｍｅｔｒｉｃ ｐｒｏｃｅｓｓ，" Ａｐｐｌ． Ｐｈｙｓ． Ｌｅｔｔ． ９３， ０２１１０６ （２００８）．
Ｈｉｒｏａｋｉ Ｍｉｎａｍｉｄｅ， Ｊｕｎ Ｚｈａｎｇ， Ｒｕｉｘｉａｎｇ Ｇｕｏ， Ｋａｔｓｕｈｉｋｏ Ｍｉｙａｍｏｔｏ， Ｓｅｉｇｏ Ｏｈｎｏ， ａｎｄ Ｈｉｒｏｍａｓａ Ｉｔｏ， "Ｈｉｇｈ‐ｓｅｎｓｉｔｉｖｉｔｙ ｄｅｔｅｃｔｉｏｎ ｏｆ ｔｅｒａｈｅｒｔｚ ｗａｖｅｓ ｕｓｉｎｇ ｎｏｎｌｉｎｅａｒ ｕｐ‐ｃｏｎｖｅｒｓｉｏｎ ｉｎ ａｎ ｏｒｇａｎｉｃ ４‐ｄｉｍｅｔｈｙｌａｍｉｎｏ‐Ｎ‐ｍｅｔｈｙｌ‐４‐ｓｔｉｌｂａｚｏｌｉｕｍ ｔｏｓｙｌａｔｅ ｃｒｙｓｔａｌ，" Ａｐｐｌ． Ｐｈｙｓ． Ｌｅｔｔ． ９７， １２１１０６ （２０１０）．
非線形光学効果を用いた従来のテラヘルツ波検出手段は、非線形光学結晶として例えばニオブ酸リチウム（ＬｉＮｂＯ３）結晶を用いる。励起光とテラヘルツ波が、非線形光学結晶に入射すると、テラヘルツ波周波数に応じた異なる波長の検出信号光が発生するので、発生した検出信号光の強度や波長を検出することによって、テラヘルツ波の強度や周波数を測定するものである。
図１は、非線形光学効果を用いた従来のテラヘルツ波検出装置の模式図であり、（Ａ）は従来のテラヘルツ波検出光学系、（Ｂ）は位相整合条件を示している。
図１（Ａ）に示すように、非線形光学効果を有するバルク結晶４（この例でＬｉＮｂＯ３）を用いてテラヘルツ波３を検出するには、バルク結晶４の内部で位相整合条件（図１（Ｂ）参照）という波長変換のために必要なエネルギー保存則（ω１＝ω２＋ω３）および運動量保存則（ｋ１＝ｋ２＋ｋ３）を同時に満たさなくてはならない。
ここで、ω１は励起光１、ω２は検出信号光２、ω３はテラヘルツ波３のそれぞれ角周波数（すなわちエネルギー）、ｋ１は励起光１、ｋ２は検出信号光２、ｋ３はテラヘルツ波３のそれぞれ波数（すなわち運動量）である。
図１（Ａ）に示すように、非線形光学効果を有するバルク結晶４（この例でＬｉＮｂＯ３）を用いてテラヘルツ波３を検出するには、バルク結晶４の内部で位相整合条件（図１（Ｂ）参照）という波長変換のために必要なエネルギー保存則（ω１＝ω２＋ω３）および運動量保存則（ｋ１＝ｋ２＋ｋ３）を同時に満たさなくてはならない。
ここで、ω１は励起光１、ω２は検出信号光２、ω３はテラヘルツ波３のそれぞれ角周波数（すなわちエネルギー）、ｋ１は励起光１、ｋ２は検出信号光２、ｋ３はテラヘルツ波３のそれぞれ波数（すなわち運動量）である。
（１）その結果、図１（Ｂ）から明らかなように、検出信号光２は励起光１と違う角度θで進む。ここで角度θは、検出信号光２と励起光１がなす角度である。
この角度θがあることにより、励起光１と検出信号光２がバルク結晶４の内部を進行するに従って空間的に離れるため、２つの光波（励起光１と検出信号光２）が相互作用する領域（相互作用領域）は小さい。そのため従来のテラヘルツ波検出装置は、テラヘルツ波を検出容易な検出信号光へ変換する変換効率が低いという問題点があった。
（２）また、テラヘルツ波３の周波数に応じて、検出信号光２の出力位置と角度θが変化する。そのため、テラヘルツ波３の周波数が変化した場合でも発生した検出信号光２を光検出器５に導くためには、バルク結晶４の後段に補正光学系６が必要であり、系全体が大型化する。
（３）さらに、テラヘルツ波３は、シリコンプリズムなどの入力結合素子７を介して、バルク結晶４の側面から斜めに導入される。したがって、テラヘルツ波３を相互作用領域に導入させるアライメントは、屈折率の高い入力結合素子７を介して行われるため困難を伴い、かつ入力結合素子７でのテラヘルツ波導入損失（フレネル損）が大きかった。
この角度θがあることにより、励起光１と検出信号光２がバルク結晶４の内部を進行するに従って空間的に離れるため、２つの光波（励起光１と検出信号光２）が相互作用する領域（相互作用領域）は小さい。そのため従来のテラヘルツ波検出装置は、テラヘルツ波を検出容易な検出信号光へ変換する変換効率が低いという問題点があった。
（２）また、テラヘルツ波３の周波数に応じて、検出信号光２の出力位置と角度θが変化する。そのため、テラヘルツ波３の周波数が変化した場合でも発生した検出信号光２を光検出器５に導くためには、バルク結晶４の後段に補正光学系６が必要であり、系全体が大型化する。
（３）さらに、テラヘルツ波３は、シリコンプリズムなどの入力結合素子７を介して、バルク結晶４の側面から斜めに導入される。したがって、テラヘルツ波３を相互作用領域に導入させるアライメントは、屈折率の高い入力結合素子７を介して行われるため困難を伴い、かつ入力結合素子７でのテラヘルツ波導入損失（フレネル損）が大きかった。
本発明は、上述した問題点を解決するために創案されたものである。すなわち本発明の目的は、テラヘルツ波導入損失を低減し、かつ励起光と検出信号光の相互作用領域を大きくして、微弱なテラヘルツ波を検出容易な検出信号光へ光パラメトリック増幅により高効率に変換することができるテラヘルツ波検出装置と方法を提供することにある。
本発明によれば、互いに平行な入射面及び出射面と、前記入射面及び前記出射面に対し直交する受光側面とを有する斜周期分極反転素子と、
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を備え、
前記斜周期分極反転素子は、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光を前記励起光と同一方向に平行に発生するように構成されている、ことを特徴とするテラヘルツ波検出装置が提供される。
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を備え、
前記斜周期分極反転素子は、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光を前記励起光と同一方向に平行に発生するように構成されている、ことを特徴とするテラヘルツ波検出装置が提供される。
前記疑似位相整合は、前記励起光に対する分極反転の角度と分極反転の周期とによって位相整合条件に運動量を与えるものである。
前記斜周期分極反転素子は、周期分極反転ニオブ酸リチウム結晶である、ことが好ましい。
前記光分離素子は、誘電体多層膜フィルター、又はグレイティングである。
前記斜周期分極反転素子と光分離素子の間に、前記検出信号光又は前記励起光の周波数を２倍にする第二高調波発生装置を備える、ことが好ましい。
前記光検出器は、テラヘルツ波の強度又は周波数を計測するフォトダイオード又はスペクトルアナライザである。
さらに、テラヘルツ波を試料に照射して通過したサンプル波を出力する試料照射光学系を備え、
前記テラヘルツ波は、前記サンプル波と試料を通過しないリファレンス波とからなり、前記光検出器によりテラヘルツ波の位相と振幅を計測する。
前記テラヘルツ波は、前記サンプル波と試料を通過しないリファレンス波とからなり、前記光検出器によりテラヘルツ波の位相と振幅を計測する。
さらに、テラヘルツ波を遅延させて遅延したサンプル波を出力する時間遅延光学系を備え、
前記テラヘルツ波は、前記サンプル波と遅延しないリファレンス波とからなり、前記光検出器によりテラヘルツ波の位相と振幅を計測する。
前記テラヘルツ波は、前記サンプル波と遅延しないリファレンス波とからなり、前記光検出器によりテラヘルツ波の位相と振幅を計測する。
前記斜周期分極反転素子は、
前記励起光又は検出信号光を内部に閉じ込めるバルク型、
該バルク型に基板を付加したリッジ型、或いは、
前記バルク型又はリッジ型の受光側面とこれに対向する側面とにそれぞれ一体化されその間より屈折率が小さい低屈折率部を有するスラブ型の導波路構造である。
前記励起光又は検出信号光を内部に閉じ込めるバルク型、
該バルク型に基板を付加したリッジ型、或いは、
前記バルク型又はリッジ型の受光側面とこれに対向する側面とにそれぞれ一体化されその間より屈折率が小さい低屈折率部を有するスラブ型の導波路構造である。
前記斜周期分極反転素子は、並列又は直列に配置された複数の分割素子からなり、
各分割素子の分極反転の角度と分極反転の周期は、互いに異なる周波数のテラヘルツ波に対応する。
各分割素子の分極反転の角度と分極反転の周期は、互いに異なる周波数のテラヘルツ波に対応する。
前記斜周期分極反転素子の入射面と励起光源を連結する光ファイバーからなる励起光ファイバーと、
前記斜周期分極反転素子の出射面に一端が結合された光ファイバーからなる検出信号光ファイバーと、
検出信号光ファイバーの他端に結合され前記検出信号光のみを前記光検出器に出力するファイバー周波数フィルターと、を備える。
前記斜周期分極反転素子の出射面に一端が結合された光ファイバーからなる検出信号光ファイバーと、
検出信号光ファイバーの他端に結合され前記検出信号光のみを前記光検出器に出力するファイバー周波数フィルターと、を備える。
さらに、前記斜周期分極反転素子の受光側面とテラヘルツ波光源を連結する光ファイバーからなるテラヘルツ波ファイバーを備える、ことが好ましい。
また本発明によれば、互いに平行な入射面及び出射面と、入射面及び出射面に対し直交する受光側面とを有する斜周期分極反転素子と、
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を準備し、
前記斜周期分極反転素子を、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光が前記励起光と同一方向に平行に発生するように構成する、ことを特徴とするテラヘルツ波検出方法が提供される。
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を準備し、
前記斜周期分極反転素子を、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光が前記励起光と同一方向に平行に発生するように構成する、ことを特徴とするテラヘルツ波検出方法が提供される。
上記本発明の装置と方法によれば、斜周期分極反転素子が、受光側面から垂直に入射するテラヘルツ波と入射面から垂直に入射する励起光との疑似位相整合により、検出信号光を励起光と同一方向に平行に発生するように構成されているので、入力結合素子を省略してテラヘルツ波導入損失を低減することができ、かつ検出信号光を励起光と同一方向に発生かつ伝搬させることができる。
斜周期分極反転素子の内部（光導波路内）で励起光と同一方向に伝搬する検出信号光は、パラメトリック位相整合を満たしているので励起光により光パラメトリック増幅される。この光パラメトリック増幅は、テラヘルツ波の入射位置から斜周期分極反転素子の出射面まで継続するので、励起光と検出信号光の相互作用領域を大きくでき、テラヘルツ波を検出信号光に変換する変換効率を飛躍的に高めることができる。
従って、本発明によれば、テラヘルツ波導入損失を低減し、かつ励起光と検出信号光の相互作用領域を大きくして、微弱なテラヘルツ波を検出容易な検出信号光へ光パラメトリック増幅により高効率に変換することができる。
以下、本発明の好ましい実施形態を図面を参照して説明する。なお、各図において、共通する部分には同一の符号を付し重複した説明を省略する。
図２は、本発明によるテラヘルツ波検出装置１０の第１実施形態を示す図であり、（Ａ）は本発明のテラヘルツ波検出光学系、（Ｂ）は本発明の位相整合条件を示している。
図２において、本発明のテラヘルツ波検出装置１０は、斜周期分極反転素子１２、励起光源１４、及び光検出器１６を備える。
斜周期分極反転素子１２は、互いに平行な入射面１２ａ及び出射面１２ｂと、入射面１２ａ及び出射面１２ｂに対し直交する受光側面１２ｃとを有する。
励起光源１４は、斜周期分極反転素子１２の入射面１２ａに垂直に励起光１を入射する。
光検出器１６は、斜周期分極反転素子１２の出射面１２ｂから垂直に出射した検出信号光Ａからテラヘルツ波３を検出する。
斜周期分極反転素子１２は、互いに平行な入射面１２ａ及び出射面１２ｂと、入射面１２ａ及び出射面１２ｂに対し直交する受光側面１２ｃとを有する。
励起光源１４は、斜周期分極反転素子１２の入射面１２ａに垂直に励起光１を入射する。
光検出器１６は、斜周期分極反転素子１２の出射面１２ｂから垂直に出射した検出信号光Ａからテラヘルツ波３を検出する。
励起光１は、波長１〜１０μｍの赤外レーザ光、又は可視レーザ光であることが好ましい。以下、波長１〜１０μｍの赤外レーザ光、又は可視レーザ光を「光波」と呼ぶ。
斜周期分極反転素子１２は、この例では、周期分極反転ニオブ酸リチウム結晶である。
斜周期分極反転素子１２の受光側面１２ｃは、この受光側面１２ｃから垂直に入射するテラヘルツ波３と励起光１との疑似位相整合により、検出信号光Ａを励起光１と同一方向に平行に発生するように構成されている。
疑似位相整合は、励起光１に対する分極反転の角度α（°）と分極反転の周期Λ（μｍ）とによって位相整合条件に運動量ｋΛを与えるものである。
疑似位相整合は、励起光１に対する分極反転の角度α（°）と分極反転の周期Λ（μｍ）とによって位相整合条件に運動量ｋΛを与えるものである。
疑似位相整合は、以下の式で示すことができる。
ｋΛｓｉｎα＝ｋ１−ｋ２・・・（１）
ｋΛｃｏｓα＝ｋ３ ・・・（２）
ｋΛ＝２π／Λ ・・・（３）
ｋｊ＝ωｊｎｊ／ｃ ・・・（４）
ここで、ｋΛは格子波数（Ｇｒａｔｉｎｇ ｗａｖｅ ｎｕｍｂｅｒ）、ｊは１，２，３の整数、ωは角周波数、ｎは角周波数に対する素子の屈折率、ｃは光速である。
ｋΛｓｉｎα＝ｋ１−ｋ２・・・（１）
ｋΛｃｏｓα＝ｋ３ ・・・（２）
ｋΛ＝２π／Λ ・・・（３）
ｋｊ＝ωｊｎｊ／ｃ ・・・（４）
ここで、ｋΛは格子波数（Ｇｒａｔｉｎｇ ｗａｖｅ ｎｕｍｂｅｒ）、ｊは１，２，３の整数、ωは角周波数、ｎは角周波数に対する素子の屈折率、ｃは光速である。
また本発明によるテラヘルツ波検出方法は、図２に示すように、互いに平行な入射面１２ａ及び出射面１２ｂと、入射面１２ａ及び出射面１２ｂに対し直交する受光側面１２ｃとを有する斜周期分極反転素子１２と、入射面１２ａに垂直に励起光１を入射する励起光源１４と、出射面１２ｂから垂直に出射した検出信号光Ａからテラヘルツ波３を検出する光検出器１６と、を準備する。
また、斜周期分極反転素子１２を、受光側面１２ｃから垂直に入射するテラヘルツ波３と励起光１との疑似位相整合により、検出信号光Ａが励起光１と同一方向に平行に発生するように構成する。
また、斜周期分極反転素子１２を、受光側面１２ｃから垂直に入射するテラヘルツ波３と励起光１との疑似位相整合により、検出信号光Ａが励起光１と同一方向に平行に発生するように構成する。
励起光１が光波（波長１〜１０μｍの赤外レーザ光、又は可視レーザ光）である場合、検出信号光Ａも光波となる。
図２において、本発明のテラヘルツ波検出装置１０は、さらに、斜周期分極反転素子１２と光検出器１６の間に、検出信号光Ａと励起光１の２光波から検出信号光Ａのみを分離する光分離素子１８を備える。
光分離素子１８は、この例では励起光１を透過させ、検出信号光Ａを反射する半透過ミラーである。
光分離素子１８は、この例では励起光１を透過させ、検出信号光Ａを反射する半透過ミラーである。
図２において、斜周期分極反転素子１２は、分極反転した部分を斜線で示している。
斜周期分極反転素子１２は、分極反転の角度α（°）と周期Λ（μｍ）によって、位相整合条件に運動量ｋΛを与えることができる（図２（Ｂ）の位相整合条件を参照）。この疑似位相整合により、検出信号光Ａは励起光１とコリニア（平行）に発生可能になり、大きな相互作用領域を実現できる。その結果、光パラメトリック増幅が、テラヘルツ波３の入射位置から斜周期分極反転素子１２の出射面１２ｂまで継続するので、テラヘルツ波３から検出信号光Ａへの高い変換効率が得られる。
ここで「光パラメトリック増幅」とは、斜周期分極反転素子１２の内部（光導波路内）で励起光１と同一方向に検出信号光Ａが伝搬するので、励起光１のエネルギーにより検出信号光Ａが増幅される現象を意味する。
また、テラヘルツ波３を斜周期分極反転素子１２の受光側面１２ｃから垂直に入射できるので、入力結合素子７（図１参照）が不要になり、入力結合素子７によるテラヘルツ波導入損失（フレネル損）をなくすことができる。
斜周期分極反転素子１２は、分極反転の角度α（°）と周期Λ（μｍ）によって、位相整合条件に運動量ｋΛを与えることができる（図２（Ｂ）の位相整合条件を参照）。この疑似位相整合により、検出信号光Ａは励起光１とコリニア（平行）に発生可能になり、大きな相互作用領域を実現できる。その結果、光パラメトリック増幅が、テラヘルツ波３の入射位置から斜周期分極反転素子１２の出射面１２ｂまで継続するので、テラヘルツ波３から検出信号光Ａへの高い変換効率が得られる。
ここで「光パラメトリック増幅」とは、斜周期分極反転素子１２の内部（光導波路内）で励起光１と同一方向に検出信号光Ａが伝搬するので、励起光１のエネルギーにより検出信号光Ａが増幅される現象を意味する。
また、テラヘルツ波３を斜周期分極反転素子１２の受光側面１２ｃから垂直に入射できるので、入力結合素子７（図１参照）が不要になり、入力結合素子７によるテラヘルツ波導入損失（フレネル損）をなくすことができる。
図３は、本発明によるテラヘルツ波検出装置１０の第２実施形態を示す図である。
この図において、本発明のテラヘルツ波検出装置１０は、さらに、テラヘルツ波３を試料Ｓに照射して通過したサンプル波３ａを出力する試料照射光学系２０を備える。
試料照射光学系２０は、この例では、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３の一部を試料Ｓに向けて反射する第１ミラー２０ａと、試料Ｓを透過したテラヘルツ波３（サンプル波３ａと呼ぶ）を斜周期分極反転素子１２の受光側面１２ｃに向けて反射する第２ミラー２０ｂとを有する。
この構成により、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３は、サンプル波３ａと試料Ｓを通過しないリファレンス波３ｂとからなる。
また、光検出器１６は、テラヘルツ波３（サンプル波３ａとリファレンス波３ｂ）の強度又は周波数を計測するフォトダイオード又はスペクトルアナライザである。
図３のテラヘルツ波検出装置１０の構成により、テラヘルツ波３を分光して、テラヘルツ波３の位相と振幅を計測することができる。
またこの構成により、テラヘルツ波検出装置１０をテラヘルツ波による光変調器として用いることができる。
この図において、本発明のテラヘルツ波検出装置１０は、さらに、テラヘルツ波３を試料Ｓに照射して通過したサンプル波３ａを出力する試料照射光学系２０を備える。
試料照射光学系２０は、この例では、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３の一部を試料Ｓに向けて反射する第１ミラー２０ａと、試料Ｓを透過したテラヘルツ波３（サンプル波３ａと呼ぶ）を斜周期分極反転素子１２の受光側面１２ｃに向けて反射する第２ミラー２０ｂとを有する。
この構成により、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３は、サンプル波３ａと試料Ｓを通過しないリファレンス波３ｂとからなる。
また、光検出器１６は、テラヘルツ波３（サンプル波３ａとリファレンス波３ｂ）の強度又は周波数を計測するフォトダイオード又はスペクトルアナライザである。
図３のテラヘルツ波検出装置１０の構成により、テラヘルツ波３を分光して、テラヘルツ波３の位相と振幅を計測することができる。
またこの構成により、テラヘルツ波検出装置１０をテラヘルツ波による光変調器として用いることができる。
図４は、本発明によるテラヘルツ波検出装置１０の第３実施形態を示す図である。
この図において、本発明のテラヘルツ波検出装置１０は、さらに、テラヘルツ波３を遅延させて遅延したサンプル波３ａを出力する時間遅延光学系２２を備える。
時間遅延光学系２２は、この例では、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３の一部を受光側面１２ｃ以外に向けて反射する第１ミラー２２ａと、第１ミラー２２ａで反射したテラヘルツ波３（サンプル波３ａと呼ぶ）を斜周期分極反転素子１２の受光側面１２ｃに向けて反射する第２ミラー２２ｂとを有する。
第２ミラー２２ｂは、第１ミラー２２ａに対して相対的に移動可能であり、その移動によりサンプル波３ａを受光側面１２ｃに直接入射するテラヘルツ波３（リファレンス波３ｂと呼ぶ）に対して遅延させるようになっている。
図４のテラヘルツ波検出装置１０の構成により、時間遅延光学系２２を用いることで、テラヘルツ波３の位相と振幅を計測することができる。
またこの構成によっても、テラヘルツ波検出装置１０をテラヘルツ波による光変調器として用いることができる。
この図において、本発明のテラヘルツ波検出装置１０は、さらに、テラヘルツ波３を遅延させて遅延したサンプル波３ａを出力する時間遅延光学系２２を備える。
時間遅延光学系２２は、この例では、斜周期分極反転素子１２の受光側面１２ｃに垂直に入射するテラヘルツ波３の一部を受光側面１２ｃ以外に向けて反射する第１ミラー２２ａと、第１ミラー２２ａで反射したテラヘルツ波３（サンプル波３ａと呼ぶ）を斜周期分極反転素子１２の受光側面１２ｃに向けて反射する第２ミラー２２ｂとを有する。
第２ミラー２２ｂは、第１ミラー２２ａに対して相対的に移動可能であり、その移動によりサンプル波３ａを受光側面１２ｃに直接入射するテラヘルツ波３（リファレンス波３ｂと呼ぶ）に対して遅延させるようになっている。
図４のテラヘルツ波検出装置１０の構成により、時間遅延光学系２２を用いることで、テラヘルツ波３の位相と振幅を計測することができる。
またこの構成によっても、テラヘルツ波検出装置１０をテラヘルツ波による光変調器として用いることができる。
図５は、本発明によるテラヘルツ波検出装置１０の第４実施形態を示す図である。
この図において、斜周期分極反転素子１２は、バルク型（Ａ）、リッジ型（Ｂ）、又はスラブ型（Ｃ）の導波路構造である。
この図において、斜周期分極反転素子１２は、バルク型（Ａ）、リッジ型（Ｂ）、又はスラブ型（Ｃ）の導波路構造である。
バルク型（Ａ）の導波路構造は、励起光１又は検出信号光Ａを全反射させてその内部に閉じ込めるようになっている。
リッジ型（Ｂ）の導波路構造は、バルク型（Ａ）に基板１３ａを付加したものであり、バルク型（Ａ）と同様に励起光１又は検出信号光Ａを全反射させてその内部に閉じ込めるようになっている。
スラブ型（Ｃ）の導波路構造は、バルク型（Ａ）又はリッジ型（Ｂ）の受光側面１２ｃとこれに対向する側面とにそれぞれ一体化されその間の屈折率ｎａより屈折率ｎｂが小さい低屈折率部１３ｂを有する。
リッジ型（Ｂ）の導波路構造は、バルク型（Ａ）に基板１３ａを付加したものであり、バルク型（Ａ）と同様に励起光１又は検出信号光Ａを全反射させてその内部に閉じ込めるようになっている。
スラブ型（Ｃ）の導波路構造は、バルク型（Ａ）又はリッジ型（Ｂ）の受光側面１２ｃとこれに対向する側面とにそれぞれ一体化されその間の屈折率ｎａより屈折率ｎｂが小さい低屈折率部１３ｂを有する。
リッジ型（Ｂ）又はスラブ型（Ｃ）のような導波路構造の場合、光波（励起光１と検出信号光Ａ）は分斜周期分極反転素子１２の内部に閉じ込められる。この構成により、狭い領域に閉じ込められた光波は、パラメトリック位相整合を満たしているので励起光１により光パラメトリック増幅される非線形光学効果を誘起し、高効率に検出信号光Ａを発生させることができる。
図６は、本発明によるテラヘルツ波検出装置１０の第５実施形態を示す図であり、（Ａ）は斜周期分極反転素子１２の構成図、（Ｂ）はテラヘルツ波３の周波数ｆと分極反転角度α及び周期Λとの関係図である。
図６（Ａ）において、斜周期分極反転素子１２は、並列に配置された複数の分割素子１２−１、１２−２、１２−３からなる。
各分割素子１２−１、１２−２、１２−３の分極反転の角度αと分極反転の周期Λは、互いに異なる周波数のテラヘルツ波３に対応している。
各分割素子１２−１、１２−２、１２−３の分極反転の角度αと分極反転の周期Λは、互いに異なる周波数のテラヘルツ波３に対応している。
図６（Ｂ）は、０．５ＴＨｚ〜３ＴＨｚのテラヘルツ波３に対して、位相整合を得られる分極反転角度α（°）と周期Λ（μｍ）を計算した結果である。入力テラヘルツ波周波数がｆ１，ｆ２，ｆ３のとき、それぞれ分極反転の角度（α１，α２，α３）と周期（Λ１，Λ２，Λ３）が図６（Ｂ）からわかる。従って、図６（Ａ）に示すように、それぞれの周波数（ｆ１，ｆ２，ｆ３）に対応する３つの分割素子１２−１、１２−２、１２−３を重ねた斜周期分極反転素子１２を作成することによって、テラヘルツ波３の周波数ｆ１，ｆ２，ｆ３を検出できる。
また、分割素子は３つに限定されず、重ねる分割素子をさらに増やすことで、広帯域のテラヘルツ波３を検出することができる。
また、分割素子は３つに限定されず、重ねる分割素子をさらに増やすことで、広帯域のテラヘルツ波３を検出することができる。
図７は、本発明によるテラヘルツ波検出装置１０の第６実施形態を示す図であり、（Ａ）は斜周期分極反転素子１２の構成図、（Ｂ）はテラヘルツ波３の周波数ｆと分極反転角度α及び周期Λとの関係図である。
図７（Ａ）において、斜周期分極反転素子１２は、直列に配置された複数の分割素子１２−１、１２−２、１２−３からなる。
各分割素子１２−１、１２−２、１２−３の分極反転角度α（°）と周期Λ（μｍ）は、互いに異なる周波数のテラヘルツ波３に対応している。
その他の構成は、図６と同様である。
この構成により、分割素子１２−１、１２−２、１２−３を直列に配置した場合においても、図６と同様に広帯域のテラヘルツ波３を検出することができる。
各分割素子１２−１、１２−２、１２−３の分極反転角度α（°）と周期Λ（μｍ）は、互いに異なる周波数のテラヘルツ波３に対応している。
その他の構成は、図６と同様である。
この構成により、分割素子１２−１、１２−２、１２−３を直列に配置した場合においても、図６と同様に広帯域のテラヘルツ波３を検出することができる。
図８は、本発明によるテラヘルツ波検出装置１０の第７実施形態を示す図である。
この図において、光分離素子１８は、（Ａ）誘電体多層膜フィルター１８Ａ又は（Ｂ）グレイティング１８Ｂである。
この図において、光分離素子１８は、（Ａ）誘電体多層膜フィルター１８Ａ又は（Ｂ）グレイティング１８Ｂである。
本発明では、励起光１と検出信号光Ａがコリニアに伝搬するため、励起光１と検出信号光Ａの２光波を空間的、あるいは周波数的に分離したのち、光検出器１６へ導入することが好ましい。
図８に示すように、励起光１と検出信号光Ａの２光波を分離するための光分離素子１８として斜周期分極反転素子１２の後段に、誘電体多層膜フィルター１８Ａ又はグレイティング１８Ｂを設けることが好ましい。分離された検出信号光Ａは光検出器１６で検出される。
図８に示すように、励起光１と検出信号光Ａの２光波を分離するための光分離素子１８として斜周期分極反転素子１２の後段に、誘電体多層膜フィルター１８Ａ又はグレイティング１８Ｂを設けることが好ましい。分離された検出信号光Ａは光検出器１６で検出される。
また、図８（Ｃ）において符号１９は第二高調波発生装置である。
第二高調波発生装置１９は、斜周期分極反転素子１２と光分離素子１８（図示せず）の間に設けられ、検出信号光Ａ又は励起光１の周波数を２倍にすることが好ましい。
この構成により、検出信号光Ａ又は励起光１の周波数を倍加して、光検出器１６で容易に検出しやすい波長（例えば可視光）に変換することができる。
第二高調波発生装置１９は、斜周期分極反転素子１２と光分離素子１８（図示せず）の間に設けられ、検出信号光Ａ又は励起光１の周波数を２倍にすることが好ましい。
この構成により、検出信号光Ａ又は励起光１の周波数を倍加して、光検出器１６で容易に検出しやすい波長（例えば可視光）に変換することができる。
図９は、本発明によるテラヘルツ波検出装置１０の第８実施形態を示す図である。
この図において、本発明のテラヘルツ波検出装置１０は、さらに励起光ファイバー２４、テラヘルツ波ファイバー２６、検出信号光ファイバー２８、及びファイバー周波数フィルター３０を備える。
この図において、本発明のテラヘルツ波検出装置１０は、さらに励起光ファイバー２４、テラヘルツ波ファイバー２６、検出信号光ファイバー２８、及びファイバー周波数フィルター３０を備える。
励起光ファイバー２４は、斜周期分極反転素子１２の入射面１２ａと励起光源１４を連結する光ファイバーからなる。
テラヘルツ波ファイバー２６は、斜周期分極反転素子１２の受光側面１２ｃとテラヘルツ波光源を連結する光ファイバーからなる。
なお、テラヘルツ波ファイバー２６は、必須ではなく、これを省略してもよい。
検出信号光ファイバー２８は、斜周期分極反転素子１２の出射面１２ｂに一端が結合された光ファイバーからなる。
ファイバー周波数フィルター３０は、検出信号光ファイバー２８の他端に結合され検出信号光Ａのみを光検出器１６に出力する。
テラヘルツ波ファイバー２６は、斜周期分極反転素子１２の受光側面１２ｃとテラヘルツ波光源を連結する光ファイバーからなる。
なお、テラヘルツ波ファイバー２６は、必須ではなく、これを省略してもよい。
検出信号光ファイバー２８は、斜周期分極反転素子１２の出射面１２ｂに一端が結合された光ファイバーからなる。
ファイバー周波数フィルター３０は、検出信号光ファイバー２８の他端に結合され検出信号光Ａのみを光検出器１６に出力する。
図９において、励起光１およびテラヘルツ波３は、それぞれ励起光ファイバー２４、テラヘルツ波ファイバー２６を介して斜周期分極反転素子１２に導入される。斜周期分極反転素子１２から出力された２光波（励起光１と検出信号光Ａ）は、検出信号光ファイバー２８へ結合され、ファイバー周波数フィルター３０へ導入される。ファイバー周波数フィルター３０で検出信号光Ａのみを取り出し、光検出器１６で検出信号光Ａを検出する。この時、ファイバー増幅器３２によって検出信号光Ａを増幅したのち、光検出器１６で測定することによって、微弱な検出信号光Ａを検出でき、結果として、テラヘルツ波３の検出感度の向上が可能である。
また、図９の構成において、励起光ファイバー２４と検出信号光ファイバー２８は光波の損失が少ないため、伝送距離を長距離（例えば数十ｍ〜数ｋｍ）にできるので、テラヘルツ波３を斜周期分極反転素子１２で増幅又は変調することによりテラヘルツ波３の位相と振幅を遠隔から計測することができる。
この構成により、本発明によるテラヘルツ波検出装置１０を光ファイバーネットワークに適用することができる。またその際に、２光波（励起光１と検出信号光Ａ）を伝送波、テラヘルツ波３を変調波として用いることができる。
この構成により、本発明によるテラヘルツ波検出装置１０を光ファイバーネットワークに適用することができる。またその際に、２光波（励起光１と検出信号光Ａ）を伝送波、テラヘルツ波３を変調波として用いることができる。
上述した本発明の装置と方法によれば、斜周期分極反転素子１２が、受光側面１２ｃから垂直に入射するテラヘルツ波３と入射面１２ａから垂直に入射する励起光１との疑似位相整合により、検出信号光Ａを励起光１と同一方向に平行に発生するように構成されているので、入力結合素子７を省略してテラヘルツ波導入損失を低減することができ、かつ検出信号光Ａを励起光１と同一方向に発生かつ伝搬させることができる。
斜周期分極反転素子１２の内部（光導波路内）で励起光１と同一方向に伝搬する検出信号光Ａは、パラメトリック位相整合を満たしているので励起光１により光パラメトリック増幅される。この光パラメトリック増幅は、テラヘルツ波３の入射位置から斜周期分極反転素子１２の出射面１２ｂまで継続するので、励起光１と検出信号光Ａの相互作用領域を大きくでき、テラヘルツ波３を検出信号光Ａに変換する変換効率を飛躍的に高めることができる。
従って、本発明によれば、テラヘルツ波導入損失を低減し、かつ励起光１と検出信号光Ａの相互作用領域を大きくして、光パラメトリック増幅により微弱なテラヘルツ波３を検出容易な検出信号光Ａへ高効率に変換することができる。
また、上述した本発明の実施形態によれば、以下の付随した効果も得ることもできる。
図３に示した構成により、テラヘルツ波３を分光して、テラヘルツ波３の強度や周波数を計測することができる。
図４に示した構成により、テラヘルツ波３の位相計測が可能である。
図６、図７に示した構成により、異なる周波数のテラヘルツ波３を検出することができる。
図９に示した構成により、光ファイバー技術と融合することで小型化および安定化できる。また光導波路技術を用いて、より高効率化できる。
図３に示した構成により、テラヘルツ波３を分光して、テラヘルツ波３の強度や周波数を計測することができる。
図４に示した構成により、テラヘルツ波３の位相計測が可能である。
図６、図７に示した構成により、異なる周波数のテラヘルツ波３を検出することができる。
図９に示した構成により、光ファイバー技術と融合することで小型化および安定化できる。また光導波路技術を用いて、より高効率化できる。
なお、本発明は上述した実施例及び実施形態に限定されず、本発明の要旨を逸脱しない範囲で種々変更できることは勿論である。
α（α１，α２，α３） 分極反転の角度
Λ（Λ１，Λ２，Λ３） 分極反転の周期
ｋΛ 運動量
Ａ 検出信号光
Ｓ 試料
ｆ（ｆ１，ｆ２，ｆ３） 周波数
１ 励起光
２ 検出信号光
３ テラヘルツ波
３ａ サンプル波
３ｂ リファレンス波
４ バルク結晶
５ 光検出器
６ 補正光学系
７ 入力結合素子
１０ テラヘルツ波検出装置
１２ 斜周期分極反転素子（周期分極反転ニオブ酸リチウム結晶）
１２−１、１２−２、１２−３ 分割素子
１２ａ 入射面
１２ｂ 出射面
１２ｃ 受光側面
１３ａ 基板
１３ｂ 低屈折率部
１４ 励起光源
１６ 光検出器（フォトダイオード、スペクトルアナライザ）
１８ 光分離素子
１８Ａ 誘電体多層膜フィルター
１８Ｂ グレイティング
１９ 第二高調波発生装置
２０ 試料照射光学系
２０ａ 第１ミラー
２０ｂ 第２ミラー
２２ 時間遅延光学系
２２ａ 第１ミラー
２２ｂ 第２ミラー
２４ 励起光ファイバー
２６ テラヘルツ波ファイバー
２８ 検出信号光ファイバー
３０ ファイバー周波数フィルター
Λ（Λ１，Λ２，Λ３） 分極反転の周期
ｋΛ 運動量
Ａ 検出信号光
Ｓ 試料
ｆ（ｆ１，ｆ２，ｆ３） 周波数
１ 励起光
２ 検出信号光
３ テラヘルツ波
３ａ サンプル波
３ｂ リファレンス波
４ バルク結晶
５ 光検出器
６ 補正光学系
７ 入力結合素子
１０ テラヘルツ波検出装置
１２ 斜周期分極反転素子（周期分極反転ニオブ酸リチウム結晶）
１２−１、１２−２、１２−３ 分割素子
１２ａ 入射面
１２ｂ 出射面
１２ｃ 受光側面
１３ａ 基板
１３ｂ 低屈折率部
１４ 励起光源
１６ 光検出器（フォトダイオード、スペクトルアナライザ）
１８ 光分離素子
１８Ａ 誘電体多層膜フィルター
１８Ｂ グレイティング
１９ 第二高調波発生装置
２０ 試料照射光学系
２０ａ 第１ミラー
２０ｂ 第２ミラー
２２ 時間遅延光学系
２２ａ 第１ミラー
２２ｂ 第２ミラー
２４ 励起光ファイバー
２６ テラヘルツ波ファイバー
２８ 検出信号光ファイバー
３０ ファイバー周波数フィルター
Claims (13)
- 互いに平行な入射面及び出射面と、前記入射面及び前記出射面に対し直交する受光側面とを有する斜周期分極反転素子と、
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を備え、
前記斜周期分極反転素子は、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光を前記励起光と同一方向に平行に発生するように構成されている、ことを特徴とするテラヘルツ波検出装置。 - 前記疑似位相整合は、前記励起光に対する分極反転の角度と分極反転の周期とによって位相整合条件に運動量を与えるものである、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。
- 前記斜周期分極反転素子は、周期分極反転ニオブ酸リチウム結晶である、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。
- 前記光分離素子は、誘電体多層膜フィルター、又はグレイティングである、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。
- 前記斜周期分極反転素子と前記光分離素子の間に、前記検出信号光又は前記励起光の周波数を２倍にする第二高調波発生装置を備える、ことを特徴とする請求項４に記載のテラヘルツ波検出装置。
- 前記光検出器は、前記テラヘルツ波の強度又は周波数を計測するフォトダイオード又はスペクトルアナライザである、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。
- さらに、前記テラヘルツ波を試料に照射して通過したサンプル波を出力する試料照射光学系を備え、
前記テラヘルツ波は、前記サンプル波と試料を通過しないリファレンス波とからなり、前記光検出器により前記テラヘルツ波の位相と振幅を計測する、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。 - さらに、前記テラヘルツ波を遅延させて遅延したサンプル波を出力する時間遅延光学系を備え、
前記テラヘルツ波は、前記サンプル波と遅延しないリファレンス波とからなり、前記光検出器により前記テラヘルツ波の位相と振幅を計測する、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。 - 前記斜周期分極反転素子は、
前記励起光又は前記検出信号光を内部に閉じ込めるバルク型、
該バルク型に基板を付加したリッジ型、或いは、
前記バルク型又は前記リッジ型の前記受光側面とこれに対向する側面と一体化されその間より屈折率が小さい低屈折率部を有するスラブ型の導波路構造である、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。 - 前記斜周期分極反転素子は、並列又は直列に配置された複数の分割素子からなり、
各分割素子の分極反転の角度と分極反転の周期は、互いに異なる周波数の前記テラヘルツ波に対応する、ことを特徴とする請求項２に記載のテラヘルツ波検出装置。 - 前記斜周期分極反転素子の前記入射面と励起光源を連結する光ファイバーからなる励起光ファイバーと、
前記斜周期分極反転素子の前記出射面に一端が結合された光ファイバーからなる検出信号光ファイバーと、
検出信号光ファイバーの他端に結合され前記検出信号光のみを前記光検出器に出力するファイバー周波数フィルターと、を備える、ことを特徴とする請求項１に記載のテラヘルツ波検出装置。 - さらに、前記斜周期分極反転素子の前記受光側面とテラヘルツ波光源を連結する光ファイバーからなるテラヘルツ波ファイバーを備える、ことを特徴とする請求項１１に記載のテラヘルツ波検出装置。
- 互いに平行な入射面及び出射面と、入射面及び出射面に対し直交する受光側面とを有する斜周期分極反転素子と、
前記入射面に垂直に励起光を入射する励起光源と、
前記出射面から垂直に出射した検出信号光からテラヘルツ波を検出する光検出器と、
前記斜周期分極反転素子と前記光検出器の間に、前記検出信号光と前記励起光の２光波から前記検出信号光のみを分離する光分離素子と、を準備し、
前記斜周期分極反転素子を、前記受光側面から垂直に入射する前記テラヘルツ波と前記励起光との疑似位相整合により、前記検出信号光が前記励起光と同一方向に平行に発生するように構成する、ことを特徴とするテラヘルツ波検出方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013081055A JP6229923B2 (ja) | 2013-04-09 | 2013-04-09 | テラヘルツ波検出装置と方法 |
| US14/249,046 US9040918B2 (en) | 2013-04-09 | 2014-04-09 | Terahertz wave detection device and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013081055A JP6229923B2 (ja) | 2013-04-09 | 2013-04-09 | テラヘルツ波検出装置と方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2014203025A JP2014203025A (ja) | 2014-10-27 |
| JP6229923B2 true JP6229923B2 (ja) | 2017-11-15 |
Family
ID=51653812
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013081055A Active JP6229923B2 (ja) | 2013-04-09 | 2013-04-09 | テラヘルツ波検出装置と方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US9040918B2 (ja) |
| JP (1) | JP6229923B2 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9557267B2 (en) * | 2014-10-22 | 2017-01-31 | The Boeing Company | Terahertz imaging via simultaneous surface and sub-surface evaluation via non-linear optical response |
| WO2016147253A1 (ja) * | 2015-03-13 | 2016-09-22 | 株式会社日立ハイテクノロジーズ | テラヘルツ波測定装置 |
| JP6439575B2 (ja) * | 2015-05-08 | 2018-12-19 | 株式会社島津製作所 | 波長変換素子及びレーザ装置 |
| JP6810954B2 (ja) | 2016-09-30 | 2021-01-13 | 国立研究開発法人理化学研究所 | テラヘルツ波生成装置、光パラメトリック増幅器、テラヘルツ波検出器、および非線形光学素子 |
| JP2019144435A (ja) * | 2018-02-21 | 2019-08-29 | 沖電気工業株式会社 | テラヘルツ波検出素子 |
| JP7232498B2 (ja) * | 2018-08-31 | 2023-03-03 | 国立研究開発法人理化学研究所 | テラヘルツ波を用いた検査装置と検査方法 |
| CN110223898B (zh) * | 2019-06-18 | 2021-07-30 | 天津中安信业集团有限公司 | 一种多用途太赫兹波电子发生器 |
| CN110380326B (zh) * | 2019-07-29 | 2020-10-23 | 武汉电信器件有限公司 | 一种光信号输出装置及方法、存储介质 |
| WO2024050518A1 (en) * | 2022-09-01 | 2024-03-07 | California Institute Of Technology | Optical parametric amplification protocols for quantum nondemolition measurement |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8040929B2 (en) * | 2004-03-25 | 2011-10-18 | Imra America, Inc. | Optical parametric amplification, optical parametric generation, and optical pumping in optical fibers systems |
| GB0416673D0 (en) * | 2004-07-27 | 2004-08-25 | Univ St Andrews | Parametric generation with lateral beam coupling |
| DE102008021791A1 (de) * | 2008-04-30 | 2009-11-26 | ARIZONA BOARD OF REGENTS, on behalf of THE UNIVERSITY OF ARIZONA, Tucson | Laserbasierte Quelle für Terahertz- und Millimeterwellen |
| JP5618119B2 (ja) | 2009-03-11 | 2014-11-05 | 国立大学法人香川大学 | テラヘルツ光検出素子および光学設備 |
| JP5240858B2 (ja) | 2009-09-03 | 2013-07-17 | 独立行政法人理化学研究所 | 単色波長可変型テラヘルツ波発生／検出システム及び方法 |
| GB201008073D0 (en) * | 2010-05-14 | 2010-06-30 | Univ St Andrews | Parametric generation with optimised lateral beam output coupling |
| JP5967867B2 (ja) * | 2010-06-03 | 2016-08-10 | キヤノン株式会社 | テラヘルツ波発生素子、テラヘルツ波検出素子、及びテラヘルツ時間領域分光装置 |
-
2013
- 2013-04-09 JP JP2013081055A patent/JP6229923B2/ja active Active
-
2014
- 2014-04-09 US US14/249,046 patent/US9040918B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| US9040918B2 (en) | 2015-05-26 |
| US20140299773A1 (en) | 2014-10-09 |
| JP2014203025A (ja) | 2014-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6229923B2 (ja) | テラヘルツ波検出装置と方法 | |
| US10331010B2 (en) | Terahertz-wave generating element terahertz-wave detecting element and terahertz time-domain spectroscopy device | |
| US9091869B2 (en) | Terahertz-wave element, terahertz-wave detecting device, and terahertz time-domain spectroscopy system | |
| US8446580B2 (en) | System for generating Raman vibrational analysis signals | |
| Amiot et al. | Cavity enhanced absorption spectroscopy in the mid-infrared using a supercontinuum source | |
| JP5240858B2 (ja) | 単色波長可変型テラヘルツ波発生／検出システム及び方法 | |
| US10725359B2 (en) | Terahertz wave generating device, optical parametric amplifier, terahertz wave detector, and nonlinear optical element | |
| Kuyken et al. | Nonlinear optical interactions in silicon waveguides | |
| EP2577393A1 (en) | Terahertz wave generation element, terahertz wave detection element, and terahertz time domain spectroscope device | |
| JP2015222414A (ja) | テラヘルツ波発生装置、及びこれを用いた測定装置 | |
| Zentgraf et al. | Ultrabroadband 50-130 THz pulses generated via phase-matched difference frequency mixing in LiIO 3 | |
| Pfeiffer et al. | Terahertz detection by upconversion to the near-infrared using picosecond pulses | |
| JP2012068621A (ja) | テラヘルツ波発生素子、テラヘルツ波検出素子、及びテラヘルツ時間領域分光装置 | |
| US8901477B2 (en) | Electromagnetic wave detection device | |
| JP4838111B2 (ja) | 電磁波検出器及び電磁波検出システム | |
| JP2019132595A (ja) | テラヘルツ波分光計測装置 | |
| Khan et al. | Ultrasensitive, room temperature detection of THz radiation using nonlinear parametric conversion | |
| Lorenc et al. | Mid-infrared Kerr index evaluation via cross-phase modulation with a near-infrared probe beam | |
| Minamide et al. | Frequency-agile terahertz-wave generation and detection using a nonlinear optical conversion, and their applications for imaging | |
| JP5963080B2 (ja) | 電磁波検出方法及び電磁波検出装置 | |
| JP2013174513A (ja) | 電磁波検出方法及び電磁波検出装置 | |
| JP2010181193A (ja) | 光検出装置及び光検出方法 | |
| Lee et al. | Coherent detection of multicycle terahertz pulses generated in periodically inverted GaAs structures | |
| Patokoski et al. | Broadband mid-infrared cavity ringdown spectroscopy with frequency upconversion detector | |
| JP2019144435A (ja) | テラヘルツ波検出素子 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160408 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20161221 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170206 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170405 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20170929 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20171005 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6229923 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |