统计物理 · 第一部分 预备知识 · 预备章 P3

电磁学与物质的电、磁性质要点

Electromagnetism and the Electric & Magnetic Properties of Matter
已完成讲义更新于 2026.10.08统计物理讲义 v1.0
本章目标

准备后面各部分要用的电磁学结论:偶极矩与极化、磁化与居里定律、物质磁性的分类(以及"经典统计中没有磁性"这一定理)、分子间作用力、电磁波与腔模、光散射与密度涨落、电路中的热噪声。

物理输入

库仑定律与静电势;均匀电场中 φ(r)=−E⋅r+\varphi(\mathbf r) = -\mathbf E\cdot\mathbf r + 常数;磁矩 μ\boldsymbol\mu 在磁场中的能量为 −μ⋅B-\boldsymbol\mu\cdot\mathbf B;两个相距 r\mathbf r 的电偶极子之间的相互作用能为 U=14πε0r3[d1⋅d2−3(d1⋅r^)(d2⋅r^)]U = \frac{1}{4\pi\varepsilon_0r^3}\left[\mathbf d_1\cdot\mathbf d_2 - 3(\mathbf d_1\cdot\hat{\mathbf r})(\mathbf d_2\cdot\hat{\mathbf r})\right](磁偶极子把 1/4πε01/4\pi\varepsilon_0 换成 μ0/4π\mu_0/4\pi)。

§P3.1偶极矩、极化与磁化

电偶极矩。一组电荷的偶极矩 d=∑iqiri\mathbf d = \sum_iq_i\mathbf r_i。若体系总电荷为零,把原点平移 a\mathbf a,d\mathbf d 变为 ∑iqi(ri−a)=d−a∑iqi=d\sum_iq_i(\mathbf r_i - \mathbf a) = \mathbf d - \mathbf a\sum_iq_i = \mathbf d:偶极矩与原点的选取无关。在均匀电场中,体系的静电能为

U=∑iqiφ(ri)=−E⋅∑iqiri=−d⋅E(P3.1)U = \sum_iq_i\varphi(\mathbf r_i) = -\mathbf E\cdot\sum_iq_i\mathbf r_i = -\mathbf d\cdot\mathbf E \tag{P3.1}

感生偶极矩。把原子中的电子云看成电荷 qq、质量 mm、固有频率 ω0\omega_0 的谐振子(§P1.4)。静电场中的平衡位移为 x=qE/(mω02)x = qE/(m\omega_0^2),所以

d=αE,α=q2mω02(P3.2)\mathbf d = \alpha\mathbf E,\qquad \alpha = \frac{q^2}{m\omega_0^2} \tag{P3.2}

α\alpha 称为极化率;在频率为 ω\omega 的电场中,由 (P1.12),α(ω)=q2χ(ω)\alpha(\omega) = q^2\chi(\omega)。

高温下的线性响应。设 H^=H^0−AF\hat H = \hat H_0 - AF,FF 是外加的"力"(电场、磁场等),AA 是与之耦合的量(偶极矩分量、磁矩分量等),且无外场时 ⟨A⟩0=0\langle A\rangle_0 = 0。若 AF≪kBTAF\ll k_{\mathrm B}T,把玻尔兹曼因子中的 eβAFe^{\beta AF} 展开到一阶:

⟨A⟩=⟨A eβAF⟩0⟨eβAF⟩0≈⟨A⟩0+βF⟨A2⟩01+βF⟨A⟩0=⟨A2⟩0kBT F(P3.3)\langle A\rangle = \frac{\langle A\,e^{\beta AF}\rangle_0}{\langle e^{\beta AF}\rangle_0}\approx\frac{\langle A\rangle_0 + \beta F\langle A^2\rangle_0}{1 + \beta F\langle A\rangle_0} = \frac{\langle A^2\rangle_0}{k_{\mathrm B}T}\,F \tag{P3.3}

(⟨⋅⟩0\langle\cdot\rangle_0 表示无外场时的平均。)响应系数等于无外场时涨落的平方平均除以 kBTk_{\mathrm B}T——这是第24章涨落–耗散定理的最简单形式。(对量子系统,当 AA 与 H^0\hat H_0 对易时上式严格成立,例如下面的磁矩问题。)

