tecktom.com

出てくるもの (1):1ループの再帰式と Wilson–Fisher 不動点

2026年10月15日12分で読めますAI生成教科書『虎の巻』第5部 第11章 くりこみ群
目次
  1. 殻の積分
  2. 流れの方程式
  3. 2つの固定点
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第5部 第11章 §5.11.3 ★★★ 殻の積分を λ\lambda の2次まで行うと rr と uu の流れが閉じる。d<4d<4 で Gauss の固定点は不安定になり、u∗=2ϵ/3u^*=2\epsilon/3 の新しい固定点で ν=12+ϵ12\nu=\frac12+\frac{\epsilon}{12}

殻の積分

相互作用 Sint=λ4!∫ϕ4S_{\mathrm{int}}=\frac{\lambda}{4!}\int\phi^4 を λ\lambda について展開して式 (5.11.6) を計算する。ϕ>\phi_> だけの Gauss 積分の期待値を⟨⋅⟩>\langle \cdot \rangle _> と書くと、キュムラント展開で

SE′[ϕ<]=S0[ϕ<]+⟨Sint⟩>−12(⟨Sint2⟩>−⟨Sint⟩>2)+O(λ3)+(定数)(5.11.10)S'_{\mathrm E}[\phi_<]=S_0[\phi_<]+\langle S_{\mathrm{int}} \rangle _>-\frac12\bigl(\langle S_{\mathrm{int}}^2 \rangle _>-\langle S_{\mathrm{int}} \rangle _>^2\bigr)+\mathcal{O}(\lambda^3)+(\text{定数})\tag{5.11.10}

である。ϕ>\phi_> の伝播関数は ⟨ϕ~>(k)ϕ~>(k′)⟩>=(2π)dδ(d)(k+k′)/(k2+r)\langle \tilde\phi_>(k)\tilde\phi_>(k') \rangle _>=(2\pi)^d\delta^{(d)}(k+k')/(k^2+r) で、第5部第7章の Wick の定理で期待値を組み合わせに分ける。

ϕ4=(ϕ<+ϕ>)4\phi^4=(\phi_<+\phi_>)^4 のうち、ϕ>\phi_> について2次の項は 6ϕ<2ϕ>26\phi_<^2\phi_>^2 である。1次の項 ⟨Sint⟩>\langle S_{\mathrm{int}} \rangle _> からはλ4!⋅6 ϕ<2⟨ϕ>2⟩>\frac{\lambda}{4!}\cdot6\,\phi_<^2\langle \phi_>^2 \rangle _> が出て、これは ϕ<2\phi_<^2 の係数を変える。⟨ϕ>(x)2⟩>=∫shellddk(2π)d1k2+r≕I\langle \phi_>(x)^2 \rangle _>=\int_{\mathrm{shell}}\frac{\mathrm{d}^dk}{(2\pi)^d}\frac{1}{k^2+r}\eqqcolon I と書くと

12 Δr ϕ<2=λ4I ϕ<2,Δr=λ2I.\frac12\,\Delta r\,\phi_<^2=\frac{\lambda}{4}I\,\phi_<^2,\qquad \Delta r=\frac{\lambda}{2}I.

2次の項からは、xx と yy の2点の ϕ<2ϕ>2\phi_<^2\phi_>^2 どうしの ϕ>\phi_> を2本ずつつないだ寄与 −12(λ4)2⋅2∫ ⁣ ⁣∫ϕ<2(x)ϕ<2(y)G>(x−y)2-\frac12(\frac{\lambda}{4})^2\cdot2\int\!\!\int\phi_<^2(x)\phi_<^2(y)G_>(x-y)^2 が出る。G>G_> は殻の波数だけでできているので、∣x−y∣∼1/Λ\lvert x-y\rvert \sim1/\Lambda の範囲で減衰する。この範囲で ϕ<\phi_< はほとんど変わらないので ϕ<2(y)≈ϕ<2(x)\phi_<^2(y)\approx\phi_<^2(x) と置き(局所近似)、∫G>(z)2ddz=∫shellddk(2π)d1(k2+r)2≕J\int G_>(z)^2\mathrm{d}^dz=\int_{\mathrm{shell}}\frac{\mathrm{d}^dk}{(2\pi)^d}\frac{1}{(k^2+r)^2}\eqqcolon J を使うと、ϕ<4\phi_<^4 の係数の変化は Δλ/4!=−λ216J\Delta\lambda/4!=-\frac{\lambda^2}{16}J、すなわち

