JP2003014623A - Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation - Google Patents

Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation

Info

Publication number
JP2003014623A
JP2003014623A JP2001202984A JP2001202984A JP2003014623A JP 2003014623 A JP2003014623 A JP 2003014623A JP 2001202984 A JP2001202984 A JP 2001202984A JP 2001202984 A JP2001202984 A JP 2001202984A JP 2003014623 A JP2003014623 A JP 2003014623A
Authority
JP
Japan
Prior art keywords
surface plasmon
plasmon resonance
spr
curve
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2001202984A
Other languages
Japanese (ja)
Inventor
Koji Suzuki
鈴木  孝治
Kazuyoshi Kurihara
一嘉 栗原
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.)
Kanagawa Academy of Science and Technology
Japan Science and Technology Agency
Original Assignee
Kanagawa Academy of Science and Technology
Japan Science and Technology Corp
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 Kanagawa Academy of Science and Technology, Japan Science and Technology Corp filed Critical Kanagawa Academy of Science and Technology
Priority to JP2001202984A priority Critical patent/JP2003014623A/en
Publication of JP2003014623A publication Critical patent/JP2003014623A/en
Pending legal-status Critical Current

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a method by which an asymmetric SPR curve usable for measuring the SPR of gold more appropriately than the equation of the Kretchmann's SPR curve can be decided. SOLUTION: This method is constituted to accurately decide a surface plasmon resonance angle based on the fitting equation of an asymmetric surface plasmon resonance curve derived from the three-layer Fresnel equation for P-polarization at the time of measuring the state of a material, by utilizing a surface plasmon resonance phenomenon by using a surface plasmon resonance sensor system composed of three layers of a prism, a thin metallic film, and a sensing layer.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【発明の属する技術分野】本発明は、表面プラズモン共
鳴現象を利用したセンシングにおける表面プラズモン共
鳴カーブの非対称表面プラズモン共鳴カーブ方程式によ
る決定方法に関するものである。
TECHNICAL FIELD The present invention relates to a method for determining a surface plasmon resonance curve in sensing using the surface plasmon resonance phenomenon by an asymmetric surface plasmon resonance curve equation.

【0002】[0002]

【従来の技術】一般に知られているように、表面プラズ
モン共鳴(SPR)センサーは、特異的な吸着や金薄膜上
に作製したホストのリカンドとゲストの目的検出物質と
の複合体によって生じる屈折率変化の測定に基づいた有
用な分析手段である。最もよく知られているSPRセンサ
ーとしては、抗原−抗体相互作用に関連したタンパク質
のような特異的相互作用を有する高分子を検出するため
に用いられるセンサーがある。SPR測定において、反射
率と入射角の関係を示すSPRカーブは、屈折率変化によ
り起こる信号を解析する基本的な概念である。理論的SP
Rカーブの式を最初に導いたのはクレッチマンであり、
そのカーブは共鳴角と通常呼ばれている角度を最小反射
率の所に有している。クレッチマンによる理論式によ
り、SPRの信号は、共鳴角の変化を観測することによっ
て得られている。しかし、クレッチマンの理論式は、銀
を利用したSPRにより通常得られる対称のSPRカーブにの
み適用でき、金を利用したSPRによる非対称のSPRカーブ
には適用できない。金が銀よりも酸化において有利な化
学的性質を持っていることから、金のSPRが銀のSPRより
も実際に広く用いられているにもかかわらず、本件発明
者らの知る限り、金のSPRに最適な理論や概念は未だに
十分に検討されていない。
BACKGROUND ART As is generally known, a surface plasmon resonance (SPR) sensor is a refractive index generated by specific adsorption or a complex of a host licand formed on a gold thin film and a guest target detection substance. It is a useful analytical tool based on measuring changes. The most well-known SPR sensors are those used to detect macromolecules with specific interactions such as proteins associated with antigen-antibody interactions. In SPR measurement, the SPR curve showing the relationship between the reflectance and the incident angle is a basic concept for analyzing a signal caused by a change in the refractive index. Theoretical SP
Kretschmann was the first to lead the R curve formula,
The curve has what is commonly called the resonance angle at its minimum reflectance. According to the Kretschmann theoretical formula, the SPR signal is obtained by observing the change in the resonance angle. However, Kretschmann's theoretical formula can be applied only to the symmetric SPR curve normally obtained by SPR using silver, and cannot be applied to the asymmetric SPR curve of SPR using gold. To the best knowledge of the inventors, gold SPR is actually more widely used than silver SPR because gold has more favorable chemistry in oxidation than silver, despite the fact that gold SPR is widely used. The optimal theory and concept for SPR has not been fully examined yet.

【0003】従来、金のSPRカーブを実際にデータ解析
する場合、共鳴角(RA)を見つける単純な方法として、
二次関数によるフィッティング、高次関数によるフィッ
ティング、重心による方法、スムージングをして1階微
分で零となる点を求める方法などが適用されている。さ
らに進んだ方法として部分的に比重をかけてフィッティ
ングする方法や、最適線形法などが提案されている。
Conventionally, when actually analyzing data of a gold SPR curve, as a simple method for finding a resonance angle (RA),
The fitting by a quadratic function, the fitting by a higher-order function, the method by the center of gravity, the method by smoothing, and the method by which the point which becomes zero in the first derivative is obtained are applied. Further advanced methods have been proposed, such as a method of partially fitting specific gravity and an optimal linear method.

【0004】[0004]

【発明が解決しようとする課題】しかしながら、残念な
ことに、上記で述べた従来の方法はいずれも、SPRの理
論に基づいたものでなく、単に最小値を求める範中の方
法であって、最小値から離れた有効で正確なSPRのデー
タは、使われていない。そこで、発明者等はSPR理論に
基づく適切なフィッティング関数を見出して、金のSPR
センシングをもっと高感度にすることができた。従っ
て、本発明は、クレッチマンのSPRカーブの式より一層
適当な、金のSPR測定に使える非対称SPRカーブを決定す
る方法を提供することを目的としている。
However, unfortunately, none of the conventional methods described above is based on the theory of SPR, and is merely a method for obtaining the minimum value. Valid and accurate SPR data away from the minimum are not used. Therefore, the inventors found an appropriate fitting function based on the SPR theory and found that the gold SPR
We were able to make sensing more sensitive. It is therefore an object of the present invention to provide a method for determining an asymmetric SPR curve that can be used for gold SPR measurement, which is more suitable than the Kretschmann SPR curve equation.

【0005】[0005]

【課題を解決するための手段】上記の目的を達成するた
めに、本発明によれば、プリズム−金属薄膜−センシン
グ層の三層から成る表面プラズモン共鳴センサーシステ
ムを用い、表面プラズモン共鳴現象を利用したセンシン
グにおける表面プラズモン共鳴カーブの非対称表面プラ
ズモン共鳴カーブ方程式による決定方法において、A、
B、C、D、Eの5つの実数のうち少なくとも1つ以上を実
数パラメータとし、sを非対称表面プラズモン共鳴カー
ブの特性因子とし、Rを表面プラズモン共鳴センサーシ
ステムにおける光の反射率、rpmsをプリズム−金属薄
膜−センシング層の境界における振幅反射率とし、kx
を境界面に平行な入射光の波数ベクトルとし、koをプ
リズムの存在しない時の金属薄膜とセンシング層との界
面に理論的に生成される表面プラズモンの複素波数と
し、kRをプリズムの存在する時のkoの補正項とすると
き、P偏光の三層フレネルの式から導いた式 で表される非対称表面プラズモン共鳴カーブのフィッテ
ィング方程式に基き表面プラズモン共鳴角を決定するこ
とを特徴としている。実験から得られた表面プラズモン
共鳴カーブのフィッティングにおいて、反射率データを
R、入射角データをXとするとき、上記フィッティング
方程式は で表わすことができ、このフィッティング方程式に基き
表面プラズモン共鳴角は決定され得る。また、最小反射
率に対する角度は で定義され、これを表面プラズモン共鳴角とすることが
できるパラメータDは修正共鳴角として定義され得る。
さらに、観測した表面プラズモン共鳴カーブはコンピュ
ータ上でソフトウエア的に処理され、共鳴角を求めるよ
うにされ得る。
In order to achieve the above object, according to the present invention, a surface plasmon resonance sensor system comprising three layers of a prism, a thin metal film and a sensing layer is used, and the surface plasmon resonance phenomenon is utilized. In the method of determining the surface plasmon resonance curve by the asymmetric surface plasmon resonance curve equation in
At least one of the five real numbers B, C, D, and E is a real parameter, s is a characteristic factor of the asymmetric surface plasmon resonance curve, R is the reflectance of light in the surface plasmon resonance sensor system, and r pms is The amplitude reflectance at the boundary between the prism, the metal thin film, and the sensing layer is represented by k x
Is the wave vector of the incident light parallel to the boundary surface, k o is the complex wave number of the surface plasmon theoretically generated at the interface between the metal thin film and the sensing layer when there is no prism, and k R is the existence of the prism. When the correction term for k o is set , the formula derived from the P-polarized three-layer Fresnel formula It is characterized in that the surface plasmon resonance angle is determined based on the fitting equation of the asymmetric surface plasmon resonance curve represented by. In the fitting of the surface plasmon resonance curve obtained from the experiment, when the reflectance data is R and the incident angle data is X, the above fitting equation is The surface plasmon resonance angle can be determined based on this fitting equation. Also, the angle for the minimum reflectance is And the parameter D, which can be the surface plasmon resonance angle, can be defined as the modified resonance angle.
Further, the observed surface plasmon resonance curve can be processed by software on a computer to determine the resonance angle.

