室蘭工業大学応用物理コース応用光学講義ノート...2 幾何光学...

1

室蘭工業大学応用物理コース応用光学講義ノート 2013/0325版

------------------- 古典電磁気学による、光の性質と物質の光学的性質 -------------------

目次

幾何光学---2

フェルマーの原理---2

レンズ---5

波の性質---7

波の干渉---8

ホイヘンスの原理干渉および回折現象---9

反射型回折格子---13

マクスウェル（Maxwell）方程式---15

マクスウェル方程式の導出 ---16

波動方程式---22

位相速度---23

群速度---25

§§§§数学的準備・複素数表示---26

ポインティング（Poynting）ベクトル・電磁波のエネルギ

ー---28

物質中の電磁波の基本的性質

(i) 電気伝導度 0=σ の場合---35

電磁波は横波---36

光の偏光特性---37

(ii) 電気伝導度 0≠σ の場合---40

吸収係数---42

(iii) 境界面での反射率---45

物質の光学的性質の概要---49

電子振動子模型---52

ドルーデ模型---55

ドルーデ模型の光学的性質---56

プラズマ振動---58

ドルーデ模型の光学的性質 ( κ,n の周波数依存性の詳

細)---60

局所場補正：ローレンツ（Lorentz）の局所電場モデル---64

クラウジウス-モソッティ（Clausius-Mossotti）の式---70

電子振動子模型における、局所電場の効果---71

