JP2008002960A - Method for measuring solid carrier mobility - Google Patents

Method for measuring solid carrier mobility Download PDF

Info

Publication number
JP2008002960A
JP2008002960A JP2006172865A JP2006172865A JP2008002960A JP 2008002960 A JP2008002960 A JP 2008002960A JP 2006172865 A JP2006172865 A JP 2006172865A JP 2006172865 A JP2006172865 A JP 2006172865A JP 2008002960 A JP2008002960 A JP 2008002960A
Authority
JP
Japan
Prior art keywords
plasmon
time
phonon
mode
longitudinal
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.)
Granted
Application number
JP2006172865A
Other languages
Japanese (ja)
Other versions
JP4831482B2 (en
Inventor
Muneaki Hase
宗明 長谷
Masahiro Kitajima
正弘 北島
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Institute for Materials Science
Original Assignee
National Institute for Materials Science
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Institute for Materials Science filed Critical National Institute for Materials Science
Priority to JP2006172865A priority Critical patent/JP4831482B2/en
Publication of JP2008002960A publication Critical patent/JP2008002960A/en
Application granted granted Critical
Publication of JP4831482B2 publication Critical patent/JP4831482B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method capable of measuring a solid carrier mobility in the dynamic state accurately in a short time. <P>SOLUTION: The method is characterized as follows: an actual time vibration structure by plasmon in a reflectivity change or in a transmissivity change is measured by using a pump-probe spectroscopy, relative to a solid which is a measuring object wherein phonon is not bonded to the plasmon; fitting of time domain data is performed by using expression (1) to the measured data, to thereby determine a relaxation time τ of the plasmon; and the value of the determined relaxation time τ is substituted for expression (2), to thereby determine the solid carrier mobility μ. In the expressions, A is the amplitude of vibration by the plasmon, τ is the relaxation time of the plasmon, t is a time, ω<SB>p</SB>is a frequency of the plasmon, ϕ is an initial phase of the vibration by the plasmon, e is a charge of the electron, and m* is an effective mass of the electron. There are two other procedures. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

本発明は、固体のキャリア移動度測定方法に関するものである。より詳しくは、本発明は、固体のレーザーパルス励起により生じた非平衡プラズモン又はプラズモンライクな縦波光学(LO)フォノン−プラズモン結合モードの振動による固体の屈折率変化を実時間領域で測定し、その振動の減衰時定数から固体の移動度を直接決定する方法に関するもので、特に半導体、金属などの材料からなるナノ構造デバイスの電気的特性評価法として応用できるものである。   The present invention relates to a method for measuring solid carrier mobility. More specifically, the present invention measures the change in the refractive index of a solid in the real-time domain due to vibrations of a nonequilibrium plasmon or plasmon-like longitudinal wave optical (LO) phonon-plasmon coupled mode caused by laser pulse excitation of the solid, The present invention relates to a method for directly determining the mobility of a solid from its vibration damping time constant, and is particularly applicable as a method for evaluating electrical characteristics of nanostructure devices made of materials such as semiconductors and metals.

近年、デバイスの極微細化が進み、ナノメートルサイズのデバイスが実現しつつある。ここで電子の平均自由行程(電子が散乱されるまでに進む平均距離;Siでは10〜50nm程度)がデバイスサイズ以下になると、電子が格子振動や格子欠陥に散乱されない領域(すなわちバリスティック伝導領域)になる。この領域が実現できれば、電界効果トランジスタ(FET)等のデバイス中で、電子はソースからドレインまで散乱されることなく進み、スイッチング速度が劇的に速いトランジスタが実現できる。   In recent years, the miniaturization of devices has progressed, and nanometer-sized devices are being realized. Here, when the average free path of electrons (average distance traveled until electrons are scattered; about 10 to 50 nm for Si) becomes smaller than the device size, the region where electrons are not scattered by lattice vibration or lattice defects (that is, ballistic conduction region) )become. If this region can be realized, in a device such as a field effect transistor (FET), electrons proceed without being scattered from the source to the drain, and a transistor with a dramatically high switching speed can be realized.

現在、バリスティック伝導に関する研究は、ほとんどが電子線リソグラフィー等の微細加工技術に依存したデバイスサイズを小さくする手法を取っている。一方、デバイスサイズを小さくする代わりに、電子の平均自由行程を長くする手法もある。例えば極低温にしてフォノン散乱の影響を無くすか、あるいは分子線エピタキシー(MBE)等の技術を駆使して高純度で欠陥の少ない結晶を生成することである。   Currently, most researches on ballistic conduction take a method of reducing the device size depending on microfabrication techniques such as electron beam lithography. On the other hand, there is also a method of increasing the mean free path of electrons instead of decreasing the device size. For example, the effect of phonon scattering can be eliminated by making the temperature extremely low, or a crystal having high purity and few defects can be produced by using a technique such as molecular beam epitaxy (MBE).

しかし、デバイスを室温で使用するためには前者の極低温という選択は難しく、また後者の結晶性についても、欠陥密度をゼロにすることはほとんど不可能である。このような理由で、微細加工技術によるデバイスサイズのナノ構造化が進んできた訳である。   However, in order to use the device at room temperature, it is difficult to select the former extremely low temperature, and it is almost impossible to make the defect density zero even for the latter crystallinity. For this reason, nano-structuring of device size has been advanced by microfabrication technology.

以上のような背景の中で重要になるのが、半導体、金属などのナノ構造における移動度の評価である。通常、半導体の移動度測定は、ホール効果を利用することが一般的である(例えば、非特許文献1など)。このホール効果は、1879年にE. H. Hallによって発見された現象であり、半導体や金属に電流が流れている時、電流に直角に磁場を印加すると両者に直角な方向に起電力が生じるというものである。このホール効果の実験では、生じた起電力と印加した磁場と電流密度との関係からホール係数が求まり、このホール係数から移動度を導き出すことができる。しかし、ナノ構造体については、ナノメートル(10−9m)オーダーの空間分解能でのホール測定は困難であると考えられる。 What is important in the above background is the evaluation of mobility in nanostructures such as semiconductors and metals. In general, the mobility measurement of a semiconductor generally uses the Hall effect (for example, Non-Patent Document 1). This Hall effect is a phenomenon discovered by EH Hall in 1879. When a current flows through a semiconductor or metal, an electromotive force is generated in a direction perpendicular to both when a magnetic field is applied perpendicular to the current. is there. In the Hall effect experiment, the Hall coefficient is obtained from the relationship between the generated electromotive force, the applied magnetic field, and the current density, and the mobility can be derived from the Hall coefficient. However, for nanostructures, it is considered difficult to measure holes with a spatial resolution on the order of nanometers (10 −9 m).

そのような背景で近年、レーザーを用いたラマン散乱を利用した光学的手法によっても、半導体の移動度を求めることができることが分かってきた(例えば、非特許文献2、3など)。空間分解能約1μmを達成できる顕微ラマン分光(例えば、非特許文献4など)や、空間分解能1μm以下を達成できる近接場ラマン分光(例えば、非特許文献5など)を利用すれば、ナノ構造体において非常に局所的な(ナノ構造体1個または数個の)キャリア移動度を評価できる可能性を持っている。   Against this background, it has recently been found that the mobility of a semiconductor can also be obtained by an optical method using Raman scattering using a laser (for example, Non-Patent Documents 2 and 3). By using microscopic Raman spectroscopy (for example, Non-Patent Document 4) that can achieve a spatial resolution of about 1 μm and near-field Raman spectroscopy (for example, Non-Patent Document 5) that can achieve a spatial resolution of 1 μm or less, It has the potential to evaluate very local carrier mobility (one or several nanostructures).

このラマン散乱を用いた方法では、まず、極性のある半導体(SiCなど)中に特有である、縦波光学(LO)フォノン(以下、LOフォノンとも称する)とプラズモン(電子の集団振動)とが結合して生じたLOフォノン−プラズモン結合モードの上の分枝(L)と下の分枝(L)のうち、プラズモンライクな分枝(L又はL)に相当するスペクトルを、カーブフィッティングすることによって、プラズモンの減衰定数γを求める。
ここで、プラズモンライクなモードとは、図1に示すように、結合モードの上の分枝(L)と下の分枝(L)のうち、裸の(フォノンと結合していない)プラズモンの周波数を示す直線Eに近い状態(図1中で、AやBの領域)を指す。CやDの領域は結合モードの減衰がフォノンの緩和に支配される「フォノンライク」な領域である。
In this method using Raman scattering, first, there are longitudinal wave optical (LO) phonons (hereinafter also referred to as LO phonons) and plasmons (collective vibrations of electrons), which are peculiar in polar semiconductors (such as SiC). The spectrum corresponding to the plasmon-like branch (L + or L ) of the upper branch (L + ) and the lower branch (L ) of the combined LO phonon-plasmon coupled mode is A plasmon damping constant γ is obtained by curve fitting.
Here, as shown in FIG. 1, the plasmon-like mode is a bare (not coupled to a phonon) of the upper branch (L + ) and the lower branch (L ) of the coupled mode. This indicates a state close to a straight line E indicating the frequency of plasmons (A and B regions in FIG. 1). Regions C and D are “phonon-like” regions where the coupling mode attenuation is governed by phonon relaxation.