【0006】上記のフィッティング方程式は、R(x)=1-
α/{(X-β)22}(ここで、Rは入射光の反射率、α、
β、γは実数パラメータ、又は入射角である。)で与え
られるクレッチマンによる式よりも適切なSPRカーブ方
程式として利用できる。また本発明における上記の方程
式は金を利用したSPRのSPRカーブをフィッティングする
関数として極めて有効である。
The above fitting equation is R (x) = 1-
α / {(X-β) 2 + γ 2 } (where R is the reflectance of the incident light, α,
β and γ are real parameters or incident angles. ) Can be used as a more appropriate SPR curve equation than the Kretschmann equation given in. The above equation in the present invention is extremely effective as a function for fitting the SPR curve of SPR using gold.

【0007】[0007]

【発明の実施の形態】以下添付図面を参照して本発明の
実施の形態について説明する。図1には、クレッチマン
減衰全反射(ATR)形態のプリズム−金属−センシン
グ層構造の三層から成るSPRセンサーシステムの簡単
な例を示しており、ガラスプリズムp、金属薄膜m及び
センシング層sから成ると仮定する。なお図1において
θは入射角、dは金属薄膜mの厚さである。P偏光の光
が入射角θ(単位ラジアン)で入射する時、光の反射率
Rは、次のようにP偏光三層フレネルの方程式で与えら
れる。 上記式において、γpmは、ガラスプリズムp−金属薄膜
mの境界における振幅反射率であり、そしてγpm及びγ
msは、P偏光フレネルの式で与えられる。εjは媒質jの
誘電率である。kjzは境界面に垂直な成分を持つ媒質jの
波数ベクトルである。kxは、境界面に平行な入射光の波
数ベクトル成分であり、実数として与えられる。dは金
属薄膜mの厚さである。ωは入射光の角周波数であり、
λは波長であり、またcは光速である。
BEST MODE FOR CARRYING OUT THE INVENTION Embodiments of the present invention will be described below with reference to the accompanying drawings. FIG. 1 shows a simple example of a SPR sensor system consisting of three layers of prism-metal-sensing layer structure of Kretschmann Attenuated Total Reflection (ATR) type, which includes a glass prism p, a metal thin film m and a sensing layer s. Suppose that In FIG. 1, θ is the incident angle and d is the thickness of the metal thin film m. When P-polarized light is incident at an incident angle θ (unit radian), the reflectance R of the light is given by the P-polarized three-layer Fresnel equation as follows. In the above equation, γ pm is the amplitude reflectance at the boundary of the glass prism p-metal thin film m, and γ pm and γ pm
ms is given by the P-polarized Fresnel equation. ε j is the permittivity of medium j. k jz is the wave number vector of medium j that has a component perpendicular to the boundary surface. k x is a wave vector component of incident light parallel to the boundary surface, and is given as a real number. d is the thickness of the metal thin film m. ω is the angular frequency of the incident light,
λ is the wavelength and c is the speed of light.

【0008】プリズムpが無い場合の金属薄膜mとセン
シング層sとの境界に発生する表面プラズモンの複素波
数koは、kxを未知変数とする次の方程式より導かれる。
The complex wave number k o of the surface plasmon generated at the boundary between the metal thin film m and the sensing layer s without the prism p is derived from the following equation with k x as an unknown variable.

【0009】最も単純な近似として、変数kxがkoに等し
い時にSPR現象が強く起きると仮定することができる。k
x=koの時、式(3)のγmsの値が(無限大となり)特異
的になるので、式(2)は次のように変形できる。
As the simplest approximation, it can be assumed that the SPR phenomenon occurs strongly when the variable k x equals k o . k
When x = k o , the value of γ ms in the equation (3) becomes peculiar (becomes infinite), so the equation (2) can be modified as follows.

