tecktom.com

この章のまとめ

2026年10月15日11分で読めますAI生成教科書『虎の巻』第3部 第7章 平均場

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第7章 §3.7.7 ★★★ 自由エネルギーの最小化を積の分布に制限すると平均場の1体の HH が出る。Ising では Weiss、気体では van der Waals と Vlasov

相互作用する多体系では e−βH\mathrm{e}^{-\beta H} が1体の因子に分かれず、分配関数は一般に計算することができない。この章では、第3部第3章の自由エネルギーの最小性 F[ρ]≥FF[\rho]\ge F から Gibbs–Bogoliubov の不等式 F≤F0+⟨H−H0⟩0F\le F_0+\langle H-H_0 \rangle _0 を出し、試す分布を粒子ごとの積に限って最小にした。停留条件は、各粒子の分布が「ほかの粒子について平均した HH」を1体の Hamiltonian heffh^{\mathrm{eff}} とするカノニカル分布になる、という自己無撞着方程式になる。得られる自由エネルギーは真の値の上界で、捨てたのは粒子どうしの相関である。

Ising 模型ではこれが Weiss の分子場理論になり、m=tanh⁡β(zJm+b)m=\tanh\beta(zJm+b)、kBTc=zJk_{\mathrm B} T_{\mathrm c}=zJ が出る。TcT_{\mathrm c} の近くで FMFF_{\mathrm{MF}} を展開するとLandau の形になり、臨界指数は 1/21/2、11、33、比熱は 32NkB\frac32N k_{\mathrm B} の跳びで、次元によらない。全部のスピンが全部と J/NJ/N で結合する模型では、Stirling の公式と最大の項の評価により平均場の自由エネルギーが N→∞N\to\infty で厳密になる。一方、短距離の相互作用では、1次元に相転移は無く、2次元の TcT_{\mathrm c} は平均場の近似の 1/1.761/1.76 で、臨界指数も違う。

相空間の気体では、1体の Hamiltonian heff=∣p∣2/(2m)+U[n](r)h^{\mathrm{eff}}=\lvert \bm{p}\rvert ^2/(2m)+U[n](\bm{r}) の Liouville 方程式が Vlasov 方程式になり、平衡の平均場の分布はその定常解である。引力の尾を平均場、剛体の芯を排除体積で扱うと van der Waals の状態方程式が出る。授業で習う分子場 λM\lambda M と Curie–Weiss の θ\theta はλ=zJ/(NμB2)\lambda=zJ/(N\mu_{\mathrm B}^2)、θ=Tc\theta=T_{\mathrm c} と対応する。

この章の結果

  • Gibbs–Bogoliubov:F≤F0+⟨H−H0⟩0F\le F_0+\langle H-H_0 \rangle _0。積の分布での停留条件は ρi∝e−βhieff\rho_i\propto\mathrm{e}^{-\beta h_i^{\mathrm{eff}}}、hieffh_i^{\mathrm{eff}} はほかの粒子で平均した HH。
  • Weiss:m=tanh⁡β(zJm+b)m=\tanh\beta(zJm+b)、kBTc=zJk_{\mathrm B} T_{\mathrm c}=zJ、FMF/N=−zJ2m2−bm−kBT hb(1+m2)F_{\mathrm{MF}}/N=-\frac{zJ}{2}m^2-bm- k_{\mathrm B} T\,h_{\mathrm b}\bigl(\frac{1+m}{2}\bigr)。
  • 臨界:m∗≈3(1−T/Tc)1/2m_*\approx\sqrt3(1-T/T_{\mathrm c})^{1/2}、χ=1/(kB(T−Tc))\chi=1/( k_{\mathrm B}(T-T_{\mathrm c}))、b=kBTc3m3b=\frac{ k_{\mathrm B} T_{\mathrm c}}{3}m^3、比熱の跳び 32NkB\frac32N k_{\mathrm B}。
  • 到達距離が無限大(Curie–Weiss 模型)なら F/N→min⁡mFMF/NF/N\to\min_mF_{\mathrm{MF}}/N が厳密。
  • 気体:heff=∣p∣2/(2m)+U[n]h^{\mathrm{eff}}=\lvert \bm{p}\rvert ^2/(2m)+U[n]、Vlasov 方程式、P=NkBT/(V−Nbv)−aN2/V2P=N k_{\mathrm B} T/(V-Nb_{\mathrm v})-aN^2/V^2。

次の章(第3部第8章)では、Vlasov 方程式に衝突を足して Boltzmann 方程式にする。磁化率と磁化の揺らぎの関係χ=β(⟨M2⟩−⟨M⟩2)/N\chi=\beta(\langle M^2 \rangle -\langle M \rangle ^2)/N は 第3部第9章の揺らぎと線形応答で、式 (3.7.14) を一般化した Landau 理論は第3部第10章で扱う。量子力学の平均場(Hartree–Fock)と超伝導の平均場は 第5部第1章・第5部第6章とハミルトニアンの選び方で固体を考える 第5部第4章で、平均場の近似の厳密化の状況は 近似の数理 第13章でまとめる。