また、カーブフィッティングでは、一般に非対称的なスペクトル波形になるLOフォノン−プラズモン結合モードを対象とするため、下記のような特殊な関数を用いる(例えば、前記非特許文献2、3など)。   In curve fitting, since the LO phonon-plasmon coupled mode, which generally has an asymmetric spectral waveform, is targeted, the following special functions are used (for example, Non-Patent Documents 2 and 3).

ここで、Γはフォノンの減衰定数、γはプラズモンの減衰定数、ω,ω,ωはそれぞれTOフォノン周波数、LOフォノン周波数、プラズマ周波数、eは電子電荷、Nは電子密度、mは電子の有効質量、εは高周波数極限の誘電率、εは真空の誘電率である。このようにして求められたプラズモンの減衰定数γを次式に代入することによって、移動度が計算される。 Here, Γ is a phonon attenuation constant, γ is a plasmon attenuation constant, ω T , ω L , and ω P are a TO phonon frequency, a LO phonon frequency, and a plasma frequency, respectively, e is an electron charge, N is an electron density, and m *. Is the effective mass of electrons, ε is the dielectric constant at the high frequency limit, and ε 0 is the dielectric constant of the vacuum. The mobility is calculated by substituting the plasmon attenuation constant γ thus obtained into the following equation.

このラマンを用いた方法では、測定試料に電極を付けたり、磁場を印加したりせずに、非接触に光学的手法で移動度を求めることができるので、極性がある半導体については有効であると思われる。
「半導体の物理」培風館, 西沢潤一編, 御子柴宣夫著, 1995年 Journal of Applied Physics, Vol. 78, No.3, pp.1996, 1995年 Journal of Applied Physics, Vol. 90, No.10, pp.5211, 2001年 Journal of Applied Physics, Vol. 78, No.5, pp.3357, 1995年 Journal of Chemical Physics, Vol. 117, No.3, pp.1296, 2002年
This method using Raman is effective for semiconductors with polarity, because the mobility can be obtained in a non-contact optical manner without attaching an electrode to the measurement sample or applying a magnetic field. I think that the.
"Semiconductor Physics" Baifukan, edited by Junichi Nishizawa, Nobuo Mikoshiba, 1995 Journal of Applied Physics, Vol. 78, No.3, pp.1996, 1995 Journal of Applied Physics, Vol. 90, No.10, pp.5211, 2001 Journal of Applied Physics, Vol. 78, No.5, pp.3357, 1995 Journal of Chemical Physics, Vol. 117, No.3, pp.1296, 2002

しかし、上記の方法は、ラマンスペクトルという周波数領域の測定手法で得られたスペクトルをかなり大がかりな数式群(式(8)〜式(12))でフィットし、元々時間領域の定数であるプラズモンの緩和時間τの逆数に相当する減衰定数γ=1/τを求めるので、非常に複雑なスペクトルのカーブフィットを必要とし、データ処理に時間がかかるという問題点があった。また、上記の方法は、静的な状態のキャリア移動度を測定しているが、動的な状態、すなわち、例えば光励起によって生じたキャリアが励起された表面から試料の深さ方向に伝導していく状態のキャリア移動度を感度良く測定することが望まれていた。   However, the above method fits the spectrum obtained by the frequency domain measurement technique called Raman spectrum with a fairly large mathematical group (formulas (8) to (12)), and the plasmon constant, which is originally a constant in the time domain. Since the attenuation constant γ = 1 / τ corresponding to the reciprocal of the relaxation time τ is obtained, a very complicated curve fitting of the spectrum is required, and there is a problem that the data processing takes time. In the above method, the carrier mobility in a static state is measured. However, in a dynamic state, for example, carriers generated by photoexcitation are conducted from the excited surface in the depth direction of the sample. It has been desired to measure the carrier mobility in various states with high sensitivity.

本発明は、以上の通りの事情に鑑みてなされたものであり、動的な状態のキャリア移動度を短時間で且つ正確に測定することができる固体のキャリア移動度測定方法を提供することを課題としている。   The present invention has been made in view of the circumstances as described above, and provides a solid carrier mobility measurement method capable of accurately measuring carrier mobility in a dynamic state in a short time. It is an issue.

本発明の固体のキャリア移動度測定方法は、上記の課題を解決するために、以下のことを特徴としている。   The solid carrier mobility measuring method of the present invention is characterized by the following in order to solve the above problems.

第1に、フォノンとプラズモンが結合しない測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンによる実時間振動構造を測定し、その測定データに対し、下記式(1)を用いて時間領域データのフィッティングを行うことにより、プラズモンの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めること。   First, using a pump-probe spectroscopy method, the real-time vibration structure due to plasmon in the change in reflectance or transmittance is measured for a solid that is a measurement target in which phonons and plasmons do not bind. The plasmon relaxation time τ is obtained by fitting time domain data using the following equation (1), and the value of the obtained relaxation time τ is substituted into the following equation (2) to obtain the solid carrier mobility. Find μ.