取向极化。具有固定偶极矩 d0d_0 的极性分子在热运动下取向随机。由 (P3.3),并用各向同性平均 ⟨dz2⟩0=d02⟨cos⁡2θ⟩=d02/3\langle d_z^2\rangle_0 = d_0^2\langle\cos^2\theta\rangle = d_0^2/3,分子的总极化率为

α总=α+d023kBT(P3.4)\alpha_{\text{总}} = \alpha + \frac{d_0^2}{3k_{\mathrm B}T} \tag{P3.4}

(德拜公式。)测量介电常数随温度的变化,就可以把两部分分开,求出分子的偶极矩。

磁化与居里定律。磁化强度 MM 是单位体积的平均磁矩,M=n⟨μz⟩M = n\langle\mu_z\rangle;磁场强度 HH 由 B=μ0(H+M)B = \mu_0(H + M) 定义,磁化率 χm=M/H\chi_m = M/H。对稀薄的顺磁体 M≪HM\ll H,B≈μ0HB\approx\mu_0H。总角动量为 JJ 的磁性离子,μz=−gJμBmJ\mu_z = -g_J\mu_{\mathrm B}m_J(§P2.5),它与 H^0\hat H_0 对易,由 (P3.3),⟨μz⟩=⟨μz2⟩0B/kBT\langle\mu_z\rangle = \langle\mu_z^2\rangle_0B/k_{\mathrm B}T。无外场时 2J+12J+1 个 mJm_J 等概率,由 ∑m=−JJm2=J(J+1)(2J+1)/3\sum_{m=-J}^Jm^2 = J(J+1)(2J+1)/3(可用归纳法证明),⟨mJ2⟩0=J(J+1)/3\langle m_J^2\rangle_0 = J(J+1)/3。于是

χm=μ0n μeff23kBT,μeff2=gJ2μB2J(J+1)(P3.5)\chi_m = \frac{\mu_0n\,\mu_{\mathrm{eff}}^2}{3k_{\mathrm B}T},\qquad \mu_{\mathrm{eff}}^2 = g_J^2\mu_{\mathrm B}^2J(J+1) \tag{P3.5}

这就是居里定律。对 J=1/2J = 1/2、gJ=2g_J = 2,μeff2/3=μB2\mu_{\mathrm{eff}}^2/3 = \mu_{\mathrm B}^2,与 (6.7)(μm=μB\mu_{\mathrm m} = \mu_{\mathrm B})一致。测得的 μeff\mu_{\mathrm{eff}} 可以用来确定离子的 JJ。

§P3.2物质磁性的唯象图像

  • 抗磁性:χm<0\chi_m<0,量级约 −10−5-10^{-5},几乎与温度无关;所有物质都有(外场使电子轨道运动产生反向的感应磁矩)。
  • 顺磁性:χm>0\chi_m>0,服从居里定律 (P3.5);来自原子或离子的固有磁矩。
  • 铁磁性:低于居里温度 TcT_{\mathrm c} 时出现自发磁化;高于 TcT_{\mathrm c} 时实验上满足居里–外斯定律 χm=C/(T−Tc)\chi_m = C/(T - T_{\mathrm c})(第18章将用平均场理论推导)。
  • 反铁磁性:低于奈尔温度 TNT_{\mathrm N} 时相邻磁矩反向排列;高温下 χm=C/(T+θ)\chi_m = C/(T + \theta)。

玻尔–范莱文定理:经典统计中没有磁性。

物理输入

磁场中带电粒子的哈密顿量为 H=∑i(pi−qiA(ri))22mi+V(r1,…,rN)\mathcal H = \sum_i\frac{(\mathbf p_i - q_i\mathbf A(\mathbf r_i))^2}{2m_i} + V(\mathbf r_1,\dots,\mathbf r_N),其中 ∇×A=B\nabla\times\mathbf A = \mathbf B。

经典配分函数 Z=1N!h3N∫e−βHd3Nr d3NpZ = \frac{1}{N!h^{3N}}\int e^{-\beta\mathcal H}d^{3N}r\,d^{3N}p。对固定的坐标,把每个动量积分变量平移为 pi′=pi−qiA(ri)\mathbf p_i' = \mathbf p_i - q_i\mathbf A(\mathbf r_i)(平移的雅可比行列式为 1,积分限仍为全空间),A\mathbf A 就从被积函数中完全消失。所以 ZZ 与 BB 无关,

M∝−∂F∂B=kBT∂ln⁡Z∂B=0(P3.6)M\propto-\frac{\partial F}{\partial B} = k_{\mathrm B}T\frac{\partial\ln Z}{\partial B} = 0 \tag{P3.6}

经典统计力学预言任何物质在热平衡时都没有磁化。磁性本质上是量子效应:自旋与量子化的轨道角动量不能用经典动量平移消去。

为什么铁磁性不是磁偶极相互作用造成的。两个相距 a≈0.25a\approx0.25 nm、磁矩为 μB\mu_{\mathrm B} 的偶极子之间的相互作用能约为 μ04πμB2a3≈5×10−25\frac{\mu_0}{4\pi}\frac{\mu_{\mathrm B}^2}{a^3}\approx5\times10^{-25} J,相当于 kB×0.04k_{\mathrm B}\times0.04 K;而铁的居里温度约为 1043 K。使磁矩平行排列的是强得多的交换作用(§P2.7),它来自静电相互作用与泡利原理。

§P3.3分子间作用力

取向力(基索姆力)。两个固有偶极矩为 d1d_1、d2d_2 的极性分子,相距 rr(取连线为 zz 轴)。由开头的偶极相互作用公式,

U=d1d24πε0r3(sin⁡θ1sin⁡θ2cos⁡ϕ−2cos⁡θ1cos⁡θ2),ϕ=φ1−φ2U = \frac{d_1d_2}{4\pi\varepsilon_0r^3}\left(\sin\theta_1\sin\theta_2\cos\phi - 2\cos\theta_1\cos\theta_2\right),\qquad \phi = \varphi_1 - \varphi_2

对取向作热平均。与 (P3.3) 同理,⟨U⟩=⟨Ue−βU⟩0/⟨e−βU⟩0≈⟨U⟩0−β(⟨U2⟩0−⟨U⟩02)\langle U\rangle = \langle Ue^{-\beta U}\rangle_0/\langle e^{-\beta U}\rangle_0\approx\langle U\rangle_0 - \beta\left(\langle U^2\rangle_0 - \langle U\rangle_0^2\right),而随机取向下 ⟨U⟩0=0\langle U\rangle_0 = 0。计算 ⟨U2⟩0\langle U^2\rangle_0:交叉项含 cos⁡ϕ\cos\phi 的一次方,平均为零;再用 ⟨sin⁡2θ⟩=2/3\langle\sin^2\theta\rangle = 2/3、⟨cos⁡2θ⟩=1/3\langle\cos^2\theta\rangle = 1/3、⟨cos⁡2ϕ⟩=1/2\langle\cos^2\phi\rangle = 1/2,括号的平方平均为 23⋅23⋅12+4⋅13⋅13=23\frac23\cdot\frac23\cdot\frac12 + 4\cdot\frac13\cdot\frac13 = \frac23。于是

⟨U⟩≈−23kBT(d1d24πε0)21r6(P3.7)\langle U\rangle\approx-\frac{2}{3k_{\mathrm B}T}\left(\frac{d_1d_2}{4\pi\varepsilon_0}\right)^2\frac{1}{r^6} \tag{P3.7}

热平均使得吸引的取向更常出现,净效果是吸引,且随温度升高而减弱。

色散力(伦敦力)。非极性分子之间也有吸引:它是一个量子力学的零点能效应。用 (P3.2) 的模型,把两个原子各看成一个三维谐振子(电荷 qq、质量 mm、频率 ω0\omega_0),电子的位移分别为 u1\mathbf u_1、u2\mathbf u_2,偶极矩 di=qui\mathbf d_i = q\mathbf u_i。由偶极相互作用公式(连线为 zz 轴),耦合能为 q24πε0r3(x1x2+y1y2−2z1z2)\frac{q^2}{4\pi\varepsilon_0r^3}(x_1x_2 + y_1y_2 - 2z_1z_2)。三个方向互不耦合。对每个方向,两个振子的耦合能写成 c x1x2c\,x_1x_2(xx、yy 方向 c=q2/4πε0r3c = q^2/4\pi\varepsilon_0r^3,zz 方向 c=−2q2/4πε0r3c = -2q^2/4\pi\varepsilon_0r^3)。由 §P1.3,简正坐标 (x1±x2)/2(x_1\pm x_2)/\sqrt2 的频率为 ω±=ω01±ϵ\omega_\pm = \omega_0\sqrt{1\pm\epsilon},ϵ=c/(mω02)\epsilon = c/(m\omega_0^2)。两个模的零点能之和,用 1+ϵ≈1+ϵ2−ϵ28\sqrt{1+\epsilon}\approx1 + \frac\epsilon2 - \frac{\epsilon^2}{8}(提示 A3):

ℏ2(ω++ω−)≈ℏω0(1−ϵ28)\frac\hbar2(\omega_+ + \omega_-)\approx\hbar\omega_0\left(1 - \frac{\epsilon^2}{8}\right)

能量降低了 ℏω0ϵ2/8\hbar\omega_0\epsilon^2/8。对三个方向求和,∑c2=(1+1+4)(q24πε0r3)2\sum c^2 = (1 + 1 + 4)\left(\frac{q^2}{4\pi\varepsilon_0r^3}\right)^2,再用 α=q2/(mω02)\alpha = q^2/(m\omega_0^2):

ΔE=−34 ℏω0 α2(4πε0)2r6(P3.8)\Delta E = -\frac34\,\frac{\hbar\omega_0\,\alpha^2}{(4\pi\varepsilon_0)^2r^6} \tag{P3.8}

两种吸引都正比于 r−6r^{-6}。

伦纳德–琼斯势。距离很近时,电子云的交叠受泡利原理阻碍,产生强烈的排斥;实用上用 r−12r^{-12} 描述:

V(r)=4ε[(σr)12−(σr)6](P3.9)V(r) = 4\varepsilon\left[\left(\frac\sigma r\right)^{12} - \left(\frac\sigma r\right)^6\right] \tag{P3.9}

由 dV/dr=0dV/dr = 0 得 (σ/r)6=1/2(\sigma/r)^6 = 1/2,即极小值位于 rmin⁡=21/6σr_{\min} = 2^{1/6}\sigma,深度为 V(rmin⁡)=4ε(14−12)=−εV(r_{\min}) = 4\varepsilon(\frac14 - \frac12) = -\varepsilon。对氩,ε/kB≈120\varepsilon/k_{\mathrm B}\approx120 K,σ≈0.34\sigma\approx0.34 nm。第17章计算实际气体的位力系数、第32章的分子模拟都将用到它。

§P3.4电磁波与谐振腔中的模式

物理输入

真空中(无电荷、无电流)的麦克斯韦方程组:∇⋅E=0\nabla\cdot\mathbf E = 0,∇⋅B=0\nabla\cdot\mathbf B = 0,∇×E=−∂B/∂t\nabla\times\mathbf E = -\partial\mathbf B/\partial t,∇×B=μ0ε0 ∂E/∂t\nabla\times\mathbf B = \mu_0\varepsilon_0\,\partial\mathbf E/\partial t。

矢量分析

梯度 ∇f\nabla f、散度 ∇⋅v\nabla\cdot\mathbf v、旋度 ∇×v\nabla\times\mathbf v;恒等式 ∇×(∇×v)=∇(∇⋅v)−∇2v\nabla\times(\nabla\times\mathbf v) = \nabla(\nabla\cdot\mathbf v) - \nabla^2\mathbf v。(详见附录 A11,恒等式的验证见 §A11.6。)

对第三式取旋度,用上述恒等式与 ∇⋅E=0\nabla\cdot\mathbf E = 0,再代入第四式:

∇2E=1c2∂2E∂t2,c=1μ0ε0(P3.10)\nabla^2\mathbf E = \frac{1}{c^2}\frac{\partial^2\mathbf E}{\partial t^2},\qquad c = \frac{1}{\sqrt{\mu_0\varepsilon_0}} \tag{P3.10}

这是波动方程。平面波 E=E0ei(k⋅r−ωt)\mathbf E = \mathbf E_0e^{i(\mathbf k\cdot\mathbf r - \omega t)} 满足它的条件是 ω=ck\omega = ck;∇⋅E=0\nabla\cdot\mathbf E = 0 要求 k⋅E0=0\mathbf k\cdot\mathbf E_0 = 0:电磁波是横波,每个 k\mathbf k 有两个独立的偏振方向。