【0010】式(13)から(16)までの4つの式に
より式(10)は次のように簡単に表すことができる。 ここで、kRはガラスプリズムpが存在する時のk0の補正
項である。クレッチマン配置における表面プラズモン
は、次式で与えられる複素波数kSPを持つ。 ksp=k0+kR (21) 式(17)は、SPRカーブを議論する時に基本となる式
である。クレッチマンによる対称SPRカーブの式を導く
近似では、2つの条件が仮定されている。つまり、 ここで、φは実数である。式(22)と式(23)の2
つの条件では、項ρσ-1は次のように表すことができ
る。 ここで、χ'とχ''の数学的記述は、複素数χの実部と
虚部をそれぞれ示す。式(25)のSPRカーブは、(k
x-k'SP)の2次の項しか含まないので、kx=k'SPの軸に対
称である。
The equation (10) can be simply expressed as follows by the four equations (13) to (16). Here, k R is a correction term of k 0 when the glass prism p exists. The surface plasmon in the Kretschmann configuration has a complex wave number k SP given by the following equation. k sp = k 0 + k R (21) Expression (17) is a basic expression when discussing the SPR curve. Two conditions are assumed in the approximation by Kretschmann to derive the symmetric SPR curve equation. That is, Here, φ is a real number. 2 of Expression (22) and Expression (23)
Under one condition, the term ρσ −1 can be expressed as Here, the mathematical description of χ ′ and χ ″ indicates the real and imaginary parts of the complex number χ, respectively. The SPR curve of equation (25) is (k
It is symmetric about the axis of k x = k ' SP because it contains only the second-order term of x -k' SP ).

【0011】本発明では、式(17)に対して何の条件
も課さずに、5つのパラメータA、B、C、D、Eを使っ
て、式(17)を次のように厳密に変形できる。 ここで、値sは、非対称SPRカーブの特性因子として定義
される。式(29)におけるパラメータのSPRカーブに
対する意味は次の通りである。AはSPRカーブの最大値を
与える。B、CはそれぞれSPRカーブの対称性及び非対称
性に寄与する。Dは、後で説明するように、屈折率変化
を観測する最も適した変数である修正共鳴角である。E
はクレッチマン配置における表面プラズモンの複素波数
の虚部であり、SPRカーブの深と幅に関係している。式
(25)のクレッチマンの式と比較すると、本発明にお
ける式(29)は、kx-k'SPの一次の項を含んでいるた
め、kx=k'SP軸に対称でない。式(29)は、基本式で
ある式(17)に何の条件も課していないので、クレッ
チマンの式よりも一般的で適切である。その結果、フィ
ッティング関数としてより適当である。
In the present invention, the expression (17) is strictly modified as follows using the five parameters A, B, C, D and E without imposing any condition on the expression (17). it can. Here, the value s is defined as a characteristic factor of the asymmetric SPR curve. The meanings of the parameters in equation (29) for the SPR curve are as follows. A gives the maximum of the SPR curve. B and C contribute to the symmetry and asymmetry of the SPR curve, respectively. D is the modified resonance angle, which is the most suitable variable for observing the refractive index change, as will be explained later. E
Is the imaginary part of the complex wave number of the surface plasmon in the Kretschmann arrangement, and is related to the depth and width of the SPR curve. When compared with the Kretschmann's equation of Equation (25), Equation (29) in the present invention includes a first-order term of k x -k ' SP , and is therefore not symmetric about the k x = k' SP axis. The expression (29) is more general and appropriate than the Kretschmann's expression because it does not impose any condition on the basic expression (17). As a result, it is more suitable as a fitting function.

【0012】式(29)において、反射率は従来の共鳴
角θminと最小反射率Rminとして次式で表される。 して示された近似はB>>CE(>0)の条件で得られる。
In the equation (29), the reflectance is represented by the following equation as the conventional resonance angle θ min and the minimum reflectance R min . The approximation shown by is obtained under the condition of B >> CE (> 0).

【0013】SPR分析の原理は、エバネッセントと表面
プラズモンの実部における運動量保存則に基づいてお
り、次式で表される。 ここで、kxは入射光により生成されたエバネッセント波
の波数であり、k'SPは表面プラズモン波の複素波数の実
部である。式(39)は次式で与えられる特別な角度θ
SPRを持つ。 従って、式(40)の近似式である下記の有効な式を基
本として、センシング層sの屈折率変化を測定できる。 上記式においてnSはセンシング層sの屈折率である。
式(25)のクレッチマンの近似式では、角度θ
SPRは、式(25)により反射率が最小になるので、従
来の共鳴角である。しかし、式(29)の非対称方程式
では、最小反射率を示す従来の共鳴角が式(36)で与
えられるので、角度θSPRは従来の共鳴角ではない。そ
こで、本明細書では最小反射率を与える従来の共鳴角と
区別してθSPRは修正共鳴角と呼ぶことにする。物理的
現象では、θSPR及びθminは、エバネッセント波と表面
プラズモンの運動量保存則がそれぞれ実部と複素数で成
り立つことによって決められる。屈折率変化を線形に決
定する際において、修正共鳴角すなわちθSPRは従来の
共鳴角よりも優れている。その理由は式(41)で与え
られる有効な式がθminよりもθSPRの良い近似であるか
らである。このことは、長い間、最小反射率を与える共
鳴角こそがセンシング層sの屈折率変化を観測する上で
最も適していると考えられてきたので、新しい概念であ
る。本発明の理論は、従来の共鳴角によって屈折率の線
形な変化を正確に決定するのに理論的限界があることを
示すものである。
The principle of SPR analysis is based on the law of conservation of momentum in the real part of evanescent and surface plasmons, and is expressed by the following equation. Here, k x is the wave number of the evanescent wave generated by the incident light, and k ′ SP is the real part of the complex wave number of the surface plasmon wave. Equation (39) is a special angle θ given by
Have an SPR . Therefore, the change in the refractive index of the sensing layer s can be measured based on the following effective formula, which is an approximate formula of the formula (40). In the above equation, n S is the refractive index of the sensing layer s.
In the Kretschmann approximation formula (25), the angle θ
The SPR is the conventional resonance angle because the reflectance is minimized by the equation (25). However, in the asymmetric equation of equation (29), the conventional resonance angle showing the minimum reflectance is given by equation (36), and therefore the angle θ SPR is not the conventional resonance angle. Therefore, in this specification, θ SPR will be referred to as a corrected resonance angle in distinction from the conventional resonance angle that gives the minimum reflectance. In a physical phenomenon, θ SPR and θ min are determined by the realization of the momentum conservation law of the evanescent wave and the surface plasmon and the complex number, respectively. The modified resonance angle, θ SPR, is superior to conventional resonance angles in determining the refractive index change linearly. The reason is that the effective equation given by equation (41) is a better approximation of θ SPR than θ min . This is a new concept because it has long been considered that the resonance angle that gives the minimum reflectance is the most suitable for observing the refractive index change of the sensing layer s. The theory of the present invention shows that there is a theoretical limit to accurately determining the linear change in the refractive index by the conventional resonance angle.

【0014】本発明における非対称SPRカーブ方程式の
数値的考察 数値計算は、マセマティカ4(Wolf共鳴角m Research社
製)及びイゴールプル(WaveMetrics社製)の2つのソ
フトウェアを使って行った。マセマティカを用いての計
算によって得られた理論的SPRカーブの数値データを、
イゴールプロのカーブフィッティングルーチンを使っ
て、本発明における非対称SPRカーブ方程式に合わせ
た。イゴールプルを使用した場合のカーブフィッティン
グは、フィッティング関数や、各点の重み、各係数に関
する偏微分を計算するためのイプシロン値をユーザーが
設定でいるように設計されている。イゴールプロは、非
線形最小二乗法で係数値を見つける問題に対して標準的
解法となっているレベンバーグ・マーカート法を使って
いる。本発明で行ったすべてのフィッティングに関し
て、重みとイプシロン値は自動選択とした。イゴールプ
ロの取扱説明書によると、重みは全ての生データに対し
てであり、各係数Kに対するイプシロン値は、10- 10×
Kである。
Numerical Consideration of Asymmetric SPR Curve Equation in the Present Invention Numerical calculation was carried out using two softwares, Mathematica 4 (manufactured by Wolf resonance angle m Research) and Igor Pull (manufactured by WaveMetrics). Numerical data of theoretical SPR curve obtained by calculation using Mathematica
Igor Pro's curve fitting routine was used to fit the asymmetric SPR curve equation in the present invention. The curve fitting using Igorpur is designed so that the user can set the fitting function, the weight of each point, and the epsilon value for calculating the partial derivative of each coefficient. Igor Pro uses the Levenberg-Marquardt method, which is the standard solution for the problem of finding coefficient values by the nonlinear least-squares method. Weights and epsilon values were automatically selected for all fittings performed in the present invention. According to Igor Pro's instructions weights be to all raw data, epsilon values for each coefficient K, 10 - 10 ×
K.

【0015】[0015]

【実施例】以下、0、1、4、10、20%の蔗糖水溶
液の屈折率を測定した実験例について説明する。蔗糖
は、和光純薬より購入した。実験は、市販のSPR装置(D
KK社製)を使って行った。この装置では、プリズムや光
源のような光学素子は自由に変えることができ、光学系
は典型的なクレッチマンATR配置である。プリズムは、
587nmの時1.80の高い屈折率を持つSFL6ガラス
(OHA共鳴角社製)で作製した半円柱プリズムであ
る。光源としては半値幅20nm中心波長630nmの発光
ダイオード(LED)と、半値幅40nm中心波長88nmの
発光ダイオードとの2つを使用した。光は円柱(シリン
ドリカル)レンズを介して集束ビームに集光され、それ
によりセンサー表面で入射角の範囲を同時に観察できる
ようにされる。電荷結合素子(CCD)(SONY社製のXC
−77)を用いて入射角度の関数として反射光強度を測
定した。直線偏光子をプリズムとCCDの間に入れて、p又
はsの偏光を持つ反射光を選択的に観測した。SPRカーブ
は、光強度の不均一分布を除くため、p偏光の反射光強
度をs偏光の反射光強度で割ることによって得た。SPR
実験に用いた金基板は、51×20mm、厚さ1mmのSFL6
ガラス基板上に20×20mmの範囲で45nm厚の金薄膜
を形成したものである。金とガラスとの固着を良くする
ために、金薄膜とガラス基板の間に厚さ5nmのクロム層
を設けた。金基板を半円柱プリズムへ固定する時ガラス
表面の反射率を減裏するために、マッチングオイルを使
用した。
EXAMPLES Hereinafter, experimental examples in which the refractive index of 0, 1, 4, 10, 20% sucrose aqueous solution was measured will be described. Sucrose was purchased from Wako Pure Chemical Industries. The experiment was conducted using a commercially available SPR device (D
KK). In this device, optical elements such as prisms and light sources can be freely changed, and the optical system is a typical Kretschmann ATR arrangement. The prism is
It is a semi-cylindrical prism made of SFL6 glass (manufactured by OHA Resonance Co.) having a high refractive index of 1.80 at 587 nm. As the light source, a light emitting diode (LED) having a half width of 20 nm and a center wavelength of 630 nm and a light emitting diode having a half width of 40 nm and a center wavelength of 88 nm were used. The light is focused through a cylindrical lens into a focused beam, which makes it possible to simultaneously observe a range of angles of incidence on the sensor surface. Charge Coupled Device (CCD) (XC made by SONY
-77) was used to measure the reflected light intensity as a function of incident angle. A linear polarizer was inserted between the prism and CCD to selectively observe the reflected light with p or s polarization. The SPR curve was obtained by dividing the reflected light intensity of p-polarized light by the reflected light intensity of s-polarized light in order to eliminate the uneven distribution of light intensity. SPR
The gold substrate used in the experiment was 51 × 20 mm, 1 mm thick SFL6.
A gold thin film having a thickness of 45 nm is formed on a glass substrate within a range of 20 × 20 mm. In order to improve the adhesion between gold and glass, a chromium layer having a thickness of 5 nm was provided between the gold thin film and the glass substrate. Matching oil was used to reduce the reflectance of the glass surface when fixing the gold substrate to the semi-cylindrical prism.

【0016】次に、式(29)で与えられる非対称SPR
カーブ方程式について、近似式としての有効性とフィッ
ティング関数としての有益性について説明する。SPR実
験の計算シミュレーションにおいて、np=1.80の高屈折
率プリズムを使った金のSPRセンサーを用いて水溶液の
屈折率を測定したことを仮定した。水の典型的な屈折率
は、波長589.3nm、温度20゜Cでn=1.3330である。
分散の効果は、金では考慮しているが、プリズムと水溶
液では考慮していなかった。金の分散を示す光学パラメ
ータとして、次式で表される金属層の誘電率εに対し
てドルーデモデルを仮定した。
Next, the asymmetric SPR given by equation (29)
Regarding the curve equation, the effectiveness as an approximate expression and the usefulness as a fitting function will be described. In the computational simulation of the SPR experiment, it was assumed that the refractive index of the aqueous solution was measured using a gold SPR sensor with a high index prism with n p = 1.80. The typical refractive index n w of water is n w = 1.3330 at a wavelength of 589.3 nm and a temperature of 20 ° C.
The effect of dispersion was taken into account for gold, but not for prisms and aqueous solutions. As an optical parameter indicating the dispersion of gold, a Drude model was assumed for the dielectric constant ε m of the metal layer represented by the following equation.

【0017】図2及び図3には、式(29)の非対称SP
Rカーブ方程式を波長600nmと800nmについて調べた結果
を示す。図2は、励起波長600nmにおける非対称SPRカー
ブ方程式の数値計算を例示し、屈折率1.3330の水を高屈
折率1.80のプリズムで金のSPR測定することを仮定して
いる。図2の(a)のグラフは、式(29)で与えられ
る非対称SPRカーブの式(破線)とp偏光三層フレネル
式による完全な式(実線)とを、金薄膜の厚さ30、40、
50、60、70nmにおける5つのSPRカーブについて比較し
たものである。(b)のグラフは、完全な式で得られた
SPRカーブ(実線)を未知変数A、B、C、D、Eを持つ式
(29)の非対称SPRカーブ方程式でフィッティングす
る例を示している。(c)のグラフは5つのパラメータ
A、B、C、D、Eについて、フィッティングで得られたも
のを非対称方程式の式(30)〜(35)で得られるも
ので割った比率を示している。
2 and 3, the asymmetric SP of equation (29) is shown.
The result of having investigated the R curve equation about wavelength 600nm and 800nm is shown. FIG. 2 illustrates the numerical calculation of the asymmetric SPR curve equation at an excitation wavelength of 600 nm, and assumes that water with a refractive index of 1.3330 is used to measure the SPR of gold with a prism of high refractive index of 1.80. The graph of FIG. 2A shows the asymmetric SPR curve formula (broken line) given by formula (29) and the complete p-polarization three-layer Fresnel formula (solid line) for the thicknesses of the gold thin films 30, 40. ,
It is a comparison of 5 SPR curves at 50, 60 and 70 nm. The graph in (b) was obtained with a complete formula
An example of fitting the SPR curve (solid line) with the asymmetric SPR curve equation of the equation (29) having unknown variables A, B, C, D and E is shown. The graph in (c) has five parameters
For A, B, C, D, and E, the ratios obtained by fitting are divided by those obtained by equations (30) to (35) of the asymmetric equation.

【0018】図3は、励起波長800nmにおける非対称SPR
カーブ方程式の数値計算を例示し、励起波長以外は、図
2で用いた数値と同じ数値を計算で用いている。図2及
び図3において、実線はp偏光三層フレネル方程式から
得た原型SPRカーブを示しており、また破線は、式(2
9)〜式(35)を用いて非対称SPRカーブ方程式を近
似して求めたSPRカーブである。30、40、50、60、70nm
と異なる金膜厚のSPRカーブについて考察した。破線のS
PRカーブと実線の原型SPRカーブが良く一致しているの
は、図2においては膜厚が50、60、70nmの時であり、図
3においては膜厚が40、50、60、70nmの時である。近似
SPRカーブは原型SPRカーブとほぼ合っているが、ぴった
りとは合っていない。
FIG. 3 shows an asymmetric SPR at an excitation wavelength of 800 nm.
The numerical calculation of the curve equation is illustrated, and the same numerical values as those used in FIG. 2 are used in the calculation except the excitation wavelength. 2 and 3, the solid line shows the prototype SPR curve obtained from the p-polarized three-layer Fresnel equation, and the broken line shows the equation (2
9) is an SPR curve obtained by approximating an asymmetric SPR curve equation using Equations (35) to (35). 30, 40, 50, 60, 70nm
We examined the SPR curve with different gold film thickness. Dashed S
The PR curve and the solid prototype SPR curve are in good agreement when the film thickness is 50, 60, 70 nm in FIG. 2 and when the film thickness is 40, 50, 60, 70 nm in FIG. Is. Approximation
The SPR curve is close to the original SPR curve, but not exactly.

【0019】図2及び図3の(b)において実線で示し
た原型SPRカーブは未知変形A、B、C、D、Eを持つ式29の
非対称SPRカーブでフィッティングをした。フィッティ
ングは図2の(b)では52-62°の10°の範囲で行い、
図3の(b)では48-58°の10°の範囲で行った。破線
で示したフィッティングカーブが見えないのは、フィッ
ティングカーブと原型SPRカーブが良く一致しているか
らである。このことは、非対称SPRカーブ方程式がSPRカ
ーブのフィッティング関数として極めて有効であること
を意味している。
The prototype SPR curves shown by the solid lines in FIGS. 2 and 3 (b) were fitted with an asymmetric SPR curve of Equation 29 having unknown deformations A, B, C, D and E. The fitting is performed in the range of 10 ° of 52-62 ° in FIG.
In FIG. 3 (b), it was performed in the range of 48 ° to 58 °, 10 °. The fitting curve shown by the broken line is not visible because the fitting curve and the prototype SPR curve are in good agreement. This means that the asymmetric SPR curve equation is extremely effective as a fitting function for the SPR curve.

【0020】図2及び図3の(c)には、フィッティイ
ングにより得たA〜Eの5つのパラメータと、非対称SPR
カーブ方程式の式(30)〜(34)で与えられるA〜E
の5つのパラメータとの比が示されている。5つのパラ
メータの中で、パラメータDだけが0.2%以内で1となって
いる。パラメータDはセンシング層の屈折率の線型変化
により正しい情報を与えるため、特に重要である。理論
とフィッティングとにおけるパラメータDの正確な一致
は、近似が非常に効果的であり、式(29)の非対称SP
Rカーブ方程式がフィッティング関数として有益である
ことを示している。
2 and 3 (c), the five parameters A to E obtained by fitting and the asymmetric SPR are shown.
A to E given by equations (30) to (34) of the curve equation
The ratios of the five parameters are shown. Of the five parameters, only parameter D is 1 within 0.2%. The parameter D is particularly important because it gives correct information due to the linear variation of the refractive index of the sensing layer. The exact agreement of the parameter D between theory and fitting is very effective in the approximation, and the asymmetric SP of equation (29)
It shows that the R-curve equation is useful as a fitting function.

【0021】SPRカーブの反射の谷の実時間変化を観測
する実際のSPRカーブでは、角度の絶対的測定でなく、
相対的測定が通常行われる。その場合、式(29)で与
えられる非対称SPRカーブ方程式では、入射角θに関す
る絶対値の情報が必要となるので、データ処理における
数値計算には不適当である。この問題は、関数kx(θ)
をθ=θSPRにおいて一次までのテイラー展開に適用する
ことによって解決できる。 ここでXは入射角であり、ここでは置き換えについての
説明はは省略する。これは実際の実験に適したクレッチ
マンSPR方程式のフィッティング関数である。
In the actual SPR curve for observing the real-time change of the reflection valley of the SPR curve, not the absolute measurement of the angle but the
Relative measurements are usually made. In that case, the asymmetric SPR curve equation given by the equation (29) requires information on the absolute value of the incident angle θ, and is not suitable for numerical calculation in data processing. This problem is a function k x (θ)
Can be solved by applying the Taylor expansion up to the first order in θ = θ SPR . Here, X is the incident angle, and a description of replacement is omitted here. This is a fitting function of the Kretschmann SPR equation suitable for actual experiments.

【0022】図4及び図5には、SPRカーブのフィッテ
ィング関数としての式(45)の非対称SPRカーブ方程
式の有益性の評価を示す。図4は、元のSPRを励起波長6
00nmでp偏光三層フレネル式を用いて計算した場合の、
非対称SPRカーブ方程式でカーブフィッティングしたも
のを示し、1.333から1.368まで0.005間隔の屈折率を持
つ水溶液を50nm金薄膜、1.80屈折率プリズムの金のSPR
で測定すると仮定する。(a)のグラフは、実線で示さ
れた計算により得たSPRカーブを非対称SPRカーブ方程式
でフィッティングした例を示し、フィッティングしたSP
Rカーブは破線で示したが重なっているためほとんど見
えない。(b)のグラフは、フィッティングパラメータ
A、B、C、D、Eのセンシング層屈折率の依存性を示し、
フィッティングパラメータは、表1にしたセンシング層
の屈折率が1.333の時に得られるパラメータでそれぞれ
規格化されている。また(c)のグラフは、パラメータ
Dとして定義される修正共鳴角及び反射で最小値を与え
る従来の共鳴角の屈折率変化に対する応答を示し、従来
の共鳴角は式(47)で示されているように
4 and 5 show an evaluation of the usefulness of the asymmetric SPR curve equation of equation (45) as a fitting function for the SPR curve. Figure 4 shows the original SPR with excitation wavelength 6
When calculated using the p-polarized three-layer Fresnel equation at 00 nm,
The curve fitting by the asymmetric SPR curve equation is shown. A 50 nm gold thin film of an aqueous solution with a refractive index of 0.005 intervals from 1.333 to 1.368 is used for the SPR of gold with a 1.80 refractive index prism.
Suppose you measure at. The graph of (a) shows an example of fitting the SPR curve obtained by the calculation indicated by the solid line with the asymmetric SPR curve equation, and the fitted SP
The R-curves are shown by broken lines, but they are almost invisible because they overlap. The graph in (b) shows the fitting parameters.
A, B, C, D, E showing the dependence of the refractive index of the sensing layer,
The fitting parameters are standardized by the parameters shown in Table 1 obtained when the refractive index of the sensing layer is 1.333. The graph in (c) shows the parameters
The modified resonance angle defined as D and the response to the change in refractive index of the conventional resonance angle which gives the minimum value in reflection are shown in equation (47).

【0023】図5は、元のSPRを励起波長800nmでp偏光
三層フレネル式を用いて計算した場合の、非対称SPRカ
ーブ方程式でカーブフィッティングしたものを示し、数
値計算は、励起波長以外図5で用いたものと同じであ
る。
FIG. 5 shows the curve fitting by the asymmetric SPR curve equation when the original SPR was calculated using the p-polarization three-layer Fresnel equation at the excitation wavelength of 800 nm, and the numerical calculation was performed except for the excitation wavelength. The same as that used in.

【0024】非対称SPRカーブ方程式にフィッティング
する元となるSPRカーブとして、図4及び図5で示すよ
うに波長600nmと800nm励起におけるp偏光フレネル方程
式により計算したSPRカーブを使用した。計算では、厚
さ50nmの金薄膜及び屈折率1.80のプリズムを備えた金の
SPRを用いて1.333から1.368まで0.005間隔の屈折率を持
つ水溶液を測定したと仮定した。図4及び図5の(a)
では、元になるSPRカーブが実線で示されており、フィ
ッティングSPRカーブは破線で示されている。フィッテ
ィングはまた、図4の(a)では52〜62°、図5の
(a)では48〜58°のそれぞれ10°の範囲で行った。フ
ィッティングカーブは、元となるカーブとフィッティン
グカーブとがほとんど重なっているため、ほとんど見え
ない。図4及び図5の(b)は、フィッティングパラメ
ータA〜Eのセンシング層の屈折率への依存性を示して
いる。これら5つのフィッティングパラメータA〜Eは
表 1に示すように、屈折率が1.333の時に得られるフィ
ッティングパラメータでそれぞれ規格化されている。
As the SPR curve which is the basis for fitting the asymmetric SPR curve equation, the SPR curve calculated by the p-polarization Fresnel equation at the wavelengths of 600 nm and 800 nm as shown in FIGS. 4 and 5 was used. Calculations show that a gold film with a thickness of 50 nm and a prism with a refractive index of 1.80 is used.
It was assumed that an aqueous solution having a refractive index of 0.005 intervals from 1.333 to 1.368 was measured using SPR. (A) of FIG. 4 and FIG.
In, the original SPR curve is shown by a solid line, and the fitting SPR curve is shown by a broken line. The fitting was also carried out in the range of 52 ° to 62 ° in FIG. 4A and in the range of 10 ° to 48 ° in FIG. 5A. The fitting curve is almost invisible because the original curve and the fitting curve almost overlap. 4 and 5B show the dependence of the fitting parameters A to E on the refractive index of the sensing layer. As shown in Table 1, these five fitting parameters A to E are standardized by the fitting parameters obtained when the refractive index is 1.333.

【0025】 [0025]

【0026】1−|γ|<10-3となる線形相関係数γを
持つパラメータは、図4の(b)ではA、D、E、図5の
(b)ではDである。図4及び図5の(b)で用いたA、
B、C、D、Eの線形相関係数を表2に示す。
Parameters having a linear correlation coefficient γ such that 1− | γ | <10 −3 are A, D and E in FIG. 4B and D in FIG. 5B. A used in FIGS. 4 and 5B,
Table 2 shows the linear correlation coefficients of B, C, D and E.

【0027】 [0027]

【0028】2つの場合において、パラメータDだけが
条件1−|γ|<10-4を満たす。すなわち、パラメータD
は他のパラメータより屈折率に関して直線関係を持って
いることを意味している。
In both cases, only the parameter D satisfies the condition 1- | γ | <10 -4 . That is, the parameter D
Means that it has a linear relationship with respect to the refractive index than other parameters.

【0029】図4及び図5の(c)には、パラメータD
に全く等しい修正共鳴角を、数学的表式が式(47)で
与えられる従来の共鳴角と比較した例を示している。従
来の共鳴角と修正共鳴角の2つの直線は、お互いに平行
な関係を持っている。図4及び図5の(c)における4
つの直線全てが条件1−|γ|<10-4を満たす線形相関
係数γを持つ。図4及び図5の結果は、式(45)の非
対称SPRカーブ式がSPRカーブのフィッティング関数とし
て極めて有益であることを意味している。
The parameter D is shown in FIGS. 4 and 5C.
The modified resonance angle exactly equal to is compared with the conventional resonance angle given by the mathematical expression (47). The two lines of the conventional resonance angle and the corrected resonance angle have a parallel relationship with each other. 4 in (c) of FIGS. 4 and 5
All three straight lines have a linear correlation coefficient γ that satisfies the condition 1− | γ | <10 −4 . The results in FIGS. 4 and 5 mean that the asymmetric SPR curve equation of equation (45) is extremely useful as a fitting function for the SPR curve.

【0030】一方、本発明の理論は、修正共鳴角が従来
の共鳴角より線形性が良いことを予測している。図4の
(c)では、修正共鳴角が従来の共鳴角よりも、線形相
関係数で1.9×10-5だけ線形性が良いことを示してい
る。図5の(c)では、修正共鳴角と従来の共鳴角で線
形性は同じである。この結果は本発明の理論を支持する
ものである。
On the other hand, the theory of the present invention predicts that the modified resonance angle has better linearity than the conventional resonance angle. FIG. 4C shows that the corrected resonance angle is better than the conventional resonance angle by linear correlation coefficient of 1.9 × 10 −5 . In FIG. 5C, the linearity is the same between the modified resonance angle and the conventional resonance angle. This result supports the theory of the present invention.

【0031】さらに、元のSPRカーブに対する装置関数
の効果を、実際のSPR実験について説明する。図6及び
図7において、装置関数を二次式で与えられる箱型関数
のコンボルーション(たたみ込み)として含んでいる。
Further, the effect of the instrumental function on the original SPR curve will be described for an actual SPR experiment. 6 and 7, the device function is included as a convolution of a box-shaped function given by a quadratic expression.

【0032】図6は、装置関数を励起波長600nmの元のS
PRカーブに含んだ時、非対称SPRカーブ方程式にフィッ
ティングしたものを示し、1.333から1.368まで0.005間
隔の屈折率を持つ水溶液を金薄膜50nm、プリズム屈折率
1.80の金SPRで測定すると仮定している。図6の(a)
のグラフにおいて、実線で示されたSPRカーブは装置関
数を含んだ元のSPRカーブであり、破線で示すSPRカーブ
は式(45)で与えられる非対称SPRカーブ方程式でフ
ィッティングしたものであり、破線で示したSPRがほと
んど見ないのは、元のSPRカーブと重なっているためで
ある。(b)のグラフは修正共鳴角及び従来の共鳴角の
屈折率に対する応答を示している。
FIG. 6 shows the device function as the original S at an excitation wavelength of 600 nm.
When included in the PR curve, the one fitted to the asymmetric SPR curve equation is shown. An aqueous solution having a refractive index of 0.005 intervals from 1.333 to 1.368 is formed with a gold thin film of 50 nm and a prism refractive index.
It is assumed to measure with a gold SPR of 1.80. FIG. 6A
In the graph of, the SPR curve shown by the solid line is the original SPR curve including the device function, and the SPR curve shown by the broken line is the one fitted by the asymmetric SPR curve equation given by equation (45), The SPR shown is almost invisible because it overlaps with the original SPR curve. The graph in (b) shows the response of the modified resonance angle and the conventional resonance angle to the refractive index.

【0033】図7は、装置関数を励起波長800nmの元のS
PRカーブに含んだ時、非対称SPRカーブ方程式でフィッ
ティングしたものを示し、計算に用いた数値は励起波長
以外図6で用いたものと同じである。図7の(a)にお
いて実線で示されたSPRカーブは装置関数を含む元のSPR
カーブであり、一方破線で示されたSPRカーブは、非対
称SPRカーブ方程式にフィッティングしたものを示して
いる。また図7の(b)には修正共鳴角及び従来の共鳴
角の屈折率に対する応答を示している。
FIG. 7 shows the instrumental function as the original S at an excitation wavelength of 800 nm.
When included in the PR curve, the one fitted by the asymmetric SPR curve equation is shown, and the numerical values used for the calculation are the same as those used in FIG. 6 except the excitation wavelength. The SPR curve shown by the solid line in FIG. 7A is the original SPR including the device function.
The curve, while the SPR curve shown by the dashed line, shows the fitted asymmetric SPR curve equation. Further, FIG. 7B shows the responses of the modified resonance angle and the conventional resonance angle to the refractive index.

【0034】実際には、ソフトウェアのイゴールプロの
箱型スムージングアルゴリズムを使用した。図6の実線
で示した元のSPRカーブは、図4の(a)のものから箱
型スムージングアルゴリズムを用いて得られた。一方、
図7の(a)の元となるSPRカーブは、図5の(a)の
ものから得られた。図6の(a)において、フィッティ
ングカーブは元になるSPRカーブとぴったりと一致して
いる。それに対して、図7の(a)では、元となるカー
ブとフィッティングカーブは僅かに重ね合っているに過
ぎない。装置関数の効果により、図7の(a)において
は、元になるカーブとフィッティング関数の違いを生じ
させているが、図6の(a)では生じさせてない。この
違いが大きくなるのは、SPRカーブの幅が装置関数の幅
より小さくなる場合である。すなわち、非対称SPRカー
ブ方程式が装置関数を含んだSPRカーブに対して有効に
なるのは、SPRカーブの幅が装置関数よりもずっと大き
い場合である。図7において、元のカーブとフィッティ
ングのカーブとが僅かに重ね合っているにすぎないが、
図7の(b)は、非対称SPRカーブ方程式が共鳴角の決
定において依然として有効であることを示している。こ
のことは、元になるカーブとフィッティングカーブの多
少の違いは、フィッティング関数において共鳴角を決定
するには許容されるということである。図6の(b)及
び図7の(b)において、修正共鳴角は、従来の共鳴角
より、線形性がそれぞれ2.8×10-5、3×10-6だけ良かっ
た。この結果は、修正共鳴角が、従来の共鳴角よりも屈
折率変化を正確に決めるのに適していることを示してい
る。
In practice, IgorPro's box-type smoothing algorithm of the software was used. The original SPR curve shown by the solid line in FIG. 6 was obtained from the one in FIG. 4 (a) using the box smoothing algorithm. on the other hand,
The SPR curve which is the basis of (a) of FIG. 7 was obtained from that of (a) of FIG. In FIG. 6A, the fitting curve exactly matches the original SPR curve. On the other hand, in FIG. 7A, the original curve and the fitting curve are only slightly overlapped. Due to the effect of the device function, the difference between the original curve and the fitting function is generated in FIG. 7A, but it is not generated in FIG. 6A. This difference becomes large when the width of the SPR curve is smaller than the width of the device function. That is, the asymmetric SPR curve equation is valid for an SPR curve including a device function when the width of the SPR curve is much larger than the device function. In FIG. 7, the original curve and the fitting curve are only slightly overlapped,
FIG. 7 (b) shows that the asymmetric SPR curve equation is still valid in determining the resonance angle. This means that some difference between the original curve and the fitting curve is acceptable for determining the resonance angle in the fitting function. In (b) of FIG. 6 and (b) of FIG. 7, the corrected resonance angle was better than the conventional resonance angle by linearity of 2.8 × 10 −5 and 3 × 10 −6 , respectively. This result shows that the modified resonance angle is more suitable for determining the refractive index change more accurately than the conventional resonance angle.

【0035】最後に、図8及び図9を参照して非対称SP
Rカーブ方程式について、0、1、4、10、20重量
%の5つの濃度の蔗糖水溶液を用いた実際の実験に基き
説明する。図8は、励起波長630nmで得られるSPRカーブ
をフィッティングすることによる非対称SPRカーブ方程
式の実験的考察を、濃度0、1、4、10、20重量%
の蔗糖水溶液の測定によって示し、(a)のグラフにお
いて、実験のSPRカーブは実線で示され、破線は式(4
5)の非対称SPRカーブ方程式でのフィッティングしたSP
Rカーブを示している。また(b)のグラフは修正共鳴
角及び従来の共鳴角の屈折率応答を示し、修正共鳴角
は、非対称SPR方程式のフィッティングでパラメータDと
して与えられる。従来の共鳴角は、
Finally, referring to FIGS. 8 and 9, the asymmetric SP
The R-curve equation will be described based on actual experiments using five concentrations of sucrose aqueous solutions of 0, 1, 4, 10, and 20% by weight. FIG. 8 shows the experimental consideration of the asymmetric SPR curve equation by fitting the SPR curve obtained at the excitation wavelength of 630 nm at concentrations of 0, 1, 4, 10, 20% by weight.
In the graph of (a), the SPR curve of the experiment is shown by the solid line, and the broken line is shown by the formula (4).
5) SP fitted with asymmetric SPR curve equation
The R curve is shown. The graph of (b) shows the refractive index response of the modified resonance angle and the conventional resonance angle, and the modified resonance angle is given as the parameter D in the fitting of the asymmetric SPR equation. The conventional resonance angle is