(上記式中、Aはプラズモンによる振動の振幅、τはプラズモンの緩和時間、tは時間、ωPはプラズモンの振動数、φはプラズモンの振動の初期位相、eは電子の電荷、m
電子の有効質量である)
(Where A is the amplitude of plasmon oscillation, τ is the plasmon relaxation time, t is the time, ω P is the plasmon frequency, φ is the initial phase of the plasmon oscillation, e is the charge of the electron, m * is (Effective mass of electrons)

第2に、フォノンとプラズモンが結合する測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォ
ノン−プラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(3)を用いて時間領域データのフィッティングを行うことにより、プラズモンライクな縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めること。
Secondly, real-time vibration of a solid, which is a measurement object in which phonon and plasmon are combined, by a plasmon-like longitudinal wave optical phonon-plasmon coupled mode in reflectance change or transmittance change using pump-probe spectroscopy. By measuring the structure and fitting the time domain data to the measured data using the following equation (3), the relaxation time τ of the plasmon-like longitudinal wave optical phonon-plasmon coupled mode is obtained, and the obtained relaxation is obtained. The solid carrier mobility μ is obtained by substituting the value of time τ into the following equation (2).

(上記式中、tは式(1)と同様、AL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードによる振動の振幅、τ(=τL+)は縦波光学フォノン−プラズモ
ン結合モードのプラズモンライクモードの緩和時間、ΩL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードの振動数、φL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードの振動の初期位相、Bは縦波光学フォノンによる振動の振幅、τLOは縦波光学フォノンの緩和時間、ΩLOは縦波光学フォノンの振動数、φLOは縦波光学フォノンの振動の初期位相、Cは縦波光学フォノン−プラズモン結合モードのフォノンライクモードによる振動の振幅、τL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの緩和時間、ΩL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの振動数、φL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの振動の初期位相、eは電子の電荷、mは電子の有効質量である)
(In the above equation, t is the same as in equation (1), A L + is the amplitude of vibration in the plasmon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, and τ (= τ L + ) is the longitudinal wave optical phonon-plasmon coupled mode. Plasmon-like mode relaxation time, Ω L + is the frequency of the plasmon-like mode in the longitudinal wave optical phonon-plasmon coupled mode, φ L + is the initial phase of the plasmon-like mode vibration in the longitudinal wave optical phonon-plasmon coupled mode, and B is the longitudinal phase Amplitude of vibration by wave optical phonon, τ LO is relaxation time of longitudinal optical phonon, Ω LO is frequency of longitudinal optical phonon, φ LO is initial phase of longitudinal optical phonon vibration, C is longitudinal optical phonon − The amplitude of vibration due to the phonon-like mode of the plasmon coupled mode, τ L− is the phonon line of the longitudinal optical phonon-plasmon coupled mode. Mode relaxation time, Ω L− is the frequency of the phonon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, φ L− is the initial phase of the phonon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, and e is the electron Charge, m * is the effective mass of electrons)

第3に、フォノンとプラズモンが結合する測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォノン−プラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(4)を用いてウェーブレット変換を行うことにより、プラズモンライクな縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めること。   Third, real-time vibrations in a plasmon-like longitudinal-wave optical phonon-plasmon coupled mode in a change in reflectance or transmittance using a pump-probe spectroscopy for a solid that is a measurement target in which phonons and plasmons are combined. The structure is measured, and the measurement data is subjected to wavelet transform using the following equation (4) to obtain the relaxation time τ of the plasmon-like longitudinal wave optical phonon-plasmon coupled mode, and the obtained relaxation time τ Obtain the solid carrier mobility μ by substituting the value into the following formula (2).

(上記式中、Ψ(x)はマザーウェーブレットと称される、波束を表す関数、xは入力信号で時間の関数、bはこの波束の時間遅延を表す変数、aはこの波束の時間軸に対する伸縮度、tは時間、f(t)は被変換関数である) (In the above equation, Ψ (x) is a function that represents a wave packet called a mother wavelet, x is a function of time as an input signal, b is a variable that represents a time delay of the wave packet, and a is a time axis of the wave packet. Degree of elasticity, t is time, f (t) is the transformed function)

第4に、上記第3の発明において、ウェーブレット変換のマザーウェーブレットとして下記式(5)で示す、ガウス関数を変形したGaborのマザーウェーブレットを用いること
Fourth, in the third invention, a Gabor mother wavelet obtained by transforming a Gaussian function expressed by the following equation (5) is used as a mother wavelet for wavelet transformation.

(上記式中、σはガウス関数の半値全幅、xは入力信号で時間の関数である)
を用いること。
(In the above formula, σ is the full width at half maximum of the Gaussian function, x is the input signal and a function of time)
Use.

第5に、フォノンとプラズモンが結合している測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォノン−プラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(6)を用いて短時間フーリエ変換を行うことにより、プラズモンライクな縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めること。   Fifth, by using pump-probe spectroscopy for a solid, which is a measurement object in which phonon and plasmon are combined, an actual measurement using a plasmon-like longitudinal wave optical phonon-plasmon coupled mode in reflectance change or transmittance change is performed. The time oscillation structure was measured and the relaxation time τ of the plasmon-like longitudinal wave optical phonon-plasmon coupled mode was determined by performing a short-time Fourier transform on the measured data using the following formula (6). By substituting the value of the relaxation time τ into the following formula (2), the carrier mobility μ of the solid is obtained.

(上記式中、Ψ(x)は窓関数、xは入力信号で時間の関数、bは時間遅延、ωは周波数、tは時間、f(t)は被変換関数である) (In the above formula, Ψ (x) is a window function, x is an input signal and a function of time, b is a time delay, ω is a frequency, t is time, and f (t) is a transformed function)

第6に、上記第5の発明において、短時間フーリエ変換の窓関数としてガウス関数を用い、下記式(7)により短時間フーリエ変換を行うこと。   Sixth, in the fifth invention, a Gaussian function is used as the window function of the short-time Fourier transform, and the short-time Fourier transform is performed by the following equation (7).

(上記式中、σはガウス関数の半値全幅である) (In the above formula, σ is the full width at half maximum of the Gaussian function)

本発明によれば、電子又は正孔の集団振動であるプラズモン又はプラズモンライクなLOフォノン−プラズモン結合モードの実時間振動波形を測定することにより、緩和時間を直接求め、それによって固体のキャリア移動度を求めるようにしたので、動的な状態のキャリア移動度を短時間で且つ正確に測定することができる固体のキャリア移動度測定方法を提供することが可能となる。   According to the present invention, the relaxation time is directly obtained by measuring the real-time oscillation waveform of the plasmon or plasmon-like LO phonon-plasmon coupled mode, which is the collective oscillation of electrons or holes, and thereby the carrier mobility of the solid. Therefore, it is possible to provide a solid carrier mobility measuring method capable of accurately measuring carrier mobility in a dynamic state in a short time.

本発明は、特に半導体、金属などの固体材料からなるナノ構造デバイスの電気的特性評
価法として応用することができる。固体、特にナノ構造の固体中の動的な状態あるいは過程は、例えば、半導体ナノ構造を用いたキャリアによる超高速スイッチを開発する際に起きている現象であり、半導体ナノデバイスの電気的特性評価には移動度の測定が必須であり、これからの超高速デバイスの開発に非常に有用になると期待される。本発明によれば、様々な条件(温度や欠陥密度など)下におけるキャリア移動度の評価を行うことができ、HEMT(High Electron Mobility Transistor)などの高速電子移動によるテラヘルツ
(1012Hz)電子デバイスの開発へ利用、超高速データ通信や超省電力機器の実現など、大容量データ送受信や省エネルギー問題に寄与することができる。
The present invention can be applied particularly as a method for evaluating electrical characteristics of nanostructure devices made of solid materials such as semiconductors and metals. The dynamic state or process in solids, especially nanostructured solids, is a phenomenon that occurs, for example, when developing ultrafast switches with carriers using semiconductor nanostructures, and the electrical characterization of semiconductor nanodevices Therefore, mobility measurement is essential, and it is expected to be very useful for the development of ultra high-speed devices in the future. According to the present invention, carrier mobility can be evaluated under various conditions (temperature, defect density, etc.), and a terahertz (10 12 Hz) electronic device using high-speed electron transfer such as HEMT (High Electron Mobility Transistor). It can contribute to the transmission / reception of large-capacity data and energy saving problems such as the development of high-speed data communication and the realization of ultra-power-saving equipment.

以下、本発明の実施の形態について詳細に説明する。   Hereinafter, embodiments of the present invention will be described in detail.

本発明では、従来法のように周波数領域でプラズモン又はプラズモンライクなLOフォノン−プラズモン結合モードの減衰定数γを求めるのではなく、時間領域で直接的にプラズモン又はプラズモンライクなLOフォノン−プラズモン結合モードの緩和時間τを求めるという全く新しい発想で、金属や半導体(極性、非極性を問わない)バルク結晶や薄膜、および金属や半導体などを含む固体のナノ構造体などにおける移動度を測定するものである。   In the present invention, the plasmon or plasmon-like LO phonon-plasmon coupled mode attenuation constant γ is not obtained in the frequency domain as in the conventional method, but the plasmon or plasmon-like LO phonon-plasmon coupled mode directly in the time domain. It is a completely new idea to determine the relaxation time τ of metal and semiconductors (both polar and non-polar) bulk crystals and thin films, and the mobility in solid nanostructures including metals and semiconductors. is there.

本発明において、プラズモン又はプラズモンライクなLOフォノン−プラズモン結合モードの緩和時間τを求めるには、ポンプ−プローブ分光法と呼ばれる一般的な時間分解測定法(例えば、時間分解反射率変化測定法や時間分解透過率変化測定法:特開2002-214137号公報などを参照)を用いて、反射率(あるいは透過率)変化においてプラズモン又は
プラズモンライクなLOフォノン−プラズモン結合モードによる実時間振動構造を測定することにより行う。ただし、ここで、プラズモン又はプラズモンライクなLOフォノン−プラズモン結合モードの緩和時間τを求める観測手段は大きく分けて2種類存在する。
In the present invention, in order to obtain the relaxation time τ of the plasmon or plasmon-like LO phonon-plasmon coupled mode, a general time-resolved measurement method called pump-probe spectroscopy (for example, time-resolved reflectance change measurement method or time Measurement of real-time vibration structure by plasmon or plasmon-like LO phonon-plasmon coupled mode in reflectivity (or transmissivity) change using a method for measuring a change in decomposed transmissivity (see JP 2002-214137 A). By doing. However, there are roughly two types of observation means for obtaining the relaxation time τ of the plasmon or plasmon-like LO phonon-plasmon coupled mode.

第1は、反射率(あるいは透過率)変化測定により、裸のプラズモンによる信号を取得する方法である。この場合、プラズモンの周波数はよく知られた関係式(前記式(13))で与えられるので、試料のキャリア密度に比例してプラズモンの周波数は増加することになる。従って、高キャリア密度になればなるほど、プラズモンの周波数は高くなり、高い時間分解能が必要になる。この手法は、プラズモンが他の素励起(特にLOフォノン)と結合しないような状態で有効である。このようなものとしては、非極性物質でありプラズモンがLOフォノンと結合しない金属、非極性半導体(Si、Geなど)が考えられる。   The first is a method of acquiring a signal by naked plasmon by measuring a change in reflectance (or transmittance). In this case, since the plasmon frequency is given by a well-known relational expression (formula (13)), the plasmon frequency increases in proportion to the carrier density of the sample. Accordingly, the higher the carrier density, the higher the plasmon frequency and the higher the time resolution. This technique is effective in a state where plasmons do not couple with other elementary excitations (especially LO phonons). Examples of such materials include metals and nonpolar semiconductors (Si, Ge, etc.) that are nonpolar substances and plasmons do not bind to LO phonons.

第2は、極性物質(GaAs、InP、SiCなど)において、積極的にプラズモンとLOフォノンの結合を利用する手法である。この場合は、ラマン散乱の場合と同じように、LOフォノン−プラズモン結合モードに相当する実時間信号を取得することになる。ここで注意すべき点は、LOフォノン−プラズモン結合モードにもフォノンライクモード(モードの減衰時間はフォノンの減衰時間で決まる)とプラズモンライクモード(モードの減衰時間はプラズモンの減衰時間で決まる)が存在するので、プラズモンライクなモードを観測することにより、その減衰時間を求めれば、近似的にプラズモンの減衰時間と見なすことができる(τ=τL+)。ここで、プラズモンライクなモードとは、ラマンの場合
と同様に図1に示すように、結合モードの上の分枝(L)と下の文枝(L)のうち、裸の(フォノンと結合していない)プラズモンの周波数を示す直線Eに近い状態を指す。すなわち、図1中で、AやBの領域を指す。CやDの領域は結合モードの減衰がフォノンの緩和に支配される「フォノンライク」な領域である。
The second is a method that positively utilizes the coupling between plasmons and LO phonons in polar substances (GaAs, InP, SiC, etc.). In this case, as in the case of Raman scattering, a real-time signal corresponding to the LO phonon-plasmon coupled mode is acquired. It should be noted that the LO phonon-plasmon coupled mode also has a phonon-like mode (mode decay time is determined by phonon decay time) and plasmon-like mode (mode decay time is determined by plasmon decay time). Therefore, if the decay time is obtained by observing a plasmon-like mode, it can be approximately regarded as the decay time of plasmon (τ = τ L + ). Here, the plasmon-like mode is a bare (phonon and phonon) of the upper branch (L + ) and the lower branch (L ) of the coupled mode as shown in FIG. This indicates a state close to a straight line E indicating the frequency of plasmons (not coupled). That is, it refers to the area A or B in FIG. Regions C and D are “phonon-like” regions where the coupling mode attenuation is governed by phonon relaxation.

いずれの場合も、実時間領域で観測されたプラズモンによる振動の減衰時間を求めなけ
ればならない。上記第1の場合は、単純な減衰振動になることから、データを
In any case, the decay time of vibration due to plasmons observed in the real-time region must be obtained. In the first case, the data is simply a damped vibration.

のような減衰調和振動の式でフィットすれば、プラズモンの緩和時間τが求まる。ここで、Aはプラズモンによる振動の振幅、τはプラズモンの緩和時間、ωPはプラズモンの振
動数、φはプラズモンの振動の初期位相、tは時間である。
The plasmon relaxation time τ can be obtained by fitting with a damped harmonic vibration equation such as Here, A is the amplitude of the plasmon vibration, τ is the plasmon relaxation time, ω P is the plasmon frequency, φ is the initial phase of the plasmon vibration, and t is the time.

しかし、上記第2の場合には、プラズモンライクなLOフォノン−プラズモン結合モードの信号は、同時に観測されるLOフォノンの信号(このLOフォノンは主に表面近傍の空乏層に存在している)に重畳することになり、上記式(1)のような単純な式ではフィットできない場合がある。   However, in the second case, the plasmon-like LO phonon-plasmon coupled mode signal is an LO phonon signal observed at the same time (this LO phonon is mainly present in the depletion layer near the surface). In some cases, a simple expression such as the above expression (1) cannot be fitted.

そこで、解決手段として下記の3つの手法をとる。まず、上記式(1)を、観測されうる全てのモードについて加算した次式、   Therefore, the following three methods are taken as solution means. First, the following equation obtained by adding the above equation (1) for all modes that can be observed:

で時間領域の振動波形をフィットすることである。ここで、第一項のAL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモード(例えばn型GaAsにおいてN>1×1018cm−3では、図1中AのL+に相当)による振動の振幅、τ(=τL+
はL+モードの緩和時間、ΩL+はL+モードの振動数、φL+はL+モードの振動の初期
位相、第二項のBはLOフォノンによる振動の振幅、τLOはLOフォノンの緩和時間、ΩLOはLOフォノンの振動数、φLOはLOフォノンの振動の初期位相である。また、第三項のCはLOフォノン−プラズモン結合モードのフォノンライクモード(例えばn型GaAsにおいてN>1×1018cm−3では、図1中DのLに相当)による振動の振幅、τL−はLモードの緩和時間、ΩL−はLモードの振動数、φL−はLモードの振動の初期位相である。
To fit the time domain vibration waveform. Here, A L + in the first term is based on a plasmon-like mode of the longitudinal wave optical phonon-plasmon coupling mode (for example, in N-type GaAs, N> 1 × 10 18 cm −3 corresponds to L + of A in FIG. 1). Amplitude of vibration, τ (= τ L + )
Is the relaxation time of the L + mode, Ω L + is the frequency of the L + mode, φ L + is the initial phase of the vibration of the L + mode, B in the second term is the amplitude of the vibration due to the LO phonon, and τ LO is the relaxation of the LO phonon. Time, Ω LO is the LO phonon frequency, and φ LO is the initial phase of the LO phonon vibration. C in the third term is the amplitude of vibration due to the phonon-like mode of the LO phonon-plasmon coupled mode (for example, N> 1 × 10 18 cm −3 in n-type GaAs corresponds to L in FIG. 1), τ L− is the relaxation time of the L mode, Ω L− is the frequency of the L mode, and φ L− is the initial phase of the vibration of the L mode.

もし、上記式(3)で振動波形を上手くフィットできない場合は、次にウェーブレット変換と呼ばれる独特な手法を用いる(「ウェーブレットビギナーズガイド」榊原進著(東京電機大学出版局,1995年)及び特開2003-296301号公報などを参照)。これは、数学的
には確立されて来ている手法であるが、物理・計測分野ではほとんど利用されていない手法である。具体的には、ある被変換関数f(t)に対して、次式
If the vibration waveform cannot be fit well with the above equation (3), then a unique method called wavelet transform is used ("Wavelet Beginners Guide" Susumu Sugawara (Tokyo Denki University Press, 1995) and JP (See 2003-296301). This is a method that has been established mathematically, but is rarely used in the fields of physics and measurement. Specifically, for a given function f (t),

で定義される。ここで、Ψ(x)はマザーウェーブレット(xは入力信号で時間の関数)と呼ばれ、波束を表す関数である。bはこの波束の時間遅延を表す変数、aはこの波束の時間軸に対する伸縮度、すなわち波束の包絡線の時間幅を表す変数である。変換後の信号
強度は、[b,1/a]平面すなわち[時間,周波数]平面で二次元の等高線図になる。このよ
うな変換を施すことにより、同じ時間軸に別々の周波数で埋もれていた信号成分を別々に取り出すことができる。マザーウェーブレットの関数型としては、例えばガウス関数などがある。これによって、LOフォノンとは周波数の異なるプラズモンライクなLOフォノン−プラズモン結合モードの信号の時間変化を得ることができる。このようにして得られたプラズモンライクモードの緩和時間τは、ウェーブレット信号上でプラズモンライクモードのピーク周波数を見つけ、この周波数成分の時間遅延方向の緩和時間を指数関数e−t/τを用いてフィットし求めることができる。
Defined by Here, Ψ (x) is called a mother wavelet (x is a function of time with an input signal), and is a function representing a wave packet. b is a variable representing the time delay of the wave packet, and a is a variable representing the degree of stretch of the wave packet with respect to the time axis, that is, the time width of the envelope of the wave packet. The signal intensity after conversion becomes a two-dimensional contour map in the [b, 1 / a] plane, that is, the [time, frequency] plane. By performing such conversion, signal components buried at different frequencies on the same time axis can be extracted separately. Examples of the function type of the mother wavelet include a Gaussian function. Thereby, it is possible to obtain a time change of a signal in a plasmon-like LO phonon-plasmon coupled mode having a frequency different from that of the LO phonon. The relaxation time τ of the plasmon-like mode obtained in this way is obtained by finding the peak frequency of the plasmon-like mode on the wavelet signal and using the exponential function e −t / τ to determine the relaxation time of this frequency component in the time delay direction. Fit and ask.

また一方、ウェーブレット変換を用いず、ある窓関数Ψ(x)を用いた短時間フーリエ変換(「ウェーブレットビギナーズガイド」榊原進著(東京電機大学出版局,1995年)などを参照)を使って、プラズモンライクモードの緩和時間τの取得も可能である。この場合、次式   On the other hand, without using the wavelet transform, using a short-time Fourier transform using a window function Ψ (x) (see “Wavelet Beginners Guide” Susumu Sugawara (Tokyo Denki University Press, 1995) etc.) It is also possible to acquire the relaxation time τ of the plasmon-like mode. In this case,

で表される変換を被変換関数f(t)に対して行い、様々な時間遅延bに対する短時間フーリエスペクトルを得る。このスペクトルは遅延時間と共に減衰していくので、プラズモンライクモードの緩和時間τがその減衰から求まる。なお、上記式中、ωは周波数である。 Is applied to the function f (t) to be converted to obtain short-time Fourier spectra for various time delays b. Since this spectrum attenuates with the delay time, the relaxation time τ of the plasmon-like mode is obtained from the attenuation. In the above formula, ω is a frequency.

以上のようにして得られた緩和時間τを、前記式(14)ではなく、より一般的な下記式(2)(「半導体の物理」培風館,西沢潤一編,御子柴宣夫著,1995年などを参照)   The relaxation time τ obtained as described above is not the above formula (14), but the following general formula (2) (“Semiconductor Physics” Baifukan, edited by Junichi Nishizawa, Nobuo Mikoshiba, 1995) reference)

に代入して移動度μを直接求めることができる。 The mobility μ can be directly obtained by substituting into.

このように、周波数領域の分光法では、非常に複雑なデータフィットを要した移動度の導出が、本発明によれば、時間領域で直接プラズモンもしくはプラズモンライクなLOフォノン−プラズモン結合モードの緩和時間τを求めることにより、非常に簡潔に行うことができるわけである。   Thus, in frequency domain spectroscopy, the derivation of mobility requiring a very complex data fit, according to the present invention, is the relaxation time τ of a direct plasmon or plasmon-like LO phonon-plasmon coupled mode in the time domain. Can be done very simply.

次に、本発明を実施例によりさらに詳細に説明する。
<実施例>
図2は、本発明による固体のキャリア移動度測定方法の原理を示したものである。フェムト秒パルスレーザーを光源としたポンプ−プローブ分光法による時間分解反射率(透過率)測定装置Aにおいて測定された試料のプラズモン信号の種類によって、3種類のデータ処理(図2のB、C、D)が考えられる。測定装置Aには反射型、透過型などがある。
Next, the present invention will be described in more detail with reference to examples.
<Example>
FIG. 2 shows the principle of the solid carrier mobility measuring method according to the present invention. Depending on the type of plasmon signal of the sample measured in the time-resolved reflectance (transmittance) measuring device A by pump-probe spectroscopy using a femtosecond pulse laser as a light source, three types of data processing (B, C, D) is conceivable. The measuring apparatus A includes a reflection type and a transmission type.

まず第1は、上記式(1)により時間領域の振動波形が簡単にフィットできる場合が図2のBである。このように単純な減衰調和振動は上記式(1)により良く再現できるため、プラズモン振動の緩和時間τをこのフィッティングで求めることができる。求めた緩和
時間τは直ちに上記式(2)に代入され、移動度μが求められる。
First, FIG. 2B shows a case where the vibration waveform in the time domain can be easily fitted by the above equation (1). Since simple damped harmonic vibration can be well reproduced by the above equation (1), the relaxation time τ of plasmon vibration can be obtained by this fitting. The obtained relaxation time τ is immediately substituted into the above equation (2), and the mobility μ is obtained.

第2は、図2のCに示すように、時間分解反射率(透過率)の信号が単純な減衰調和振動ではないが、プラズモンライクなLOフォノン−プラズモン結合モードとLOフォノンなど複数の減衰調和振動の和(上記式(3))で良く再現できる場合である。この場合、得られたプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間τをこのフィッティングで求めることができる。求めた緩和時間τは直ちに上記式(2)に代入され、移動度μが求められる。 Second, as shown in FIG. 2C, the signal of time-resolved reflectance (transmittance) is not a simple damped harmonic oscillation, but a plurality of damped harmonics such as a plasmon-like LO phonon-plasmon coupled mode and LO phonon. This is a case where the vibration can be well reproduced by the sum of vibrations (formula (3) above). In this case, the relaxation time τ of the obtained plasmon-like LO phonon-plasmon coupled mode (L + ) can be obtained by this fitting. The obtained relaxation time τ is immediately substituted into the above equation (2), and the mobility μ is obtained.

第3は、図2のDのように非常に複雑な振動構造で、上記式(1)でも上記式(3)でもデータをフィットできない場合である。この場合は、ウェーブレット変換(上記式(4))もしくは短時間フーリエ変換(上記式(6))を用いることになる。   The third is a case where the vibration structure is very complicated as shown in FIG. 2D, and the data cannot be fitted by the above formula (1) or the above formula (3). In this case, wavelet transform (the above formula (4)) or short-time Fourier transform (the above formula (6)) is used.

図3は、非常に複雑な振動構造であり上記式(1)でも上記式(3)でもデータをフィットできない場合で、ウェーブレット変換(上記式(4))を用いてプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間を求めた実施例の解析結果を示す図である。図3(A)は、図2のDを拡大したもので時間領域での振動波形を表す図である。この時間波形をウェーブレット変換する訳であるが、ここでマザーウェーブレットとしてガウス関数を変形したGaborのマザーウェーブレット: FIG. 3 shows a very complicated vibration structure in which data cannot be fitted by either the above formula (1) or the above formula (3), and a plasmon-like LO phonon-plasmon using the wavelet transform (the above formula (4)). It is a figure which shows the analysis result of the Example which calculated | required relaxation time of coupling mode (L <+> ). FIG. 3A is an enlarged view of D in FIG. 2 and shows a vibration waveform in the time domain. This time waveform is wavelet transformed, but here Gabor's mother wavelet with a modified Gaussian function as the mother wavelet:

を用いることでMeyerのマザーウェーブレットなど他のマザーウェーブレットを用いる場
合よりもずっと滑らかな[時間, 周波数]平面でのスペクトルを得ることができ、従って、より明確にLOフォノン−プラズモン結合モードとLOフォノンを区別できる。ただし、式(5)でσはガウス関数の半値全幅(Full width at half maximum = FWHM)である。こ
のGaborのマザーウェーブレットを用いて変換した信号の時間−周波数プロットが図3(
B)である。図3(B)のウェーブレット変換スペクトルでは、約26THz付近にピークを持つプラズモンライクなLOフォノン−プラズモン結合モード(L)と約7THz付近にピークを持つLOフォノンとが時間−周波数領域で非常に良く分離観測されていることが分かる。このような分離は、図3(A)のような実時間領域では大変難しい。図3(B)から直ちにプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間τがおおよそ0.2ps(ピコ秒)であることが分かる。より正確にτを求めるためには、図3(B)でプラズモンライクなLOフォノン−プラズモン結合モード(L)のピーク周波数位置の強度の時間変化を抽出すればよい。その結果を図3(C)に示す。図3(C)はウェーブレット変換スペクトル上で、プラズモンライクモード(L)のピーク強度を遅延時間(b)に対してプロットした図である。
Provides a much smoother spectrum in the [time, frequency] plane than with other mother wavelets, such as Meyer's mother wavelet, and therefore more clearly LO phonon-plasmon coupled modes and LO phonons. Can be distinguished. In Equation (5), σ is the full width at half maximum (FWHM) of the Gaussian function. A time-frequency plot of the signal converted using the Gabor mother wavelet is shown in FIG.
B). In the wavelet transform spectrum of FIG. 3B, the plasmon-like LO phonon-plasmon coupled mode (L + ) having a peak near about 26 THz and the LO phonon having a peak near about 7 THz are very much in the time-frequency domain. It can be seen that they are well separated. Such separation is very difficult in the real-time region as shown in FIG. It can be seen from FIG. 3B that the relaxation time τ of the plasmon-like LO phonon-plasmon coupled mode (L + ) is approximately 0.2 ps (picoseconds). In order to obtain τ more accurately, it is only necessary to extract the time change of the intensity at the peak frequency position in the plasmon-like LO phonon-plasmon coupled mode (L + ) in FIG. The result is shown in FIG. FIG. 3C is a diagram in which the peak intensity of the plasmon-like mode (L + ) is plotted against the delay time (b) on the wavelet transform spectrum.

これから直ちにプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間τが約0.2ps(ピコ秒)であることが分かる。従って、上記式(2)から移動度μが約5250cm/Vsと求まる。ただし、ここで電子の電荷e=1.60219×10−19C、電子の有効質量m=0.067m(m=9.10956×10−31kg)を用いた。 From this, it is immediately understood that the relaxation time τ of the plasmon-like LO phonon-plasmon coupled mode (L + ) is about 0.2 ps (picosecond). Therefore, the mobility μ is determined to be about 5250 cm 2 / Vs from the above formula (2). However, here, an electron charge e = 1.60219 × 10 −19 C and an effective mass of electron m * = 0.067 m (m = 9.110956 × 10 −31 kg) were used.

図4は、非常に複雑な振動構造であり上記式(4)でも上記式(6)でもデータをフィットできない場合(図2のDもしくは図3(A)の場合)で、短時間フーリエ変換(上記式(6))を用いてプラズモンライクなLOフォノン−プラズモン結合モード(L)の
緩和時間を求めた実施例の解析結果を示す図である。
FIG. 4 shows a very complicated vibration structure, and when the data cannot be fitted by the above formula (4) or the above formula (6) (in the case of D in FIG. 2 or FIG. 3A), short-time Fourier transform ( It is a figure which shows the analysis result of the Example which calculated | required the relaxation time of plasmon-like LO phonon-plasmon coupling mode (L <+> ) using the said Formula (6).

短時間フーリエ変換(上記式(6))における窓関数はガウス関数とした。何故なら矩形関数など他の窓関数を用いる場合よりもずっと滑らかな短時間スペクトルを得ることができ、従って、より明確に縦波光学(LO)フォノン−プラズモン結合モードとLOフォノンを区別できるからである。すなわち、短時間フーリエ変換(上記式(6))で窓関数をガウス関数にすると、次式   The window function in the short-time Fourier transform (the above formula (6)) was a Gaussian function. This is because a much smoother short-time spectrum can be obtained than when using other window functions such as a rectangular function, and therefore, the longitudinal wave optical (LO) phonon-plasmon coupling mode and LO phonon can be more clearly distinguished. is there. That is, when the window function is made a Gaussian function by the short-time Fourier transform (the above formula (6)),

となる。ただしここで、σはガウス関数の半値全幅(Full width at half maximum = FWHM)である。図4(A)には図2のDの時間領域振動波形データに短時間フーリエ変換を施
した0.036ps毎のスペクトルを示した。約26THz付近にピークを持つプラズモンライクなLOフォノン−プラズモン結合モード(L)と約7THz付近にピークを持つLOフォノンとが周波数領域で非常に良く分離観測されていることが分かる。このような分離は、図3(A)のような実時間領域では大変難しい。図4(B)にはプラズモンライクなLOフォノン−プラズモン結合モード(L)のピーク周波数位置の強度の時間変化を遅延時間bに対してプロットした。これから直ちにプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間τが約0.2ps(ピコ秒)であることが分かる。従って、上記式(2)から移動度μが約5250cm/Vsと求まる。ただし、ここで電子の電荷e=1.60219×10−19C、電子の有効質量m=0.067m(m=9.10956×10−31kg)を用いた。
It becomes. Where σ is the full width at half maximum of the Gaussian function (Full width at half maximum = FWHM). FIG. 4A shows a spectrum every 0.036 ps obtained by performing a short-time Fourier transform on the time-domain vibration waveform data of D of FIG. It can be seen that the plasmon-like LO phonon-plasmon coupled mode (L + ) having a peak in the vicinity of about 26 THz and the LO phonon having a peak in the vicinity of about 7 THz are very well separated and observed in the frequency domain. Such separation is very difficult in the real-time region as shown in FIG. In FIG. 4B, the time change of the intensity at the peak frequency position in the plasmon-like LO phonon-plasmon coupled mode (L + ) is plotted with respect to the delay time b. From this, it is immediately understood that the relaxation time τ of the plasmon-like LO phonon-plasmon coupled mode (L + ) is about 0.2 ps (picosecond). Therefore, the mobility μ is determined to be about 5250 cm 2 / Vs from the above formula (2). However, here, an electron charge e = 1.60219 × 10 −19 C and an effective mass of electron m * = 0.067 m (m = 9.110956 × 10 −31 kg) were used.

最後に図5は、実際の試料n型GaAs単結晶で測定されたプラズモンライクなLOフォノン−プラズモン結合モード(L)での実施例の測定結果である。時間領域振動の解析には短時間フーリエ変換(上記式(5))を用いた。図5(A)は、時間分解反射率変化測定で実際に得られた信号を示す。図5(B)は(A)の時間領域振動波形を短時間フーリエ変換したスペクトルを示す図であり、スペクトル毎の遅延時間bの間隔は可変(0.05〜1.0ps)である。図5(C)は、短時間フーリエ変換したスペクトル図5(B)上で、プラズモンライクモード(L)のピーク強度を遅延時間(b)に対してプロットした図である。図5(C)のプラズモンライクなLOフォノン−プラズモン結合モード(L)の強度時間依存性から、Lモードの緩和時間τがおおよそ0.2ps(ピコ秒)であることが分かる。従って、上記式(2)から移動度μが約5250cm/Vsと求まる。ただし、ここで電子の電荷e=1.60219×10−19C、電子の有効質量m=0.067m(m=9.10956×10−31kg)を用いた。得られた移動度μの値は、室温においてホール測定で得られている実験値μ=8800cm/Vsにほぼ近い値である。 Finally, FIG. 5 shows the measurement results of the example in the plasmon-like LO phonon-plasmon coupling mode (L + ) measured with an actual sample n-type GaAs single crystal. A short-time Fourier transform (the above formula (5)) was used for the analysis of the time domain vibration. FIG. 5A shows a signal actually obtained by the time-resolved reflectance change measurement. FIG. 5B is a diagram showing a spectrum obtained by subjecting the time domain vibration waveform of FIG. 5A to a short-time Fourier transform, and the interval of the delay time b for each spectrum is variable (0.05 to 1.0 ps). FIG. 5C is a diagram in which the peak intensity of the plasmon-like mode (L + ) is plotted against the delay time (b) on the spectrum diagram 5B subjected to the short-time Fourier transform. From the intensity time dependence of the plasmon-like LO phonon-plasmon coupled mode (L + ) in FIG. 5C, it can be seen that the relaxation time τ of the L + mode is approximately 0.2 ps (picoseconds). Therefore, the mobility μ is determined to be about 5250 cm 2 / Vs from the above formula (2). However, here, an electron charge e = 1.60219 × 10 −19 C and an effective mass of electron m * = 0.067 m (m = 9.110956 × 10 −31 kg) were used. The obtained value of the mobility μ is a value substantially close to the experimental value μ = 8800 cm 2 / Vs obtained by the Hall measurement at room temperature.

以上図2〜図5記載の実施例の結果は、プラズモンライクなLOフォノン−プラズモン結合モードとしてLに限定したものではなく、プラズモンライクなLOフォノン−プラズモン結合モードとしてLモードを観測しても同様な結果が得られる。このように、本発明によって、これまで周波数領域で得られたラマンスペクトルを複雑な数式でフィットして求めていた半導体の移動度μが、実時間領域で簡単に求めることが可能であると示された。 Or 2-5 results of Example described the plasmon-like LO phonon - not for limiting the L + as plasmon coupling mode, the plasmon-like LO phonon - as plasmon coupling mode L - mode observed in the Gives similar results. Thus, according to the present invention, it has been shown that the semiconductor mobility μ, which has been obtained by fitting the Raman spectrum obtained in the frequency domain so far with a complicated mathematical formula, can be easily obtained in the real time domain. It was done.

また、本発明の測定対象の固体試料としては、図5の測定結果を得た実施例にあるような半導体単結晶に限るものではなく、金属単結晶や薄膜、半導体超格子や半導体量子井戸
構造、金属・半導体ナノ結晶(量子ドットなど)及びナノワイヤーなどプラズモンを有する全ての固体とすることができ、本発明はこれらの電気的特性(移動度)の評価法として有効である。
Further, the solid sample to be measured of the present invention is not limited to the semiconductor single crystal as in the embodiment where the measurement results of FIG. 5 are obtained, but is a metal single crystal, a thin film, a semiconductor superlattice, or a semiconductor quantum well structure. All solids having plasmons such as metal / semiconductor nanocrystals (such as quantum dots) and nanowires can be used, and the present invention is effective as a method for evaluating these electrical characteristics (mobility).

LOフォノン−プラズモン結合モードの上の分枝(L)と下の分枝(L)の周波数のキャリア密度依存性を示す図である。LO phonon - shows the carrier density dependence of the frequency of - branching on the plasmon coupling mode (L +) and branched lower (L). 本発明による固体のキャリア移動度測定方法の原理の説明図である。It is explanatory drawing of the principle of the solid carrier mobility measuring method by this invention. 本発明の手法のうち、ウェーブレット変換を用いてプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間を求めた実施例の解析結果を示す図であり、(A)は図2のDを拡大したもので時間領域での振動波形を示す図、(B)はウェーブレット変換スペクトルを示す図、(C)はウェーブレット変換スペクトル上で、プラズモンライクモード(L)のピーク強度を遅延時間(b)に対してプロットした図である。It is a figure which shows the analysis result of the Example which calculated | required the relaxation time of the plasmon-like LO phonon-plasmon coupling mode (L <+> ) using the wavelet transformation among the methods of this invention, (A) is D of FIG. (B) is a diagram showing a wavelet transform spectrum, (C) is a wavelet transform spectrum, and the peak intensity of plasmon-like mode (L + ) is expressed as a delay time (L + ). It is the figure plotted with respect to b). 本発明の手法のうち、短時間フーリエ変換を用いてプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間を求めた実施例の解析結果を示す図であり、(A)は図2のDの時間領域振動波形を短時間フーリエ変換したスペクトルを示す図、(B)は短時間フーリエ変換したスペクトル図4(A)上で、プラズモンライクモード(L)のピーク強度を遅延時間(b)に対してプロットした図である。It is a figure which shows the analysis result of the Example which calculated | required the relaxation time of the plasmon-like LO phonon-plasmon coupling mode (L <+> ) using short-time Fourier transform among the methods of this invention, (A) is FIG. The figure which shows the spectrum which carried out the short-time Fourier transform of the time-domain vibration waveform of D of (D), (B) is the spectrum which carried out the short-time Fourier transform, and the peak intensity of plasmon-like mode (L <+> ) is shown as delay time ( It is the figure plotted with respect to b). 本発明の手法のうち、短時間フーリエ変換を用いて、n型GaAsにおいて実験的にプラズモンライクなLOフォノン−プラズモン結合モード(L)の緩和時間を求めた実施例の測定結果を示す図であり、(A)は時間分解反射率測定で得られた信号、(B)は(A)の時間領域振動波形を短時間フーリエ変換したスペクトルを示す図、(C)は短時間フーリエ変換したスペクトルを示す図5(B)上で、プラズモンライクモード(L)のピーク強度を遅延時間(b)に対してプロットした図である。The figure which shows the measurement result of the Example which calculated | required the relaxation time of plasmon-like LO phonon-plasmon coupling mode (L <+> ) experimentally in n-type GaAs using short-time Fourier transform among the methods of this invention. Yes, (A) is a signal obtained by time-resolved reflectance measurement, (B) is a diagram showing a spectrum obtained by short-time Fourier transform of the time-domain vibration waveform of (A), and (C) is a spectrum obtained by short-time Fourier transform. FIG. 5B is a diagram in which the peak intensity of the plasmon-like mode (L + ) is plotted against the delay time (b).

Claims (6)

フォノンとプラズモンが結合しない測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンによる実時間振動構造を測定し、その測定データに対し、下記式(1)を用いて時間領域データのフィッティングを行うことにより、プラズモンの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めることを特徴とする固体のキャリア移動度測定方法。
(上記式中、Aはプラズモンによる振動の振幅、τはプラズモンの緩和時間、tは時間、ωPはプラズモンの振動数、φはプラズモンの振動の初期位相、eは電子の電荷、m
電子の有効質量である)
Using a pump-probe spectroscopy method, the real-time vibration structure due to plasmons in the change in reflectance or transmittance is measured for a solid that is a measurement target in which phonons and plasmons do not bind. 1) is used to obtain the plasmon relaxation time τ by fitting time domain data, and the value of the obtained relaxation time τ is substituted into the following equation (2) to obtain the solid carrier mobility μ. A solid carrier mobility measuring method characterized by the above.
(Where A is the amplitude of plasmon oscillation, τ is the plasmon relaxation time, t is the time, ω P is the plasmon frequency, φ is the initial phase of the plasmon oscillation, e is the charge of the electron, m * is (Effective mass of electrons)
フォノンとプラズモンが結合する測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォノン−プラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(3)を用いて時間領域データのフィッティングを行うことにより、プラズモンライクな縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めることを特徴とする固体のキャリア移動度測定方法。
(上記式中、tは式(1)と同様、AL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードによる振動の振幅、τ(=τL+)は縦波光学フォノン−プラズモ
ン結合モードのプラズモンライクモードの緩和時間、ΩL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードの振動数、φL+は縦波光学フォノン−プラズモン結合モードのプラズモンライクモードの振動の初期位相、Bは縦波光学フォノンによる振動の振幅、τLOは縦波光学フォノンの緩和時間、ΩLOは縦波光学フォノンの振動数、φLOは縦波光学フォノンの振動の初期位相、Cは縦波光学フォノン−プラズモン結合モードのフォノンライクモードによる振動の振幅、τL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの緩和時間、ΩL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの振動数、φL−は縦波光学フォノン−プラズモン結合モードのフォノンライクモードの振動の初期位相、eは電子の電荷、mは電子の有効質量である)
Using pump-probe spectroscopy, we measure the real-time vibration structure of the plasmon-like longitudinal-wave optical phonon-plasmon coupling mode in the reflectance change or transmittance change for the solid that is the object of measurement where phonon and plasmon are combined. Then, the relaxation time τ of the plasmon-like longitudinal wave optical phonon-plasmon coupled mode is obtained by fitting time domain data to the measured data using the following equation (3), and the value of the obtained relaxation time τ The solid carrier mobility μ is obtained by substituting into the following formula (2).
(In the above equation, t is the same as in equation (1), A L + is the amplitude of vibration in the plasmon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, and τ (= τ L + ) is the longitudinal wave optical phonon-plasmon coupled mode. Plasmon-like mode relaxation time, Ω L + is the frequency of plasmon-like mode in longitudinal optical phonon-plasmon coupled mode, φ L + is the initial phase of plasmon-like mode vibration in longitudinal wave optical phonon-plasmon coupled mode, and B is longitudinal Amplitude of vibration by wave optical phonon, τ LO is relaxation time of longitudinal optical phonon, Ω LO is frequency of longitudinal optical phonon, φ LO is initial phase of longitudinal optical phonon vibration, C is longitudinal optical phonon − The amplitude of vibration due to the phonon-like mode of the plasmon coupled mode, τ L− is the phonon line of the longitudinal optical phonon-plasmon coupled mode. Mode relaxation time, Ω L− is the frequency of the phonon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, φ L− is the initial phase of the phonon-like mode of the longitudinal wave optical phonon-plasmon coupled mode, and e is the electron Charge, m * is the effective mass of electrons)
フォノンとプラズモンが結合する測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォノン−プ
ラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(4)を用いてウェーブレット変換を行うことにより、プラズモンライクな縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めることを特徴とする固体のキャリア移動度測定方法。
(上記式中、Ψ(x)はマザーウェーブレットと称される、波束を表す関数、xは入力信号で時間の関数、bはこの波束の時間遅延を表す変数、aはこの波束の時間軸に対する伸縮度、tは時間、f(t)は被変換関数である)
Using pump-probe spectroscopy, we measure the real-time vibration structure of the plasmon-like longitudinal-wave optical phonon-plasmon coupling mode in the reflectance change or transmittance change for the solid that is the object of measurement where phonon and plasmon are combined. The measurement data is subjected to wavelet transform using the following equation (4) to obtain the relaxation time τ of the plasmon-like longitudinal wave optical phonon-plasmon coupled mode, and the value of the obtained relaxation time τ is expressed by the following equation: A solid carrier mobility measuring method, wherein the solid carrier mobility μ is obtained by substituting in (2).
(In the above equation, Ψ (x) is a function that represents a wave packet called a mother wavelet, x is a function of time as an input signal, b is a variable that represents a time delay of the wave packet, and a is a time axis of the wave packet. Degree of elasticity, t is time, f (t) is the transformed function)
ウェーブレット変換のマザーウェーブレットとして下記式(5)で示す、ガウス関数を変形したGaborのマザーウェーブレットを用いることを特徴とする請求項3に記載の固体
のキャリア移動度測定方法。
(上記式中、σはガウス関数の半値全幅、xは入力信号で時間の関数である)
を用いることを特徴とする固体のキャリア移動度μを求める方法。
4. The solid carrier mobility measuring method according to claim 3, wherein a Gabor mother wavelet obtained by transforming a Gaussian function represented by the following equation (5) is used as a mother wavelet for wavelet transformation.
(In the above formula, σ is the full width at half maximum of the Gaussian function, x is the input signal and a function of time)
A method for obtaining a carrier mobility μ of a solid, characterized in that
フォノンとプラズモンが結合する測定対象である固体に対し、ポンプ−プローブ分光法を用いて、反射率変化又は透過率変化におけるプラズモンライクな縦波光学フォノン−プラズモン結合モードによる実時間振動構造を測定し、その測定データに対し、下記式(6)を用いて短時間フーリエ変換を行うことにより、縦波光学フォノン−プラズモン結合モードの緩和時間τを求め、求めた緩和時間τの値を下記式(2)に代入することにより固体のキャリア移動度μを求めることを特徴とする固体のキャリア移動度測定方法。
(上記式中、Ψ(x)は窓関数、xは入力信号で時間の関数、bは時間遅延、ωは周波数、tは時間、f(t)は被変換関数である)
Using pump-probe spectroscopy, we measure the real-time vibration structure of the plasmon-like longitudinal-wave optical phonon-plasmon coupling mode in the reflectance change or transmittance change for the solid that is the object of measurement where phonon and plasmon are combined. The measurement data is subjected to short-time Fourier transform using the following equation (6) to obtain the relaxation time τ of the longitudinal optical phonon-plasmon coupled mode, and the value of the obtained relaxation time τ is represented by the following equation ( A solid carrier mobility measuring method, wherein the solid carrier mobility μ is obtained by substituting in 2).
(In the above formula, Ψ (x) is a window function, x is an input signal and a function of time, b is a time delay, ω is a frequency, t is time, and f (t) is a transformed function)
短時間フーリエ変換の窓関数としてガウス関数を用い、下記式(7)により短時間フーリエ変換を行うことを特徴とする請求項5に記載の固体のキャリア移動度測定方法。
(上記式中、σはガウス関数の半値全幅である)
6. The solid carrier mobility measuring method according to claim 5, wherein a Gaussian function is used as a window function of the short-time Fourier transform, and the short-time Fourier transform is performed by the following equation (7).
(In the above formula, σ is the full width at half maximum of the Gaussian function)
JP2006172865A 2006-06-22 2006-06-22 Method for measuring carrier mobility of individuals Expired - Fee Related JP4831482B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006172865A JP4831482B2 (en) 2006-06-22 2006-06-22 Method for measuring carrier mobility of individuals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006172865A JP4831482B2 (en) 2006-06-22 2006-06-22 Method for measuring carrier mobility of individuals

Publications (2)

Publication Number Publication Date
JP2008002960A true JP2008002960A (en) 2008-01-10
JP4831482B2 JP4831482B2 (en) 2011-12-07

Family

ID=39007467

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006172865A Expired - Fee Related JP4831482B2 (en) 2006-06-22 2006-06-22 Method for measuring carrier mobility of individuals

Country Status (1)

Country Link
JP (1) JP4831482B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313849A (en) * 2010-06-30 2012-01-11 株式会社神户制钢所 The evaluation method of oxide semiconductor thin-film and the method for quality control of oxide semiconductor thin-film
KR101667150B1 (en) * 2015-05-10 2016-10-17 고려대학교 산학협력단 Surface plasmon resonance curve fitting method and analyzing method of target substance using spr sensor
WO2018021259A1 (en) * 2016-07-23 2018-02-01 国立大学法人千葉大学 Infrared light element

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1060510A (en) * 1996-08-16 1998-03-03 Nkk Corp Method for predicting abnormal furnace condition of blast furnace
JP2002214137A (en) * 2001-01-12 2002-07-31 National Institute For Materials Science Method and device for optically measuring characteristic of metal sample
JP2002296251A (en) * 2001-03-30 2002-10-09 Railway Technical Res Inst Technique of evaluating concrete by impact sound
JP2003043021A (en) * 2001-07-31 2003-02-13 Railway Technical Res Inst Concrete hammering inspection method and apparatus therefor
JP2003506675A (en) * 1999-07-30 2003-02-18 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Improved optical method for measuring metal thin films
JP2004183140A (en) * 2002-12-03 2004-07-02 Murata Mach Ltd Apparatus and method for analyzing yarn tension for textile machine
JP2004349582A (en) * 2003-05-23 2004-12-09 National Institute For Materials Science Method and device of measuring defect in solid
JP2005092147A (en) * 2003-09-19 2005-04-07 Japan Science & Technology Agency Far infrared optical switch based on coherent control of organic charge transfer complex
JP2005235831A (en) * 2004-02-17 2005-09-02 Toyota Central Res & Dev Lab Inc P-type semiconductor and semiconductor device
JP2005233683A (en) * 2004-02-17 2005-09-02 Aisin Seiki Co Ltd Multichannel terahertz wave spectrum measuring method and measuring device

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1060510A (en) * 1996-08-16 1998-03-03 Nkk Corp Method for predicting abnormal furnace condition of blast furnace
JP2003506675A (en) * 1999-07-30 2003-02-18 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Improved optical method for measuring metal thin films
JP2002214137A (en) * 2001-01-12 2002-07-31 National Institute For Materials Science Method and device for optically measuring characteristic of metal sample
JP2002296251A (en) * 2001-03-30 2002-10-09 Railway Technical Res Inst Technique of evaluating concrete by impact sound
JP2003043021A (en) * 2001-07-31 2003-02-13 Railway Technical Res Inst Concrete hammering inspection method and apparatus therefor
JP2004183140A (en) * 2002-12-03 2004-07-02 Murata Mach Ltd Apparatus and method for analyzing yarn tension for textile machine
JP2004349582A (en) * 2003-05-23 2004-12-09 National Institute For Materials Science Method and device of measuring defect in solid
JP2005092147A (en) * 2003-09-19 2005-04-07 Japan Science & Technology Agency Far infrared optical switch based on coherent control of organic charge transfer complex
JP2005235831A (en) * 2004-02-17 2005-09-02 Toyota Central Res & Dev Lab Inc P-type semiconductor and semiconductor device
JP2005233683A (en) * 2004-02-17 2005-09-02 Aisin Seiki Co Ltd Multichannel terahertz wave spectrum measuring method and measuring device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313849A (en) * 2010-06-30 2012-01-11 株式会社神户制钢所 The evaluation method of oxide semiconductor thin-film and the method for quality control of oxide semiconductor thin-film
KR101667150B1 (en) * 2015-05-10 2016-10-17 고려대학교 산학협력단 Surface plasmon resonance curve fitting method and analyzing method of target substance using spr sensor
WO2018021259A1 (en) * 2016-07-23 2018-02-01 国立大学法人千葉大学 Infrared light element
JPWO2018021259A1 (en) * 2016-07-23 2019-05-16 国立大学法人千葉大学 Infrared light element

Also Published As

Publication number Publication date
JP4831482B2 (en) 2011-12-07

Similar Documents

Publication Publication Date Title
Huang et al. Surface optical rectification from layered MoS2 crystal by THz time-domain surface emission spectroscopy
Ge et al. Coherent longitudinal acoustic phonon approaching THz frequency in multilayer molybdenum disulphide
Suess et al. Power-dependent Raman analysis of highly strained Si nanobridges
Bensmann et al. Near-field imaging and spectroscopy of locally strained GaN using an IR broadband laser
Mante et al. THz acoustic phonon spectroscopy and nanoscopy by using piezoelectric semiconductor heterostructures
WO2016007950A1 (en) Apparatus and methods for probing a material as a function of depth using depth-dependent second harmonic generation
Lanzillotti-Kimura et al. Towards GHz–THz cavity optomechanics in DBR-based semiconductor resonators
Cunningham et al. On-chip terahertz systems for spectroscopy and imaging
Pashnev et al. Terahertz time-domain spectroscopy of two-dimensional plasmons in AlGaN/GaN heterostructures
US11635376B2 (en) Fiber-optic based material property measurement system and related methods
Nakanishi et al. Comparison between laser terahertz emission microscope and conventional methods for analysis of polycrystalline silicon solar cell
Beardsley et al. Nanomechanical probing of the layer/substrate interface of an exfoliated InSe sheet on sapphire
JP4831482B2 (en) Method for measuring carrier mobility of individuals
Cinquanta et al. Charge carrier dynamics in 2D materials probed by ultrafast THzspectroscopy
Gorbatova et al. Terahertz emission from a monolayer tungsten diselenide surface
Yang et al. Generation and coherent control of terahertz acoustic phonons in superlattices of perovskite oxides
Molis et al. Intervalley separation in the conduction band of InGaAs measured by terahertz excitation spectroscopy
Poyser et al. Coherent phonon optics in a chip with an electrically controlled active device
Chatzakis et al. Broadband terahertz modulation in electrostatically-doped artificial trilayer graphene
Ahn et al. Dynamic near-field optical interaction between oscillating nanomechanical structures
Armstrong Visualizing Carrier Dynamics in Transition Metal Dichalcogenide Nanoflakes using Femtosecond Pump-Probe Microscopy
Yang et al. Non-destructive plasma frequency measurement for a semiconductor thin film using broadband surface plasmon polaritons
Beleckaitė et al. Determination of the terahertz pulse emitting dipole orientation by terahertz emission measurements
Li et al. Characterization of ultrathin films by laser-induced sub-picosecond photoacoustics with coherent extreme ultraviolet detection
Van Goethem Imaging Charge Carrier and Acoustic Phonon Dynamics in Semiconductor Nanomaterials using Ultrafast Pump-Probe Microscopy

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090611

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110601

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110614

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110805

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: 20110906

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110908

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140930

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140930

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees