跳转至

留数定理与实积分计算

本节侧重留数定理与积分应用。第一次学习时,先读求留数:方法选择与计算,掌握孤立奇点分类、实际极点阶数与六种计算方法,再进入专项习题和 44 题综合题集。

物理问题引入:如何用复数「绕道」解决算不出来的实积分?

在理论物理中,从经典电动力学辐射场计算、阻尼系统受迫共振响应,到量子散射截面与量子场论单圈传播子积分,物理学家无时无刻不在面对形如下式的实数反常定积分:

I=∫−∞+∞cos⁡(kx)x2+a2dx,∫0∞sin⁡xxdx,∫0∞sin2⁡xx2dxI = int_-inf^+inf cos(k x)x^2 + a^2 dx, quad int_0^inf sin xx dx, quad int_0^inf sin^2 xx^2dx

若使用大一微积分的初等换元或分部积分,根本无法找到初等原函数进行牛顿 - 莱布尼茨公式计算. 然而,复分析提供了一种「降维打击」式的数学策略:

  1. 复化延拓:将实数轴上的被积函数延拓到复平面 z=x+iyz = x + iy;
  2. 闭合围道:给实轴积分补上一段半圆弧、扇形或矩形路径,使其成为封闭围道;
  3. 留数定理:围道内的积分完全不依赖于整条路径的繁琐细节,仅仅由被包围奇点处的局部代数系数(留数)之和乘以 2πi2pi i 直接决定!

实数领域中几乎不可逾越的积分天堑,在复平面上只需做简单的代数除法与求导即可迎刃而解.


1. 留数的诞生哲学与计算工具箱

1.1 留数(Residue) 诞生的代数直觉

设 z0z_0 为函数 f(z)f(z) 的孤立奇点.在以 z0z_0 为中心的去心圆环域内将函数作洛朗级数展开:

f(z)=⋯+b2(z−z0)2+b1z−z0+a0+a1(z−z0)+a2(z−z0)2+…f(z) = dots + b_2(z - z_0)^2 + b_1z - z_0 + a_0 + a_1(z - z_0) + a_2(z - z_0)^2 + dots

现在,我们用一个包含 z0z_0 的逆时针简单闭合曲线 CC 对整个级数进行逐项积分.回顾基础整幂积分公式:

∮C(z−z0)ndz={2πi,n=−10,n≠−1oint_C (z - z_0)^n dz = cases 2pi i, & n = -1 0, & n != -1 cases
  • 所有正整数次幂项 (z−z0)n(z-z_0)^n(n≥0n >= 0)沿闭合回路的积分全部为零(柯西积分定理);
  • 所有负二次及更高负幂次项 (z−z0)−k(z-z_0)^-k(k≥2k >= 2)沿闭回路的积分也全部为零(由高阶导数推广公式);
  • 全级数无穷多项积分之后,唯独只有负一次幂项 b1z−z0b_1/z - z_0 留下了非零的数值贡献 2πib12pi i b_1!

因为 b1b_1(习惯上也记作 a−1a_-1 或 c−1c_-1)是围道积分中唯一「遗留」下来的非零核心,物理学家与数学家称之为 留数(Residue),记为:

Res⁡[f(z),z0]=b1=a−1Res[f(z), z_0] = b_1 = a_-1

1.2 极点留数的极速计算法则

在物理实战中,无须每次都求出完整的洛朗级数,而是利用如下求导与极限法则快速求得留数:

法则 1:一阶单极点极限法与「商式法则」

若 z0z_0 为一阶单极点,直接乘以 (z−z0)(z - z_0) 消除分母后取极限:

Res⁡[f,z0]=limz→z0(z−z0)f(z)Res[f, z_0] = lim_z -> z_0 (z - z_0) f(z)
商式法则(物理中应用最广的留数公式)

若被积函数为分式 f(z)=P(z)Q(z)f(z) = P(z)Q(z),其中 P(z0)≠0P(z_0) != 0,z0z_0 为分母 Q(z)Q(z) 的一阶单零点(即 Q(z0)=0Q(z_0) = 0 且 Q′(z0)≠0Q'(z_0) != 0),利用洛必达法则:

Res⁡[f,z0]=limz→z0(z−z0)P(z)Q(z)=limz→z0P(z)⋅limz→z0z−z0Q(z)−Q(z0)=P(z0)Q′(z0)Res[f, z_0] = lim_z -> z_0 (z - z_0)P(z)Q(z) = lim_z -> z_0 P(z) * lim_z -> z_0 z - z_0Q(z) - Q(z_0) = P(z_0)Q'(z_0)

只需将极点代入分子,除以分母的导数,即可在数秒内算出留数!

法则 2:mm 阶高阶极点法则

若 z0z_0 为 mm 阶极点(m≥2m >= 2),两端乘 (z−z0)m(z-z_0)^m 并求 m−1m-1 阶导数:

Res⁡[f,z0]=1(m−1)!limz→z0dm−1dzm−1[(z−z0)mf(z)]Res[f, z_0] = 1(m - 1)! lim_z -> z_0 d^m-1dz^m-1left[(z - z_0)^m f(z)right]

法则 3:无穷远点留数与全平面留数和定理

定义扩充复平面中无穷远点 z=∞z = inf 的留数(注意围道方向关于无穷远点为 顺时针,负号来源于此外侧积分取向):

Res⁡[f(z),∞]=−12πi∮CRf(z)dz=−Res⁡[1z2f(1z),0]Res[f(z), inf] = -12pi ioint_C_R f(z)dz = -Resleft[1z^2fleft(1zright), 0right]

全平面留数和定理:若函数 f(z)f(z) 在扩充复平面 C¯C 上除有限个孤立奇点 z1,…,znz_1, dots, z_n 外处处解析,则 其在全平面所有有限奇点的留数与无穷远点留数之和严格为零:

∑k=1nRes⁡[f(z),zk]+Res⁡[f(z),∞]=0sum_k=1^n Res[f(z), z_k] + Res[f(z), inf] = 0

当围道内部包含大量繁琐奇点,而围道外部奇点极少(或只有无穷远点)时,利用此定理从外部反求可极大简化运算.


2. 留数定理(Residue Theorem)

定理(柯西留数定理):设简单闭曲线 CC 为正向闭围道,f(z)f(z) 在 CC 上全纯解析,在 CC 内部除有限个孤立奇点 z1,z2,…,znz_1, z_2, dots, z_n 外处处解析,则:

∮Cf(z)dz=2πi∑k=1nRes⁡[f(z),zk]oint_C f(z)dz = 2pi i sum_k=1^n Res[f(z), z_k]
例题:Fibonacci 生成函数在不同半径圆周上的围道积分

计算积分 IR=∮|z|=Rz1−z−z2dzI_R = oint_|z|=R z1 - z - z^2dz.
分母零点为两个单极点:z1=5−12≈0.618z_1 = 5-12 ~ 0.618,z2=−5+12≈−1.618z_2 = -5+12 ~ -1.618.
由商式法则计算两奇点的留数(P(z)=z,Q(z)=1−z−z2,Q′(z)=−1−2zP(z) = z, Q(z) = 1 - z - z^2, Q'(z) = -1 - 2z):

Res⁡[F,z1]=z−1−2z|z1=5−12−1−(5−1)=−5−125=5−510Res[F, z_1] = left.z-1 - 2zright|_z_1 = 5-12-1 - (5-1) = -5-125 = 5 - 510 Res⁡[F,z2]=z−1−2z|z2=−5+12−1+(5+1)=−5+125=−5+510Res[F, z_2] = left.z-1 - 2zright|_z_2 = -5+12-1 + (5+1) = -5+125 = -5 + 510
  1. 当 R=0.5R = 0.5 时:圆周内无任何奇点,由柯西定理:I0.5=0I_0.5 = 0;

  2. 当 R=1.0R = 1.0 时:圆周内仅包含奇点 z1z_1:

    I1.0=2πiRes⁡[F,z1]=2πi(5−510)=πi5−55I_1.0 = 2pi i Res[F, z_1] = 2pi i left(5 - 510right) = pi i 5 - 55
  3. 当 R=2.0R = 2.0 时:圆周同时包含 z1z_1 与 z2z_2:

    I2.0=2πi(Res⁡[F,z1]+Res⁡[F,z2])=2πi(5−510−5+510)=−2πiI_2.0 = 2pi i (Res[F, z_1] + Res[F, z_2]) = 2pi i left(5 - 510 - 5 + 510right) = -2pi i

可用无穷远展开 F(z)=−1/z+O(z−2)F(z)=-1/z+O(z^-2) 核对有限留数之和为 −1-1。讲义中的半径 0.6,0.7,1.6,1.70.6,0.7,1.6,1.7 分别包围零个、一个、一个、两个极点,对应计算见综合题集。圆周经过极点时,普通留数定理不适用。


3. 理论物理实积分的五大围道计算范式

根据被积函数的代数与解析特征,求解实数定积分有五大经典围道范式:

flowchart TD
    A["实积分计算需求"] --> B["三角周期积分 \int_0^{2\pi} R(\cos\theta,\sin\theta)\mathrm{d}\theta<br/>【范式一:单位圆代换】"]
    A --> C["实轴无穷有理分式 \int_{-\infty}^\infty \frac{P(x)}{Q(x)}\mathrm{d}x<br/>【范式二:上半平面大半圆弧】"]
    A --> D["傅里叶振荡型积分 \int_{-\infty}^\infty R(x)e^{ikx}\mathrm{d}x<br/>【范式三:若当引理半圆围道】"]
    A --> E["实轴含一阶/二阶奇点(Dirichlet型)<br/>【范式四:凹陷围道法 Indented Contours】"]
    A --> F["高斯傅里叶变换型积分<br/>【范式五:矩形闭合围道】"]

范式一:单位圆代换求解三角函数周期积分

适用于形如 ∫02πR(cos⁡θ,sin⁡θ)dθint_0^2pi R(costheta, sintheta)dtheta 的有理三角积分.
令 z=eiθz = e^itheta(θtheta 从 00 变到 2π2pi 正向绕单位圆周 C:|z|=1C: |z| = 1 一周):

cos⁡θ=z+z−12,sin⁡θ=z−z−12i,dθ=dzizcostheta = z + z^-12, quad sintheta = z - z^-12i, quad dtheta = dziz

将三角定积分完全化为单位圆周上的复有理分式围道积分,通过留数定理求出落在 |z|<1|z| < 1 内部的极点留数之和.

例题:计算 ∫02πdθa+cos⁡θint_0^2pi da + costheta(a>1a > 1)

代入参数化变换:

I=∮|z|=11a+z+z−12dziz=∮|z|=12i(z2+2az+1)dzI = oint_|z|=1 1a + z + z^-12 dziz = oint_|z|=1 2i(z^2 + 2az + 1)dz

分母零点为 z±=−a±a2−1z_+/- = -a +/- a^2 - 1.由于 a>1a > 1:

  • |z−|=a+a2−1>1|z_-| = a + a^2 - 1 > 1(位于圆外);
  • |z+|=a−a2−1<1|z_+| = a - a^2 - 1 < 1(唯一位于单位圆内的单极点).
    由商式法则计算 z+z_+ 处的留数:
Res=2/i2z+2a|z=z+=2/i2a2−1=1ia2−1Res = left.2/i2z + 2aright|_z = z_+ = 2/i2a^2 - 1 = 1ia^2 - 1

代入留数定理:

I=2πi⋅1ia2−1=2πa2−1I = 2pi i * 1ia^2 - 1 = 2a^2 - 1

范式二:无穷实轴有理分式积分

适用于 ∫−∞+∞P(x)Q(x)dxint_-inf^+inf P(x)/Q(x)dx,其中分母在实轴上无实根,且分母次数比分子至少高 2 次(deg⁡Q≥deg⁡P+2deg Q >= deg P + 2).
补上半平面大半圆弧 CR:z=ReiθC_R: z = R e^itheta(θ∈[0,π]theta in [0, pi]).当 R→∞R -> inf 时大圆弧积分由于 ∼1/R~ 1/R 衰减为 0.

∫−∞+∞P(x)Q(x)dx=2πi∑Im⁡zk>0Res⁡[P(z)Q(z),zk]int_-inf^+inf P(x)Q(x)dx = 2pi i sum_Im z_k > 0 Resleft[P(z)Q(z), z_kright]

范式三:傅里叶振荡积分与约当引理(Jordan's Lemma)

适用于求解频域谱响应反演积分:

I=∫−∞+∞R(x)eikxdx(k>0)I = int_-inf^+inf R(x) e^ikx ,dx quad (k > 0)

引理(约当引理/Jordan's Lemma):
设 CRC_R 为上半平面大半圆弧 z=Reiθz = R e^itheta(θ∈[0,π]theta in [0, pi]).若在此半圆弧上函数 R(z)R(z) 一致地满足当 R→∞R -> inf 时 MR=maxz∈CR|R(z)|→0M_R = max_z in C_R|R(z)| -> 0(仅需一次衰减,不要求二次衰减!),且 k>0k > 0,则:

limR→∞∫CRR(z)eikzdz=0lim_R -> inf int_C_R R(z) e^ikzdz = 0
补充说明

证明要点:在 CRC_R 上 |eikz|=|eikR(cos⁡θ+isin⁡θ)|=e−kRsin⁡θ|e^ikz| = |e^ik R(costheta + isintheta)| = e^-k Rsintheta.

|∫CRR(z)eikzdz|≤MRR∫0πe−kRsin⁡θdθ=2MRR∫0π/2e−kRsin⁡θdθleft|int_C_R R(z)e^ikzdzright| <= M_R R int_0^pi e^-k Rsinthetadtheta = 2 M_R R int_0^pi/2 e^-k Rsinthetadtheta

利用约当不等式 sin⁡θ≥2πθsintheta >= 2/pitheta:

≤2MRR∫0π/2e−kR2θπdθ=2MRR⋅π2kR(1−e−kR)=πkMR(1−e−kR)→MR→00<= 2 M_R R int_0^pi/2 e^-k R2pidtheta = 2 M_R R * 2kRleft(1 - e^-kRright) = k M_R (1 - e^-kR) M_R -> 0 0
例题:求解狄拉克滤波实积分 ∫−∞∞cos⁡kxx2+a2dxint_-inf^inf cos kx/x^2 + a^2dx(k>0,a>0k > 0, a > 0)

构造辅助积分 ∮Ceikzz2+a2dzoint_C e^ikzz^2 + a^2dz.围道由实轴 [−R,R][-R, R] 和上半平面半圆弧 CRC_R 构成.

  • 分母极点为 z=±iaz = +/- ia;仅有 z=iaz = ia 位于上半平面闭围道内部;

  • 计算单极点留数(商式法则):

    Res⁡[f,ia]=eikz2z|z=ia=eik(ia)2ia=e−ka2iaRes[f, ia] = left.e^ikz2zright|_z = ia = e^ik(ia)2ia = e^-ka2ia
  • 由约当引理,半圆弧积分趋于 0.取极限 R→∞R -> inf:

    ∫−∞+∞eikxx2+a2dx=2πi⋅e−ka2ia=πae−kaint_-inf^+inf e^ikxx^2 + a^2dx = 2pi i * e^-ka2ia = ae^-ka

取其实部即得:

∫−∞+∞cos⁡kxx2+a2dx=πae−kaint_-inf^+inf cos kxx^2 + a^2dx = ae^-ka

范式四:实轴奇点与凹陷围道法(Indented Contours)

当被积函数的奇点恰好落在实数积分路径上(例如原点 x=0x = 0)时,我们必须构造 绕过该实数奇点的微小圆弧(顺时针或逆时针小半圆凹陷).

小圆弧引理(分数留数引理):
若 z0z_0 为一阶单极点,圆弧 CεC_varepsilon 的半径为 εvarepsilon,圆心角为 αalpha(顺时针方向),则当 ε→0varepsilon -> 0 时:

limε→0∫Cεf(z)dz=−αiRes⁡[f(z),z0]lim_varepsilon -> 0 int_C_varepsilon f(z)dz = -alpha i Res[f(z), z_0]

若为顺时针小半圆弧(α=πalpha = pi),其积分贡献严格等于 −πiRes⁡[f,z0]-pi i Res[f, z_0].

例题 1:Dirichlet 积分 ∫0∞sin⁡xxdx=π2int_0^inf sin x/xdx = pi/2

辅助复函数 F(z)=eizzF(z) = e^izz,在原点具有一阶极点.
设计上半平面闭合凹陷围道 ΓGamma(由 [−R,−ε][-R, -varepsilon]、小半圆弧 CεC_varepsilon(从 −ε-varepsilon 顺时针绕至 εvarepsilon)、[ε,R][varepsilon, R] 及大半圆弧 CRC_R 构成).
围道内部处处全纯,由柯西定理:

∮Γeizzdz=0oint_Gamma e^izzdz = 0
  1. 大圆弧 CRC_R 由约当引理衰减趋于 0;

  2. 小半圆弧 CεC_varepsilon 在原点展开:eizz=1+iz+…z=1z+i+…e^izz = 1 + iz + dots/z = 1/z + i + dots.参数化 z=εeiθz = varepsilon e^itheta(θtheta 从 πpi 变到 00):

    ∫Cεeizzdz=∫π01εeiθiεeiθdθ=i∫π0dθ=−iπint_C_varepsilon e^izzdz = int_pi^0 1varepsilon e^itheta ivarepsilon e^ithetadtheta = i int_pi^0 dtheta = -ipi
  3. 实轴两段合并(对负半轴作代换 x→−tx -> -t):

    ∫−R−εeixxdx+∫εReixxdx=∫εReix−e−ixxdx=2i∫εRsin⁡xxdxint_-R^-varepsilon e^ixxdx + int_varepsilon^R e^ixxdx = int_varepsilon^R e^ix - e^-ixxdx = 2iint_varepsilon^R sin xxdx

取极限 ε→0,R→∞varepsilon -> 0, R -> inf 代入柯西等式:

2i∫0∞sin⁡xxdx−iπ=0⟹∫0∞sin⁡xxdx=π22i int_0^inf sin xxdx - ipi = 0 implies int_0^inf sin xxdx = 2
例题 2:先判实际阶数,再用凹陷围道求 ∫0∞sin2⁡xx2dxint_0^inf sin^2 x/x^2dx

由偶函数对称性与二倍角降幂:

I=∫0∞sin2⁡xx2dx=14∫−∞+∞1−cos⁡2xx2dxI = int_0^inf sin^2 xx^2dx = 14int_-inf^+inf 1 - cos 2xx^2dx

构造复变辅助函数 f(z)=1−e2izz2f(z) = 1 - e^2izz^2,由于分子在 z=0z = 0 处有一阶零点,原点为 一阶极点.
在相同的上半平面凹陷围道上积分,由柯西定理 ∮Γf(z)dz=0oint_Gamma f(z)dz = 0.

  • 实轴两段实部给出 4I4I;
  • 大圆弧衰减为 0;
  • 小圆弧上展开分子:1−(1+2iz−2z2+…)z2=−2iz+2+…1 - (1 + 2iz - 2z^2 + dots)/z^2 = -2i/z + 2 + dots.
    顺时针积分为 −2i(−iπ)=−2π-2i(-ipi) = -2pi.
    代回柯西方程:
4I−2π=0⟹I=∫0∞sin2⁡xx2dx=π24I - 2pi = 0 implies I = int_0^inf sin^2 xx^2dx = 2

范式五:矩形闭合围道与高斯傅里叶变换

例题:证明高斯函数傅里叶变换的自对偶性 ∫−∞+∞e−πx2e−2πixξdx=e−πξ2int_-inf^+inf e^-pi x^2 e^-2pi i x xidx = e^-pi xi^2

设 ξ>0xi > 0.考察全纯整函数 f(z)=e−πz2f(z) = e^-pi z^2.
构造顶点分别为 −R,R,R+iξ,−R+iξ-R, R, R + ixi, -R + ixi 的矩形闭合围道 γRgamma_R.由柯西积分定理:

∮γRe−πz2dz=0oint_gamma_R e^-pi z^2dz = 0
  1. 底边(y=0y = 0):∫−RRe−πx2dx→R→∞1int_-R^R e^-pi x^2dx R -> inf 1(标准高斯积分);
  2. 左右两侧竖直边(x=±Rx = +/- R):由于被积函数含有模长衰减因子 e−π(R2−y2)≤e−πR2eπξ2e^-pi(R^2 - y^2) <= e^-pi R^2 e^pi xi^2,当 R→∞R -> inf 时指数衰减至 0;
  3. 顶边(z=x+iξ,xz = x + ixi, x 从 RR 变到 −R-R):

    ∫R−Re−π(x+iξ)2dx=−∫−RRe−π(x2−ξ2+2ixξ)dx=−eπξ2∫−RRe−πx2e−2πixξdxint_R^-R e^-pi(x + ixi)^2dx = -int_-R^R e^-pi(x^2 - xi^2 + 2ixxi)dx = -e^pixi^2int_-R^R e^-pi x^2 e^-2pi i xxidx

四边求和令 R→∞R -> inf:

1+0−eπξ2∫−∞+∞e−πx2e−2πixξdx+0=0⟹∫−∞+∞e−πx2e−2πixξdx=e−πξ21 + 0 - e^pixi^2int_-inf^+inf e^-pi x^2 e^-2pi i xxidx + 0 = 0 implies int_-inf^+inf e^-pi x^2 e^-2pi i xxidx = e^-pixi^2

讲义 3.7 的计算链:留数、边界项、实积分

下面三个例子将求留数步骤与讲义的应用连接起来。相应训练见综合题集。

二阶极点:有理积分与余弦积分

对 f(z)=z2/(z2+1)2f(z)=z^2/(z^2+1)^2,上半平面只有二阶极点 ii,且

Res⁡(f,i)=ddzz2(z+i)2|z=i=−i4.Res(f,i)=left.mathrm dmathrm dz z^2(z+i)^2right|_z=i=-dfrac i4.

大圆上 |f(z)|≤R2/(R2−1)2|f(z)|<= R^2/(R^2-1)^2,弧长为 πRpi R,故圆弧积分为 O(R−1)O(R^-1)。因此

∫−∞∞x2(x2+1)2dx=2πi(−i4)=π2.int_-inf^infx^2(x^2+1)^2 ,mathrm dx =2pi ileft(-dfrac i4right)=dfracpi2.

对余弦积分,改用 g(z)=eiz/(z2+1)2g(z)=e^iz/(z^2+1)^2。必须对指数与剩余分母一起求导:

Res⁡(g,i)=[ieiz(z+i)2−2eiz(z+i)3]z=i=−i2e.Res(g,i) =left[ie^iz(z+i)^2-2e^iz(z+i)^3right]_z=i =-i2e.

大圆弧由约当引理消失;留数定理得到复实轴积分 π/epi/e,取实部即有

∫−∞∞cos⁡x(x2+1)2dx=πe.int_-inf^infcos x(x^2+1)^2 ,mathrm dx=dfracpi e.

凹陷围道:辅助函数有极点,原正弦积分可以收敛

设 a,m>0a,m>0,考虑 I=∫0∞sin⁡(mx)/[x(x2+a2)]dxI=int_0^infsin(mx)/[x(x^2+a^2)] ,mathrm dx。原实被积函数在零点趋于 m/a2m/a^2,但辅助函数 h(z)=eimz/[z(z2+a2)]h(z)=e^imz/[z(z^2+a^2)] 在原点有单极点,留数为 1/a21/a^2。

在上半平面用从 −ε-varepsilon 到 εvarepsilon 的顺时针小半圆绕开原点,小弧贡献为 −iπ/a2-ipi/a^2。围道内只含 iaia,其留数为 −e−ma/(2a2)-e^-ma/(2a^2)。实轴两段合为 2iI2iI,大圆弧由约当引理消失,所以

2iI−iπa2=−πia2e−ma,I=π2a2(1−e−ma).2iI-ia^2=-pi ia^2e^-ma, qquad I=2a^2(1-e^-ma).

当 m→0+m->0^+ 时结果趋零;当 mm 很大时趋于 π/(2a2)pi/(2a^2)。这些极限也可用来核对小圆弧的符号。

傅里叶积分与 Yukawa 核

对实数 a≠0ane0、p∈Rpinmathbb R,应分别处理 pp 的正、负、零三种情形。p>0p>0 向上闭合,p<0p<0 向下闭合并计入顺时针负号,p=0p=0 用普通有理函数的弧长估计;三者合并为

∫−∞∞eipxx2+a2dx=π|a|e−|ap|.int_-inf^infe^ipxx^2+a^2 ,mathrm dx =dfrac|a|e^-|ap|.
三维动量积分的径向核(进阶)

讲义采用没有 (2π)−3(2pi)^-3 前因子的傅里叶约定。设 r=|r|>0,μ>0r=|mathbf r|>0,mu>0,将极轴沿 rmathbf r 取向,先完成角积分,再按球形截断取径向反常积分:

I(r)=∫R3eik⋅rk2+μ2d3k=4πr∫0∞ksin⁡(kr)k2+μ2dk=2πir∫−∞∞keikrk2+μ2dk.I(r)=int_mathbb R^3e^imathbf k*mathbf rk^2+mu^2 ,mathrm d^3k =4rint_0^infksin(kr)k^2+mu^2 ,mathrm dk =2irint_-inf^infk e^ikrk^2+mu^2 ,mathrm dk.

最后一步的全实轴积分按对称截断理解。上半平面唯一极点为 iμimu,留数为 e−μr/2e^-mu r/2;有理因子 k/(k2+μ2)=O(1/k)k/(k^2+mu^2)=O(1/k) 满足约当引理条件。因此

I(r)=2πir(2πie−μr2)=2π2re−μr.I(r)=2irleft(2pi ie^-mu r2right) =2pi^2re^-mu r.

原三维积分并非绝对收敛,不能任意交换积分顺序;这里采用上述球形截断解释,也可采用相容的正则化或分布解释。若傅里叶逆变换含 (2π)−3(2pi)^-3,则径向核为 e−μr/(4πr)e^-mu r/(4pi r)。这个指数衰减给出了 Yukawa 势的短程径向形式。


4. 理论物理中的深度连接:因果律与极点处方

在相对论与量子物理中,「结果不能先于原因发生」是物理世界的基本铁律——因果律(Causality).这一物理原理在复分析中表现为极其深刻的数学结构:

4.1 柯西主值(P.V.)与索霍茨基 - 普莱姆利(Sokhotski–Plemelj) 公式

当极点落在实轴上时,物理学家在虚轴方向加上微小的扰动 x±iϵx +/- iepsilon(ϵ→0+epsilon -> 0^+)给出绕过实轴极点的处方:

limϵ→0+1x±iϵ=P1x∓iπδ(x)lim_epsilon -> 0^+ 1x +/- iepsilon = P1x mp ipi delta(x)

其中 PP 代表柯西主值(Cauchy Principal Value),δ(x)delta(x) 是狄拉克 Delta 函数。等式在分布意义下成立,不能当作 x=0x=0 处的普通函数等式。用于具体响应函数时,主值项与 Delta 项的物理意义还取决于傅里叶约定和响应模型。

4.2 克拉默斯 - 克勒尼希(Kramers–Kronig) 色散关系

在光学与凝聚态物理中,物质的复介电常数 ε(ω)=ε1(ω)+iε2(ω)varepsilon(omega) = varepsilon_1(omega) + ivarepsilon_2(omega)。采用 χ(ω)=∫0∞χ(t)eiωtdtchi(omega)=int_0^infchi(t)e^iomega t ,mathrm dt 的约定,对稳定且满足相应收敛条件的因果响应,响应函数在上半复平面解析。若高频极限为 11,且减去该极限后的大圆弧贡献消失,再结合实时间响应的对称性,可由柯西积分公式得到:

ε1(ω)−1=2πP∫0∞ω′ε2(ω′)ω′2−ω2dω′varepsilon_1(omega) - 1 = 2piPint_0^inf omega' varepsilon_2(omega')omega'^2 - omega^2domega'

这说明实部与虚部由全频段的积分关系联系起来。某一频段内可以有很小的吸收而仍有色散,不能把该关系解释为每个频率上必须同时出现非零吸收。

4.3 量子场论费曼(Feynman)iϵiepsilon 处方

在求解相对论 Klein-Gordon 波动方程自由传播子(格林函数)时,奇点落在动量实轴上(p2−m2=0p^2 - m^2 = 0).为了确保粒子向前传播、反粒子向后传播的时间因果律,费曼规定将分母微扰为:

ΔF(p)=1p2−m2+iϵDelta_F(p) = 1p^2 - m^2 + iepsilon

在度规约定 p2=(p0)2−p2p^2=(p^0)^2-mathbf p^2 下,令 Ep=p2+m2>0E_mathbf p=mathbf p^2+m^2>0,两能量极点位于 p0≈Ep−iϵ/(2Ep)p^0~ E_mathbf p-iepsilon/(2E_mathbf p) 和 p0≈−Ep+iϵ/(2Ep)p^0~-E_mathbf p+iepsilon/(2E_mathbf p)。这一处方定义时间排序的费曼传播子;推迟与超前格林函数采用各自的极点处方,不能与费曼传播子混同。


5. 学习衔接


知识小测与巩固练习

补充说明

课后核心技能自测点

  1. 用商式法则计算 f(z)=ezz4−1f(z) = e^z/z^4 - 1 在 z=iz = i 处的留数;
  2. 计算实定积分 ∫−∞∞x2(x2+1)(x2+4)dxint_-inf^inf x^2/(x^2+1)(x^2+4)dx;
  3. 说明为什么约当引理中圆弧半径仅需一次衰减(MR→0M_R -> 0),而普通有理分式需要二次衰减(MR∼1/R2M_R ~ 1/R^2)?

参考资料

本节奇点、留数定义与计算例题对照用户提供的《3.5 奇点分类、3.6 留数定理发现之旅与 3.7 留数定理》,即仓库 References/数学物理方法课件/mmplecture 3.5-3.7.md。求留数的局部方法和扩充复平面的符号约定见方法选择与计算,讲义思考题的训练与解析见综合题集。

本页小测与复习



评论