【0036】図9は、励起波長880nmで得られるSPRカー
ブをフィッティングすることによる非対称SPRカーブ方
程式の実験的考察を示し、測定には濃度0、1、4、1
0、20重量%の蔗糖水溶液を用いた。
FIG. 9 shows an experimental consideration of the asymmetric SPR curve equation by fitting the SPR curve obtained at the excitation wavelength of 880 nm, and the concentrations of 0, 1, 4, 1 were measured.
An aqueous sucrose solution of 0 or 20% by weight was used.

【0037】図8及び図9の(a)における実線は、63
0と880nm励起波長の実験で得られたものであり、破線は
非対称SPRカーブ方程式によるフィッティングカーブで
ある。表3には、図8及び図9の(a)における蔗糖0
重量%のSPRカーブをフィッティングする際に得られた
5つのパラメータを示す。
The solid line in FIGS. 8 and 9 (a) is 63.
It was obtained in the experiment of 0 and 880 nm excitation wavelength, and the broken line is the fitting curve by the asymmetric SPR curve equation. Table 3 shows that sucrose 0 in FIG.
The five parameters obtained when fitting the wt% SPR curve are shown.

【0038】 [0038]

【0039】注目すべきは、半円柱プリズムや円柱レン
ズによる光学歪みや、光源のスペクトル広がりのため
に、理論上のSPRカーブと実験のSPRカーブが大きく違っ
ているにもかかわらず、フィッティングカーブが実験の
カーブと良く一致していることである。実験結果は、本
発明における非対称SPRフィッティングカーブ方程式
が、金のSPRカーブをフィッティングするのに最も適当
であることを示している。つまり、5つよりも少ないパ
ラメータのフィッティングでは、実験のSPRカーブをぴ
ったりとフィッティングすることは不可能であり、5つ
よりも多いパラメータによるフィッティングは、SPRカ
ーブの無意味な情報を与えてしまうということである。
It should be noted that, although the theoretical SPR curve and the experimental SPR curve are greatly different due to the optical distortion due to the semi-cylindrical prism or the cylindrical lens and the spectral spread of the light source, the fitting curve is It is in good agreement with the experimental curve. Experimental results show that the asymmetric SPR fitting curve equation in the present invention is most suitable for fitting the gold SPR curve. In other words, it is impossible to exactly fit the experimental SPR curve with fitting of less than 5 parameters, and fitting with more than 5 parameters gives meaningless information on the SPR curve. That is.