腔中的模式数。边长为 LL 的立方腔,用周期性边界条件(§P2.2),k=2πL(nx,ny,nz)\mathbf k = \frac{2\pi}{L}(n_x,n_y,n_z),每个 k\mathbf k 在 kk 空间占体积 (2π)3/V(2\pi)^3/V。频率不超过 ω\omega 的模式数为 2×V(2π)3⋅4π3(ωc)3=Vω33π2c32\times\frac{V}{(2\pi)^3}\cdot\frac{4\pi}{3}\left(\frac\omega c\right)^3 = \frac{V\omega^3}{3\pi^2c^3}(因子 2 来自偏振),对 ω\omega 求导得模式密度:

g(ω)=Vω2π2c3(P3.11)g(\omega) = \frac{V\omega^2}{\pi^2c^3} \tag{P3.11}

每个模式是一个谐振子。把场按模式展开,每个模式的振幅满足谐振子方程,其能量具有 12(Q˙2+ω2Q2)\frac12(\dot Q^2 + \omega^2Q^2) 的形式(这是麦克斯韦方程的结构性结论,此处不展开)。量子化后,模式的能级为 (n+12)ℏω(n + \frac12)\hbar\omega,nn 就是该模式中的光子数,每个光子带能量 ℏω\hbar\omega、动量 ℏk\hbar\mathbf k。

经典的困难。若按经典均分,每个模式(两个平方项)的平均能量为 kBTk_{\mathrm B}T,则单位体积、单位频率间隔的辐射能量为

u(ω)=ω2π2c3 kBT(P3.12)u(\omega) = \frac{\omega^2}{\pi^2c^3}\,k_{\mathrm B}T \tag{P3.12}

(瑞利–金斯公式),对全部频率积分发散——"紫外灾难"(§9.10)。第15章用量子谐振子的平均能量 (6.14) 代替 kBTk_{\mathrm B}T,得到普朗克公式。

§P3.5光的散射与密度涨落

物理输入

偶极矩随时间变化的偶极子向外辐射电磁波,辐射功率为 P=d¨2/(6πε0c3)P = \ddot d^2/(6\pi\varepsilon_0c^3)(拉莫尔公式)。

单个分子的散射。频率为 ω\omega 的光照射分子,感生偶极矩 d=αE0cos⁡ωtd = \alpha E_0\cos\omega t,d¨=−ω2αE0cos⁡ωt\ddot d = -\omega^2\alpha E_0\cos\omega t。对时间平均(⟨cos⁡2⟩=12\langle\cos^2\rangle = \frac12):

⟨P⟩=α2E02 ω412πε0c3(P3.13)\langle P\rangle = \frac{\alpha^2E_0^2\,\omega^4}{12\pi\varepsilon_0c^3} \tag{P3.13}

散射正比于 ω4\omega^4(瑞利散射):蓝光比红光散射得强得多,所以天空是蓝色的,而日落时的阳光偏红。

介质的散射:只有涨落才散射(简化论证)。把介质分成许多比光波长小得多的小格子,第 jj 个格子含 NjN_j 个分子,它向某个侧向方向散射的电场正比于 NjN_j,相位由格子的位置决定。若每个格子的分子数都等于平均值 Nˉ\bar N,各格子的散射波在侧向上相干抵消——完全均匀的介质不向侧面散射光(光在玻璃中直线传播就是这个原因)。写 Nj=Nˉ+δNjN_j = \bar N + \delta N_j,平均部分抵消后只剩涨落部分;不同格子的涨落相互独立(格子大于涨落的关联长度时),它们的散射强度直接相加:

I侧向∝∑j⟨δNj2⟩,⟨δN2⟩=Nˉ n kBTκT(P3.14)I_{\text{侧向}}\propto\sum_j\langle\delta N_j^2\rangle,\qquad \langle\delta N^2\rangle = \bar N\,n\,k_{\mathrm B}T\kappa_T \tag{P3.14}

第二式由 (8.8) 得到。对理想气体 κT=1/p\kappa_T = 1/p,⟨δN2⟩=Nˉ\langle\delta N^2\rangle = \bar N,散射强度等于各分子独立散射之和;在临界点附近 κT→∞\kappa_T\to\infty,散射急剧增强,介质变得浑浊——这就是临界乳光(爱因斯坦 1910 年的解释;第19章)。光散射因此是测量密度涨落与压缩系数的实验手段。