クラマース・クロニッヒの関係式---72

)(),( ωεωχ の総和則---77

結晶中での電磁波の伝播光学的異方性・複屈折---79

索引 85

注意：このファイルの著作権は、室蘭工業大学矢野隆治にあります。ファイルの内容を紙に印刷した後、鉛筆・

ボールペンなどで書き込みしてもかまいません。しかし、著者の許諾なく、ファイルの内容そのものに対して

勝手に書き換え・改変する事は、しないでください。また、間違い・ミスプリントなどが見つかった場合は、

著者（矢野隆治）まで連絡をお願いします。

2

幾何光学光の波長は、我々が目にする物体の大きさよりもはるかに小さいため、光の波としての性質を無視して、直

線的に進む光線として扱う。このような扱いをする分野を、幾何光学という。

フェルマーの原理

光線がある位置から別の位置へ伝播する時、光線はその時間が極値（極小値）をとるような道筋をたどる．

これをフェルマーの原理という。この原理から、2つの媒質の境界面での、反射、屈折の関係式を導く。

光の反射・屈折

x 軸を境界面として 0>y の屈折率 1n 、 0<y の屈折率 2n である、2 つの物質の境界面での、光の反射お

よび屈折を考える。図のように入射角度θ 、反射角度 'θ 、屈折角度φ を決めた時、それぞれの角度の間に成

り立つ関係式を求める。その際の指針となるのが、「時間を最短にするように光は伝播し、その結果として、反

射角度、屈折角度が決まる」事である。点Aから点B へ到る最短経路、および。点Aから点C へ到る最短経

路を求める。

屈折：屈折率 jn での光の速さを jj ncc /= とする。 ),0(: 1hA から点 )0,(: xD 、さらに点 ),(: 2haB − へ

進むのに掛かる時間 )(xf は、長さを光の速さで割り、

2

22

2

1

22

1

1

)()(

c

xah

c

xhxf

−++

+=

で与えられる。 )(1 xf を微分すると、

22

2222

11

1

)(

11)('

xah

xa

cxh

x

cxf

−+

−−

+=

22

1

1 )(xh

xxg

+= は、 x の増加に伴い値が増加し、

0)(,0)0( 11 >= agg である。

一方22

2

2

)()(

xah

xaxg

−+

−= はx の増加に伴い値が減少し、 0)(,0)0( 22 => agg である。

よって、 )('1 xf は、 ax << 00 となる点 0x においてのみ、 0)(' 01 =xf となる。これから、 00 xx << で

0)('1 <xf であり、 0)(' 01 =xf 、さらに axx <<0 で 0)('1 >xf なので、 )(1 xf は 0xx = で極小値を取る。

その時、

22

0

2

2

0

212

0

2

1

0

1 /

sin

)(

1,

/

sin1

ncxah

xa

cncxh

x

c

φθ=

−+

−=

+

とすれば、 0)(' 01 =xf の条件から、入射角θ と屈折角φ に関する、次の関係式を得る。

1

22101

sin

sinsinsin0)('

n

nnnxf =→=→=

φθ

φθ

θ θ’

φ

a

h1

h3

A

D x

－h2

b

C

B

0

n1

n2

3

反射：点 ),(: 3hbC とすれば、点Aから点C へ光が伝播するのにかかる時間は、次の式で与えられる。

1

22

3

1

22

1

2

)()(

c

xbh

c

xhxf

−++

+=

先ほどと同様の考察で、 )(2 xf が極小値を取る条件を求めると、

''sin)(

sin

)(0

)(

11)('

22

3

22

1

22

3

22

1

22

3122

11

2

θθθθ =→=−+

−==

+→

−+

−=

+→=

−+

−−

+=

xbh

xb

xh

x

xbh

xb

xh

x

xbh

xb

cxh

x

cxf

となる。すなわち、入射角θ と反射角 'θ は等しい。

反射、屈折角の別の導出方（幾何学的な考察）

反射

右図のように、点 ),(: 3hbC を、x 軸に対して折り返す。点C

を折り返した点を、点 ),(:' 3hbC − をする。また、光線と 0=x と

の交点を点 D とする。点 CDA →→ と進む経路と、点

'CDA →→ と進む経路は、距離が等しく、かかる時間も同じ

である。

さて、 'CDA →→ と進む経路の中で、最も距離が短い、従

って最も時間が短い経路は、3つの点 ',, CDA が直線上にある時

である。この時、入射角θ と反射角 'θ は等しくなっている。屈折

速さcの光から見ると、点Aから点 )0,(: xD まで進むには、縦に 1h 、横にx 進む必要がある。一方点D か

ら点B まで進むには、縦に 2h 、横に xa − 進む必要がある。反射の場合と同様に考えて、光線が点 BA, 間を

まっすぐに進む事を要求すればよいように見える。しかしそうすると、入射角と屈折角が同じになってしまう。

実はこの方法では、屈折率の違い（ 21 nn ≠ ）による光学的な距離（屈折率と距離を掛けた値）の違いが全

く考慮されていない。光が伝播する時間は、前ページで見たよう

に、光学的距離を真空中での光の速度cで割った値で与えられる。

屈折率の違いを考慮するため、任意の2つの点の間の距離を光学

的距離で表わした場合、点Aから点 )0,(: xD まで進むには、縦

に 11hn 、横に xn1 進むと考えると良いように思える。しかしその

場合、原点0と点D 間の水平方向の光学的な距離 D0 が、一方で

xn1 、他方で xn2 となるため、 21 nn = でなければ矛盾する事に

なる。

よって、この方法では、入射角と屈折角が一般に異なる事を説

明できない。

θ

A

C

C’

b

h3

0

－h3

D

n1

n1

θ'’

θ

A

B

a

h1

0

－h2

D

n1

n2

x

φ

4

乱反射

普通、物体の表面には。細かな凸凹があり、入射した光線は全てが同

じ方向に反射されるわけではない。それぞれの光線で入射光に対する反

射角度が決まり、反射角度がいろいろの値をとるため、反射光は、あら

ゆる方向に反射される。このような反射を乱反射という。

逆進性

光がある点P から点Q へ進む時、光は、点Q から点P へ、P からQに到達した際に通ったのと同じ光路を

逆向きに進む事も出来る。これを光の逆進性という。

見かけの深さ

水中の物体を空気中から見ると、実際の深さよりも浅く見える。これにつ

いて考察する。

右図のように、水中にある物体P から出た光線が、水と空気との境界面と

交差する点をQ とする。点Qを通る光線が我々の眼に入る。しかし、光線は

点 'P からやって来たと我々は考える。

物体P を通る境界面に対する垂線が、境界面と交差する点をO とする。物

体P から出た光線が水と空気の境界面に入射する角度をθ、空気中で屈折す

る角度を 'θ とする。空気の屈折率＝1、水の屈折率＝n とする。

図から、角度 '', θθ =∠=∠ QOPOPQ である。 '', hOPhOP == とし、OQ の長さをx とする。

'tan

tan

tan/

'tan/''tan

',tan

θθ

θθ

θθ ==→==x

x

h

h

h

x

h

x

',θθ が十分小さければ、 'tan~'sin,tan~sin θθθθ が成り立つことを利用し、

'sinsin:,1

'sin

sin~

'tan

tan'θθ

θθ

θθ

=== nnotenh

h　　（p.2の最後の式を見よ）

を得る。 3.1~n では、 3.1/~/' hnhh = となり、水中の物体は実際よりも浅く見える。

h’

h

P

Q O

P’

θ

θ’

水 n

空気 1

5

レンズ

光軸：レンズの中心部を通り、レンズの面と垂直な線を光軸という。

凸レンズ：中央部が厚く、周辺にいくほど厚さが薄くなっているレンズ。

凹レンズ：中央部が薄く、周辺にいくほど厚さが厚くなっているレンズ。

凸レンズの性質

① 凸レンズの軸に平行な光線は、レンズを通過後、軸上の一点（焦点）を通る。焦点はレンズの両側にあり、

レンズからの距離は同じである。レンズと焦点との距離を、焦点距離 f と言う。

② 焦点を通過して凸レンズに当たる光線は、凸レンズを通過後、軸と平行に進む。

③ 凸レンズの中心を通過する光線は、レンズ通過後も、その向きを変えない。

これらの性質を利用すると、凸レンズによる像の作図を行う事ができる。

レンズの公式

fba

111=+

物体がレンズから距離a の場所にあると

しよう物体がレンズから距離a の位置に

あるとして、 fa > なら実像ができ、

fa < なら、虚像ができる。

像の作図

1）実像の場合

物体のレンズからの距離をa とする。凸レンズの焦点距離は 0>f であり， 0>a である．光線 ①，③で

作図すると，図において fa > であれば倒立像が得られる。

物体に対する像の大きさ・倍率は、以下のようにして求まる。図において、三角形OAB と '' BOA は相似で

ある。よって、物体の大きさを

h 、像の大きさを 'h とすれば、

大きさの hh /' 比は、

a

b

h

h=

'

となる。これが、物体に対する

像の倍率である。

f f

①

②

③ 光軸

焦点焦点

a b

f f

a b

f f

h

h’

物体

像

B

A

A’

B’

O

6

2）虚像の場合

fa <<0 であれば、正立像が得られる。レンズを覗き込むと像が見えるので、その像を虚像という。この

時b の値は、レンズの公式

fba

111=+

より求まり、 0<b である。この場合、像（虚

像）はレンズをはさんで物体と同じ側に出来る。

図で、右の焦点の位置をF 、レンズで点Aと

同じ高さの所を点P などとする。図において三

角形 '' BFA とFPO は、相似にある。よって、

物体に対する像の大きさの比（倍率）は、

1||

1||'

>+=+

=f

b

f

fb

h

h

で与えられる。これが、虫めがね（ルーペ）で物体が拡大される（像が大きく見える）理由である。

凹レンズの性質

① 凹レンズに平行な光線が、凹レンズを通過した後の経路を逆向きに延長すると、軸上の一点（焦点）を通

る。焦点は、レンズの両側にあり、レンズからの距離は等しい。

② 凹レンズの向こう側の焦点に向かって進む光線は、凹レンズを通過後、光軸に平行に進む。

③ 凹レンズの中心を進む光線は、その向きを変えない。

像の作図

物体の位置が焦点距離よりも（a）短い時、(b)長い時の光線の結像を図に示す。像は焦点を結ばないので、

全て虚像である。(a)(b)共に像は正立であるが、像の大きさは物体よりも小さくなっている。

a |b|

f f

物体

A’

A

B B’

O

P

F

h’

h

虚像

f f

①

②

③

f f

物体虚像

(a) f f

物体虚像

(b)

7

波の性質

水や空気のように、波・振動を伝える物質を媒質という。媒質中のある場所（波源）に振動が生じると、そ

の振動に少し送れて隣の場所が振動する。これが隣から隣へと、次々と振動が伝わる。この現象を、波・波動

という。連続して振動している波のうち、波形が正弦曲線（サイン関数）で表される波を正弦波という。正弦

波が伝わる時、媒質の各点は波源の振動と同じ振動をする。

一直線上を波が伝わる時、媒質の各点は、波源と同じ振動をする。媒質の振動の周期（1 回振動する時間）

を波の周期 )s(T 、媒質が振動する速さ（1秒間に何回振動するか）を振動数または周波数 )1/s()Hz( ff = と

いう。 fT /1= の関係がある。角振動数または角周波数ω は、 fπω 2= で定義されている。また、媒質の振

動の大きさを、波の振幅 )m(A という。また、正弦波において、隣り合う山と山の頂点の間の距離 )m(λ を波

長という。なお、以下では、 f を周波数と呼ぶ事にする。ω は角周波数と呼ぶが、省略して周波数と頻繁に呼

ぶ。 f とω が記号で区別できるのでOK。

時間の経過と共に、波は移動する。波が1回振動

する間に、波は1波長だけ移動する。したがって正

弦波の進む速さ )m/s(v は、 )m/s()m/s( λfv = で

ある。

波がx 軸を正の方向に伝播しているとする。波の速さが )m/s(v なら、原点から距離 )m(x まで振動が伝わ

るのに )s(/ vx 時間が掛かるので、位置 )m(x における変位は、時刻 )s(/ vxt − での原点の変位に等しい。し

たがって原点での波を )sin(),0( tAtxy ω==

とすれば（3 角関数は、周期 π2 で元の値に戻

る。周期T の時間が経過すれば、元の値に戻る

必要あり）、点 x まで波が伝わるまでに、時間

vx / かかるので、点x での波は、原点での波の

振動に対して、時間 vx / だけ遅れる。よって、

点x での波は、 vk /ω= とおいて、

( )kxtAxv

tAv

xtAtxy −=

−=

−= ωω

ωω sinsinsin),(

となる。ここで、 λπλπω /2)/(2/ === ffvk を波数と言う。また、3角関数の中の項

kxtxv

tv

xt −−

− ωω

ωω ,,

を位置 )m(x における位相という。位相が π2 変わるごとに、振動の変位 ),( txy が元の値に戻る。

波長 λ 振幅 A

波

の

高

さ位置

波

の

高

さ位置

8

波の干渉

2 つ以上の波を重ね合わせた時に、波が強めあったり、弱めあったりする現象を干渉という。今、振幅が同

じ2つの波 )sin()(1 tAtf ω= 、 )sin()(2 φω += tAtf の波の合成（干渉）を考える。「2つの波による媒質

の変位は、それぞれの波による変位の合計に等しい」という、重ね合わせの原理が成り立つ。よって、観測点

における波の位相がそろっていれば（ 0=→+= φφωω tt ）、合成された波 )()()( 21 tftftf += は、

)sin(2)sin()sin()( tAtAtAtf ωωω =+= (1)

となり、波の振幅はさらに大きくなる。一方、位相差が180度異なれば（ πφ = ）、

0)sin()sin()( =++= πωω tAtAtf (2)

となり、下の図のように振幅はゼロとなる。

また、位相差φ が0、π 以外の場合は

( ) ( ) ( ) ( )αωφφωω ++=++= tAtAtAtf sincos22sinsin (3)

となり、振幅は位相差φ に依存して、 φcos22 +A となる。ただし、

+= −

φφ

αcos1

sintan 1 (4)

である。[(3)(4)の計算： ( )φω +tsin を展開し、 αφαφ sinsin,coscos1 BB ==+ 　とおいてみよ。]

)(tf が光の波である場合、その周波数 πω 2/ が非常に大きい（可視光で1310 Hz程度）ため、人間の目で感じる

ことが出来るのは光の強度である。光強度 I は、 )(tf の振幅の2乗に比例する｡位相がそろっている2つの波が干

渉した場合、(1)式より、光強度は24AI ∝ となり、光波が1つのときの4倍となり、位相が180度異なる波の干

渉では、光強度はゼロになる。それ以外の場合は(4)式より

( ) 2cos22 AI φ+∝ (5)

となり、2光波の位相差φ に依存して光強度が変化する。なお、光を含む電磁波の強度が振幅の2乗に比例する事

は、ポインティングベクトルの項(p.30)で説明する。

位相が180度異なる波の合成

9

ホイヘンスの原理干渉および回折現象

ホイヘンス（Huygens）の原理

波動の伝播とは、何らかの振動(水面の振動、空気密度の振動、弦の振動、電磁場の振動など)が、媒質を次々

に伝わる現象である。振動する場所を波源という。波動の伝播をホイヘンスの原理に従って説明する。下の図

において，波面 321 ,, ΣΣΣ は点線で示されている。

ある時刻に波面が波面 1Σ の所まで進んできたとすると、今度はその波面上の各点が点源となり無数の新し

い球面波（2 次球面波）が放出される。この 2 次球面波は、互いに重なった点での振幅が同じ向きの場合に、

初めて波となる。したがって 2 次波の包絡面が新しい波面，波面 2Σ となる。すぐ近くの位相の等しい波の点

を結んだ時、できた等位相面が直線になっていれば、その波を平面波という。

この説明では，点源から出る波の重ね合わせでできる波面が存在するため、元に戻る波面も生じる事になる。

しかし実際は、元に戻る波面は存在しない。（キルヒホッフ（Kirchhoff）は、後退する波動も考慮し、元に戻る

波面が存在しない事を示した。この説明は、例えば、砂川重信理論電磁気学 (紀伊国屋書店) 第 8章電磁波 §6

電磁波の回折、Huygens の原理にある。）波源から十分に遠方の波動では、波面は狭い範囲内では平面とみな

せる。したがって，十分遠方では広い範囲にわたり、近似的に平面波とみなせる。（地球の表面も丸いが、私た

ちの体の大きさからすれば、十分平面である。これと同じ。）

ヤング(Young)の干渉実験

2 つ以上の波が同時にある点で重なり、山同志あるいは谷同志の波が重なると互いに強めあい、山と谷とが

重なると弱めあう。この現象を干渉という。隣り合う山同志あるいは谷同志の間隔が波長λである。

ヤング(Young)の干渉実験では、下の図にあるように、点光源S が、光が少し透過するような小さな隙間（ス

リット） 21, SS （スリット同士の間隔a ）から等距離にある場合の、スクリーン上の点P における光の強度を

考える。

波長λ の単色の点光源S により照射される 2点 21, SS では、光

の振動は同位相であるとする。 21, SS からの球面波がスクリーン

上で重なりあう。距離 PS1 と PS2 との差が光の波長の整数倍の

とき、2つの光は強く干渉し、スクリーン上で明るい部分になる。

P 点での時間平均した光強度 )(xI は、 La << および Lx << を

満足すれば、

ホイヘンスの原理による波の伝播

(a)球面波の場合,(b)平面波の場合 (a) 球面波

2次波の合成

波面

Σ３ Σ１ Σ２

2次波の合成波面

(b) 平面波

Σ１ Σ２ Σ３

S S1

S2

L

P(x)

x θ

10

( )

+∝ xL

axI

λπ2

cos1

で与えられる。よって、スクリーン上での光強度の明暗の周期 x∆ は、次式で与えられる。

a

Lxx

L

a λ∆π∆

λπ

=→= 22

電場の可干渉性

2つの光電波が互いに干渉するか否か（可干渉であるかどうか）は、以下の例で考えよう。Youngの実験と

類似の光源などの配置を考える。2 つの点光源 21,SS から出伝播した光電場がスクリーン上で干渉するとしよ

う。2つの電場の振幅、（角）周波数、波数が同じで、それぞれ λπω /2,,0 =kE であるとし、 21,SS からスク

リーン上のある点P までの距離をそれぞれ 21, xx とする。その時、スクリーン上での光電場 )(),( 21 tEtE がそ

れぞれ次式で与えられるとする。（ ..cc は、複素共役を意味する。x を実数として、 )sin()cos()exp( xixix +=

が成り立つ。）

..)](exp[)(

..)](exp[)(

2202

1101

cckxtiEtE

cckxtiEtE

+−−−=

+−−−=

φω

φω

その時、この2つの電場の干渉によるスクリーン上の光電場の強度は、光強度が光電場の2乗で与えられる事

（p.30を参照。ここでは、とりあえず認めよう。）から、

|)])(cos(1[2

)])((exp[)])((exp[2

)])(exp[)](exp[(

)](exp[)](exp[(

))()(()(

)()()()()),((

2121

2

0

212102121

2

0

220110

220110

2

21

2

21

φφ

φφφφ

φωφω

φωφω

−+−+=

−+−+−+−−+=

−−+−−×

−−−+−−−=

+∝

∝→+=

xxkE

xxkiExxkiE

kxtiEkxtiE

kxtiEkxtiE

tEtEtI

tEtItEtEtxPE

　　　

　　　

　　　　　　

　　　

さて、電場の位相の中の定数項 21,φφ が時間的に安定なら、スクリーン上の光強度の明暗は、差 21 xx − が周

期 λπ =k/2 で明暗となる干渉縞を示す事がわかる（図は、 21 φφ = の場合を示す。）。つまり、Young の実験

において、異なる2つの光源 )(),( 21 tEtE を用いた場合でも、同じ波長（周波数）の光源、かつ =− 21 φφ 定数

であれば、スクリーン上にきれいな干渉縞（光強度の明暗）を観測する事が可能である。このように干渉縞が

きれいに見える場合、2つの光電場は、互いに可干渉である（コヒーレンスがある）という。しかし、 21,φφ が

全く相関なくすばやく時間変動する場合、干渉縞の明暗の模様はスク

リーン上で激しく移動するため、干渉縞が観測されない。このような

時、それらの光電場は、互いに干渉しない（コヒーレンスがない）と

いう。

S1

)(1 tE

)(2 tE

2x

1x

S2

P

11

単一スリットによる回折

波が障害物に当たったとき､波が障害物の裏側に回りこむ現象を回折という。右図の

ように、小さな隙間AB(スリットという)を通過した光は障害物の裏側にも伝播し、ス

リットの端Aのそばを通った光と、Bのそばを通った光が干渉することで、P点の光の

強さは、干渉の条件(位相差)により変化する。

今、下の図のように、幅a のスリットを通過して、十分遠方（ aL >> ）にあるス

クリーンでの光強度の分布を考える。スクリーン上での点P における光強度は、幅a

のスリット上の任意の点から点P に来る全ての光の干渉の結果、与えられる。下の図

で、点P 、スリットの中心、およびスクリーンの中心で作られる角度をθ とする。この時 Lx << として、スク

リーン上の点P での光強度は、以下の式で与えられる。

( )2

2sin

x

xL

a

xIλπ

∝

スクリーン上での光強度分布

スリットによる光電場の干渉計算

平面波である光電場が、スリットに垂直に入射する場合を考えよう。また、 1/,1/ <<<< LxLa が成り立つ

とする。パラメーターは、単一スリットのすぐ上の図と

右図を参照のこと。

スリット上での電場は、 )sin()( tAtE ω= の時間依存

性を持つとしよう。スリットAB 上に、スリットの端A

から距離ξ 離れた点 )(ξQ をとる。

点 )0( == ξQA からスクリーン上の点P までの距

離をR とすれば、点 )(ξQ からスクリーンまでの距離は、

θξ sin+R で与えられる。よって、スリット上の任意

の点 )(ξQ での光電場が )sin( tE ω であるなら、スクリ

ーン上の点P での光電場は、スリット上の光電場より、

電場の伝播時間 cR /)sin( θξ+ だけ遅れる。ここでc

は光電場の速度である。

以上の事柄を考慮すると、スリットを通過後、スリットに垂直な方向からθ だけ傾いた方向に進む光電場は、

A

B

P

スリットによる回折

0 2λL a

4λL a

-4λL a

-2λL a

A

B

a

L

P(x)

x θ

単スリットでの回折

A

B

a

L

θ

Q(ξ) ξsinθ

R

12

次の式で与えられる。（すぐ右上の、単一スリットの図を参照）

∫

∫∞+

∞−

+∞

∞−

−−−=

+−=

ξξρωξθω

ξξρθξ

ωθ

dc

Rt

cA

dc

RtAE

)(sin

sin

)(sin

sin)(

　　　

(1)

ここで、 )(ξρ はスリットの分布関数であり、以下のようになる。

(i)単一スリット（スリット幅a）

<−<

≤≤−=

)2/,2/(0

)2/2/(1)(

ξξξ

ξρaa

aa

　　

　 (2)

(ii)二重スリット（スリット間隔a。スリット自体の幅は十分小さいとする。）

)2/()2/()( aa ++−= ξδξδξρ )(xδ はδ関数 (3)

δ 関数の性質：∫+∞

∞−=− )()()( 00 xfdxxxxf δ

それぞれの場合について、 )(θE と光強度分布 )(xI を求めよう。

(i)単―スリットの場合

微小な幅 ξd のスリットが、 2/a−=ξ から 2/a+=ξ まで並んでいると考え、積分する。(2)式を(1)式に代

入して積分を計算する。 1/ <<Lx として、 Lx≈≈ θθsin の近似、および ωπλ /2 c= を用いると、

−−−−

−−=

−−−−=⋅⋅

−−−=−−

∫

c

Rt

c

a

c

Rt

c

aAc

c

Rt

c

Acd

c

Rt

ccAE

a

a

a

a

ωθω

ωθω

θω

ωξθω

θωξ

ωξ

θωθ

2

sincos

2

sincos

sin

sincos

sin1

sinsin)(

2/

2/

2/

2/

　　　

ここで、ここで、ここで、ここで、 bababababababa sinsin2)sinsincos(cos)sinsincos(cos)cos()cos( =−−+=+−− を用い、を用い、を用い、を用い、

)4(sin2

sinsin

sin

2)( 　　　＝

−

c

Rt

c

acAE

ωω

θωθω

θ

となる。光強度は電場の2乗の1周期分の時間平均に比例するから、 ωπ /2=T として、

( ) ( )

)5(

sin

2

sin1

2

sinsin

sin

21

2

2

22

0 0

22

2

2

−−−

=

−

=∝ ∫ ∫

　　

xL

a

xL

a

aA

dtc

Rt

Tc

acAdtE

TI

T T

λπ

λπ

ωω

θωθω

θθ

となる。 (5)式で、比例係数を省略すると(6)式が得られる。

( )2

2sin

x

xL

a

xIλπ

∝ (6)

13

(ii)二重スリットの場合

(i)と同様に(3)式を(1)式に代入して、積分を行う。

)7()]/(sin[2

sincos2

2

sincos)]/(sin[2cos)'sin(2

)2/(sin,/':)'sin()'sin(

sin)2/(sin

sin)2/(sin

)]2/()2/([sin

sin)(

−−−−⋅

=

−==

=−=−++=

−−+

+−=

++−

+−= ∫

+∞

∞−

　　　　　

　　　

　　　　

　　　

cRtc

aA

c

acRtAtA

cacRtttAtA

c

aRtA

c

aRtA

daac

RtAE

ωθω

θωωφω

θωφφωφω

θω

θω

ξξδξδθξ

ωθ

光強度分布は(7)式を2乗して時間平均をとれば得られる。 1/ <<Lx として、 Lx≈≈ θθsin の近似、および

ωπλ /2 c= を用いると、

( )

+∝

∝ xL

ax

L

aAxI

λπ

λπ 2

cos1cos22 (8)

となる。

反射型回折格子（Grating）

反射型回折格子は非常に幅の狭い多数の溝が刻まれた鏡である。光が回折格子に入射すると、回折格子に刻ま

れた溝によって、光はさまざまな方向に回折される。しかし、ある条件を満足する方向に、特に強く回折される。

ある波長λ の平面波の光が、格子間隔d の回折格子に、入射角度 inθ で入射し、出射角度 outθ で回折格子から

出てくる場合を考えよう。

光が回折格子に入射する時は、②の光線の光路が①の光線よりも ind θsin （青い線分）だけ長い距離を移動する。

一方回折格子からの出射時は、①の光線は、②の光線よりも outd θsin （赤い線分）だけ長い距離を移動する。

よって、入射・出射時の距離の差をとると、

λθθ ⋅=− md inout )sin(sin （m は整数） (1)

θin

θout

d

θin

θout

反射型回折格子の模式図

①

①

②

②

θin、θoutと dの関係

14

を満たす方向 outθ には、異なる溝からの光が干渉して強め合い、強い回折光が観測される。逆に(1)を満たさな

いような方向では干渉の結果弱め合い、特に右＝ λ)2/1( +m の時、光強度は非常に弱くなる。 0=m の時は波

長にかかわらず inout θθ = であり、これは鏡面反射に等しい。分光において特に重要になるのは 1±=m に相当

する回折光（1次の回折光とよばれる）である。いま 1+=m とすると、

λθθ =− )sin(sin inoutd (2)

であり、 inout θθ sinsin − を変化させることによって、異なる波長λ を選択することができる。

この原理を利用したのが回折格子型分光器である。下の図に分光器の構造を示す｡入射スリットを通過した光

はコリメートミラーで平行光となり、回折格子に入射する｡回折角は波長によって違うため、回折光は射出面で

は異なる位置に集光される｡射出面にスリットをおくことで一つの波長成分のみを取り出すことができる｡

入射スリット

射出スリット

回折格子

集光ミラー

コリメートミラー

回折格子型分光器の構造

15

マクスウェル（Maxwell）方程式

光は電磁波の一種である。以下で、電磁波のマクスウェル方程式の導出と説明を行う。マクスウェル方程式は、

次の(1)~(4)までの式をいう。(5)～(7)は、補助的な関係式。（ iBHED P,,,, , は3次元ベクトル）

(7)

(6)

(5)

)4(0div

)3(div

)2(rot

)1(rot

Ei

HB

ED

B

D

iD

H

BE

σµ

ε

ρ =

=

=

=

=

=∂∂

−

∂∂

−=

　　t

t

電気伝導度透磁率、誘電率、

電流密度電荷密度、

磁束密度磁場、電場、電束密度、

:::

: :

::::

σµε

ρ i

BHED

真空中での値： m/A)(T(H/m),104m),(C/V(F/m),10854.8 7

0

12

0 ⋅×=⋅×= −− πµε

PED += 0ε

P ：分極ベクトル物質中に生じた単位体積あたりの電気双極子（微小な領域での平均をとろう）

EPEP χε 0=→∝ 、 )1(0 χεε += 、χ ：電気感受率

（異なる定義もある。 χε 0 ：電気感受率 χ ：比誘電率）

IHB += 0µ

I ：磁化ベクトル物質中に生じた単位体積あたりの磁気双極子

HIHI mχµ0=→∝ 、 )1(0 mχµµ += 、 mχ ：磁気感受率、磁化率

（異なる定義。 mχµ0 ：磁気感受率、磁化率 mχ ：比磁化率）

IP, の起源：外部から物質に電場0E を加えると、それに応じる形で、原子内の原子核と電子の位置の平衡点

からのずれや、正イオンと負イオンとの位置の平衡点からのずれなどが生じる。このような正負のイオンや原

子核・電子の位置のずれを電気双極子といい、 pで表す。電荷 q− から電荷 q+ へ向かう位置ベクトルがrの

時、 rp q= で与えられる。単位体積あたりの電気双極子 pを全て集めた量（ベクトル）を、分極P で表す。

磁気双極子も同様に考える。外部から掛けられた磁場0B により誘起される単位体積あたりの磁気双極子を I

で表わす。

補足：分極のとらえ方

単位体積[ 3m ]あたりの電気双極子[ mc ⋅ ]は、その単位で見ると、23 m/cm/mc =⋅ となる。つまり、分極

は、単位体積あたりの電気双極子ベクトルの和として定義もできるし、外部電場によって電気双極子が出来た

時の単位面積当たりの電荷の量としても定義もできる。

),(),()],(rot[

),(),()],(rot[),,(),()],(rot[

),(),(),(),(div

rotdiv,)),(),,(),,((),(3

tAy

tAx

t

tAx

tAz

ttAz

tAy

t

z

tA

y

tA

x

tAt

tAtAtAt

xyz

zxyyzx

zyx

zyx

xxxA

xxxAxxxA

xxxxA

xxxxA

∂∂

−∂∂

=

∂∂

−∂∂

=∂∂

−∂∂

=

∂∂

+∂

∂+

∂∂

=

= は以下の演算に対して、次元ベクトル

16

マクスウェル方程式の導出

以下で、(1)~(4)のマクスウェル方程式の導出を行う。

ガウス（Gauss）の法則

真空中で、原点Oに電荷 )0(>q の点電荷がある。それから距離r 離れた位置での電場の大きさE は

2

04

1

r

qE

πε=

であり、電場の向きは、原点から外向きである。この電場を、半径r の球で囲む。球の面積は24 rπ なので、

qrEq

rr

qrE =⋅→=⋅=⋅ 2

0

0

2

2

0

2 444

14 πε

επ

πεπ

電場Eの向きは常に外向きなので、半径r の球の表面での、面に垂直な外向きの単位ベクトルをnとすれば、

その点の球の微小面積をdS として、

qdSqdSr

qdS

r

qdS

SSSS

=⋅→==⋅=⋅ ∫∫∫∫ nEnrnE 02

0

03

0

004

1

4

1ε

πεε

πεεε

とかける。この結果を一般化すると、電荷分布 )(xρ で与えられる閉曲面S

で囲まれた領域V において、

∫∫ =⋅VS

dVdS0

)()(

ερ x

nxE

を得る。さて、閉局面S で囲まれた領域V でガウスの定理： dVdS VS ∫∫ =⋅ )(div)()( xAxnxA が成り立つ。

ここで、 )(xA は任意のベクトルである。関係式 )()( 0 xExD ε= を用いて電束密度 )(xD の式に書き換え、

さらにV は任意に取れる事から、次のガウスの法則の微分形を得る。

)()(div)()(div)](div[)( 00 xxDxxDxEnxE ρρεε =→===⋅ ∫∫∫∫VVVS

dVdVdVdS

真空中で、原点Oに磁荷 )0(>mq の点磁荷がある。それから距離r 離れた位置での磁場の大きさH は

2

04

1

r

qH m

πµ=

であり、磁場H の向きは、原点から外向きである。この磁場を、半径r の球で囲む。球の面積は24 rπ なので、

mmm qrH

qr

r

qrH =⋅→=⋅=⋅ 2

0

0

2

2

0

2 444

14 πµ

µπ

πµπ

磁場の向きは常に外向きなので、半径r の球の表面での、面に垂直な外向きの単位ベクトルをnとすれば、そ

の点の球の小面積をdS 、磁束密度 )()( 0 xHxB µ= の関係式を用いると、電荷の場合と同様の計算で、

m

SS

m

S

m

S

m

S

qdSdSqdSr

qdS

r

qdS =⋅=⋅→==⋅=⋅ ∫∫∫∫∫ nxBnxHnrnxH )()(

4

1

4

1)( 023

0

00 µππµ

µµ

を得る。さて、単磁極のみ（N極または、S極のみ）は見つかっていない。よって 0=mq であり、

0)()(0 =⋅=⋅ ∫∫SS

dSdS nxBnxHµ

E

n

x

O

17

である。これは非常に単純な場合であるが、これがガウスの法則の積分形である。微分形は、次のとおり。

0)(div =xB

まとめると、

0)(div),()(div == xBxxD ρ

以下では、電磁場に時間依存性を認め、 ),( txE 、 ),( txB 等と書いて議論する。

ファラデー（Faraday）の電磁誘導の法則

閉じた導線回路の近くにある磁石が動くと、回路内に電流が流れる。回路の抵抗R 、電流の大きさ I 、起電

力をφ とし、閉回路で囲まれる任意の曲面S を貫く磁束をN とした時、

dt

dNRI −==φ （1）

で与えられる。この式は、実験より得られた式である。磁束の変化のあ

る空間に導線回路があれば、磁束の変化により誘起された電場により、

回路内に電流が流れる。

真空中の閉曲線 C によって囲まれる磁束の時間的変化により電場

),( trE が発生する。この電場 ),( trE により、点電荷eには、 EF e= の力が作用する。rは、閉曲線C上の

点。この力の作用に抗って、点電荷eが閉曲線Cを1周すれば、外部にする仕事の量は

∫ ⋅C

dte rrE ),(

である。ここで積分は閉回路を反時計回りに1回りする積分。起電力φ に電荷eを掛けた値が、電場が電荷に

する仕事量なので、

∫ ⋅=C

dt rrE ),(φ （２）

が成り立つ。

磁束N は、磁束密度 ),( txB の、閉曲線 Cで囲まれる任意の曲面 S上の面積分で表される。 xは、曲面 S

上の点。

∫ ⋅=S

dStN nxB ),( （３）

（1）が成り立つので、

∫∫ ⋅∂∂

−=⋅SC

dStt

dt nxBrrE ),(),( （４）

となる。ストークス（Stokes）の定理： ∫∫ ⋅=⋅SC

dSd nxArrA )(rot)( を )(xE に対して用いると、

∫∫∫ ⋅−=⋅=⋅SSC

dStdt

ddStdt nxBnxErrE ),(),(rot),( （５）

N

B n

dr

C

18

が成り立つ。曲面S を固定すると、積分領域は時間と共に変わる事はないので、右辺の時間変化は ),( txB の

中のtで与えられる。（注意： xは単なる積分変数。時間に依存する変数ではない。）よって、以下のように時

間微分を積分の中に入れる事ができる。

∫∫ ⋅∂

∂−=⋅−

SS

dSt

tdSt

dt

dn

xBnxB

),(),( （6）

（5）～（6）から、

0),(

),(rot),(

),(rot =⋅

∂

∂+→⋅

∂∂

−=⋅ ∫∫∫SSS

dSt

ttdS

t

tdSt n

xBxEn

xBnxE （7）

積分領域S は任意なので、この式が常に成り立つためには、被積分関数がゼロである事すなわち、

0),(

),(rot =∂

∂+

t

tt

xBxE （8）

を得る。

アンペール（Ampere）の法則

定常電流のまわりの真空中に、静磁場が発生する。実験によると、電流 I とその

周りに出来る磁束密度Bとの間には、

IdC

0)( µ=⋅∫ rrB （１）

の関係式（定常電流なので、 )(xB に時間依存性は無い。）が成り立つ。ここで線積分は、電流の周りを一回り

する任意の閉曲線C 上の線積分であり、rは径路C 上の点である。

閉曲線Cで囲まれる任意の曲面S での電流密度を )(xi とすると、 xを曲面S 上の点として、

∫ ⋅=S

dSI nxi )( （２）

（1）（2）より、

∫∫ ⋅=⋅SC

dSd nxirrB )()( 0µ （３）

この関係式にストークスの定理 ∫∫ ⋅=⋅SC

dSd nxBrrB )(rot)( を適用すると、

( ) 0)()(rot)()(rot)( 00 =⋅−→⋅=⋅=⋅ ∫∫∫∫SSSC

dSdSdSd nxixBnxinxBrrB µµ （４）

をえる。曲面Sの取り方は任意なので、この式が常に成り立つためには、被積分関数がゼロである事すなわち、

)()(rot 0 xixB µ= （５）

をえる。これが時間変化する電流 ),( txi に対しても成立する事を要請すると、以下の式を得る。

),(),(rot 0 tt xixB µ= （６）

ｄｒｒｒｒ

n(x)

i(x)

dr

S

19

電荷の保存則

電荷の総量はいかなる物理的過程においても一定不変である。この事を式で表そう。空間内に閉曲面Sで囲

まれた領域Vを考える。閉曲面Sの外向きの法線方向の単位ベクトルを )(xn とする。Sを通って、外からこ

の領域に流れ込む電荷の量は、電流密度を ),( txi とすると、単位時間当たり、

dStS

∫ ⋅− )(),( xnxi （１）

である。 )(xn の方向を外向きにとったから、この領域に流れ

込む量には（マイナス）の符号をつける。

電荷の総量が不変であるという事は、V内に流れ込んだ電気

がどこにも失われずに、そのままV内における電荷の量の単位時間当たりの増加量（あるいは、増加の割合）

に等しくなる事を意味する。この関係は、次式のとおりである。

dStdVtdt

d

SV

∫∫ ⋅−= )(),(),( xnxixρ （２）

閉曲面を固定（時間的に変化しない閉曲面を仮定）すると、右辺の時間変化は ),( txρ の中のtで与えられる。

よって、時間微分を積分の中に入れる事ができるので、

dVt

tdVt

dt

d

VV

∫∫ ∂∂

=),(

),(x

xρ

ρ （３）

となる。Gaussの定理 dVdS V∫∫ =⋅ )(div)()( xAxnxAS を ),( txi に対して用いると、

dVtdStVS

∫∫ =⋅ ),(div)(),( xixnxi （４）

となるので、

0),(div),(

=

+∂

∂∫ dVt

t

t

V

xixρ

（５）

をえる。領域Vは任意にとれるので、この式が常に成り立つためには、被積分関数がゼロである事すなわち、

0),(div),(

=+∂

∂t

t

txi

xρ （６）

これが、電荷の保存則である。

アンペール-マクスウェル（Ampere- Maxwell）の法則

ここで、アンペールの法則と電荷保存の法則をひとつにまとめる。定常電流により静磁場が出来る時のアン

ペールの法則を、電流が時間依存性を持っている場合にも拡張し、我々は

),(),(rot1

0

tt xixB =µ

（７）

を要請した。両辺の発散（div）をとると、

n

i x S

20

),(div)],(div[rot1

0

tt xixB =µ

となる。ベクトル解析の公式により、左辺は恒等的にゼロ（任意のベクトル )(xA について、 0)(divrot =xA

が成り立つ。次ページの補足参照）である。したがって、

0),(div =txi

となるべきである。しかしこの結果は、先ほど導いた電荷の保存則と矛盾する。この困難は、マクスウェルに

よりアンペールの法則を次のように拡張することで、解決された。

t

ttt

∂∂

+=),(

),(),(rot1

0

0

xExixB ε

µ、（８）

t

t

∂∂ ),(

0

xEε ：変位電流

（８）式の発散（div）をとると、

t

tt

t

tt

t

ttt

∂∂

+=

∂∂

+=∂

∂+=

),(),(div

)),(div(),(div

)),(div(),(div),(divrot

1 00

0

xxi

xExi

xExixB

ρ

εε

µ

　　　　　　　

となるべきである。左辺は、任意のベクトル )(xA で 0)(divrot =xA が成り立つ事から、ゼロになることが

分かる。右辺は、電荷の保存則そのものであり、ゼロである。このようにして、アンペールの法則の矛盾が解

消された。

さて、もし 0),( =txi なら、

t

tt

t

ttt

∂∂

+=→∂

∂+=

),(0),(rot

1),(),(),(rot

10

0

0

0

xExB

xExixB ε

µε

µ

となる。この式は、 0),( =txi の時、「電場の時間変化に伴い、そのまわりの真空中に磁場を生じる」事を示

す。すなわち、ファラデーの電磁誘導の法則と逆の現象の存在を示している。変位電流 tt ∂∂ /),(0 xEε は磁場

を生じさせる一種の電流である。それで、“変位電流”の名前が付いた。なお、マクスウェルの方程式の第2式（2）（p.15）は、

),(),(),,(),( 00 tttt xExDxHxB εµ ==

の関係式を用いると、次式に示すようにして得られる。

t

ttt

t

ttt

∂∂

+=→∂

∂+=

),(),(),(rot

),(),(),(rot

10

0

xDxixH

xExixB ε

µ

21

補足： ),,( zyx=x とし、3次元ベクトル ( ))(),(),()( xxxxA zyx AAA= の時、 0))(div(rot =xA の説明

0)()()(

)]()([)]()([)]()([)](rot[div

)()()()()](rot[

)()()](rot[),()()](rot[),()()](rot[

)()()()(div

222222

222222

22

=

∂∂∂

−∂∂

∂+

∂∂∂

−∂∂

∂+

∂∂∂

−∂∂

∂=

∂∂∂

−∂∂

∂+

∂∂∂

−∂∂

∂+

∂∂∂

−∂∂

∂=

∂∂∂

−∂∂

∂=

∂∂

−∂∂

∂∂

=∂∂

∂∂

−∂∂

=∂∂

−∂∂

=∂∂

−∂∂

=

∂∂

+∂∂

+∂∂

=

xxx

xxxxxxxA

xxxxxA

xxxAxxxAxxxA

xxxxA

zyx

xyzxyz

yzyzx

xyzzxyyzx

zyx

Axyyx

Azxxz

Ayzzy

Ayz

Axz

Axy

Azy

Azx

Ayx

Azx

Ayx

Az

Ayxx

Ay

Ax

Ax

Az

Az

Ay

Az

Ay

Ax

　　　　　

よって

Maxwell方程式のまとめ電荷がある場合の電磁波の方程式

ガウスの法則

マクスウェルの法則アンペール

ファラデーの法則

　　　　　

　

　　　

=

=

=∂

∂−

=∂

∂+

),(),(div

0),(div

),(),(

),(rot

0),(

),(rot

-

tt

t

tt

tt

t

tt

xxD

xB

xixD

xH

xBxE

ρ

真空中および空気中では、次の関係式を要請する。

EDHB 00 , εµ == 。

空間的に一様な媒質中では、次の関係が成り立つとする。

EDHB εµ == , 。

さらにガラスなどの非磁性材料では、 0µµ = である。

22

波動方程式

1次元で、 z 方向に伝播する波を表わす式（波動方程式）は、

2

2

22

2 ),(1),(

t

tzU

uz

tzU

∂∂

=∂

∂

で与えられる。この解は、（x を実数、 1−=i として、 xixix sincos)exp( += の関係がある。）

)](exp[),cos( 00 kztiUkztU −−− ωω 　

などで与えられる。ただし、 uk=ω の関係がある。（この関係は、波動方程式に解を代入して確認せよ。）

この波が伝播する速さを求めるため、指数の中の因子（位相）が時間的に一定と仮定して時間微分をする。

ukdt

dz

dt

dzk

t

kzt==→=→=

∂−∂ ω

ωω

0)(

この速さu を位相速度という。 ω,k を、それぞれ波数、（角）周波数という。ω は、単位時間当たりに波が振

動する回数、すなわち周波数 f と fπω 2= の関係にある。波の波長をλ とすれば、波が 1回振動すると 1波

長λ だけ波は進む。波が1秒で f 回振動すると、 λf だけ進む。よって波の速さu は、 λfu = で与えられる。

また、波数k と波長λ には、 kfukfk /2/2/ πλλπω =→=== の関係がある。

1次元の波は、

)](exp[),cos( 00 kztiUkztU −−− ωω 　

などで与えられた。3次元では

)](exp[0 xk ⋅−− tiU ω

と書く事が出来る。kは波数ベクトル、ω は電磁波の角周波数である。

ここで、 =⋅− xktω 定数の状況を考える。簡単のため、 0=⋅− xktω とする。

tt ωω =⋅→=⋅− xkxk 0

なので、xベクトルとkベクトルの内積が一定の値 tω をとる。これは、xベクトルのkベクトルと平行な方

向の成分が、一定値 tω を取る事を意味する。そのような xベクトルの終点の集合は、kベクトルと垂直な平

面である。

補足： xベクトルを、波数ベクトル //kkek = と平行および垂直な

成分に分ける。 ⊥kk ee ,// は、大きさ1の単位ベクトルである。

1||||,0 ////

////

===⋅

+=

⊥⊥

⊥⊥

kkkk

kkkk xx

eeee

eex

　

すると、

k

tx

tkxtxxt

k

kkkkk

ω

ωωω

=→

=→=+⋅→=⋅ ⊥⊥

//

////// )( eekxk

となり、 xベクトルのk ベクトルと平行な成分は一定の値をとる。

k x

k

Y

Z

23

位相速度

一様な空間（電磁波の伝播方向などによらず、ε 誘電率、µ 透磁率が定数である空間）を伝播する電磁波は、

次の波動方程式を満足する。

2

2

2

2

2

2

2

2

0,zyxt ∂

∂+

∂∂

+∂∂

≡∆=

∂∂

−∆ ラプラシアン Eεµ

§ 波動方程式の導出：まず、　　0rot =∂

∂+

t

HE µ でrot をとり、 0rot

EH =

∂∂

−t

ε を代入する。

0rotrotrotrotrotrotrotrotrot02

2

=∂∂

+=∂∂

∂∂

+=∂∂

+=

∂

∂+=

ttttt

EE

EEHE

HE µεεµµµ

恒等式： AAA ∆−= graddivrotrot 、および　0div =E を用いると、波動方程式を得る。

00graddiv2

2

2

2

2

2

EE

EE

EE =

∂∂

−∆→=∂∂

+∆−=∂∂

+∆−ttt

εµµεµε

§ AAA ∆−= graddivrotrot の確認

),()],(grad[div

),(),(div

),(),(),(),(

),(),(),(

),(),(),(

),(),(),(),()],([rotrot

),(),()],(rot[

),,(),()],(rot[),,(),()],(rot[

)(,

)(,

)()(grad,

),(),(),(),(div

2

2

2

2

2

2

2

2

2

2

2

2

2

2

tAt

tAtx

tAzyz

tAz

tAy

tAxx

tAzy

tAz

tAyx

tAzy

tAxz

tAxy

tAx

tAzz

tAy

tAxy

t

tAy

tAx

t

tAx

tAz

ttAz

tAy

t

z

f

y

f

x

ff

z

tA

y

tA

x

tAt

xx

x

xzyx

xzy

xzy

zxxyx

xyz

zxyyzx

zyx

xxA

xxA

xxxx

xxx

xxx

xxxxxA

xxxA

xxxAxxxA

xxxx

xxxxA

∆−=

∆−∂∂

=

∂∂

+∂∂

+∂∂

−∂∂

+∂∂

+∂∂

∂∂

=

∂∂

+∂∂

−∂∂

+∂∂

∂∂

=

∂∂

+∂∂

−∂∂

∂∂

+∂∂

∂∂

=

∂∂

−∂∂

∂∂

−∂∂

−∂∂

∂∂

=

∂∂

−∂∂

=

∂∂

−∂∂

=∂∂

−∂∂

=

∂∂

∂∂

∂∂

=∂

∂+

∂

∂+

∂∂

=

　　　　　　　

　　　　　　

　　　　　　

　　　　　　

　　　　　　

　　

他の成分も同様にして求める事ができ、 AAA ∆−= divgradrotrot が成り立つ事を示す事ができる。

そこで、一様な空間（ µε , が一定の値である）で z+ 方向に伝播する単一周波数ω の電磁波を考える。電磁

波は横波であるため、電場はx 成分 ),( tzEx を持つとする（ zy, 成分＝ゼロ。）。（一様な空間を伝播する電磁

波が横波である事は、p.35で示す。電磁場の振幅ベクトルは互いに直交するため、磁場は y 成分のみを持つ。）

)](exp[),( 0 kztiEtzE xx −−= ω

もし観測者が常に電場のある一定値のところを観測しながら移動すれば、観測者は、一定=− tkz ω の条件を

満たすように移動しなければならない。一定=− kztω を時間微分すると、観測者は z+ の向きに

kdt

dzc

dt

kztd ωω==→=

−0

)(

24

の速さcで移動する。 k/ω をこの波の位相速度と呼ぶ。（この時点では、k の値は確定していない。）また、

波長 k/2πλ = である。

波動方程式に電場 ),( tzEx の式を代入し、 εµω ,,,k に関する次の式を得る。

00)](exp[)()(0 22

0

22

2

2

=−→=−−−−→=

∂∂

−∆ εµωωωεµεµ kkztiEiikt

x E

よって、k として、次の値を得る。（ガラスなどの非磁性材料を想定しているので、透磁率 0µµ = である。）

n

cc

kck 0

0

00

0000

0

/

11, =======

εεεε

µεµεεε

εµω

εµω 　

ここで、 0c は真空中（ 00 , µµεε == ）での位相速度 )m/s(10998.2/1 8

000 ×== µεc である。0/εε=n は

物質の屈折率である。媒質中での位相速度cはnを用いて ncc /0= となり、真空中での位相速度の n/1 になる。

一方、 z− 方向へ伝播する電場は

)](exp[),( 0 kztiEtzE xx +−= ω

で与えられる。一般には、波は z± 方向へ伝播する波の重ね合わせで書ける。

)](exp[)](exp[),( 00 kztiEkztiEtzE xxx +−+−−= −+ ωω

磁場の場合も、電場の場合と同様の計算で、一様な空間（ µε , が定数）での波動方程式を得る。

0rot E

H =∂∂

−t

ε でrot をとり、　　0rot =∂

∂+

t

HE µ を代入する。

00divgrad

0rotrot)(rotrotrotrotrotrotrot0

2

2

2

2

2

2

HH

HH

HH

HHEH

EH

=

∂∂

−∆→=∂

∂+∆−→

=∂

∂+=

∂∂

−∂∂

−=∂∂

−=

∂∂

−=

tt

ttttt

εµµε

µεµεεε

　　

一様な空間で z± 方向に伝播する電場の振幅がx 成分を持てば、磁場の振幅は、 y 成分のみゼロでない。（p.36

を見よ。）そこで電場と同じ周波数ω を持つ磁場の y 成分を

)](exp[),( 0 kztiHtzH yy −−= ω

とおく。 ),( tzH y を波動方程式に代入する事で、磁場も電場と同様、位相速度 εµω /1/ == kc で空間を伝播

する事が分かる。

次に、電場と磁場の振幅 00 , yx HE の間の関係式を求める。マクスウェル方程式の 1つである次式において、

x 成分に注目すると（ εµω /1/ == kc ）、

00000//0rot xxyxyxy EEk

HEiikHtEzHt µε

εεωεωεε ==→−=−→∂∂+=∂−∂→=

∂∂

− E

H

となる。電場と磁場の振幅 00 , yx HE の間には、物質中のインピーダンス εµ /=Z とおけば、

Z

EEEH x

xxy0

000/

1===

εµµεε

の関係がある。真空中では、 )(377/ 000 Ω=== εµZZ となる。 0Z を、自由空間（真空）のインピーダ

ンスという。

25

群速度

周波数および波数が近い2つの電場（平面波）の重ね合わせを考える。図では、赤色（破線）と青色（実線）

の2つの波が該当する。

( ) ( )])()[(exp])()[(exp 00 xkkExkkEE ⋅−−−−+⋅+−+−= ∆ω∆ω∆ω∆ω titi

これは、

( ) ( ) ( )( ) ( )][expcos2

]][exp][[exp][exp

0

0

xkxkE

xkxkxkEE

⋅−−⋅−=

⋅−++⋅−−⋅−−=

tit

tititi

ω∆ω∆

∆ω∆∆ω∆ω

　

となる。2つの波の重ね合わせ（緑色の波）は、包絡線（緑の曲線を囲む 2つの黒色破線）を持つ波の伝播で

ある。なお黒色の破線は関数 ( )xk ⋅∆−∆± tωcos を表わす）。この包

絡線（黒色の破線）の移動する速度は、位相速度 k/ω では進まない。

簡単のため、 k∆ が z 成分のみを持つとすれば、cosの中の因子が時

間的に変化しない条件から、2つの波の重ね合わせは、

ku

dt

dzzkt

dt

dg ∆

∆==→=⋅∆−∆

ωω 0)(

で進む。 gu を群速度という。

電磁波の場合、真空での位相速度 0c を用いて、 kn

c0=ω から、

−=−=+=

==−

dk

dn

n

k

n

c

dk

dn

n

kc

n

c

dk

dnkc

n

ck

n

c

dk

d

dk

dug 10

2

00

1

000ω

をえる。たいていの光学的物質（光が伝播するような物質）では、屈折率は、周波数の増加（波数の増加）に

伴い大きくなる。すなわち、 0/ >dkdn 。よって群速度 gu は位相速度 ncc /0= よりも遅い。光信号は、波

の塊（光パルスの組合せ）で送らなくてはならないので、信号は、位相速度よりも速く送ることはできない。

真空中を伝播する電磁波では、 0/ =dkdn である。また、空気中を伝播する可視光の電磁波では、空気中に

可視光の波長を吸収する物質がほとんど存在しないため、 0/ ≈dkdn と見なせる。この時、電磁波の群速度と

位相速度は等しくなる。さらに、最初にパルス波形を保っていた電磁波の周波数ごとの位相関係が伝播中も保

たれるため、パルス波形の形状は物質中を伝播しても、壊れる事無くそのままである。可視領域における真空

中や空気中での屈折率 1=n である事から、光パルスの群速度は、 (m/s)10998.2 8

0 ×=c で与えられる。

一方 0/ ≠dkdn である物質では、電磁波の位相速度が周波数ごとに異なる。これが原因で、電磁波のパルス

波形を保つために必要な周波数ごとの位相関係が崩れるため、パルス幅が拡がる。よって、最初に十分短いパ

ルス電磁波を用いても、 0/ ≠dkdn の物質中では、伝播に伴い、パルス幅が拡がる。

26

群速度に関する、論理的な考察

z方向に進む、中心周波数 0ω の平面波の波束を考える。波の周波数スペクトルは、 δωωωδωω +<<− 00

に分布し、次のように書き表せるとする。

∫∫+

−

+

−

−−=

−−=

δωω

δωω

δωω

δωω

ωωωωωωω

0

0

0

0

)exp)()(exp)(),(00

zc

ntiAdz

c

ntiAdtzu

ここで ω/0kcn = は、波束が伝播している媒質の屈折率で、 λπ /2=k は波数である。もし、 )(ωn のω 依

存性が小さいなら、次のような展開が可能である。

ωωω )()( nf = と置けば（テーラー展開の1次まで）、

　　　　　　　　

　　00

)(]))(([

)()()(

)(~)(],)([)(

00000ωωωω

ωωω

ωωωωωω

ωω

ωωωωωω

ω==

−+=−+∂∂

=d

ndn

d

dfffn

d

df

これより波束は、

( )

[ ]

∫

∫

∫

∫

∫∫

+

− =

+

− =

+

− =

+

− =

+

−

+

−

−−−=

−−=

−−−

−−=

−+−−

+−=

−++−−−

+−=

−−=

δωω

δωω ωω

δωω

δωω ωω

δωω

δωω ωω

δωω

δωω ωω

δωω

δωω

δωω

δωω

ωωω

ωωωω

ωω

ωωω

ωωω

ωωω

ωωω

ωωωω

ωωωωω

ωωω

ωωωωωωωωω

ωωωω

ωωωω

0

0 0

0

00

0

00

0

00

0

0

0

0

])([)(exp)(

)(exp

])([)(exp

)(exp)(

)(])([1

)(exp)(

exp)(

)(])([1)(

exp)(exp)(~

expexp)()exp)(),(

0

00

0

000

0

0

0

00

0

0

0

0

000

0

00

0000

00

d

nd

c

ztiAdu

zc

ntiu

d

nd

c

ztiz

c

ntiAd

zd

nd

citiz

c

nitiAd

zd

nd

cc

nititiAd

zc

nitiAdz

c

ntiAdtzu

　　　　

　　　

　　　

　　　

　　　

　　　　

　　　　

0u は、この波の包絡線を表す。この波が伝播する速さ gu は、 0u の指数関数の中の位相項＝一定とし、位相項

の時間微分を行う事で、得られる。

)(

0

0

00

0

0

0

00

)(

])([

])([110

])([

ω

ωω

ωω

ωωωω

ω

ωω

ω

ωωωω

ωωω

ωω

k

g

g

dk

d

c

nd

du

d

nd

cu

dt

dz

dt

dz

d

nd

cd

nd

c

zt

t

=

=→

=≡→⋅=→=

−

∂∂

=

=

==

こうして、いわゆる群速度 gu が求まる。

注）0

0

c

nkk

n

c ωω =→=

27

§§§§数学的準備：複素関数表示

電磁気学や光学では、正弦関数的に時間変化する電場や磁場を、 ..cc をその複素共役として、

)](exp[),(

)](exp[),(

0

0

xkHxH

xkExE

⋅−−=

⋅−−=

tit

tit

ω

ω または

..)](exp[),(

..)](exp[),(

0

0

cctit

cctit

+⋅−−=

+⋅−−=

xkHxH

xkExE

ω

ω

のように、指数関数を用いて書く事がある。その理由は、この記法だと、微積分の計算が簡単になるからである。

さて、この講義ノートでは、

)](exp[),()],(exp[),( 00 xkHxHxkExE ⋅−−=⋅−−= tittit ωω

と書き、その実数部分を考える、という立場で、さまざまな物理量の計算や議論を行っていく。このような取り

扱いをする際の注意点や、計算が便利になる結果について、以下で述べる。

例として次の関数を考える。 0, 00 >BA として、

)cos()(),cos()( 00 ba tBtbtAta φωφω +=+= （１）

複素振幅 )exp()exp( 00 ba iBBiAA φφ −=−= 、と定義すると、

)]exp(Re[)()],exp(Re[)( tiBtbtiAta ωω −=−= （２）

と書き直せる。Reは、実数部分を示す。（１）の代わりに

)exp()(),exp()( tiBtbtiAta ωω −=−=

と書く場合もあるが、この実数部分を考えよ、という事を意味するものと考える。この講義ノートでも、この方

式で物理量を記述する。する。たいていの計算では、この置き換え（読み替え）は問題ない。例外は正弦関数の

積または冪が関与する場合である。

まず例外を見る。 )(),( tbta は以下のように表わされる。

)](exp)([exp2

)cos()](Re[

)](exp)([exp2

)cos()](Re[

00

00

bbb

aaa

titiB

tBtb

titiA

tAta

φωφωφω

φωφωφω

++++−=+=

++++−=+=

)(),( tbta の実数表示 )](Re[)],(Re[ tbta の積を計算すると、 bababa sinsincoscos)cos( m=± を利用し、

)]cos()2[cos(2

)](Re[)](Re[ 00babat

BAtbta φφφφω −+++=

となる。一方複素数表示から計算し、その後実数表示を取る（ )]()(Re[)()( tbtatbta → ）と、

)2cos()]()(Re[)]2(exp[)2exp()()( 0000 baba tBAtbtatiBAtiABtbta φφωφφωω ++=→++−=−=

となる。 )]()(Re[ tbta の実数部分を見ると、時間依存の関数の係数が異なる以外に、 )](Re[)](Re[ tbta の時間

によらない項 2/)cos(|| baAB φφ − が脱落している事がわかる。つまり、 )(),( tbta が複素数表示の関数の時、

)]()(Re[)](Re[)](Re[ tbtatbta ≠

である。この講義ノートでは、「指数関数で表わした物理量 )(),( tbta は、本当はその実数部分を扱うという約

束」なので、2つの物理量 )(),( tbta の正しい積は、その実数部分同士の積を行わなければならない。

28

次に、複素表示が便利な例を示す。

微分の比較:

)sin(

)]cos([)](Re[

0

0

a

a

tA

tAdt

dta

dt

d

φωω

φω

+−=

+=

　　　　　

一方

)sin(

)sin()()(Re

)exp()exp()](exp[)(

0

0

00

a

a

aa

tA

tiAitadt

d

tiiAitiAdt

dta

dt

d

φωω

φωω

ωφωφω

+−=

+−⋅−=

→

−−−=+−=

　　　　　　

このように、 )(ta が複素数表示の場合でも、複素数表示のままで微分した結果の実数部分は、 )(ta の実数部分

を微分した結果と等しい事が分かる：

)](Re[)(Re tadt

dta

dt

d=

さらに、以下のように、 )(),( tata IR を実数として、 )()()( tiatata IR += と置いてみれば、n階の微分でも、

その関係は成り立つ事が分かる。

)](Re[)(Re

)()]()(Re[

)()()(Re))()((Re

tadt

dta

dt

d

tadt

dtiata

dt

d

tadt

dta

dt

dita

dt

dtiata

dt

d

n

n

n

n

Rn

n

IRn

n

Rn

n

In

n

Rn

n

IRn

n

=

→

=+

=

+=

+

別の有用な例：正弦関数の積の時間平均（大変便利）

周波数が等しい2つの正弦関数 )cos()(),cos()( 00 ba tBtbtAta φωφω +=+= の積：（ πω 2=T ）

)cos(2

)cos()cos(1

)](Re[)](Re[ 00

000 ba

T

ba

BAdttBtA

Ttbta φφφωφω −=++= ∫

ここで横棒は時間平均を意味する。

さて、 )exp(),exp( 00 ba iBBiAA φφ −=−= であるから（ *B はB の複素共役）、

)cos(*]Re[

)](exp[)exp()exp(*

00

0000

bb

bbba

BAAB

iBAiBiAAB

φφ

φφφφ

−=→

−−=+−=

　　

よって、複素数表示された同じ時間依存性を持つ2つの関数 )](Re[)],(Re[ tbta の時間平均 )](Re[)](Re[ tbta は、

複素振幅AおよびB を用いて次のように表わされる：

*]Re[2

1)](Re[)](Re[ ABtbta = 。

この関係式は、電磁波のエネルギーの計算などで頻繁に利用する。光の分野では、複素数表示の関数 )(),( tbta な

どを扱うことが多い。物理量は実数であるから、複素数表示の物理量は、その実数部分を取り扱うと解釈するの

が一般である。複素数表示で記述した2つの物理量の積（内積、外積など）の時間平均では、上記関係を用いれ

ば、計算時間が短縮される。この講義ノートでは、もっぱらこの関係式を用いる。

29

ポインティング (Poynting) ベクトル電磁波のエネルギー

物質におけるMaxwell方程式を出発点にとり、空間を伝播する電磁波のエネルギーに関する関係式を導く。

物質中電磁波のエネルギーの議論では、非磁性の物質（ 0µµ = ）を想定し、電磁波のエネルギーを計算し、

真空中の電磁場のエネルギーとの関連を議論する。

Maxwell方程式のうちの2つの式

0rot,rot =∂∂

+=∂∂

−tt

BEi

DH

において、Eと Hrot とのスカラー積およびH と Erot とのスカラー積（内積 ][⋅ ）をとると、

tttt ∂∂

⋅−=⋅→∂∂

−=⋅+∂∂

⋅=⋅→+∂∂

=B

HEHB

EiED

EHEiD

H rotrot,rotrot 。

EHHE rotrot ⋅⋅ ---- を公式 BAABBA rotrot)(div ⋅−⋅=× を用いて計算すると、上の式から、次式を得る。

tt ∂∂

⋅+⋅+∂∂

⋅=⋅+⋅−=×B

HiED

EHEEHHE rotrot)(div----

BAAB

BA

rotrot

)()()()()()(

)()()(

)()()(

)()()()(div

⋅−⋅=

∂∂

−∂

∂−

∂∂

−∂

∂−

∂

∂−

∂∂

−∂

∂−

∂∂

+∂∂

−∂

∂+

∂

∂−

∂∂

=

∂∂

−∂

∂+

∂∂

−∂

∂+

∂

∂−

∂∂

+

∂

∂−

∂∂

+∂∂

−∂∂

+∂∂

−∂

∂=

−∂∂

+−∂∂

+−∂∂

=×

　　　　

　　　　　

　　　　　　　　

　　　　　

y

B

x

BA

x

B

y

BA

y

B

y

BA

y

A

y

AB

x

A

z

AB

z

A

y

AB

y

BA

y

BA

y

BA

y

BA

x

BA

x

BA

z

AB

z

AB

y

AB

y

AB

x

AB

x

AB

BABAz

BABAy

BABAx

xy

zzx

y

yzx

yxz

zxy

yxx

xy

y

xz

xx

z

y

zz

y

y

xx

yx

zx

xz

y

y

z

xyyxzxxzyzzy

ここで特に物質中の電磁波である事を明確にするため、 HBED µε == , （物質中において、 HBED ,,, の

間に成り立つ関係式）を代入する（誘電率ε と透磁率µ を用いる）と、

iEHEHE ⋅+

+∂∂

=× 22

22)(div

µεt

----

を得る。なおこの関係式を導くにあたり、次の定義および計算を行った。 ][⋅ は、ベクトルの内積を意味する。

2

,,

2

,,

,,

22

22

),,(),,(

EE

ED

E

EEEED

εεεε

ε

tE

tt

EE

tt

EEEEEEE

zyxj

j

zyxj

j

j

zyxj

jzyxzyx

∂∂

=∂∂

=∂

∂=

∂∂

⋅=∂∂

⋅→

=⋅=⋅≡=

∑∑

∑

==

=

　　　　　　

　　

2H の項についても、同様の計算で、2

2H

BH

µtt ∂

∂=

∂∂

⋅ を用いる。これを利用し、上の関係式を得る。

ここで、次の物理量

22

22HE

µε+

30

を、物質に蓄えられている電磁波のエネルギー密度（単位体積あたりの電磁波のエネルギー）という。係数を含

めて、この形が電磁波のエネルギーになる事は、以下に簡単に示す。

1個の荷電粒子の運動を考えよう（複数個の荷電粒子の場合でも、同じ結果になる）。質量m 、電荷q の荷電

粒子の速度をv とすれば、電場E、磁束密度B における運動方程式は、

BvEv

×+= qqdt

dm

で与えられる。右辺第2項は、荷電粒子に加わるローレンツ力である。v との内積 ][⋅ を取り、

)()(2

)( 2 BvvEvvBvEvv

v ×⋅+⋅=

→×+⋅=⋅ qqm

dt

dqq

dt

dm

外積の定義から、 Bvv ×, は互いに直交するので、 0)( =×⋅ Bvv となる。これを用いると、荷電粒子の電流を

vi q= とおいて、

EiEvv ⋅=+⋅=

0)(

2

2 qm

dt

d

を得る（複数個の荷電粒子の場合は、k 番目の荷電粒子の質量を km 、電流を kkk q vi = とおけばよい。）。この式

を、既に導いた式に代入すると、

+

+∂∂

=×→⋅+

+∂∂

=× 22222

222)(div

22)(div vHEHEiEHEHE

m

dt

d

tt

µεµε--------

となり、22

22HE

µε+ は、エネルギーの次元を持つ。また、その形（係数）も、これでよい事が分かる。

大雑把な話になるが、例えば水道ホースを流れる水との対比で、電磁波のエネルギーとパワーを考えよう。こ

の時、ホースの体積 Vδ にある密度ρ の水の質量は、 Vδρ ⋅ で与えられる。よって、単位体積あたりの水の質

量は、 ρδδρ =⋅ VV / で与えられる（正しくは、これが密度の定義である。）。この単位体積あたりの水の質

量に対応する量が、電磁波のエネルギー密度である。また、ホースの中から時間 t∆ 当たり出てくる水の質量は、

「流速v の水が時間 t∆ でホースを移動する長さ tv ∆⋅ 」と「ホースの断面積 S∆ 」の積で与えられる体積

Stv ∆×∆⋅ に、水の密度ρ をかけたもので与えられる。よって、単位時間・単位断面積当たり、ホースから流

れ出る水の質量は、 vStStv ρρ =∆⋅∆∆×∆⋅× )/( で与えられる。この量に対応するのが、単位面積あたりの

電磁波のパワーである。

さて、この式を3次元空間に渡って積分し、ガウスの定理 dSdVSV ∫∫ ⋅= nAAdiv を

適用する。ここにAは任意のベクトル、n はV を囲む表面S の単位法線ベクトルである。

結果は、以下のとおり。

∫∫∫ ⋅+

+∂∂

=⋅×−VVS

dVdVt

dS iEHEnHE22

22)(

µε

ここで、左辺は、S で囲まれた空間へ流入するエネルギーを意味する。 HES ×= をポインティングベクト

ルという。右辺第1項の体積積分は、「物質に蓄えられる電磁波のエネルギーの、単位時間あたりの変化量」

を表す。第2項の体積積分は、「運動する電荷（電流）によって電磁場が失う単位時間あたりのエネルギー（ジ

ュール熱の発生）」である。外部から流入した電磁波のエネルギーは、体積V の空間に電磁波のエネルギーと

n

V

31

して蓄えられる。また、運動する電子による熱の発生として消費される。

今度は、物質中にある分極（単位体積あたりの電気双極子）P の存在を明確にした計算を行う。先ほどと同様

の計算を HBPED 00 , µε =+= の関係式を用いて行う。用いる式は、以下のとおり。

∂∂

−=

+∂

+∂=

→

=∂∂

+

=∂∂

−

t

t

t

t

)(rot

)(rot

0rot

rot

0

0

HE

iPE

H

BE

iD

H

µ

ε

Eと Hrot とのスカラー積（内積 ][⋅ ）をとると、

iEP

EEiEP

EE

EHEiPE

H ⋅+∂∂

⋅+

∂∂

=⋅+∂∂

⋅+∂∂

⋅=⋅+∂

+∂=

ttttt

200

0

2rot:

)(rot

εε

ε。

となる。同様にH と Erot とのスカラー積（内積）をとると

∂∂

−=∂

∂⋅−=⋅

∂∂

−= 200

0

2rot:

)(rot H

HHEH

HE

µµ

µttt

。

これから、ベクトル公式 BAABBA rotrot)(div ⋅−⋅=× （p.29）を利用して、

iEP

EHEHE ⋅+∂∂

⋅+

+∂∂

=×tt

2020

22)(div

µε----

を得る。同じ物質に対して、 BD, の記述の仕方が異なるだけなので、先の計算と比較し、次の等式が成り立つ。

ttt ∂∂

⋅+

+∂∂

=

+∂∂ P

EHEHE202022

2222

µεµε

これは、単位体積あたりで考えた場合、「物質中に蓄えられる電磁場のエネルギーの単位時間当たりの変化量」

が、「(物質の部分を真空にしたとして）その真空の部分での電磁場のエネルギーの単位時間当たりの変化量」

と「単位時間当たり、電気双極子（分極）に吸収されるエネルギー」の和に等しい事を示す。

先ほどと同様にして、体積積分の計算をすると、

∫∫∫∫ ⋅+∂∂

⋅+

+∂∂

=⋅×−VVV

SdVdV

tdV

tdS iE

PEHEnHE 2020

22)(

µε

となる。右辺第1項 ( )2/2/)/( 2

0

2

0 HE µε +∂∂ t の体積積分は、「空間V （その空間の部分は真空と考える）に

蓄えられる電磁場のエネルギーの、単位時間あたりの変化量」である。第2項 t∂∂⋅ P/E の体積積分は、「体積V

の中の電気双極子（分極）の運動を電磁波が誘起する事で、電気双極子（分極）が吸収するエネルギーの単位時

間当たりの変化量」である。第3項 iE ⋅ の体積積分は、「運動する電荷によって電磁場が失う単位時間あたりの

エネルギー（ジュール熱）」である。左辺の被積分関数中のnは体積V の境界面で外向きである単位長さの垂線

ベクトルを示す。よってこの式は、V の外部から流入する電磁波のエネルギーが、それら3つのエネルギーに変

わる事を示す。

32

電磁波の単位面積当たりの平均パワー

一様平面電磁波により伝播方向に運ばれる単位面積当たりの平均パワーを求める。水道水をホースで流す場合

との比較で話しよう。ホースの中にある単位体積あたりの水の質量に対応するのが、電磁波のエネルギー密度で

ある。そして、ホースの中から単位時間当たり出てくる水の質量に相当する量が、電磁波のパワーになる。単位

時間当たりに出てくる水の質量をホースの断面積で割った量に対応するのが、単位面積あたりの電磁波のパワー

である。

次のような電磁場を考える。繰り返すが、以下の物理量（電場と磁場）の実数部分を考えよう。

)](exp[)0,,0()](exp[))(),(),((),(

))](exp[)0,0,()](exp[))(),(),((),(

00

00

kztiHkztitHtHtHtz

kztiEkztitEtEtEtz

yzyx

xzyx

−−=−−==

−−=−−==

ωω

ωω

HH

EE

真空中での電磁波のエネルギーの流れ

ここでは、まず真空中で考える。電場と磁場の振幅 00 ,HE の振幅の間に成り立つ関係式を求める。

Maxwell方程式の中の1つの式（次式）に、上の電場および磁場の式を代入する。

0),(

),(rot),(

),(rot 0 =∂

∂+=

∂∂

+t

tt

t

tt

xHxE

xBxE µ

その後、 y 成分の比較を行う事で、電場および磁場間の振幅の関係式を求める。（ 000 /1 εµ=c ）

( ) ( )

( ) kccZEZ

H

EEEkc

kE

kHHiikE

kztiHikztiikEt

tH

x

tE

z

tE

xy

xxxxyyx

yx

yzx

0000

2/1

0000

0

0

02/1

00

02/1

000

0

00

0

0

0000

0000

,/1,/,1

/

1

/

10)(

0)](exp[)(0)](exp[0)()()(

====→

====→=−+→

=−−−+−−−→=∂

∂+

∂∂

−∂

∂

ωεµεµ

εµεµµµωµωµ

ωωµωµ

　　　

　　

)(3770 Ω=Z は、自由空間（真空）のインピーダンスである。この関係式を電磁波のエネルギーの式に代入

すると（時間平均を取るための計算）、

2

00

2

00

2

0

00

02

002

002

00

000

00000

)/(2222)(

2)(

2*

2*

2

E

HHEEHHEE

εε

εµµεµεµεµε

==

+=+=⋅+⋅=⋅+⋅

x

xxyx

E

EEHE

　　　　　　　　　

よってその時間平均は、（ )exp()(),( titt ω−∝HE の時間依存性がある場合、p.28を見よ）

2

00002020

2*

2*

2Re

2

1

22EHHEEHE

εµεµε=

⋅+⋅=+

で与えられる。一様な平面電磁波により伝播方向に運ばれる単位面積当たりの平均パワーは、真空中では電磁

波の速さが 0c なので、以下の結果を得る。

2

0002

00

022

EEc

cεε

=×

今度は、ポインティング(Poynting)ベクトルから電磁波のエネルギーを計算しよう。

)](exp[)0,,0()](exp[),(

)](exp[)0,0,()](exp[),(

00

00

kztiHkztitz

kztiEkztitz

y

x

−−=−−=

−−=−−=

ωω

ωω

HH

EE

33

ポインティングベクトルのz成分を計算する（注意： yx, 成分＝ゼロ）。

[ ]

2

000

2

000

2

02/1

00

02

02/1

00

00

2

1

2

1

)(2

1

)/(

1

2

1

2

1*)()(Re

2

1

E

HEHES z

cEc

EEHE

x

xxyxyx

εε

εµε

εµ

==

====×=

　　

よって、真空中での単位時間、単位面積当たりの電磁波のエネルギーの流れ（電磁波のパワー）は、

2

0002

1Ecε

で与えられる。なお、ポインティングベクトル HES ×= の yx, 成分は、計算すればわかるように、ゼロで

ある。このように HES ×= は、電磁波のエネルギーの流れる方向と電磁波のパワー(単位断面積あたり単位

時間当たりの、電磁波のエネルギーの流れ)を表わす。

物質中での電磁波のエネルギーの流れ

電気伝導度 0=σ の物質（誘電率ε 、透磁率 0µ とする）を考え、屈折率を実数nとする。電場と磁場の振幅

との間には、次の関係がある。（真空中での電場と磁場の振幅の間に成り立つ関係式を求めたのと同様の計算で、

導く事ができる。前ページの計算を見よ。）

2/1

02/1

000

)/(,)/(

εεεµ

=== nEn

EZ

nH

この関係式を電磁波のエネルギーの式に代入する（時間平均を取るため、前もって必要な計算（p.28を見よ））。

なお、 *E （ *H ）はE（H ）の複素共役なベクトルである。物質は非磁性で 0µµ = とする。

2

0

2

0

2

0

00

002

0

2

002

0000

00)/(

/

2222)(

2)(

2*

2*

2

E

HHEEHHEE

εε

εµεεµεµεµεµε

==

+=+=⋅+⋅=⋅+⋅

x

xxyx

E

EEHE

　　　　　　　　

　

電磁波のエネルギー密度の時間平均を取ると、

2

0022

2*

2*

2Re

2

1

22EHHEEHE

εµεµε=

⋅+⋅=+

となる。一様な平面電磁波により、物質内で伝播方向に運ばれる単位面積当たりの平均パワーは、電磁波の速

さ（位相速度）が ncc /0= なので、（ n=2/1

0)/( εε にも注意）

2

000

2

000

2

00

2/1

00

2

00

2/1

0

2

02/1

0

0

2

002

002

0

2

1

2

1)/(

2

1)(

2

1

)/(2

222

E

E

ncEncEcEcEc

n

cE

n

cEc

xxxx

xx

εεεεεεεεε

ε

εεε

=====

==×

　　

となる。よって、物質中での、単位時間、単位面積当たりの電磁波のエネルギーの流れ（単位面積あたりの電

磁波のパワー）は、以下の式で与えられる。

2

000

2

00

2

1,

2

1EE nc

n

cε

ε

34

ポインティングベクトルを用いた電磁波のパワーの計算

今度は、ポインティングベクトルから電磁波のエネルギーを計算し、上と同じ結果を得よう。ポインティン

グベクトルの z 成分（電磁波のエネルギーの流れる方向は、 z 方向）を計算する。 )0,0),(( tEx=E 、

)0),(,0( tH y=H であたえられるので、 HE × は、 z 成分のみを持つ。 000 /1 εµ=c 、 0/εε=n の関

係式を用いると、

[ ]

( )( )

==⋅

==

====

=

===×=

−

2

002

002

02/1

000

2

02/12

0

2

02/1

0

2

000

2

000

2

00

2/1

00

2

02/1

0

2/1

00

02

02/1

000

2

02/1

0

0000

2

1

2

1

))/((2

1

)/(

1

2

1

)/(

1

2

1

2

1

2

1/

2

1

/)(2

1

))/((

1

2

1

)/(

1

2

1

2

1*)()(Re

2

1)(

E

E

HEHES

n

cE

n

cEEE

ncEncEcEE

EHE

xxxx

xxxx

xyxyxzz

εεεεεµ

εεεµεµ

εεεεεεεεµ

εεεεµ

εµ

　　

物質中での単位時間、単位面積当たりの電磁波のエネルギーの流れ（単位面積あたりの電磁波のパワー）は、

次式で与えられる。

2

000

2

00

2

1,

2

1EE nc

n

cε

ε

注意：一様な空間を伝播する電磁波では、平面波の伝播する方向（波数k ）とポインティングベクトルS の方

向は等しいが、空間的に一様でない物質（電場の感じる屈折率が、電場の振幅の向きによって異なる物質）中

では、一般に同じ方向ではない（ kS ∝ が成り立たない。）。「結晶中での電磁波の伝播」を参照。

分極による電磁波エネルギーの損失

電磁波に運動が誘起されることで、電気双極子（分極）が蓄える単位体積あたりのエネルギーの単位時間当

たりの変化量は、

t∂∂

⋅P

E

で与えられる。簡単のため、E//P と仮定し、ベクトル表示をしないで、スカラー表示 )(),( tPtE にする。そ

れらが複素表示で

)exp()exp()(),exp()( 0000 tiEtiPtPtiEtE ωχεωω −=−=−=

で与えられているとする。ここで、 χ は電気感受率である。さらに IR χχ , を実数として、電気感受率

IR iχχχ += とおくと、

( )[ ] [ ] [ ]

[ ] [ ]2

00

2

00

2

00

2

00000

||2

1

||)(Re2

1||)(Re

2

1

||)*(Re2

1**)(Re

2

1*)()(Re

2

1

E

EiEii

EiiEiEtPitEt

I

IRIR

IR

χωε

χχωεχχωε

χχωεχωεω

=

+=−=

+=+⋅=−=∂∂

⋅

　　　

　　　

PE

となる。 Iχ の値が大きいほど、電気双極子（分極）に電磁波が吸収されることが分かる。

35

物質中の電磁波の基本的性質

ここでは、電磁波の伝播方向や偏光方向に関係なく、誘電率ε および透磁率µ が一定の値を持つ媒質（等方

性媒質）での電磁波の伝播の様子をみる。媒質の電気伝導度 0=σ または 0≠σ であり、電荷密度 0)( =xρ は

共通とする。電磁波は横波であり、電場、磁場と伝播方向は、互いに直交することを理解する。

(i) 電気伝導度 0=σ の場合

電磁波の伝播方向や偏光方向に関係なく、パラメーターε およびµ が一定の値を持つ媒質（等方性媒質）での

電磁波の伝播を考えよう。また、電気伝導度 0=σ かつ、真電荷は存在しない（自由電子がない誘電体などの場

合：電荷密度 0)( =xρ ）とする。すなわち、ε とµ が定数、かつ電気伝導度 0=σ （ 0== Ei σ ）である等

方性の媒質中での電磁波の波動方程式を導き、電磁波の特性などを議論する。なお、光の周波数で振動する磁場

に対する物質の透磁率は，真空における値 0µ とほぼ等しい（例外は磁性体）。

マクスウェル方程式は、

0),(div

0),(div

0rot

0rot

=

=

=∂∂

−

=∂

∂+

t

t

t

t

xB

xE

E

H

HE

ε

µ

注意：

Ei

HB

ED

σµε

=

=

=

まず、 0rot =∂

∂+

t

HE µ でrot をとり、 0rot

EH =

∂∂

−t

ε を代入する。ラプラシアン2

2

2

2

2

2

zyx ∂∂

+∂∂

+∂∂

≡∆

0rotrotrotrotrotrotrotrotrot02

2

=∂∂

+=∂∂

∂∂

+=∂∂

+=

∂

∂+=

ttttt

EE

EEHE

HE µεεµµµ

恒等式： AAA ∆−= graddivrotrot （p.23を見よ）、および 0div =E を用いると、上式は

00graddiv2

2

2

2

2

2

EE

EE

EE =

∂∂

−∆→=∂∂

+∆−=∂∂

+∆−ttt

εµµεµε

となり、電場 Eの満たす方程式（波動方程式）が得られる。この式は、電場が等方性媒質中を伝播する事を示

す。同様にして、磁場 H においても、同様の計算で、波動方程式

0)(rotrotrotrotrotrotrot02

2

2

2

HH

EEHE

H =

∂∂

−∆→∂

∂−−=

∂∂

−=

∂∂

−=tttt

εµµεεε

が得られる。なお、これら電磁波の波動方程式において、変位電流tt ∂

∂=

∂∂ ED

ε の存在は重要である。もし、変位

電流が存在しないなら、上の2つの式はラプラス方程式（ 0 E =∆ 、 0 H =∆ ）になってしまう。すなわち、電磁

波が発生しない（電磁波の時間依存性が無い）。

36

電磁波は横波である事

一様な空間をk方向に伝播する、平面波の電磁波を考える。一定の周波数を持つ交流電磁波の場合は、電磁波

は横波である事を簡単に示す事ができる。電磁波をそれぞれ、

)](exp[),(

)](exp[),(

0

0

xkBxB

xkExE

⋅−−=

⋅−−=

tit

tit

ω

ω 、

),,(

),,(

0000

0000

zyx

zyx

BBB

EEE

=

=

B

E：定ベクトル

と書く。kは波数ベクトル、ω は電磁波の角周波数である。真電荷なし（自由電子がない誘電体などの場合：

0)( =xρ ）では、

0),(div0)()],(div[),(div 0 =→=== ttt xExxExD ρε

0

0)](exp[

)](exp[)(),(),(),(),(div

)](exp[)](exp[)](exp[)(),(

)](exp[div),(div

)),(),,(),,((),(

0

0

000

000

0

=⋅→

=⋅−−⋅=

⋅−−⋅+⋅+⋅=∂∂

+∂∂

+∂∂

=

⋅−−⋅=⋅−−∂∂

=⋅−−∂∂

=∂∂

⋅−−=

=

kE

xkkE

xkxxxxE

xkxkxkEx

xkExE

xxxxE

　　　

　　　　　　

　

　　

　

tii

tikEkEkEitEz

tEy

tEx

t

tiikEtiEx

tix

tEx

tit

tEtEtEt

zzyyxxzyx

xxxxx

zyx

ω

ω

ωωω

ω

同様にして、

00),(div 0 =⋅→= kBxB t

を得る。さらに、 0),(

),(rot =∂

∂+

t

tt

xBxE を利用して kBE ,, 00 の間に成り立つ関係式を求める。それぞれ

の項の計算をすると、以下のとおり。

)](exp[)(),(rot 0 tiit ω−⋅×= xkEkxE

)](exp[][

)](exp[)(

)](exp[)](exp[

),(),()],(rot[

0

00

00

xkEk

xk

xkxk

xxxE

⋅−−×=

⋅−−⋅−⋅=

⋅−−⋅−⋅−−⋅=

∂∂

−∂∂

=

tii

tikEkEi

tikiEtikiE

tEz

tEy

t

x

zyyz

zyyz

yzx

ω

ω

ωω

　　　　　

　　　　　

　　　　　

)](exp[)](exp[),(

00 xkBxkBxB

⋅−−+=⋅−−∂∂

−=∂

∂− tiiti

tt

tωωω

よって、

00

00 )](exp[)](exp[)(0),(

),(rot

BEk

xkBxkEkxB

xE

ω

ωωω

=×→

⋅−−+=⋅−−×→=∂

∂+

　　　

tiitiit

tt

は互いに直交するkBE

BEk

kB

kE

,,0

0

00

00

0

0

→

=×

=⋅

=⋅

ω

電磁波の伝播方向はkなので、電場、磁場の向き 00 ,BE は、伝播方向k と

直交する、横波である。

E0

B0 = k×E0/ω

k

37

光の偏光特性直線偏光、円偏光

自由空間や空間的に一様な媒質中を伝播する電磁波は、電磁波の振幅 00 ,BE が電磁波の伝播方向kと直交

する横波である事を示した。ここでは、電磁波の偏光に関する定義や、波数や屈折率との関係について述べる。

直線偏光

電磁波を、

)](exp[),(

)](exp[),(

0

0

xkBxB

xkExE

⋅−−=

⋅−−=

tit

tit

ω

ω 、一定:, 00 BE

とし、 00 ,BE を実数のベクトルとすれば、この電磁波は直線偏光しているという。光学の分野では、電場の方

向（ベクトル 0E の向き）を偏光方向とするのが、慣習になっている。

自然光（無偏光の光）

自然光は、 00 ,BE が時間的に一定でない、ランダムな偏光を

持つ。なお、偏光方向や振幅に何らかの空間的時間的な規則性

があれば、自然光（無偏光）とは言わない。

直線偏光板：linear polarizer

直線偏光板は、右図のように、偏光がランダムな向きの光が

ある場合、その光の、ある向きの直線偏光の光のみを通す。

直線偏光の光電場の偏光方向が、偏光板の向き（直線偏光の光を通す向き）に対して、図のように、角度θ

をなすとする。偏光板の「直線偏光の光を通す向き」と同じ偏光方向のベクトルを xe （ 1|| =xe ）、それと直交

する方向のベクトルを ye （ 1|| =ye ）とすると、電場ベクトル 0E は、次のように分解される：

|,sincos)()( 0000000 E|eeEeeEeeE =+=⋅+⋅= EEE yxyyxx θθ 。

光電場の強度 I は電場の振幅E の2乗（2E ）に比例する（p.30電磁波のパワー

を参照）。よって、偏光板を透過する光強度は、入射光強度を 0I とすれば、

θθθ 2

0

22

0 coscos)cos( IIEI =→∝∝

となる。

円偏光

0=z の場所で、

yx

yx

tEtEtz

EkztEkztEt

eeE

eexE

)sin()cos(),0(

0,)sin()()cos(),(

00

000

ωω

ωω

−==→

>−−+−=

で与えられる電場を考えよう。この電場の向きが時間t に対してどのよう

に変化するか、次の表を利用して調べる。

tω 0 2/π π 2/3π π2 2/5π π3

)cos( tω 1 0 －1 0 1 0 －1

－ )sin( tω 0 －1 0 +1 0 －1 0

自然光

直線偏光

θ

E0

ex

ey

ωt=0、2π

π/2,5π/2

π

3π/2

ey

ex

38

すると、図に示すように、電場の向きは、その大きさを変えることなく回転する。図では、場所を固定した場

合を示したが、電場が伝播する場合、螺旋（らせん）回転しながら進む。

円偏光：右周りの円偏光と左周りの円偏光の定義

例えば、コマが回転しているとしよう。ある人は、それを見て、右回りに回転していると言うかも知れない。

でも、反対（裏）側から見ると、左回りに回転している。同じ回転でも、見る方向により、回転方向は違う。

よって回転方向の明確な定義が必要になる。我々は、空間伝播しなが

らその電場ベクトルが回転する光電場がある事を、知った。その回転

方向の一般的な定義を見る。先ほど考察した電場を考えよう。 0=z の場所で、

yx

yx

tEtEtz

EkztEkztEt

eeE

eexE

)sin()cos(),0(

0,)sin()()cos(),(

00

000

ωω

ωω

−==→

>−−+−=

で与えられる電場を考えよう。この電場の向きは、先ほど調べたよう

に、表を利用すると、簡単に分かる。

tω 0 2/π π 2/3π π2 2/5π π3

)cos( tω 1 0 －1 0 1 0 －1

－ )sin( tω 0 －1 0 +1 0 －1 0

電場の伝播方向は、位相項 kzt −ω が一定になる事を要求すると、

00)(

>=→=−

kdt

dz

dt

kztd ωω

となり、 z+ の方向に伝播する事がわかる。 −∞=z から 0=z の方向を見よう。すると、 0=z での電場ベ

クトルの時間的な動きは、左回り（反時計周り）の回転になる。しかし、 +∞=z から 0=z の方向を見ると、

0=z での電場ベクトルの時間的な動きは、右回り（時計周り）の回転になる。

一般的な定義では、位置を固定（例では、 0=z の位置）して電場の伝播する向きと逆向き（例では +∞→z

から 0=z の方向）から電場の時間発展を見る時、電場ベクトルが右回りの回転であれば、右回りの円偏光と

いう。

円偏光の電場の様子を見る際、時間をある値に固定して、その位置依存性を見る方法もある。先ほどと同じ

電場で、簡単のため 0=t の場合を見よう。

yx

yx

kzEkzEtz

EkztEkztEt

eeE

eexE

)sin()cos()0,(

0,)sin()()cos(),(

00

000

+==→

>−−+−= ωω

先ほどと同様の表をkz に対して作ると、下のようになる。

kz 0 2/π π 2/3π π2 2/5π π3

)cos(kz 1 0 －1 0 1 0 －1

)sin(kz 0 1 0 －1 0 1 0

この場合、電場の伝播方向に向かって（この例では +∞=→−∞= zz の方向に進む）電場ベクトルの回転方

向を見る（時間t は一定とする）と、右回りで電場ベクトルが回転しながら進む。当然ながら、時間と位置の

いずれを固定してその様子を見ようとも、同じ電場であるなら、同じ呼び方をされなければならない。よって、

ωt=0、2π

π/2,5π/2

π

3π/2

ex

ey

39

時間を固定した場合、電場が伝播する方向に向かって電場ベクトルが

右回りに回転すれば、その電場は、右回りの円偏光という。右周りの

円偏光と逆向きの回転をする電場は左回りの円偏光という。右図は、

ある時間における、 z+ 方向へ伝播する右回り偏の円光の振幅ベクト

ルの変化する様子を示す。

波数k と周波数ω

以下では、次式で表わされる、周波数ω の単色光の平面波である電磁波を考える。物質（媒質）の誘電率は

ε で透磁率はµ とする。電場および磁場を平面波で次のように記述する。なお、 0=σ より、 0== Ei σ 。

)](exp[),()](exp[),( 00 xkBxB,xkExE ⋅−−=⋅−−= tittit ωω

Maxwell方程式から、 µεω ,,,k の間に成り立つ関係式を求める。

0rot =∂∂

+t

BE でさらにrotをとり、 0rot0rot

EB

EH =

∂∂

−→=∂∂

−tt

εµε を代入する。

0rotrotrotrotrotrotrot02

2

E

EBEBE =∂∂

+=∂∂

+=

∂∂

+=ttt

εµ

さて、

)](exp[)]([),(rotrot 0 xkEkkxE ⋅−−××−= tit ω

)](exp[)]([)](exp[)]([),(rotrot

)](exp[)(),(rot

)](exp[][

)](exp[)()](exp[)](exp[

),(),()],(rot[

00

0

0

0000

xkEkkxkEkkxE

xkEkxE

xkEk

xkxkxk

xxxE

⋅−−××−=⋅−−××=→

⋅−−×=

⋅−−×=

⋅−−⋅−⋅=⋅−−⋅−⋅−−⋅=

∂∂

−∂∂

=

titiiit

tiit

tii

tikEkEitikiEtikiE

tEz

tEy

t

x

zyyzzyyz

yzx

ωω

ω

ω

ωωω

　　

　　　

　　　

一方、

)]](exp[)](exp[)( 0

2

0

2

2

2

xkE xkE E

⋅−−−=⋅−−−=∂∂

titiit

ωεµωωωεµεµ

これより、

0

2

0

0

2

0

)(

0)](exp[)](exp[)]([

E Ekk

xkE xkEkk

εµω

ωεµωω

−=××→

=⋅−−−+⋅−−××−

　　

titi

が成り立つ。次にベクトルの外積の公式： )()()( BACCABCBA ⋅−⋅=×× （説明は、次ページ）およ

び 0),(div =txE から得られる式 00 =⋅ kE を利用し、 )( 0Ekk ×× を簡単にする。

00)](exp[

),(),(),()](exp[div),(div

00

0

=⋅→=⋅−−⋅=

∂∂

+∂∂

+∂∂

=⋅−−=

kExkkE

xxxxkExE

tii

tEz

tEy

tEx

tit zyx

ω

ω

　　　　　

すると、

x

y

z

伝播方向

（時間固定）

振幅ベクトルのz方向依存性

右回り円偏光

40

2

0

2

00

2

000

2

0

)(

)()()()(

εµω

εµωεµωεµω

=⋅→

−=⋅−→−=⋅−⋅×→−=××

kk

E EkkE kkEEkkE Ekk

この式で、右辺 02 >εµω であるから、 kk ⋅ は実数である。そこで 1||,// =kk eek となる単位ベクトル ke を

用いて kkek = と書けば、波数k と周波数ω の間の関係を表わす、次の式が得られる。

22 εµω kk ==⋅ k 。

波数k と屈折率n

媒質が非磁性体であるとして、 0µµ = とする。真空での電磁波の速度000 /1 µε=c を用いると、

0000

2

0

2

0

2

0

00

22

0

222 ,1

εεωω

εε

ωεε

ωεε

µεωεµεµω ==⋅=→==→=→= nc

nc

kc

kkk

ここで定義した0/εε=n を、物質の屈折率という。電磁場は、p.24 で示したように、 k/ω の（位相）速度

で伝播する。上の結果から、 nck // 0=ω である。

)()()( BACCABCBA ⋅−⋅=×× の説明

x

zzyyxxxzzyyxxx

zzyyxzzyyx

zxxzzxyyxyyzzyx

BACABACCACACAB

BACACCACAB

CBCBACBCBAAA

)]()([

)()(

)()(

)()()()()]([

BACCAB

CBCBCBA

⋅−⋅=

++−++=

+−+=

−−−=×−×=××

　　　　　　

　　　　　　

　　　　　　

他の成分も、同様に成り立つことが示される。

(ii) 電気伝導度 0≠σ の場合

0≠σ の場合、物質は電磁波を吸収する性質を持つ。ここでは、複素屈折率や吸収係数を説明する。

0=σ の場合と同様にして、電場が平面波： )](exp[0 xkEE ⋅−−= ti ω で 0≠σ の場合、 Ei σ= 、

ED ε= の関係から、

ωσ

εεεωεωσ

εωσωε it

tiiii

tt +≡

∂∂

=−=+−=+−=+∂∂

= ',),(

'')()(),(rotxE

EEEEiD

xH

を得る。 'ε を、複素誘電率という。こうする事で、物質の電磁波に対する応答が 'ε に含まれた。この式を用い

て、物質中の電磁波の伝播を考察する。注意：真電荷は存在しない（ 0)( =xρ ）とする。

0)( =xρ の理由

電気伝導度σ 、誘電率ε の物質中での電荷分布を考える。

電荷保存則 0),(div),(

=+∂

∂t

t

txi

xρに、 ρεσ === DEDEi div,, を代入すると、

0),(),(

),(div1),(

),(div),(

),(div),(

=+∂

∂=⋅+

∂∂

=+∂

∂=+

∂∂

tt

tt

t

tt

t

tt

t

tx

xxD

xxE

xxi

xρ

εσρ

εσ

ρσ

ρρ

これから、

41

−= ttεσ

ρρ exp)0,(),( xx

を得る。この式は、いかなる電荷分布があっても、 σε / 程度の時間で、電荷が散逸してしまう事を意味する。

そして、何らかの方法で電荷を集めたり、電荷を注入したり（電流 0),( ≠txi ）しなければ、電荷は散逸した

ままである。ここでは 0),( =txρ が常に成り立つものとし、

0),(div =txD

の扱いをする。そこで、Gauss の法則から得られる式は、 0),(div,0),(div == tt xDxB とする事にする。

この時マクスウェル方程式は以下のように、 0=σ の場合とほぼ同じ式になる。物質の透磁率 0µµ = とした。

ωσ

εεµε it

t

tt

t

+=

=

=

=∂∂

−

=∂∂

+

→

=

=

=∂∂

−

=∂∂

+

',

0div

0div

0'rot

0rot

0div

0div

),(rot

0rot

0

B

E

E

B

BE

B

E

iD

xH

BE

これから、 0=σ の場合と同様の計算で、次の波動方程式を得る。

2

2

2

2

2

2

2

2

02

2

0 0,'0,'zyxtt ∂

∂+

∂∂

+∂∂

≡∆=

∂∂

−∆=

∂∂

−∆ ラプラシアン　　　 B E µεµε

また、 0=σ の場合と同様の計算をする事で、 00 ,, BEk が互いに直行することがわかる。具体的には、次の3つ

の式から、 00 ,, BEk 間の直交性が導かれる。

は互いに直交するkBE

BEk

kB

kEB

EBE ,,:0

0

0rot,0div,0div 00

00

0

0

=×

=⋅

=⋅

→=∂∂

+==

ωt

よって、 0)( =xρ である吸収媒質（ 0≠σ ）では、電磁波は横波である。

波数 "i' kkk += と周波数ω

以下では、次式のような、周波数ω の単色光の平面波である電磁波を考える。

)](exp[),()](exp[),( 00 xkBxB,xkExE ⋅−−=⋅−−= tittit ωω （2-24

Maxwell方程式 0'rot,0rot 0 E

BB

E =∂∂

−=∂∂

+tt

µε でさらにrotをとる。

0'rotrotrotrotrotrotrot02

2

0 E

EBEBE =∂∂

+=∂∂

+=

∂∂

+=ttt

µε

ここで、 0=σ の場合に行った計算と同様の計算をする。

)](exp[)]([)](exp[)]([),(rotrot

)](exp[)(),(rot

)](exp[][),(),()],(rot[

00

0

0

xkEkkxkEkkxE

xkEkxE

xkEkxxxE

⋅−−××−=⋅−−××=

⋅−−×=

⋅−−×=∂∂

−∂∂

=

titiiit

tiit

tiitEz

tEy

t xyzx

ωω

ω

ω

)](exp[')](exp[)('' 0

2

00

2

02

2

0 xkE xkE E

⋅−−−=⋅−−−=∂∂

titiit

ωωµεωωµεµε

42

これより、

0

2

00

0

2

00

')(

0)](exp[')](exp[)]([

E Ekk

xkE xkEkk

ωµε

ωωµεω

−=××→

=⋅−−−+⋅−−××−

　　

titi

ベクトルの外積の公式： )()()( BACCABCBA ⋅−⋅=×× （p.40に説明）と 00),(div 0 =⋅→= kExE t を

用いると、

0

2

000

2

0000

2

00 ')(')()(')( E EkkE kkEEkkE Ekk ωµεωµεωµε −=⋅−→−=⋅−⋅×→−=××

すなわち、

2

0' ωµε kk =⋅

となる。右辺は複素数（ ωσεε /' i+= ）なので、左辺のkも複素数である。よって、実数のベクトル ",' kk

を用いて "i' kkk += と表わすと、

)xkxkEE ⋅−⋅−−= "'ti exp()](exp[0 ω

となる。 0"≠k の場合、電磁波は伝播に伴い、振幅が変化（減衰または増大）する。

波数kと複素屈折率n~

吸収のない（ 0=σ ）物質の場合に、実数の波数k ベクトルを単位ベクトル ke （ 1|| =ke ）を用いて kkek =

と書いた。そして、誘電率ε 、透磁率 0µ の媒質を伝播する電磁波の波数（大きさk ）が、電磁波の周波数ω と

次のような関係を持っている事を示した。ここで、 000 /1 µε=c 、屈折率 0/εε=n である。

000

2

0

2

0

2

00

0

22

0 )(c

nc

kc

kωω

εεω

εε

ωµεεε

ωεµ =⋅=→==→=⋅ kk

この関係式をk が複素数のベクトルの場合に拡張する。簡単のため "//' kk とし、 1||,// =kk eek となる単位

ベクトル ke を用いて kikk ek )"'( += と書こう。そうすれば、波数 kk ikkk ekek )"'( +=→= の変化に伴い、

