「電荷密度」の版間の差分

2009年5月12日 (火) 21:36時点における版

2025年7月20日 (日) 03:38時点における最新版

(24人の利用者による、間の28版が非表示)

1行目:

{{出典の明記|date=2015年7月}}

'''電荷密度'''（でんかみつど、charge density）は、単位[[体積]]当たりの[[電荷]]の分布量（体積[[密度]]）。電荷を担うものとしては[[電子]]や[[原子核]]、[[イオン]]のような粒子（素粒子や[[正孔]]などを含む）であったり、仮想的に一様に分布する電荷のような場合（→参照：[[ジェリウムモデル]]）もある。

{{物理量

| 名称 =

| 英語 = charge density

| 画像 =

| 記号 =''ρ''

| 次元 = [[時間|T]] [[電流|I]] [[長さ|L]]{{sup-|3}}

| 階 =[[スカラー]]

| SI =クーロン毎立方メートル

| CGS =クーロン毎立方センチメートル

| MTS =

| FPS =

| MKSG =

| CGSG =

| FPSG =

| プランク =

| 原子 =

}}

'''電荷密度'''（でんかみつど、{{lang-en-short|charge density}}）とは、ある特定の空間領域において、単位体積あたりに存在する電荷の量を示す物理量である。この概念を理解する上で極めて重要なのは、どの大きさのスケール（体積）で空間を定義するかによって、電荷密度の値が全く異なる意味を持つという点である。

物質は、正の電荷を持つ原子核と、負の電荷を持つ電子から構成される原子（または分子）の集合体である。この基本構造に基づき、電荷密度を異なるスケールで考察する。

実験的にはX線回折実験による構造解析から得られた結果を最大エントロピー法などを使って実空間での電子の電荷分布（→電子密度に相当）が求まる。また中性子回折実験の結果から同様な手法により原子核の密度が求まる（？）。

'''ミクロスケール'''（原子内部のスケール）

原子1個よりもはるかに小さい、極めて微小な体積で空間を定義する場合、電荷密度は場所によって著しく変動する。

*原子核の位置では、正電荷が密集しているため、電荷密度は非常に大きな正の値をとる。

*その周囲を運動する電子が存在する領域では、負電荷が存在するため、電荷密度は負の値となる。

*原子核と電子の間の真空の空間では、電荷密度はゼロである。

このように、ミクロスケールにおける電荷密度は、極めて不均一に分布している。

'''マクロスケール'''（原子サイズ以上のスケール）

一方、原子数個分、あるいはそれ以上の大きさの体積で空間を定義する場合、状況は本質的に異なる。

*通常、物質は電気的に中性であり、定義された体積内に含まれる原子核の正電荷の総量と、電子の負電荷の総量は互いに相殺される。

*その結果、このスケールにおける電荷密度は、空間的に平均化され、実質的にゼロと見なされる。

結論として、「電荷密度」という物理量を扱う際は、それが原子内部のようなミクロな電荷分布を論じているのか、物質全体としてのマクロな電気的性質を論じているのか、その定義スケールを明確に意識することが不可欠である。

（注：原子核が持つ正電荷の根源は陽子にあるが、陽子自体もクォークなどの素粒子から構成される。電荷の起源をどのレベルで議論するかは、物理学で扱うエネルギーや空間のスケールに依存する。）

物性物理学の文脈では、物質全体の電荷密度（通常はゼロ）とは別に、特定の役割を担う電荷の密度だけを意図的に取り出して議論することが頻繁にある。その代表例が、金属中の伝導電子である。

'''特定の担い手（キャリア）に注目した電荷密度'''

金属の銅（Cu）を例にとると、これは規則正しく並んだ「銅イオン（Cu+）」と、その周りを自由に動き回る「伝導電子（e-）」の集合体としてモデル化される。この場合、電荷密度は以下のように区別して扱われる。

*銅イオンによる電荷密度: 各格子点に局在した、正の電荷密度。

*伝導電子による電荷密度: 金属全体に広がった、負の電荷密度。

このように、電流や熱輸送といった現象を議論する際には、物質の性質を支配する「伝導電子の電荷密度」という言葉が、文脈上、単に「電荷密度」として用いられることがある。これは、静的なイオンを背景として捉え、動的な電子の振る舞いに焦点を当てているためである。

'''電荷密度の実験的な測定手法'''

理論モデルだけでなく、電荷密度の空間的な分布は実験によって直接観測することも可能である。この測定には、粒子線と物質との相互作用の違いが利用される。

[[X線回折]]

*X線は、原子内の電子によって強く散乱される性質を持つ。そのため、結晶にX線を照射し、その回折パターンを解析することで、結晶内における電子の空間分布、すなわち負の電荷密度の分布を精密に決定することができる。