Δr=λ2 I,Δλ=−32λ2J(5.11.11)\Delta r=\frac{\lambda}{2}\,I,\qquad \Delta\lambda=-\frac32\lambda^2J\tag{5.11.11}

を得る。Δr\Delta r は第5部第8章の「おたまじゃくし」の図、Δλ\Delta\lambda は3つのチャネルの「泡」の図を、殻の運動量だけで計算したものである。同じ組み合わせから ϕ<6\phi_<^6 の項や ϕ<2(∇ϕ<)2\phi_<^2(\nabla\phi_<)^2 の項も現れるが、それらは式 (5.11.7) で y<0y<0 の係数で、次の節の固定点では ϵ2\epsilon^2 の大きさにとどまるので捨てる。1ループでは勾配の項の係数は変わらない(II が外の運動量によらない)ので、ζ\zeta は相互作用の無いときと同じ e(d−2)Δℓ/2\mathrm{e}^{(d-2)\Delta\ell/2} でよい。

流れの方程式

殻が薄いとき ∫shellddk(2π)df(∣k∣)=KdΛdf(Λ) Δℓ\int_{\mathrm{shell}}\frac{\mathrm{d}^dk}{(2\pi)^d}f(\lvert k\rvert )=K_d\Lambda^df(\Lambda)\,\Delta\ell、Kd≔Sd(2π)d=2(4π)d/2Γ(d/2)K_d\coloneqq\frac{S_d}{(2\pi)^d}=\frac{2}{(4\pi)^{d/2}\Gamma(d/2)}(SdS_d は単位球面の面積。K4=1/(8π2)K_4=1/(8\pi^2)、K3=1/(2π2)K_3=1/(2\pi^2))である。式 (5.11.11) に尺度の取り替え(式 (5.11.7) の e2Δℓ\mathrm{e}^{2\Delta\ell} と eϵ Δℓ\mathrm{e}^{\epsilon\,\Delta\ell})を合わせ、Λ\Lambda を単位にした無次元の量rˉ≔r/Λ2\bar r\coloneqq r/\Lambda^2、u≔KdλΛ−ϵu\coloneqq K_d\lambda\Lambda^{-\epsilon} で書くと

drˉdℓ=2rˉ+u2(1+rˉ),dudℓ=ϵu−3u22(1+rˉ)2(5.11.12)\frac{\mathrm{d} \bar r}{\mathrm{d} \ell}=2\bar r+\frac{u}{2(1+\bar r)},\qquad \frac{\mathrm{d} u}{\mathrm{d} \ell}=\epsilon u-\frac{3u^2}{2(1+\bar r)^2}\tag{5.11.12}

になる。これが1ループの再帰式である。これ以降は rˉ\bar r の上の棒を省いて rr と書く。

2つの固定点

式 (5.11.12) の固定点は2つある。1つは (r∗,u∗)=(0,0)(r^*,u^*)=(0,0) で、相互作用の無い理論である(Gauss 固定点)。ここでの行列 MM は上三角で、固有値は 22 と ϵ\epsilon である。d>4d>4(ϵ<0\epsilon<0)では uu の方向は固定点に近づき、Gauss 固定点が臨界点を支配する。そこでの指数は平均場の値になり、第3部第10章で Ginzburg の判定条件が d>4d>4 で Landau 理論を支持したことと一致する。d<4d<4 では uu の方向も外へ流れ出し、Gauss 固定点は臨界点を支配することができない。