屈折率も κinnn +=→ ~ のように変わると予想される。すなわち、以下のような関係式になるだろう。

ωσ

εεκωω

εεω

iinc

ncc

ikkk +≡+===+= '),(~'"'

0000

---（1）

ここでn~は、複素屈折率と呼ばれる。また、nは屈折率であり、κ は消衰係数である。屈折率nは、 0=σ の

媒質と同様、電場の位相速度の遅くなる程度を表わす定数であり、消衰係数κ は、電場の伝播に伴って、電場

の振幅が変化する度合いを表わす定数である。 ωσεε /' i+≡ の場合の、 κ,n の具体的な式は、次の「吸収係

数α 」の項目で、与えられている。

吸収係数α ：Beerの法則

光の伝播方向に沿っての、物質による光吸収を考える。物質の吸

収係数α は、単位長さあたり物質により光が吸収される割合として

定義される。すると、長さL の物質により光が吸収される割合は、

Lα になる。

光が z 軸方向に伝播し、物質の表面から深さ z の場所での光強度

を )(zI とする。光の位置が z から dzz + だけ変化する時、光が吸 z=0 z z+dz

I(z) I(z+dz)

43

収される割合は dzα で与えられる。よって、光の位置 dzzz +, での光強度の関係は、次のようになる。ここ

で、 )0(I は入射光強度である。

)exp()0()(

)()(

)()()(

zIzI

zIdz

zdIdzzIzIdzzI

α

αα

−=→

−=→−=−+

　

この関係をBeerの法則という。

先ほど得られた関係式 )(0

κω

inc

k += を用いると、物質中の電場ベクトルは、 κ,n を用いて、

( )

−−−=

+−−=−−= z

cz

c

ntizin

ctikztitz κ

ωωκ

ωωω

00

0

0

00 )(exp)(exp)(exp),( EEEE

となる。n は光電場の位相速度を決め、一方κ は振幅の大きさが電磁波の伝播とともに減衰（増大）する割合

を決めることがわかる。n およびκ をそれぞれ屈折率および消衰係数と呼ぶ。さらに光強度は電場の 2乗に比

例するので、光強度 )(zI は

)exp()0()( zIzI α−=

のように変化する。ここで0

2

c

ωκα = を吸収係数と呼ぶ。

κinn +=~ を ωσε ,, を用いて表わす事

さて、既に示した式（p.42の（1）式）を用いると、

=

−=→+−=+=+=→+=

κεε

κεεκκκ

εε

εε

εε

κεε

n

ninniniin

I

RIR

2/'

/'2)(

''''

0

22

0222

0000

の関係式が出る。これとωσ

εεεε ii IR +=+= ''' の関係式から、 κ,n を σωε ,, で表わすと、次のようになる。

2/12

0

2/12

0

112

,112

++−=

++=εωσ

εε

κεωσ

εε

n 。

ただし 0, >κn とした。また、波数ベクトル )("'0

κω

inc

ikkk +=+= は、ωσ

εεεε ii IR +=+= ''' の関係から

2/12

00

2/12

00

112

",112

'

++−=

++=εωσ

εεω

εωσ

εεω

ck

ck 　

で与えられる。計算の詳細は、以下。

0, >κn の値の求め方（ 0',' >IR εε としての計算） 22 , κ−n は、以下の2次方程式の解である。

44

04

''

)())2/('(

)(/'

2/'

/'2

0

2

0

2

222

0

22

0

0

22

0 =−−→

−⋅=−

−+=→

=

−=

εε

εε

κεε

κεε

κεε

κεεIR

I

R

I

R xxn

n

n

n

この方程式を解く。

2

0

2

2

0

2

0

2

0

2

2

0

22

0

2

0

2

2

0

22

0

2

0

2

0

2

''

2

1

2

'

''

4

1

2

'0

4

'

4

'

2

'0

4

''

εε

εε

εε

εε

εε

εε

εε

εε

εε

εε

εε

IRR

IRRIRRIR

x

xxxx

+±=→

+=

−→=−−

−→=−−

　　

よって、 0,0 >> κn として、

2/12

0

2

00

2

0

2

2

0

2

0

2

2

0

2

2

0

2

0

2

2/12

0

2

00

2

0

2

2

0

2

0

2

'''

2

1''

2

1

2

'''

2

1

2

'

'''

2

1''

2

1

2

'

+

+−=→++−=→+−=−

+

+=→++=

εε

εε

εε

κεε

εε

εε

κεε

εε

εε

κ

εε

εε

εε

εε

εε

εε

IRRIRRIRR

IRRIRR nn

ωσ

εεε == ',' IR を代入して、 κ,n 等の値を得る。

45

(iii) 空気と物質（ 0≠σ ）との境界面での反射率R

空気中から垂直入射で電磁波が物体へ入射する時、電磁波の反射率R は、物体

の複素屈折率をn~ として、次式で与えられる。

22

222

)1(

)1(

1~1~

κκ

+++−

=+−

=n

n

n

nR

以下、R の求め方

0=z を境界面とし、 0<z で空気中または真空中、 0>z で物質中とする。

電荷密度 0)( =xρ とする。また、 0µµ = の物質とする。電磁波は、境界面に

垂直入射する。

入射電磁波：incident z+ 方向、波数はk ベクトルの向きの、計算での確認

)](exp[)0,,0()](exp[),(

)](exp[)0,0,()](exp[),(

00

00

kztiHkztitz

kztiEkztitz

i

y

i

i

x

i

−−=−−=

−−=−−=

ωω

ωω

HH

EE

i

i

)0,/,0(/

)0,,0(

),0,0(),0,0,(

000

000

00

ωω

ω

kE

kE

kE

i

ii

iii

=×=→

==×

==

EkB

BEk

kE

　

反射電磁波：reflected z− 方向、波数は k−

( )[ ]( )[ ]kztiHzktitz

kztiEzktitz

r

y

r

r

x

r

+−=−−−=

+−=−−−=

ωω

ωω

exp)0,,0()])((exp[),(

exp)0,0,()])((exp[),(

00

00

HH

EE

r

r

)0,/,0(/

)0,,0(

),0,0(),0,0,(

000

000

00

ωω

ωrrrr

rr

rr

kE

kE

kE

−=×=→

=−=×

−==

EkB

BEk

kE

r

r

　

透過電磁波：transmitted z+ 方向、波数は 'k

)]'(exp[)0,,0()]'(exp[),(

)]'(exp[)0,0,()]'(exp[),(

00

00

zktiHzktitz

zktiEzktitz

t

y

t

t

x

t

−−=−−=

−−=−−=

ωω

ωω

HH

EE

t

t

)0,/',0(/

)0,',0(

)',0,0(),0,0,(

