物理工学科博士前期課程

H21 前期光計測工学(谷）１．非線形光学（I）(注 1a)
１．１

はじめに

本章では非線形光学現象を取り扱うための基礎的な事柄を学ぶ。
光の強度が低い場合は、物質は光に対して線形な応答を示し、媒質の誘電率は一定であり、媒質を透過した光の周波数は透過する前と同じものが得られる。しかし、光の強度が強くなると、物質と光による相互作用が強まり、媒質の誘電率は光強度に依存して変化する。物質の誘電率は電磁波の電界によって誘起される誘電分極によって決まるが、誘電率が変化するのはこの誘電分極が電磁波の電界 E に比例(線形応答)しなくなるためである。
電磁波は振動電界であるため、誘起された誘電分極 P（線形分極＋非線形分極）も振動する。分極 P が振動すると分極 P に比例(十分遠方では 2 次の時間微分に比例）した電磁波が放射されるが、P の非線形性のためにその周波数成分には入射する電磁波の周波数と一致するものだけではなく、入射波の 2 倍の周波数のものや、２つの入射波がある場合はそれらの和および差周波数のものなどが含まれる。新たに発生した複数の電磁波成分は媒質のなかで入射波と非線形な相互作用を通じてエネルギーのやりとりをしたあと、媒質を透過してくるが、どの周波数成分がどの程度強く得られるかは、入射波の強度、媒質の非線形応答の強さ、後に述べる位相整合条件などで決まる。
光は主として物質中の電子と相互作用する(注 1b)が、非線形な効果を得るためには、光の振動電界が物質中にある電子をその平衡点からポテンシャルが調和型でなくなる領域、すなわち平衡点からの距離を x としたとき、ポテンシャルが px2（p はポテンシャルで決まる定数）で近似できなくなる領域にまで変位させる必要がある。このためには非常に大きな電界が必要で、通常自然界にある光では、強度が十分ではなく非線形光学効果を観測するまでには至らない。しかし、レーザーの出現により、非常に強度が強く、しかもコヒーレンス(注 1c)性の良い（波の位相がそろっているという意味）光源を我々は利用することができるようになったので、非線形光学効果を用いて、１つ（または２つ）の波長から別の波
長の光を生成することが可能になった。
現在、半導体レーザーをはじめとしたレーザー技術の進歩が著しいが、あらゆる波長領域で発振が得られているわけではない。したがって、非線形光学効果による波長変換技術
はレーザー発振が容易でない波長域のコヒーレント光を得るために重要な技術である。また非線形光学効果は電気光学効果による光の変調、ソリトン波の伝播、Photorefractive 効果による光メモリーなど光通信や光エレクトロニクスにおいてさまざまな形で応用されている。
１．２

非線形感受率とその起源

媒質の誘電率は線形光学（通常の電磁気学）では
注１a：本テキストの非線形光学の記述は Robert W. Boyd の“Nonlinear Optics,”Seconde Edition（Academice Press）の 1 章、2 章、11 章にほぼ相当する。ただしガウス単位系ではなく MKS 有
理単位系を用いる。ちなみにガウス単位系では真空の誘電率ε0 と真空の透磁率μ0 が１なので式にはあらわに現れない。したがって、どちらの単位系を用いているかはε0 が式中に現れてくるかこないかでほぼ見分けがつく。
注１b：電子が光と相互作用しやすいのは、電子が原子核に比べてずっと軽くて、動きやすいからである。

1

物理工学科博士前期課程

H21 前期光計測工学(谷）

~
~ ~
~
~
~
D ≡ ε 0 E + P = ε 0 E + χ e ε 0 E = (1 + χ e )ε 0 E (線形な電束密度)(1.2.1)で定義される電束密度 D から

ε ≡ (1 + χ e )ε 0 = κ e ε 0(線形な誘電率)(1.2.2)で定義されるのであった。ここでε0 は真空の誘電率である。また各ベクトル量の頭にチル
ダ ” ～ ” を 付 け 、 そ れ が 時 間 的 に 振 動 す る 位 相 項 exp(−iωt) を 含 む も の （ 例 え ば~E(t ) = E exp(−iωt ) = A exp(−iωt + ikr ) など）であることを示すようにし、時間的に振動する位相項を含まない振幅成分（先の例では E に対応）と区別して書くことにする。ここ~~で(1.2.1)式における誘電分極 P はその点における電場 E に比例するとした(注 2)。

~~
P(t ) = χ e ε 0 E(t )(線形な誘電分極)(1.2.3)~~

しかし、光強度が強くなり、したがって電磁波の電界強度 E が大きくなると誘電分極 P は(1.2.3)式の比例関係で与えられるのではなく、~~P = χ ( E )ε 0

(1.2.4)~という式で一般に与えられる（ここで χ は電気感受率というよりは、単に E に依存する関~~数である）。とりあえず簡単のため、 P および E をスカラー量とし、この式を E についてのべきで展開すると、~~~~P (t ) = ε 0 ( χ (1) E (t ) + χ ( 2 ) E 2 (t ) + χ ( 3) E 3 (t ) + ")~~~≡ P (1) (t ) + P ( 2 ) (t ) + P ( 3) (t ) + "(1.2.5)と表すことができる。(1.2.5)式と(1.2.3)式を比べれば分かるように線形な感受率(linear
susceptibility)は(1.2.5)式の第 1 項の係数 χ
とχ

( 3)

(1)

に対応している。第 2 項と第 3 項の係数 χ( 2)はそれぞれ、2 次の非線形感受率(second-order nonlinear susceptibility)、3 次の非線形感受率(third-order nonlinear susceptibility)と呼ばれる。これらの電気感受率の係数は実際にはスカラーではなく、一般にテンソル量として表される。すなわち、直交座標系~~

において線形な感受率は分極 P の各成分に対応する座標係数 i と電場 E の各成分に対応する座標係数 j により、2 階のテンソルとして χ ij(1)などと表される（誘電率の異方性に対応）。

~

同様に 2 次および 3 次の非線形感受率は各分極成分 Pi に対して、それぞれ２つの電場成分（Ej と Ek）及び３つの電場成分（Ej, Ek,及び El）との積で与えられるため、 χ ijk

( 2)

および

注１c：
“コヒーレンス”は「可干渉性」と訳される。波長のまったく異なる波は干渉しないので、コヒーレント光と言うとき、単色性の良い光を意味する場合もあるが、その本質は波の位相がそろっていることである。
注 2a：物質の電磁波に対する線形応答は「電磁気学」のおさらいである。

2

物理工学科博士前期課程

H21 前期光計測工学(谷）electronicresonanceχ acoustic
nonlinear susceptibility

phonon resonanceabsorption

opticalphonon resonanceωfrequency
図 1.2.1 非線形感受率の周波数依存性と共鳴周波数

U
U=px2+qx3
U=px2
P=-ex
E(ω)

U=px2-rx4-epositionx

図 1.2.2 電荷-e が物質中のポテンシャル U の井戸にあり、振動電界 E･exp(-iωt)で強制振動する古典モデル。ポテンシャルは調和型の U=px2 からわずかにずれている。ポテンシャ
ルが反対称成分を含む U= px2+ qx3 と中心対称な U= px2－rx4 に対応して分極 P(t)=－ex(t)の振動にはそれぞれ２次及び３次の非線形成分が現れる。

3

物理工学科博士前期課程

H21 前期光計測工学(谷）χ ijkl (3) などと、それぞれ 3 階及び４階のテンソルで表される。2 次及び 3 次の非線形感受~ ( 2)率に対応して、 P

非線形分極と呼ぶ。

= ε 0 χ ijk

~

( 2)

~
~
( 3) ~E (t ) 2 を 2 次の非線形分極、 P ( 3) = ε 0 χ ijkl E (t ) 3 を 3 次の~(1.2.5)式では分極 P (t ) は同時刻の電界 E (t ) で表されると仮定している。これは媒質が~瞬時に電界 E (t ) に応答することを意味するが、これは媒質に吸収や分散がないと仮定していることと同等である（Kramers-Kronig の関係(注 4a)から導かれる）。実際には媒質には分散や吸収があり、非線形感受率（線形感受率も含め）は波長または周波数に依存する。図
1.2.1 に示すように非線形感受率には低周波数域では電子的な寄与とフォノン的な寄与があり、光学フォノンによる共鳴分散域より上では電子的な寄与のみが重要となる(注 4b)。電子的な共鳴、すなわち光吸収端より上の周波数では電子的な寄与も小さくなっていくが、それよりも実際には吸収による影響のほうが大きくなり非線形光学効果を利用することが困難になる。
いま、フォノンによる影響が無視できるほど十分振動数が高く、電子励起による吸収が起こらない領域、すなわち非線形感受率は電子的寄与のみを考慮すればよい周波数領域を考える(注 4c)。誘電体中の束縛電子（電荷を−e とする）は図 1.2.2 のようにあるポテンシャ
ルの平衡点にあると考える。このときある振動数の光が入射したとすると、電子は平衡点を中心に揺さぶられるが、その程度が小さいときは、ポテンシャルは U=px2（p は定数）で近似されるので、電子は強制振動を伴う調和振動子として振舞う。したがって分極 P(t)=－ex(t)は入射電場の振動数と同じ振動数で振動し、非線形な効果は生じない。光の電界強度が強くなり、x(t)の振幅も大きくなってくると、電子の感じる実際のポテンシャルは
U=px2 からずれはじめ、最初の補正項として qx3（q は定数）の項が加わるが、媒質の結晶格子が中心対称性を持つときは、ポテンシャルが左右対称なためこの項はゼロになる（q=0）。
したがって、中心対称性をもつ結晶の非線形光学効果は rx4（r は定数）の項、すなわち 3次の非線形性が一番低次数の非線形光学効果となる。
Si や Ge の単結晶や、アモルファス半導体などのガラス状物質は中心対称（均一）なので、２次の非線形性は持たない。また、気体や液体もガラス状物質同様、特定の軸方位がなく平均的に中心対称といるので、通常は 2 次の非線形性は持たない(注 4d)。一方、中心対
称でない結晶、たとえば GaAs などの化合物半導体やその他の非線形結晶は 3 次の非線形光学効果と同時に 2 次の非線形光学効果を持つ。
ここで非線形感受率の大きさを大雑把に見積もってみよう。2 次の非線形分極 ε 0 χ

( 2)

~
E2

は結晶格子中の束縛電子が感じる電界 Eat と外部電界 E が等しくなったときに 1 次の分極~

ε 0 χ (1) E と同程度になると考えられるので、 χ ( 2 ) ≅ χ (1) / E at となる。 χ (1) は１に近いので（無次元量） χ( 2)≅ 1 / E at となり E at = e /( 4πε 0 a 0 )（a0 はボーア半径、e は電子の素電荷）2から求めると
注 4a：吸収と分散(誘電率または感受率の周波数依存性)との間には因果関係があり、Kramers-Kronigの関係式を用いると吸収スペクトルから分散スペクトルが求められ、また逆も同様に求められる。
注 4b：共鳴領域では線形、非線形感受率ともに大きな値を持つが、同時に吸収も強くなるため、非線形
光学効果は利用しにくい。また、共鳴境域の低周波数側では感受率が高周波数側よりも大きくなっていることに注意しよう。このため低周波数側では共鳴領域があるたびに感受率が加算的に大きくなる。したがって静的な誘電率が大きい物質はどこかの周波数で強い共鳴吸収を持つ。

4

物理工学科博士前期課程

H21 前期光計測工学(谷）χ ( 2 ) ≈ 1 / E at ≅ 2 x 10-12 (m/V)(1.2.6)

同様に 3 次の非線形感受率は

χ (3) ≈ 1 / E at 2 ≅ 4 x 10-24 (m2/V2)(1.2.7)と見積もられる。これらのことから、3 次の非線形光学効果は 2 次の非線形光学効果より
一般に何桁も小さいことが分かる。また、非線形光学効果は光強度が非常に強くなくては観測されないことが分かる。
非線形光学では非線形分極が重要な役割を果たす。たとえば、非線形媒質中での電磁波の振る舞いはしばしば次の形の波動方程式で表される（実電流、実電荷のない場合）(注 5a)

~
~
∂2E
∂2P
~
∇ E − εμ 2 = μ 2∂t∂t2(1.2.8)~右辺の電気分極の時間についての 2 次微分 ∂ P / ∂t は電子が加速されることを意味するの22で、この項がゼロでない場合は、分極電荷が振動し、電磁波を放射することを示している。
したがって、非線形媒質中の電磁波の伝播は非線形分極によって発生した電磁波と入射電
磁波との非線形相互作用（媒質を介した）によって記述される。

１．３

各種非線形光学過程

この節では、代表的ないくつかの非線形光学過程について定性的な説明を試み、非線形光学で用いられる表式及び記号についてのルールや、非線形光学における基本概念を紹介する。
（１） 第 2 高調波発生（Second-Harmonic Generation, SHG）
まず非線形光学過程として第 2 高調波発生（SHG）、すなわち周波数ωの入射波がその倍の周波数の光に変換される過程について考えてみよう（図 1.3.1）。まず、電磁波の電界強度をつぎのように表す。

~
E (t ) = E exp(−iωt ) + c.c.

(1.3.1)

ここで c.c.は複素共役（Complex conjugate）を表す(注 5b)。この波が 2 次の非線形感受率を持つ結晶に入射すると、非線形分極は~~P ( 2) (t ) = ε 0 χ ( 2 ) E 2 (t )= 2ε 0 χ ( 2) EE * +(ε 0 χ ( 2 ) E 2 exp(−2iω t ) + c.c.)(1.3.2)

注 4c：物質にもよるが通常可視から中赤外域までの波長に対応している。紫外域では電子励起による吸収が、遠赤外域では格子振動（フォノン）による吸収が強くなる。
注 4d：中心対称性を持つ物質でも表面や外部電場・磁場がかかると、対称性が崩れ 2 次の非線形光学効果が観測される。

5

物理工学科博士前期課程

H21 前期光計測工学(谷）

(b)

(a)

ω

χ(2)

ω

ω

2ω

ω

図 1.3.1 (a) 第 2 次高調波発生の様子。

2ω(b) 第 2 次高調波発生におけるエネルギーレ
ベルで与えられる。この式から 2 次の非線形分極は周波数がゼロの DC 成分（第 1 項）と周波数が 2 ω（第 2 項）の成分からなることが分かる。2 ωの項は(2.2.8)式で示されるように、
2 ωの周波数の電磁波を発生する SHG 過程に対応している。DC の項は振動しないので、電磁波は発生しない。このような静電分極を引き起こす過程を Optical Rectification（略してOR、和訳：光整流）という。
適当な条件の下で SHG 過程は非常に効率が良く、入射した光のほとんど 100%を 2ωに変換することも可能である。代表的な例を挙げれば、Nd:YAG レーザーの基本波 1.06μmの近赤外光（目には見えない）は LBO や CLBO などといった非線形結晶を用い、0.53μmのグリーンの光に効率良く変換される。SHG は図 1.3.1(b)に図示されるように、光子変換の観点からも説明される。すなわち、=ω のエネルギーの光子が２つ消滅し、新たに 2=ω のエネルギーの光子１つが生成される量子過程として理解される( 2=ω = =ω + =ω )。図1.3.1(b)では下側の実線が原子の基底状態に対応し、点線が仮想励起状態（virtual states）に対応する。ここでいう仮想励起状態は物理的に準位が存在しないという意味ではなく、定常的な原子の固有状態ではないが、原子と光子（複数）が相互作用した瞬間だけ存在している過渡的な状態を意味している。
（２） 和 周 波 発 生 (Sum Frequency Generation, SFG) と 差 周 波 発 生 (Difference
Frequency Generation, DFG)
つぎに周波数の異なる２つの光が 2 次の非線形感受率を持つ結晶に入射する場合を考えよう。電界振幅は~E (t ) = E 1 exp(−iω 1 t ) + E 2 exp(−iω 2 t ) + c.c.

(1.3.3)で表され、結晶に誘起される非線形分極は
注 5a：(1.2.8)式の右辺がゼロのときは真空中を伝播する電磁波の波動方程式（電界成分）になっている。
~
したがって右辺の ∂ 2 P / ∂t 2 は媒質の分極による放射に対応している。
注 5b：ここで電界の振幅を表す E は phasor 表示で振動電界を E ph exp(−iω t ) と書いた場合の E

phと比べて 1/2 倍になっていることに注意。phasor 表示では物理量としての電界振幅の大きさは~Re( E ph exp(−iω t )) を意味しているので、E (t ) = Re( Ephexp(−iω t )) = 12 E ph exp(−iω t ) + c.c.

である。これと(1.3.1)式を見比べると E= Eph/2 であることが分かる。

6

物理工学科博士前期課程

H21 前期光計測工学(谷）~~P ( 2) (t ) = ε 0 χ ( 2) E 2 (t )= ε 0 χ ( 2) [ E1 exp(−2iω1 t ) + E 2 exp(−2iω 2 t )22+ 2 E1 E 2 exp(−i (ω1 + ω 2 ) t ) + 2 E1 E 2 * exp(−i (ω1 − ω 2 ) t ) + c.c.](1.3.4)+ 2ε 0 χ ( 2 ) [ E1 E1 * + E 2 E 2 *]で与えられる。これを次のように表現すると便利である。

~
P ( 2) (t ) = ∑ P(ω n ) exp(−iω n t )(1.3.5)nここで和はすべての正と負の周波数ωn についてとるものとする。ωn に対応する各周波数成分は次で与えられる。

P (2ω1 ) = ε 0 χ ( 2 ) E12 （SHG1）
P(2ω 2 ) = ε 0 χ ( 2) E 22 （SHG2）
P (ω1 + ω 2 ) = 2ε 0 χ ( 2 ) E1 E 2 （SFG）(1.3.6)P (ω1 − ω 2 ) = 2ε 0 χ ( 2 ) E1 E 2 * （DFG）
P (0) = 2ε 0 χ ( 2 ) ( E1 E1 * + E 2 E 2 *) （OR）
さらに負の周波数に対応してP (−2ω1 ) = ε 0 χ ( 2 ) E1 *2 （－SHG1），

P(−2ω 2 ) = ε 0 χ ( 2 ) E2 *2 （－SHG2）

P (−ω1 − ω 2 ) = 2ε 0 χ ( 2 ) E1 * E 2 * （－SFG）， P (−ω1 + ω 2 ) = 2ε 0 χ ( 2 ) E1 * E 2 （－DFG）(1.3.7)の項が現れるが、これらの量は(1.3.6)式のそれぞれの項の複素共役になっており、どちらか一方についてのみ考慮（我々は正の周波数を考慮）すれば十分である(注 7)。
(1.3.6)式から５つの分極成分（周波数ゼロでない成分は４つ）が得られるが、通常位相整合条件(後述)は２つの周波数成分に対して同時に満たされることがないので、実際に観測されるのはこられのうち位相整合条件を満たす１つの成分のみである。周波数を選ぶのに、実用上は入射光の偏光方法を選んだり、結晶の角度を調整したりする。

注 7：マイナスの周波数がでてきたのは、電磁波の振動電界を E exp(−iω t ) とその複素共役で表すという数学的な表現方法に起因している。マイナスの周波数を持った電磁波は物理的には意味はない
ように思えるが、直線偏光が右円偏光と左円偏光成分の重ね合わせであることから、プラスとマ
イナスの周波数はそれぞれ右円偏光と左円偏光に対応している（あるいはその逆）と考えてもよい。

7

物理工学科博士前期課程

A.

H21 前期光計測工学(谷）和周波発生（SFG）
(b)

(a)

ω1

χ(2)

ω2

ω2

ω3=ω1+ω2

図 1.3.2 (a) 和周波発生の様子。

ω1ω3(b) 和周波発生におけるエネルギーレ
ベル

図 1.3.2 に描かれているような、和周波発生（Sum Frequency Generation, SFG）、すな
わち周波数ω1とω2の２つの入射波が 2 次の非線形感受率をもつ結晶に入射し、それらの和の周波数ω3=ω1+ω2の光に変換される過程について考えてみる。その非線形分極は(1.3.6)式からP (ω1 + ω 2 ) = 2ε 0 χ ( 2 ) E1 E 2 （SFG）(1.3.8)で与えられる。SFG 過程は入射波が２つの異なる周波数を持つことを除き多くの点で SHG
過 程 と 似 て い る 。 SFG は エ ネ ル ギ ー が =ω1 と =ω 2 の ２ つ 光 子 が 消 滅 し 、 新 た に=ω 3 = =(ω1 + ω 2 ) のエネルギーの光子１つが生成される量子過程として理解される。SFG
は波長の異なる可視光を２つ重ね合わせ、波長可変な紫外光を発生させたりするのに用いられる。
B.

差周波発生（DFG）
差周波発生（Difference Frequency Generation, DFG）は次の形の非線形分極で記述される(図 1.3.3)。

P (ω1 − ω 2 ) = 2ε 0 χ ( 2 ) E1 E 2 * （DFG）(1.3.9)

この場合は発生する波の周波数は入射波ω1とω2の差ω3=ω1－ω2 (ω1>ω2)である(注 8)。DFG は
２つの異なる波長の可視光から波長可変の赤外光を発生させたりするのに用いられる。
表面的には DFG と SFG は似ているようにみえる。しかし図 1.3.3(b)にみられるように、量子過程としての DFG は SFG と様相を異にする。すなわち DFG 過程はω1の光子が消滅
し、ω2とω3の光子が新たに発生する過程である（ω3だけでなくω2の光子が作られている！）。
これは周波数の低いほうの入射波ω2が増幅されていることを意味する。このため DFG 過程
は光パラメトリック増幅（Optical Parametric Amplification, OPA）としても知られている。
注 8：周波数ω2の前のマイナス符号と振幅Ε2の複素共役が対応していることに注意。

8

物理工学科博士前期課程

H21 前期光計測工学(谷）

(b)

(a)

ω1
ω2

ω2
χ(2)

ω1

ω3=ω1−ω2

図 1.3.3 (a) 差周波発生の様子。

ω3(b) 差周波発生におけるエネルギーレベル

図 1.3.3(b)のエネルギー準位図では基底状態の原子がまずω1の光子を吸収し、高いほうの仮想励起状態へ遷移する。そして、この励起状態は入射波のω2の光子に刺激されて（Stimulated)、
ω2とω3の光子を同時に放出して基底状態へ緩和する。この2光子放出過程はω2の光子が最初なかったとしても起こる。この場合は2光子の仮想励起状態からの自然放出過程(Spontaneous two photon emission)なので、発生するω3とω2の波の強度は弱い。この過程はパラメトリック蛍光（Parametric Fluorescence）(注 9a)と知られている。
C.

光パラメトリック発振（Optical Parametric Oscillation, OPO）
DFG 過程ではω2またはω3の光子の存在が他方の光子の放出を誘発することを述べた。もし、このような過程で用いられる非線形結晶を光共振器(注 9b)の中にいれると(図 1.3.4)、ω2
またはω3の波は共振器のなかで成長・蓄積され非常に大きな強度になりうる。このような
光 共 振 器 を 用 い た 非 線 形 光 学 装 置 を 光 パ ラ メ ト リ ッ ク 発 振 器 (Optical Parametric
Oscillator)という。光パラメトリック発振器は波長可変の光源が得にくい赤外域で利用されることが多い。入射波のω1 はポンプ(Pump)周波数、得たいほうの周波数、たとえばω2 をシグナル(Signal)周波数、もう一方の必要でないほうの周波数ω3 をアイドラー(Idler、車のアイドリングと同じで、idle は“無駄な”とか“遊びの”と言う意味)周波数という。シグナルとアイドラーの区別はどちらを必要としているかで決まるが、物理的には両者は等価である。

ω1=ω2+ω3pumpω2(signal)
ω3(idler)

χ(2)図 1.3.4 光パラメトリック発振器
注 9a：ω1の光子のみを入射させ、ω1の強度を上げていくと、共振器のない状態でも、ある閾値を超えると自然放出過程の Parametric Fluorescence から誘導放射の 1 種である光パラメトリック発生(Optical Parametric Generation, OPG)になり、発生するω3とω2の波の強度は非常に強くなる。
注 9b：共振器の波長は通常シグナル、アイドラー光のどちらかに設定されることが多い。またシグナル
光の波長を変化(微調)させるためには結晶の角度を変えて、位相整合条件を変化させるなどする。

9物理工学科博士前期課程(3)

H21 前期光計測工学(谷）3 次の非線形分極（単一周波数入射の場合）

つぎの形の 3 次の非線形分極を考えることにする。

~~
P ( 3) (t ) = ε 0 χ ( 3) E 3 (t )(3 次の非線形分極)(1.3.10)~

ここで入射波 E (t ) は一般的には複数の周波数成分の重ね合わせで与えられるが、複雑になるので、単一周波数ωの電磁波が入射した簡単な場合をまず考えてみる。このとき入射波は~E (t ) = E exp(−iω t ) + c.c. = E exp(−iω t ) + E * exp(iω t )(1.3.11)と表せるので、(1.3.10)式は~~P ( 3) (t ) = ε 0 χ ( 3) E 3 (t ) = ε 0 χ ( 3) ( E exp(−iω t ) + c.c.) 3(1.3.12)2= ε 0 χ ( 3) E 3 exp(−3iω t )+ 3ε 0 χ (3) E E exp(−iω t )+ c.c.
Third Harmonic Generation

Optical Kerr Effectとなる。ここで EE*=|E|2 を用いた。この式の最初の項は 3 倍の周波数 3ωで振動する分極を表している。これは図 1.3.5 に示されるように第 3 次高調波(Third Harmonic
Generation, THG)発生を表している。図 1.3.5(b)に示されるように、これは３つの周波数ωの光子が消滅し、周波数 3ωの光子ひとつが作られる過程に対応している。
一方、第2項の分極成分は周波数が入射波と同じ（線形分極の周波数）で、|E|2E に比例している。光強度 I は|E|2 に比例するので、これは次に示すように実効的な分極率が光強度に比例して変化すること示している。分極率の変化は、すなわち誘電率、屈折率の変化であるので、第2項は周波数ωの光（入射波）が感じる媒質の屈折率が非線形光学効果（3次）によって変化することを示している(注10)。

(b)

(a)

ω

χ(3)

ω

ω

ω

3ω

3ω

ω
図 1.3.5(a) 第 3 次高調波発生の様子。

(b) 第 3 高調波発生におけるエネルギーレベル。

注 10：屈折率の虚部が光強度に比例して変化するのが 2 光子吸収過程である。したがって 2 光子吸収は
3 次の非線形光学効果である。

10

H21 前期

物理工学科博士前期課程

光計測工学(谷）3次の非線形分極による屈折率の変化は現象論的に

n = n0 + n2 I

(1.3.13)と書ける。ここで n0 は通常（非線形光学効果がないとき）の屈折率、I は光強度、ｎ2 は非線形光学効果の強さを表す係数である。ｎ2 は次のようにして求められる。
いま2次の非線形分極が無く（中心対称性の結晶の場合など）、周波数ωについての非線形分極を考えると、それは線形な分極成分と(1.3.12)式の第 2 項の和で書けるので、

P (ω ) = ε 0 χ (1) E (ω ) + 3ε 0 χ ( 3) E (ω )2

E (ω ) ≡ ε 0 χ eff E (ω )(1.3.14)ここで実効的電気感受率をχ eff = χ (1) + 3χ (3) E (ω )2

(1.3.15)で定義した。誘電率と電気感受率の関係から[

ε2= n 2 = 1 + χ eff = 1 + χ (1) + 3χ ( 3) E (ω )ε0](1.3.16)となる。そこで(1.3.13)式の 2 乗をとり、ｎ2 は非常に小さいのでｎ22I2 の項を省略して、(1.3.16)式と等しいとおくとn 2 = (n0 + n2 I ) 2 ≅ n0 + 2n0 n2 I = n0 + 4n0 n2 ε 0 / μ 0 E (ω )
2

2

2

2

= (1.3.16) = 1 + χ (1) + 3χ ( 3) E (ω ) = n0 + 3χ ( 3) E (ω )ここで光強度 I = 2n 0

ε 0 / μ0 E

22 注( 11)、

2(1.3.17)2

1 + χ (1) = ε (1) / ε 0 = n0 であることを用いた。
2

この式から

n2 =

3
4n0

2

μ 0 ( 3)χ
ε0(1.3.18)が得られる。3 次の非線形分極によるこのような屈折率の変化は光 Kerr 効果(Optical Kerr
Effect)と呼ばれている。
・Self-Focusing（自己集束）
光 Kerr 効果すなわち、光強度依存の屈折率の変化は レーザービームの Self-Focusing（自己集束）を引き起こす。すなわち(1.3.13)式の係数 n2 が正のとき、ビームの中心部分
注 11：光強度は単位面積あたりのポインティングベクトルの時間平均値で与えられる。電磁気学の教科書でよく用いられる phasor 表示の電界振幅を用いると係数は 1/4 になる。また真空中の光速度

c = 1 / ε 0 μ 0 を用いて書き直すと， I = (c / n0 ) 12 ε E
1
2

2

2

となる。これはエネルギー密度がε E の電磁波が屈折率 n0 = ε / ε 0 の媒質中を光速度 c / n0 で運ぶエネルギーに等しい。
11

物理工学科博士前期課程

H21 前期光計測工学(谷）の光強度が高いため、屈折率は中心部分で高くなり、正のレンズ(凸レンズ)として作用する。
この現象は応用上重要である。たとえば、媒質中で（光学ガラスなど）レーザービームがSelf-Focusing により絞りこまれ、強度がその媒質の閾値を超えると、ダメージを与えることになる。またチタンサファイアレーザーに代表される超短パルスレーザーでは、この
Self-Focusing を利用して光パルスが光共振器内で自動的に圧縮されるようになっている（Kerr レンズモード同期と呼ばれる手法）。
(4)

3 次の非線形分極（一般の場合）

さてここで、３次の非線形分極の一般的な場合を見てみよう。いま３次の非線形分極~~P ( 3) (t ) = ε 0 χ ( 3) E 3 (t )(3 次の非線形分極)(1.3.10)における入射波が３つの周波数成分をもつとする。すなわち~
E (t ) = E1 exp(−iω1 t ) + E 2 exp(−iω 2 t ) + E3 exp(−iω 3 t ) + c.c.

(1.3.19)~
E 3 (t ) を計算すると周波数の負の成分も含めると 44 個の周波数成分が得られる(注 12)。これ
らの周波数を具体的に書き出してみると

ω1 , ω 2 , ω 3 , 3ω1 , 3ω 2 , 3ω 3 ,(ω1 + ω 2 + ω 3 ), (ω1 + ω 2 − ω 3 ), (ω1 − ω 2 + ω 3 ), (−ω1 + ω 2 + ω 3 ),(2ω1 ± ω 2 ), (2ω1 ± ω 3 ), (2ω 2 ± ω1 ), (2ω 2 ± ω 3 ), (2ω 3 ± ω1 ), (2ω 3 ± ω 2 ),(1.3.20)

及びこれらの符号を反転させたものである。3 次の非線形分極についても 2 次の場合同様、~
P ( 3) (t ) = ∑ P(ω n ) exp(−iω n t )(1.3.21)nと書くと、それぞれの周波数ωn に対応する非線形分極は（正の周波数のものだけ書くと）P (ω 1 ) = ε 0 χ ( 3) (3E 1 E 1 * +6 E 2 E 2 * +6 E 3 E 3 *) E 1 ,P (ω 2 ) = ε 0 χ ( 3) (6 E 1 E 1 * +3E 2 E 2 * +6 E 3 E 3 *) E 2 , （Kerr Effect）(1.3.22a)
P (ω 3 ) = ε 0 χ (3) (6 E 1 E 1 * +6 E 2 E 2 * +3E 3 E 3 *) E 3 ,P(3ω1 ) = ε 0 χ ( 3) E13 , P(3ω 2 ) = ε 0 χ ( 3) E 23 , P(3ω 3 ) = ε 0 χ (3) E33 , （THG）

P(ω1 + ω 2 + ω3 ) = 6ε 0 χ (3) E1 E 2 E3 , P(ω1 + ω 2 − ω3 ) = 6ε 0 χ (3) E1 E 2 E3 *,P(ω1 − ω 2 + ω3 ) = 6ε 0 χ (3) E1 E 2 * E3 , P(−ω1 + ω 2 + ω3 ) = 6ε 0 χ (3) E1 * E 2 E3 ,(1. 3.22b)(1.3.22c)

注 12：3 つの波の重ね合わせである E(t)は複素共役のマイナス周波数の項を含め 6 つの項で表されている。E3(t)を計算すると 6x6x6=228 個の項が現れるが，同じ周波数を与える組み合わせがあり，周波数としては 44 個が現れる。このうち半分はマイナスの周波数なので，実際にはその半分の正の
周波数 22 個を考慮すればよい。

12

H21 前期

物理工学科博士前期課程

光計測工学(谷）P (2ω1 + ω 2 ) = 3ε 0 χ (3) E1 E 2 , P(2ω1 + ω 3 ) = 3ε 0 χ ( 3) E1 E3 ,22P (2ω 2 + ω1 ) = 3ε 0 χ (3) E1 E 2 , P(2ω 2 + ω 3 ) = 3ε 0 χ (3) E 2 E3 ,22

(1.3.22d)

P (2ω 3 + ω1 ) = 3ε 0 χ ( 3) E3 E1 , P(2ω 3 + ω 2 ) = 3ε 0 χ ( 3) E3 E 2 ,22

P (2ω1 − ω 2 ) = 3ε 0 χ ( 3) E1 E 2 *, P(2ω1 − ω 3 ) = 3ε 0 χ ( 3) E1 E3 *,22P (2ω 2 − ω1 ) = 3ε 0 χ ( 3) E1 * E 2 , P(2ω 2 − ω 3 ) = 3ε 0 χ ( 3) E 2 E3 *,22(1.3.22e)P (2ω 3 − ω1 ) = 3ε 0 χ ( 3) E3 E1 *, P(2ω 3 − ω 2 ) = 3ε 0 χ ( 3) E3 E 2 *,22

これらの式をよくみると次のことが見て取れる。すなわち、右辺に現れる電界振幅に対応した周波数の和が左辺の非線形分極の周波数に対応している。たとえば一番の最後のP(2ω3−ω2)では、右辺の電界振幅の積は E32E2*であるが、これは E3 が 2 回と E2 の複素共
役 E2*（したがってその位相部分である周波数は符号が反転）の積なので、対応する周波数はω3+ω3－ω2=2ω3－ω2 となっている。また右辺の項の数因子（１，３，６など）はその周波数を作る３つ周波数の可能な組み合わせ（Permutation）の数に対応している。たとえば先の P(2ω3 －ω2)例では数因子３は (E3 E3 E2*)→(ω3+ω3 －ω2), (E3 E2*E3)→(ω3 －ω2+ω3), (E2*E3 E3) →(－ω2+ω3+ω3)の３つの組み合わせに対応している。光 Kerr 効果を表す非線形感受率 χ( 3)の前の係数３は、同様に非線形分極 P(ω=ω+ω－ω)に対して、(EEE*)→(ω+ω－ω), (EE*E)→(ω－ω+ω), (E*E E) →(－ω+ω+ω)という３つの組み合わせがあることに対応している(注 13)。
図 1.3.6 に代表的な３次の非線形光学過程の様子と対応するエネルギーレベルを模式的に示した。
・

パラメトリック過程と非パラメトリック過程
以上で述べた非線形光学過程はすべてパラメトリック過程として分類される。パラメト

リック過程は光が相互作用する媒質（あるいは原子）の量子力学的な初状態と終状態が同じになるものを指す。すなわち、原子状態は非常に短い時間、仮想励起状態へ一時的に遷移し、もとの状態へ復帰する。仮想励起状態へとどまる時間 Δt は、基底状態からのエネル
ギー差をΔE とすると、不確定性原理より、 Δt = = / ΔE 程度である。パラメトリック過程
では電子やその他の素励起（フォノンなど）の実励起は伴わない。これに対し電子などの実励起をともなう過程が非パラメトリック過程である。パラメトリック過程は線形及び実数の非線形感受率で記述されるが、非パラメトリック過程は複素数(虚部を持つ)の感受率で記述される。また同時にパラメトリック過程では光子のエネルギーは保存されるが、非パ
ラメトリック過程では光子のエネルギーの一部または全部が物質系（媒質）へ散逸（または物質から流入）し、光子系のエネルギーは一般に保存されない。
注 13：これら同じ周波数を作る異なる組み合わせ(電界成分の掛け算の順序が違うだけ)に対応する各項
は物理的にはまったく等価である。これは数学的な表現が対応する物理過程に対して冗長性を持っているために，3 つの項に分かれているだけであり，後に述べる intrinsic permutation
symmetry に対応している。

13

物理工学科博士前期課程

H21 前期光計測工学(谷）(b)(a) Optical Kerr Effect

ω1
ω2
ω3

χ(2)

ω4=ω1=ω1−ω1+ω1=ω2−ω2+ω1=ω3−ω3+ω1

χ(2)ωκ (κ=1,2,3)

ω1

ω4=ω1
=ωκ−ωκ+ω1(d)(c) Sum frequency of 3 waves

ω1
ω2
ω3

ωκ

ω3

ω4=ω1+ω2+ω3

ω2

ω4

ω1
(f)

(e)

ω1
ω2
ω3

χ(2)

ω4=ω1+ω2−ω3

(g)

ω2
ω3

ω2

ω3

ω1

ω4=ω1+ω2−ω3

ω3

ω1

ω3

ω4=2ω3−ω2

(h)

χ(2)

ω4=ω3+ω3−ω2

図 1.3.6 3 次の非線形光学過程：(a),(b) 光 Kerr 効果のの様子とそのエネルギーレベル。
(c),(d) ３つの波の和周波発生の様子とそのエネルギーレベル。(e),(f) ３つの波の周波数混合（ω1+ ω2−ω3）の様子とそのエネルギーレベル。(g),(h) 2 つの波の周波数混合（2 ω2−ω3）の様子とそのエネルギーレベル。
パラメトリック過程と非パラメトリック過程の区別を示す簡単な例として、複素屈折率がある。複素屈折率の実部はパラメトリック過程(線形な)な結果であり、一方吸収を表す虚部は非パラメトリック過程、すなわち電子の実励起その他（たとえばフォノン）による光子の吸収の結果である。以下で非線形な非パラメトリック過程の例として、可飽和吸収14物理工学科博士前期課程

H21 前期光計測工学(谷）(Saturable Absorption)、光双安定(Optical Bistability)、２光子吸収、誘導ラマン散乱ついて簡単に説明する。
・ 可飽和吸収(Saturable Absorption)
非パラメトリック過程のひとつの例として可飽和吸収(Saturable Absorption)があげられる(注 15a)。多くの物質ではレーザー強度が強いと吸収係数が小さくなる。可飽和吸収が生じるときの吸収係数はしばしば次のような形で表される。

α=α0(1.3.23)

1+ I / Is

ここでα0 は光強度 I が弱いときの吸収係数、Is は吸収係数がもとの値の 1/2 に達するときの光強度で（物質に依存した定数）で飽和強度と呼ばれる。
・ 光双安定(Optical Bistability)
可飽和吸収の結果として観測される現象の一つが光双安定(Optical Bistability)である。
いま図 1.3.7 のような Fabry-Perot 光共振器のなかに可飽和吸収体が置かれている場合を考える。
入射レーザービームの強度を低いレベルから増加させていくと共振器内の強度もそれにつれて増加する（状態１）。光強度が増加すると可飽和吸収体の吸収係数が非線形に低下し、ある強度以上で急激に透過率が上がる（遷移 2）。一度共振器内の強度が上がってしまうと、過飽和吸収のため、透過率がよい状態にとどまり続ける（状態３）。しかし、入射レーザーの強度をかなり落とすと、可飽和吸収の効果が非線形に失われ、透過率が急激に低下し、もとの透過率が低い状態へ戻る（遷移４）。図 1.3.7 から分かるように、ある範囲の入射レ
ーザー強度に対して、２つの透過強度を取ることが可能である。またどちらの透過強度を取るかは、それまでのレーザー強度の履歴で決まる(注 15b)。
(b)

(a)

Iout

Iin

saturable
absorber

3

Iout

2
4
1
Iin

図 1.3.7 (a) Fabry-Perot 共振器のなかの過飽(b) 光双安定の様子。Iin の値に対して

和吸収体にレーザービームを入射する場合。

Iout はその履歴に依存して２つの値をとりうる。

注 15a：可飽和吸収( saturable absorption)は特別な現象ではなく，吸収過程そのものがもともと非線形な過程であり，基底準位のポピュレーションが励起準位よりもずっと大きな場合だけ線形であると近似できる。入射する光強度を強くしていき基底準位と励起準位ポピュレーション比が 1 対 1になったとき，吸収はなくなり物質は透明になる。なんらかの方法で励起準位のポピュレーショ
ンが基底準位のポピュレーションを超える(反転分布が生じる)と吸収係数はマイナスになる。すな
わ ち 誘 導 放 射が 起 こ る 。 これ がレーザー(Light Amplification by Stimulated Emission ofRadiation, LASER)の原理である。

15

物理工学科博士前期課程

H21 前期光計測工学(谷）・ ２光子吸収2 光子吸収過程(図 1.3.8)では原子（または媒質の分子）は基底状態から光子を 2 個同時に吸収し、上位の固有状態（励起状態）へ遷移する。2 光子吸収における吸収断面積σは入射波の強度に比例する。

ω

real excitedstatevirtual stateωground state
図 1.3.8 2 光子吸収過程
すなわち、

σ = σ ( 2) I

(1.3.24)

ここでσ(2)は 2 光子吸収を記述する係数である（線形なシステムでは吸収は光強度に依存しないことに注意）。原子の遷移確率は吸収断面積に単位時間、単位面積当たりに入射する光子数をかけたものなので、したがって、2 光子吸収による原子の遷移確率ＲはR = σ I / =ω = σ ( 2) I 2 / =ω

(1.3.25)で与えられる。したがって、２光子吸収による吸収は光強度の 2 乗に比例して強くなる。
・ 誘導ラマン散乱(Stimulated Raman Scattering)
ラマン散乱(Raman Scattering)は光が物質と相互作用し、周波数ωの光子が消滅し、周波数がωv だけずれた周波数ωS=ω− ωv の光が生成し、物質の原子または分子がエネルギー =ω vの励起状態へ遷移する過程である。このときの周波数のずれを Stokes シフトという（光が原子分子からエネルギーをもらって高い周波数へシフトする場合は Anti-Stokes シフトという）。誘導ラマン散乱(Stimulated Raman Scattering)(注 16)では散乱によって励起された周波数 ω v の分子振動(光学フォノン)が入射波と結合して、さらにラマン散乱を誘導するので、通常の自然放出的なラマン散乱より効率がかなり大きくなる。入射光に対するエネルギー
散乱効率が 10%あるいはそれ以上になる場合もある。同様な誘導散乱過程には誘導
Brillouin 散乱（Brillouin 散乱はラマン散乱の 1 種でとくに散乱が音響フォノンによるものをいう）、誘導 Rayleigh 散乱（Rayleigh 散乱は波長より非常に小さい散乱体による散乱で波長はシフトしない）などがある。
1.4.

非線形電気感受率の一般的な定義（テンソル表現）

これまでは非線形光学過程を吸収や分散のない場合の非線形分極で記述した。本節では、媒質が吸収や分散を持ち、入射波と分極をベクトルとして扱う、より一般的な記述につい注 15b：このような光双安定の性質を利用すると光フリップ・フロップなどの論理スイッチが実現できる。また光バッファメモリ(履歴の性質を利用)としての応用も検討されている。
注 16：誘導ラマン散乱は 2 光子吸収同様，3 次の非線形光学過程である。誘導ラマン散乱の場合，光強度に依存して分極率(Polarizability)が変化する。一方，2 光子吸収の場合は，吸収係数(電気感受率の虚部に対応)が変化する。

16

物理工学科博士前期課程

H21 前期光計測工学(谷）て考えよう。吸収のある一般の場合には非線形感受率はこれまでのように実数ではなく、複素振幅で表される電界と分極に対応して複素数として取り扱う必要がある。

~

まず、任意の振動電界ベクトル E(r, t) をいくつかの周波数成分の和として記述する。

~~
E(r, t) = ∑ ' E n (r, t)(1.4.1)n~

ここでプライム“ ’ ”は和は正の周波数のみについてとることを表す。 E n (r, t) は正の周波数部分と負の周波数部分の和として表すと便利である。すなわち

~
~ (+)
~ (−)
E n (r, t) = E n (r, t) + E n (r, t)(1.4.2a)

ここで

~ (+)
E n (r, t) = E n (r )exp(-iωn t)
および
~ (−)
~ (+)
E n (r, t) = E n * (r ) = E n * (r )exp(iω n t)
~ (−)

ここで E n(1.4.2b)(1.4.2c)~ (+)~~(r ) を E n (r ) の複素共役とすることで、 E n (r, t) すなわち E(r, t) が実数の物理量であることが保障されている。さらに空間的にゆっくりと変動する振幅 An(ベクトル量)を次のように定義しておく。

E n = A n exp(ik n ⋅ r )(1.4.3)

したがって、全振動電界ベクトルはこれらの振幅を用いて~E(r, t ) = ∑ ' E n exp(−iω n t ) + c.c.
n(1.4.4a)= ∑ ' A n exp{i (k n ⋅ r − ω n t )} + c.c.
n

しばしばこれらの振幅はつぎのように書くことがある。

E n = E(ω n ),A n = A(ω n )(1.4.5)

この表記法では周波数ωn は単にパラメーターであり周波数ωn で振動する関数を意味しているのではないが、これを使うと便利である。例えば実数条件(1.4.2) の最後の式はE(−ω n ) = E(ω n )*, or A(−ω n ) = A(ω n ) *と表すことができる。この表記法を用いて全振動電界ベクトルは~
E(r, t ) = ∑ E(ω n ) exp(−iω n t )n(1.4.4b)= ∑ A(ω n ) exp{i (k n ⋅ r − ω n t )}nと簡潔に表すことができる（プライムと複素共役の項が省けた）。ここで和はすべての正と負の周波数についてとる。

17

物理工学科博士前期課程

H21 前期光計測工学(谷）同様な形式で非線形分極も表すことができ~P(r, t ) = ∑ P(ω n ) exp(−iω n t )(1.4.6)nとなる。
さてここで 2 次の非線形分極テンソルを定義しておこう。2 次の非線形感受率は分極の方向と、２つの電界ベクトルの方向で決まるので、3 階のテンソル量であり、

χ ijk ( 2 ) (ω n + ω m , ω n , ω m )（2 次の非線形感受率テンソル）(1.4.7)と表記される。ここで i,j,k は直交座標系でそれぞれの方向を表す添え字、ωn とωm は入射波の周波数、ωn＋ωm は誘起される非線形分極の周波数である。２次の非線形分極はこの非線形感受率を用いて2 次の非線形分極のベクトル表現

Pi (ω n + ω m ) = ε 0 ∑ ∑ χ ijk
jk

( 2)

(ω n + ω m , ω n , ω m ) E j (ω n ) E k (ω m )(1.4.8)( nm )と表現される。ここで i j k は直交座標系での各方向成分を表す。2 番目のサンメンションの(nm)は周波数の和ωn+ ωm が同じになるような組み合わせについて和を取ることを意味する。上の式では非線形感受率 χ( 2)(ω 3 , ω 2 , ω1 ) は３つ周波数をパラメーターとしてもつが、最初の周波数ω3 は必ず他の２つの周波数の和になっている（ω3 =ω2+ω1）。このため最初の周波数パラメーターを他の２つと区別して、

χ ijk ( 2 ) (ω 3 ; ω 2 , ω1 ) や χ ijk ( 2 ) (ω 3 = ω 2 , ω1 ) ←(他の表記法)(1.4.9)などと書いたりする場合もある。
次にテンソル表現された一般の電気感受率で記述される非線形光学過程についていくつかの例を挙げて検討してみる。
(1)

和周波発生 SFG

２つの入射波の周波数がω1 とω2 である和周波発生(SFG)の場合についてみてみよう。和周波数はω3 =ω1+ω2 であるから、(1.4.8)式より SFG の非線形分極（の各ベクトル成分）は18

H21 前期

物理工学科博士前期課程

Pi (ω 3 ) = ε 0 ∑ [ χ ijk( 2)

光計測工学(谷）(ω 3 , ω1 , ω 2 ) E j (ω1 ) E k (ω 2 ) + χ ijk( 2)(ω 3 , ω 2 , ω1 ) E j (ω 2 ) E k (ω1 )]jk(1.4.10)で与えられる。後により詳しく議論するが非線形感受率の対称性から

χ ijk ( 2 ) (ω m + ω n , ω m , ω n ) = χ ikj ( 2) (ω m + ω n , ω n , ω m )(1.4.11)であるので、(1.4.10)はPi (ω 3 ) = 2ε 0 ∑ χ ijk( 2)(ω 3 , ω1 , ω 2 ) E j (ω1 ) E k (ω 2 )(1.4.12)jkと書ける。例えば２つの入射波の偏光がともに x 軸方向を向いている場合では

Pi (ω 3 ) = 2ε 0 χ ixx (ω 3 , ω1 , ω 2 ) E x (ω1 ) E x (ω 2 )( 2)(1.4.13)となる。
(2)第 2 高調波発生 SHG

入射波の周波数がω1 でω3 =2ω1 の周波数の非線形分極（光）が発生する第 2 高調波発生(SHG)の場合を考えよう。やはり(1.4.8)式を用いて、Pi (ω 3 = 2ω 1 ) = ε 0 ∑ χ ijk( 2)(ω 3 , ω 1 , ω 1 ) E j (ω 1 ) E k (ω 1 )(1.4.14)jkが得られる（SFG の場合のように係数に２がついていないことに注意）。入射波の x 軸方向を向いている特殊な場合では

Pi (ω 3 ) = ε 0 χ ixx( 2)(ω 3 , ω1 , ω1 ) E x (ω1 ) 2(1.4.15)となる。ここで(1.4.13)と比べると係数の２が付いていないが、これは SHG の場合、入射波が１つで SFG の場合と比べて、入射波の全強度が小さくなっていることに対応している。
一般に(1.4.8)式における周波数に関するサンメンション（∑ ）を実行すると(nm )Pi (ω n + ω m ) = ε 0 D ∑ χ ijk( 2)(ω n + ω m , ω n , ω m ) E j (ω n ) E k (ω m )jk

19(1.4.16)

物理工学科博士前期課程

H21 前期光計測工学(谷）という形が得られる。ここで現れる数係数 D は縮退因子(Degeneracy Factor)と呼ばれる。
縮退因子は同じ和周波数を与える入射周波数ωn とωm の（交換可能な）組み合わせの数である。2 次の非線形分極(ベクトル)を与える(1.4.8)式は、より高次の非線形分極の場合にも一般化できる。3 次の非線形分極の場合は3 次の非線形分極のベクトル表現

Pi (ω o + ω n + ω m ) = ε 0 ∑
jkl

∑χ

( 3)
ijkl

(ω o + ω n + ω m , ω o , ω n , ω m )( nmo )

× E j (ω o ) E k (ω n ) El (ω m )(1.4.17)= ε 0 D ∑ χ ijkl (ω o + ω n + ω m , ω o , ω n , ω m )E j (ω o ) E k (ω n ) El (ω m )
( 3)

jkl

ここで D =

∑ は入射周波数ω ,ω 及びω に関する縮退因子である
onm(注 20)。

(nm )

注 20：例えば SFG の場合の(1.4.12)式では右辺第 1 項と第 2 項は周波数成分 E(ω1)と E(ω2)の掛け算の順序が違っているだけで同じ周波数の SFG に寄与する。したがって縮退因子 D == 2 とな∑( nm )

る。縮退因子は先に述べた数学的表現の冗長性のために現れる。

20

物理工学科博士前期課程

１．５

H21 前期光計測工学(谷）古典的非調和振動子の非線形感受率

Lorentz モデル、すなわち原子を調和振動子として扱うモデルは、原子分子気体と非金属固体の光学線形応答を非常によく説明できることが知られている。この節では、このLorentz モデルを電子の強制振動において復元力に非線形応答が含まれる場合に拡張する。
媒質が中心対称性を持つか持たないかで、まったくその応答特性が異なってくる。我々はまず、中心対称性を持たない媒質の場合を検討するが、その場合には 2 次の非線形応答が現れることが示される。次に中心対称性を持つ媒質の場合を検討し、その場合には 2 次の非線形応答は現れず、最低次数の非線形応答は 3 次の非線形応答であることが示される。
古典モデルでは個々の原子または分子は 1 個の共鳴振動数しか持たないが、量子論的扱い
では一個の原子は複数の共鳴振動数を同時に持っている。したがって、古典論ではある種の非線形感受率の共鳴特性（例えば 1 光子もしくは 2 光子の同時共鳴など）をうまく説明できず、より一般的には量子論的なモデルを用いねばならない。しかし、古典的モデルは直感的で分かりやすいという利点があり、また電子的共鳴周波数(図 1.2.1)より十分低い周波数領域では、古典的モデルは非常にうまく光学応答を記述することが可能である。
（１） 中心対称性も持たない媒質の場合
いま電子(電荷−e、質量 m)が図 1.2.2 のように、ある物質中のポテンシャルの井戸 U の平衡点にあるとしよう。調和型ポテンシャルからのずれが小さいとして、x 軸方向へのテーラー展開を行い、U を x の 3 次の項までで近似すると、

U = px 2 + qx 3

(1.5.1)と書ける(p>0 と q はそれぞれ定数)。物質の結晶配置が中心対称性を持つ場合は x の奇数次の qx3 の項はゼロになるが、中心対称性を持たない場合はゼロではない。このポテンシャ
ルにおける復元力はFrestore = −∂U= −2 px − 3qx 2∂x(1.5.2)で与えられる。また速度 dx/dt に比例したダンピング力（または摩擦力）をFdamp = −Γ∂x
≡ −Γx∂t（Γはダンピング係数）(1.5.3)とする（マイナスの符号は力が速度の向きと逆であることを示す）。このとき周波数ωの電磁波（偏光方向は x 軸に平行）が入射したとすると、電子はこのポテンシャルのもとで強制振動を行い、その振動の運動方程式は慣性力 Finert = mx と外力＝｛ダンピング Fdamp ＋復21

物理工学科博士前期課程

H21 前期光計測工学(谷）~~元力 Frestore ＋外部電場 E (t ) によるクーロン力（ − eE (t ) ）｝を釣り合わせた形で書けるので~~Finert = Fdamp + Frestore + (−eE ) Æ Finert − Fdamp − Frestore = −eE~Γ ~ 2 p ~ 3q ~ 2
eE (t )
2~
2
~
~
~
~





Æ x+
x+
x+
x = x + 2γx + ω0 x + ax = −mmmminertia damping linear
nonlinear
restoring

restoring(1.5.4)

external
field~で与えられる（x は速く振動する E (t ) によって、やはり速く振動するので頭にチルダをつけた）。ここで、ω02=2p/m、a=3q/m である。またダンピング係数をΓ=2γm で表したが、2γは固有振動（ω=ω0)の線幅（半値全幅＝FWHM）に対応している。
入射波の電界として次の形で与えられるものを仮定しよう。

~
E (t ) = E1 exp(−iω1 t ) + E 2 exp(−iω 2 t ) + c.c.

(1.5.5)このような入射波に対する運動方程式(1.5.4)の一般的な解は知られていないが、入射波のx に対して非線形項 a~x
強度が十分弱く、それによって引き起こされる変位 ~x が線形項 ω 0 ~22~に比べて非常に小さければ摂動論の手法で解くことができる。すなわち、まず E (t ) を~

λE (t ) で置き換える。ここで係数λは摂動の強さを表し、ゼロから１の間で連続な値をとる
が、計算の最後には１とおく。この置き換えで(1.5.4)式は~
λeE (t )2~~
x + 2γ~
x + ω 2 ~0 x + ax = −m(1.5.6)~x = λ~
x (1) + λ2 ~x ( 2 ) + λ3 ~x ( 3) + ⋅ ⋅ ⋅(1.5.7)x をλで展開する。
と書かれる。そして摂動論の手続きに従い、 ~x が、λのあらゆる値で方程式(1.5.6)の解であるためには、λについて同じ次数の
この形の ~
項同士が独立に方程式を満たす必要がある。したがってλの3次までの項については、それ
ぞれ次の式が成立する((1.5.7)を(1.5.6)に代入してλの各次数について整理）。

~
eE (t )
2 ~ (1)
(1)
(1)
~
~


x + 2γx + ω 0 x = −
m
( 2)
~
~
x ( 2 ) + 2γ~
x ( 2 ) + ω 2 ~+ a[ x (1) ] 2 = 00 x(λについて１次の項をまとめた)(1.5.8a)(λについて 2 次の項をまとめた)(1.5.8b)2 ( 3)~x ( 3) + 2γ~
x (3) + ω 0 ~x + 2a~x (1) ~x ( 2) = 0(λについて 3 次の項をまとめた)(1.5.8c)~(t ) についての方程式(1.5.8a)は通常の調和振動子（その強制力が − eE (t ) ）のx (1) (t ) は線形光学応答に対応している。
ものと同じになっていることがわかる。すなわち ~最低次の ~x(1)

その定常解は22

H21 前期

物理工学科博士前期課程

光計測工学(谷）~x (1) (t ) = x (1) (ω1 ) exp(−iω1 t ) + x (1) (ω 2 ) exp(−iω 2 t ) + c.c.
で与えられる。ここで各周波数に対する振幅 x

x (1) (ω j ) = −(1)(1.5.9)(ω j ) は

Ej
e
e Ej
≡−
2
2
m (ω0 − ω j − 2iω j γ )m D(ω j )(1.5.10)の形で与えられる(注 23)。ここで、分母に現れる複素関数をD(ω ) = ω 0 − ω 2 − 2iωγ
2(1.5.11)と定義した。
この ~x(1)(t ) についての解を次の次数の方程式(1.5.8b)に代入して解くことにより、( 2)~
x (t ) の 解 を 得 る こ と が で き る 。 ~x ( 2 ) (t ) は ２ つ の 入 射 周 波 数 に 対 応 し て 、± 2ω1 , ± 2ω 2 , ± (ω1 + ω 2 ), ± (ω1 − ω 2 ), および 0 の周波数成分を含むので、それぞれの周
波数について方程式(1.5.8b)を解く。各周波数についての解が得られれば、一般解はそれらの重ね合わせで書ける。例えば、周波数 2ω1 での応答を求めるために次の方程式を解く。
2⎛ eE ⎞ exp(−2iω1t )2 ( 2)
~x ( 2 ) + 2γ~
x ( 2 ) + ω 0 ~x = − a⎜ 1 ⎟
D 2 (ω1 )⎝ m ⎠(1.5.12)

これは強制力の周波数が 2ω1 でその振幅が－a(eE1/m)2/D2(ω1)で与えられる調和振動子の運動方程式に相当している。この方程式に対する定常解として~
x ( 2 ) (t ) = x ( 2) (2ω1 ) exp(−2iω1t )(1.5.13)の形のものを考え、これを(1.5.12)式に代入すると
2

2

E1
⎛e⎞
x (2ω1 ) = −a⎜ ⎟2⎝ m ⎠ D(2ω1 ) D (ω1 )( 2)（SHG1）(1.5.14a)（SHG2）(1.5.14b)（SFG）(1.5.14c)が得られる。同様に他の周波数についても計算すると
2

2

E2
⎛e⎞
x (2ω 2 ) = −a⎜ ⎟2⎝ m ⎠ D(2ω 2 ) D (ω 2 )( 2)

2

E1 E 2⎛e⎞x ( 2) (ω1 + ω 2 ) = −2a⎜ ⎟⎝ m ⎠ D(ω1 + ω 2 ) D(ω1 ) D(ω 2 )

注 23：実際には(1.5.9)式の解を仮定し，(1.5.8a)式に代入して振幅 x(1)を求める。

23

H21 前期

物理工学科博士前期課程

光計測工学(谷）2

E1 E 2 *⎛e⎞x (ω1 − ω 2 ) = −2a⎜ ⎟⎝ m ⎠ D(ω1 − ω 2 ) D(ω1 ) D(−ω 2 )
( 2)

(DFG)

(1.5.14d)2⎞
E1 E1 *
E2 E2 *⎛e⎞ ⎛⎟⎟+x (0) = −2a⎜ ⎟ ⎜⎜⎝ m ⎠ ⎝ D(0) D(ω1 ) D(−ω1 ) D(0) D(ω 2 ) D(−ω 2 ) ⎠( 2)（OR） (1.5.14e)

これらを次のようにまとめて表すことができる。

Ei (ωi ) E j (ω j )⎛e⎞x (ω k = ωi + ω j ) = − ∑ a⎜ ⎟i , j:ωk =ωi +ω j ⎝ m ⎠ D (ω k ) D (ω i ) D (ω j )2( 2)(1.5.14f)

こ こで ωi , ω j = ±ω1 ,±ω 2 で ある 。また Ei * (ωi ) = Ei ( −ωi ) で あ る。サ ンメンションはω k = ωi + ω j となるすべての ωi と ω j の組み合わせについてとる(注 24)。
(1)
( 2)
つぎにこれらの振動振幅に対する結果を線形感受率 χ と 2 次の非線形感受率 χ で表してみよう。線形感受率は線形分極と電界振幅の関係から

P (1) (ω i ) = ε 0 χ (1) (ω i ) E (ω i )(1.5.15)と定義されていた。いま考えている古典的振動子モデルでは、原子の数密度を N とすると
線形分極は

⎛ e E (ω i ) ⎞ e 2 N E (ω i )⎟⎟ =P (1) (ω i ) = − Nex (1) (ω i ) = (− Ne)⎜⎜ −m D(ω i )⎝ m D(ω i ) ⎠(1.5.16)

ε'
ε''

2γ
ε', ε''ε0

ω0Frequency
図 1.5.1 振動子モデルによる線形の感受率 x(1)(ω ) の周波数依存特性。

注 24：たとえば，ω k = ω1 − ω 2 = −ω2 + ω1 のとき組み合わせは (ω1 ,−ω 2 ) と (−ω2 , ω1 ) であり，これらは周波数の組み合わせ順序を変えただけで同じものであるので係数の２を得る。また
E2 (−ω 2 ) = E2 * (ω 2 ) であるので(1.5.14d)式の右辺の形の式を得る。

24

H21 前期

物理工学科博士前期課程

光計測工学(谷）で与えられるので(1.5.15)式との比較から、線形の感受率 χ(1)はe2 N 1e2 N
1χ (ωi ) ==2
ε 0 m D(ωi ) ε 0 m (ω0 − ωi 2 − 2iγωi )(1)(1.5.17)で与えられる。
誘電率と線形の感受率 χ(1)はε ≡ ε '+iε ' ' = ε 0 (1 + χ (1) (ω )) = ε 0 (1 + χ '(1) (ω ) + iχ ' '(1) (ω ))(1.5.18a)の関係があるので、誘電率の実部ε’と虚部ε’’はそれぞれ、

ε ' = ε 0 (1 + x ′ (1) (ω )) = ε 0 +ω0 − ω 2e2 N,
m (ω 0 2 − ω 2 ) 2 + (2γω ) 2
2(1.5.18b)e2 N2γω
ε ' ' = ε 0 x' ' (1) (ω ) =2m (ω 0 − ω 2 ) 2 + (2γω ) 2と表すことができる。この式は変位分極の場合の誘電分散を良く表している（配向分極の場合は別の形の分散を示す）。(1.5.18b)で表される、複素誘電率の固有振動数ω0 付近での周波数依存特性を図 1.5.1 に示す(注 25)。
2 次の非線形感受率 χ( 2)

も線形の感受率 χ(1)と同様な方法で計算される。例えば、SHG

に対する非線形分極はP ( 2 ) (2ω1 ) = ε 0 χ ( 2) (2ω1 , ω1 , ω1 ) E (ω1 ) 2で与えられ（ここで P( 2)(1.5.19)(2ω1 ) は 2ω1 で振動する非線形分極成分の振幅を表す）、古典的振動子モデルでは対応する分極は(1.5.14a)を用いて

P

( 2)

⎛ ⎛ e ⎞2⎞E1 (ω1 ) 2⎟(2ω1 ) = − Nex (2ω1 ) = (− Ne)⎜ − a⎜ ⎟⎜ ⎝ m ⎠ D (2ω ) D 2 (ω ) ⎟
1
1 ⎠
⎝
( 2)

(1.5.20)⎛ e3 ⎞
E1 (ω1 ) 2= aN ⎜⎜ 2 ⎟⎟2
⎝ m ⎠ D (2ω1 ) D (ω1 )で与えられるので、SHG に対する非線形感受率

χ ( 2) (2ω1 , ω1 , ω1 ) =aN e 3
1
2
ε 0 m D(2ω1 ) D 2 (ω1 )(SHG)を得る。この式は(1.5.17)の関係を用いて線形感受率 χχ ( 2 ) (2ω1 , ω1 , ω1 ) =ε0 a m2

N 2 e3

χ (1) (2ω1 )[ χ (1) (ω1 )] 2

(1)(1.5.21a)で表すことができて、(SHG)(1.5.21b)

注 25：Lorentz モデルにおいて(1.5.18b)式で表される誘電率虚部(図 1.5.1 の破線)は共鳴周波数ω0 を中心とする，いわゆる Lorentz 型をしている。また図 1.5.1 に示されているように誘電率実部(実線)は共鳴周波数ω0 付近でマイナスの値を持つこともある。

25

物理工学科博士前期課程

H21 前期光計測工学(谷）となる。周波数ω2に対する SHG の非線形感受率はω1 →ω2と置き換えるだけでよい。その他の周波数に対する非線形感受率も同様な手続きで得られる。

χ ( 2 ) (ω1 + ω 2 , ω1 , ω 2 ) =aN e 3
1
2
ε 0 m D(ω1 + ω 2 ) D (ω1 ) D(ω 2 )ε 02a m(1)=

ε 02a m(1)=

ε 02a mχ (ω1 + ω 2 ) χ (ω1 ) χ (ω 2 )N 2 e3aN e 31
χ ( 2 ) (ω1 − ω 2 , ω1 ,−ω 2 ) =2
ε 0 m D(ω1 − ω 2 ) D(ω1 ) D(−ω 2 )=(1)

χ (ω1 − ω 2 ) χ (ω1 ) χ (−ω 2 )N 2 e31aN e 3
χ ( 2 ) (ω1 − ω1 , ω1 ,−ω1 ) =2
ε 0 m D(ω1 − ω1 ) D(ω1 ) D(−ω1 )N 2 e3(SFG)(1.5.22)(DFG)(1.5.23)(OR)(1.5.24)

(1)

(1)

(1)

χ (ω1 − ω1 ) χ (ω1 ) χ (−ω1 )
(1)

(1)

(1)

以上のことから、中心対称性を持たない媒質では非線形光学効果の最低次数は 2 次の非線形応答(分極)であることが分かる（電界振幅の 2 乗に比例）。
より高次の非線形応答についても 3 次の摂動についての式(1.5.8c)、およびより高い次数のλn についての同様な方程式を順次解いていくことで χ

(n)

を得る。

・Miller の法則
Miller は 2 次の非線形感受率と線形感受率の間の次のような比が中心対称性を持たない物質ではほぼ一定になることを経験則として見出した（Miller の法則）。

χ ( 2) (ω1 + ω 2 , ω1 , ω 2 )
χ (1) (ω1 + ω 2 ) χ (1) (ω1 ) χ (1) (ω 2 )(1.5.25)

先に導出した関係式によれば、定数

ε 0 2 ma

(1.5.26)

N 2e3が物質によらずほぼ一定であれば、Miller の法則が成立しているのが分かる。実際、固体の原子の数密度 N は物質によらずほぼ一定で（~1022 cm-3）、真空の誘電率ε0、電子の質量m 及び電荷 e は基本定数である（注 26）。ここで非線形定数 a（我々のモデルではポテンシャ
ルの調和関数からのずれとして導入した x3 項の係数に対応していた）は物質の結晶構造、
格子定数、電子状態のバンド構造などで決まるが、物質によってそれほど大きな違いはな
いとすれば、(1.5.25)はほぼ物質によらずほぼ一定となり、Miller の法則が成り立つ。
Miller の法則は言い換えると、各周波数(ω1+ω2, ω1,ω2)での線形感受率あるいは線形誘電
注 26：固体中では有効質量 m*、有効電荷 e*、でしばしば置き換えられる。

26

物理工学科博士前期課程

率ε

(1)

H21 前期光計測工学(谷）= ε 0 ( χ (1) + 1) が分かれば、2 次（およびより高次）の非線形感受率はある程度推定できることを意味している。
（２） 中心対称性を持つ媒質の場合
つぎに中心対称性を持つ媒質の場合の古典振動子モデルについて考えよう。
やはり図 1.2.2 のように、電子がある物質中のポテンシャルの井戸 U の平衡点にあるとしよう。同様に調和型ポテンシャルからのずれが小さいとして、x 軸方向へのテーラー展開を行うが、中心対称性をもつ媒質中では U は x の偶関数で表されるので３次の項はゼロになり、ポテンシャルを x についての 4 次の項まで含めて近似する。したがってポテンシャル
は、

U = px 2 − rx 4

(1.5.27)と書ける(p と r はそれぞれ正の定数)。ここで非線形定数 r は実際の物質では通常プラスになる(注 27)。このポテンシャルにおける復元力はFrestore = −∂U= −2 px + 4rx 3∂x(1.5.28)で与えられる。したがって(1.5.4)式は~
eE (t )
Γ  2 p ~ 4r ~ 32~x + 2γ~x−x = ~x + ω 0 ~x + b~
x3 = −x + ~
x+mm
m
m
nonlinearinertia damping linear
restoring

restoring(1.5.29)

external
fieldと書き直すことができる。ここで、b = 4r/m である。この節の後半ではっきり示されるように、この場合の応答、すなわち中心対称性をもつ物質中での振動子の最低次数の非線形応答（分極）は 3 次の電気感受率 χ( 3)で記述される。中心対称を持たない媒質の場合もそうであるが、非線形感受率のテンソル特性は媒質固有の対称性(結晶の対称性)が分かっていな
いと記述できない。中心対称性を持つ媒質でもっとも重要な例の一つは、ガラス状物質や液体などの、均一物質の場合である。以下に均一物質の場合の 3 次の非線形感受率のテンソル特性について検討してみよう。
均一媒質の場合の復元力のベクトル表現は次のように書ける((1.5.2)に対応）。

~
2
Frestore = −mω 0 ~r + mb(~
r ⋅~
r )~
r

(1.5.30)r は平衡点からの変位ベクトルである。したがって、電子の運動方程式は
ここで ~~eE(t )2~
r + 2γ ~r + ω 0 ~r − b(~
r ⋅~
r )~
r =−
m

(1.5.31)と表せる。ここで入射波の電界（ベクトル）を次のように３つの周波数成分の和として与える。
注 27：結晶中では原子あるいは分子の規則正しい配列のなかで，電子は周期ポテンシャルを感じている。
このポテンシャルは原子・分子どうしのちょうど中間あたりで極大をとる（無限に大きくはならない）とすれば，係数 r は正である必要がある( 図 1.2.2 参照)。

27

物理工学科博士前期課程

H21 前期光計測工学(谷）~E(t ) = E1 exp(−iω1t ) + E 2 exp(−iω 2 t ) + E 3 exp(−iω 3t ) + c.c.

(1.5.32)

この表記は 3 次の非線形応答を考える場合の最も一般的な形であるが、これをそのまま用いて計算すると非常に表式が煩雑になる。そこで、これを次のように表記することにする。

~
E(t ) = ∑ E(ω n ) exp(−iω n t ).

(1.5.33)n

ここで、周波数についてのサンメンション（∑ ）は複素共役な項（マイナス周波数に対n応）を含めた 3ｘ２＝６個の周波数について取る。この入射波に対して、中心対称性を持たない媒質について行ったのと同様のやり方で、電子の運動方程式を摂動論を用いて解く。

~~

すなわち、 E(t ) を λ E (t ) で置き換え、 ~r (t ) をλ について展開し、λ, λ2, λ3,･･･のそれぞれの次数についての次の方程式を解く。

~
eE(t )~~~~
r + 2γ ~
r + ω 2 ~0 r − b( r ⋅ r ) r = −λm(1.5.34)に~r (t ) = λ ~r (1) (t ) + λ2 ~r ( 2 ) (t ) + λ3 ~r ( 3) (t ) + "(1.5.35)を代入し、λ, λ2, 及びλ3の各項ごとに方程式を作ると~
eE(t )2 (1)~r (1) + 2γ ~r (1) + ω 0 ~r =−m(1.5.36a)2 ( 2)~r ( 2 ) + 2γ ~r ( 2) + ω 0 ~r =0(1.5.36b)

2 ( 3)~r ( 3) + 2γ ~r (3) + ω 0 ~r − b( ~r (1) ⋅ ~r (1) )~
r (1) = 0

(1.5.36c)(1.5.36a)式は良く知られた調和振動子の運動方程式のベクトル表現となっている。その定常解は~r (1) (t ) = ∑ r (1) (ω n ) exp(−iω n t )(1.5.37)nとしたとき

eE(ω n ) / m~
r (1) (ω n ) = −D (ω n )(1.5.38)で与えられる。ここで D (ω n ) は以前と同様、 D (ω n ) = ω 0 − ω n − 2iω n γ で定義される。
2

したがって周波数ωn の線形な分極成分は

28

2

物理工学科博士前期課程

H21 前期

P (1) (ω n ) = − Ner (1) (ω n ) =Ne 2 E (ω n ) / m
D(ω n )

光計測工学(谷）(1.5.39a)で与えられ、一方、線形な感受率テンソルは直交座標系における線形分極の成分と電界成分の関係からPi (ω n ) = ε 0 ∑ χ ij (ω n ) E j (ω n ), χ ij (ω n ) = χ (1) (ω n )δ ij
(1)

(1)

(1)

(1.5.39b)jで定義される（均一媒質の場合）ので、(1.5.39a)と(1.5.39b)を比べることにより、線形な
感受率は

χ (1) (ω n ) =Ne 2 / m
ε 0 D(ω n )(線形な感受率)(1.5.39c)で与えられる。
つぎに 2 次の非線形応答であるが、対応する方程式の(1.5.36b)の右辺がゼロなので、強r
制力が働いていない。ダンピングの項 2γ ~( 2)が左辺にあるので、この場合の定常解はゼロになる。

~
r ( 2) = 0(1.5.40)

したがって、2 次の非線形感受率もゼロになる。これは媒質が均質（中心対称）でそのポテンシャルの形が偶関数であることの結果である。

r
3 次の非線形応答を計算するには(1.5.38)で得られた ~(1)(ω n ) を(1.5.36c)式へ代入する。

すると次の式が得られる。

be [E(ω m ) ⋅ E(ω n )]E(ω p )2 ( 3)
~r ( 3) + 2γ ~r ( 3) + ω 0 ~r = −∑ 3
mnp m D (ω m ) D (ω n ) D (ω p )3(1.5.41)× exp{−i (ω m + ω n + ω p )t}m,n,p についてのサンメンションが右辺に現れるが、このためにこの式は３つの周波数± ω1 ,±ω 2 ,±ω 3 の組み合わせによる多くの周波数成分を含んでいる。したがって、これらの周波数成分について個々に調和振動の方程式を解けばよい。一般解はそれらの解の組み合わせで与えられる。そこでこれらの周波数のうちのひとつを ω q = ω m + ω n + ω p と表すと、(1.5.41)の解は~
r ( 3) = ∑ r ( 3) (ω q ) exp(−iω q t )(1.5.42)qという形で表すことができる。(1.5.42)式を(1.5.41)式に代入すると周波数ωq についての解r ( 3) (ω q ) は次で与えられる。

(− ω q − iω q 2γ + ω 0 )r ( 3) (ω q ) = − ∑22( mnp )

be 3 [E(ω m ) ⋅ E(ω n )]E(ω p )m 3 D(ω m ) D(ω n ) D(ω p )29(1.5.43)

H21 前期

物理工学科博士前期課程

光計測工学(谷）あるいは D (ω q ) = −ω q − iω q 2γ + ω 0 -なので22

be3 [E(ωm ) ⋅ E(ωn )]E(ω p )r (ωq ) = − ∑( 3)(1.5.44)m 3 D(ωq ) D(ωm ) D(ωn ) D(ω p )( mnp)

ここでサンメンション∑ は ω = ω + ω + ω となるすべての組み合わせについての
q

m

n

p

( mnp )

和を取ることを意味している。さて、古典モデルにおいて周波数ωq で振動する 3 次の非線
形分極の成分は P

P

( 3)

( 3)

(ω q ) = − Ne r

(ω q ) = − Ner ( 3) (ω q ) で与えられるので( 3)(ω q ) = ∑( mnp )

= ∑Nbe 4 [E (ω m ) ⋅ E (ω n )]E (ω p )m 3 D (ω q ) D (ω m ) D (ω n ) D (ω p )Nbe 4 [ E x (ω m ) E x (ω n ) + E y (ω m ) E y (ω n ) + E z (ω m ) E z (ω n )]E(ω p )m 3 D (ω q ) D (ω m ) D (ω n ) D (ω p )( mnp )

= ∑∑Nbe 4 [ E j (ω m ) E k (ω n )δ jk ]E(ω p )m 3 D (ω q ) D (ω m ) D (ω n ) D (ω p )( mnp ) jk

(1.5.45a)で表される(注 30)。成分表示にすると

Pi (ω q ) = ∑

∑

( 3)

Nbe 4 [ E j (ω m ) E k (ω n )δ jk ]El (ω p )δ il(1.5.45b)m 3 D(ω q ) D(ω m ) D(ω n ) D(ω p )( mnp )

jklと表現できる。一方 3 次の非線形感受率テンソルの定義は

Pi (ω q ) = ε 0 ∑∑ χ (ω ,ω ,ω ,ω )E (ω ) E (ω ) E (ω )

( 3)

( 3)
ijkl

q

m

n

p

j

m

k

n

l

p

(1.5.46)( mnp )

jkl

で与えられる。この式はサンメンションのなかに m, n,及び p というダミー変数（添え字1,2,3 に対する変数）を含んでいる。したがって、非線形感受率の表現の仕方に任意性がある。例えば、(1.5.45b)と(1.5.46)を比較して、

χ (ω q , ω mω n , ω p ) =
( 3)
ijkl

Nbe 4δ jk δ il

ε 0 m 3 D(ω q ) D(ω m ) D(ω n ) D(ω p )(1.5.47)と書くことができる。しかし、ここではどれを Ej(ωm)とし、どれを Ek(ωn)とし、どれを El(ωp)とするかに関して任意性がある。すなわち Ej(ωm)、Ek(ωn)、El(ωp)を並べ替えて同じものを表現できる。この順序のとり方には、6 種類あり、それをクロネッカーのデルタの組み合わせで書くと[δ ij δ kl , δ kl δ ij ], [δ ik δ jl , δ jl δ ik ], [δ il δ jk , δ jk δ il ]注 30： クロ ネ ッ カ ー の デ ル タ は δ ij = 0(i ≠ j ) , δ ij = 1 (i = j ) で定義される。したがって∑ E E δ = E E + E E + E E などとなる。
i

j

ij

1

1

2

2

3

(1.5.48)

3

ij

30

物理工学科博士前期課程

H21 前期光計測工学(谷）である((1.5.47)式は最後の組み合わせに対応)。括弧[]で囲んだ組み合わせはクロネッカ

ーのデルタの順番を入れ替えただけなので同じものである。これらの組み合わせについての対称性（固有置換対称性、intrinsic permutation symmetry と呼ばれる）を表現するために、これら 6 つの各表現に 1/6 をかけてすべて足し合わせた、次のような表現を用いることにする。

χ (ω q , ω mω n , ω p ) =
( 3)
ijkl

Nbe 4 [δ ij δ kl + δ ik δ jl + δ il δ jk ]3ε 0 m 3 D(ω q ) D (ω m ) D (ω n ) D (ω p )(1.5.49)

先に求めた(1.5.39c)による線形の感受率を用いて各周波数に対応する D(ω)を消去すると( 3)(ω q , ω mω n , ω p ) =χ ijklε 0 3bm[ χ (1) (ω q ) χ (1) (ω m ) χ (1) (ω n ) χ (1) (ω p )][δ ij δ kl + δ ik δ jl + δ il δ jk ]3N 3e 4と表すことができる(注 31)。
(1.5.50)

注 31 ： こ の 式 か ら 均 一 媒 質 中 の χ( 3)

χ (3)1233 な ど は ゼ ロ に な る こ と が 分 か る 。 一 方 ，
χ (3)1212 = χ (3) 2323 = χ (3) 3131 ， χ (3)1221 = χ (3) 23332

1234 やχ (3)1122 = χ (3) 2233 = χ (3) 3311 ，= χ ( 3) 3113 および χ ( 3)1111 などはゼロでない値を持つ。またら均一媒質中では方向に区別はない
( 3)
( 3)
( 3)
ので χ 1111 = χ 2222 = χ 3333 である。これらのことから液体など等方的，均質な媒質中では独立な 3 次の非線形感受率は４つしかないことが分かる。

31

物理工学科博士前期課程

２．６

H21 前期光計測工学(谷）非線形感受率の対称性

この節では非線形感受率の対称性について検討してみよう。例えば３つの周波数ω1、ω2 、
ω3=ω1+ω2 の波が 2 次の非線形分極を通して相互作用する場合、非線形分極 P(ωl=ωn+ω m)は一般に次の形で表される。

Pi (ω n + ω m ) = ε 0 ∑ ∑ χ ijk( 2) (ω n + ω m , ω n , ω m )E j (ω n ) Ek (ω m )jk(1.6.1)( nm )

したがって、3 つの波の相互作用を記述するには次の 12 個の周波数の非線形感受率テンソルを求める必要がある。

χ ijk( 2) (ω1 , ω3 ,−ω 2 ),

χ ijk( 2 ) (ω1 ,−ω 2 , ω3 ),χ ijk( 2 ) (ω 2 , ω3 ,−ω1 ),χ ijk( 2 ) (ω 2 ,−ω1 , ω3 ),χ ijk( 2 ) (ω3 , ω1 , ω 2 ),

χ ijk( 2) (ω3 , ω 2 , ω1 ),(1.6.2)およびこれらの負の周波数に対応した 6 つの非線形感受率
これらの表式で添え字 i, j, k はそれぞれ独立に x, y, z の値をとるので 3x3x3=27 個の組み合わせがある。したがって 12 個の各項に 27 個の直交座標成分があるので 324 個（=12 個
ｘ27 個）のテンソル成分（複素数）が存在する。
幸い対称性に基づくいくつもの制約があるので、非線形相互作用を記述するために実際に取り扱わなければならない χ( 2)のテンソル成分は 324 個よりずっと少ない数で済む。以下で 2 次の非線形感受率テンソルの場合について対称性による特性について説明するが、同様の議論は 3 次の非線形感受率テンソル（及びより高い次数の非線形感受率テンソル）にもそのまま当てはまる。
(1)

物理量の実数性に基づく対称性

例えば SFG の場合の非線形分極は、正の和周波数ωn+ωm に対応する分極成分と負の和周波数−(ωn+ωm)に対応する分極成分の和で表されるのであった。

~
Pi (r, t ) = Pi (ω n + ω m ) exp[−i (ω n + ω m )t ] + Pi (−ω n − ω m ) exp[i (ω n + ω m )t ](1.6.3)~Pi (r, t ) は物理量なので実数でなければならない。したがって、Pi (−ω n − ω m ) = Pi (ω n + ω m ) *(1.6.4)が成り立つ(注 32)。電界成分についても同様で、その複素振幅はE j (−ω n ) = E j (ω n )*,(1.6.5)Ek (−ω m ) = Ek (ω m )*,の関係を満たす。

注 32：これは負の周波数成分は正の周波数成分の複素共役と定義したからで、定義を確認したにすぎない。

32

物理工学科博士前期課程

H21 前期光計測工学(谷）電界と 2 次の非線形分極は 2 次の非線形感受率を通じて関係付けられるので、正と負の周波数に対応する非線形感受率はしたがって次の関係を満たす(注 33a)。

χ ijk( 2) (−ω n − ω m , − ω n , − ω m ) = χ ijk( 2 ) (ω n + ω m , ω n , ω m ) *(2)(1.6.6)固有置換対称性（Intrinsic Permutation Symmetry）

先に古典的な非調和振動子の非線形感受率の表現である(1.5.47)式を固有置換対称性を用いて (1.5.49) 式に書き直した。ここでは固有置換対称性についてもう少し一般的な観点から捉えなおしてみよう。
(1.6.1)式を見ると分かるように、非線形分極 Pi (ω n + ω m ) の成分は

χ ijk( 2) (ω n + ω m , ω n , ω m ) E j (ω n ) E k (ω m ) という形をしている。この積において E j (ω n ) とEk (ωm ) の順序を変えて χ ikj( 2) (ω n + ω m , ω m , ω n ) E k (ω m ) E j (ω n ) と表しても同じ非線形分極成分を表現しているので、非線形感受率にχ ijk( 2) (ω n + ω m , ω n , ω m ) = χ ikj( 2 ) (ω n + ω m , ω m , ω n )(固有置換対称性) (注 33b)(1.6.7)という対称性の規則を導入すると、

χ ijk( 2) (ω n + ω m , ω n , ω m ) E j (ω n ) E k (ω m ) = χ ikj( 2) (ω n + ω m , ω m , ω n ) E k (ω m ) E j (ω n )(1.6.8)という関係式が得られる。これは非線形分極を表す表式において 2 番目と 3 番目の周波数に対するダミー変数を入れ替えたとき、すなわちｎと m を交換したとき、同時に 2 番目と3 番目の座標に対するダミー変数 j と k を同時に交換したものは同じものを表していることを意味している。この交換関係あるいは対称性は便宜的なもので、(1.6.8)の関係を別の形で定義することも可能である。例えば(1.6.8)式の組み合わせの一方をゼロと定義しておき、和をとるときにはもう一方の項を 2 倍するように決めておいても、結果として得られる物
理量 Pi (ω n + ω m ) の意味は変わらない。
(3)

吸収の無い媒質の場合の２つの対称性

吸収の無い媒質の場合には２種類の対称性が得られる。１つは吸収の無い媒質では非線形感受率 χ ijk (ω n + ω m , ω n , ω m ) のすべての成分は実数になるということである。すなわ( 2)ちχ ijk( 2) (ω n + ω m , ω n , ω m ) = χ ijk( 2 ) (ω n + ω m , ω n , ω m ) *(実数の条件)(1.6.9)この厳密な証明は量子力学的な χ ijk (ω n + ω m , ω n , ω m ) の表式を用いる必要があるが、古( 2)

典的な非調和振動子に対して得られた非線形感受率の表式からすぐに理解される。古典モ
注 33a：これで独立な非線形感受率テンソルの成分は半分に減る。
注 33b：これでの独立な非線形感受率テンソルの成分はさらに半分に減り、81 個になる。

33

物理工学科博士前期課程

H21 前期光計測工学(谷）デルでは 2 次および 3 次の非線形感受率は実数の係数と相互作用に含まれる周波数に対応した線形感受率の積で表された。また 1 次元の線形感受率は(1.5.39c)式で

χ (1) (ω ) =Ne 2 / m
Ne 2 / m=
ε 0 D(ω ) ε 0 (ω 0 2 − ω 2 − 2iω γ )(線形な感受率)(1.6.10)で表されていた。この式の分母の周波数の関数において、虚数部である 2iω γの項は振動子のダンピング、すなわち吸収を表している（複素誘電率または複素屈折率の虚部に相当）。
したがって、吸収のない媒質ではこの虚数項をゼロとおくことができ、線形な感受率は実数で与えられ、その積で与えられる非線形感受率も実数になる。2iωγの項が無視できるた
めには、ω02－ω2 の項が相対的に非常に大きい、すなわち入射波の周波数が電子の共鳴周波数ω0 よりも十分離れていればよい（ω0 >>ωまたはω0 <<ω）(注 34a)。これは物理的には物質の光吸収端（電子励起を伴う共鳴周波数に相当）より入射波の周波数が遠ざかるほど物質の吸収が少なくなるという事実に一致している。
吸収の無い媒質の場合の２つめの感受率の対称性は「完全置換対称性(Full permutation
symmetry)」と呼ばれている。これは「感受率の周波数パラメーターの順序は(ωp =ωn+ωmの関係は維持しつつ)対応する座標軸の添え字(i,j,k など)と同時に置換するのであれば、自由に変更しても良い（同じ値を与える）」というものである。例えばこの対称性によれば、

χ ijk( 2) (ω 3 = ω1 + ω 2 ) = χ (jki2 ) (−ω1 = ω 2 − ω 3 )(1.6.11)などとなる。この例ではω1 ，ω1，ω3 の位置が循環的にずらされ（符号も変化していることに注意）、対応する座標の添え字が周波数の順序とともに入れ替わっている。さらに実数条
件を用いるとχ ijk( 2) (ω 3 = ω1 + ω 2 ) = χ (jki2 ) (−ω1 = ω 2 − ω 3 ) = χ (jki2) (ω1 = −ω 2 + ω 3 ) *= χ (jki2) (ω1 = −ω 2 + ω 3 )(1.6.12)が得られる。同様にして

χ ijk( 2) (ω 3 = ω1 + ω 2 ) = χ kij( 2 ) (ω 2 = ω 3 − ω1 )(1.6.13)が得られる(注 34b)。
(4)非線形媒質中のエネルギー密度

一般に吸収の無い媒質中で完全置換対称性が成立することを証明するためには量子力学的な非線形感受率の表式を用いる必要があるが、電磁場の非線形媒質内でのエネルギー密度についての考察から導くことが可能である。
注 34a：これは電子的な共鳴吸収端から下の透明な周波数領域においては、(1.6.11)式が良い近似として成立することを意味している。
( 2)
注 34b：この対象性の導入で、独立な非線形感受率の成分は 1/3 に減る。2 次の非線形感受率 χ ijkの場
合独立な成分は 81 個から 27 個に減る。

34

物理工学科博士前期課程

H21 前期光計測工学(谷）線形媒質の場合、~
Ei (t ) = ∑ Ei (ω n ) exp(−ω n t )(1.6.14)nで表される入射波のエネルギー密度は~ ~~ ~U = D ⋅ E = ∑ Di Ei(1.6.15)iで表される。ここで上についているバーは時間平均をとることを意味する。電束密度ベク
トルは~~Di (t ) = ∑ ε ij E j (t ) =∑ ∑ ε ij E j (ω n ) exp(−iω n t )jj(1.6.16)nで与えられる。ここで誘電率テンソルεij を

ε ij (ω n ) = ε 0δ ij + ε 0 χ ij (1) (ω n )(1.6.17)の形で与えると、エネルギー密度は次の式で与えられる。

U = 12 ε 0 ∑ ∑ E i (ω n ) E i (ω n ) + 12ε 0 ∑ ∑ E i (ω n ) χ ij(1) E j (ω n )
*

i

*

n

ij

(1.6.18)n

この式の第 1 項は真空中での電磁場のエネルギー密度を表し、第 2 項は媒質中に分極として蓄えられているエネルギーの密度を表している（時間平均を取ったために各項の前に 1/2の係数が付くことに注意）。
非線形相互作用のある媒質中での分極によるエネルギー密度は、(1.6.18)式をより一般的にした次のような式で表される。

U = 12 ε 0 ∑ ∑ χ ij(1) (ω n ) Ei (ω n ) E j (ω n )
*

ij

n

+ ε 0 ∑ ∑ χ ijk( 2) ' (−ω n − ω m , ω m , ω n ) Ei (ω m + ω n ) E j (ω m ) E k (ω n )
*

1
3

ij k

mn

+ 14 ε 0 ∑ ∑ χ ijk( 3) ' (−ω o − ω n − ω m , ω m , ω n , ω o ) Ei (ω m + ω n + ω o ) E j (ω m ) E k (ω n ) El (ω o )
*

ij kl

mno

+"
(1.6.19)．．はとりあえず非線形感受率ではなく分極エネルギー密度 U を
ここで係数 χ ijk ' ， χ ijkl ' ，．
( 2)

( 3)

入射波の電界で展開したときの係数と考えておく。この U の表現では、電界を掛け合わせる順序は、係数 χ ijk ' ， χ ijkl ' ，．．．の座標の添え字の順序と周波数パラメーターの順序をそ
( 2)

( 3)

35

H21 前期

物理工学科博士前期課程

光計測工学(谷）れに合わせておきさえすれば、どのような順序でも構わない。したがって、完全置換対称を満たしている。
さて、この表式と非線形分極、したがって非線形感受率を関係付けるために、次の関係式を用いる。

Pi (ω n ) =∂U∂Ei* (ω n )(1.6.20)

したがって、(1.6.19)式で表される U を電界で変微分することにより、まず線形分極についてはPi (ω m ) = ε 0 ∑ χ ij(1) E j (ω m )(1)(1.6.21a)jが得られ、非線形分極に対してはPi ( 2 ) (ω m + ω n ) = ε 0 ∑jk∑ χ ' (−ω − ω , ω , ω ) E (ω ) E (ω )
( 2)
ijk

m

n

m

n

j

m

k

(1.6.21b)n( mn )Pi ( 3) (ω m + ω n + ωo ) = ε 0 ∑j kl∑ χ ' (−ω − ω − ω ,ω ,ω ,ω ) E (ω ) E (ω ) E (ω )
( 3)
ijkl

m

n

o

m

n

o

j

m

k

n

l

o

( mno )

(1.6.21c)が得られる。この表式は (1.4.8)式や(1.4.17)式で表された非線形分極の式と一致している(注36)。したがって、

χ (n ) ' が完全置換対称性を持つので χ (n ) も同様に完全置換対称性を持つと結論できる（証明終わり）。
(5)

Kleinman の対称性（Kleinman’s Symmetry）

媒質の持つ最低の共鳴周波数ω0 よりも、入射波の周波数ωi がずっと低い場合、非線形感
受率は周波数にほとんど依存しなくなる。このような低周波数領域では媒質は入射波に対して瞬間的に応答すると考えることができ、非線形分極は時間領域で~~P (t ) = χ ( 2) E 2 (t )(1.6.22)と書くことができる。ここで χ( 2)は定数である。ところで入射波の周波数ωi が共鳴周波数ω0 よりもずっと小さい場合は、吸収は無視できるので完全置換対称性が成り立ち、

χ ijk( 2) (ω 3 = ω1 + ω 2 ) = χ (jki2 ) (ω1 = −ω 2 + ω 3 ) = χ kij( 2 ) (ω 2 = ω 3 − ω1 )= χ ikj( 2) (ω 3 = ω 2 + ω1 ) = χ (jik2 ) (ω1 = ω 3 − ω 2 )( 2)= χ kji(ω 2 = −ω1 + ω 3 )

注 36： χ(n )' では最初の周波数を表す部分が χ (n ) とは逆符号のものを用いているが、これは定義の違いによるもので、 χ ijk ' (−ω m − ω n , ω m , ω n ) = χ ijk (ω m + ω n , ω m , ω n ) と定義しなおせば(1.4.8)
( 2)

( 2)

式の表式と一致する。

36

物理工学科博士前期課程

H21 前期光計測工学(谷）(Full Permutation Symmetry)
一方、低い周波数領域では χ( 2)(1.6.23)は周波数には依存しない定数とみなせるので、i, j, k の添え字を入れ替えることなく(1.6.23）式の周波数の順序を変えることができて、

χ ijk( 2) (ω 3 = ω1 + ω 2 ) = χ (jki2 ) (ω 3 = ω1 + ω 2 ) = χ kij( 2) (ω 3 = ω1 + ω 2 )= χ ikj( 2) (ω 3 = ω1 + ω 2 ) = χ (jik2 ) (ω 3 = ω1 + ω 2 )(1.6.24)( 2)= χ kji(ω 3 = ω1 + ω 2 )となる。したがって吸収がなく χ( 2)は周波数には依存しないという条件のもとでは非線形感受率テンソルの座標に関する添え字を入れ替えたものは互いに等しい ( 注 37) 。これをKleinman の対称性と呼ぶ。Kleinman の対称性が成り立つと仮定したとき、χ ijk は周波数
( 2)には依存しないので、周波数の指定を省いて非線形感受率を書いてよい。
(6)

非線形感受率の縮約表現（Contracted Notation）dil

Kleinman の対称性が成り立つときに、2 次の非線形感受率テンソルを次のような記号を用いて表す。

1d ijk = χ ijk( 2)2(1.6.25)

ここで Kleinman の対称性が成り立つためには、周波数依存性がないことが条件だったので、周波数の指定は省いた。この記号を用いると 2 次の非線形分極はPi (ω n + ω m ) = ∑ ∑ 2ε 0 d ijk E j (ω n ) Ek (ω m )jk(1.6.26)( nm )と表される。ここで dijk には周波数依存性がないとしたので 2 番目と 3 番目の添え字ｊとｋはいつでも入れ替えても構わない。そこで、dijk の表現を単純にするために次のような対応関係を定義し、添え字を書き直す。

jk : 11, 22, 33, (23, 32), (31, 13), (12, 21)
l:
1, 2,
3,
4,
5,
6

(1.6.27)

ちなみに、この添え字ｊとｋの入れ替えは、SHG の場合には Kleinman の対称性が成り立たない一般の場合でも有効である（周波数ωn=ωm であり、入れ替え可能なので）。
この対応関係を用いて非線形感受率テンソルを３ｘ６の行列として次のように書くことができる。

注 37：この対称性の導入により 2 次の非線形感受率テンソルの独立な成分は 27 個から 10 個に減る。

37

物理工学科博士前期課程

H21 前期光計測工学(谷）⎡d11 d 12 d 13 d14 d15 d 16 ⎤d il = ⎢⎢d 21 d 22 d 23 d 24 d 25 d 26 ⎥⎥⎢⎣d 31 d 32 d 33 d 34 d 35 d 36 ⎥⎦(1.6.28)

これを非線形感受率の縮約表現という（Contracted Notation）。さらに Kleinman の対象性を用いると dijk の添え字の順序は自由に入れ替えることができたので、例えば

d12 ≡ d122 = d 212 ≡ d 26

(1.6.29a)や

d14 ≡ d123 = d 213 ≡ d 25 = d 321 ≡ d 36

(1.6.29b)などである。このことから dil は Kleinman の対称性が有効であれば、実際には 10 個の独立した要素で記述できる(注 38)。したがって(1.6.28)は⎡d11
d il = ⎢⎢d16
⎢⎣d15

d12

d13

d 14

d 15

d 22

d 23

d 24

d14

d 24

d 33

d 23

d 13

d16 ⎤
d12 ⎥⎥
d14 ⎥⎦(1.6.30)と書くことができる。
SHG の非線形分極を dil の行列を用いて書くと⎡ Px (2ω ) ⎤⎡d11 d12 d13 d14 d15⎢⎥⎢⎢ Py (2ω )⎥ = 2ε 0 ⎢d 21 d 22 d 23 d 24 d 25
⎢ P (2ω ) ⎥⎢⎣d 31 d 32 d 33 d 34 d 35⎣ z⎦

⎡ E x (ω ) 2⎤⎢⎥
2⎢ E y (ω )⎥d16 ⎤ ⎢⎥
E z (ω ) 2⎥⎢⎥
d 26 ⎥⎢2 E y (ω ) E z (ω ) ⎥d 36 ⎥⎦ ⎢⎥⎢2 E x (ω ) E z (ω ) ⎥⎢2 E (ω ) E (ω )⎥y⎣ x⎦（SHG の縮約行列表現）(1.6.31)

この式は Kleinman の対称性が成り立たない場合でも有効である。SFG の非線形分極(ω3 =ω1＋ω2)は Kleinman の対称性が成り立つ場合、dil の行列を用いて表すと

注 38：dil の成分は 3x6=18 個あるが、そのうち独立なものは(1.6.30)式で示されるように以下の 10 個である。
d11, d22, d33 ,d12=d26,d13=d36,d14=d25=d36,d15=d31,d16=d21,d23=d34,d24=d32

38

物理工学科博士前期課程

H21 前期光計測工学(谷）⎡ E x (ω1 ) E x (ω 2 )⎤⎢ E (ω ) E (ω )
⎥
1
2
y
y
⎢
⎥
d16 ⎤ ⎢
⎥
(
)
(
)
E
ω
E
ω
1
2
z
z
⎥
d 26 ⎥⎥ ⎢⎢ E y (ω1 ) E z (ω 2 ) + E z (ω1 ) E y (ω 2 ) ⎥d 36 ⎥⎦ ⎢⎥⎢ E x (ω1 ) E z (ω 2 ) + E z (ω1 ) E x (ω 2 ) ⎥⎢⎣ E x (ω1 ) E y (ω 2 ) + E y (ω1 ) E x (ω 2 )⎥⎦⎡ Px (ω3 ) ⎤⎡d11 d12 d13 d14 d15⎢⎥⎢⎢ Py (ω3 )⎥ = 4ε 0 ⎢d 21 d 22 d 23 d 24 d 25
⎢ P (ω ) ⎥⎢⎣d 31 d 32 d 33 d 34 d 35
⎣ z 3 ⎦

（SFG の縮約行列表現）(1.6.32)

SFG の表現で SHG の場合と比べて係数が 2 倍になっているのは周波数の(nm)に関するサンメンションで縮退が２つあるためである。
(7)ｄの実効値：deff

入射波の方向と偏光が固定されているときに、2 次の非線形分極はスカラー量の関係として記述することが可能である。たとえは SFG の場合、

P (ω3 ) = 4ε 0 d eff E (ω1 ) E (ω 2 )(1.6.33)

また SHG の場合

P (2ω ) = 2ε 0 d eff E (ω ) 2(1.6.34)と書ける（係数が４から２になったのは入射波が１つであることに対応）。ここで deff は(1.6.26)式におけるサンメンション∑を具体的に実行することで得られる。そのために

jk

は非線形結晶の対称性と方向を知る必要がある。例えば対称性が 3m という結晶点群（3方晶系の結晶に見られる対称性）に属する結晶の場合で、２つの入射波の偏光方向が一致
xz(c-axis)φθy

図 1.6.1 結晶座標 x, y, z 軸と光の伝播方向の関係。図は正の 1 軸性結晶(正常屈折率
nx=ny=no が異常屈折率 ne より小さい)中を 光学軸(c 軸= z 軸)に対して角度θ(極角)，x 軸に対する方位角φで伝播するとして描かれている。

39

物理工学科博士前期課程

H21 前期光計測工学(谷）している場合（Type I として知られる条件）はd eff = d 31 sin θ − d 22 cos θ sin 3φ(1.6.35a)で、２つの入射波の偏光が直交している場合（（Type II として知られる条件）、はd eff = d 22 cos 2 θ cos 3φ(1.6.35b)で与えられる。ここでθは光の進行方向と結晶の z 軸(注 39)が作る角度、φは進行方向と結晶の x-z 面の作る方位角(azimuthual angle)である(図 1.6.1 参照)。
(a)
3斜晶系（Triclinic)(b)
単斜晶系（Monoclinic)(b)
斜方晶系(Orthorhombic)

c

c

c

b

b

b

a
222, mm2, mmm

a2, m, 2/m

a
1, 1
(d)
3方晶系(Trigonal)(e)
4方晶系(Tetragonal)

a

c

(f)

6方晶系 (Hexagonal)

c
a

a

a
a
3, 32, 3m, 3 , 3 m(g)
立方晶系(Cubic)aaa432, 43m , 23,m3, m3m60o

a4, 4 , 422, 4mm
4 2m , 4/m, 4/mmm

a6, 6 , 622, 6mm
6m 2 , 6/m, 6/mmm

図 1.6.2 ７つの結晶構造と点群。群論から導かれる結晶点群は全部で 32 個存在している。
例：・点群”4”に属する結晶は対称軸に対して 4 重の回転対称性を持っており，n/4 回転(n は整数)させるともとの結晶配置と一致する。
・点群“ 4 ”は 4 重の回転反射(回映)対称性を持つ。すなわち，n/4 回転させ，回転軸に垂直な面で反射させる
ともとの結晶配置と一致する。
・点群”3m”に属する結晶は対称軸に対して 3 重の回転対称性を持っていると同時に、対称軸に平行な面について
鏡面対称性(鏡映)を持っている(m は鏡面対称を表す)。
・点群”4/m” に属する結晶は対称軸に対して 4 重の回転対称性を持っていると同時に、対称軸に垂直な面について鏡面対称性(鏡映)を持っている。

注 39：光学軸、あるいは c 軸とも呼ばれる。一軸性結晶ではこの方向沿って光が進む場合は複屈折がゼ
ロになる。

40

物理工学科博士前期課程

H21 前期光計測工学(谷）(8) 結晶の対称性
結晶は原子またはイオンが規則的な配列をした 3 次元周期構造を持っている。結晶は単位格子(ユニットセル)を 3 次元的に積み重ねたものとみなすことができる。ユニットセルの３つの辺に対応した並進ベクトル a, b, c の整数倍だけ結晶をずらすと，元の結晶配列と一致する（単位格子と並進ベクトルの定義については固体物理の教科書を参照）。ユニットセ
ルを表すこれらの並進ベクトル a, b, c の長さが互いに等しいかどうか，それらの成す角度が直角かどうかに応じて，結晶構造は図 1.6.2 のように７つに分類される。また結晶は回転，鏡映，回映（回転と鏡映を組み合わせた対称操作）などの対称操作及びその組み合わせに対応した対称性を持ち，全部で 32 個の点群に分類される。結晶の対称性から結晶の線形，及び非線形感受率には制限が課せられ，ゼロでないテンソル成分の数が限られる。対称性の高い結晶（図 1.6.2 で記号の多い結晶）ほどゼロでない独立なテンソル成分の数は少なくなる。線形な感受率あるいは誘電率のテンソル（行列）成分はそれぞれの結晶構造に対応
して表 1.１のようになる。
表 1.1

結晶構造と線形な電気感受率のゼロでないテンソル成分 χ ijの関係。液体など等方性(Isotropic)の媒質の線形な電気(1)感受率は cubic な結晶と同じ対称性を示す。
結晶構造
3 斜晶系(Triclinic)

χ ij (1) または ε ij (1)⎡ xx, xy, xz ⎤⎢ yx, yy, yz ⎥⎢⎥⎢⎣ zx, zy, zz ⎥⎦⎡ xx, 0, xz ⎤⎢ 0, yy, 0⎥⎢⎥⎢⎣ zx, 0, zz ⎥⎦⎡ xx, 0, 0 ⎤⎢ 0, yy, 0⎥⎢⎥⎢⎣ 0, 0, zz ⎥⎦⎡ xx, 0, 0 ⎤⎢ 0, xx, 0⎥⎢⎥⎢⎣ 0, 0, zz ⎥⎦⎡ xx, 0, 0 ⎤⎢ 0, xx, 0⎥⎢⎥⎢⎣ 0, 0, xx ⎥⎦単斜晶系(Monoclinic)

斜方晶系(Orthorhombic)

Trigonal(3 方晶系)
Tetragonal(4 方晶系)
Hexagonal(6 方晶系)
Cubic(立方晶系)
Isotropic(等方性媒質)

41

物理工学科博士前期課程(9)

H21 前期光計測工学(谷）２次の非線形感受率と結晶の空間反転対称性

非線形感受率も線形な感受率同様，非線形媒質あるいは非線形結晶が持つ対称性によって制限される。例えば z 軸方向に４回の回転対称である場合，x 軸方向を向いた場合と，y軸方向を向いた場合とでは同じ結晶構造が見えている。したがって， χ zxx( 2)と χ zyy( 2)

は等

しくなることが分かる。非線形光学において最も重要な対称性は反転対称性(Inversion
Symmetry)である。反転対称性を持つ結晶の 2 次の非線形感受率はゼロになる。32 個の点群に属する結晶のうち 11 個の点群に属する結晶は反転対称性を持っている。したがって、
これらの結晶は 2 次の非線形感受率がゼロになることはあらかじめ分かるので、このルー
ルは非常に便利である。
反転対称性を持つ媒質において，2 次の非線形感受率はゼロになることは古典振動子モデルを用いた非線形感受率の導出ですでに明らかになっているが，以下のような考察からも定性的に理解することができる。ここでは吸収がなく，電界に対して瞬時に応答する媒質における SHG を考える。このとき SHG の非線形分極は~~~P (t ) = ε 0 χ ( 2 ) E 2 (t ) for E (t ) = E0 cos ω t(a)(1.6.36)(b)~
E (t )
Input field

Linear response

t

t

(d)
P (t )(c)~
P (t )

Nonlinear, noncentrosymmetric medium

Nonlinear, centrosymmetric mediumtt

図 1.6.3 非線形応答に関係した波形。(a)入射波の波形、(b)線形分極の波形、(c)中心対称な媒質における非線形分極の波形、(d) 中心対称な媒質における非線形分極の波形。

42

H21 前期

物理工学科博士前期課程

光計測工学(谷）~で与えられる。入射電界 E (t ) の符号を反転させると，媒質は反転対称性を持つので非線形~~

分極 P (t ) の符号も反転するはずである。したがって，反転した入射電界－ E (t ) に対して~~~
− P (t ) = ε 0 χ ( 2) [− E (t )]2 = ε 0 χ ( 2 ) E (t ) 2(1.6.37)~~が成り立つはずである。式(1.6.36)と式(1.6.37)の比較より P (t ) =− P (t ) であり，これが常に成り立つためには χ( 2)= 0 でなければならない。

非線形分極の様子はその時間依存を見たほうが直感的に分かりやすい。図 1.6.3(a)に入射~~

波 E (t ) = E exp(−iωt ) + c.c. の時間依存の波形、図 1.6.3(b)に線形分極 P (t ) = ε 0 χ

~

(1)

~
E (t ) の波形、図 1.6.3(c)に反転対称性を持つ媒質の非線形分極 P (t ) （線形成分も含む）の波形、~

図 1.6.3(d)に中心対称性を持たない媒質の非線形分極 P (t ) の波形、をそれぞれ示した。図1.6.3(c)では線形な分極成分とともに 3ωの THG 成分が重なっており、時間的な振幅の平均はゼロになっている（Optical Rectification 成分が無い）。図 1.6.3(d)では線形な分極成分
とともに 2ωの SHG 成分と DC 成分（Optical Rectification 成分）が含まれている。
(10)

その他の対称性

媒質に結晶固有の対称性があれば、さらに非線形感受率に対する制限が増えるので、結晶の持つ対称性に応じて感受率テンソルの成分の数は少なくなる。32 個の結晶の対称性に対してゼロにならないテンソル成分の一覧を表 1.2 に示す。表 1.2 の結果を dil の縮約表現を用いて視覚的に図示したものが図 1.6.4 である。
例えば 3m の結晶点群に対して、dil 行列は⎡ 0
d il = ⎢⎢− d 22⎢⎣ d 31

0

0

0

d 22

0

d 31

d 31

d 33

0

d 31 − d 22 ⎤00 ⎥⎥ =
00 ⎥⎦(1.6.38)で与えられる。ここで右側の図は縮約行列を模式的に示したもので、小さい黒丸（・）はゼロになる成分を表し、大きい黒丸（●）と白丸（○）はゼロにならない成分を表す。白丸（○）と黒丸（●）の成分は符号が互いに逆である。また同じ値の成分同士は実線で結ばれている。点線で結ばれている成分は Kleinman の対称性が有効なときのみ同じ値をとる。さらにこの例では示されていないが、Kleinman の対称性が有効なときにはゼロになるが、Kleinman の対称性が有効でないときゼロにならない成分は四角（■または□）で表される（図 1.6.4 参照）。
3 次の非線形感受率のテンソル成分についても同様に各結晶の対称性に応じて一覧を作ることができる（例えば Nonlinear Optics by Boyd の Table1.5.2 を参照）。

43

物理工学科博士前期課程

H21 前期光計測工学(谷）表 1.3 に代表的な非線形光学結晶の 2 次の非線形感受率を示した。
(11)

非線形感受率テンソル χ ijk (ω 3 , ω 2 , ω1 ) の独立な成分数について( 2)

先に 2 次の非線形感受率のテンソル成分は全部で 324 個あると述べたが、実際にはゼロでない独立な成分は非常に数が少なくなる。まず、実数条件から、独立な成分は 324 個の半分，すわわち 162 個になる。さらに固有対称性で、独立なパラメーターは 81 個に減る。
吸収のない媒質では χ ijk は実数になり、また完全置換対称性により、27 個の成分が独立成( 2)分になる。SHG 発生では縮約表現が適用でき、この場合は独立な成分は 18 個になる。さらに Kleinman の対称性（ χ ijk が周波数に依存しないこと）が成り立つ場合は、10 個の成( 2)

分のみが独立成分となる。この上に結晶固有の対称性が加わると、ゼロで無い独立な成分
はさらに少なくなる。実際、図 1.6.4 に見られるように対称性が高い結晶（対称性の記号がたくさんついた結晶群）ほど、ゼロでない独立の dil の数が少ないのが分かる。

44

物理工学科博士前期課程

表 1.2

H21 前期光計測工学(谷）結 晶 の 対 称 性に よ る テ ンソ ル 成 分 の一 覧 （ Nonlinear Optics by Boyd[1] のTable1.5.1 より転載）
結晶構造
3 斜晶
(Triclinic)
単斜晶(Monoclinic)

斜方晶(Orthorhombic)
４方晶(Tetragonal)

結晶点群

ゼロでないテンソル成分1
1
2 (y 軸についての2 重の回転対称)m (y 軸に垂直な鏡
面対称)
2/m
222
mn2
mmm
4

全ての成分が独立でゼロでない。
全てゼロ。
xyz, xzy, xxy, xyx, yxx, yyy, yzz, yzx, yxz, zyz,
zzy, zxy, zyxxxx, xyy, xzz, xzx, xxz, yyz, yzy, yxy, yyx, zxx,
zyy, zzz, zzx, zxz
全てゼロ。
xyz, xzy, yzx, yxz, zxy, zyx
xzx, xxz, yyz, yzy, zxx, zyy, zzz
全てゼロ。
xyz= － yxz, xzy= － yzx, xzx=yzy, xxz=yyz,
zxx=zyy, zzz, zxy=－zyx
xyz=yxz, xzy=yzx, xzx= － yzy, xxz= － yyz,
zxx=－zyy, zxy=zyx
xyz=－yxz, xzy=－yzx, zxy=－zyx
xzx=yzy, xxz=yyz, zxx=zyy, zzz
xyz=yxz, xzy=yzx, zxy=zyx
全てゼロ。
xyz=－xzy=yzx=－yxz=zxy=－zyx
xyz=xzy=yzx=yxz=zxy=zyx
xyz=yzx=zxy, xzy=yxz=zyx
全てゼロ。
xxx=－xyy=－yyz=－yxy, xyz=－yxz, xzy=－yzx, xzx=yzy, xxz=yyz, yyy=－yxx=－xxy=－
xyx, zxx=zyy, zzz, zxy=－zyx
xxx=－xyy=－yyx=－yxy, xyz=－yxz, xzy=－yzx, zxy=－zyx
xzx=yzy, xxz=yyz, zxx=zyy, zzz, yyy=－yxx=－xxy=－xyx
全てゼロ。
xyz= － yxz, xzy= － yzx, xzx=yxy, xxz=yyz,
zxx=zyy, zzz, zxy=－zyx
xxx= － xyy= － yxy= － yyx, yyy= － yxx= －xyx=－xxyxyz=－yxz, xzy=－yxz, zxy=－zyx
xzx=yzy, xxz=yyz, zxx=zyy, zzz
yyy=－yxx=－xxy=－xyx
全てゼロ。

4

立方晶
(Cubic)

３方晶(Trigonal)

422
4mm
42m
4/m, 4mmm
432
43m
23
m3, m3m
332

６方晶(Hexagonal)

3m(x 軸 に 垂 直 な
鏡面対称)
3 , 3m
6

6
622
6mm

6m2
6/m, 6/mmm

注 42：この表では点群の表示に Hermann-Mauguin (H-M)記号を用いている。H-M 記号では数字は回転対称性の次数を表す。H-M 記号は結晶学でよく用いられるが，分子科学では Schonflies（S）記号というのがよく用いられる。S 記号では H-M 記号の”4”は C4 などと表記される。

45

物理工学科博士前期課程

H21 前期光計測工学(谷）2 軸性結晶 (Biaxial crystal classes)
Class 1

Class 2

Class mClass 222Class mm2

1 軸性結晶(Uniaxial cryslat classes)
Class 3

Class 3m

Class 6Class 6m 246

物理工学科博士前期課程

H21 前期光計測工学(谷）Classes 6 and 4

Classes 6mm and 4mm

Classes 622 and 422

Class 4

Class 32Class 4 2m

等方結晶（Isotoropic crystal classes）
Classes 4 3m and 23

Class 432

図 1.6.4 中心対称性も持たない結晶に対する dil 行列の形。小さい黒丸（・）はゼロになる成分を表し、大きい黒丸（●）はゼロにならない成分を表す。ゼロにならない成分のうち黒丸で示された成分と反対符号をとるものは白丸（○）で表されている。また同じ値の成分同士は実線で結ばれている。点線で結ばれている成分は Kleinman の対称性が有効なときのみ同じ値をとる。さらにここでは示されていないが、Kleinman の対称性が有効なときゼロのみにならない成分は四角（■または□）で表される

47

物理工学科博士前期課程

H21 前期光計測工学(谷）表 1.3 代表的な非線形結晶の 2 次の非線形係数
結晶点群

dil (10-9 esu)

Quartz(結晶石英)

32 = D3

LiNbO3（LN）

3m = C3v

LiTaO33m = C3v

BaTiO34mm = C4v

NH4H2PO4 (ADP)

4 2 m =D2d

KH2PO4 (KDP)

42m =D2d

AgGaS2AgGaSe2
ZnGeP2LiIO342m = D2d
4 2m = D2d
4 2m = D2d
6 = C6CdSe

6mm = C6v

ZnO

6mm = C6v

CdTe
Te
ZnTe

43m = Td
32 = D3
43m = Td

ZnSe

43m = Td

β-ZnS

43m = Td

GaAs

43m = Td

GaP

43m = Td

InP
InAs
InSb
GaSb
β- BaB2O4(BBO)
CsLiB6O10 (CLBO)

43m = Td
4 3m = Td
43m = Td
43m = Td

d11 = 0.96
d14 = 0.02
d22 = 7.4
d31 = 14d33 = 98
(d22 = 4.96)(d31 = 3.06)(d33 = 46.3)
d15 = -41
d31 = -43
d33 = -16
d14= 1.2
d36 = 1.2
d14= 1.2
d36 = 1.1
(d14=139)
d36 =81(d14= 265)
d35 = -13
d36 = -10
d15 = 74
d31 = 68
d33 = 130(d15 = 5.5)(d31 = 5)(d33 = 16.7)(d14= 265)(d11 = 2201)(d36 = 255)
(d14= 226)(d36 = 75.6)(d14= 187)
(d36 = 58.7)(d14= 73)(d36 =594)(d14= 654)
(d36 = 98)((d14= 204)(d14= 398)
(d14= 997)
(d14= 1241)(d14= 1499)
d11 = 4.6
deff (Type II)= 0.95

結晶・dil =12dil(10-12m/V)(d11 = 0.402)(d14 = 0.0 08)(d22 = 3.1)
(d31 = 5.86)(d33 = 41)
d22 = 2.08
d31 = 1.28
d33 = 19.4
(d15 = -17.2)(d31 = -18)
(d33 = -6.7)(d14= 0.5)(d36 = 0.5)
(d14= 0.5)(d36 = 0.46)
d14= 56.5(d14=33.9)
d14= 111(d35 = -5.4)(d36 = -4.2)(d15 = 31)(d31 = 28)(d33 = 54)
d15 = 2.32
d31 = 2.1
d33 = 7.0
d14= 111
d11 = 922
d36 = 106.7
d14= 94.6
d36 = 31.7
d14= 78.3
d36 = 24.6
d14= 30.6
d36 = 249
d14= 274
d36 = 41.2
d14= 85.6
d14= 167d14= 418
d14= 520
d14= 628(d11 = 1.926)(deff = 0.398)

出典
[1]
[1]
[2]
[1]
[1]
[1]
[2]
[2]
[2]
[1]
[1]
[2]
[2]
[3]
[2]
[2]
[2]
[2]
[2]
[2]
[2]
[2]
[2]
[1]

χ ( 2) の値は主として SHG 発生の効率から得られたもの。括弧で括られた値は、他方の単位系から換算した値。(MKS)=0.4189x10-3 (esu).
データの出典は[1] Nonlinear Optics by Boyd, Table 1.5.3[2] CRC Handbook of Lasers with Selected Data on Optical Technology, pp.497-505.
[3] C. C. Davis, Lasers and Electro-Optis, Table 20.1.

48

H21 前期

物理工学科博士前期課程

1.7

光計測工学(谷）非線形光学応答の時間領域表示

これまでは非線形分極を２つまたは３つの周波数成分を含む入射波に対して記述した。
すなわち非線形分極を入射波の和周波数，差周波数などそれぞれの周波数に対応する周波数成分 P(ω)に分解し，周波数ωを与える入射波の周波数成分 E(ω'), E(ω'')・・・の組を考え，それらの積が P(ω)に比例するとした(注 49a)。そして成分 E(ω')の組による P(ω)への寄与の大きさを表す係数として非線形感受率χ(ω=ω'＋ω''+･･･)が定義された。これは周波数領域で非線形光学応答を記述していることになる。

~

時間的に変動する任意の入射電界 E (t ) に対する場合はどのように記述できるであろう~か？離散的な周波数成分を持つ場合と同様に， E (t ) を周波数分解し，特定の周波数成分ω
に対する非線形分極 P(ω)に寄与する電場成分の組み合わせの和（積分で表される）につい~~て同様な記述をすることは可能である。しかし入射波 E (t ) に応答する非線形分極 P (t ) を直接時間領域で記述することも可能である。周波数領域での記述と時間領域での記述は実はフーリエ変換を通じて 1 対１に対応しており，同じ非線形光学応答を表現している。時間領域での表現は，例えば入射波が短パルス（したがって多数の周波数成分を持つ）である
場合などの過渡的な非線形光学応答を直感的に記述するのに便利であり，周波数領域での記述は入射波が単色か離散的な周波数成分を持つときに便利である(注 49b)。
時間領域で非線形光学応答を記述するために，まず時間領域で線形な光学応答を考えて~みよう。任意の入射波 E (t ) に対する媒質の線形な応答，すなわち線形分極は次のように表すことができる。
∞~~
P (1) (t ) = ∫ R (1) (τ ) E (t − τ ) dτ

(1.7.1)0

ここで R (τ ) は媒質の線形な光学応答を表す関数である。すなわち時刻 t よりτ だけ前の
(1)

~

~ (1)

時刻(t−τ)における電界 E (t − τ ) が時刻 t の線形分極 P (t ) にどれだけ寄与するかを表す関~ (1)数である。時刻 t の線形分極 P (t ) は時刻 t 以前の全ての時刻における入射電界からの寄与~で決まるので，各時刻の寄与 R (τ ) E (t − τ ) の積分をとる。因果律（causality）のために(1)

時刻 t 以前（τ < 0）の寄与はないので，τ についての積分はゼロから無限大までとなってい
る。

~ (1)

~

一方，（1.7.1)式における時間領域のそれぞれの物理量（ P (t ) と E (t ) ）および係数

R (1) (τ ) はフーリエ変換を用いて，周波数領域の表現に書き換えることができる。まず，入射電界は∞ ~
E (ω ) = ∫ E (t )e iω t dt

(1.7.2a)1 ∞~E (t ) =E (ω )e −iω t dω
∫
−
∞
2π

(1.7.2b)−∞(1.7.2b)式を(1.7.1)式に代入することにより，
注 49a：例えば 2 次の非線形分極を表す(1.4.8)式

Pi (ω ) = ε 0 ∑
jk

∑χ
ω ω ω

( = ' + '')

( 2)
ijk

(ω = ω '+ω ' ' ) E j (ω ' ) E k (ω ' ' )を想起しよう。
注 49b：実用上，周波数領域あるいは周波数成分に分解したほうがほとんどの場合扱いやすい。

49

H21 前期

物理工学科博士前期課程

光計測工学(谷）∞∞ dω~P (1) (t ) = ∫ dτ ∫R (1) (τ ) E (ω )e −iω ( t −τ )−∞ 2π
0∞ dω ∞=∫dτ R (1) (τ )e iω τ E (ω )e −iω t−∞ 2π ∫0(1.7.3)∞ dω
~P (1) (t ) = ε 0 ∫χ (1) (ω ;ω ) E (ω )e −iω t−∞ 2π(1.7.4)

あるいは

ここで χ

(1)

(ω ;ω ) は線形な感受率であり，次で与えられる。
∞

ε 0 χ (1) (ω ;ω ) = ∫ dτ R (1) (τ )e iω τ

(1.7.5)0(1.7.4)式をみると線形な分極の時間変化が入射波の各周波数成分の振幅と周波数に依存した線形な感受率で与えられていることが分かる。分極の時間変化はフーリエ変換を用いて∞ dω
~
P (1) (t ) = ∫P (1) (ω )e −iω t−∞ 2π(1.7.6)で表される。(1.7.4)式の右辺と(1.7.6)式の右辺を比較すると

P (1) (ω ) = ε 0 χ (1) (ω ;ω ) E (ω )(1.7.7)であることが分かる。これは周波数領域での通常の線形な分極と入射電界の関係式である。
非線形分極も同様に時間領域で記述することができる。すなわち 2 次の非線形分極は∞∞~P ( 2) (t ) = ∫ dτ 1 ∫ dτ 2 R ( 2 ) (τ 1 , τ 2 ) E (t − τ 1 ) E (t − τ 2 )0(1.7.8)0と表すことができる。 τ 1 およびτ 2 に関する積分は因果律から正の領域についてのみとる。
つまり，時刻 t より後(未来)の E(t−τ 1)E(t−τ 2)による非線形分極への寄与はゼロ，すなわちR ( 2) (τ 1 ,τ 2 ) =0 である。ここで E(t−τ 1)と E(t−τ 2)をその周波数成分に分解，すなわちフーリ
エ変換で表すと，∞~P ( 2) (t ) = ∫−∞∞dω1 ∞ dω 2 ∞dτdτ 2 R ( 2) (τ 1 , τ 2 ) E (ω1 )e −iω1 ( t −τ1 ) E (ω 2 )e −iω2 ( t −τ 2 )
1
∫
∫
∫
−
∞
0
0
2π
2π

= ε0 ∫∞−∞dω1 ∞ dω 2 ( 2 )χ (ωσ ; ω1 , ω 2 ) E (ω1 ) E (ω 2 )e −iωσ t
∫
−
∞
2π
2π

(1.7.9)

ここで ωσ = ω1 + ω 2 とし，2 次の非線形感受率を次で与えた。
∞

∞

0

0

ε 0 χ ( 2) (ωσ ; ω1 , ω 2 ) = ∫ dτ 1 ∫ dτ 2 R ( 2 ) (τ 1 , τ 2 )e i (ω τ +ω τ )
1 1

50

2 2

(1.7.10)

H21 前期

物理工学科博士前期課程

光計測工学(谷）(1.7.9)式による時間領域の非線形分極を周波数分解すると∞P ( 2 ) (ω ) = ∫ dtP ( 2 ) (t )e +iω t =(1.7.9)式を代入
−∞

∞

∞

−∞

−∞

= ε 0 ∫ dt ∫dω1 ∞ dω2 ( 2)χ (ωσ ; ω1 , ω2 ) E (ω1 ) E (ω2 )e −i (ωσ −ω )t
∫
−
∞
2π
2π

(1.7.11)上の式で t についての積分は ωσ − ω = 0 以外ではゼロになる(注 51a)ので， ω1 と ω 2 についての 2 重積分は ω1 または ω 2 についての 1 重積分に書き換えることができる。すわなち 2 次の非線形分極の周波数領域での表式は次で与えられる。

P ( 2 ) (ω ) =∞ε0dω1 χ ( 2 ) (ω ; ω1 , ω 2 = ω − ω1 ) E (ω1 ) E (ω 2 = ω − ω1 ) (注 51b) (1.7.12)
3 / 2 ∫−∞
(2π )上の式で，条件 ω = ω1 + ω 2 のもとでの積分∑ に相当している。

∞∫ dω が 1.3 節などで示したサンメンション−∞

1( nm )

同様にして，より高次の非線形分極を時間領域の応答として記述することができる。3次の非線形分極は次のように表すことができる。
∞ dω
∞ dω
∞ dω
~
3 ( 3)1
2
P ( 3) (t ) = ε 0 ∫χ (ωσ ; ω1 , ω 2 , ω3 ) E (ω1 ) E (ω 2 ) E (ω3 )e −iωσ t∫∫−∞ 2π −∞ 2π −∞ 2π(1.7.13)

ここで ωσ = ω1 + ω 2 + ω3 とし，3 次の非線形感受率を次で与えた。
∞

∞

∞

0

0

0

ε 0 χ ( 3) (ωσ ; ω1 , ω 2 , ω3 ) = ∫ dτ 1 ∫ dτ 2 ∫ dτ 3 R ( 3) (τ 1 , τ 2 , τ 3 )e i (ω τ +ω τ +ω τ )注 51a：

∫

∞

−∞

1 1

2 2

3 3

∞

(1.7.14)f (ωσ )e −i (ωσ −ω ) t dt = f (ωσ ) ∫ e −i (ωσ −ω ) t dt = f (ωσ )δ (ωσ − ω ) 2π となる。ここで−∞f (ωσ ) は ωσ についての任意の関数である。
注 51b：係数 1/(2π)3/2 は規格化の違いにによるもので，ここではあまり気にしなくて良い。

51

物理工学科博士前期課程

H21 前期光計測工学(谷）52