【0040】図8及び図9の(b)は、20゜C、589nm波
長における蔗糖水溶液の屈折率を関数として、パラメー
タDに等しい修正共鳴角と従来の共鳴角をプロットした
ものである。本発明における実際の測定条件は、参考デ
ータの波長とは明らかに異なるものの、純水に対する蔗
糖水溶液の屈折率の相対変化は波長にかかわらずほとん
ど同じなので、線形応答に関する本発明における考察は
正しいことが認められる。図8及び図9の(b)におい
て、修正共鳴角は従来の共鳴角よりも、2.0×10-4およ
び1.4×10-4だけ線形相関係数が良い。この結果は、本
発明の理論が数値計算の場合よりも実際の実験の方が有
効であることを示している。
FIG. 8 and FIG. 9 (b) are plots of the modified resonance angle equal to the parameter D and the conventional resonance angle as a function of the refractive index of the sucrose aqueous solution at 20 ° C. and 589 nm wavelength. Although the actual measurement conditions in the present invention are clearly different from the wavelength of the reference data, the relative change in the refractive index of the sucrose aqueous solution with respect to pure water is almost the same regardless of the wavelength, and therefore the consideration in the present invention regarding the linear response is correct. Is recognized. In FIGS. 8 and 9B, the corrected resonance angle has a linear correlation coefficient of 2.0 × 10 −4 and 1.4 × 10 −4 better than the conventional resonance angle. This result shows that the theory of the present invention is more effective in actual experiments than in the case of numerical calculation.