総合問題 古典の回転子の強磁性 ★★★

オリジナル格子の各点に、慣性モーメント II の直線分子の回転子が1つずつ固定されている。回転子 ii の軸の向きを単位ベクトル ni\bm{n}_i(極角 θi\theta_i、方位角 φi\varphi_i)で表し、その共役運動量を pθip_{\theta i}、pφip_{\varphi i} とする。Hamiltonian は

H=∑i=1N(pθi22I+pφi22Isin⁡2θi)−J∑⟨ij⟩ni⋅nj−b∑icos⁡θiH=\sum_{i=1}^N\Bigl(\frac{p_{\theta i}^2}{2I}+\frac{p_{\varphi i}^2}{2I\sin^2\theta_i}\Bigr)-J\sum_{\langle ij\rangle}\bm{n}_i\cdot\bm{n}_j-b\sum_i\cos\theta_i

で(J>0J>0、配位数 zz、外場は zz 方向)、温度 TT の熱浴に接している。相空間の測度は1つの回転子あたり dθ dφ dpθ dpφ/h2\mathrm{d}\theta\,\mathrm{d}\varphi\,\mathrm{d} p_\theta\,\mathrm{d} p_\varphi/h^2 とする。

  1. (1) (使う道具:3.7.2節・第1部第8章)積の分布で ⟨nj⟩=mez\langle \bm{n}_j \rangle =m\bm{e}_z とするとき、回転子 ii の1体の有効 Hamiltonian heffh^{\mathrm{eff}} を書け。
  2. (2) (使う道具:第3部第3章)heffh^{\mathrm{eff}} の分配関数 Z1Z_1 を求め、⟨cos⁡θ⟩=ℓ(x)≔coth⁡x−1/x\langle \cos\theta \rangle =\ell(x)\coloneqq\coth x-1/x(Langevin 関数)を示せ。xx は何か。
  3. (3) (使う道具:3.7.3節)自己無撞着方程式を書き、b=0b=0 で自発磁化が現れる温度 TcT_{\mathrm c} を求めよ。T=0.5TcT=0.5T_{\mathrm c} での mm を数値で求めよ。
  4. (4) (使う道具:3.7.4節)TcT_{\mathrm c} の直下の mm と、TcT_{\mathrm c} の上の磁化率 χ=∂m/∂b∣b=0\chi=\partial m/\partial b|_{b=0} を求め、Ising 模型と比べよ。
  5. (5) (使う道具:3.7.4節・第3部第6章)1つの回転子あたりの比熱を、TcT_{\mathrm c} のすぐ上とすぐ下で求めよ。
解答

(1) 相互作用 −Jni⋅nj-J\bm{n}_i\cdot\bm{n}_j を nj\bm{n}_j について平均すると −Jm ni⋅ez=−Jmcos⁡θi-Jm\,\bm{n}_i\cdot\bm{e}_z=-Jm\cos\theta_i。隣が zz 個あるので、式 (3.7.8) から

heff=pθ22I+pφ22Isin⁡2θ−(zJm+b)cos⁡θ.h^{\mathrm{eff}}=\frac{p_\theta^2}{2I}+\frac{p_\varphi^2}{2I\sin^2\theta}-(zJm+b)\cos\theta .

運動エネルギーの部分は、球面の上の自由な回転子の Hamiltonian(L=I2(θ˙2+sin⁡2θ φ˙2)L=\frac I2(\dot\theta^2+\sin^2\theta\,\dot\varphi^2) の Legendre 変換)そのもので、平均場の近似は位置エネルギーの部分だけを変える。

(2) θ\theta を止めて運動量で積分すると、∫e−βpθ2/2Idpθ=2πIkBT\int\mathrm{e}^{-\beta p_\theta^2/2I}\mathrm{d} p_\theta=\sqrt{2\pi I k_{\mathrm B} T}、∫e−βpφ2/(2Isin⁡2θ)dpφ=2πIkBTsin⁡θ\int\mathrm{e}^{-\beta p_\varphi^2/(2I\sin^2\theta)}\mathrm{d} p_\varphi=\sqrt{2\pi I k_{\mathrm B} T}\sin\theta で、積は 2πIkBTsin⁡θ2\pi I k_{\mathrm B} T\sin\theta。運動量の積分から、球面の面積要素 sin⁡θ dθ dφ\sin\theta\,\mathrm{d}\theta\,\mathrm{d}\varphi が出てくる。x=β(zJm+b)x=\beta(zJm+b) とおくと