000

00

00

ωω

ωtttt

tttt

ttt

Ek

Ek

kE

=×=→

==×

==

EkB

BEk

kE

0

　

物質中のマクスウェル方程式は、 0µµ = かつ損失がある場合（ 0≠σ ）に行った計算から、次のような式

になった。（ 0),( =txρ を仮定する。 ωσεε /' i+= 。p.41の式）

0),(div

0),(div

),('),(rot),(

),(),(rot

0),(

),(rot

=

=∂

∂==

∂∂

−

=∂

∂+

t

t

t

ttt

t

tt

t

tt

xD

xB

xExHxi

xDxH

xBxE

εまたは

また、空気中と物質中で、次の関係式が成り立つ。

),(),(,,:0

0),(,,:0

0

00

ttz

tz

xExiEDHB

xiEDHB

σεµ

εµ

===>

===<

また、誘電率ε 、透磁率 0µµ = 、電気伝導度σ （複素誘電率 ωσεε /' i+= ）である物質を伝播する電磁波

の波数 'k （複素数）および物質の複素屈折率n~ は、

ωσ

εεκεεω

κωω

εεω

εε

ωµεεε

ωµε iinnc

inc

nc

kc

k' 2 +=+==+====== ',~',)(~'

':''

'00000

2

0

2

0

2

00

0

2

0

E0

B0=k×E0/ω

k

入射

反射

透過

nn ~1　　　　=

ｚ

z=0 Exi

Exr

Ext Hy

i

Hyr

Hyt

46

で与えられる（p.42の（1）式）。ここでは、 µσε ,, の代わりに、もっぱら κ,,~ nn を用いて、入射、反射、透

過する電磁波の振幅の間の関係式を求め、反射率R を求める。

最初は、Maxwell方程式の中の1つの式 0),(

),(rot),(

),(rot 0 =∂

∂+=

∂∂

+t

tt

t

tt

xHxE

xBxE µ に、

)](exp[)0,,0()](exp[),(

)](exp[)0,0,()](exp[),(:,

00

00

kztiHkztitz

kztiEkztitz

y

x

−−=−−=

−−=−−=

ωω

ωω

HH

EEHE

を代入し、 y 成分の比較を行う事で、電場および磁場間の振幅の関係式を求める。

まず、透過する電磁波では、（ ),(),,( tztz HE で、 'kk → の読み替え） y 成分を比較すると、

( ) ( )( ) 2/1

0000

0

02/1

00

02/1

000

0

00

0

0

0

0000

/,~

/

~

/

~~'

)]'(exp[)(0)]'(exp['

εµεµεµµµωµ

ωωµωµ

======→

−−−−=−−−→∂

∂−=

∂∂

−∂

∂

ZEZ

nE

nE

nE

c

nE

kH

zktiHizktiEikt

H

x

E

z

E

t

x

t

x

t

x

t

x

t

x

t

y

t

y

t

x

t

yt

z

t

x

)(3770 Ω=Z を、自由空間（真空）のインピーダンスという。

入射電磁波i

y

i

x HE 00 , の場合は、上の関係式で、 kk →→ ',' 0εε の置き換えを行い、

i

x

i

y EZ

H 0

0

0

1=

を得る。

最後に、反射する電磁波r

y

r

x HE 00 , では、伝播方向が z− だから波数は k− となる。そこで ),(),,( tztz HE

で、波数 kk −→ の読み替えを行い、真空中での伝播である事から誘電率 0' εε = （ 0,0 == σεε ）とおき、

( ) 2/1

0000

0

0

0

0

0000

/,1

)])((exp[)(0)])((exp[)(

εµωµ

ωωµωµ

=−=−

=→

−−−−−=−−−−−→∂

∂−=

∂∂

−∂

∂

ZEZ

Ek

H

zktiHizktiEkit

H

x

E

z

E

r

x

r

x

r

y

r

y

r

x

r

yr

z

r

x

となる。まとめると、

t

x

t

y

r

x

r

y

i

x

i

y EZ

nHE

ZHE

ZH 0

0

00

0

00

0

0

~,

1,

1=−==

である。

t

x

t

y

r

x

r

y

i

x

i

y HEHEHE 000000 ,,,,, 間の関係式は、次の2つのマクスウェル方程式を利用して電磁波の境界面で

の連続の条件を求め、その結果を利用して導く事ができる。

),(),(),(

),(rot,0),(

),(rot ttt

tt

t

tt xExi

xDxH

xBxE σ==

∂∂

−=∂

∂+

最初に、

0),(

),(rot =∂

∂+

t

tt

xBxE

において、空気中と物質との境界線を含む、長さ hl δδ , 面積 hlS δδδ = の微小な長方形で、面積積分を行う。

47

その際、ストークスの定理 ∫∫ ⋅=⋅SC

dStdt nxErrE ),(rot),( を用い、電場の面積積分を線積分に変える。

∫∫

∫∫∫

⋅∂

∂−=⋅→

⋅∂

∂−=⋅=⋅→=

∂∂

+

SC

SSC

dSt

tdt

dSt

tdStdt

t

tt

nxB

rrE

nxB

nxErrExB

xE

),(),(

),(),(rot),(0

),(),(rot

　　

ベクトルt は、境界面と平行な任意方向の単位ベクトルである。nを境界面に垂直な法線方向の単位ベクトル

とし、bをnおよびt に垂直な単位ベクトル（ tnb ×= ）とする。 1E を空気中の電場、 2E を物質中の電場

とする。

Eの線積分を具体的に計算しよう。右図の線積分の長方形の経路で、境界面に垂直方向の経路では、空気中

（領域1）で 1hδ 、物質中（領域2）で 2hδ の移動の長さとし、 21 hhh δδδ += とする。境界線に沿った経路の

長さは lδ とする。境界線をまたぐ 2 つの線積分は、相殺してゼロになる

（ lh δδ << として、 hδ の線積分の寄与を無視できると考えても、同じ結果を

得る。）。 ),( txE および ),( txB の値は、 lh δδ , の大きさが十分小さいとして、

長方形の経路内の 1点 xで代表させる。電場Eは線積分、磁束密度B は面積

積分を行えば、

∑=

⋅∂

∂−=−⋅+⋅

2,1

21

),()(

j

j

jhl

t

tll δδδδ b

xBtEtE

をえる。ここで 021 →+= hhh δδδ とすれば、右辺はゼロになるので、

0)(0)( 2121 =⋅−→=−⋅+⋅ tEEtEtE ll δδ

を得る。この式は、 21 EE − のt 成分がゼロ、すなわち、tr EEEEE i =+= 21 , のt 成分が等しい事を示す。

さて、t をx 方向にとれば、 21 EE − のx 方向の成分はゼロなので、

t

x

r

x

i

x EEE 000 =+

を得る。同様にして、

),(),(

),(rot tt

tt xE

xDxH σ+

∂∂

=

において、空気中と物質との境界線を含む、長さ hl δδ , 面積 hlS δδδ = の微小な長方形における面積積分を行

う。 1H を空気中の磁場、 2H を物質中の磁場とする。先ほどと同様に考え、ストークスの定理を用いると、

∑

∫∫

=

⋅

+

∂

∂=−⋅+⋅→

⋅

+∂

∂=⋅→+

∂∂

=

2,1

21 ),),(

)(

),(),(

),(),(),(

),(rot

j

jj

j

SC

hltt

tll

dStt

tdtt

t

tt

δδσδδ

σσ

b(xExD

tHtH

nxExD

rrHxExD

xH

　　

となる。ここで 021 →+= hhh δδδ とすれば、

0)(0)( 2121 =⋅−→=−⋅+⋅ tHHtHtH ll δδ

を得る。この式は 21 HH − のt 方向の成分がゼロ、すなわち、tr HHHHH i =+= 21 , のt 成分が等しい事

を示す。そこでt を y 方向にとれば、

δl

t δh

領域1

領域2

tnb ×=

t

n

48

t

x

r

x

i

x

t

x

r

x

i

x

t

x

r

x

i

x EnEEEZ

nE

ZE

ZHHH 0000

0

0

0

0

0

000~

~11=−→=−→=+

を得る。反射率は、入射電磁波強度に対する反射電磁波強度の比で定義される。また、電磁波の強度は電磁波

の振幅の2乗に比例するので、以下のようになる。

2

0

000

000

000

1~1~

1~1~

0)1~()1~(~ +

−=→

+−

−=→=++−→

=−

=+

n

nR

n

n

E

EEnEn

EnEE

EEEi

x

r

xr

x

i

xt

x

r

x

i

x

t

x

r

x

i

x

これより反射率R は、

2

1~1~

+−

=n

nR

で与えられる。磁場の連続の式から計算しても、同じ結果が得られる。以下は、その途中計算の部分である。

1~1~

~1

~10)~1()~1(

~1

0

000

000

000

+−

=+−

−=→=++−→

=−

=+

n

n

n

n

H

HHnHn

Hn

HH

HHH

i

x

r

xr

x

i

xt

y

r

y

i

y

t

x

r

x

i

x

補足：

0),(div,0),(div == tt xDxB

この2つの式から、境界面において、 BD, の境界面に垂直な成分（法線成分）が連続である事が求まる。

図のような、境界面を含む、厚さの薄い直方体の中で 0),(div =txB 、および 0),(div =txD の体積積分を

行う。Gaussの定理 dSdVSV ∫∫ ⋅= nAAdiv を利用し、nを境界面の法線方向の単位ベクトルとすると、

0)(

0)(,0

21

21212211

=⋅−→

=⋅−→−===⋅+⋅

nDD

nDDnnnnDnD SSS δδδ

これは、 21,DD のn方向の成分が等しい事を示す。同様の計算で、

0)( 21 =⋅− nBB

を得る。この場合も、 21,BB のn方向の成分が等しい事を示す。

n2

n1 n

49

物質の光学的性質

物質の光学的性質は、誘電率や吸収係数などで特徴付けられる。そのメカニズムを古典的なモデル（模型）

により説明する。代表的な模型には、電子がポテンシャルに束縛された電子振動模型と、自由電子の運動を扱

うドルーデ模型がある。

振動の一般論・概論

原子の振動子

振動している電気双極子は、1887年にヘルツ（Hertz）により実験的に証明されたように、電磁波を放射す

る。原子を振動する双極子とみなす考えは、Hertzの実験に先立って、1878年にローレンツ（Lorentz）によ

り提案された。

原子の振動子模型では、次のように考える。電子は、原子核の周りの安定な軌道を回る。そして、ばねは

電子の平衡状態からのずれを元に戻す復元力を表わす。負に帯電した電子（図で白丸）と、正に帯電した原子

核（図で黒丸）により、電気双極子を作る。その大きさは、電子と原子核との距離に比例する。

電子の質量、原子核の質量をそれぞれ Mm, とし、復元力のばね定数を SK とすれば、周波数 0ω は、

Mmmm

K

redred

S 111,0 +==ω

で与えられる。ここで redm は換算質量という。原子には、多くの遷移周波数

が存在する（近赤外、可視、紫外の領域に、多くの吸収と放出のスペクトル

が存在）ので、我々は、多くの双極子の存在を仮定しなければならない。

双極子とスペクトルとの関係は、双極子の振動を考える事で理解できる。

原子核の質量は電子に比べはるかに大きいので、原子双極子の振動の最中、ほとんどその位置を変えない。一

方電子は周波数 0ω で原子核に近づいたり遠ざかったりする。よって、電子の振動で時間変化する双極子が生

じる。 e+ の電荷を持つ原子核と e− の電荷を持つ電子との距離を )(tx （原子核が十分重くほとんど移動しな

いと考え、 )(tx は電子の平衡状態の位置からのずれと考えてよい。）である。とすると、時間変化する電気双

極子 )(tp は以下のように書かれる：

)()()( txetp −= 。

振動する双極子は、古典的な Hertz の双極子の実験が示すように、周波数 0ω の電磁波を放出する。よって原

子は、十分なエネルギーが与えられれば、その双極子固有の周波数の電磁波を放出する。

イオン・極性分子による振動子

もし媒質がイオン結晶なら、正負のイオンを含む。これらの帯電したイオンが平衡状態の位置から振動する

なら、振動する双極子になる。これは、原子核の周りで振動する電子と同じである。よって、イオン結晶中の

イオンと電磁波との相互作用を考える時、これら双極子による光学的効果を考える必要がある。

古典的な分極した分子では、2 つのイオンがばねにより結合し、安定な状態にある。2 つのイオンは、安定

な位置の周辺で振動し、振動する電気双極子を形成する。これは、原子核に束縛された電子の振動と同じであ

る。この模型の換算質量は、先ほどの原子と電子での双極子の換算質量よりもはるかに大きいので、電子が振

動する場合と比較すると周波数が低くなり、振動の周波数は近赤外の領域の周波数になる。

50

分子の振動と光との相互作用は、分子に対して外部から加える電場

によって生じる。電場により帯電した原子（イオン）が動くので、分子

が帯電している事は重要である。分極した分子では、電子雲が複数の原

子のうちの1つに高確率で位置している。

原子による電磁波の吸収

双極子模型は、原子と周波数ω の外部電磁波との相互作用を理解する助けにもなる。交流電場は、電子と原

子核に対して力を加え、周波数ω で振動させる。もしω が原子の持つ周波数のうちの１つ 0ω に等しいなら、

共鳴現象が起こる。共鳴により、電子と原子核の振動は非常に大きくなり、外部電場から原子へとエネルギー

が移動する。つまり 0ωω = なら、原子は電磁場のエネルギーを吸収する。

媒質による電磁波の吸収

量子論によると、電磁波の吸収では、原子が電磁波を吸収し、励起状態に遷移する。これは、終状態、始状

態のエネルギーを 12 , EE とし、 12 EE −=ωh の時に限って起こる。ここでhはh をプランク定数として、

)2/( πh=h である。もし原子が励起されたら、原子は無輻射遷移（電磁波を出さない遷移）により基底状態

に戻る。この際、吸収された電磁波は熱になる。または、ある時間遅れを伴い、電磁波を放出する。原子から

放出された電磁波は、互いにコヒーレンスが無い（干渉の性質を示さない、または放出される電場の位相がラ

ンダムである。また、放出される電磁波の向きは、ある定まった方向があるわけではなく、全ての方向に電磁

波は放出される。よって、電磁波が原子により構成された媒質を通過後、電磁波のエネルギーが減少する。こ

れが媒質による電磁波の吸収である。なお、電場の可干渉性、コヒーレンスについては、簡単な記述が p.10

にある。

媒質中において光電場の速度が減少する理由

もし 0ωω ≠ なら、原子は電磁波を吸収せず、原子を含む媒質はその電磁波に対して透明である。この状況

では、電磁波は原子を非共鳴的に振動させ、その振動周波数は電磁波の周波数ω である。原子の振動は電磁波

の振動に追従するが、その振動には位相の遅れが生じる。この位相の遅れは摩擦による減衰により生じる（減

衰のある調和振動子の強制振動で、位相の遅れのあることが、大抵の本で説明されている。）。振動している原

子は全て電磁波を放射するが、電場位相の遅れは、原子を含む媒質にわたって起こるため、電磁波の波面に時

間遅延が生じる。これが、媒質中の電磁波の速度が真空中の速度よりも遅い理由である。速度の減少は、この

ように起こる。そして、媒質の屈折率は、真空中の光電場の速度に対する媒質中の光電場の速度の比により定

義される。

自由電子の集団運動

電子と原子核による双極子、イオン同士による双極子は、どちらも束縛された振動子である。ところが、金

属やキャリアをドープした半導体では、自由電子の数が非常に多い。このような自由電子は、個別の振る舞い

ではなく、集団としての運動を行う。

E

- - ++++

51

平衡状態では、自由電子が空間的に一様に分布しているとしよう。しかし何らかの原因でその分布に乱れが

生じると、ある場所での電荷密度分布が高くなる。しかし電子が集合すると、同じ符号の電荷を持つ電子にク

ーロン力が働き、集合した電子は、もとの平衡点に戻ろうとする。このようにして、自由電子の集団としての

動きが生じる。平衡点に自由電子が戻っても、いきなり止まる事はできないため、平衡点を通過し、その後再

び平衡点に戻ろうとする。自由電子の集合は、このような振動運動を行う。このような振動を、プラズマ振動

という。

自由電子が支配的な金属では、プラズマ周波数（後ほど説明する）よりも高い周波数の電磁波に対して透明

になり、プラズマ周波数よりも低い周波数の電磁波を全反射で反射する。金属光沢は、この考えに基づいた模

型で説明される。

52

電子振動子模型 --- 誘電体・絶縁体に対する模型 ---

この模型は、量子力学の発展以前に原子物理学者たちによって用いられてきた。電磁波に最もよく応答するの

は，物質内の電子である。この模型では、1個の原子の中で、調和振動子のポテンシャルの谷の中に捉えられた

電子により電気双極子が誘起される。電子振動子の共鳴周波数と吸収線の線幅は、実際の遷移にあうように選ぶ。

また、振動子と電磁場の相互作用の強さを表すパラメーターである振動子強度（oscillator strength）を導入す

る。これを用いて計算した吸収強度が実験値と合うように定める。

1次元における質量m 、電荷 e− の電子振動子の運動方程式は、

)()()()(

2

2

tEm

etx

m

k

dt

tdx

dt

txd−=++ γ

ここで )(tx は電子の平衡位置からの変位であり、γ は減衰定数、kx は復元力、 )(tE は電場の瞬時値を表す。

なお、この模型では、電子は原子核に束縛されており、巨視的な電荷の流れ、すなわち電流はない。

電場 )(tE と原子の変位 )(tx を次のように表す。（電場や変位は実数であるため、計算は )exp( tiω− のまま

行ない、最後に実部をとればよい。）

)exp()()(),exp()( 0 tiXtxtiEtE ωωω −=−=

固有（角）周波数 mk /0 =ω とおくと、上式は

ωγωωωγωω

i

EmeXEmeXiX

−−−=→−=−−

22

0

00

22

0

)/()/()(

となり、変位 )(tx は、次のようになる。

)exp()/(

)( 022

0

tiEi

metx ω

ωγωω−⋅

−−−=

われわれは電磁波の周波数 0~ ωω の近傍に興味があるので、次の変形をする。

γωωωω

γωωωωωω

ωγωωωωω

i

tiEme

i

tiEme

i

tiEmetx

−−−

−=−+−

−−≈

−+−−

−=)(2

)exp()/(

))((

)exp()/(

))((

)exp()/()(

0

00

0000

0

00

0

分極（巨視的な双極子モーメント） )(tP は、単位体積中にN個の振動子が存在する場合、1個の振動子（電

子）の双極子モーメント )()()( txetp −= のN 倍で与えられ、

)exp(2/)(

)2/(

)(2

)exp()/()()()( 0

0

0

2

0

00

2

tiEi

mNe

i

tiEmNetxeNtP ω

γωωω

γωωωω

−+−

−=−−

−=−=

と表される。電子の電気感受率 )(ωχ は、誘起された分極 )()()( 0 tEtP ωχε= の式で定義され、一般には複

素数となる。 )(")(')( ωχωχωχ i+= とおくと、

22

000

2

22

0

0

00

2

000

2

)2/()(

2/

2)(",

)2/()(2)('

2/

1

2)(

γωωγ

εωωχ

γωωωω

εωωχ

γωωεωωχ

+−⋅=

+−−

⋅−=

+−⋅−=

m

Ne

m

Ne

im

Ne

　

となる。 )(" ωχ は半値全幅（full width at half maximum:FWHM） 2/γ のローレンツ形（ )1/(1 2 +x のよう

な関数形）であり、 )(' ωχ はそれを微分したような形状である。半値全幅は、半値幅とも言われる。ここでは、

)(" ωχ がとる最大値 Max"χ の半分の値 2/"Maxχ を与える2つの周波数の差を、半値幅という。

53

)(ωχ が複素数であるため、物質の誘電率は複素数になるので、 'ε と書く。電束密度D の定義：

EEEPEED ))(1()(' 0000 ωχεωχεεεε +=+=+==

から、誘電率 'ε は、

))(1(' 0 ωχεε +=

である。よって、 1|)(| <<ωχ の場合、 κin + は次のように簡単に求まる。（p.42の（1）式参照）

)("2

1),('

2

11

)("2

)]('2

11[)(

2

11~)(1

))(1(',

'

0

0

00

ωχκωχ

ωχωχωχωχε

ωχεεε

κεε

=+=→

++=++=+

=+=

n

iin

　　　

一般には、

( ) ( )

( ) ( )

+++−−=

++++=

→++=+

2/122

2/122

0

0

)(")('1)('12

1

)(")('1)('12

1

)(")]('1[))(1(

ωχωχωχκ

ωχωχωχωχωχ

εωχε

n

i

で与えられる。

複数の周波数を持つ振動子が単位体積あたりN 個あり、そのうち周波数 jω を持つ振動子の割合がjf なら、

1,2/

1

22/

1

2)(

0

2

0

2

=−−

⋅=+−

⋅−= ∑∑∑j

j

j jjj

j

j jjj

jf

im

fNe

im

fNe　　

γωωεωγωωεωωχ

となる。ここでjf を振動子強度という。

図は 4.0/5/0 =γω の場合

00

2

0max

2

0max

2

0

0

max

)(""

12/

1

"

)(",

12/

2/

"

)('

γεωωχχ

γωωχ

ωχ

γωω

γωω

χωχ

m

Ne==

+

−=

+

−

−

−= 　　

12/

2/

2/

1

22

2

1

2)(

2

0

0

00

2

000

2

000

2

+

−

−−

⋅−=

+−

⋅⋅

−=+−

⋅−=

γωω

γωω

γεω

γωωγ

εωγ

ωωεωωχ

i

m

Ne

im

Ne

im

Ne

　　　

1|)(| <<ωχ の場合、複素屈折率 κinn +=~ は上で計算したように、きれいな式になる( 2/)(1~ ωχ+=n )

ので、屈折率nの電磁波の周波数ω 依存性に関する以下の議論は、周波数の範囲を含め厳密に正しい。

さて、 )(' ωχ は Lωγωω ≡−= 2/0 に最大値、 Rωγωω ≡+= 2/0 に最小値を持つ分散型( ]1/[ 2 +xx の形)

-1

-0.5

0

0.5

1

1.5

-4 -2 0 2 4

χ"(ω)/χ"max

χ'(ω)/χmax

(ω−ω0)/γ

54

の曲線である（このページの補足参照）。一方 )(" ωχ は、 0ωω = にピークを持つローレンツ型( ]1/[1 2 +x の

形)の曲線で、その半値全幅（ )(" ωχ が最大値の1/2の値をとる時の2つのω の差： LR ωω − ）はγ である。

Lωω << および Rωω >> の領域では、物質は電場に対してほぼ透明である（ 0~0~)(" κωχ → 。注意：

真空中と物質との境界で反射がゼロという意味ではない。吸収がない事をさして、透明と言っている。）。また

Lωω < および ωω <R の領域では、屈折率は、周波数ω が大きくなるに従い大きくなる。一方 0ω 近傍

（ RL ωωω << ）では、周波数ω が大きくなるに従い屈折率は小さくなる。屈折率が周波数に依存する現象

を分散とよび、 0/ >ωddn の場合（ Lωω < および ωω <R の領域）を正常分散、 0/ <ωddn の場合

（ RL ωωω << の領域）を異常分散とよぶ。

1|)(| >ωχ の場合は、上の議論は定性的に正しい。周波数の範囲が少し変わる。

複数の電子振動子がある場合（図では

3つ）、右図のような振る舞いをする。

三角プリズムによる光の分散現象：波長ごとに屈折率が異なる

ため、異なる波長成分が空間的に分離される。

補足：

12/

2/

"

)('2

0

0

max+

−

−

−=

γωω

γωω

χωχ の振る舞い

0

00

11

0|)(

|)('

|

)1(

1

)1(

2)1()1()2()1(1)('

)1(1

)(,2/

22

2

22

222212

12

2

0

増加

　　

極小値減少極大値

　　

増加

　　　　　　

+

∞+

−

+−

+

∞−

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

+−

=+

++−=+⋅−⋅−+⋅−=

+⋅−=+

−=−

≡

−−

−

xf

xf

x

x

x

x

xxxxxxxf

xxx

xxfx

γωω

よって )(' ωχ は、 12/

0 ±=−

γωω で極値を取る。

赤緑紫青

白色光

ω

n

κ

ω

55

ドルーデ（Drude）模型

重い金属イオンと軽い電子との混じった電気的に中性の集合体を、プラズマと言う。金属や、不純物をドー

プした半導体は、同じ数の位置が固定されたイオンと自由電子からなるので、プラズマとしての扱いが出来る。

自由電子は、外部から力が加わっても、復元力（元の位置に戻ろうとする力：ばねによる力など）が働かない。 1次元の自由電子の運動方程式は、

)()()(

2

2

teEdt

tdxm

dt

txdm −=+

τ

ここで )(tx は電子の位置であり、τ は電子同士の衝突までの時間、または電子がイオンに衝突するまでの時間

の平均値、 )(tE は電場の瞬時値、 )( e− は電子の電荷を表す。電場 )(tE と原子の変位 )(tx を次のように表す。

（電場や変位は実数であるため、計算は )exp( tiω− のまま行ない、最後に実部をとればよい。）

)exp()()(),exp()( 0 tiXtxtiEtE ωωω −=−=

すると、運動方程式から

02

0

2 1)/( E

im

eXEmeX

iX ⋅

+⋅=→−=−−

τω

ωτω

ω


)exp(1

)( 02

tiEim

etx ω

τω

ω−⋅

+⋅=

1個の電子の双極子モーメント )(tp は )()()( txetp −= であるから、単位体積中にN個の振動子が存在する

場合、巨視的な双極子モーメントすなわち分極 )(tP は、 )(tp のN倍となり、

)exp(1

)exp(1

)()()()( 02

2

02

tiEim

NetiE

im

eeNtxeNtP ω

τω

ωω

τω

ω−⋅

+⋅−=−⋅

+⋅⋅−=−=

と表される。分極 )(tP と電場 )(tE との関係式： )()()( 0 tEtP ωχε= から電気感受率 )(ωχ が求まる。

m

Ne

iim

Ne

im

NeP

P

0

22

2

2

20

2

2

2

0 ,1

)(1

)(ε

ω

τω

ω

ω

τω

ωεωχ

τω

ωωχε =

+−=

+⋅−=→

+⋅−= 　。

Pω をプラズマ周波数という。実部と虚部に分けて書くと

22

2

22

2

22

2

2

2

2

2

2

2

2

20

2

/1

)/(

/1

1/1

)(

11

1

1

11)(

τωωτω

τωω

ωττωω

ωτ

ωτωω

ωτωω

τω

ω

ω

τω

ωεωχ

++

+−=

−⋅+

−=+

−⋅−=

+⋅−=

+−=

+⋅−=

PP

PPPP

i

i

i

iiim

Ne

　　　

となり、 )(ωχ は一般には複素数となる。

56

ドルーデ模型の光学的性質（緩和定数 ∞=τ の場合の電気感受率 )(ωχ ）

まず、運動がほとんど抑制されない（ ∞=τ ）電子の系を考える。すると、

2

2

2

0

2

)(/

1)(

ωω

ωχτωωε

ωχ P

im

Ne−=→

+⋅−=

となる。誘電率 ))(1(0 ωχεε += であり、誘電率ε と κ,n の関係式は次のとおり。（p.42の（1）式）（注意：

誘電率が実数でも、屈折率は複素数になる。）

2

2

0

22

0

00

1)/1('

,~

ωω

εωωε

εε

κεε PPinn −=

−=+==

空気と物体との境界面での反射率R は、垂直入射の場合、

22

222

)1(

)1(

1~1~

κκ

+++−

=+−

=n

n

n

nR

で与えられる（既に示した。p.45）。

物体の厚さが十分厚いとしよう。 Pωω < の時 n~/ 0 =εε は純虚数である。よって、電磁波の反射率 1=R

となる。 Pωω > では反射率が下がり、 ∞→ω では 1=n となり、反射率 0→R となる。金属には、自由電

子がたくさんあるため、 Pω は十分大きく、 Pωω < であるような可視光の領域で 1=R （全反射）である。こ

れが、金属が全反射して光沢を持つ理由である。

ドルーデ模型での自由電子の緩和の効果（伝導率、joule熱への影響）

電子の運動量 dttmdxtp /)()( = として、電子の運動方程式を書き換える。

)()()(1)(

)()()(

2

2

tEetpdt

tdpteE

dt

tdxm

dt

txdm −=+→−=+

ττ

この運動方程式は、電子が電場 )(tE により加速されるが、その運動量を時間τ 程度で失う事を示す。よってτは、運動量散乱時間（momentum scattering time）といわれる。電子が電場により加速されても、イオンとの衝

突、電子同士の衝突により一定速度以上には速くならない。これは、衝突によりエネルギーが失われる事を意

味し、そのエネルギーの損失は、ジュール（joule）熱となって表れる。

AC（交流）電場 )exp()( 0 tiEtE ω−= を電子に加える。解として、 )exp()( 0 tiptp ω−= を仮定する。

)exp(/1

)()()()()(

1)()()()(

1)( 0 tii

EetptEetptpitEetp

dt

tdpω

τωτω

τ−

+−−

=→−=+−→−=+

電子の速度 mtptv /)()( = は、

)exp(1

1)()exp(

/1

/)(/)()( 00 tiE

im

etiE

i

memtptv ω

ωττ

ωτω

−−

⋅−

=−+−

−==

で与えられるので、自由電子の密度（単位体積あたりの電子の個数）をN として、電流 )(tJ および伝導度 )(ωσが以下のように求まる。（オームの法則が交流でも成立するとしよう： )()()( tEtJ ωσ= から )(ωσ が求まる。）

)()()exp(1

1)()()( 0

2

tEtiEim

NetveNtJ ωσω

ωττ

=−−

⋅=−=

反

射

率

1 ω/ωP

1

0

57

m

Ne

iim

NeAC

τσ

ωτσ

ωττ

ωσ2

00

2

,11

1)(: =

−=

−⋅= 　伝導度。

ここで mNe /2

0 τσ = は、DC（直流）電場で測定した伝導度である。金属、キャリア（電子または正孔）を

ドープした半導体の運動量散乱時間 )s(10~10 1314 −−=τ であり、光(可視光程度)の周波数に相当する。よって

DC伝導からのずれを観測するためには、光の周波数に相当する電磁波、つまり光を用いる必要がある。

さて、この場合のジュール熱 )()( tEtJ を計算しよう。 )exp()(),( titbta ω−∝ の時、 *]Re[2

1)()( BAtbta ⋅=

の関係（p.28）から、時間平均したジュール熱 )()( tEtJ を求めると、

2

02

22

02

22

0

2

2

0

2

00

2

)(1

1

2

1

)(1

1Re

2

1

1

1Re

2

1)()(

1

1)exp()exp(

1

1)(*)(

Em

NeE

i

m

NeE

im

NetEtJ

Eim

NetiEtiE

im

NetEtJ

ωττ

ωτωττ

ωττ

ωττ

ωωωτ

τ

+⋅=

++

⋅=

−⋅=→

−⋅=+×−

−⋅=

　　

となる。

以下では、2つの極端な場合について考えよう。

(i) 1<<ωτ ： ωτ /1<< より、交流電場の振動よりも頻繁に電子の散乱が生じるので、電場の向きが変わら

ない間に電子散乱（電子とイオンや、電子同士の衝突）が起こると見なせる（ ω/1~ よりも短い時間の間では、

電場の向きは、ほとんど変化しない。）。よって、電子が散乱されるまでの時間では、電場の大きさ・向きが変

化しない直流電場と見なせる。その場合、電場による力 )(teE と電子の速度は同位相で振動するため、ジュー

ル熱 0)()( ≠tEtJ である。これは、電場から電子にエネルギーが与えられる事を意味する。

(ii) 1>>ωτ ： τω <</1 より、電場が振動する間に電子がほとんど散乱されず、電子が孤立して真空中にい

るかのように見なせる。この時、電場 )exp()( 0 tiEtE ω−= に対して電子の速度 )(tv は、

)2/exp()(

)exp()(

~)exp(1

1)()( 000 πω

ωττ

ωωτ

τω

ωττ

itiEm

etiE

m

eitiE

im

etv +−⋅

−=−⋅

−⋅−

−⋅

−=

となる。電子の速度変化と電場の振動は位相が90度ずれるので、電場が電子に与えるエネルギーは0になる。

（注意：速度は )(tv ではなく、 )](Re[ tv として、理解する事。）

0)(1

1

2

1)()( 2

02

2

→+

⋅= Em

NetEtJ

ωττ

よって、電子と電磁波とのエネルギーのやり取りはない。

光学的手法による )(ωχ とAC伝導度 )(ωσ の同等性

さて、1次元の自由電子の運動から、巨視的な電気双極子（分極） )(tP と電流 )(tJ を導いた。

)exp(1

1)()()(),exp(

1)()( 0

2

02

2

0 tiEim

NetEtJtiE

im

NeEtP ω

ωττ

ωσω

τω

ωωχε −

−⋅==−⋅

+⋅−==

)(ωχ を式変形すれば、以下のように、 )(ωσ と関連付けられる。

58

ωωσ

ωτωσ

ωτωτ

ωτωτ

ωτωτ

τω

ωωχε

)(

1

1

1

1/

1

111)( 0

222

2

2

0

i

ii

i

mNei

im

Nei

im

Ne

im

Ne

=

−⋅=

−⋅=

−⋅=

+⋅−=

+⋅−=

　　　　

となり、ドルーデ模型によれば、電気感受率 )(ωχ とAC伝導度 )(ωσ とは、お互いに関連することが分かる。

すなわち、光学的手法による )(ωχ の測定は、AC伝導度での )(ωσ 測定と同等である。

プラズマ振動

ドルーデ模型においてほとんど運動が抑制されない（ ∞=τ ）時、電子系の電気感受率 )(ωχ は、

2

2

22

2

22

2

)(/1

)/(

/1)(

ωω

ωχτω

ωτωτω

ωωχ PPP i −=→

++

+−=

であった。また、

2

2

0

2

2

0

00

1

)1(

,~

ωω

εωω

ε

εε

κεε P

P

inn −=−

=+==

から、 Pωω ~ で、 1~)( −ωχ であり、 0~~ κn を得る。その時電束密度 EED ))(1(0 ωχεε +== で与

えられるため、 0=D である。電場E はイオンと自由電子で作られ、分極P の向きは電場E の向きと逆向き

である。このように、 Pωω ~ の振動は、誘電率 0=ε での振動に対応する。

ここで、自由電子群に働く電場の力から、自由電子群が振動する事を定

性的に理解しよう。外部電場のない金属を考える。固定された格子の位置

にある正のイオンと、自由電子から構成されており、その電荷は全体で中

性である。簡単のため、イオンと電子からなる平板を考える。平衡状態で

は、自由電子が空間的に一様に分布しているとしよう。

図に示しように、自由電子全体の位置を右方向にu だけずらす。すると、正のイオン位置は変わらず、電子

の位置が変化するので、分極（電気双極子）ができる。分極は、電場を作る。電場の向きは、自由電子の移動

を妨げる向き（図で右向き）であり、電場による力は、負の電荷を持つ自由電子群を元の位置（左方向）に戻

そうとする。この電場の力は、電子が平衡点からずれた場合、常に元の位置に戻すような復元力にとして働く。

そのため、電子の集団は、平衡点に戻ろうとする。このようにして、

自由電子の集団としての動きが生じる。平衡点に自由電子群が戻って

も、いきなり止まる事はできない（動きを急に止める力が無い）ため、

平衡点を通過し、その後再び平衡点に戻ろうとする。この繰り返しに

より、自由電子群は、固定されたイオンの格子の周りで振動する。こ

の振動は、プラズマ振動と呼ばれている。

プラズマ振動の周波数は、以下のようにして求まる。平衡状態では、正のイオンと電子群は、ずれることな

く空間的に重なっている。しかし、電子群の位置の移動u により、分極P が生じる。電子の位置がu ずれると、

そのずれの向きと反対方向に1つの電子あたり、 ueuexxep )()0()( −=−=−= −+ の電気双極子が出来る。

ここで、 +x は正イオンの位置（位置の変化なしで0）であり、 −x は電子の位置である。すると、N を単位体

+－－－－ +－－－－ +－－－－ +－－－－ +－－－－

－－－－ +－－－－ +－－－－ +－－－－ +－－－－+

+－－－－ +－－－－ +－－－－ +－－－－ +－－－－ +

+

+

+

電場Eの向き

----

----

----

u

平衡状態

－－－－ +－－－－ +－－－－ +－－－－ +－－－－ +－－－－

+－－－－ +－－－－ +－－－－ +－－－－ +－－－－+

－－－－ +－－－－ +－－－－ +－－－－ +－－－－ +－－－－

59

積あたりの電子の数として、巨視的な電気双極子すなわち分極P は ueNNpP )(−== で与えられる。一方、

分極によって生じる電場（外部電場ではない。自由電子によって生じる電場の事。）E は、 0/εNeuE = であ

る（このページの補足参照）。この時、電束密度 00 =+== PEED εε である。

電場により自由電子群が受ける力の向きは、自由電子が元に戻るような向きである。イオンは十分重く動か

ないとする。すると、微小面積dS 、微小長さdL （体積dSdL ）の自由電子群の運動方程式は、電子の質量

をm として、

um

Ne

dt

udNeuedSdLNEedSdLN

dt

udmdSdLN

0

2

2

2

0

2

2

))(())(()(εε

−=→⋅−=⋅−=

となる。これは、単振動の運動方程式と同等であり、角周波数は )/( 0

2 mNe ε で与え

られる。この値は 1 次元の自由電子のドルーデ模型で得られたプラズマ周波数 Pω に一致する。すなわち、電

子群は、プラズマ周波数 Pω で単振動する。

まとめると、分極P は ueNP )(−= で与えられる（分極の向きは、負の電子から正のイオンヘ向かう方向）。

一方、正のイオンと電子によって作られる電場E の向きは、正のイオンから負の電子へ向かう方向である。よ

って、分極ベクトルの向きは電場の向きとは逆向きである。そのため、電束密度 00 =+== PEED εε とな

り、誘電率 0=ε となる。逆に言えば、プラズマ振動は、誘電率 0=ε の時の自由電子群の振動である。

補足：電子群のずれによって生じる電場の大きさが 0/εNeu で与えられる事

平行に並んだ 2枚の導体の板にそれぞれ単位面積当たり QQ −+ , の電荷を帯電させる。そのうち、+に帯電

した導体の表面にまたがる底面積dS の微小な円筒を考え、この円筒状でGaussの法則を適用する。対称性か

ら電場の成分は導体板の面に垂直な成分しかない事が分かる。誘電率 0ε の物体の中を電子群が移動するとして、

その電場を +E とすれば、

00 2εεQ

EQdS

dSEdSE =→+

=+ +++

となる。 +E の電場がプラスの符号であるのは、右図が示す

ように、電場の向きが導体から外向きである事を示す。

Q− で電荷で帯電した導体板についても同様の計算をす

ると、電場 −E は、

00 2

)(

εεQ

EdSQ

dSEdSE −=→−

=+ −−− 。

ここで、2 つの電場 +E と −E は、向きが反対方向で、その大きさは同じである事に注意しよう。 −E の電場

にマイナスの符号がついているのは、右上の図が示すように、電場の向きが導体に向かっての向きである事を

示す。すると、2つの導体板にはさまれた領域の電場の大きさをE とすれば、電場の向きを考慮し、電場の合

成（先ほどの図を参照）により、

0εQ

EEE =−= −+

となる。この結果を、プラズマ振動の問題に適用する。電子と正のイオンとのずれの大きさをu とすると、ず

れによって生じる単位面積あたりの電荷は、 NeuQ = である。よって、電子群のずれによって生じる電場の

大きさは、 0/εNeu で与えられる。

dL

dS

dS

E+ E+ E－ E－

+Q －Q

60

ドルーデ模型の光学的性質： κ,n の周波数依存性の詳細

電束密度： EEPEED )(000 ωχεεεε +=+== から、複素誘電率 ))(1(' 0 ωχεε += であり、既に求め

た関係式（p.55）から、 IR iεεε +=' と書けば、 IR εε , は次式で与えられる。

22

2

0

22

2

00

22

22

2

22

2

2

2

0

/1

)/(,

/11,

/1

)/(

/1)(

))(1('

τωωτω

εε

τωω

εε

εω

τωωτω

τωω

τω

ω

ωωχ

ωχεεεε

+=

+−=→=

++

+−=

+−=

+=+=

PIPRP

PPP

IR

m

Nei

i

i

となる。また、

=

−=→+=+→+=

κεε

κεεκεεεεκ

εε

n

niniin

I

R

IR2/

/)(//

'

0

22

02

00

0

の関係式が成り立つ。さらに、以下の式を既に導出した（p.44）。

2/12

0

2

00

2/12

0

2

00 2

1,

2

1

+

+−=

+

+=

εε

εε

εε

κεε

εε

εε IRRIRRn 　

例えば銅の場合、 )s(101.4),s(107.1 1131116 −−− ×=×= τωPなので、以下では

1−>> τωP として考察する。

（I） 1−>>>> τωω P （ 1,1 <<>>ωω

ωτ P ）の場合

電子が散乱される時間τ の間に、電磁波は多数回の振動をする。これは、電子があまり散乱されない（衝突

するまでに、電子は多数回の振動をする）事を意味するので、伝導電子による電磁波の吸収は小さい。

0~/,1~/1~)(1/'

0~~)/1(

~)]/(1[/

1)(

000

2

2

2

2

2

2

2

0

2

εεεεωχεε

ωω

ωω

ωτωω

τωωεωχ

IR

PPP

iiim

Ne

→+=→

−∞+

−+

−=+

⋅−=

　　

であるから、 0~,1~ κn となる。（ κ,n の具体的な表式からも、 κ,n の大きさの程度が求まる。）よって反射

率 0~R （p.45）となり、金属は電磁波に対し完全に透明となり、吸収もない。この領域は、紫外領域である。

（II） )1(1 >>=>>= − ωττωτωω PP の場合

11

/

/1

)/(/,

1

1

11

/11/

2222

2

22

2

02222

22

22

2

0 +=

+=

+=

+=

+−=

+−=

τωτω

τωωτω

τωωτω

εετωτω

τωτω

ωεε

P

P

P

PPPI

PP

PPR

近似計算をすると、 0~~ κn を得る。この時、プラズマ振動をする。

κ,n の具体的な表式からの計算（近似計算の詳細）

0~11

2

1~

1

1

1

1

2

11

1

1

1

1

2

1

11

1

1

1

2

1

2

1

2/1

22

2/1

2222

2/1

222

2222

2/12

22

2

2222

2/12

0

2

00

+

++

+=

+

++

+=

++

++

+=

+

+=

τωτωτωτωτω

τωτω

τωτω

τωτωεε

εε

εε

PPPP

P

PP

P

P

PP

IRRn

　

61

0~1

1

1

1

2

11

1

1

1

1

2

1

11

1

1

1

2

1

2

1

2/1

22222

2/1

222

2222

2

22

2

2222

2/12

0

2

00

++

+−=

+

++

+−=

++

++

+−=

+

+−=

τωτωτω

τωτω

τωτω

τωτωεε

εε

εε

κ

PP

P

PP

P

P

PP

IRR

　　

（III） 1−>>>> τωωP （ 1,1,1 >>>>>> τωωω

ωτ PP ）の場合

−

→

∞+⋅

+=

+=

−−=∞+

−+

−=+

−=

2

2

0

2

2

0

2

2

2

2

22

2

22

2

0

2

2

2

2

2

2

22

2

22

2

0

1~/

~/

1~

)/11(

1~

])/(11[

)/(

/1

)/(/

~1)/11(

1~])/(11[

1/1

1/

ωω

ωτεε

ωω

εε

ωω

ωτωω

ωτωτωωτω

τωωτω

εε

ωω

ωω

ωω

ωτωω

τωω

εε

PI

PR

PPPPI

PPPPPR

これより、多少の近似計算の後、ωω

κ Pn ~,0~ を得る。

計算の詳細

02

1~01

2

1~

11

2

1

1

2

1~

2

1

2/1

2

2

2

22/1

2

2

2

2

22/1

2

2

2

2

2

2

2/12

2

22

2

2

2

2

2/12

0

2

00

=

+−

++−

++−=

+

−+−

+

+=

ωω

ωω

ωω

ωω

ωτωω

ωω

ωω

ωτωω

ωω

εε

εε

εε

PPPPPP

PPPIRRn

　　

ωω

ωω

ωω

ωω

ωτωω

ωω

ωω

ωτωω

ωω

εε

εε

εε

κ

PPPPPP

PPPIRR

~22

1~01

2

1~

11

2

1

1

2

1~

2

1

2/1

2

22/1

2

2

2

2

22/1

2

2

2

2

2

2

2/12

2

22

2

2

2

2

2/12

0

2

00

++

++=

+

−+

−−

+

+−=

　　

ωωκ /~ Pかつ ωω >>P から、 1>>κ となる。よって、 1>>κ から電磁波の伝播は起こらず減衰する事が

わかる（p.43）。また、周波数が小さくなるにつれ、減衰が大きくなる。伝播が無く減衰が起こっている事、お

よび 0~,1 n>>κ から、電磁波が反射される割合が大きい事（反射率 1~R に近い状態）がわかる。

（IV） ωτω >>>> −1

P （ 1,1,1 >>>><< τωωω

ωτ PP ）の場合

電磁場が 1 回振動する時間（1 回の振動時間： ωπ /2~ ）よりも、電子の散乱時間τ のほうが小さい

)/1 τω >> 。すなわち、電子が多数回散乱されても、外部から加わる電場（外部電場）の大きさと向きはほと

んど変わらないため、電子から見ると、外部電場は直流電場のように見える。よって、外部電場により分極し

た電子が作る電場により、外部電場が遮断されるので、電磁波は導体の中深くには進入できない。以下の計算

で、その事を確認しよう。

電磁波の周波数ω が1−<< τω である時、次のような近似を用いる。

62

1~1

/

/1

)/(/2

~1~1

1/1

1/

2

22

2

22

2

0

2222

22

22

22

2

0

22

>>⋅=+

=+

==

−−+

−=+

−==−

τωωω

ωτω

τωωτω

τωωτω

εεκ

τωτωτω

τωτω

ωεεκ

PPPPP

I

PPPP

R

n

n

、

すると、 κ,n に対して、次のような値を得る。

2/12

2~~

ωτω

κ Pn

近似計算の詳細

( ) ( )

( ) ( )

2/12

2/122

2/122

22

2/12

22222

2/12

2222222

2/12

222222

2/12

0

2

00

2/12

2/122

2/122

22

2/12

22222

2/12

2222222

2/12

222222

2/12

0

2

00

22

1~

2

1~

11

2

1~

2

1~

2

1~

2

1

22

1~

2

1~

11

2

1~

2

1~

2

1~

2

1

=

++

+++

+++

+−++

+

+−=

=

+−

++−

++−

+−+−

+

+=

ωτω

ωττω

ωττω

τωωτ

τωτω

ωττω

τωτωω

τωτωτω

εε

εε

εε

κ

ωτω

ωττω

ωττω

τωωτ

τωτω

ωττω

τωτωω

τωτωτω

εε

εε

εε

PPPPPP

PPP

PPP

IRR

PPPPPP

PPP

PPP

IRRn

　　

　　　

表皮深さ（skin depth）（金属中を電磁波が進入する深さ）

さて、物質中の電場ベクトルは、

( )[ ]

−−−=+−−= z

cz

c

ntizikktitz κ

ωωω

00

00 )(exp)"'(exp),( EEE

