跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2014年9月実施 II-3

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

物質中の電子(質量 mm)の状態を考えよう。電子スピンは、zz 軸を量子化軸とし sz=±1/2s_z=\pm1/2 の2方向しか取れない状態にある。この電子の集まり(“電子系”とよぶ)をスピンを持つ自由粒子の系と考える。以下では電子の電荷の効果は無視する。この電子系は NN 個の電子からなり、体積 VV の容器に入っているとする。電子の1粒子状態は運動量空間に一様な密度 V/(2πℏ)3V/(2\pi\hbar)^3 で分布している。スピン自由度も考慮すると、密度は 2V/(2πℏ)32V/(2\pi\hbar)^3 となる。絶対零度では、電子は3次元運動量空間において原点を中心に描いた半径 pFp_F の球の内部の状態を満たし、pFp_F の大きさは内部の状態数が全電子数 NN に等しいという条件から求めることができる。pFp_F をフェルミ運動量、ϵF=pF2/(2m)\epsilon_F=p_F^2/(2m) をフェルミエネルギーとよぶ。

(1) フェルミ運動量 pFp_F とフェルミエネルギー ϵF\epsilon_F を N,V,m,ℏN,V,m,\hbar を用いて求めよ。また D(ϵF)=dN/dϵFD(\epsilon_F)=dN/d\epsilon_F で定義される状態密度 D(ϵF)D(\epsilon_F) を求めよ。

(2) 絶対零度における電子系の全運動エネルギー U0U_0 は N,ϵFN,\epsilon_F を用いて

U0=∫0ϵFϵD(ϵ) dϵ=a0NϵFU_0=\int_0^{\epsilon_F}\epsilon D(\epsilon)\,d\epsilon=a_0N\epsilon_F

と表される。a0a_0 の値を求めよ。

(3) フェルミ気体とみなせる電子系の圧力は P=−(dU0/dV)P=-(dU_0/dV) で定義されるとし、電子系の PVPV を U0U_0 を用いて表せ。ただし、全粒子数 NN は一定とする。

次に、フェルミ縮退をしている NN 電子系の絶対零度における磁化 MM を以下の手順で求めよう。電子スピン sz=±1/2s_z=\pm1/2 は磁気モーメント ∓μB\mp\mu_B を持ち、外部磁場 B  (>0)B\;(>0) の下でのゼーマンエネルギーは ±μBB\pm\mu_BB である。磁場 BB における sz=1/2s_z=1/2 の電子、sz=−1/2s_z=-1/2 の電子の数をそれぞれ、N+=N(1−α)/2, N−=N(1+α)/2N^+=N(1-\alpha)/2,\ N^-=N(1+\alpha)/2 とする(α>0\alpha>0)。

(4) 有限磁場での sz=1/2s_z=1/2 の電子系が持つ運動エネルギー E+E^+ を、磁場 B=0B=0 のときのフェルミエネルギー ϵF\epsilon_F、および N,αN,\alpha を用いて表せ。同様に、sz=−1/2s_z=-1/2 の電子系が持つ運動エネルギー E−E^- を求めよ。

(5) 電子系の全ゼーマンエネルギー EZE_Z を N,B,μB,αN,B,\mu_B,\alpha を用いて表せ。

(6) 電子系の全エネルギーを Etot=E++E−+EZE_{\mathrm{tot}}=E^++E^-+E_Z とする。μBB≪ϵF\mu_BB\ll\epsilon_F のとき EtotE_{\mathrm{tot}} は α≪1\alpha\ll1 で最小値をとる。EtotE_{\mathrm{tot}} を最小にする α\alpha を求め、これに対応する電子系の磁化 MM を ϵF,N,B,μB\epsilon_F,N,B,\mu_B を用いて表せ。

さらに、この NN 電子系の電子間に次のような相互作用が働くとする。すなわちスピンが平行な電子間の相互作用エネルギーを −J  (J≥0)-J\;(J\ge0) とし、スピンが反平行な電子間には相互作用は働かないとする。

(7) N≫1N\gg1 の条件下で、sz=1/2s_z=1/2 の電子間に働く全相互作用エネルギーはどのように表されるか。J,N,αJ,N,\alpha を用いて表示せよ。