もう1つの固定点は u≠0u\neq0 のもので、dudℓ=0\frac{\mathrm{d} u}{\mathrm{d} \ell}=0 から u∗=2ϵ3(1+r∗)2u^*=\frac{2\epsilon}{3}(1+r^*)^2、これを drdℓ=0\frac{\mathrm{d} r}{\mathrm{d} \ell}=0 に入れて

r∗=−ϵ6+ϵ=−ϵ6+O(ϵ2),u∗=2ϵ3+O(ϵ2)(5.11.13)r^*=-\frac{\epsilon}{6+\epsilon}=-\frac{\epsilon}{6}+\mathcal{O}(\epsilon^2),\qquad u^*=\frac{2\epsilon}{3}+\mathcal{O}(\epsilon^2)\tag{5.11.13}

を得る。これを Wilson–Fisher 固定点と呼ぶ1K. G. Wilson, M. E. Fisher (1972), Phys. Rev. Lett. 28, 240。下の脚注へ。ϵ>0\epsilon>0 で u∗>0u^*>0 となり、物理的に意味のある(ϕ4\phi^4 の係数が正で重みが積分できる)領域に入る。u∗u^* が ϵ\epsilon に比例して小さいことが、λ\lambda についての展開をϵ\epsilon についての展開として正当化する根拠になる。

固定点のまわりの行列 MM を ϵ\epsilon の1次まで求める。

M=(2−u∗2(1+r∗)212(1+r∗)3u∗2(1+r∗)3ϵ−3u∗(1+r∗)2)=(2−ϵ312O(ϵ2)−ϵ)+O(ϵ2).M=\begin{pmatrix}2-\dfrac{u^*}{2(1+r^*)^2}&\dfrac{1}{2(1+r^*)}\\[8pt]\dfrac{3u^{*2}}{(1+r^*)^3}&\epsilon-\dfrac{3u^*}{(1+r^*)^2}\end{pmatrix} =\begin{pmatrix}2-\dfrac{\epsilon}{3}&\dfrac12\\[4pt]\mathcal{O}(\epsilon^2)&-\epsilon\end{pmatrix}+\mathcal{O}(\epsilon^2).

左下が ϵ2\epsilon^2 の大きさなので、固有値は対角成分で ϵ\epsilon の1次まで決まり

yt=2−ϵ3,yu=−ϵ,ν=1yt=12+ϵ12+O(ϵ2)(5.11.14)y_t=2-\frac{\epsilon}{3},\qquad y_u=-\epsilon,\qquad \nu=\frac{1}{y_t}=\frac12+\frac{\epsilon}{12}+\mathcal{O}(\epsilon^2)\tag{5.11.14}

になる。Wilson–Fisher 固定点には relevant な方向が1つ(温度)しかなく、uu の方向は固定点に近づく。したがって、臨界温度を1つ合わせるだけで、どの u>0u>0 から出発しても同じ固定点に行き着く。図 5.11.1 にこの流れを示す。

図

図 5.11.1 d=3d=3(ϵ=1\epsilon=1)の1ループの再帰式 (5.11.12) の流れ(近似なので実線)。G は Gauss 固定点、WF は Wilson–Fisher 固定点 (r∗,u∗)=(−1/7, 0.490)(r^*,u^*)=(-1/7,\,0.490)。太線は WF に近づく臨界線(臨界面)で、その上の点はすべて WF に行き着く。臨界線の両側から出た細線は、WF の近くを通ってから yt>0y_t>0 の方向に沿って高温相(r→+∞r\to+\infty)か低温相(r→−∞r\to-\infty)へ離れていく。G から出た細線は uu が小さい所から出発しても WF の近くに引き寄せられる