であり、その時、光強度 )(zI は、2|),(|)( tzzI E∝ から、

)exp()0()( zIzI α−=

で与えられ、0/2 cωκα = は吸収係数であった（p.43）。

ωτ >>−1であるような低周波数領域での電磁波の伝播を考える。プラズマ周波数 Pω 、真空での光の速さ 0c 、

およびDC電気伝導度 0σ は、

00

2

0

2

0

0

0

2

00

22 1

,,1

µετ

στε

σττεε

ω ===== cm

Ne

m

Ne

m

NeP

で与えられた。（IV）で計算した結果を利用すると、吸収係数α は、次式で与えられる。

( ) 2/1

00

2/1

0

0

2

00

0

2/1

2

00

0

2/1

0

0

2

0

2/1

2

0

22/1

2

00

22222

2

22ωµσ

µµ

εωσ

εωσ

τεσωττωω

ωτωωωκ

α =

⋅=

=

⋅=

=

==

cccccc

PP 。

よって、吸収係数α は「DC伝導率 0σ と周波数ω の積」の平方根に比例する。この結果は、AC（交流）電場

は伝導体の非常に短い距離しか進入できない事（2/1ωα ∝ ）を示す。この現象は、表皮効果（skin effect）と

いわれる。また、

63

00

22

ωµσαδ ==

を表皮深さ（skin depth）という。表皮深さは、伝播する電場の振幅が e/1 の大きさになるときの電場の進入長

として定義されている。よって、吸収係数α の半分の値になる。

銅の場合を見よう。なお、真空での透磁率 m/A)(T(H/m),104 7

0 ⋅×= −πµ 。銅の室温での DC 伝導率は、

m)/1(105.67

0 Ω×=σ である。周波数 Hz50=f では、 4rad/s31=ω になる。この時、 mm8.8=δ である。

また MHz100=f は、 rad/s106.28 8×=ω であり、その時、 m6.2µδ = である。

さて、周波数 Pω は、金属では可視領域にある。よって Pωω < であるすべての可視領域の光は反射される事

になり、金属は金属特有の光沢面を示す。ただし、金属の厚さl が δ≤l の場合は、電磁波は、その振幅をあ

まり小さくする事なく、金属を伝播し透過する。

64

局所電場補正：ローレンツ（Lorentz）の局所電場モデル

物質内の点rにおける巨視的な電場（原子の大きさよりは十分大きな、ある大きさで平均を取った、物質内

の電場の大きさ）Eが、一般にMaxwell方程式で用いられている電場の量である。

さて、外部から物質に電場0E を加えると、それに応じる形で、原子内の原子核と電子の位置の平衡点から

のずれや、正イオンと負イオンとの位置の平衡点からのずれなどが生じる。このような正負のイオンや原子核・

電子の位置のずれを電気双極子といい、 pで表す。電気双極子の単位体積あたりの大きさを、分極P で表す。

正負のイオン、原子などがその平衡点からずれると、それらは、電場を作る。分極P により作られる電場の大

きさは、物質に依存する。分極P により作られる電場の存在により、物質内の巨視的な電場Eと外部から物質

に加えられた電場0E との違いを生じる。

では、物質を構成する個々の原子に加わる電場の大きさは、巨視的な電場Eおよび分極P と、どのような関

係があるのだろうか。物質内における巨視的な電場Eは、平均操作を施された量であるため、物質中の1個の

原子に加わる電場の大きさと異なってもよい。原子を1個選び出したとき、その原子の位置での電場を、局所

電場localE という。

気体のような希薄な物質の場合は、巨視的な電場Eと原子の位置での局所電場localE との違いは無視できる

かもしれない。しかし、液体や気体の場合には、それら物質中に数多くの原子が存在する（着目する原子の近

くにも他の原子がたくさん存在する）ため、巨視的電場Eと局所電場localE の違いが大きくなる。このように

して、物質中のある原子に加わる電場localE は、巨視的な電場Eそのものではない。

では、どのようにして局所電場localE を求めるのだろうか？まず、物質中のある原子に着目する。外部電場E

が物質に加えられている際に、ほかに影響を与えないように（どのようにするのだろうか？）しながら、着目

する原子を囲む微小な部分をくりぬいて真空にする。そのくりぬかれた位置に現れる電場が、局所的な電場

localE である。物質を連続体として扱い、微小な球をくりぬいたとすると、その球の表面には表面電荷σ が現

れる。物質の分極をP とすれば、球の中心での電場の大きさlocalE は、

03εP

EE +=local

で与えられる。（この関係式は、すぐ後で説明する。）この局所電場 localE を Lorentz（ローレンツ）の局所電

場という。実際の物質は、こんな単純な構造をしていないので、連続体での近似や球をくりぬく事は、多少な

りとも正確さに欠けると考えたほうがよい。しかし、現象の粗い把握には、便利である。

localE の求め方

既に述べたように、物質内における巨視的な電場E（「外部から加えられた電場」と「物質の分極によって

出来る、平均的な電場」の、2 つの電場を足し合わせて出来る電場）は、平均操作を施された電場である。し

かし、物質内のある点 xに着目すると、点 xの近くにある双極子がとる微視的な空間的構造に起因した電場に

よる寄与もある。そこで、局所電場localE を求めるにあたり、次のような手続きで計算する。

物質に外部から電場が加われば、それに伴い、物質中に分極P が出来る。分極P により生じる電場を計算す

るためには、電気双極子 pにより出来る電場を計算し、物質の内部の電荷密度 )(xeρ と物質の表面にできる表

面電荷密度 )(xσ により出来る電場を求める必要がある。さらに、物質内部の電気双極子 pの微視的な空間的

配置に起因した電場も考慮する必要がある。そのため、物質中の着目する点 xを中心とする径aの球の内部と

65

球の表面、それ以外の物質中の内部全て、最後に物質の外表面の 4 つの箇所に分けて電場の寄与を考えよう。

半径aの球を考えるのは、計算を簡単にするためである。

このように物質を4つの部分に分けて考えると、分極P に起因する

物体内の電荷密度 )(xeρ 、および表面電荷密度 )(xσ を求める必要が

ある事がわかる。 )(xeρ と )(xσ は分極P から生じるため、それらを

分極P の関数として求めなくてはならない。

そこで、最初に、電気双極子によってできる電場を求め、その電場

から、電荷密度 )(xeρ と表面電荷密度 )(xσ を求める。次に4つに分けた物質のそれぞれにおいて発生する電

場を求め、それらの和としての局所電場を求める。

外部電場Eにより物質内に誘起される分極P は、微小な電気双極子 pの集合体である。そこで、電気双極

子 pが1個で作る静電ポテンシャルを求める。電気双極子 pの電荷をq とし、その長さをsとする（ qq −+ ,

の2つの電荷がs離れて位置する。）と、 sp q= で与えられる。電気双極子の中間の点を原点に取り、そこか

らr離れた場所での静電ポテンシャル )(rφ は、

++−

−+=

+−

−=

θθπεπεφ

cos)2/(

1

cos)2/(

1

4)2/(

1

)2/(

1

4)(

222200 rssrrssr

qq

22srsr

r

である。 sr >> とし、上の式を rs / で展開し、 rs / の1次まで取る。 xxf += 1/1)( の1次までの近似は、

2/312/12/1 )1(2/1)1(2/1)(')1()( −−−− +⋅−=+⋅−=→+= xxxfxxf

である事を利用すると、 2/1)0(')0()( xxffxf −=+≈ から、

++−

−+=

++−

−+=

rsrsrsrsr

q

rssrrssr

q

/cos)2/(1

1

/cos)2/(1

1

4

cos)2/(

1

cos)2/(

1

4)(

220

22220

θθπε

θθπεφ

　　　

r

3

0

2

00

22

0

4

4

coscos

4]/cos)2/[(

2

11]/cos)2/[(

2

11

4

r

r

qs

r

s

r

qrsrsrsrs

r

q

πε

πεθθ

πεθθ

πε

rp ⋅=

=⋅=

+−−

−−≈

　　　

　　　

を得る。電気双極子のある点（2つの電荷の中点）を 'x 、観測点をxとすれば、 'xxr −= なので、ポテンシ

ャル )(rφ を次のように書き換える。

3

0 |'4

)'()',(

xx|

xxpxx

−−⋅

=πε

φ

この式を、 'x に関する微分演算 )'/,'/,'/(' zyxx ∂∂∂∂∂∂=∇ に関する下の関係式を用いて、式変形する。

( ) 12/1222222

222

''

)'()'()'()2/1()')(1(2)'()'()'(''

|'

)'()'()'(|')'