(8) 電子系の全エネルギー EtotE_{\mathrm{tot}} に、スピン間に働く相互作用も考慮した全エネルギーを E~tot\widetilde E_{\mathrm{tot}} とする。μBB≪ϵF, NJ≪ϵF\mu_BB\ll\epsilon_F,\ NJ\ll\epsilon_F のとき E~tot\widetilde E_{\mathrm{tot}} は α≪1\alpha\ll1 で最小値をとる。E~tot\widetilde E_{\mathrm{tot}} を最小にする α\alpha を求め、対応する磁化 MM を ϵF,N,B,μB,J\epsilon_F,N,B,\mu_B,J を用いて表せ。

题目描述​

考虑物质中电子(质量为 mm)的状态。以 zz 轴为量子化轴,电子自旋只能取 sz=±1/2s_z=\pm1/2 两个方向。将这一电子集合(称为“电子系”)视为具有自旋的自由粒子系统。以下忽略电子电荷的作用。系统由 NN 个电子组成,置于体积为 VV 的容器内。电子的单粒子态在动量空间以均匀密度 V/(2πℏ)3V/(2\pi\hbar)^3 分布;计入自旋自由度后,密度为 2V/(2πℏ)32V/(2\pi\hbar)^3。在绝对零度下,电子填满三维动量空间中以原点为球心、半径为 pFp_F 的球内部的状态。由球内状态数等于总电子数 NN 的条件可求出 pFp_F。称 pFp_F 为费米动量,ϵF=pF2/(2m)\epsilon_F=p_F^2/(2m) 为费米能量。

(1) 用 N,V,m,ℏN,V,m,\hbar 求费米动量 pFp_F 和费米能量 ϵF\epsilon_F,并求由 D(ϵF)=dN/dϵFD(\epsilon_F)=dN/d\epsilon_F 定义的态密度 D(ϵF)D(\epsilon_F)。

(2) 绝对零度下,电子系的总动能 U0U_0 可用 N,ϵFN,\epsilon_F 表示为

U0=∫0ϵFϵD(ϵ) dϵ=a0NϵF.U_0=\int_0^{\epsilon_F}\epsilon D(\epsilon)\,d\epsilon=a_0N\epsilon_F.

求 a0a_0 的值。

(3) 将电子系视为费米气体,并将压强定义为 P=−(dU0/dV)P=-(dU_0/dV),用 U0U_0 表示电子系的 PVPV。总粒子数 NN 保持不变。

接着按以下步骤求费米简并的 NN 电子系在绝对零度下的磁化 MM。电子自旋 sz=±1/2s_z=\pm1/2 具有磁矩 ∓μB\mp\mu_B,在外磁场 B  (>0)B\;(>0) 中的塞曼能为 ±μBB\pm\mu_BB。磁场 BB 下,自旋 sz=1/2s_z=1/2 和 sz=−1/2s_z=-1/2 的电子数分别为 N+=N(1−α)/2, N−=N(1+α)/2N^+=N(1-\alpha)/2,\ N^-=N(1+\alpha)/2,其中 α>0\alpha>0。

(4) 用零磁场 B=0B=0 时的费米能量 ϵF\epsilon_F 及 N,αN,\alpha,表示有限磁场中 sz=1/2s_z=1/2 电子系的动能 E+E^+。同样求出 sz=−1/2s_z=-1/2 电子系的动能 E−E^-。

(5) 用 N,B,μB,αN,B,\mu_B,\alpha 表示电子系的总塞曼能 EZE_Z。

(6) 设电子系的总能量为 Etot=E++E−+EZE_{\mathrm{tot}}=E^++E^-+E_Z。当 μBB≪ϵF\mu_BB\ll\epsilon_F 时,EtotE_{\mathrm{tot}} 在 α≪1\alpha\ll1 处取得最小值。求使 EtotE_{\mathrm{tot}} 最小的 α\alpha,并用 ϵF,N,B,μB\epsilon_F,N,B,\mu_B 表示相应的磁化 MM。

进一步设该 NN 电子系的电子间存在如下相互作用:自旋平行的电子间相互作用能为 −J  (J≥0)-J\;(J\ge0),自旋反平行的电子间没有相互作用。

(7) 在 N≫1N\gg1 的条件下,sz=1/2s_z=1/2 电子之间的总相互作用能如何表示?用 J,N,αJ,N,\alpha 写出表达式。

