ドルーデモデル (英 : Drude model )またはドルーデ模型 は、1900年 にパウル・ドルーデ により提唱された、電気伝導 についてのモデル で、物質(特に金属)内部の電子 の特性について記述する。このモデルは気体分子運動論 を応用しており、固体中の電子の微視的挙動は古典的に扱えるものとし、重く動きづらい陽イオン の間をピンボール のように電子が常に行き来しながら満たしているという仮定をおく。
ドルーデモデルにおける電子(青で示す)はより重い静的な結晶イオン(赤で示す)の間で常に反跳している。
ドルーデモデルから導かれる最も重要な結論は、電子の運動方程式
d
d
t
p
(
t
)
=
q
(
E
+
p
(
t
)
×
B
m
)
−
p
(
t
)
τ
{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\boldsymbol {p}}(t)=q\left({\boldsymbol {E}}+{\frac {{\boldsymbol {p}}(t)\times {\boldsymbol {B}}}{m}}\right)-{\frac {{\boldsymbol {p}}(t)}{\tau }}}
と、電流密度 J と電場 E との間の線形な関係式
J
=
(
n
q
2
τ
m
)
E
{\displaystyle {\boldsymbol {J}}=\left({\frac {nq^{2}\tau }{m}}\right){\boldsymbol {E}}}
の2つである。ここで t は時間、 p , q, n, m, τ はそれぞれ電子の運動量 、電荷 、数密度 、質量 、陽イオンとの衝突の間の平均自由時間 [注 1] を示す。後者の式は、電磁気学 において最も普遍的な関係式の1つであるオームの法則 が、何故成立するのかを半定量的に説明することができる点で特に重要である。
このモデルは1905年 にローレンツ により拡張された[注 2] 、古典的 なモデルである。後の1933年 に、ゾンマーフェルト とベーテ により量子論 の結果が取り込まれ、ドルーデ・ゾンマーフェルトモデル へと発展した。
ドルーデモデルでは、金属が正に帯電したイオンの集まりと、それから放出された膨大な数の「自由電子 」から構成されていると考える。このことは、原子の価電子準位が他の原子によるポテンシャルと接触することによって非局在化していると考えることもできる。
ドルーデモデルでは電子とイオン、もしくは電子同士の間に働く一切の長距離相互作用は無視される。自由電子が環境との間に持つ唯一の相互作用は衝突の一瞬のうちにのみ行なわれる。自由電子が次に衝突するまでの平均時間は τ であり、衝突する相手の性質はドルーデモデルの計算や結果には影響しない。
直流電場
ドルーデモデルによる最も単純な解析では、電場 E が一様かつ静的に印加されており、電子の熱運動速度 が十分に高く、無限小の運動量 dp が τ 秒ごとにくりかえされる衝突の間に蓄積していくものと仮定する。
このとき、時刻 t における孤立電子は最後に衝突してから平均して時間 τ だけ経過しており、従って蓄積された運動量は以下のように表わされる。
Δ
⟨
p
⟩
=
q
E
τ
{\displaystyle \Delta \langle {\boldsymbol {p}}\rangle =q{\boldsymbol {E}}\tau }
最後の衝突の際に、この電子が前向きに反跳した確率と、後ろ向きに反跳した確率とは等しいので、衝突以前の電子の運動量の寄与は無視できるものと考えられるので、電子の運動量は次式で表わされる。
⟨
p
⟩
=
q
E
τ
{\displaystyle \langle {\boldsymbol {p}}\rangle =q{\boldsymbol {E}}\tau }
この式に以下の二つの式を代入すると、前述したオームの法則が得られる。
⟨
p
⟩
=
m
⟨
v
⟩
{\displaystyle \langle {\boldsymbol {p}}\rangle =m\langle {\boldsymbol {v}}\rangle }
J
=
n
q
⟨
v
⟩
{\displaystyle {\boldsymbol {J}}=nq\langle {\boldsymbol {v}}\rangle }
J
=
(
n
q
2
τ
m
)
E
{\displaystyle {\boldsymbol {J}}=\left({\frac {nq^{2}\tau }{m}}\right){\boldsymbol {E}}}
時間変動の解析
実効抗力 を導入することによって、上と同じ特性を説明することもできる。時刻 t = t 0 + dt における電子の平均運動量は以下のように表わせる。
⟨
p
(
t
0
+
d
t
)
⟩
=
(
1
−
d
t
τ
)
(
⟨
p
(
t
0
)
⟩
+
q
E
d
t
)
{\displaystyle \langle {\boldsymbol {p}}(t_{0}+\mathrm {d} t)\rangle =\left(1-{\frac {\mathrm {d} t}{\tau }}\right)\left(\langle {\boldsymbol {p}}(t_{0})\rangle +q{\boldsymbol {E}}\mathrm {d} t\right)}
なぜなら、平均すれば 1 − dt / τ だけの電子はまだ衝突していないはずであり、既に衝突した電子は総運動量に無視できるオーダーの寄与しかもたないからである。
代数的な処理を施して dt 2 のオーダーの項を無視すると、以下の微分方程式 が結果として得られる。
d
d
t
⟨
p
(
t
)
⟩
=
q
E
−
⟨
p
(
t
)
⟩
τ
{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\langle {\boldsymbol {p}}(t)\rangle =q{\boldsymbol {E}}-{\frac {\langle {\boldsymbol {p}}(t)\rangle }{\tau }}}
ここで ⟨ p ⟩ は平均運動量を示す。この線形非斉次微分方程式は以下のような一般解を持つ。
⟨
p
(
t
)
⟩
=
q
τ
E
(
1
−
e
−
t
/
τ
)
+
⟨
p
(
0
)
⟩
e
−
t
/
τ
{\displaystyle \langle {\boldsymbol {p}}(t)\rangle =q\tau {\boldsymbol {E}}(1-e^{-t/\tau })+\langle {\boldsymbol {p(0)}}\rangle e^{-t/\tau }}
よって、定常解 (d⟨ p ⟩ / dt = 0 ) は
⟨
p
⟩
=
q
τ
E
{\displaystyle \langle {\boldsymbol {p}}\rangle =q\tau {\boldsymbol {E}}}
上述のとおり、平均運動量は平均速度と関連しており、それを通じて電流密度と関連づけることができる。
⟨
p
⟩
=
m
⟨
v
⟩
{\displaystyle \langle {\boldsymbol {p}}\rangle =m\langle {\boldsymbol {v}}\rangle }
J
=
n
q
⟨
v
⟩
{\displaystyle {\boldsymbol {J}}=nq\langle {\boldsymbol {v}}\rangle }
ここから、直流 電気伝導率 σ 0 でオームの法則を満たす物質は以下を満たすことが示せる。
J
=
(
n
q
2
τ
m
)
E
{\displaystyle {\boldsymbol {J}}=\left({\frac {nq^{2}\tau }{m}}\right){\boldsymbol {E}}}
τ = 10− 5 、σ 0 = 1 と置いたときの複素電気伝導率の周波数依存性。
ドルーデモデルにより、角周波数 ω で時間変動する電場への応答を予測することもできる。
σ
(
ω
)
=
σ
0
1
+
i
ω
τ
=
σ
0
1
+
ω
2
τ
2
+
i
ω
τ
σ
0
1
+
ω
2
τ
2
{\displaystyle \sigma (\omega )={\frac {\sigma _{0}}{1+i\omega \tau }}={\frac {\sigma _{0}}{1+\omega ^{2}\tau ^{2}}}+i\omega \tau {\frac {\sigma _{0}}{1+\omega ^{2}\tau ^{2}}}}
ここで、以下の二つを仮定している。
E
(
t
)
=
ℜ
(
E
0
e
i
ω
t
)
;
{\displaystyle E(t)=\Re (E_{0}e^{i\omega t});}
J
(
t
)
=
ℜ
(
σ
(
ω
)
E
0
e
i
ω
t
)
{\displaystyle J(t)=\Re (\sigma (\omega )E_{0}e^{i\omega t})}
i を全て − i で置き換えた表式を用いることもある。虚部は電子が変動する電場に追随して加速するまでにおよそ τ だけの時間を要することに起因する電流の電場に対する遅れを表わす。ここまで、電子についてドルーデモデルを適用してきたが、このモデルは電子と正孔 (半導体中の正に帯電した電荷担体 )のどちらにも適用することができる。 σ (ω ) の曲線をグラフに示す。
実際の物質のドルーデ応答
ドルーデ金属のふるまいの時間的または周波数的特徴、つまり時定数 τ での減衰または上記のような σ (ω ) の周波数依存性は、ドルーデ応答と呼ばれる。 典型的で単純な実在の金属(例:室温におけるナトリウム、銀、金)においては、特性周波数 τ − 1 が、ドルーデモデルでは無視されている物性(バンド構造 など)が重要な働きをする赤外領域に入っているため、ドルーデ応答を実験的に観測することはできない。しかし、他の特定の金属性物質ではドルーデモデルの予測と非常によく一致する周波数依存伝導率 σ (ω ) を示すこともある。このような物質では、緩和速度 τ − 1 が非常に低い周波数領域にある。このような物質の例として、ドープ された単結晶半導体や、高移動度二次元電子ガス 、重い電子系 を持つ金属があげられる。
歴史的には、ドルーデの式は電荷担体を理想気体 とみなす誤った仮定から導かれた。 現在では、電荷担体はフェルミ分布 に従うこと、また無視できない相互作用をもつことが知られているが、にもかかわらずドルーデモデルが有効である。1957年 にランダウ が相互作用をもつ粒子の気体は、ほとんど相互作用をしない準粒子 の系によって記述できることを示し、このために金属中の伝導電子にドルーデモデルが有効であることがわかった。
単純で古典的なドルーデモデルは室温下の金属の直流・交流 電気伝導率、ホール効果 、熱伝導率 (の電子による寄与)をよく説明する。 このモデルにより、1853年に発見されたヴィーデマン=フランツ則 も説明される。しかし、このモデルでは金属の電子による熱容量を大きく過大評価してしまう。実際には、金属と不導体は室温ではほとんど同じ熱容量を示す。このモデルが正孔にも有効であることは、ホール効果によって立証された。
この理論にまつわるトリビアとして、ドルーデは彼の原論文で根本的な間違いを犯しており、電気伝導度を古典的な正しい値の半分としてしまっていたことがあげられる。
注釈
つまり、電子が直前に陽イオンと衝突してから平均してどのくらい経っているかであり、衝突と衝突の時間間隔の平均ではない。
したがって、このモデルはドルーデ・ローレンツモデル としても知られている。
原論文
Drude, Paul (February 22, 1900). “Zur Elektronentheorie der metalle” (PDF ). Annalen der Physik 306 (3): 566-613. Bibcode : 1900AnP...306..566D . doi :10.1002/andp.19003060312 . ISSN 0003-3804 . LCCN 50-13519 . OCLC 5854993 . http://onlinelibrary.wiley.com/doi/10.1002/andp.19003060312/epdf .
Drude, Paul (September 9, 1900). “Zur Elektronentheorie der Metalle; II. Teil. Galvanomagnetische und thermomagnetische Effecte” (PDF ). Annalen der Physik 308 (11): 369-402. Bibcode : 1900AnP...308..369D . doi :10.1002/andp.19003081102 . ISSN 0003-3804 . LCCN 50-13519 . OCLC 5854993 . http://onlinelibrary.wiley.com/doi/10.1002/andp.19003081102/pdf .
Dressel, M.; Scheffler, M. (October 10, 2005). “Verifying the Drude response” (PDF ). Ann. Phys. 15 (7–8): 535–544. Bibcode : 2006AnP...518..535D . doi :10.1002/andp.200510198 . ISSN 0003-3804 . LCCN 50-13519 . OCLC 5854993 . http://onlinelibrary.wiley.com/doi/10.1002/andp.200510198/pdf .
van Exter, M.; Grischkowsky, D. (December 22, 1989). “Carrier dynamics of electrons and holes in moderately doped silicon” (PDF ). Phys. Rev. B 41 (17): 12140–12149. Bibcode : 1990PhRvB..4112140V . doi :10.1103/PhysRevB.41.12140 . ISSN 1098-0121 . OCLC 643849502 . http://journals.aps.org/prb/pdf/10.1103/PhysRevB.41.12140 .
Scheffler, M.; Dressel, M.; Jourdan, M.; Adrian, H. (July 25, 2005). “Extremely slow Drude relaxation of correlated electrons”. Nature 438 (7071): 1135–1137. Bibcode : 2005Natur.438.1135S . doi :10.1038/nature04232 . ISSN 0028-0836 . OCLC 263593080 .
Burke, P. J.; Spielman, I. B.; Eisenstein, J. P.; Pfeiffer, L. N.; West, K. W. (June 4, 1999). “High frequency conductivity of the high-mobility two-dimensional electron gas” (PDF ). Appl. Phys. Lett. 76 (6): 745–747. Bibcode : 2000ApPhL..76..745B . doi :10.1063/1.125881 . ISSN 0003-6951 . OCLC 1580952 . http://authors.library.caltech.edu/2073/1/BURapl00b.pdf .
書籍
洋書
訳書
Solid State Physics
Electricity and Magnetism