解析力学から様々な理論分野を組み立てる 第3部 第10章 §3.10.3 ★★★ f ( m ) f(m) f ( m ) の最小点が T c T_{\mathrm c} T c で 0 0 0 から2つに分かれ、m ∝ ( T c − T ) 1 / 2 m\propto(T_{\mathrm c}-T)^{1/2} m ∝ ( T c − T ) 1/2 ・比熱の跳び・χ ∝ ∣ T − T c ∣ − 1 \chi\propto\lvert T-T_{\mathrm c}\rvert ^{-1} χ ∝ ∣ T − T c ∣ − 1 が出る
最小点の分かれ方
外場 h = 0 h=0 h = 0 の Landau 関数 (3.10.2) の最小点を探す。f ′ ( m ) = m ( a + b m 2 ) = 0 f'(m)=m(a+bm^2)=0 f ′ ( m ) = m ( a + b m 2 ) = 0 の解は m = 0 m=0 m = 0 と m 2 = − a / b m^2=-a/b m 2 = − a / b で、T > T c T>T_{\mathrm c} T > T c (a > 0 a>0 a > 0 )では m = 0 m=0 m = 0 だけが最小、T < T c T<T_{\mathrm c} T < T c (a < 0 a<0 a < 0 )では m = 0 m=0 m = 0 が極大になり
m 0 = ± a 0 ( T c − T ) b (3.10.17) m_0=\pm\sqrt{\frac{a_0(T_{\mathrm c}-T)}{b}}\tag{3.10.17} m 0 = ± b a 0 ( T c − T ) ( 3.10.17 )
の2点が最小になる(図 3.10.1 )。m 0 m_0 m 0 は T c T_{\mathrm c} T c で 0 0 0 から連続に立ち上がり、その傾きは無限大である。2つの最小点のどちらが選ばれるかは f f f からは決まらず、わずかに残る外場や履歴で決まる。これが対称性の自発的な破れの、Landau 理論での姿である。
図 3.10.1 Curie–Weiss 模型の Landau 関数 f − f 0 = 1 2 T − T c T c m 2 + 1 12 m 4 f-f_0=\frac12\frac{T-T_{\mathrm c}}{T_{\mathrm c}}m^2+\frac{1}{12}m^4 f − f 0 = 2 1 T c T − T c m 2 + 12 1 m 4 (単位は k B T c k_{\mathrm B} T_{\mathrm c} k B T c 。b = k B T c / 3 b= k_{\mathrm B} T_{\mathrm c}/3 b = k B T c /3 とした)。T c T_{\mathrm c} T c より上では最小点は m = 0 m=0 m = 0 の1つ、下では ± m 0 \pm m_0 ± m 0 の2つに分かれる。
熱力学量
最小点での値は
f min = f 0 ( T ) + { 0 ( T > T c ) − a 2 4 b = − a 0 2 ( T − T c ) 2 4 b ( T < T c ) (3.10.18) f_{\min}=f_0(T)+\begin{cases}0 & (T>T_{\mathrm c})\\[2pt] -\dfrac{a^2}{4b}=-\dfrac{a_0^2(T-T_{\mathrm c})^2}{4b} & (T<T_{\mathrm c})\end{cases}\tag{3.10.18} f m i n = f 0 ( T ) + ⎩ ⎨ ⎧ 0 − 4 b a 2 = − 4 b a 0 2 ( T − T c ) 2 ( T > T c ) ( T < T c ) ( 3.10.18 )
である。f min f_{\min} f m i n そのものと、その T T T についての1階微分(エントロピー − ∂ f / ∂ T -\partial f/\partial T − ∂ f / ∂ T )は T c T_{\mathrm c} T c で連続で、2階微分が跳ぶ。比熱は C = − T ∂ 2 f / ∂ T 2 C=-T\,\partial^2 f/\partial T^2 C = − T ∂ 2 f / ∂ T 2 (第1部第7章 ・第3部第4章 )なので、T c T_{\mathrm c} T c を下から近づけたときと上から近づけたときの差は
Δ C = C ( T c − 0 ) − C ( T c + 0 ) = T c a 0 2 2 b (3.10.19) \Delta C=C(T_{\mathrm c}-0)-C(T_{\mathrm c}+0)=\frac{T_{\mathrm c}a_0^2}{2b}\tag{3.10.19} Δ C = C ( T c − 0 ) − C ( T c + 0 ) = 2 b T c a 0 2 ( 3.10.19 )
となる。低温側で比熱が大きいのは、温度を上げると秩序が壊れていく分のエントロピーが余分に吸収されるからである。
外場への応答を見る。磁化率 χ = ∂ m / ∂ h ∣ h = 0 \chi=\partial m/\partial h|_{h=0} χ = ∂ m / ∂ h ∣ h = 0 は、f ′ ( m ) = h f'(m)=h f ′ ( m ) = h を h h h について1次まで解けば得られる。T > T c T>T_{\mathrm c} T > T c では a m = h am=h am = h からχ = 1 / a \chi=1/a χ = 1/ a 、T < T c T<T_{\mathrm c} T < T c では m = m 0 + δ m m=m_0+\delta m m = m 0 + δ m と置いて f ′ ′ ( m 0 ) δ m = h f''(m_0)\,\delta m=h f ′′ ( m 0 ) δ m = h 、f ′ ′ ( m 0 ) = a + 3 b m 0 2 = 2 ∣ a ∣ f''(m_0)=a+3bm_0^2=2\lvert a\rvert f ′′ ( m 0 ) = a + 3 b m 0 2 = 2 ∣ a ∣ から
χ = { 1 a 0 ( T − T c ) ( T > T c ) 1 2 a 0 ( T c − T ) ( T < T c ) (3.10.20) \chi=\begin{cases}\dfrac{1}{a_0(T-T_{\mathrm c})} & (T>T_{\mathrm c})\\[6pt] \dfrac{1}{2a_0(T_{\mathrm c}-T)} & (T<T_{\mathrm c})\end{cases}\tag{3.10.20} χ = ⎩ ⎨ ⎧ a 0 ( T − T c ) 1 2 a 0 ( T c − T ) 1 ( T > T c ) ( T < T c ) ( 3.10.20 )
を得る。どちらの側でも ∣ T − T c ∣ − 1 \lvert T-T_{\mathrm c}\rvert ^{-1} ∣ T − T c ∣ − 1 で発散し、係数の比は2である。T = T c T=T_{\mathrm c} T = T c ちょうどでは a = 0 a=0 a = 0 で、f ′ ( m ) = h f'(m)=h f ′ ( m ) = h は b m 3 = h bm^3=h b m 3 = h 、すなわち m = ( h / b ) 1 / 3 m=(h/b)^{1/3} m = ( h / b ) 1/3 になる。Curie–Weiss 模型の係数 (3.10.16) を入れた Δ C = 3 2 k B \Delta C=\frac32 k_{\mathrm B} Δ C = 2 3 k B (1スピンあたり)と χ = 1 / ( k B ( T − T c ) ) \chi=1/( k_{\mathrm B}(T-T_{\mathrm c})) χ = 1/ ( k B ( T − T c )) は、第7章の式 (3.7.17) ・(3.7.16) と同じである(第7章の z J zJ z J をこの章の J J J と読む)。
臨界指数
T c T_{\mathrm c} T c の近くで物理量が ∣ T − T c ∣ \lvert T-T_{\mathrm c}\rvert ∣ T − T c ∣ や h h h の冪で振る舞うとき、その冪を臨界指数 と呼ぶ。Landau 理論の値(平均場の値)を表 3.10.1 にまとめる。相関長の指数 ν \nu ν と相関関数の指数 η \eta η は3.10.5 節で出す。
表 3.10.1 臨界指数。t = ( T − T c ) / T c t=(T-T_{\mathrm c})/T_{\mathrm c} t = ( T − T c ) / T c 。3次元 Ising の値は Kos ら(2016)のスケーリング次元から指数の関係式で計算した(3.10.5 節の脚注)。 指数 定義 Landau(平均場) 3次元 Ising α \alpha α C ∝ ∣ t ∣ − α C\propto\lvert t\rvert ^{-\alpha} C ∝ ∣ t ∣ − α 0 0 0 (跳び)0.110 0.110 0.110 β \beta β m 0 ∝ ( − t ) β m_0\propto(-t)^{\beta} m 0 ∝ ( − t ) β 1 / 2 1/2 1/2 0.326 0.326 0.326 γ \gamma γ χ ∝ ∣ t ∣ − γ \chi\propto\lvert t\rvert ^{-\gamma} χ ∝ ∣ t ∣ − γ 1 1 1 1.237 1.237 1.237 δ \delta δ m ∝ h 1 / δ m\propto h^{1/\delta} m ∝ h 1/ δ (t = 0 t=0 t = 0 )3 3 3 4.79 4.79 4.79 ν \nu ν ξ ∝ ∣ t ∣ − ν \xi\propto\lvert t\rvert ^{-\nu} ξ ∝ ∣ t ∣ − ν 1 / 2 1/2 1/2 0.630 0.630 0.630 η \eta η G ( r ) ∝ r − ( d − 2 + η ) G(r)\propto r^{-(d-2+\eta)} G ( r ) ∝ r − ( d − 2 + η ) (t = 0 t=0 t = 0 )0 0 0 0.0363 0.0363 0.0363
指数の値は a 0 a_0 a 0 , b b b によらない。どの物質でも f f f が m m m の解析的な偶関数であるかぎり同じ値が出るのは、Landau の展開では T c T_{\mathrm c} T c の近くで m 2 m^2 m 2 と m 4 m^4 m 4 の項しか残らず、物質の違いは係数にしか入らないからである。ただし実際の3次元の物質の指数は表の右の列のように平均場の値からずれる。このずれの原因が3.10.5 節の揺らぎである。
図 3.10.2 Curie–Weiss 模型の自発磁化。点線は自己無撞着方程式 m = tanh ( m T c / T ) m=\tanh(mT_{\mathrm c}/T) m = tanh ( m T c / T ) の厳密解(N → ∞ N\to\infty N → ∞ )、実線は Landau の展開で4次まで残した m 0 = 3 ( 1 − T / T c ) m_0=\sqrt{3(1-T/T_{\mathrm c})} m 0 = 3 ( 1 − T / T c ) 。T c T_{\mathrm c} T c の近くで2つは接し、離れると Landau の値は大きすぎる(m ≤ 1 m\le1 m ≤ 1 の制約を知らない)。
図 3.10.2 は、Curie–Weiss 模型について厳密な m 0 m_0 m 0 と Landau の値を比べたものである。Landau の展開は m m m の小さい所の展開なので、T c T_{\mathrm c} T c から離れるほど合わなくなる。
近似の裏付け
型 B(極限) | 小さいもの ∣ T − T c ∣ / T c → 0 \lvert T-T_{\mathrm c}\rvert /T_{\mathrm c}\to0 ∣ T − T c ∣ / T c → 0 (m m m が小さく、f ( m ) f(m) f ( m ) の m 6 m^6 m 6 以上の項を落とせる) | 証明 条件付き(f f f が m = 0 m=0 m = 0 で解析的で b > 0 b>0 b > 0 。Curie–Weiss 模型では式 (3.10.10) が解析的なので成り立つ。相互作用が短距離の系では揺らぎのために前提が崩れる。3.10.5 節) | 行き先 近似の数理 第2章
例題 3.10.3 Curie–Weiss 模型の比熱の跳びと自発磁化 ★★
オリジナル Curie–Weiss 模型の Landau の係数 (3.10.16) を使う。
(1) 1スピンあたりの比熱の跳び Δ C \Delta C Δ C を求めよ。スピン 1 m o l 1\,\mathrm{mol} 1 mol あたりではいくらか(気体定数 R = 8.314 J K − 1 m o l − 1 R=8.314\,\mathrm{J\,K^{-1}\,mol^{-1}} R = 8.314 J K − 1 mo l − 1 )。
(2) b = k B T c / 3 b= k_{\mathrm B} T_{\mathrm c}/3 b = k B T c /3 として、Landau の m 0 m_0 m 0 を T / T c T/T_{\mathrm c} T / T c で表せ。T = 0.95 T c T=0.95\,T_{\mathrm c} T = 0.95 T c と 0.8 T c 0.8\,T_{\mathrm c} 0.8 T c での値を、厳密な m 0 m_0 m 0 (0.379 0.379 0.379 と 0.710 0.710 0.710 )と比べよ。
(3) 厳密な Curie–Weiss 模型でも T → T c − 0 T\to T_{\mathrm c}-0 T → T c − 0 で比熱が (1) の値に近づくことを、エネルギー U / N = − J m 0 2 / 2 U/N=-Jm_0^2/2 U / N = − J m 0 2 /2 から説明せよ。
解答 (1) 式 (3.10.19) に a 0 = k B a_0= k_{\mathrm B} a 0 = k B , b = k B T c / 3 b= k_{\mathrm B} T_{\mathrm c}/3 b = k B T c /3 を入れて Δ C = T c k B 2 / ( 2 k B T c / 3 ) = 3 2 k B \Delta C=T_{\mathrm c} k_{\mathrm B}^2/(2 k_{\mathrm B} T_{\mathrm c}/3)=\frac32 k_{\mathrm B} Δ C = T c k B 2 / ( 2 k B T c /3 ) = 2 3 k B 。1 m o l 1\,\mathrm{mol} 1 mol では 3 2 R = 12.5 J K − 1 m o l − 1 \frac32R=12.5\,\mathrm{J\,K^{-1}\,mol^{-1}} 2 3 R = 12.5 J K − 1 mo l − 1 。T c T_{\mathrm c} T c より上では m = 0 m=0 m = 0 なので、この模型の比熱は 0 0 0 で、T c T_{\mathrm c} T c で 0 0 0 から 3 2 k B \frac32 k_{\mathrm B} 2 3 k B に跳ぶ。
(2) m 0 2 = a 0 ( T c − T ) / b = 3 ( 1 − T / T c ) m_0^2=a_0(T_{\mathrm c}-T)/b=3(1-T/T_{\mathrm c}) m 0 2 = a 0 ( T c − T ) / b = 3 ( 1 − T / T c ) 。T = 0.95 T c T=0.95\,T_{\mathrm c} T = 0.95 T c で m 0 = 0.15 = 0.387 m_0=\sqrt{0.15}=0.387 m 0 = 0.15 = 0.387 (厳密 0.379 0.379 0.379 、ずれ 2 % 2\% 2% )、T = 0.8 T c T=0.8\,T_{\mathrm c} T = 0.8 T c で 0.6 = 0.775 \sqrt{0.6}=0.775 0.6 = 0.775 (厳密 0.710 0.710 0.710 、ずれ 9 % 9\% 9% )。
(3) 自己無撞着方程式 m = tanh ( m T c / T ) m=\tanh(mT_{\mathrm c}/T) m = tanh ( m T c / T ) を m m m の3次まで展開すると m ≈ T c T m − 1 3 ( T c T ) 3 m 3 m\approx\frac{T_{\mathrm c}}{T}m-\frac13\bigl(\frac{T_{\mathrm c}}{T}\bigr)^3m^3 m ≈ T T c m − 3 1 ( T T c ) 3 m 3 で、T c T_{\mathrm c} T c の近くでは m 0 2 ≈ 3 ( 1 − T / T c ) m_0^2\approx3(1-T/T_{\mathrm c}) m 0 2 ≈ 3 ( 1 − T / T c ) となり (2) と一致する。U / N = − J 2 m 0 2 ≈ − 3 2 k B ( T c − T ) U/N=-\frac{J}{2}m_0^2\approx-\frac32 k_{\mathrm B}(T_{\mathrm c}-T) U / N = − 2 J m 0 2 ≈ − 2 3 k B ( T c − T ) を T T T で微分すると 3 2 k B \frac32 k_{\mathrm B} 2 3 k B 。(数値でも T = 0.998 T c T=0.998\,T_{\mathrm c} T = 0.998 T c 〜0.999 T c 0.999\,T_{\mathrm c} 0.999 T c の差分で 1.496 k B 1.496\, k_{\mathrm B} 1.496 k B になる。)
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com