(8) 在电子系总能量 EtotE_{\mathrm{tot}} 的基础上,再计入自旋间相互作用,记所得总能量为 E~tot\widetilde E_{\mathrm{tot}}。当 μBB≪ϵF, NJ≪ϵF\mu_BB\ll\epsilon_F,\ NJ\ll\epsilon_F 时,E~tot\widetilde E_{\mathrm{tot}} 在 α≪1\alpha\ll1 处取得最小值。求使 E~tot\widetilde E_{\mathrm{tot}} 最小的 α\alpha,并用 ϵF,N,B,μB,J\epsilon_F,N,B,\mu_B,J 表示相应的磁化 MM。

Kai​

(1)​

N=2V(2πℏ)34πpF33,pF=ℏ(3π2NV)1/3,N=\frac{2V}{(2\pi\hbar)^3}\frac{4\pi p_F^3}{3}, \qquad \boxed{p_F=\hbar\left(\frac{3\pi^2N}{V}\right)^{1/3}},
ϵF=ℏ22m(3π2NV)2/3,D(ϵF)=3N2ϵF.\boxed{\epsilon_F=\frac{\hbar^2}{2m}\left(\frac{3\pi^2N}{V}\right)^{2/3},\qquad D(\epsilon_F)=\frac{3N}{2\epsilon_F}}.

より一般に D(ϵ)=3N2ϵF3/2ϵD(\epsilon)=\frac{3N}{2\epsilon_F^{3/2}}\sqrt\epsilon。

(2)​

U0=3N2ϵF3/2∫0ϵFϵ3/2dϵ=35NϵF,a0=35.U_0=\frac{3N}{2\epsilon_F^{3/2}}\int_0^{\epsilon_F}\epsilon^{3/2}d\epsilon =\frac35N\epsilon_F, \qquad\boxed{a_0=\frac35}.

(3)​

NN 一定では U0∝V−2/3U_0\propto V^{-2/3} より

PV=23U0.\boxed{PV=\frac23U_0}.

(4)​

スピン一方向のフェルミ運動量は pF(1∓α)1/3p_F(1\mp\alpha)^{1/3} であるから

E±=310NϵF(1∓α)5/3.\boxed{E^\pm=\frac3{10}N\epsilon_F(1\mp\alpha)^{5/3}}.

(5)​

EZ=(N+−N−)μBB=−NαμBB.\boxed{E_Z=(N^+-N^-)\mu_BB=-N\alpha\mu_BB}.

(6)​

Etot=35NϵF+13NϵFα2−NμBBα+O(α4).E_{\rm tot}=\frac35N\epsilon_F+\frac13N\epsilon_F\alpha^2-N\mu_BB\alpha+O(\alpha^4).

微分をゼロとして

α=3μBB2ϵF,M=(N−−N+)μB=NαμB=3NμB2B2ϵF.\boxed{\alpha=\frac{3\mu_BB}{2\epsilon_F},\qquad M=(N^--N^+)\mu_B=N\alpha\mu_B=\frac{3N\mu_B^2B}{2\epsilon_F}}.

ここで MM は体積で割らない全磁気モーメントである。

(7)​

同じ電子対を二重に数えず

Eint+=−J2N+(N+−1)≃−JN28(1−α)2.\boxed{E_{\rm int}^+=-\frac J2N^+(N^+-1)\simeq-\frac{JN^2}{8}(1-\alpha)^2}.

(8)​

両スピンを加えると Eint≃−JN2(1+α2)/4E_{\rm int}\simeq-JN^2(1+\alpha^2)/4。したがって

E~tot=定数+N(ϵF3−JN4)α2−NμBBα+O(α4).\widetilde E_{\rm tot}=\text{定数} +N\left(\frac{\epsilon_F}{3}-\frac{JN}{4}\right)\alpha^2-N\mu_BB\alpha+O(\alpha^4).

仮定より二次係数は正であり、

α=3μBB2ϵF(1−3JN4ϵF)−1,M=3NμB2B2ϵF(1−3JN4ϵF)−1.\boxed{\begin{gathered} \alpha=\frac{3\mu_BB}{2\epsilon_F} \left(1-\frac{3JN}{4\epsilon_F}\right)^{-1},\\ M=\frac{3N\mu_B^2B}{2\epsilon_F} \left(1-\frac{3JN}{4\epsilon_F}\right)^{-1} \end{gathered}}.

同方向スピンを有利にする相互作用は磁化率を増大させる。

Reference​