Z1=2πIkBTh2∫02π ⁣ ⁣dφ∫0πsin⁡θ excos⁡θdθ=8π2IkBTh2⋅sinh⁡xx.Z_1=\frac{2\pi I k_{\mathrm B} T}{h^2}\int_0^{2\pi}\!\!\mathrm{d}\varphi\int_0^\pi\sin\theta\,\mathrm{e}^{x\cos\theta}\mathrm{d}\theta=\frac{8\pi^2I k_{\mathrm B} T}{h^2}\cdot\frac{\sinh x}{x}.

⟨cos⁡θ⟩=∂ln⁡Z1/∂x=ddxln⁡sinh⁡xx=coth⁡x−1x=ℓ(x)\langle \cos\theta \rangle =\partial\ln Z_1/\partial x=\frac{\mathrm{d} }{\mathrm{d} x}\ln\frac{\sinh x}{x}=\coth x-\frac1x=\ell(x)。

(3) 自己無撞着方程式は m=ℓ(β(zJm+b))m=\ell\bigl(\beta(zJm+b)\bigr)。ℓ(x)=x3−x345+⋯\ell(x)=\frac x3-\frac{x^3}{45}+\cdots で原点での傾きは 1/31/3 なので、b=0b=0 で m≠0m\neq0 の解が現れるのは βzJ/3>1\beta zJ/3>1、すなわち

kBTc=zJ3.k_{\mathrm B} T_{\mathrm c}=\frac{zJ}{3}.

Ising 模型の zJzJ の 1/31/3 である。向きが連続な n\bm{n} では、zz 成分の2乗の平均が ⟨cos⁡2θ⟩=1/3\langle \cos^2\theta \rangle =1/3 と小さいので、そろえる力が弱い。T=0.5TcT=0.5T_{\mathrm c} では βzJ=6\beta zJ=6 で、m=ℓ(6m)m=\ell(6m) を数値で解いて m=0.789m=0.789(Ising 模型の同じ T/TcT/T_{\mathrm c} では 0.9580.958)。

(4) βzJ=3Tc/T\beta zJ=3T_{\mathrm c}/T とおき、ℓ\ell を3次まで展開すると m=TcTm−2745(TcT)3m3m=\frac{T_{\mathrm c}}{T}m-\frac{27}{45}\bigl(\frac{T_{\mathrm c}}{T}\bigr)^3m^3、よってm2=53(TcT−1)(TTc)3≈53(1−TTc)m^2=\frac53\bigl(\frac{T_{\mathrm c}}{T}-1\bigr)\bigl(\frac{T}{T_{\mathrm c}}\bigr)^3\approx\frac53\bigl(1-\frac{T}{T_{\mathrm c}}\bigr)、

m≈5/3 (1−TTc)1/2.m\approx\sqrt{5/3}\,\Bigl(1-\frac{T}{T_{\mathrm c}}\Bigr)^{1/2}.

T>TcT>T_{\mathrm c} では1次の項だけで m=β3(zJm+b)m=\frac{\beta}{3}(zJm+b)、これを解いて

χ=13kB(T−Tc).\chi=\frac{1}{3 k_{\mathrm B}(T-T_{\mathrm c})}.

係数(3\sqrt3 が 5/3\sqrt{5/3} に、11 が 1/31/3 に)は Ising 模型と違うが、指数 1/21/2 と 11 は同じである。平均場の近似の臨界指数が、スピンの種類によらないことの一例になっている。

(5) 運動エネルギーは運動量の2次式が2つで、等分配則(第3部第6章)により1つの回転子あたり kBTk_{\mathrm B} T、比熱 kBk_{\mathrm B}。これは mm によらない。相互作用の部分は積の分布で Uint=−NzJ2m2U_{\mathrm{int}}=-\frac{NzJ}{2}m^2(b=0b=0)。T>TcT>T_{\mathrm c} では 00、TcT_{\mathrm c} の直下では (4) からUint≈−NzJ2⋅53(1−TTc)U_{\mathrm{int}}\approx-\frac{NzJ}{2}\cdot\frac53\bigl(1-\frac{T}{T_{\mathrm c}}\bigr) で、

dUintdT=NzJ2⋅53Tc=N⋅3kBTc2⋅53Tc=52NkB.\frac{\mathrm{d} U_{\mathrm{int}}}{\mathrm{d} T}=\frac{NzJ}{2}\cdot\frac{5}{3T_{\mathrm c}}=\frac{N\cdot3 k_{\mathrm B} T_{\mathrm c}}{2}\cdot\frac{5}{3T_{\mathrm c}}=\frac52N k_{\mathrm B} .

よって1つの回転子あたりの比熱は、TcT_{\mathrm c} のすぐ上で kBk_{\mathrm B}、すぐ下で kB+52kB=72kBk_{\mathrm B}+\frac52 k_{\mathrm B}=\frac72 k_{\mathrm B}。跳びは 52kB\frac52 k_{\mathrm B} で、Ising 模型の 32kB\frac32 k_{\mathrm B} より大きい。運動量の自由度は相互作用に関係しないので、相転移の跳びには加わらない。

木下の勉強ステータス

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

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

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