[[中性子回折法]]

*一方、電気的に中性である中性子線は、電子とはほとんど相互作用せず、主に原子核によって散乱される。したがって、中性子線回服実験を行えば、X線回折とは相補的に、結晶内における原子核の位置、すなわち正の電荷密度の中心位置を正確に特定することが可能となる。

このように、X線と中性子線という二つの異なるプローブ（探針）を使い分けることで、物質内の正負の電荷密度分布をそれぞれ分離して可視化することができる。

==バンド計算での電荷密度==

[[バンド計算]]では通常、電荷密度とは電子の密度のことを示す。従って、この場合は'''電子密度ｍｎ'''（electron density）とも言う。電子以外の電荷（例えば[[イオン (化学)|イオン]]など）に対しても "電荷密度" の表記を用いることがある；ｌｋので注意が必要。　

ｋｊｈ[[バンド計算]]では、実空間での電荷密度 ''ρ''('''''r''''') は[[波動関数]ｆｇｊｆ] ''ψ''<SUB>''i'','''''k'''''</SUB>('''''r''''') のノルムを取ることにより求められる：

:<math>\rho (\boldsymbol{r}) = \sum_{i,\boldsymbol{k}} f_{i,\boldsymbol{k}} |\psi_{i,\boldsymbol{k}} (\boldsymbol{r})|^2. </math>

i, '''k''' はそれぞれ[[バンド構造|バンド]]と[[k点]]の指標。''f''<S。ｋｊUB>i,'''ｖｎk'''</SUB> は、各 k 点上の各バンドでの電子の占有数。尚、バンド計算では普通[[原子単位]]を用いるので素電荷は、''e''<SUP>2</SUP> = 1（ハｚートリー原子ｋｂ単位の場合）としている。ここで占有数は、''N'' を系の全電子数とすると