固定点の r∗r^* が負であることにも意味がある。rr は Landau 理論の a/κ∝T−Tc(0)a/\kappa\propto T-T_{\mathrm c}^{(0)} なので、臨界点 r=r∗<0r=r^*<0 は Landau 理論の転移温度より低い。揺らぎが秩序を壊す向きに働き、本当の転移温度 TcT_{\mathrm c} は Tc(0)T_{\mathrm c}^{(0)} より下がる。これは非普遍的な量で、Λ\Lambda や uu の出発点に依存する。

近似の裏付け

型 B(極限) | 小さいもの u∗∝ϵ→0u^*\propto\epsilon\to0(ϵ=4−d\epsilon=4-d。殻の積分を λ\lambda の2次で打ち切り、ϕ6\phi^6 と微分の項を捨てる) | 証明 要確認(捨てた項が固定点で ϵ2\epsilon^2 の大きさにとどまることは Wilson–Fisher(1972)の議論による。短距離の d=3d=3 の ϕ4\phi^4 理論で非自明な固定点が存在することの数学的な証明の状況は原論文で確かめていない) | 行き先 近似の数理 第13章

例題 5.11.3 ϵ=1\epsilon=1 で再帰式を展開せずに解く ★★

オリジナル式 (5.11.12) を ϵ=1\epsilon=1(d=3d=3)とし、ϵ\epsilon で展開せずにそのまま扱う。

  1. (1) 固定点 (r∗,u∗)(r^*,u^*) を求めよ。
  2. (2) 固定点での行列 MM を数値で求め、その固有値と ν=1/yt\nu=1/y_t を求めよ。
  3. (3) 式 (5.11.14) の ν=12+ϵ12\nu=\frac12+\frac{\epsilon}{12} の ϵ=1\epsilon=1 での値、および ν=1/(2−ϵ/3)\nu=1/(2-\epsilon/3) と比べよ。3つの値が違うことは何を意味するか。
解答

(1) 式 (5.11.13) に ϵ=1\epsilon=1 を入れて r∗=−1/7=−0.1429r^*=-1/7=-0.1429、u∗=23(67)2=0.4898u^*=\frac23\bigl(\frac67\bigr)^2=0.4898。

(2) 1+r∗=6/71+r^*=6/7 を使うと M11=2−u∗2(1+r∗)2=2−13=1.667M_{11}=2-\frac{u^*}{2(1+r^*)^2}=2-\frac13=1.667、M12=12(1+r∗)=712=0.583M_{12}=\frac{1}{2(1+r^*)}=\frac{7}{12}=0.583、M21=3u∗2(1+r∗)3=87=1.143M_{21}=\frac{3u^{*2}}{(1+r^*)^3}=\frac87=1.143、M22=1−3u∗(1+r∗)2=1−2=−1M_{22}=1-\frac{3u^*}{(1+r^*)^2}=1-2=-1。固有値は 23⋅12±(43)2+712⋅87\frac{2}{3}\cdot\frac12\pm\sqrt{(\frac{4}{3})^2+\frac{7}{12}\cdot\frac{8}{7}} の形で、数値では yt=1.897y_t=1.897、yu=−1.230y_u=-1.230。したがって ν=1/1.897=0.527\nu=1/1.897=0.527。

(3) ϵ\epsilon の1次の式では ν=0.583\nu=0.583、ν=1/(2−ϵ/3)=1/(5/3)\nu=1/(2-\epsilon/3)=1/(5/3) では 0.6000.600 である。3次元 Ising の値 0.6300.630(5.11.5節)に最も近いのは1次の展開で、展開しないで解いた 0.5270.527 は最も遠い。再帰式 (5.11.12) 自体が λ\lambda の2次で打ち切ったものなので、ϵ\epsilon の2次以上の部分は捨てた項と同じ大きさで、意味を持たない。展開しないで解くと、意味の無い高次の部分が数値に混ざる。3つの値の差(0.070.07 程度)は、1ループの結果が ϵ=1\epsilon=1 で持つ不確かさの目安である。

1 K. G. Wilson, M. E. Fisher (1972), Phys. Rev. Lett. 28, 240。

木下の勉強ステータス

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

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

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