tecktom.com

出てくるもの2:Pauli 常磁性と圧縮率

2026年10月15日9分で読めますAI生成教科書『虎の巻』第4部 第1章 自由電子気体と Sommerfeld 展開
目次
  1. Pauli 常磁性
  2. 圧縮率
  3. Fermi 準位の状態密度が全部を決める
  4. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第4部 第1章 §4.1.4 ★★ 磁場でスピンを偏らせる応答も、密度を変える応答も、Fermi 準位の席の数 D(εF)D(\varepsilon _{\mathrm F}) で決まる。帯磁率は温度にほとんどよらず、比熱と磁化率の比(Wilson 比)は 11 になる

Pauli 常磁性

磁束密度 B=Bez\bm{B}=B\bm{e}_z の中では、電子のスピンの向きでエネルギーがずれる。電子の磁気モーメントはスピンと逆向きで大きさ μB\mu_{\mathrm B}(g=2g=2 とする)なので、スピンの zz 成分が上(σ=+1\sigma=+1)の電子は +μBB+\mu_{\mathrm B}B、下(σ=−1\sigma=-1)の電子は −μBB-\mu_{\mathrm B}B だけエネルギーが上がる。ここでは軌道運動への磁場の効果(Landau 準位)を考えず、スピンへの効果だけを入れる。各スピンの状態密度は 12D(ε)\frac12D(\varepsilon ) で、その原点が σμBB\sigma\mu_{\mathrm B}B だけずれるので、スピン σ\sigma の電子の密度は

nσ=∫12D(ε−σμBB)f(ε) dε≈n2−σμBB2∫D′(ε)f(ε) dε=n2−σμBB2∫D(ε)(−∂f∂ε)dεn_\sigma=\int\frac12D(\varepsilon -\sigma\mu_{\mathrm B}B)f(\varepsilon )\,\mathrm{d}\varepsilon \approx\frac{n}{2}-\frac{\sigma\mu_{\mathrm B}B}{2}\int D'(\varepsilon )f(\varepsilon )\,\mathrm{d}\varepsilon =\frac n2-\frac{\sigma\mu_{\mathrm B}B}{2}\int D(\varepsilon )\Bigl(-\frac{\partial f}{\partial \varepsilon }\Bigr)\mathrm{d}\varepsilon

である(最後は部分積分)。磁化は下向きスピンが多い分だけ生じ、M=μB(n↓−n↑)M=\mu_{\mathrm B}(n_{\downarrow}-n_{\uparrow}) である。絶対零度では −∂f/∂ε=δ(ε−εF)-\partial f/\partial\varepsilon =\delta(\varepsilon -\varepsilon _{\mathrm F}) なので

M=μB2D(εF)B,χm=μ0μB2D(εF)(4.1.15)M=\mu_{\mathrm B}^2D(\varepsilon _{\mathrm F})B,\qquad \chi_{\mathrm m}=\mu_0\mu_{\mathrm B}^2D(\varepsilon _{\mathrm F})\tag{4.1.15}

を得る(M=χmHM=\chi_{\mathrm m}H、B=μ0HB=\mu_0H。磁化率は無次元)。これが Pauli 常磁性である。

古典的な(あるいは局在した)独立なスピン NN 個なら、磁化率は Curie 則 χ∝NμB2/(kBT)\chi\propto N\mu_{\mathrm B}^2/( k_{\mathrm B} T) に従い、低温で大きくなる。伝導電子の磁化率がそうならないのは、比熱と同じ理由である。磁場でスピンを反転できるのは Fermi 準位のまわり幅 μBB\mu_{\mathrm B}B の中の電子だけで、深い電子は反転先がすでに埋まっている。数 NN の代わりに D(εF)×(エネルギーの幅)D(\varepsilon _{\mathrm F})\times(\text{エネルギーの幅}) が入り、Curie 則の 1/(kBT)1/( k_{\mathrm B} T) は D(εF)/n∼1/εFD(\varepsilon _{\mathrm F})/n\sim1/\varepsilon _{\mathrm F} に置き換わる。