【0041】図2から図7までの数値計算シミュレーシ
ョンと図8及び図9の実験解析は、非対称SPRカーブ方
程式がSPRカーブでのフィッティング関数として極めて
有益であることを示している。また、センシング層の屈
折率が線形に変化するのを観測するのには、非対称SPR
カーブ方程式にフィッティングしてパラメータDとして
得られる修正共鳴角の方が、従来の共鳴角よりも優れて
いることを示している。本発明におけるフィッティング
関数は、また、光吸収型SPRのSPRカーブに対しても有効
である。
The numerical simulations of FIGS. 2 to 7 and the experimental analyzes of FIGS. 8 and 9 show that the asymmetric SPR curve equation is extremely useful as a fitting function in the SPR curve. Also, to observe that the refractive index of the sensing layer changes linearly, an asymmetric SPR
It is shown that the modified resonance angle obtained by fitting the curve equation as the parameter D is superior to the conventional resonance angle. The fitting function in the present invention is also effective for the SPR curve of the light absorption type SPR.

【0042】[0042]

【発明の効果】以上説明してきたように、本発明によれ
ば、プリズム−金属薄膜−センシング層の三層から成る
表面プラズモン共鳴センサーシステムを用い、表面プラ
ズモン共鳴現象を利用して物質状態を測定する際にP偏
光の三層フレネルの式から導いた式 で表される非対称表面プラズモン共鳴カーブのフィッテ
ィング方程式に基き表面プラズモン共鳴角を正確に決定
するように構成していので、クレッチマンによる式より
も適切なSPRカーブ方程式として利用でき、特に金を利
用したSPRのSPRカーブをフィッティングする関数として
極めて有効であり、SPRカーブより屈折率の線形変化を
決定することができ、またフィッティングは通常の金の
SPRセンサーで得られる非対称SPRカーブに対して有効で
あり、銀のSPRで得られる深いピークの対称SPRカーブに
対しても当然有効である。つまり、フィッティング方程
式でパラメータDとして定義される修正共鳴角は屈折率
の線形変化を観測する最適な変数として用いることがで
き、従って本発明による方法は、分析物を正確に決める
SPR装置のデータ処理において極めて有用である。
As described above, according to the present invention, the surface plasmon resonance sensor system comprising three layers of prism-metal thin film-sensing layer is used to measure the substance state by utilizing the surface plasmon resonance phenomenon. Equation derived from the P-polarized three-layer Fresnel equation Since it is configured to accurately determine the surface plasmon resonance angle based on the fitting equation of the asymmetric surface plasmon resonance curve represented by, it can be used as a more appropriate SPR curve equation than the expression by Kletschmann, and especially SPR using gold. It is extremely effective as a function for fitting the SPR curve of, and the linear change of the refractive index can be determined from the SPR curve.
It is effective for the asymmetric SPR curve obtained by the SPR sensor, and naturally effective for the deep peak symmetrical SPR curve obtained by the silver SPR. That is, the modified resonance angle, defined as the parameter D in the fitting equation, can be used as an optimal variable for observing linear changes in the refractive index, and thus the method according to the present invention accurately determines the analyte.
Very useful in data processing of SPR equipment.

