15.1

Problem 15.1

schwarzChapter 15

习题 15.1

来源: 第15章, PDF第299页


15.1 Evaluate the Casimir force using the Gaussian regulator in Eq. (15.29).

Referenced Equations:

Equation (15.29):

E(r)=12nωne(ωnπΛ)2,(15.29)E(r) = \frac{1}{2} \sum_{n} \omega_n e^{-(\frac{\omega_n}{\pi \Lambda})^2}, \tag{15.29}

习题 15.1 - 解答


物理背景与题目分析

在 1 维卡西米尔效应(Casimir effect)中,假设两个平行板(或腔体)之间的距离为 rr。对于无质量标量场,在板两端施加 Dirichlet 边界条件(即波函数在边界处为零),允许的驻波模式的波数为 kn=nπrk_n = \frac{n\pi}{r},对应的本征频率为:

ωn=nπr,n=1,2,3,\omega_n = \frac{n\pi}{r}, \quad n = 1, 2, 3, \dots

系统的零点能(真空能)形式上是一个发散的级数 E(r)=12n=1ωnE(r) = \frac{1}{2} \sum_{n=1}^{\infty} \omega_n。为了得到有物理意义的有限结果,题目要求使用式 (15.29) 给出的高斯正规化(Gaussian regulator)方法:

E(r)=12n=1ωne(ωnπΛ)2E(r) = \frac{1}{2} \sum_{n=1}^{\infty} \omega_n e^{-\left(\frac{\omega_n}{\pi \Lambda}\right)^2}

其中 Λ\Lambda 是具有能量量纲的截断参数(cutoff)。当 Λ\Lambda \to \infty 时,我们将提取出与 rr 相关的有限能量部分,并据此求出卡西米尔力。

推导过程

ωn=nπr\omega_n = \frac{n\pi}{r} 代入正规化的能量表达式中:

E(r)=12n=1nπre(nπ/rπΛ)2=π2rn=1nen2r2Λ2E(r) = \frac{1}{2} \sum_{n=1}^{\infty} \frac{n\pi}{r} e^{-\left(\frac{n\pi / r}{\pi \Lambda}\right)^2} = \frac{\pi}{2r} \sum_{n=1}^{\infty} n e^{-\frac{n^2}{r^2 \Lambda^2}}

为了方便计算,定义无量纲的小参数 α=1r2Λ2\alpha = \frac{1}{r^2 \Lambda^2},则能量表达式可以简写为:

E(r)=π2rn=1f(n)E(r) = \frac{\pi}{2r} \sum_{n=1}^{\infty} f(n)

其中 f(x)=xeαx2f(x) = x e^{-\alpha x^2}

为了在 α0\alpha \to 0(即 Λ\Lambda \to \infty)的极限下对该级数求和,我们使用欧拉-麦克劳林求和公式(Euler-Maclaurin formula):

n=1f(n)=0f(x)dx+f(0)+f()2+k=1B2k(2k)!(f(2k1)()f(2k1)(0))\sum_{n=1}^{\infty} f(n) = \int_{0}^{\infty} f(x) dx + \frac{f(0) + f(\infty)}{2} + \sum_{k=1}^{\infty} \frac{B_{2k}}{(2k)!} \left( f^{(2k-1)}(\infty) - f^{(2k-1)}(0) \right)

其中 B2kB_{2k} 为伯努利数,B2=16,B4=130B_2 = \frac{1}{6}, B_4 = -\frac{1}{30} 等。