§P3.6电路、电阻与热噪声

物理输入

电容器的能量 Q2/2CQ^2/2C;电感的能量 LI2/2LI^2/2(I=dQ/dtI = dQ/dt);电阻满足 V=IRV = IR,消耗功率 I2RI^2R。

LC 电路是一个谐振子。回路方程 L dI/dt+Q/C=0L\,dI/dt + Q/C = 0,即 LQ¨+Q/C=0L\ddot Q + Q/C = 0,与 mx¨+kx=0m\ddot x + kx = 0 形式相同(Q↔xQ\leftrightarrow x,L↔mL\leftrightarrow m,1/C↔k1/C\leftrightarrow k),频率 ω0=1/LC\omega_0 = 1/\sqrt{LC},能量 12LI2+Q2/2C\frac12LI^2 + Q^2/2C。在温度 TT 下,电容器上的电荷有热涨落;由均分定理 (5.16),⟨Q2⟩/2C=12kBT\langle Q^2\rangle/2C = \frac12k_{\mathrm B}T,所以

⟨V2⟩=⟨Q2⟩C2=kBTC(P3.15)\langle V^2\rangle = \frac{\langle Q^2\rangle}{C^2} = \frac{k_{\mathrm B}T}{C} \tag{P3.15}

例如 C=1C = 1 pF、300 K 时,电压涨落的方均根约为 64 μV。这个"kBT/Ck_{\mathrm B}T/C 噪声"与电路中电阻的大小无关。

电阻的热噪声(约翰逊–奈奎斯特噪声,启发式推导)。把电阻 RR 与电容 CC 并联。把电阻看成一个理想电阻串联一个随机的噪声电压源,其噪声在各频率上均匀分布("白噪声"),单位频率间隔(f≥0f\ge0)内的方均电压记为 SVS_V。电容上的电压通过 RCRC 低通电路响应噪声源:频率为 ff 的分量在电容上的振幅比为 1/∣1+i2πfRC∣1/\lvert1 + i2\pi fRC\rvert(用复阻抗,提示 A13)。各频率分量的贡献相加(提示 A14):

⟨VC2⟩=∫0∞SV df1+(2πfRC)2=SV2πRC⋅π2=SV4RC\langle V_C^2\rangle = \int_0^\infty\frac{S_V\,df}{1 + (2\pi fRC)^2} = \frac{S_V}{2\pi RC}\cdot\frac\pi2 = \frac{S_V}{4RC}

(令 u=2πfRCu = 2\pi fRC,∫0∞du/(1+u2)=π/2\int_0^\infty du/(1+u^2) = \pi/2。)平衡时它必须等于 (P3.15) 的 kBT/Ck_{\mathrm B}T/C,而且这对任何 RR 都成立,所以

SV=4kBTR(P3.16)S_V = 4k_{\mathrm B}TR \tag{P3.16}

例如 R=1 MΩR = 1\ \mathrm{M\Omega}、带宽 10 kHz、300 K 时,噪声电压的方均根为 4kBTRΔf≈13\sqrt{4k_{\mathrm B}TR\Delta f}\approx13 μV。电阻(耗散)越大,热噪声(涨落)越大:这是涨落–耗散关系的又一个例子,第24章将给出严格推导。

自测题

  1. 用 (P3.16) 估计 R=50 ΩR = 50\ \Omega、带宽 1 GHz、300 K 时的噪声电压方均根。[答:约 29 μV]
  2. 对 J=1/2J = 1/2 逐项验证 ∑mm2=J(J+1)(2J+1)/3\sum_mm^2 = J(J+1)(2J+1)/3,并验证 (P3.5) 在 J=1/2J = 1/2、gJ=2g_J = 2 时与 (6.7) 一致。
  3. 由 (P3.9) 求氩的 rmin⁡r_{\min},并说明在 r≫σr\gg\sigma 时伦纳德–琼斯势的吸引部分与 (P3.7)、(P3.8) 中哪一种力相对应。[答:rmin⁡≈0.38r_{\min}\approx0.38 nm;氩没有固有偶极矩,对应色散力 (P3.8)]