涨落–响应关系
引理
哈密顿量 H=H0−λA(λ 与可观测量 A 共轭,在权重 e−βH 指数中),β=1/kBT。则
∂λ∂⟨A⟩=β(⟨A2⟩−⟨A⟩2).
证明. ⟨A⟩=Z1∑cAce−βH0,c+βλAc,Z=∑ce−βH0,c+βλAc。∂λeβλAc=βAceβλAc,故 ∂λZ=βZ⟨A⟩。商求导:
∂λ∂⟨A⟩=βZ∑cAc2e⋯−β⟨A⟩Z∑cAce⋯=β(⟨A2⟩−⟨A⟩2).■
磁化率 = 磁化涨落
H=H0−hΦ,Φ=∫Vφddr 总序参量。取 λ=h, A=Φ:
χtot=∂h∂⟨Φ⟩=β(⟨Φ2⟩−⟨Φ⟩2)=β∫V∫VG(r−r′)ddrddr′.
减项 −⟨Φ⟩2 自动给出连通关联 G(不是裸乘积)。平移不变性下双重积分 =V∫Gdds,单位体积
χ=kBT1∫ddrG(r).
用高斯结果 ∫G=G(k=0)=2akBT 得 χ=2a1∝∣t∣−1,即 γ=1,与朗道理论直接算出的结果一致。
热容 = 能量涨落
E 在 e−βE 中与 β 共轭。由 ⟨E⟩=−∂βlnZ 再求一次 β 导:
∂β∂⟨E⟩=−(⟨E2⟩−⟨E⟩2),CV=∂T∂⟨E⟩=kBT2⟨E2⟩−⟨E⟩2.
与磁化率同出一个引理(M 对 h、E 对 β)。⟨⋯⟩≥0 的方差给出 χ,CV≥0。
临界点发散
χ=kBT1∫G。T→Tc 时 ξ→∞,G 长程化、∫G 发散,故 χ 发散——磁化率发散的微观来源是关联涨落变长程。
(这些是静态(等时)涨落–响应;严格的涨落耗散定理是其动力学推广,把动态响应 χ(ω) 的虚部与涨落谱相连,本系列不涉及。)