我们逐项计算 f(x)f(x) 在边界处的值及其导数:

  1. 积分项
    0xeαx2dx=[12αeαx2]0=12α\int_{0}^{\infty} x e^{-\alpha x^2} dx = \left[ -\frac{1}{2\alpha} e^{-\alpha x^2} \right]_{0}^{\infty} = \frac{1}{2\alpha}
  2. 边界值: 由于 f(0)=0f(0) = 0f()=0f(\infty) = 0,函数值项为 00
  3. 导数项: 计算 f(x)f(x) 的各阶导数:
    f(x)=(12αx2)eαx2    f(0)=1f'(x) = (1 - 2\alpha x^2) e^{-\alpha x^2} \implies f'(0) = 1
    f(x)=(6αx+4α2x3)eαx2    f(0)=0f''(x) = (-6\alpha x + 4\alpha^2 x^3) e^{-\alpha x^2} \implies f''(0) = 0
    f(x)=(6α+24α2x28α3x4)eαx2    f(0)=6αf'''(x) = (-6\alpha + 24\alpha^2 x^2 - 8\alpha^3 x^4) e^{-\alpha x^2} \implies f'''(0) = -6\alpha
    xx \to \infty 时,所有阶导数均趋于 00

将上述结果代入欧拉-麦克劳林公式:

n=1neαn2=12αB22!f(0)B44!f(0)\sum_{n=1}^{\infty} n e^{-\alpha n^2} = \frac{1}{2\alpha} - \frac{B_2}{2!} f'(0) - \frac{B_4}{4!} f'''(0) - \dots
=12α1/62(1)1/3024(6α)= \frac{1}{2\alpha} - \frac{1/6}{2} (1) - \frac{-1/30}{24} (-6\alpha) - \dots
=12α112α120+O(α2)= \frac{1}{2\alpha} - \frac{1}{12} - \frac{\alpha}{120} + \mathcal{O}(\alpha^2)

α=1r2Λ2\alpha = \frac{1}{r^2 \Lambda^2} 代回能量表达式:

E(r)=π2r(r2Λ221121120r2Λ2+)E(r) = \frac{\pi}{2r} \left( \frac{r^2 \Lambda^2}{2} - \frac{1}{12} - \frac{1}{120 r^2 \Lambda^2} + \dots \right)
E(r)=πrΛ24π24rπ240r3Λ2+E(r) = \frac{\pi r \Lambda^2}{4} - \frac{\pi}{24r} - \frac{\pi}{240 r^3 \Lambda^2} + \dots

物理意义与最终答案

Λ\Lambda \to \infty 的极限下,能量表达式包含三部分:

  1. 发散项 πrΛ24\frac{\pi r \Lambda^2}{4}:该项与板间距离 rr 成正比,代表了真空的体能量密度(bulk vacuum energy density)。在计算卡西米尔力时,由于板外侧同样存在真空,移动板仅仅是将体积从内部转移到外部,总的体能量保持不变。因此,该项对可观测的卡西米尔力没有贡献(或者可以通过重整化减去)。
  2. 有限项 π24r-\frac{\pi}{24r}:这是与 rr 成反比的有限能量,即重整化后的卡西米尔能量 ECasimir(r)E_{\text{Casimir}}(r)
  3. 消失项 π240r3Λ2-\frac{\pi}{240 r^3 \Lambda^2}:当 Λ\Lambda \to \infty 时,该项趋于零。

因此,重整化后的卡西米尔能量为:

ECasimir(r)=π24rE_{\text{Casimir}}(r) = -\frac{\pi}{24r}

卡西米尔力 FF 定义为能量对距离 rr 的负导数(即板受到的吸引力):

F=ECasimir(r)r=r(π24r)=π24r2F = -\frac{\partial E_{\text{Casimir}}(r)}{\partial r} = -\frac{\partial}{\partial r} \left( -\frac{\pi}{24r} \right) = -\frac{\pi}{24r^2}

负号表示该力是吸引力。

F=π24r2\boxed{F = -\frac{\pi}{24r^2}}
15.2

Problem 15.2

schwarzChapter 15

习题 15.2

来源: 第15章, PDF第299页


15.2 Show that the Casimir force from the vacuum energy of fermions has the opposite sign than from bosons.

习题 15.2 - 解答