有限温度の補正は、∫D(−∂f/∂ε) dε\int D(-\partial f/\partial\varepsilon )\,\mathrm{d}\varepsilon に Sommerfeld 展開を当てて求められる。H=D′H=D' に (4.1.9) を当てると D(μ)+π26(kBT)2D′′(μ)D(\mu)+\frac{\pi^2}{6}( k_{\mathrm B} T)^2D''(\mu) で、3次元の自由電子では μ(T)\mu(T) の (4.1.11) と合わせて

χm(T)χm(0)=1−π212(TTF)2+⋯(4.1.16)\frac{\chi_{\mathrm m}(T)}{\chi_{\mathrm m}(0)}=1-\frac{\pi^2}{12}\Bigl(\frac{T}{T_{\mathrm F}}\Bigr)^2+\cdots\tag{4.1.16}

になる。金属の T/TF∼10−3T/T_{\mathrm F}\sim10^{-3} では、Pauli 常磁性は温度によらないと言ってよい。

圧縮率

密度を変える応答も同じ形になる。化学ポテンシャルを dμ\mathrm{d}\mu だけ上げると、絶対零度では Fermi 準位の上に幅 dμ\mathrm{d}\mu の席が開き、dn=D(εF) dμ\mathrm{d} n=D(\varepsilon _{\mathrm F})\,\mathrm{d}\mu だけ電子が増える。熱力学では、NN を固定したときの等温圧縮率 κ=−1V∂V∂P\kappa=-\frac1V\frac{\partial V}{\partial P} が κ=1n2∂n∂μ\kappa=\frac{1}{n^2}\frac{\partial n}{\partial \mu} と書ける(Gibbs–Duhem の関係 dP=n dμ\mathrm{d} P=n\,\mathrm{d}\mu を絶対零度で使う)。したがって

∂n∂μ=D(εF),B=1κ=n2D(εF)=23nεF(4.1.17)\frac{\partial n}{\partial \mu}=D(\varepsilon _{\mathrm F}),\qquad B=\frac1\kappa=\frac{n^2}{D(\varepsilon _{\mathrm F})}=\frac23n\varepsilon _{\mathrm F}\tag{4.1.17}

で、(4.1.6) の BB と一致する。

Fermi 準位の状態密度が全部を決める

比熱 (4.1.13)、帯磁率 (4.1.15)、圧縮率 (4.1.17) の3つは、どれも D(εF)D(\varepsilon _{\mathrm F}) に普遍的な定数を掛けたものである。これは、3つとも「Fermi 準位のまわりの狭い幅に電子を出し入れする」応答であり、幅の中の席の数だけが問題になるからである。そこで、比熱の係数 γV=π23D(εF)kB2\gamma_V=\frac{\pi^2}{3}D(\varepsilon _{\mathrm F}) k_{\mathrm B}^2 と帯磁率の比をとると D(εF)D(\varepsilon _{\mathrm F}) が消え、

RW=π2kB23μ0μB2 χmγV=1(4.1.18)R_{\mathrm W}=\frac{\pi^2 k_{\mathrm B}^2}{3\mu_0\mu_{\mathrm B}^2}\,\frac{\chi_{\mathrm m}}{\gamma_V}=1\tag{4.1.18}

になる。RWR_{\mathrm W} を Wilson 比と呼ぶ。自由電子気体では RW=1R_{\mathrm W}=1 が恒等的に成り立つ。実在の金属で RWR_{\mathrm W} が 11 からずれていれば、それは「比熱で見た状態密度」と「スピンで見た状態密度」が違うということで、電子間の相互作用が残っている証拠になる。相互作用を戻したときに比熱と帯磁率が別々の因子で変わることは、第5部第1章の Fermi 液体で扱う。

例題 4.1.4 銅の Pauli 常磁性と体積弾性率 ★★