,'

,'

(grad

−−+−+−⋅⋅−−⋅=−+−+−

∂∂

=∂

−∂

−+−+−=−∂∂

∂∂

∂∂

=≡∇

zzyyxxxxzzyyxxxx

zzyyxxzyx

xx

xx|

xx|,　　

r

2/s+2/s− θ'x

4

3

1 2

誘電体

66

222 )'()'()'(

)'(

'

|'

zzyyxx

xx

x −+−+−

−−=

∂−∂ xx|

これより、

|'

'|''

xx|

xxxx|

−−

−=−∇ x

よって、以下の結果を得る。

32'2'|'

'

|'

)'(

|'

1|'

|'

1

|'

1:

'

)'(

)'(

))'(())'((

' xx|

xx

xx|

xx

xx|xx|

xx|xx|

x

x

xx

−−

=−−−

⋅−−

=−∇−−

=−

∇→∂

∂∂

∂=

∂∂

xxx

u

u

ufuf

x

この関係を用いて、ポテンシャル )(rφ は、次のように書き換えられる。（注意：観測点を x、電気双極子 pの

ある点を 'x とした。）

|'

1grad

4|'

1

4|'

)'(

4|'4

)'()',( '

0

'

0

3

0

3

0 xx|

p

xx|

p

xx|

xxp

xx|

xxpxx

−⋅=

−∇⋅=

−−

⋅=−

−⋅= xx πεπεπεπε

φ

誘電体中には、このような電気双極子 pが多数含まれる。その集合体である単位体積あたりの電気双極子が分

極を表わすベクトルP に他ならない。さらに )',( xxφ を分極 )'(xP を用いて書き換えると、電気双極子 pの

ある点を 'x としたので、

|'

1grad

4

)'()',( '

0 xx|

xPxx

−⋅= xπε

φ

である。これが、点 'x での分極によって出来る単位体積あたりの静電ポテンシャルである。なお、観測点は点

xである。

次に、分極P の空間的な分布の考察から、分極P に起因する物体内の

電荷密度 )(xeρ 、および表面電荷密度σ を求める。誘電体が作る電場

は、図のように、4 つの部分からなる。1 番目は、半径aでくり抜かれ

た誘電体の部分が、くり抜かれないで残っているために作る電場である。

残りの誘電体の部分による電場の寄与（図の 2～4の部分）を )(xDφ で

表わそう。

観測点xで生じる静電ポテンシャル )(xDφ は、誘電体中の点 'x での分極 )'(xP を用いると、誘電体中の全

ての点 'x に関して体積積分を行ない、

∫∫∫ −⋅=

−⋅==

V

x

V

x

V

D ddd '|'

1grad)'(

4

1'

|'

1grad

4

)'(')',()( '