卡西米尔效应(Casimir effect)源于量子场的真空零点能随空间边界条件(如两平行极板间距 aa)的变化。要证明费米子真空能产生的卡西米尔力与玻色子符号相反,核心在于由自旋-统计定理决定的量子化条件:玻色子满足对易关系,而费米子满足反对易关系。

1. 玻色子的真空零点能

对于自由的复玻色场(如带电标量场),其哈密顿量在动量空间可以展开为粒子与反粒子的产生和湮灭算符: HB=kωk(akak+bkbk)H_B = \sum_{k} \hbar \omega_k (a_k^\dagger a_k + b_k b_k^\dagger) 根据玻色子的对易关系 [bk,bk]=δk,k[b_k, b_{k'}^\dagger] = \delta_{k,k'},我们有 bkbk=bkbk+1b_k b_k^\dagger = b_k^\dagger b_k + 1。将其代入哈密顿量得到: HB=kωk(akak+bkbk+1)H_B = \sum_{k} \hbar \omega_k \left( a_k^\dagger a_k + b_k^\dagger b_k + 1 \right) 作用于真空态 0|0\rangle(定义为 ak0=bk0=0a_k|0\rangle = b_k|0\rangle = 0),得到玻色场的真空零点能为正: E0B=0HB0=kωkE_0^B = \langle 0 | H_B | 0 \rangle = \sum_{k} \hbar \omega_k (注:若为实玻色场如电磁场,每个自由度贡献的零点能为 +12ωk+\frac{1}{2}\hbar\omega_k。无论哪种情况,玻色子的真空能均为正值)。

2. 费米子的真空零点能

对于自由的费米场(如狄拉克场),其哈密顿量同样展开为粒子与反粒子的算符。为了保证能量本征值为正,哈密顿量形式为: HF=k,sωk(ak,sak,sbk,sbk,s)H_F = \sum_{k, s} \hbar \omega_k (a_{k,s}^\dagger a_{k,s} - b_{k,s} b_{k,s}^\dagger) 根据费米子的反对易关系 {bk,s,bk,s}=δk,kδs,s\{b_{k,s}, b_{k',s'}^\dagger\} = \delta_{k,k'}\delta_{s,s'},我们有 bk,sbk,s=1bk,sbk,sb_{k,s} b_{k,s}^\dagger = 1 - b_{k,s}^\dagger b_{k,s}。将其代入哈密顿量得到: HF=k,sωk(ak,sak,s+bk,sbk,s1)H_F = \sum_{k, s} \hbar \omega_k \left( a_{k,s}^\dagger a_{k,s} + b_{k,s}^\dagger b_{k,s} - 1 \right) 作用于真空态 0|0\rangle(定义为 ak,s0=bk,s0=0a_{k,s}|0\rangle = b_{k,s}|0\rangle = 0),得到费米场的真空零点能为负: E0F=0HF0=k,sωkE_0^F = \langle 0 | H_F | 0 \rangle = -\sum_{k, s} \hbar \omega_k (注:在狄拉克海图像中,这对应于所有负能级被填满所产生的无限大负能量。按自由度平均,每个费米子自由度贡献的零点能为 12ωk-\frac{1}{2}\hbar\omega_k)。

3. 卡西米尔能量与作用力的符号比较

设两平行极板间距为 aa。由于边界条件的限制,动量 kk 成为依赖于 aa 的离散模式 kn(a)k_n(a)。系统的卡西米尔能量 E(a)\mathcal{E}(a) 定义为极板存在时的真空能与极板移至无穷远时(连续谱)的真空能之差(经过适当的紫外正规化)。

对于玻色子,其正规化后的卡西米尔能量为: EB(a)=[kn12ωkn(a)dN(k)12ωk]reg\mathcal{E}_B(a) = \left[ \sum_{k_n} \frac{1}{2} \hbar \omega_{k_n}(a) - \int dN(k) \frac{1}{2} \hbar \omega_k \right]_{\text{reg}}

对于费米子,假设其具有与玻色子相同的边界条件和能谱结构(或仅相差一个正的简并度因子 gg),由于其零点能带有整体负号,其卡西米尔能量为: EF(a)=g[kn12ωkn(a)dN(k)(12ωk)]reg=gEB(a)\mathcal{E}_F(a) = g \left[ -\sum_{k_n} \frac{1}{2} \hbar \omega_{k_n}(a) - \int dN(k) \left(-\frac{1}{2} \hbar \omega_k\right) \right]_{\text{reg}} = -g \mathcal{E}_B(a)

卡西米尔力由真空能量对极板间距 aa 的负导数给出: F=E(a)aF = -\frac{\partial \mathcal{E}(a)}{\partial a}

因此,费米子与玻色子产生的卡西米尔力之间的关系为: FF=EF(a)a=gEB(a)a=gFBF_F = -\frac{\partial \mathcal{E}_F(a)}{\partial a} = g \frac{\partial \mathcal{E}_B(a)}{\partial a} = -g F_B

由于简并度 g>0g > 0,费米子产生的卡西米尔力 FFF_F 与玻色子产生的卡西米尔力 FBF_B 符号严格相反。

FFFB\boxed{F_F \propto -F_B}

15.3

Problem 15.3

schwarzChapter 15

习题 15.3

来源: 第15章, PDF第299页


15.3 It has been proposed that geckos use the Casimir force to climb walls. It is known that geckos do not use suction (like salamanders) or capillary adhesion (like some frogs). A gecko's foot is covered in a million tiny hairs called setae, which terminate in spatula-shaped structures around 0.5 μm0.5\ \mu\text{m} wide. Use dimensional analysis and the form of the Casimir force to decide whether you think this could be possible.

习题 15.3 - 解答


题目分析与物理解释

本题要求通过量纲分析和 Casimir 力(卡西米尔力)的形式,评估壁虎利用该力攀爬墙壁的假说是否具有物理可行性。

Casimir 力是由于量子真空涨落引起的一种宏观量子效应。当两个不带电的表面靠得非常近时,它们之间会产生一种吸引力。壁虎脚上长有数以百万计的刚毛(setae),末端分叉为极小的匙突(spatulae),这种微观结构使得壁虎脚能够与墙壁表面达到极小的间距 dd,从而可能激发显著的 Casimir 力。

第一步:Casimir 力的量纲分析

设两平行表面之间的 Casimir 压强为 PP(单位面积上的力)。由于该力源于相对论性量子真空涨落,它必然依赖于约化普朗克常数 \hbar、光速 cc 以及两表面之间的距离 dd

我们写出各物理量的量纲:

  • 压强 [P]=ML1T2[P] = \text{M} \text{L}^{-1} \text{T}^{-2}
  • 约化普朗克常数 []=ML2T1[\hbar] = \text{M} \text{L}^2 \text{T}^{-1}
  • 光速 [c]=LT1[c] = \text{L} \text{T}^{-1}
  • 距离 [d]=L[d] = \text{L}

假设压强的表达式形式为 PαcβdγP \propto \hbar^\alpha c^\beta d^\gamma,代入量纲有: ML1T2=(ML2T1)α(LT1)β(L)γ=MαL2α+β+γTαβ\text{M} \text{L}^{-1} \text{T}^{-2} = (\text{M} \text{L}^2 \text{T}^{-1})^\alpha (\text{L} \text{T}^{-1})^\beta (\text{L})^\gamma = \text{M}^\alpha \text{L}^{2\alpha+\beta+\gamma} \text{T}^{-\alpha-\beta}

比较两边的指数,得到线性方程组:

  1. 质量 M\text{M}: α=1\alpha = 1
  2. 时间 T\text{T}: αβ=2    β=1-\alpha - \beta = -2 \implies \beta = 1
  3. 长度 L\text{L}: 2α+β+γ=1    2(1)+1+γ=1    γ=42\alpha + \beta + \gamma = -1 \implies 2(1) + 1 + \gamma = -1 \implies \gamma = -4

因此,Casimir 压强的形式为: Pcd4P \sim \frac{\hbar c}{d^4} (注:对于理想导电平行板,严格的理论推导给出的系数为 π2240\frac{\pi^2}{240},即 P=π2c240d4P = \frac{\pi^2 \hbar c}{240 d^4}。在量纲估算中,我们主要关注其对距离 dd 的极度敏感性 1/d41/d^4)。

总的 Casimir 力 FCF_C 为压强乘以有效接触面积 AAFCcAd4F_C \sim \frac{\hbar c A}{d^4}

第二步:参数估算

为了判断可行性,我们需要估算壁虎脚与墙壁的有效接触面积 AA 以及壁虎的重力 WW

  1. 有效接触面积 AA: 已知一只脚上有 N106N \sim 10^6 根刚毛,匙突的宽度 w0.5 μm=5×107 mw \sim 0.5\ \mu\text{m} = 5 \times 10^{-7}\ \text{m}。 假设每个匙突的接触面积为 A1w2=(5×107 m)2=2.5×1013 m2A_1 \sim w^2 = (5 \times 10^{-7}\ \text{m})^2 = 2.5 \times 10^{-13}\ \text{m}^2。 单只脚的总有效接触面积为: A=NA1106×2.5×1013 m2=2.5×107 m2A = N \cdot A_1 \sim 10^6 \times 2.5 \times 10^{-13}\ \text{m}^2 = 2.5 \times 10^{-7}\ \text{m}^2
  2. 壁虎的重力 WW: 一只典型壁虎的质量约为 m50 g=0.05 kgm \sim 50\ \text{g} = 0.05\ \text{kg}。 其重力为: W=mg0.05 kg×9.8 m/s20.5 NW = mg \approx 0.05\ \text{kg} \times 9.8\ \text{m/s}^2 \approx 0.5\ \text{N}
  3. 接触距离 dd: 分子间作用力(包括 Casimir 力和范德华力)在纳米尺度才显著。由于匙突极其微小且柔软,能极好地贴合表面粗糙度,我们合理假设接触距离在 d10 nm=108 md \sim 10\ \text{nm} = 10^{-8}\ \text{m} 的量级。

第三步:数量级计算与可行性判断

将上述参数代入量纲分析得到的 Casimir 力公式中(取 1.05×1034 Js\hbar \approx 1.05 \times 10^{-34}\ \text{J}\cdot\text{s}c3×108 m/sc \approx 3 \times 10^8\ \text{m/s}): c3.15×1026 Jm\hbar c \approx 3.15 \times 10^{-26}\ \text{J}\cdot\text{m}

单只脚产生的 Casimir 力数量级为: FC(3.15×1026 Jm)×(2.5×107 m2)(108 m)4F_C \sim \frac{(3.15 \times 10^{-26}\ \text{J}\cdot\text{m}) \times (2.5 \times 10^{-7}\ \text{m}^2)}{(10^{-8}\ \text{m})^4} FC7.875×10331032 N0.79 NF_C \sim \frac{7.875 \times 10^{-33}}{10^{-32}}\ \text{N} \approx 0.79\ \text{N}

结论分析

计算结果表明,在间距 d10 nmd \sim 10\ \text{nm} 时,单只脚产生的 Casimir 力估算值(0.79 N\sim 0.79\ \text{N})已经大于或等于整只壁虎的重力(0.5 N\sim 0.5\ \text{N})。如果考虑壁虎有四只脚,或者间距 dd 进一步缩小(例如 d5 nmd \sim 5\ \text{nm} 时,由于 1/d41/d^4 的关系,力将增大 16 倍),该吸引力将极其巨大。

虽然在真实的凝聚态物理中,极短距离(d<10 nmd < 10\ \text{nm})下起主导作用的通常是非推迟的范德华力(van der Waals force,与 1/d31/d^3 成正比),但仅从题目要求的 Casimir 力(推迟范德华力)形式和量纲分析来看,其提供的宏观力在数量级上完全足以支撑壁虎的体重。

Yes, it is possible. Dimensional analysis yields FCcAd4. Given A2.5×107 m2 and assuming a distance d10 nm, the Casimir force is 0.8 N per foot, which is sufficient to support a typical gecko’s weight of 0.5 N.\boxed{\text{Yes, it is possible. Dimensional analysis yields } F_C \sim \frac{\hbar c A}{d^4}. \text{ Given } A \sim 2.5 \times 10^{-7}\ \text{m}^2 \text{ and assuming a distance } d \sim 10\ \text{nm}, \text{ the Casimir force is } \sim 0.8\ \text{N per foot, which is sufficient to support a typical gecko's weight of } \sim 0.5\ \text{N.}}

15.4

Problem 15.4

schwarzChapter 15

习题 15.4

来源: 第15章, PDF第299页


15.4 The vacuum energy of massive particles also contributes to the Casimir force. Before doing the calculation, how do you expect the Casimir force to depend on mass? Now do the calculation and see if you are correct (use any approximations you want – this problem is challenging).

习题 15.4 - 解答


1. 物理预期 (Physical Expectation)

在进行具体计算之前,我们可以通过量子场论的物理直觉来预期卡西米尔力对质量 mm 的依赖关系。

卡西米尔效应源于真空中虚粒子的量子涨落。对于质量为 mm 的标量场,产生一对虚粒子至少需要 2m2m 的能量。根据海森堡不确定性原理 ΔEΔt1\Delta E \Delta t \sim 1,这些虚粒子在真空中的存活时间约为 Δt1/m\Delta t \sim 1/m(自然单位制 =c=1\hbar=c=1),其能够传播的特征距离即为康普顿波长 λc1/m\lambda_c \sim 1/m

当极板间距 aa 远大于康普顿波长时(即大质量极限 ma1ma \gg 1),虚粒子在存活期间几乎无法在两个极板之间完成往返传播(往返距离为 2a2a)。因此,能够“感知”到边界条件并对卡西米尔力作出贡献的真空涨落会被指数级压低。我们预期,在大质量极限下,卡西米尔力将随质量呈指数衰减,其压低因子应正比于 e2mae^{-2ma}


2. 计算过程 (Calculation)

假设我们在三维空间中考虑两个面积为 AA、间距为 aa 的平行极板(位于 z=0z=0z=az=a),场在极板上满足 Dirichlet 边界条件 ϕ=0\phi=0

允许的动量模式为 kz=nπak_z = \frac{n\pi}{a}n=1,2,3,n=1, 2, 3, \dots),横向动量 k\mathbf{k}_\perp 是连续的。系统的零点能(真空能)为: E0=A2d2k(2π)2n=1k2+(nπa)2+m2E_0 = \frac{A}{2} \int \frac{d^2k_\perp}{(2\pi)^2} \sum_{n=1}^\infty \sqrt{k_\perp^2 + \left(\frac{n\pi}{a}\right)^2 + m^2}

步骤 2.1:横向动量积分的维度正规化 (Dimensional Regularization)

为了处理横向动量的紫外发散,我们将其推广到 dd 维空间进行维度正规化。利用标准积分公式: ddk(2π)d(k2+M2)s=1(4π)d/2Γ(sd/2)Γ(s)(M2)s+d/2\int \frac{d^d k_\perp}{(2\pi)^d} (k_\perp^2 + M^2)^s = \frac{1}{(4\pi)^{d/2}} \frac{\Gamma(-s - d/2)}{\Gamma(-s)} (M^2)^{s + d/2} 代入 d=2d=2s=1/2s=1/2,并令 Mn2=(nπ/a)2+m2M_n^2 = (n\pi/a)^2 + m^2。由于 Γ(1/2)=2π\Gamma(-1/2) = -2\sqrt{\pi}Γ(3/2)=43π\Gamma(-3/2) = \frac{4}{3}\sqrt{\pi},我们得到: d2k(2π)2(k2+Mn2)1/2=14πΓ(3/2)Γ(1/2)Mn3=16πMn3\int \frac{d^2 k_\perp}{(2\pi)^2} (k_\perp^2 + M_n^2)^{1/2} = \frac{1}{4\pi} \frac{\Gamma(-3/2)}{\Gamma(-1/2)} M_n^3 = -\frac{1}{6\pi} M_n^3 因此,单位面积的真空能为: E=E0A=112πn=1[(nπa)2+m2]3/2\mathcal{E} = \frac{E_0}{A} = -\frac{1}{12\pi} \sum_{n=1}^\infty \left[ \left(\frac{n\pi}{a}\right)^2 + m^2 \right]^{3/2}

步骤 2.2:使用 Abel-Plana 公式提取有限部分

上述级数仍然发散,我们需要提取依赖于极板间距 aa 的物理有限部分。应用 Abel-Plana 公式: n=1f(n)=0f(x)dx12f(0)+i0f(it)f(it)e2πt1dt\sum_{n=1}^\infty f(n) = \int_0^\infty f(x) dx - \frac{1}{2}f(0) + i \int_0^\infty \frac{f(it) - f(-it)}{e^{2\pi t} - 1} dt 对于 f(z)=[(zπa)2+m2]3/2f(z) = \left[ \left(\frac{z\pi}{a}\right)^2 + m^2 \right]^{3/2}

  1. 第一项 0f(x)dx\int_0^\infty f(x) dx 对应无极板时的连续真空能,通过重整化扣除。
  2. 第二项 12f(0)=12m3-\frac{1}{2}f(0) = -\frac{1}{2}m^3 与间距 aa 无关,对卡西米尔力没有贡献。
  3. 第三项给出了有限的卡西米尔能。

考虑虚轴上的取值 z=±itz = \pm it。当 t>maπt > \frac{ma}{\pi} 时,根号内部为负。根据解析延拓的 iϵi\epsilon 处方(从右半平面逼近虚轴),我们有: f(it)=(m2(tπa)2+iϵ)3/2=i((tπa)2m2)3/2f(it) = \left( m^2 - \left(\frac{t\pi}{a}\right)^2 + i\epsilon \right)^{3/2} = -i \left( \left(\frac{t\pi}{a}\right)^2 - m^2 \right)^{3/2} f(it)=(m2(tπa)2iϵ)3/2=i((tπa)2m2)3/2f(-it) = \left( m^2 - \left(\frac{t\pi}{a}\right)^2 - i\epsilon \right)^{3/2} = i \left( \left(\frac{t\pi}{a}\right)^2 - m^2 \right)^{3/2} 因此,当 t>maπt > \frac{ma}{\pi} 时,f(it)f(it)=2i((tπa)2m2)3/2f(it) - f(-it) = -2i \left( \left(\frac{t\pi}{a}\right)^2 - m^2 \right)^{3/2};当 t<maπt < \frac{ma}{\pi} 时,差值为零。

代入积分并作变量代换 u=tπau = \frac{t\pi}{a}dt=aπdudt = \frac{a}{\pi} du),单位面积的卡西米尔能为: ΔE=112π(ima/π2i((tπ/a)2m2)3/2e2πt1dt)=a6π2m(u2m2)3/2e2au1du\Delta \mathcal{E} = -\frac{1}{12\pi} \left( i \int_{ma/\pi}^\infty \frac{-2i \left( (t\pi/a)^2 - m^2 \right)^{3/2}}{e^{2\pi t} - 1} dt \right) = -\frac{a}{6\pi^2} \int_m^\infty \frac{(u^2 - m^2)^{3/2}}{e^{2au} - 1} du

步骤 2.3:大质量极限下的渐近展开 (ma1ma \gg 1)

为了计算积分,作变量代换 x=2aux = 2auΔE=196π2a32am(x24a2m2)3/2ex1dx\Delta \mathcal{E} = -\frac{1}{96\pi^2 a^3} \int_{2am}^\infty \frac{(x^2 - 4a^2m^2)^{3/2}}{e^x - 1} dx 在大质量极限 ma1ma \gg 1 下,积分主要由下限 x2amx \approx 2am 附近贡献。令 x=2am+yx = 2am + y,其中 yO(1)y \sim \mathcal{O}(1)。我们可以作如下近似: x24a2m2=(2am+y)24a2m24amyx^2 - 4a^2m^2 = (2am+y)^2 - 4a^2m^2 \approx 4amy ex1e2am+ye^x - 1 \approx e^{2am+y} 将近似代入积分中(这等价于拉普拉斯方法): ΔE196π2a3e2am0(4amy)3/2eydy=8a3/2m3/296π2a3e2am0y3/2eydy\Delta \mathcal{E} \approx -\frac{1}{96\pi^2 a^3} e^{-2am} \int_0^\infty (4amy)^{3/2} e^{-y} dy = -\frac{8 a^{3/2} m^{3/2}}{96\pi^2 a^3} e^{-2am} \int_0^\infty y^{3/2} e^{-y} dy 利用 Gamma 函数 0y3/2eydy=Γ(5/2)=3π4\int_0^\infty y^{3/2} e^{-y} dy = \Gamma(5/2) = \frac{3\sqrt{\pi}}{4},得到: ΔE112π2a3/2m3/2e2am(3π4)=116π3/2m3/2a3/2e2am\Delta \mathcal{E} \approx -\frac{1}{12\pi^2 a^{3/2}} m^{3/2} e^{-2am} \left( \frac{3\sqrt{\pi}}{4} \right) = -\frac{1}{16\pi^{3/2}} \frac{m^{3/2}}{a^{3/2}} e^{-2am}

步骤 2.4:计算卡西米尔力

单位面积的卡西米尔力(压强)为能量对间距的负导数 F=ΔEaF = -\frac{\partial \Delta \mathcal{E}}{\partial a}。 对上述结果求导: F=a(116π3/2m3/2a3/2e2am)F = -\frac{\partial}{\partial a} \left( -\frac{1}{16\pi^{3/2}} m^{3/2} a^{-3/2} e^{-2am} \right) F=116π3/2m3/2[32a5/2e2am2ma3/2e2am]F = \frac{1}{16\pi^{3/2}} m^{3/2} \left[ -\frac{3}{2} a^{-5/2} e^{-2am} - 2m a^{-3/2} e^{-2am} \right]ma1ma \gg 1 的极限下,方括号中的第二项(来自指数的导数)占据绝对主导地位。保留领头项,我们得到最终的卡西米尔力: F116π3/2m3/2(2ma3/2e2am)=m5/28π3/2a3/2e2amF \approx \frac{1}{16\pi^{3/2}} m^{3/2} \left( -2m a^{-3/2} e^{-2am} \right) = -\frac{m^{5/2}}{8\pi^{3/2} a^{3/2}} e^{-2am}

负号表示该力是吸引力。计算结果完美印证了我们在第1部分中的物理预期:大质量粒子的卡西米尔力确实受到 e2mae^{-2ma} 的指数压低。

Fm5/28π3/2a3/2e2am(for ma1)\boxed{ F \approx -\frac{m^{5/2}}{8\pi^{3/2} a^{3/2}} e^{-2am} \quad (\text{for } ma \gg 1) }