オリジナル例題 4.1.1 の銅を自由電子気体とみなす。μ0=1.2566×10−6 N A−2\mu_0=1.2566\times10^{-6}\,\mathrm{N\,A^{-2}}、μB=9.274×10−24 J T−1\mu_{\mathrm B}=9.274\times10^{-24}\,\mathrm{J\,T^{-1}} とする。

  1. (1) D(εF)D(\varepsilon _{\mathrm F}) を eV−1m−3\mathrm{eV^{-1}m^{-3}} と「1原子あたり eV−1\mathrm{eV^{-1}}」で求めよ。
  2. (2) Pauli の磁化率 χm\chi_{\mathrm m} を求めよ。軌道の反磁性(Landau 反磁性)は自由電子では −13-\frac13 倍の大きさで、これを足すといくらか。
  3. (3) B=10 TB=10\,\mathrm T で、全電子のうちスピンが向きを変えた電子の割合 (n↓−n↑)/n(n_\downarrow-n_\uparrow)/n を求めよ。
  4. (4) (4.1.18) の RWR_{\mathrm W} が自由電子で 11 になることを、γV\gamma_V と χm\chi_{\mathrm m} の式を代入して確かめよ。
解答

(1) D(εF)=3n/(2εF)=3×8.47×1028/(2×7.03 eV)=1.81×1028 eV−1m−3D(\varepsilon _{\mathrm F})=3n/(2\varepsilon _{\mathrm F})=3\times8.47\times10^{28}/(2\times7.03\,\mathrm{eV})=1.81\times10^{28}\,\mathrm{eV^{-1}m^{-3}}。nn で割ると1原子あたり 0.213 eV−10.213\,\mathrm{eV^{-1}} である。

(2) D(εF)D(\varepsilon _{\mathrm F}) を J−1m−3\mathrm{J^{-1}m^{-3}} に直して(1.13×10471.13\times10^{47})、χm=μ0μB2D(εF)=1.22×10−5\chi_{\mathrm m}=\mu_0\mu_{\mathrm B}^2D(\varepsilon _{\mathrm F})=1.22\times10^{-5}。Landau 反磁性を足すと 23×1.22×10−5=8.12×10−6\frac23\times1.22\times10^{-5}=8.12\times10^{-6} である。Landau 反磁性の導出は磁場中の軌道運動を量子化する必要があり、第4部第4章で扱う。実在の銅では、内殻の電子の反磁性がこれより大きく、全体としては反磁性を示す。伝導電子の寄与だけを取り出すには、スピンだけに反応する測り方が要る(4.1.5節)。

(3) n↓−n↑=D(εF)μBBn_\downarrow-n_\uparrow=D(\varepsilon _{\mathrm F})\mu_{\mathrm B}B なので、割合は D(εF)μBB/n=32μBB/εF=1.24×10−4D(\varepsilon _{\mathrm F})\mu_{\mathrm B}B/n=\frac32\mu_{\mathrm B}B/\varepsilon _{\mathrm F}=1.24\times10^{-4} である。10 T10\,\mathrm T という強い磁場でも、向きを変えるのは1万個に1個程度である。Zeeman エネルギー μBB=0.58 meV\mu_{\mathrm B}B=0.58\,\mathrm{meV} が εF\varepsilon _{\mathrm F} に比べて小さいからである。

(4) χm/γV=μ0μB2D(εF)/(π23D(εF)kB2)=3μ0μB2/(π2kB2)\chi_{\mathrm m}/\gamma_V=\mu_0\mu_{\mathrm B}^2D(\varepsilon _{\mathrm F})/\bigl(\frac{\pi^2}{3}D(\varepsilon _{\mathrm F}) k_{\mathrm B}^2\bigr)=3\mu_0\mu_{\mathrm B}^2/(\pi^2 k_{\mathrm B}^2) なので、前の因子と掛けて RW=1R_{\mathrm W}=1。D(εF)D(\varepsilon _{\mathrm F}) が分子と分母で消えるので、nn や有効質量によらず 11 である。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