0

'

0

xxx|

xPxxx|

xPxxxx

πεπεφφ

で与えられる。 '''' dzdydxd =x 。ここで、次の関係式を用いて、 )(xDφ を書き直す。

|'

1grad)'(

|'

)'(div

|'

)'(div

|'

1

')'(

|'

1)]'([

'|'

)]'([

'|'

)'(

':

|'

)'(div

''

'

''

'

'

xx|xP

xx|

xP

xx|

xP

xx|xP

xx|xP

xx|

xP

xx|

xP

xx|

xP

−+

−=

−

→

−∂∂

+−∂

∂=

−∂∂

=

−∂

∂

−

xx

x

xx

x

xxxxx

　

　　

この関係式を用いると、 )(xDφ は、以下の式で与えられる。

4

3

1 2

67

∫∫

∫∫

−−

+

−=

−−

−=

−⋅=

V

x

V

x

V

xx

V

xD

dd

dd

'|'

)'(div

4

1'

|'

)'(div

4

1

'|'

)'(div

|'

)'(div

4

1'

|'

1grad)'(

4

1)(

'

0

'

0

''

0

'

0

xxx|

xPx

xx|

xP

xxx|

xP

xx|

xPx

xx|xPx

πεπε

πεπεφ

　　　

ここで、右辺1番目の式にガウスの定理を用いて、体積積分を面積積分に変える。また、誘電体とそれを取り

囲む真空中との境界面においては、誘電体から真空中へ向かう方向に法線ベクトルの向きを取り、単位長さの

法線ベクトルを )(xn として計算する。すると、静電ポテンシャル )(xDφ は以下のように書き換えられる。

∫∫∫

∫∫

−⋅

+−

−+

−⋅

=

−−

+

−

=

)(0

'

0)(0

'

0

'

0

|'

)'()'(

4

1'

|'

)'(div

4

1

|'

)'()'(

4

1

'|'

)'(div

4

1'

|'

)'(div

4

1)(

outSV

x

inS

V

x

V

xD

dSddS

dd

xx|

xnxPx

xx|

xP

xx|

xnxP

xxx|

xPx

xx|

xPx

πεπεπε

πεπεφ

　　　

右辺第1項は、くり抜かれた空洞の部分の表面 )(inS の表面電荷による寄与（法線ベクトルnの向きは、誘電体

からくり抜かれた部分へ向かう方向）であり、第2項は、誘電体中の、分極による電荷の寄与である。第3項

は、誘電体の表面に誘起された表面 )(outS の表面電荷による寄与（法線ベクトルnの向きは、誘電体から外へ

向かう方向）である。

さて、空間内に電荷分布 )'(xeρ があり、誘電体表面に電荷密度 )'(xσ の電荷が分布する場合、xの場所に

おける静電ポテンシャル )(xφ は、

∫∫ −+

−=

V

e

S

ddS '|'

)'(

4

1

|'

)'(

4

1)(

00

xxx|

x

xx|

xx

ρπε

σπε

φ

である。 )(xDφ と )(xφ を比較すると、誘電体の表面 )()( , outin SS での電荷密度σ 、および分極による物体内の

電荷密度 )(xeρ は、法線ベクトルnの向きは、誘電体から外、またはくり抜かれた部分へ向かう方向とし、

)(div)(),()()( xPxxnxPx xe −=⋅= ρσ

で与えられる。これで、表面での電荷密度 )(xσ および物体内の電荷密度 )(xeρ の一般式が求まった。

注意：半径aの球の内部についても、同様の手続きで )(div)( xPx xe −=ρ が求まると考えるかもしれない。

しかし、分極 )(xP は巨視的な物理量なので、物質中の電気双極子の微視的な空間構造に起因する電荷密度を

求める事は、適当ではない。

以上の結果を元に、空洞にする前の誘電体中のある点 xにおける電場（局所電場）を求める。誘電体が作る

電場は、図のように、4つの部分からなる。番号 1の部分は、半径aでくり抜かれた誘電体の部分が、くり抜

かれないで残っているために作る電場である。残りの誘電体の部分による電場の寄与（図の番号2～4の部分）

は、 )(xDφ で表わされる。

最初に、半径aでくり抜かれた誘電体の部分が作る電場（図の1の部分）を考える。くりぬいた球の内部の

原子・分子が作る電場は、球を構成する結晶構造により異なるが、単純な構造の場合は、次の計算で示すよう

に、ゼロになる。一般の結晶構造の場合もほぼゼロになるとしよう。

球内の球の中心から位置ベクトルr離れた場所にある双極子 pが、球の中心に作る電場 dip-inE は、

⋅+−=

53

0

)(3

4

1)(

rrdip-in

rprprE

πε


68

電場 )(rE dip-inの求め方：電気双極子のある点を原点0 にとり、を観測点を ),,( zyx=r とする。電気双極子の静電

ポテンシャル )(rφ （p.65）のgradをとれば、双極子による電場が求まる。

∂

∂∂

∂∂

∂=

⋅−=−=→

⋅=

z

f

y

f

x

ff

rrdip-in

)(,

)(,

)()(grad,

4grad)(grad)(

4)(

3

0

3

0

rrrr

rprrE

rpr

πεφ

πεφ

電場の計算

( )

⋅+−=

⋅+−=→

⋅+−=

⋅+−=

∂∂⋅

−−−=

=⋅++=∂

++∂=

∂∂

∂∂

−⋅

−−=

∂∂⋅

−⋅∂∂

−=

⋅∂∂

−=

+++−

−

+−

5

0

3

0

5

0

3

0

5

0

3

0

113

0

3

0

)13(

0

3

0

12/12222/1222

)13(

0

3

0

3

0

3

0

3

0

4

)(3

4

1

4

)(3

4)(

43

443

4

1

4)3(

4

1)]([

2)(2

1)(

)3(44

11

44

1][

4)]([

rrrr

r

x

r

p

r

x

rp

x

rr

rp

r

xxzyx

x

zyx

x

r

x

rr

rp

rxrxrx

dip-in

jjj

j

j

jjdip-in

j

j

jj

j

j

jjj

jdip-in

πεπεπεπε

πεπεπεπεπεπε

πεπεπεπεπε

rprprrpprE

rprprprE

rprprp

rprE

　　　

球の対称性を考慮すると、電場 dip-inE の、外部電場に垂直な成分はゼロになる。電

場 dip-inE の外部電場と同じ方向の成分 Externaldip-in //)(E は、右図の記号を用いると、

3

2

0

53

0

//

)1cos3(

4

1)cos(cos3

4

1)(

r

p

r

prr

r

pExternaldip-in

−=

+−=θ

πεθθ

πεE

となる。電気双極子が半径a の球の空間（くりぬいた球の内部）に一様に分布していると仮定して、球内の全

ての双極子 pによる寄与を求める。次の式がその計算。 φθ ,,r が微小量 φθ dddr ,, 変わることで作られる微

小な体積は、 φθθ drrddr sin×× で与えられる（p.69の補足参照）ので、

[ ] 0)13())(13(sin)1cos3(:cos

0sin)1cos3(4

sin)1cos3(

4

1

1

1

3

1

1

2

1

1

2

0

2

2

00

2

0

3

2

0

3

2

0

=−=−=−−=−=

=−=⋅⋅−

−

+

−

−

+

+

∫∫∫

∫∫∫∫

xxdxxdxxdx

dddrr

rpdrrddr

r

pa

π

ππ

θθθθ

φθθθπε

φθθθ

πε

変数変換　　

すなわち、 0)( // =Externaldip-inE 。以上より、くりぬいた球の内部の原子・分子が作る電場は、ゼロである。

次に、誘電体中の分極による電荷密度 )(xeρ （図の3の部分）の電場への寄与を求める。電荷密度 )(xeρ は、

)(div)( xPx −=eρ

で与えられる。その電荷密度 )(xeρ で作られる静電ポテンシャルは、既に求めたように、

∫∫ −−

=−

V

x

V

e dd '|'

)'(div

4

1'

|'

)'(

4

1 '

00

xxx|

xPx

xx|

x

περ

πε

であった。外部電場は空間的に一様であり、分極ベクトル )(xP の空間依存性がない場合、電荷密度

0)(div)( =−= xPxeρ 。半径a の球の外側の部分での、分極に起因した電荷密度 )(xeρ で作られるポテン

シャルはゼロ（空間依存性がない）である。よって、この静電ポテンシャルを空間微分する事で得られる電場

＝ゼロである。

続いて、誘電体の外側の表面 )(outS に存在する表面電荷σ （図 4の部分）による、静電ポテンシャルへの寄

与を考える。これによる寄与は、

P

θ

E

'x

x

'xxr −=

69

∫ −⋅

)(0 |'

)'()'(

4

1

outS

dSxx|

xnxP

πε

で与えられる。この遠方の表面電荷による静電ポテンシャルは、着目している空洞の場所から見れば、誘電体

の結晶構造の詳細には関係ない、巨視的静電場である。誘起される分極P の向きは外部電場0E と同方向であ

るため、分極P によって出来る表面電荷で作られる電場 1E の向きは外部電場0E と逆向きになる。この電場 1E

を、反分極場という。誘電体内部の巨視的電場Eは、外部電場0E と反分極場

1E の和で与えられる。

ここで、誘電体が図に示すような長方形の場合に、誘電体内部の巨視的電場Eを求めよう。誘電体が長方形

の場合、反分極場1E の大きさと向きは、プラズマ振動の箇所で求めた

（p.59の補足）。それによれば、反分極場の大きさ 1E は、

00

1 εεPNp

E ==

で与えられた。ここで p は外部電場によって誘起された電気双極子 1個の大きさであり、N は単位体積あた

りの電気双極子の数（電気双極子の数密度）である。この電場1E は

0E と反対向きの電場ベクトルである。分

極の大きさP は NpP = で与えられ、P は 0E と同じ方向を向く。よって、長方形の形をした誘電体内部の巨

視的電場Eは、

0

010 εP

EEEE −=+=


最後に、切り抜かれた部分の表面電荷によって作られる電場を求めよう。一様な外部電場0E が加えられて

いる誘電体を連続体として扱い、微小な球をくりぬくと、その球の表面には、分極P に起因した表面電荷σ が

現れる。次は、この表面電荷σ により作られる電場PE を求める。

外部から加えられた一様な電場0E により作られる分極P は、外部

電場0E と向きが同じであるとする。球の表面電荷密度σ は、右図の

ように、誘電体から半径a の球がくりぬかれた状況を考えると、

θθπσ cos)cos( PP −=−=⋅= nP

で与えられる。図の(＋）（－）の記号は、分極P による表面電荷の符

号を示す。球座標表示をとれば、微小な面積は、 φθθ daaddS sin⋅= （すぐ下の補足参照）で与えられる。

補足：積分での微小面積と微小体積の求め方

図で、右方向に z 軸を取り、 z 軸を周回する方向に角度φ を取り、

z 軸の+方向から－方向に移動する向きに角度θ をとる。極座標で、

半径aの球で角度 φθ , の微小な変化で出来る微小面積dS を考える。

半径aの変化量 θd の回転で、 θad の幅を持つ。また、半径 θsina の

角度 φd 回転により、幅 φθdasin になる。これら2つの長さで出来る

微小面積は、 θφθ addadS ×= sin である。微小体積dV の場合は、

ra → （動径）と読み換える。さらに動径r がdr だけ変わると、

φθ dddr ,, の変化のため、微小体積は、 dSdrdV ×= すなわち、

+ -

+ -

+ -

1E 0E

nP

θ

E

+ －

+ －

+ －

+ －

+ －

θ

a

adθ

z φ

adθ

dθ

0

asinθ

asinθdφ

70

θφθ rddrdrdV ××= sin になる。

この微小な面積dS における電荷の大きさは、 )sin( φθθσσ daaddSdq ⋅== となり、この微小電荷の、球

の中心との距離はa である。この微小電荷で出来る電場の、分極P に垂直な成分は、くりぬかれた球の対称性

を考えるとゼロになる。そこで、微小部分で出来る電場の、分極P の方向の電場成分のみを考慮する。球の中

心から見ると、図の角度θ をなす部分で出来る電場の向きは、法線ベクトルnの向きになる。そこで、分極P

と平行な電場成分の大きさは、電場の大きさに )cos( θπ − を掛けたもので与えられる。

以上の考察から、分極P に起因する表面電荷σ によって作られる電場PE の、P 方向の大きさ PE は、次式


00

3

00

22

0

2

2

0

2

0

2

0)(

2

0

33

cos2

4

1cossin

4

1

cos)sin)(cos(4

1cos

4

1

4

)cos(

εθ

ππε

θθθφπε

θφθθθπε

θσπεπε

θπ

π

θ

ππ PPddPa

a

daadPa

dSaa

dqE

inS

P

+=

−⋅+=⋅+=

⋅−−=−=−

=

=

ΩΩ

∫∫

∫∫∫

　　

よって、この表面電荷による電場は

03εP

となる。

以上の結果をまとめると、局所電場localE は、誘電体内部の巨視的電場をEとすれば、ベクトル表示で、

03εP

EE +=local

となる。この結果は、仮定した球の半径a によらない。この局所電場 localE を、ローレンツの局所電場という。

クラウジウス-モソッティ（Clausius-Mossotti）の式

物質に加わる巨視的電場Eと分極P の間に、 EPPEED )( 00 εεεε −=→+== の関係があるので、

EEEP

EE0

00

00 3

2)(

3

1

3 εεε

εεεε

+=−+=+=local

を得る。一方、局所電場localE が誘電体内の1個の原子・分子に加わると、原子・分子は分極して電気双極子 p

を作る。 pは localE に比例するので、次のように書ける。

localEp α=

α を原子・分子の分極能率という。分極P は pと単位体積あたりの原子（分子）数N の積で与えられる。

localNN EpP α==

よって、次の下線を引いた関係式（クラウジウス-モソッティの式）を得る。

00

0

0

00

0

0

0

323

2)(

)(

3

2

εα

εεεε

εεε

αεεεε

εεε

αα NN

NN local=

+−

→+

=−→

−=

+==

EE

EP

EEP

71

電子振動子模型における、局所電場の効果

ここでは、電子振動子模型において、局所電場の効果を見よう。1 次元の電子振動子の運動方程式は、外部

電場を )(tE として、

)()()()(

2

2

tEm

etx

m

k

dt

tdx

dt

txd−=++ γ

で与えられた。このような、原子群の集団としての運動を考えよう。ここで )(tx は電子の平行位置からの変位

であり、γ は減衰定数、kx は復元力、 )(tE は電場の瞬時値、 )( e− は電子の電荷を表す。外部電場 )(tE に対

して、全ての電子が同じ振幅、同位相で振動すると仮定する。簡単のため、原子群の巨視的な電場を外部電場

)(tE に等しいと見なす。電場 )(tE を、局所電場に置き換える。原子数密度をN とすると、

)()(3

1)()(

3 00

txeNtEtElocal −+→↔+=εε

PEE

すると、電子の運動方程式は、次のように書き換えられる。

)()(3

)()(

)]()(3

1)([)(

)()()()(

)()(

0

2

2

2

0

2

2

2

2

tEm

etx

m

Ne

m

k

dt

tdx

dt

txd

txeNtEm

etx

m

k

dt

tdx

dt

txdtE

m

etx

m

k

dt

tdx

dt

txd

−=

−++→

−+−=++→−=++

εγ

εγγ

局所電場の効果を考慮しない場合との違いは、運動方程式において、電子の束縛（ばね定数）に関する項が

−=

−=−→=

0

22

0

0

2

0

22

03

13

13 ε

ωεε

ωk

Ne

k

Ne

m

k

m

Ne

m

k

m

k

となり、電気双極子の存在によりばね定数が実効的に小さくなる事である。電気双極子の密度N に依存した電

子の共鳴周波数は、次式で与えられる。なお、N が小さな希薄な気体では、 00 ' ωω ≈ である。

0

2

003

1'ε

ωωk

Ne−=

なお、電子がばらばらな位相で振動している場合、巨視的な電気振動子（分極）の統計平均としての値 )(tP は、

0)()()()(11

=−== ∑∑ ==

N

j j

N

j j txetptP となるため、周波数の低減は見られない。

72

クラマース・クロニッヒ（Kramers-Kronig）の関係式

誘電体に電場E を加えると、分極P を生じるが、E に対するP の応答には、一般に時間の遅れが生じる。

よってE が時間的に変動する場合には、 )(tE と )(tP との間に単純な比例関係が成り立たない。そこで、 )(tE

を入力、 )(tP を応答と考え、その間の一般的な関係を求める。まず、原理的な考察により、

(i) 入力 )(),( 21 tEtE に対し、それぞれ )(),( 21 tPtP が生じるなら、入力 )()( 21 tbEtaE + に対し、

)()( 21 tbPtaP + が生じる。（重ね合わせの原理）。

(ii) 入力 )(tE に対し )(tP が生じるなら、 )( attE − に対し )( attP − が生じる（時間に対する一様性）。

(iii) att < で 0)( =tE （入力信号＝ゼロ）なら、 att < で 0)( =tP （因果律）。

（i）の重ね合わせの性質から、入力 )(tE に対する応答 )(tP は、 )(tE に関して線形性がある。よって、

∫+∞

∞−= ')'()':()( 0 dttEttGtP ε

と置く事ができる。、（ii）（iii）により、

0,0)(

0'),()':(

<=

≥−==

ττττ

　G

ttGttG

となる。 )(τG を応答関数という。 )()( ttE δ= とすると、これに対する応答は、 )()( 0 tGtP ε= である。

ここで、 )exp()( 0 tiEtE ω−= とすると、

)()(

',)exp()()exp(

')exp())'(exp()'(')'exp()'()(

0

0

00

-

00

-

00

tE

ttdiGtiE

dttittittGEdttittGEtP

ωχε

ττωττωε

ωωεωε

=

−=−=

−−+−=−−=

∫

∫∫∞+

∞∞

　　　

　　　　　

ｔｔ

ただし、複素数で定義された電気感受率 )(ωχ は

∫+∞

=0

)exp()()( τωττωχ diG

である。 )("),(' ωχωχ を実数として )(")(')( ωχωχωχ i+= とおくと、すぐ上の式から、

∫∫+∞+∞

+==+=00

)]sin())[cos(()exp()()(")(')( τωτωτττωττωχωχωχ diGdiGi

となるので、実数部分と虚数部分同士を比較し、

∫∫+∞+∞

==00

)sin()()(",)cos()()(' τωττωχτωττωχ dGdG 　　

である。また上の式で ωω −→ とおけば、 )sin()sin(),cos()cos( ωτωτωτωτ −=−=− の性質から、

)(")("),(')(' ωχωχωχωχ −=−=−

が成立する。

フーリエ（Fourier）変換の公式

よく用いられる係数を含めた公式

73

∫∫+∞

∞−

+∞

∞−

−== ωωτωπ

ττωττπ

ω diFfdifF )exp()(2

1)(,)exp()(

2

1)( 　

において、 )(2)( ωπω FR = とおくと、

∫∫∫

∫∫∞+

∞−

∞+

∞−

∞+

∞−

+∞

∞−

+∞

∞−

−=−=−=

==→=

ωωτωπ

ωωτωππ

ωωτωπ

τ

τωττωπωτωττπ

ω

diRdiFdiFf

difFRdifF

)exp()(2

1)exp()(2

2

1)exp()(

2

1)(

)exp()()(2)()exp()(2

1)(

これにより )(,)( τω fR 　間でフーリエ変換を再定義すると、

∫∫+∞

∞−

+∞

∞−

−== ωωτωπ

ττωττω diRfdifR )exp()(2

1)(,)exp()()( 　

となる。

続いて、 )(),( τω fR の偶奇性の考察をする。

もし )()( ωω RR =− （偶関数）なら、

∫

∫∫∫∞+

+∞

∞−

+∞

∞−

+∞

∞−

=

=+=−=

0

)cos()(1

)cos()(2

1)]sin())[cos((

2

1)exp()(

2

1)(

ωωτωπ

ωωτωπ

ωωτωτωπ

ωωτωπ

τ

dR

dRdiRdiRf

　　　

となる（フーリエコサイン変換）。その時、

∫∫+∞+∞

∞−

==0

)cos()(1

)cos()(2

1)( ωωτω

πωωτω

πτ dRdRf

において ττ −→ とおけば、 )()( ττ −= ff 、すなわち )(τf は偶関数である事が分かる。

一方、 )()( ωω RR −=− （奇関数）の時、

∫

∫∫∫∞+

+∞

∞−

+∞

∞−

+∞

∞−

−=

−×=−=−=

0

)sin()(

)sin()()(2

1)]sin())[cos((

2

1)exp()(

2

1)(

ωωτωπ

ωωτωπ

ωωτωτωπ

ωωτωπ

τ

dRi

diRdiRdiRf

　　　

となる（フーリエサイン変換）。ここで ττ −→ とおけば、 )()( ττ −−= ff 、すなわち )(τf は奇関数である

事が分かる。

以上をまとめると、

)(ωR と )(τf のフーリエ変換の関係は、

∫∫+∞

∞−

+∞

∞−

−== ωωτωπ

ττωττω diRfdifR )exp()(2

1)(,)exp()()( 　

)(ωR が偶関数の時、 )(τf も偶関数であり、 ∫+∞

=0

)cos()(1

)( ωωτωπ

τ dRf 。

)(ωR が奇関数の時、 )(τf も奇関数であり、 ∫+∞−

=0

)sin()()( ωωτωπ

τ dRi

f 。

74

これらの性質を用いて、 )(")(')( ωχωχωχ i+= （ )("),(' ωχωχ は実数）と )(τG に関する考察をする。

さて、

∫∫

∫∞+∞+

+∞

==

=+=

00

0

)sin()()(",)cos()()('

,)exp()()(")(')(

τωττωχτωττωχ

τωττωχωχωχ

dGdG

diGi

　　

　

において、積分区間を +∞<<∞− τ に変更する。その際、フーリエコサイン変換、サイン変換では、 )(τG の

範囲を +∞<< τ0 から +∞<<∞− τ に拡張し、それぞれ偶関数、奇関数と解釈しなおした計算になる。

まず、 )(τG を仮に偶関数と見なして、 )(' ωχ と )(τG との関係を拡張する。

∫+∞

=0

)cos()()(' 　τωττωχ dG

に対し、仮に )(τG が 0<τ でも定義され、 )()( ττ GG =− の性質を持つと仮定すれば、偶関数どうしのフー

リエ逆変換を用いると（注意： )(")("),(')(' ωχωχωχωχ −=−=− を既に示した。）、

∫∫+∞+∞

=↔=00

)cos()('1

)()cos()(1

)( 　τωτωχπ

τωωτωπ

τ dGdRf

が得られる。 )(" ωχ についても同様である。 )(" ωχ の場合は、仮に )()( ττ GG −=− の奇関数の性質がある

として、奇関数のフーリエ変換の計算をすればよい。これらをまとめると、

∫∫+∞+∞

==00

)sin()("1

)cos()('1

)( 　　 τωτωχπ

τωτωχπ

τ ddG

である。

次に、 )(' ωχ と )(" ωχ との関係式を求める。これ以降は、複素関数論の本を参照する事。複素積分を行うた

めにω を複素数 'ω でおき換え、 )'(ωε を複素平面上で定義された関数と考えよう。

∫∫+∞+∞

→=00

)'exp()()exp()()( ττωττωττωχ diGdiG

',' 21 ωω を実数として ''' 21 ωωω i+= で置き換えると、次の式を得る。

( )

∫

∫∫∞+

+∞+∞

−=

+==

0

21

0

21

0

)'exp()'exp()(

)''(exp)()'exp()()'(

ττωτωτ

ττωωτττωτωχ

diG

diiGdiG

　　　

)(τG は、 +∞<< τ0 で定義されている。もし 0'2 <ω なら、 +∞→τ の時、 +∞→− τω '2 となるため、

+∞→− )'exp( 2 τω になる。よって、 )(τG が発散しない自然な関数であるなら、 )'(ωχ が発散する異常な振

る舞いをするのは、 0'2 <ω の場合に限られる。逆にいえば、「 )'(ωχ は 0'2 >ω の複素平面の上半分で正則」

である。この結論は、後ほど、総和側を導く際に用いる。

さて、コーシー（Cauchy）の積分定理により、

∫ ∫∫∫∫ =−++=− ∞+−∞

+∞

−∞0'

'

)'(

DCCA

A

dω

ωωω

ωχ

75

ここで +∞∞− ~~ D の部分の経路積分で、円の半径を十分大きくとると、

0=∫ ∞+∞− D

になる。続いて、次の公式（一周する円の経路でのθ積分は、反時計回り（左回り）が＋になる）を用いる。

)(2)(

lim)()(

)()(0

2

0

0

000

0

0

zifidee

ezfed

zez

ezfdz

zz

zf i

i

i

C

i

i

i

Cπθε

εε

εε

ε πθ

θ

θ

ε

θθ

θ

=+

=−+

+=

− ∫∫∫ →

すると、経路 CAω での積分は、右回り（ )1(−× ）半円（ 2/1× ）の積分のため、 2/)(2)1( ωχπi×− で与え

られる。これから、

)(''

)'(

)()1(''

)'(0

ωχπωωω

ωχ

ωπχωωω

ωχω

id

id

C

A

C

A

CCA

A

=−

=+→

+×−+=++=−

=

∫∫∫

∫∫∫∫∫∫∞+

−∞

∞+

−∞

+∞

−∞

+∞

−∞

P　　

をえる。ここで、PPPPはコーシー（Cauchy）の積分の主値である。 )(")(')( ωχωχωχ i+= とおくと、

))(")('(''

)'(")'(')('

'

)'(ωχωχπω

ωωωχωχ

ωχπωωω

ωχiid

iid +=

−+

→=− ∫∫

+∞

−∞

+∞

−∞PP

実数部、虚数部を比較して、

''

)'('1)(",'

'

)'("1)(' ω

ωωωχ

πωχω

ωωωχ

πωχ dd ∫∫

+∞

−∞

+∞

−∞ −−=

−= PP 　

を得る。この関係式は応答関数の実部と虚部を結びつけるので、分散公式（クラマース-クロニッヒ

Kramers-Kronigの式）とよぶ。この関係式は、 )("),(' ωχωχ の一方を知っている場合に、他方を導くのに

用いられる。

さて、 )(")("),(')(' ωχωχωχωχ −=−=− の性質を用いると、以下のように変形できる。

∫∫

∫∫∫

∫∫∫∫

∞+∞+

∞+∞+∞+

+∞

∞

+∞

∞

−=

−=

−

++

=−

++

−=

−+−

−−−

=−

+−−−

−=

022

022

000

4

0

0

0

0

''

)'("'2'

'

)'("'21

')'("'

1

'

11'

'

)'("1'

'

)'("1)1(

''

)'("1)'(

'

)'("1'

'

)'("1)'(

)'(

)'("1)('

　　

　　　

　　　

　　　

　　　

ωωωωχω

πω

ωωωχω

π

ωωχωωωωπ

ωωω

ωχπ

ωωω

ωχπ

ωωω

ωχπ

ωωω

ωχπ

ωωω

ωχπ

ωωω

ωχπ

ωχ

dd

ddd

dddd

PP

PPP

PPPP--

∫∫

∫∫∫∫

∫∫∫

∞+∞+

∞+∞+∞+

∞+

+∞

∞−

+∞

∞

−−=

−

++

−=

−−

+−

−=−

−−−−

−−=

−−+

−−=

−−=

022

0

00

4

0

0

0

0

''

)'(''2')'('

'

1

'

11