【図面の簡単な説明】[Brief description of drawings]

【図1】プリズム−金属−センシング層の界面から構成
されるクレッチマン配置を示す概略断面図。
FIG. 1 is a schematic cross-sectional view showing a Kretschmann arrangement composed of a prism-metal-sensing layer interface.

【図2】励起波長600nmにおける非対称SPRカーブ方程式
の数値計算例を示し、(a)は非対称SPRカーブの式
(破線)とp偏光三層フレネル式による完全な式(実
線)とのSPRカーブ比較を示すグラフ、(b)は完全な
式で得られたSPRカーブ(実線)を未知変数A、B、C、
D、Eを持つ非対称SPRカーブ方程式でフィッティングす
る場合を例示するグラフ、また(c)は5つのパラメー
タA、B、C、D、Eについて、フィッティングで得られた
ものを非対称方程式で得られるもので割った比率を示す
グラフ。
FIG. 2 shows an example of numerical calculation of an asymmetric SPR curve equation at an excitation wavelength of 600 nm, (a) shows an SPR curve comparison between the asymmetric SPR curve equation (dashed line) and the complete p-polarization three-layer Fresnel equation (solid line). (B) is a graph showing the SPR curve (solid line) obtained by a complete equation with unknown variables A, B, C,
A graph exemplifying the case of fitting with an asymmetric SPR curve equation having D and E, and (c) is the one obtained by the fitting with respect to the five parameters A, B, C, D and E. Graph showing the ratio divided by.

【図3】励起波長800nmにおける非対称SPRカーブ方程式
の数値計算例を示し、(a)は非対称SPRカーブの式
(破線)とp偏光三層フレネル式による完全な式(実
線)とのSPRカーブ比較を示すグラフ、(b)は完全な
式で得られたSPRカーブ(実線)を未知変数A、B、C、
D、Eを持つ非対称SPRカーブ方程式でフィッティングす
る場合を例示するグラフ、また(c)は5つのパラメー
タA、B、C、D、Eについて、フィッティングで得られた
ものを非対称方程式で得られるもので割った比率を示す
グラフ。
FIG. 3 shows a numerical calculation example of an asymmetric SPR curve equation at an excitation wavelength of 800 nm, (a) shows an SPR curve comparison between an asymmetric SPR curve equation (dashed line) and a complete p-polarization three-layer Fresnel equation (solid line). (B) is a graph showing the SPR curve (solid line) obtained by a complete equation with unknown variables A, B, C,
A graph exemplifying the case of fitting with an asymmetric SPR curve equation having D and E, and (c) is the one obtained by the fitting with respect to the five parameters A, B, C, D and E. Graph showing the ratio divided by.

【図4】元のSPRを励起波長600nmでp偏光三層フレネル
式を用いて計算した場合の、非対称SPRカーブ方程式で
カーブフィッティングしたものを示し、(a)は、実線
で示された計算により得たSPRカーブを非対称SPRカーブ
方程式でフィッティングした例を示すグラフ、(b)
は、フィッティングパラメータA、B、C、D、Eのセンシ
ング層屈折率の依存性を示すグラフ、また(c)は、パ
ラメータDとして定義される修正共鳴角及び反射で最小
値を与える従来の共鳴角の屈折率変化に対する応答を示
すグラフである。
FIG. 4 shows curve fitting by an asymmetric SPR curve equation when the original SPR was calculated using a p-polarized three-layer Fresnel equation at an excitation wavelength of 600 nm, and (a) is calculated by a solid line. A graph showing an example of fitting the obtained SPR curve with an asymmetric SPR curve equation, (b)
Is a graph showing the dependency of the fitting parameters A, B, C, D, and E on the refractive index of the sensing layer, and (c) is the conventional resonance that gives the minimum value at the modified resonance angle and reflection defined as the parameter D. It is a graph which shows the response with respect to the refractive index change of a corner.

【図5】元のSPRを励起波長800nmでp偏光三層フレネル
式を用いて計算した場合の、非対称SPRカーブ方程式で
カーブフィッティングしたものを示し、(a)は、実線
で示された計算により得たSPRカーブを非対称SPRカーブ
方程式でフィッティングした例を示すグラフであり、
(b)は、フィッティングパラメータA、B、C、D、Eの
センシング層屈折率の依存性を示すグラフであり、また
(c)は、パラメータDとして定義される修正共鳴角及
び反射で最小値を与える従来の共鳴角の屈折率変化に対
する応答を示すグラフである。
FIG. 5 shows curve fitting by an asymmetric SPR curve equation when the original SPR is calculated using a p-polarized three-layer Fresnel equation at an excitation wavelength of 800 nm, and (a) shows the calculation by the solid line. It is a graph showing an example of fitting the obtained SPR curve with an asymmetric SPR curve equation,
(B) is a graph showing the dependence of the fitting parameters A, B, C, D, and E on the refractive index of the sensing layer, and (c) is the modified resonance angle defined as the parameter D and the minimum value at reflection. 5 is a graph showing a response to a change in the refractive index of the conventional resonance angle that gives

【図6】装置関数を励起波長600nmの元のSPRカーブに含
んだ時、非対称SPRカーブ方程式にフィッティングした
ものを示し、(a)は装置関数を含んだ元のSPRカーブ
である実線で示されたSPRカーブ及び非対称SPRカーブ方
程式でフィッティングした破線で示すSPRカーブを示す
グラフであり、また(b)は修正共鳴角及び従来の共鳴
角の屈折率に対する応答を示すグラフである。
FIG. 6 shows the fitting to an asymmetric SPR curve equation when the instrument function was included in the original SPR curve with an excitation wavelength of 600 nm, and (a) is the original SPR curve including the instrument function, which is shown by a solid line. 2B is a graph showing the SPR curve shown by a broken line fitted with the SPR curve and the asymmetric SPR curve equation, and (b) is a graph showing the response to the refractive index of the modified resonance angle and the conventional resonance angle.

【図7】装置関数を励起波長800nmの元のSPRカーブに含
んだ時、非対称SPRカーブ方程式でフィッティングした
ものを示し、(a)は装置関数を含んだ元のSPRカーブ
である実線で示されたSPRカーブ及び非対称SPRカーブ方
程式でフィッティングした破線で示すSPRカーブを示す
グラフであり、また(b)は修正共鳴角及び従来の共鳴
角の屈折率に対する応答を示すグラフである。
FIG. 7 shows the fitting by an asymmetric SPR curve equation when the instrument function was included in the original SPR curve with the excitation wavelength of 800 nm, and (a) is shown by the solid line which is the original SPR curve including the instrument function. 2B is a graph showing the SPR curve shown by a broken line fitted with the SPR curve and the asymmetric SPR curve equation, and (b) is a graph showing the response to the refractive index of the modified resonance angle and the conventional resonance angle.

【図8】励起波長630nmで得られるSPRカーブをフィッテ
ィングすることによる非対称SPRカーブ方程式を、濃度
0、1、4、10、20重量%の蔗糖水溶液の測定によ
る実験的考察を示し、(a)は実験のSPRカーブ(実
線)及び非対称SPRカーブ方程式でのフィッティングし
たSPRカーブ(破線)を示すグラフであり、また(b)
は修正共鳴角及び従来の共鳴角の屈折率応答を示すグラ
フである。
FIG. 8 shows an asymmetric SPR curve equation obtained by fitting an SPR curve obtained at an excitation wavelength of 630 nm, and experimental consideration by measurement of a sucrose aqueous solution having a concentration of 0, 1, 4, 10, or 20% by weight, (a). Is a graph showing the experimental SPR curve (solid line) and the fitted SPR curve (dashed line) in the asymmetric SPR curve equation, and (b)
6 is a graph showing the refractive index response of a modified resonance angle and a conventional resonance angle.

【図9】励起波長880nmで得られるSPRカーブをフィッテ
ィングすることによる非対称SPRカーブ方程式の実験的
を示し、(a)は実験のSPRカーブ(実線)及び非対称S
PRカーブ方程式でのフィッティングしたSPRカーブ(破
線)を示すグラフであり、また(b)は修正共鳴角及び
従来の共鳴角の屈折率応答を示すグラフである。
FIG. 9 shows experimentally an asymmetric SPR curve equation by fitting an SPR curve obtained at an excitation wavelength of 880 nm, (a) shows an experimental SPR curve (solid line) and an asymmetric S
FIG. 6 is a graph showing a fitted SPR curve (dashed line) in the PR curve equation, and FIG. 6B is a graph showing the refractive index response of the modified resonance angle and the conventional resonance angle.

【符号の説明】[Explanation of symbols]

p プリズム m 金属層 s センシング層 θ 入射角 d 金属膜厚 p prism m Metal layer s Sensing layer θ incident angle d Metal film thickness