''i'', '''''k''''' はそれぞれ[[バンド構造|バンド]]と[[ブリュアンゾーン#k点|k点]]の指標。''f''<SUB>''i'','''''k'''''</SUB> は、各 ''k'' 点上の各バンドでの電子の占有数。なお、バンド計算では普通[[原子単位]]を用いるので素電荷は、''e'' = 1（[[ハートリー原子単位|ハートリー原子単位系]]の場合）としている。ここで占有数は、''N'' を系の全電子数とすると

ｈｊｈ

となる。バンド計算において波動関数は規格化されており、占有数 はロード非整数となる場合がある。

:<math> \sum_{i,\boldsymbol{k}} f_{i,\boldsymbol{k}} = N </math>

となる。バンド計算において波動関数は規格化されており、占有数 ''f''<SUB>''i'','''''k'''''</SUB> は非整数となる場合がある。

（''i'' は[[虚数単位]]）であり、上式左辺の ρ('''G''') は構造因子とガトリング言われるが、このことを逆空間表示での電荷密度と言う場合もある。

実空間の電荷密度を[[フーリエ変換]]したものは、

:<math> \rho (\boldsymbol{G}) = \frac{1}{V} \int \rho (\boldsymbol{r}) e^{-i\boldsymbol{G}\cdot\boldsymbol{r}} d\boldsymbol{r}

</math>

（''i'' は[[虚数単位]]）であり、上式左辺の ''ρ''('''''G''''') は構造因子と言われるが、このことを逆空間表示での電荷密度と言う場合もある。

==運動量密度==

実空間の波動関数をフーリエ変換して（指標 ''i'', '''''k''''' は省略、''V''：系の体積）、

:<math> \psi(\boldsymbol{G}) = \frac{1}{V} \int \psi(\boldsymbol{r}) e^{-i\boldsymbol{G}\cdot\boldsymbol{r}} d\boldsymbol{r}

:<math> \psi(\mathbf{G}) =

{1 \over V} \int

\psi(\mathbf{r}) e^{-i\mathbf{G}\cdot\mathbf{r}}

d\mathbf{r}

</math>

を得る。''ψ''('''''G''''') は[[逆格子空間]]（[[運動量空間]]）での波動関数であり、これの[[ノルム]]をとると、

:<math> P(\mathbfboldsymbol{G}) = |\psi(\mathbfboldsymbol{G})|^2 </math>

となり、上式左辺の ''P''('''''G''''') は逆格子空間での電荷密度と言えるが、通常は'''運動量密度'''(momentum density)と呼ばれる。

運動量密度は、[[コンプトン散乱]]や電子‐[[陽電子]]消滅実験などの実験によって観測される量で、対象が[[金属]]（含む[[半金属]]）の場合、[[フェルミ面]]の情報を含んでいる。

[[自由電子]]の場合の運動量密度 ''ρ''('''''P''''') は、自由電子の[[実空間（3]]（[[3次元]]）での波動関数 ''ψ'' が[[平面波]] <math> e^{-i\boldsymbol{k} \cdot \boldsymbol{r}}</math>であるから、

:<math> e^{-i\mathbf{k} \cdot \mathbf{r}}</math>

であるから、

:<math>\iint e^{-i\boldsymbol{G} \cdot (\boldsymbol{r}-\boldsymbol{r}')}

:<math> \int

e^{-i\mathbfboldsymbol{Gk} \cdot (\mathbfboldsymbol{r}-\mathbfboldsymbol{r}') } \cdot

d\boldsymbol{r}d\boldsymbol{r'}

e^{i\mathbf{k} \cdot (\mathbf{r}-\mathbf{r}') }

= \iint e^{-i(\boldsymbol{G}-\boldsymbol{k})\cdot(\boldsymbol{r}-\boldsymbol{r}')} d\boldsymbol{r} d\boldsymbol{r}'

d\mathbf{r}d\mathbf{r'} =

= \delta (\boldsymbol{G}-\boldsymbol{k})

\int

= \rho_{\boldsymbol{k}}(\boldsymbol{G})</math>

e^{-i(\mathbf{G}-\mathbf{k})\cdot(\mathbf{r}-\mathbf{r}')}

d\mathbf{r}d\mathbf{r}'</math>

:<math>= \delta (\mathbf{G}-\mathbf{k})

= \rho_{\mathbf{k}}(\mathbf{G})

</math>

となり（[[体積]]は省略）、

:<math> \rho(\mathbfboldsymbol{G}) =

= \sum_{\mathbfboldsymbol{k}} \rho_{\mathbfboldsymbol{k}} (\mathbfboldsymbol{G}) =

= \sum_{|\mathbfboldsymbol{k}| \le k_F} f_{\boldsymbol{k}} \rho_{\boldsymbol{k}} (\boldsymbol{G})

f_{\mathbf{k}} \rho_{\mathbf{k}} (\mathbf{G})

</math>

を得る（''f''<SUB>'''''k'''''</SUB> は[[フェルミ分布関数]]←電荷密度での占有数と表記が類似するが異なるものである）。実際は、[[2次元]]ないし [[1次元]]表示したものが[[実験]]による観測結果と比較される。

;2次元表示

:<math> \rho(G_y,G_z)

= \int \rho (\mathbfboldsymbol{G}) dG_x

= \int_{G_x^2 \le G_F^2 - G_y^2 - G_z^2} dG_x

</math>

= 2 \sqrt{G_F^2 - G_y^2 - G_z^2}.

:<math>

= \int_{G_x^2 \le G_F^2 - G_y^2 - G_z^2} dG_x

= 2 \sqrt{G_F^2 - G_y^2 - G_z^2}.

</math>

;1次元表示

:<math> \rho(G_z)

= \iint \rho (\mathbfboldsymbol{G}) d G_x dG_y

= \iint_{G_x^2 + G_y^2 + G_z^2 \le G_F^2} dG_x dG_y

</math>

= G_F^2 - G_z^2.

:<math>

= \iint_{G_x^2 + G_y^2 + G_z^2 \le G_F^2} dG_x dG_y

= G_F^2 - G_z^2.

</math>

以上から、3次元での[[自由電子]]の運動量密度の2次元表示は[[半球]]状、1次元表示は[[放物線]]となる。実際に観測されるものは、[[アルカリ金属]]のような[[価電子]]が自由電子的であるような場合を除いて自由電子のものとは大分異なった形状になることが多い。

==関連項目==

* [[物性物理学]]

* [[第一原理バンド計算]]

** [[電荷密度の混合の仕方]]

* [[分極電荷密度 ]]

* [[電荷・電流密度]]

[[Category{{DEFAULTSORT:固体物理学|てんかみつと]]}}

[[Category:密度]]

[[Category:固体物理学]]

[[ca:Densitat de càrrega]]

[[de:Ladungsdichte]]

[[en:Charge density]]

[[es:Carga eléctrica#Densidad de carga eléctrica]]

[[fa:چگالی بار الکتریکی]]

[[fi:Varaustiheys]]

[[fr:Densité de charge]]

[[he:צפיפות מטען]]

[[mn:Цэнэгийн нягт]]

[[nl:Ladingsdichtheid]]

[[pl:Gęstość ładunku]]

[[pt:Densidade de carga]]

[[ru:Плотность заряда]]

[[sl:Gostota naboja]]

[[uk:Густина заряду]]