''

)'('1'

'

)'('1)1('

'

)'('1)'(

)'(

)'('1

''

)'('1'

'

)'('1'

'

)'('1)("

　　

　　　　

　　

　　　

　　　

ωωωωχω

πωωχ

ωωωωπ

ωωω

ωχπ

ωωωωχ

πω

ωωωχ

πω

ωωωχ

π

ωωω

ωχπ

ωωω

ωχπ

ωωω

ωχπ

ωχ

dd

dddd

ddd

PP

PPPP

PPP-

すなわち、

∫∫+∞+∞

−−=

−=

022

022

''

)'(''2)(",'

'

)'("'2)('

　　

ωωωωχω

πωχω

ωωωχω

πωχ dd PP

76

となる。

さて、 ))(")('()()( 000 ωχωχεωχεεωε i+==− の関係と上で得られた式から、 )("),(' ωεωε の間に成

り立つ関係式を導く。この際、 )()( 00 ωχεεωε =− を利用し、（注意： )(")(')( ωεωεωε i+= ）

))(()(")]("))([(''

)'("))'('(

))((''

)'()('

'

)'(()('

'

)'(

000

00

00

0

εωεπωπεωεεωεπωωω

ωεεωε

εωεπωωω

εωεωχεπω

ωωωχε

εωχπωωω

ωχ

−+−=+−=−

+−→

−=−

−→=

−→×−=

−

∫

∫∫∫∞+

∞

+∞

∞

+∞

∞

+∞

∞

iiidi

ididid

-

---

P

PP)P

　

　　　

　　

実数部、虚数部を比較して、

''

)'('1)(",'

'

)'("1)(' 0

0 ωωω

εωεπ

ωεωωω

ωεπ

εωε dd ∫∫+∞

∞

+∞

∞ −−

−=−

=−--

PP　　

　

を得る。

1次元の電子振動子での、 )(")("),(')(' ωχωχωχωχ −=−=− の性質の確認をしよう。1次元の電子振動子

の運動方程式は、 )exp()(/ 00 tiEtEmk ωω −== 、で、

)()()()(

2

2

tEm

etx

m

k

dt

tdx

dt

txd−=++ γ


ωγωωω

i

tiEmetx

−−

−−=

22

0

0 )exp()/()(

電気双極子 )(tp は、

)exp()/(

)()()()()( 022

0

0 tiEi

metEtxetp ω

ωγωωωχε −

−−−==−=

で与えられる。よって )(")(')( ωχωχωχ i+= は、

2222

00

2

2222

0

22

0

0

2

2222

0

22

0

0

2

22

00

2

22

0

0

2

)()()()(

)()(

1)/()(

ωγωωωγ

εωγωωωω

ε

ωγωωωγωω

εωγωωεωγωωε

ωχ

+−+

+−

−=

+−

+−=

−−=

−−=

m

ei

m

e

i

m

e

im

e

i

me

　　　　

これより、 )(")("),(')(' ωχωχωχωχ −=−=− の関係が成り立っている事が確かめられた。なお、複数の振

動子が存在する場合は、単位体積あたりの振動子の数密度N 、k 番目の振動子の共鳴周波数 kω と振動子強度

kf を用いて

　　　　 1,)/(

)(22

0

2

=−−

= ∑∑k

k

k k

k fi

fmeN

ωγωωε

ωχ

と書く事ができる。

屈折率 )(ωn と消衰係数 )(ωκ との間の関係式

続いて、屈折率 )(ωn と消衰係数 )(ωκ との間の関係式を導く。われわれは、次の関係式を、損失がある場

合の電磁波の伝播の箇所で、既に導いた。

77

00

2),()()(1

'

cin

ωκαωκωωχ

εε

=+=+=

さて、「 )(τG が発散しない自然な関数であるなら、 ''' 21 ωωω i+= とおいて、 )'(ωχ が発散するような異常な

振る舞いをするのは、 0'2 <ω の場合に限られる。逆にいえば、 )'(ωχ は 0'2 >ω である複素平面の上半分で

正則」である」事を議論した（p.74下）。 )'(ωχ がこのような特性であるなら、 )(1)(1)(1 ωκωωχ in +−=−+

も、 0'2 >ω である複素平面の上半分で正則になる事がわかる。また、 1|| >>ω の領域では、物質が共鳴する

周波数から十分離れているので、 0)( →ωχ すなわち 0)(1)( →+− ωκω in である。よって、 )("),(' ωχωχの関係式を導いた手続きを、 )(,1)( ωκω −n でもそのまま使える。対応関係は、 1)()(' −↔ ωωχ n 、および

)()(" ωκωχ ↔ であるから、以下の関係式を得る。

∫

∫∫∫∞+

∞+∞+∞+

−−

−=

−=

−=

−=−

0

22

0

22

0

0

22

0

0

22

''

]1)'(['2)(

''

)'('

'

'2

)'('

2'

'

)'('21)(

ωωω

ωωπ

ωκ

ωωω

ωαπ

ωωωωωα

ω

πω

ωωωκω

πω

dn

dc

d

c

dn

P

PPP

)(),( ωεωχ の総和則

電子振動子模型に限らず、電子（分極）の振動は、外部から掛けられた電場の周波数 ∞→ω の時、追従で

きないと考えられる。

さて、複数の電子振子が存在するとし、 0ω の最大値を Mω としよう。 )(" ωχ は Mωω >> の時、電磁波の

吸収がないため、ゼロになる。この関係を用いて、次の )(' ωχ の積分表示において、 ',ωω の大小関係で次の

ように分けて考える。

∫+∞

−=

0

22'

'

)'("'2)(' ω

ωωωχω

πωχ dP

'ωω >> の時は、

222

)'("'~

'

)'("':'

ωωχω

ωωωχω

ωω−−

>>

とできる。一方、 Mωωω >>>' の時、222 '' ωωω ≈− と見なす事が出来る。また、 'ω は共鳴周波数の最大値

から大きく外れているので、 0)'(" →ωχ である。その時、以下のように大きさを評価することが出来る。

0'

0'~

'

)'("':,'

222=

−⋅

−>>

ωω

ωωωχω

ωωω M

よって、被積分関数を次のように書き換えても問題はないだろう。

222

)'("'

'

)'("':,'

ωωχω

ωωωχω

ωωω−

→−

>> M

次に、 'ωω >> 、 Mωωω >>>' における不等号を[ >>→> ]のように解釈しなおすと、 Mωω >> の領域

で、次の近似式を得る。なお、この式変形の最後の式は、発散する事がないので、主値を取ることを意味するP

を除いて書いた。

78

∫

∫∫∫

∫∫∫

∞+

∞+∞+

+∞+∞

−=

−=

−+

−≈

−+

−=

−=

02

022

02

220

220

22

')'("'2

')'("'2

')'("'2

')'("'2

''

)'("'2'

'

)'("'2'

'

)'("'2)('

ωωχωπω

ωω

ωχωπ

ωω

ωχωπ

ωω

ωχωπ

ωωωωχω

πω

ωωωχω

πω

ωωωχω

πωχ

ω

ω

ω

ω

d

ddd

ddd

　　　

　　　 PPP

PPP

注意：「 )(τG が発散しない自然な関数であるなら、 )'(ωχ が発散する異常な振る舞いをするのは、 0'2 <ω の

場合に限られる。逆にいえば、 )'(ωχ は 0'2 >ω の複素平面の上半分で正則である」事を、既に議論している。

次に、 )("),(' ωχωχ の関係式 ∫+∞

−=0

2')'("'

2)(' ωωχω

πωωχ d を用いて、 )(),( ωαωn 間の関係式を求める。

まず、既に得られた関係式

00

2,)(1

'

cin

ωκακωχ

εε

=+=+=

に対して、 1|)(| <<ωχ が成り立つ場合は、 )()(" ωαωχ 、の間に

)(2

)("2

1,1)()('

2

1

2

)(1)()(

2

1

2

)()()()()(

2

11~)(1

'

0

00

0

ωαω

ωχωωχ

ωωα

ωωχωωα

ωωκωωχωχεε

cn

cin

cinin

=−=→

+−=→+=+=++=

　　

の関係式を得る。この結果を、先ほどの )("),(' ωχωχ の関係式に代入すると、次式を得る。

∫∫∫∫+∞+∞+∞+∞

=→=⋅==00

2

0

0

0

0

0

2

')'(2

')'(')'('

'')'("'2

ωωαπω

ωωαωωαω

ωωωχωπω

dc

dcdc

d PP

このようにして得られた関係式：

∫∫+∞+∞

==00

2

0

2

')'(2

,')'("'2

ωωαπω

ωωχωπω

dc

d PP

は、一般に総和則と呼ばれる。これら関係式は、吸収スペクトル測定に際して、未測定の波長領域において大

きな吸収があるかどうか判断する材料として用いられる事がある。

電子振動子模型での総和則を求めよう。共鳴周波数の最大値を Mω として、 Mωω >> が成り立つ時、

m

NemeNf

meNfmeNP

P

k

k

k

k

0

22

2

2

2

0

2

2

0

2

2

0

2

,1)/()/()/(

)(ε

ωωω

ωε

ωε

ωε

ωχ =−=⋅−=−

=−

≈ ∑∑

をえる。両者の比較を行うと、

∫∫∫ ∞+∞+

+∞

=→−=−→

−=≈

−=

0

2

0

22

2

2

2

0

2

')'("'2

')'("'2

)(')(

')'("'2

)('

ωωχωπω

ωωχωπωω

ω

ωω

ωχωχ

ωωχωπω

ωχdd

dPP

P

が、電子振動子模型での総和則である。

79

結晶中での電磁波の伝播光学的異方性・複屈折

電磁波の伝播では、物質は、等方的と仮定した。したがって、誘起された分極P の方向は、電場E の方向

と一致し、分極と電場とは、方向に依存しない定数の係数で結び付けられた（ EP χε 0= ）。しかし誘電体結

晶では、一般にこのような状況にはない。原子間の距離は、全ての方向で同じとは限らない。また、原子また

はイオン間の力の大きさも変わるだろう。すると、原子に固有な振動の周波数が異なり、それに伴い屈折率の

方向依存性が出てくるだろう。このような方向依存性は、「電子振動子模型」で説明した事柄から、予想される

事である。この状況は、異方性がない気体・液体と異なる。

光学的異方性に関連した現象のうち、光学素子に応用されている簡単な現象について、まず考察する。その

後、複屈折における簡単な場合について述べる。結晶は原子・イオンの周期的配列で出来ているので、誘起さ

れる分極の向きと大きさは、加えられた電場の向きと結晶の配置に依存すると予想される。そこで、

=

=

z

y

x

z

y

x

E

E

E

P

P

P

333231

232221

131211

00 ,

χχχχχχχχχ

εε EP χχχχ

の式が成立する。電気感受率 )( ijχ=χχχχ は一般に複素数の成分を持つので、以下では )( ijχ=χχχχ を複素感受率

テンソルと呼ぶ。 ijχ の値は、一般に結晶格子に対する zyx ,, 軸の取り方に依存する。そこで簡単のため、（立

方体を含めた）直方体の形の結晶を考えよう。結晶軸にそろえて座標軸を取れば、 zyx ,, の3つの軸が互いに

直交するような軸になる。また、非対角成分はすべてゼロになり、

=

z

y

x

z

y

x

E

E

E

P

P

P

33

22

11

0

00

00

00

χχ

χε

となる。一般には、任意の構造の結晶では、3 つの独立なベクトルを複素感受率テンソルの固有ベクトルとす

れば、複素感受率テンソルの固有値が複素感受率テンソルの対角成分になる。このあたりの話は、行列の対角

化・固有ベクトルを扱った線形代数などの本を参照する事。

)( ijχ=χχχχ ：複素感受率テンソルの要素について、結晶構造が簡単な場合について、説明する。

・立方体（全ての辺の長さが等しい結晶構造）

この場合、 zyx ,, 軸を区別する事は難しい。言い換えれば、どの結晶軸を zyx ,, 軸にとってもよい。よって

332211 χχχ == となる。

・単軸結晶、1軸性結晶(uniaxial crystal)

Tetragonal、hexagonal、trigonalな対称性の結晶の総称。2つの辺の長さが等しく、残りの辺の長さが異なる。長

さが異なる軸を一般にはz軸とする。光学的性質は、 yx, 方向では同じでz方向では異なるので、

332211 χχχ ≠=

である。 z 軸を光学軸、c軸（optic axis）という。

80

以下では、簡単のため、各軸が互いに直交する単軸結晶（1 軸性結晶）の場合について、光学的異方性と、

それに関連した現象を説明する。それぞれの軸方向の誘電率または電気感受率テンソル ijχ と屈折率との関係

は、以下の関係式から導かれる。

=

33

22

11

00

00

00

χχ

χχχχχ

の時、電束密度Dは、

+

+

+

=

+

=+=

=

+=+=+==

33

22

11

33

22

11

033

022

011

0

000

100

010

001

00

00

00

100

010

001

/

/

/

)(

χ

χχ

χ

χχ

ε

εε

ε

εεεε

χχχχ

εεεε

εεεεεεεε

εεεε

χχχχχχχχ

I

EIEEPEED 0

となる。さらに結晶は光電場を吸収しない（ 0=κ ）とする。それぞれの軸方向の屈折率（実数）の2乗が「物

質の誘電率を真空の誘電率で割った値： 0/ εε 」と等しい（ 22

0 )(~/' κεε inn +== 、 0=κ ）。よって、1,2

軸方向（光電場が直線偏光しているとして、偏光方向が1軸または2軸に平衡）の屈折率を 0n 、3軸方向（光

電場が直線偏光しているとして、偏光方向が3軸に平行）の屈折率を en とすれば、

2

33

2

2211 1,11 eo nn =+=+=+ χχχ

が成り立つ。以下では、単軸性結晶において、光電場が x 軸方向に伝播する場合を考える。

光電場の偏光の変化

x+ 方向へ伝播する平面波の光電波を考え、入射電場が次の式で表される場合を考える。真空中では、何れ

の偏光方向でも波数は同じだから、 0/ ck ω= として、 0<x で

)](exp[),0,0()](exp[)0,,0(

)](exp[)](exp[),(

00

00

kxtiEkxtiE

kxtiEkxtiEtx

zy

zzyy

−−+−−=

−−+−−=

ωω

ωω

　　　

eeE

である。 zyx eee ,, は、結晶の軸方向を示す、互いに直交する単位ベクトルである。電場は、x 軸方向に伝播し、

単軸性結晶内を進む。また、電場は、 zy, 成分を持つ。

電場が結晶面に垂直に入射する位置を 0=x とすると、それ

ぞれの偏光方向での波数ベクトルは、 0/ cnk ω= の一般的

な関係から、単軸性結晶において、以下のようになる。

00

0 ,c

nkc

nk ezy

ωω== 。

)( zy kk は、光電場の偏光方向が )(zy 方向の波数の値である。ここでは、 0nne > とする。すると、結晶に

ez

ex ey

k // ex

k // ex

x=0 L

81

入射後の電場は、

)](exp[)](exp[),( 00 xktiExktiEtx zzzyyy −−+−−= ωω eeE

ここで全体を )](exp[0

xc

nti o

ωω −− でくくると、電場 ),( txE は

])(exp[)()],(exp[)(),(0

000

0

0 xc

nniEExxc

ntixtx oezzyyo

ωωω −+=−−= eeEEE

で与えられる。すると、結晶の長さL 伝播後の電場は、以下のようになる。

])(exp[)()],(exp[)(),(0

000

0

0 Lc

nniEELLc

ntiLtL oezzyyo

ωωω −+=−−= eeEEE 　 ---（議論で用いる式）

結晶から出射した光電場の偏光特性は、光電場の z成分が )]/)((exp[ 0cLnni oe ω− の因子を持つため、一般

に入射電場の偏光特性とは異なる（必ずしも直線偏光のままで結晶から出射するとは限らない。）。以下では、

すぐ上の式を用いて、光電場の偏光特性が結晶伝播に伴い変化する様子を、 ππω ,2/)/)(( 0 =− cLnn oeの 2

つの場合で見る。簡単のため、　0000 >== EEE zy として、結晶の zy, 軸と45度の角度をなす直線偏光の入

射光電場を考える。

πω

=−0

)(c

Lnn oe

の場合

)](exp[)(

)](exp[)exp(

)](exp[)(),(

1)exp()]/)((exp[)/)((

0

00

0

00

0

0

00

Lc

ntiEE

Lc

ntiiEE

Lc

ntiLtL

icLnnicLnn

ozy

ozy

o

oeoe

ωω

ωωπ

ωω

πωπω

−−−+=

−−+=

−−=

−==−→=−

ee

ee

EE

　　　

　　　

となる。この場合、z方向の光電場の向きが入射時と逆向きになるので、ベクトルの足し算を考えると、結晶

から出てくる光電場の偏光方向は、入射電場と比べ 90 度回転（直交）している。図では、最初点線の向きの

偏光方向が、結晶を通過後、実線の向きに変わる事を示す。

結晶から出てくる光電場の時間変化を見よう。そのため、 0)/( 0 =− Lo tcLn ωω となるような時間 Lt を考

え、その後の時間経過 Lttt −=' として、時刻 't での電場の大きさを見る。電場の表示は実数表示にとる。

（注： )sin()cos()exp( xixix += ）

( )[ ]

)'cos()()'cos(

'exp)(Re)](exp[])(exp[Re),(

00

00

00

00

tEtE

tiEELc

ntiLc

nniEEtL

zy

zyooezy

ωω

ωω

ωω

−+=

−−+=

−−

−+=

ee

eeeeE

　　　

　

ez

ey E

ez

ey E

偏光方向のシフト

82

直線偏光の光電場の偏光方向を任意に回転させる方法

上の説明では、入射光電場が zy, 軸と 45度の角度をなす偏光方向の場合を考えたが、直線偏光の光電場の

偏光方向を結晶に対して適当な角度にとる事で、光電場の偏光方向を任意に回転させることが出来る。この事

を示すため、入射光電場として

)](exp[cossin),( 00 kxtiEEtx zy −−+= ωθθ eeE

を仮定する。 πω =− )/)(( 0cLnn oe の時、

)](exp[cos)(sin),(0

00 Lc

ntiEEtL ozy

ωωθθ −−−+= eeE

となり、入射時に z 軸（+）方向に対してθ の角度をなしていた直線偏光の光電場の向きが、出射時には z 軸

（－）方向にθ の角度をなす。すなわち、直線偏光の光電場と単軸結晶の光学軸とのなす角度を適当に選ぶ事

で、光電場の直線偏光はそのままで、任意の偏光方向に変える事が出来る。1/2 波長板は、この原理を利用す

る。

2/)(0

πω

=−c

Lnn oe

の場合

分かりやすくするため、光電場の実数表示をする。 zy, 方向の光電場の振幅が 00 >E である光電場が、x 方

向に伝播する場合を考える。結晶の長さをL とし、 0)/( 0 =− Lo tcLn ωω となるような時間 Lt を考え、その

後の時間経過 Lttt −=' での光電場の大きさを見る。（注： )sin()cos()exp( xixix += ）

( )[ ] ( )[ ])2/'cos()'cos(

2/'(exp)'exp(Re'exp)2/exp(Re

)'(exp)2/exp(Re),(

00

0000

0

00

πωω

πωωωπ

ωωπ

−+=

−−+−=−+=

−+−+=

tEtE

tiEtiEtiiEE

Lc

nttiiEEtL

zy

zyzy

oLzy

ee

eeee

eeE

　　　

　　　

この電場の向きがどのように変化するか、次の表を利用して調べよう。

'tω 0 2/π π 2/3π π2 2/5π π3

)'cos( tω 1 0 －1 0 1 0 －1

)2/'cos( πω −t 0 1 0 －1 0 1 0

また、電場の振幅の2乗は、

1)'(sin)'(cos)2/'(cos)'(cos 2222 =+=−+ tttt ωωπωω

となるため、時間に関係なく一定の大きさの振幅である。右図に示すよ

うに、電場は、その大きさを変えることなく偏光方向が電場の伝播に伴

い回転する。このような光電場を、円偏光という。（ 4/1 波長板の原理）

ωt=0、2π

π/2,5π/2

π

3π/2

ey

ez

83

複屈折

光学的異方性が原因で起こる特徴的な現象に、複屈折(double refraction)がある。これは、無偏光の光が、例え

ば右図のように、空間的に2の光に分離する現象である。この現

象は、結晶が2つの異なる屈折率を持つ場合に可能である。この

現象で大事な所は、境界面（結晶面）に垂直に入射しても、光線

が空間的に2つに分かれることである。

最初は、簡単な場合で複屈折を理解する。紙面に垂直方向

に直線偏光する光電場の屈折率が en 、それ以外の2方向に直

線偏光する光電場の屈折率が on とする。ここで、 oe nn > と

しよう。ここに無偏光の光電場が入射すると、紙面と平行な

偏光方向の電場（屈折率 on ）の屈折角よりも、紙面と垂直な

偏光方向の電場（屈折率 en ）の屈折角のほうが小さく、図の

ように、光電場が結晶内部で 2 つに分かれる。（屈折角の大

きさについては、p.2 の最後の式を参照。）。これが複屈折の

基本現象である。

光電場が結晶に垂直入射する場合の複屈折現象は以下のように、直感的に理解することが出来る。ただし、

以下の説明は、現象を理解する事を優先にしており、正確な説明ではない事を念頭に置こう。注意：先ほどの

例とは、 )(zc 軸の向きが異なる。

平面波である無偏光の光電場の伝播方向は、x 軸に対して垂直で、z軸（c軸）と角度θ で結晶に入射する

とする。

x 軸方向に平行な偏光の向きを持つ光電場（紙面に対して垂直方向）では、光電場の感じる結晶の屈折率は

on である。その波数ベクトルの向きは入射時の伝播方向と同じで、大きさは 0/ cnk oω= である。この場合、

光電場は結晶に垂直に入射した後、結晶中をまっすぐに進む。

一方、 yz 面に平行な偏光方向の光電場（紙面に対して水平方向）では、光電場の感じる結晶の屈折率は、

光電場の偏光方向が y 軸に平行（光電場の伝播方向が z 軸と平行）な時 on であり、光電場の偏光方向が z 軸

に平行（光電場の伝播方向が y 軸と平行）な時 en である。真空中での電磁波の光速度を 0c 、波長を 0λ 、波数

を 0k とすると、 00kc=ω 、 )//(2/2 000 ck ωππλ == の関係式が成り立つ。

入射する光電場の波数ベクトルk を結晶の zy, 軸の2つの方向に分けると、

)cossincossin)()( 000 θθθθ zyzyzzyy kkk e(eeeeekeekk +=+=⋅+⋅=

となる。簡単のため（ほんとうに簡単のためである事を断わっておく。）、屈折率nの結晶に入射後の波数 'k の

無偏光の光

θθθθ

ey

c-軸

無偏光の光

zzyy eekeekk )()( ⋅+⋅=

無偏光の光

入射角

屈折角

θθθθ

ey

eｚ

ex

eｚ

84

各成分が、 00 )//( nknck == ω の関係で変わるとする。すると、結晶に入射後の波数 'k は、次のようになる。

)cossin()cossin'0

00 θθθθ ze

yoozeyn

nnknnk eee(ek +=+=

波数の変化 'kk → は、

)cossin(')cossin0

00 θθθθ ze

yozyn

nnkk eeke(ek +=→+=

となる。ここで、波数 ',kk の zy, 成分の相対的な大きさを比較しよう。 oe nn > の場合、波数k と比較して、

波数 'k の y 成分が z 成分よりも大きくなっている事がわかる。つまり、結晶の伝播に伴い、波数の向きが右

下方に向く事が分かる。

上の説明は波数ベクトルの変化から見た場合であるが、波面の向きの変化からも、結晶伝播に伴う光線の向

きの変化を理解する事ができる。光電場（光線）が結晶内を伝播する際、1 秒間に ye 方向には enc /0 、 ze 方

向には onc /0 だけ波面が進む。もし oe nn > なら、

eo n

c

n

c 00 >

となり、光電場の波面は、右図で ze 方向により進むた

め、光電場の進行方向は右下方向を向く。この様子を

右図に模式的に示した。図に示すように、光の波面が

伝播する方向は、入射方向（水平方向）から斜め下方

向にずれる。（注意：波面の伝播する方向が波数ベクト

ルで表わされるので、波数と波面の一見異なる説明は、同じ説明を、単に言葉を変えて行っているにすぎない。）

以上の説明は、光電場が結晶中を伝播する時の、波面の伝播する方向についての説明であり、光線（光電場

のエネルギー）の伝播方向についての説明ではない。しかし、およそこのような理由で、無偏光の光電場（光

線）が光学的異方性を持つ結晶に入射した場合は、たとえ垂直入射でも、c軸が結晶に入射時の光電場の伝播

方向と異なれば、結晶内部で2つの光電場（光線）に分かれる。この性質を利用したさまざまな偏光ビームス

プリッターが作成されている。

より厳密な理解のためには、上で行ったような電磁波の波面の伝播する向きを示す波数ベクトルkについて

の議論ではなく、電磁波のエネルギーの流れる方向、すなわちポインティング（Poynting）ベクトルの伝播方

向についての理解・議論が必要になる。また、結晶中における光電磁波の伝播方向の波数ベクトルを求めるた

めには、ここでは説明しなかったが、フレネル（Fresnel）の光線方程式を用いる必要がある。

on

c00c

0c

波面の様子

en

c0

真空中結晶中

85

索引

あ行あ行あ行あ行

アンペールの法則 18

アンペール-マクスウェルの法則19

位相 7

位相速度 22、24

一様な空間23

インピーダンス 24

真空の---- 24

円偏光 37

円偏光右回り、左回り 38

か行か行か行か行

回折 11

ガウスの法則 16

可干渉性 10

逆進性 4

干渉 8

換算質量 49

光学的な距離 3

吸収係数 42、43

局所電場 64

虚像 6

屈折率 24、42

群速度 25

クラウジウス-モソッテティの式 70

光学軸 79

コヒーレンス 10

さ行さ行さ行さ行

磁化ベクトル 15

実像 5

周波数7

ジュール熱 30

消衰係数 42

総和則 77

振動子強度 53

た行た行た行た行

直線偏光 37

直線偏光板 37

電荷の保存則 19

電気感受率 15、35、52

電気双極子 15、49

電気伝導度 15

電気双極子が吸収するエネルギー 31

電束密度 15、53

電磁波が横波である事 36

電磁波のエネルギー密度 30

電磁波の吸収 50

透磁率 15、23

等方性媒質 35

な行な行な行な行

は行は行は行は行

反射型回折格子 13

反分極場 69

波長 7

波源 7

反射率 45

半値全幅 52

表皮深さ 63

表面電荷密度 64

Beerの法則 43

フーリエ変換 72

複素感受率テンソル 79

複屈折 83

ファラデーの電磁誘導の法則 17

フェルマーの原理 2

プラズマ 55

プラズマ周波数 55

プラズマ振動 58

分極ベクトル 15

分極による電磁波の

エネルギーの損失 34

分散 54

正常---- 異常---- 54

分散公式 75

複素屈折率 42

平均パワー（電磁波の） 32

変位電流 20

ホイヘンスの原理 9

ポインティングベクトル 30

ま行ま行ま行ま行

マクスウェル方程式 15

や行や行や行や行

ヤングの干渉実験 9

ら行ら行ら行ら行

ラプラス方程式 35

ラプラシアン 23、35

乱反射 4

レンズの公式 5

ローレンツ形 52

ローレンツの局所電場 64

ベクトル演算の定義・公式ベクトル演算の定義・公式ベクトル演算の定義・公式ベクトル演算の定義・公式の説明の説明の説明の説明

),(div txA 15

),(rot txA 15

),(rotrot txA 23

)(div BA× 29

)(divrot xA 21

)( CBA ×× 40

室蘭工業大学応用物理コース応用光学講義ノート...2 幾何光学...

Documents

Transcript of 室蘭工業大学応用物理コース応用光学講義ノート...2 幾何光学...

室蘭工業大学 応用物理コース 応用光学 講義ノート...2 幾何光学...

Documents

Transcript of 室蘭工業大学 応用物理コース 応用光学 講義ノート...2 幾何光学...

室蘭工業大学応用物理コース応用光学講義ノート...2 幾何光学...

Transcript of 室蘭工業大学応用物理コース応用光学講義ノート...2 幾何光学...