───────────────────────────────────────────────────── フロントページの続き (72)発明者 栗原 一嘉 神奈川県川崎市中原区井田杉山町4−1− 305クレールメゾン大瀬戸 Fターム(参考) 2G059 AA05 BB12 CC16 EE02 EE04 EE05 GG02 GG04 HH01 HH02 HH06 JJ11 JJ12 JJ19 KK04 MM01    ─────────────────────────────────────────────────── ─── Continued front page    (72) Inventor Kazuyoshi Kurihara             4-1 Idasugiyama-cho, Nakahara-ku, Kawasaki-shi, Kanagawa             305 Claire Maison Oseto F term (reference) 2G059 AA05 BB12 CC16 EE02 EE04                       EE05 GG02 GG04 HH01 HH02                       HH06 JJ11 JJ12 JJ19 KK04                       MM01

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】プリズム−金属薄膜−センシング層の三層
から成る表面プラズモン共鳴センサーシステムを用い、
表面プラズモン共鳴現象を利用したセンシングにおける
表面プラズモン共鳴カーブの非対称表面プラズモン共鳴
カーブ方程式による決定方法において、A、B、C、D、E
の5つの実数のうち少なくとも1つ以上を実数パラメー
タとし、sを非対称表面プラズモン共鳴カーブの特性因
子とし、Rを表面プラズモン共鳴センサーシステムにお
ける光の反射率、rpmsをプリズム−金属薄膜−センシ
ング層の境界における振幅反射率とし、kxを境界面に
平行な入射光の波数ベクトルとし、koをプリズムの存
在しない時の金属薄膜とセンシング層との界面に理論的
に生成される表面プラズモンの複素波数とし、kRをプ
リズムの存在する時のkoの補正項とするとき、P偏光の
三層フレネルの式から導いた式 で表される非対称表面プラズモン共鳴カーブのフィッテ
ィング方程式に基き表面プラズモン共鳴角を決定するこ
とを特徴とする方法。
1. A surface plasmon resonance sensor system comprising three layers of a prism, a metal thin film and a sensing layer is used.
A, B, C, D, E in the method of determining the surface plasmon resonance curve in the sensing using the surface plasmon resonance phenomenon by the asymmetric surface plasmon resonance curve equation
Of at least one of the five real numbers of ## EQU1 ## as a real number parameter, s as a characteristic factor of the asymmetric surface plasmon resonance curve, R as the light reflectance in the surface plasmon resonance sensor system, and r pms as the prism-metal thin film-sensing layer. of the amplitude reflectance at the boundary, the k x and wave vector parallel incident light on the boundary surface, k o the theoretically generated surface plasmons at the interface between the metal thin film and the sensing layer in the absence of the prism An equation derived from the P-polarized three-layer Fresnel equation, where complex wave number and k R are the correction terms for k o when a prism is present. A method of determining a surface plasmon resonance angle based on a fitting equation of an asymmetric surface plasmon resonance curve represented by.
【請求項2】実験から得られた表面プラズモン共鳴カー
ブのフィッティングにおいて、反射率データをR、入射
角データをXとするとき、上記フィッティング方程式が で表され、このフィッティング方程式に基き表面プラズ
モン共鳴角を決定する請求項1に記載の方法。
2. In the fitting of a surface plasmon resonance curve obtained from an experiment, when the reflectance data is R and the incident angle data is X, the fitting equation is The method according to claim 1, wherein the surface plasmon resonance angle is determined based on this fitting equation.
【請求項3】最小反射率に対する角度を で定義し、これを表面プラズモン共鳴角とする請求項1
又は2に記載の方法。
3. The angle for minimum reflectance The surface plasmon resonance angle is defined as
Or the method described in 2.
【請求項4】パラメータDが修正共鳴角として定義され
る請求項1に記載の方法。
4. The method according to claim 1, wherein the parameter D is defined as a modified resonance angle.
【請求項5】観測した表面プラズモン共鳴カーブをコン
ピュータ上でソフトウエア的に処理して共鳴角を求める
請求項1〜4のいずれか一項に記載の方法。
5. The method according to claim 1, wherein the observed surface plasmon resonance curve is processed by software on a computer to obtain a resonance angle.
JP2001202984A 2001-07-04 2001-07-04 Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation Pending JP2003014623A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001202984A JP2003014623A (en) 2001-07-04 2001-07-04 Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001202984A JP2003014623A (en) 2001-07-04 2001-07-04 Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation

Publications (1)

Publication Number Publication Date
JP2003014623A true JP2003014623A (en) 2003-01-15

Family

ID=19039691

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001202984A Pending JP2003014623A (en) 2001-07-04 2001-07-04 Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation

Country Status (1)

Country Link
JP (1) JP2003014623A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006125860A (en) * 2004-10-26 2006-05-18 Fujikura Ltd Surface plasmon sensor, and apparatus for measuring surface plasmon
JP2009115623A (en) * 2007-11-06 2009-05-28 Sharp Corp Surface plasmon sensor
JP2009236690A (en) * 2008-03-27 2009-10-15 Sharp Corp Surface plasmon sensor
US7830521B2 (en) 2006-09-29 2010-11-09 Fujifilm Corporation Detection apparatus, detection method, and optically transparent member
JP2013093596A (en) * 2007-12-21 2013-05-16 Asml Holding Nv Beam characterization monitor for sensing properties of optical beam
KR101667150B1 (en) * 2015-05-10 2016-10-17 고려대학교 산학협력단 Surface plasmon resonance curve fitting method and analyzing method of target substance using spr sensor
CN115218796A (en) * 2022-08-29 2022-10-21 重庆市计量质量检测研究院 DEABC method for detecting parameters of multilayer film based on in-situ common-angle SPR

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006125860A (en) * 2004-10-26 2006-05-18 Fujikura Ltd Surface plasmon sensor, and apparatus for measuring surface plasmon
US7830521B2 (en) 2006-09-29 2010-11-09 Fujifilm Corporation Detection apparatus, detection method, and optically transparent member
JP2009115623A (en) * 2007-11-06 2009-05-28 Sharp Corp Surface plasmon sensor
JP2013093596A (en) * 2007-12-21 2013-05-16 Asml Holding Nv Beam characterization monitor for sensing properties of optical beam
JP2009236690A (en) * 2008-03-27 2009-10-15 Sharp Corp Surface plasmon sensor
KR101667150B1 (en) * 2015-05-10 2016-10-17 고려대학교 산학협력단 Surface plasmon resonance curve fitting method and analyzing method of target substance using spr sensor
CN115218796A (en) * 2022-08-29 2022-10-21 重庆市计量质量检测研究院 DEABC method for detecting parameters of multilayer film based on in-situ common-angle SPR

Similar Documents

Publication Publication Date Title
Kurihara et al. Asymmetric SPR sensor response curve-fitting equation for the accurate determination of SPR resonance angle
Nenninger et al. Data analysis for optical sensors based on spectroscopy of surface plasmons
US9304046B2 (en) Sensor system and method for characterizing a coated body
Suzuki et al. Effects of gold film thickness on spectrum profile and sensitivity of a multimode-optical-fiber SPR sensor
Johansen et al. Imaging surface plasmon resonance sensor based on multiple wavelengths: Sensitivity considerations
US20120038920A2 (en) Method and Apparatus for Ellipsometry Measurement
US9423791B2 (en) Sensor system and method for determining paper sheet quality parameters
CA2409486A1 (en) Coupled plasmon-waveguide resonance spectroscopic device
CN103743349B (en) Method and device for measuring nano film
Jääskeläinen et al. On reflectometric measurement of a refractive index of milk
CN109470154B (en) Film thickness initial value measuring method suitable for spectrum ellipsometer
CN112461787B (en) Lithium niobate optical sensor and method based on Bloch surface wave
Uddin et al. A theoretical model for determination of optimum metal thickness in kretschmann configuration based surface plasmon resonance biosensors
Kalas et al. Bloch surface waves biosensing in the ultraviolet wavelength range–Bragg structure design for investigating protein adsorption by in situ Kretschmann-Raether ellipsometry
Yue et al. Determination of thin metal film’s thickness and optical constants based on SPR phase detection by simulated annealing particle swarm optimization
Kalas et al. Protein adsorption monitored by plasmon-enhanced semi-cylindrical Kretschmann ellipsometry
JP2003014623A (en) Method of deciding surface plasmon resonance curve for sensing utilizing surface plasmon resonance phenomenon based on asymmetric surface plasmon resonance curve equation
Luna-Moreno et al. Virtual instrumentation in LabVIEW for multiple optical characterizations on the same opto-mechanical system
US7362435B1 (en) Method of determining initial thickness value for sample surface layer, and use of splines for fitting ellipsometric data
CN103674892B (en) A kind of method carrying out monitoring film growth based on total internal reflection polarization phasic difference measurement
Luo et al. Sensitivity enhancement of surface plasmon resonance sensor based on wavelength and angular combined modulations
Urban III et al. Numerical ellipsometry: Analysis of thin metal layers using n–k plane methods with multiple incidence angles
Meyzonnette et al. Refractive index of optical materials
Karabchevsky et al. Techniques for signal analysis in surface plasmon resonance sensors
Goldina Calculation of the reflection coefficient of metal-dielectric structures in frustrated total internal reflection

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20031031

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20040129