9.1

Problem 9.1

peskinChapter 9

习题 9.1

来源: 第9章, PDF第312页


9.1 Scalar QED. This problem concerns the theory of a complex scalar field ϕ\phi interacting with the electromagnetic field AμA^\mu. The Lagrangian is

L=14Fμν2+(Dμϕ)(Dμϕ)m2ϕϕ,\mathcal{L} = -\frac{1}{4} F_{\mu\nu}^2 + (D_\mu \phi)^* (D^\mu \phi) - m^2 \phi^* \phi,

where Dμ=μ+ieAμD_\mu = \partial _\mu + ieA_\mu is the usual gauge-covariant derivative.

(a) Use the functional method of Section 9.2 to show that the propagator of the complex scalar field is the same as that of a real field:

Feynman propagator diagram for a complex scalar field with momentum p.
=ip2m2+iϵ.= \frac{i}{p^2 - m^2 + i\epsilon}.

Also derive the Feynman rules for the interactions between photons and scalar particles; you should find

Feynman diagrams for scalar QED interactions: a three-point vertex representing a photon interacting with a scalar and an anti-scalar, and a four-point vertex representing a two-photon two-scalar interaction.

(b) Compute, to lowest order, the differential cross section for e+eϕϕe^+ e^- \rightarrow \phi \phi^*. Ignore the electron mass (but not the scalar particle's mass), and average over the electron and positron polarizations. Find the asymptotic angular dependence and total cross section. Compare your results to the corresponding formulae for e+eμ+μe^+ e^- \rightarrow \mu^+ \mu^-.

(c) Compute the contribution of the charged scalar to the photon vacuum polarization, using dimensional regularization. Note that there are two diagrams. To put the answer into the expected form,

Πμν(q2)=(gμνq2qμqν)Π(q2),\Pi^{\mu\nu}(q^2) = (g^{\mu\nu} q^2 - q^\mu q^\nu) \Pi(q^2),

it is useful to add the two diagrams at the beginning, putting both terms over a common denominator before introducing a Feynman parameter. Show that, for q2m2-q^2 \gg m^2, the charged boson contribution to Π(q2)\Pi(q^2) is exactly 1/4 that of a virtual electron-positron pair.

习题 9.1 - 解答


习题 9.1 分析与解答

(a) 标量场的传播子与费曼规则

标量 QED 的拉格朗日量为: L=14Fμν2+(Dμϕ)(Dμϕ)m2ϕϕ\mathcal{L} = -\frac{1}{4} F_{\mu\nu}^2 + (D_\mu \phi)^* (D^\mu \phi) - m^2 \phi^* \phi 其中协变导数 Dμ=μ+ieAμD_\mu = \partial_\mu + ieA_\mu。展开标量场的动能项: (Dμϕ)(Dμϕ)=(μieAμ)ϕ(μ+ieAμ)ϕ=μϕμϕ+ieAμ(ϕμϕϕμϕ)+e2AμAμϕϕ(D_\mu \phi)^* (D^\mu \phi) = (\partial_\mu - ieA_\mu)\phi^* (\partial^\mu + ieA^\mu)\phi = \partial_\mu \phi^* \partial^\mu \phi + ieA_\mu (\phi^* \partial^\mu \phi - \phi \partial^\mu \phi^*) + e^2 A_\mu A^\mu \phi^* \phi 自由标量场的拉格朗日量为 L0=μϕμϕm2ϕϕ\mathcal{L}_0 = \partial_\mu \phi^* \partial^\mu \phi - m^2 \phi^* \phi。在作用量中分部积分,可将其写为: S0=d4xϕ(2m2)ϕS_0 = \int d^4x \, \phi^* (-\partial^2 - m^2) \phi 复标量场的传播子是算符 (2m2)(-\partial^2 - m^2) 的逆。在动量空间中,μipμ\partial_\mu \to -ip_\mu,算符变为 p2m2p^2 - m^2。引入 Feynman 处方 +iϵ+i\epsilon 以处理极点,得到复标量场的传播子(与实标量场相同,但带有表示电荷流向的箭头): ip2m2+iϵ\boxed{ \frac{i}{p^2 - m^2 + i\epsilon} }

相互作用拉格朗日量为 Lint=ieAμ(ϕμϕϕμϕ)+e2AμAμϕϕ\mathcal{L}_{int} = ieA_\mu (\phi^* \partial^\mu \phi - \phi \partial^\mu \phi^*) + e^2 A_\mu A^\mu \phi^* \phi。 对于三点顶点(一个光子,一个标量,一个反标量),考虑作用量 iSintid4xieAμ(ϕμϕϕμϕ)iS_{int} \supset i \int d^4x \, ieA_\mu (\phi^* \partial^\mu \phi - \phi \partial^\mu \phi^*)。设入射标量粒子动量为 pp,出射标量粒子动量为 pp'(等价于入射反粒子动量为 p-p')。在动量空间中,导数作用于场产生动量因子:μϕipμϕ\partial^\mu \phi \to -ip^\mu \phiμϕipμϕ\partial^\mu \phi^* \to ip'^\mu \phi^*。提取 AμA_\mu 的系数并乘以 ii,得到顶点规则: i[ie(ipμ)ie(ipμ)]=ie(p+p)μi \left[ ie(-ip^\mu) - ie(ip'^\mu) \right] = \boxed{ -ie(p + p')^\mu } 对于四点顶点(两个光子,一个标量,一个反标量),考虑作用量 iSintid4xe2gμνAμAνϕϕiS_{int} \supset i \int d^4x \, e^2 g^{\mu\nu} A_\mu A_\nu \phi^* \phi。由于有两个相同的光子场 AμA_\muAνA_\nu,在泛函求导或 Wick 展开时会产生一个对称因子 2。顶点规则为: i×2×e2gμν=2ie2gμνi \times 2 \times e^2 g^{\mu\nu} = \boxed{ 2ie^2 g^{\mu\nu} }

(b) e+eϕϕe^+ e^- \rightarrow \phi \phi^* 的微分散射截面

在最低阶(树图阶),该过程通过 ss 沟道光子交换进行。设电子和正电子的动量分别为 k1,k2k_1, k_2,出射标量和反标量的动量分别为 p1,p2p_1, p_2。忽略电子质量。散射振幅为: iM=[vˉ(k2)(ieγμ)u(k1)](igμνq2)[ie(p1p2)ν]i\mathcal{M} = \left[ \bar{v}(k_2) (-ie\gamma^\mu) u(k_1) \right] \left( \frac{-ig_{\mu\nu}}{q^2} \right) \left[ -ie(p_1 - p_2)^\nu \right] 其中 q=k1+k2=p1+p2q = k_1 + k_2 = p_1 + p_2,且 q2=sq^2 = s。注意标量顶点中,出射 ϕ(p1)\phi(p_1) 对应动量 p1p_1,出射 ϕ(p2)\phi^*(p_2) 对应入射动量 p2-p_2,故顶点因子为 ie(p1+(p2))ν-ie(p_1 + (-p_2))^\nu。 对初始自旋求平均并对末态求和,振幅平方为: 14spinsM2=e44s2Tr[γμ̸k1γν̸k2](p1p2)μ(p1p2)ν\frac{1}{4} \sum_{\text{spins}} |\mathcal{M}|^2 = \frac{e^4}{4s^2} \text{Tr}[\gamma^\mu \not{k}_1 \gamma^\nu \not{k}_2] (p_1 - p_2)_\mu (p_1 - p_2)_\nu 计算迹 Tr[γμ̸k1γν̸k2]=4(k1μk2ν+k1νk2μgμνk1k2)\text{Tr}[\gamma^\mu \not{k}_1 \gamma^\nu \not{k}_2] = 4(k_1^\mu k_2^\nu + k_1^\nu k_2^\mu - g^{\mu\nu} k_1 \cdot k_2),并令 Δp=p1p2\Delta p = p_1 - p_2,得到: 14M2=e4s2[2(k1Δp)(k2Δp)(k1k2)(Δp)2]\frac{1}{4} \sum |\mathcal{M}|^2 = \frac{e^4}{s^2} \left[ 2(k_1 \cdot \Delta p)(k_2 \cdot \Delta p) - (k_1 \cdot k_2)(\Delta p)^2 \right] 在质心系中,设束流能量为 EEs=4E2s=4E^2),标量粒子速度为 β=1m2/E2\beta = \sqrt{1 - m^2/E^2}。运动学变量为: k1=(E,0,0,E)k_1 = (E, 0, 0, E), k2=(E,0,0,E)k_2 = (E, 0, 0, -E) p1=(E,p)p_1 = (E, \vec{p}), p2=(E,p)p_2 = (E, -\vec{p}),其中 p=Eβ|\vec{p}| = E\beta Δp=(0,2p)\Delta p = (0, 2\vec{p})(Δp)2=4E2β2(\Delta p)^2 = -4E^2\beta^2θ\theta 为出射标量粒子与入射电子束的夹角,则 k1Δp=2Epcosθ=2E2βcosθk_1 \cdot \Delta p = -2E|\vec{p}|\cos\theta = -2E^2\beta\cos\thetak2Δp=2E2βcosθk_2 \cdot \Delta p = 2E^2\beta\cos\thetak1k2=2E2k_1 \cdot k_2 = 2E^2。代入上式: 14M2=e416E4[2(2E2βcosθ)(2E2βcosθ)(2E2)(4E2β2)]=e416E4[8E4β2(1cos2θ)]=12e4β2sin2θ\frac{1}{4} \sum |\mathcal{M}|^2 = \frac{e^4}{16E^4} \left[ 2(-2E^2\beta\cos\theta)(2E^2\beta\cos\theta) - (2E^2)(-4E^2\beta^2) \right] = \frac{e^4}{16E^4} \left[ 8E^4\beta^2(1 - \cos^2\theta) \right] = \frac{1}{2} e^4 \beta^2 \sin^2\theta 微分散射截面为: dσdΩ=164π2spk(14M2)=β64π2s(12e4β2sin2θ)=α2β38ssin2θ\frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 s} \frac{|\vec{p}|}{|\vec{k}|} \left( \frac{1}{4} \sum |\mathcal{M}|^2 \right) = \frac{\beta}{64\pi^2 s} \left( \frac{1}{2} e^4 \beta^2 \sin^2\theta \right) = \boxed{ \frac{\alpha^2 \beta^3}{8s} \sin^2\theta } 总截面为: σ=dσdΩdΩ=2π11d(cosθ)α2β38s(1cos2θ)=πα2β33s\sigma = \int \frac{d\sigma}{d\Omega} d\Omega = 2\pi \int_{-1}^1 d(\cos\theta) \frac{\alpha^2 \beta^3}{8s} (1-\cos^2\theta) = \boxed{ \frac{\pi \alpha^2 \beta^3}{3s} } 比较: 对于 e+eμ+μe^+ e^- \rightarrow \mu^+ \mu^-,渐近微分散射截面(β1\beta \to 1)正比于 (1+cos2θ)(1+\cos^2\theta),而标量 QED 中正比于 sin2θ\sin^2\theta。这是因为标量对产生时沿束流方向的角动量投影为零,为了守恒光子的自旋角动量,必须有轨道角动量贡献,导致前向和后向散射为零。渐近总截面 σ(e+eμ+μ)4πα23s\sigma(e^+ e^- \rightarrow \mu^+ \mu^-) \to \frac{4\pi\alpha^2}{3s},标量 QED 的结果 πα23s\frac{\pi\alpha^2}{3s} 恰好是其 1/41/4

(c) 光子真空极化

标量场对光子真空极化 Πμν(q)\Pi^{\mu\nu}(q) 有两个单圈图贡献:

  1. 包含两个三点顶点的圈图: iΠ1μν(q)=d4k(2π)4(ie)(2k+q)μ(ie)(2k+q)ν((k+q)2m2)(k2m2)=e2d4k(2π)4(2k+q)μ(2k+q)ν((k+q)2m2)(k2m2)i\Pi_1^{\mu\nu}(q) = \int \frac{d^4k}{(2\pi)^4} \frac{(-ie)(2k+q)^\mu (-ie)(2k+q)^\nu}{((k+q)^2-m^2)(k^2-m^2)} = -e^2 \int \frac{d^4k}{(2\pi)^4} \frac{(2k+q)^\mu (2k+q)^\nu}{((k+q)^2-m^2)(k^2-m^2)}
  2. 包含一个四点顶点的圈图: iΠ2μν(q)=d4k(2π)42ie2gμνk2m2=2ie2gμνd4k(2π)41k2m2i\Pi_2^{\mu\nu}(q) = \int \frac{d^4k}{(2\pi)^4} \frac{2ie^2 g^{\mu\nu}}{k^2-m^2} = 2ie^2 g^{\mu\nu} \int \frac{d^4k}{(2\pi)^4} \frac{1}{k^2-m^2} 为了合并分母,利用积分变量的平移不变性,将 Π2\Pi_2 写为对称形式: 1k2m2=12(1k2m2+1(k+q)2m2)=12(k+q)2m2+k2m2((k+q)2m2)(k2m2)\int \frac{1}{k^2-m^2} = \frac{1}{2} \int \left( \frac{1}{k^2-m^2} + \frac{1}{(k+q)^2-m^2} \right) = \frac{1}{2} \int \frac{(k+q)^2-m^2 + k^2-m^2}{((k+q)^2-m^2)(k^2-m^2)} 因此 iΠ2μν(q)=ie2gμνd4k(2π)42k2+2kq+q22m2((k+q)2m2)(k2m2)i\Pi_2^{\mu\nu}(q) = ie^2 g^{\mu\nu} \int \frac{d^4k}{(2\pi)^4} \frac{2k^2+2k\cdot q+q^2-2m^2}{((k+q)^2-m^2)(k^2-m^2)}。将两项相加: iΠμν(q)=e2d4k(2π)4(2k+q)μ(2k+q)ν+gμν(2k2+2kq+q22m2)((k+q)2m2)(k2m2)i\Pi^{\mu\nu}(q) = e^2 \int \frac{d^4k}{(2\pi)^4} \frac{-(2k+q)^\mu (2k+q)^\nu + g^{\mu\nu}(2k^2+2k\cdot q+q^2-2m^2)}{((k+q)^2-m^2)(k^2-m^2)} 引入 Feynman 参数 xx,分母变为 D=(k+xq)2+x(1x)q2m2D = (k+xq)^2 + x(1-x)q^2 - m^2。令 l=k+xql = k + xq,则 D=l2ΔD = l^2 - \Delta,其中 Δ=m2x(1x)q2\Delta = m^2 - x(1-x)q^2。 将分子用 ll 表达并丢弃 ll 的奇数次项: 分子 =(2l+(12x)q)μ(2l+(12x)q)ν+gμν[2(lxq)2+2(lxq)q+q22m2]= - (2l + (1-2x)q)^\mu (2l + (1-2x)q)^\nu + g^{\mu\nu} \left[ 2(l-xq)^2 + 2(l-xq)\cdot q + q^2 - 2m^2 \right] 4lμlν(12x)2qμqν+gμν[2l2+(12x+2x2)q22m2]\to -4l^\mu l^\nu - (1-2x)^2 q^\mu q^\nu + g^{\mu\nu} \left[ 2l^2 + (1-2x+2x^2)q^2 - 2m^2 \right] 在维数正规化下,用 1dl2gμν\frac{1}{d} l^2 g^{\mu\nu} 替换 lμlνl^\mu l^\nuiΠμν(q)=e2μ4d01dxddl(2π)d1(l2Δ)2{gμν[(24d)l2+(12x+2x2)q22m2](12x)2qμqν}i\Pi^{\mu\nu}(q) = e^2 \mu^{4-d} \int_0^1 dx \int \frac{d^dl}{(2\pi)^d} \frac{1}{(l^2-\Delta)^2} \left\{ g^{\mu\nu} \left[ \left(2-\frac{4}{d}\right)l^2 + (1-2x+2x^2)q^2 - 2m^2 \right] - (1-2x)^2 q^\mu q^\nu \right\} 计算 gμνg^{\mu\nu} 项中包含 l2l^2 的积分: ddl(2π)d(24/d)l2(l2Δ)2=(24d)i(4π)d/2d2Γ(1d/2)Δd/21=i(4π)d/2(d2)Γ(1d/2)Δd/21\int \frac{d^dl}{(2\pi)^d} \frac{(2-4/d)l^2}{(l^2-\Delta)^2} = \left(2-\frac{4}{d}\right) \frac{-i}{(4\pi)^{d/2}} \frac{d}{2} \Gamma(1-d/2) \Delta^{d/2-1} = \frac{i}{(4\pi)^{d/2}} (d-2) \Gamma(1-d/2) \Delta^{d/2-1} 利用 (d2)Γ(1d/2)=2Γ(2d/2)(d-2)\Gamma(1-d/2) = -2\Gamma(2-d/2),该项变为 2Δi(4π)d/2Γ(2d/2)Δd/22-2\Delta \frac{i}{(4\pi)^{d/2}} \Gamma(2-d/2) \Delta^{d/2-2}。 将所有 gμνg^{\mu\nu} 的系数合并,提取公因子 i(4π)d/2Γ(2d/2)Δd/22\frac{i}{(4\pi)^{d/2}} \Gamma(2-d/2) \Delta^{d/2-2},括号内剩余: 2Δ+(12x+2x2)q22m2=2(m2x(1x)q2)+(12x+2x2)q22m2=(14x+4x2)q2=(12x)2q2-2\Delta + (1-2x+2x^2)q^2 - 2m^2 = -2(m^2 - x(1-x)q^2) + (1-2x+2x^2)q^2 - 2m^2 = (1-4x+4x^2)q^2 = (1-2x)^2 q^2 因此,积分结果自然给出了横向结构: iΠμν(q)=i(gμνq2qμqν)e2(4π)d/201dx(12x)2Γ(2d/2)Δ2d/2i\Pi^{\mu\nu}(q) = i (g^{\mu\nu}q^2 - q^\mu q^\nu) \frac{e^2}{(4\pi)^{d/2}} \int_0^1 dx \, (1-2x)^2 \frac{\Gamma(2-d/2)}{\Delta^{2-d/2}} 由此得到 Π(q2)\Pi(q^2)Π(q2)=α4π01dx(12x)2[2ϵln(m2x(1x)q2μ~2)]\boxed{ \Pi(q^2) = \frac{\alpha}{4\pi} \int_0^1 dx \, (1-2x)^2 \left[ \frac{2}{\epsilon} - \ln\left(\frac{m^2 - x(1-x)q^2}{\tilde{\mu}^2}\right) \right] } 对于 q2=Q2m2-q^2 = Q^2 \gg m^2,重整化后的真空极化函数 Π^(q2)=Π(q2)Π(0)\hat{\Pi}(q^2) = \Pi(q^2) - \Pi(0) 的渐近行为由对数项主导: Π^scalar(Q2)α4πln(Q2m2)01(12x)2dx=α4πln(Q2m2)[x2x2+43x3]01=α12πln(Q2m2)\hat{\Pi}_{\text{scalar}}(Q^2) \approx \frac{\alpha}{4\pi} \ln\left(\frac{Q^2}{m^2}\right) \int_0^1 (1-2x)^2 dx = \frac{\alpha}{4\pi} \ln\left(\frac{Q^2}{m^2}\right) \left[ x - 2x^2 + \frac{4}{3}x^3 \right]_0^1 = \frac{\alpha}{12\pi} \ln\left(\frac{Q^2}{m^2}\right) 对于旋量 QED(虚电子-正电子对),相应的渐近结果为: Π^spinor(Q2)2απln(Q2m2)01x(1x)dx=2απln(Q2m2)(16)=α3πln(Q2m2)\hat{\Pi}_{\text{spinor}}(Q^2) \approx \frac{2\alpha}{\pi} \ln\left(\frac{Q^2}{m^2}\right) \int_0^1 x(1-x) dx = \frac{2\alpha}{\pi} \ln\left(\frac{Q^2}{m^2}\right) \left( \frac{1}{6} \right) = \frac{\alpha}{3\pi} \ln\left(\frac{Q^2}{m^2}\right) 比较两者比例: Π^scalar(Q2)Π^spinor(Q2)=α/12πα/3π=14\boxed{ \frac{\hat{\Pi}_{\text{scalar}}(Q^2)}{\hat{\Pi}_{\text{spinor}}(Q^2)} = \frac{\alpha / 12\pi}{\alpha / 3\pi} = \frac{1}{4} } 这证明了带电标量玻色子的贡献恰好是虚电子-正电子对贡献的 1/4。
9.2

Problem 9.2

peskinChapter 9

习题 9.2

来源: 第9章, PDF第312,313,314页


9.2 Quantum statistical mechanics.

(a) Evaluate the quantum statistical partition function

Z=tr[eβH]Z = \text{tr}[e^{-\beta H}]

(where β=1/kT\beta = 1/kT) using the strategy of Section 9.1 for evaluating the matrix elements of eiHte^{-iHt} in terms of functional integrals. Show that one again finds a functional integral, over functions defined on a domain that is of length β\beta and periodically connected in the time direction. Note that the Euclidean form of the Lagrangian appears in the weight.

(b) Evaluate this integral for a simple harmonic oscillator,

LE=12x˙2+12ω2x2,L_E = \frac{1}{2}\dot{x}^2 + \frac{1}{2}\omega^2 x^2,

by introducing a Fourier decomposition of x(t)x(t):

x(t)=nxn1βe2πint/β.x(t) = \sum_{n} x_n \cdot \frac{1}{\sqrt{\beta}} e^{2\pi int/\beta}.

The dependence of the result on β\beta is a bit subtle to obtain explicitly, since the measure for the integral over x(t)x(t) depends on β\beta in any discretization. However, the dependence on ω\omega should be unambiguous. Show that, up to a (possibly divergent and β\beta-dependent) constant, the integral reproduces exactly the familiar expression for the quantum partition function of an oscillator. [You may find the identity

sinhz=zn=1(1+z2(nπ)2)\sinh z = z \cdot \prod_{n=1}^{\infty} \left( 1 + \frac{z^2}{(n\pi)^2} \right)

useful.]

(c) Generalize this construction to field theory. Show that the quantum statistical partition function for a free scalar field can be written in terms of a functional integral. The value of this integral is given formally by

[det(2+m2)]1/2,\left[ \det(-\partial^2 + m^2) \right]^{-1/2},

where the operator acts on functions on Euclidean space that are periodic in the time direction with periodicity β\beta. As before, the β\beta dependence of this expression is difficult to compute directly. However, the dependence on m2m^2 is unambiguous. (More generally, one can usually evaluate the variation of a functional determinant with respect to any explicit parameter in the Lagrangian.) Show that the determinant indeed reproduces the partition function for relativistic scalar particles.

(d) Now let ψ(t),ψˉ(t)\psi(t), \bar{\psi}(t) be two Grassmann-valued coordinates, and define a fermionic oscillator by writing the Lagrangian

LE=ψˉψ˙+ωψˉψ.L_E = \bar{\psi}\dot{\psi} + \omega\bar{\psi}\psi.

This Lagrangian corresponds to the Hamiltonian

H=ωψˉψ,with {ψˉ,ψ}=1;H = \omega\bar{\psi}\psi, \quad \text{with } \{\bar{\psi}, \psi\} = 1;

that is, to a simple two-level system. Evaluate the functional integral, assuming that the fermions obey antiperiodic boundary conditions: ψ(t+β)=ψ(t)\psi(t + \beta) = -\psi(t). (Why is this reasonable?) Show that the result reproduces the partition function of a quantum-mechanical two-level system, that is, of a quantum state with Fermi statistics.

(e) Define the partition function for the photon field as the gauge-invariant functional integral

Z=DAexp(d4xE[14(Fμν)2])Z = \int \mathcal{D}A \exp \left( - \int d^4x_E \left[ \frac{1}{4} (F_{\mu\nu})^2 \right] \right)

over vector fields AμA_\mu that are periodic in the time direction with period β\beta. Apply the gauge-fixing procedure discussed in Section 9.4 (working, for example, in Feynman gauge). Evaluate the functional determinants using the result of part (c) and show that the functional integral does give the correct quantum statistical result (including the correct counting of polarization states).

习题 9.2 - 解答


(a) 在量子力学中,演化算符 eiHte^{-iHt} 的矩阵元可以通过路径积分表示为: xfeiHtxi=Dxexp(i0tdtL(x,x˙))\langle x_f | e^{-iHt} | x_i \rangle = \int \mathcal{D}x \exp\left( i \int_0^t dt' L(x, \dot{x}) \right) 为了计算配分函数 Z=tr[eβH]Z = \text{tr}[e^{-\beta H}],我们进行 Wick 转动,将时间解析延拓到虚时间 τ=it\tau = it。令 t=iβt = -i\beta,演化算符变为 eβHe^{-\beta H}。此时作用量变为: i0iβdt[12m(dxdt)2V(x)]=0βdτ[12m(dxdτ)2V(x)]=0βdτLEi \int_0^{-i\beta} dt' \left[ \frac{1}{2}m\left(\frac{dx}{dt'}\right)^2 - V(x) \right] = \int_0^\beta d\tau \left[ -\frac{1}{2}m\left(\frac{dx}{d\tau}\right)^2 - V(x) \right] = -\int_0^\beta d\tau L_E 其中 LE=12mx˙2+V(x)L_E = \frac{1}{2}m\dot{x}^2 + V(x) 是欧几里得拉格朗日量。因此,虚时间演化的矩阵元为: xfeβHxi=x(0)=xix(β)=xfDxexp(0βdτLE)\langle x_f | e^{-\beta H} | x_i \rangle = \int_{x(0)=x_i}^{x(\beta)=x_f} \mathcal{D}x \exp\left( - \int_0^\beta d\tau L_E \right) 配分函数要求对该矩阵元求迹,即令 xf=xi=xx_f = x_i = x 并对 xx 积分。这等价于对所有满足周期性边界条件 x(0)=x(β)x(0) = x(\beta) 的路径求泛函积分: Z=Dxexp(0βdτLE)\boxed{ Z = \oint \mathcal{D}x \exp\left( - \int_0^\beta d\tau L_E \right) } 积分域为在长度为 β\beta 的时间区间上定义的周期函数。

(b) 对于简谐振子,欧几里得拉格朗日量为 LE=12x˙2+12ω2x2L_E = \frac{1}{2}\dot{x}^2 + \frac{1}{2}\omega^2 x^2。引入傅里叶展开: x(τ)=n=xn1βe2πinτ/βx(\tau) = \sum_{n=-\infty}^{\infty} x_n \frac{1}{\sqrt{\beta}} e^{2\pi i n \tau / \beta} 由于 x(τ)x(\tau) 是实函数,有 xn=xnx_{-n} = x_n^*。代入欧几里得作用量 SE=0βdτLES_E = \int_0^\beta d\tau L_E 中,利用正交性可得: SE=12n=xn2[(2πnβ)2+ω2]S_E = \frac{1}{2} \sum_{n=-\infty}^{\infty} |x_n|^2 \left[ \left(\frac{2\pi n}{\beta}\right)^2 + \omega^2 \right] 泛函积分测度变为 Dx=N(β)ndxn\mathcal{D}x = \mathcal{N}(\beta) \prod_n dx_n。执行高斯积分: Z=N(β)n=[(2πnβ)2+ω2]1/2Z = \mathcal{N}'(\beta) \prod_{n=-\infty}^{\infty} \left[ \left(\frac{2\pi n}{\beta}\right)^2 + \omega^2 \right]^{-1/2}n=0n=0nnn-n 的项配对合并,提取出与 ω\omega 相关的部分: Z=C(β)1ωn=1[(2πnβ)2+ω2]1=C(β)1ωn=1[1+ω2β2(2πn)2]1Z = C(\beta) \frac{1}{\omega} \prod_{n=1}^{\infty} \left[ \left(\frac{2\pi n}{\beta}\right)^2 + \omega^2 \right]^{-1} = C'(\beta) \frac{1}{\omega} \prod_{n=1}^{\infty} \left[ 1 + \frac{\omega^2 \beta^2}{(2\pi n)^2} \right]^{-1} 利用提示中的恒等式 sinhz=zn=1(1+z2n2π2)\sinh z = z \prod_{n=1}^{\infty} \left( 1 + \frac{z^2}{n^2\pi^2} \right),令 z=ωβ2z = \frac{\omega\beta}{2},可得: n=1[1+ω2β2(2πn)2]1=ωβ/2sinh(ωβ/2)\prod_{n=1}^{\infty} \left[ 1 + \frac{\omega^2 \beta^2}{(2\pi n)^2} \right]^{-1} = \frac{\omega\beta/2}{\sinh(\omega\beta/2)} 因此,配分函数为: Z1ωωβ/2sinh(ωβ/2)1sinh(ωβ/2)=2eβω/21eβωZ \propto \frac{1}{\omega} \frac{\omega\beta/2}{\sinh(\omega\beta/2)} \propto \frac{1}{\sinh(\omega\beta/2)} = \frac{2 e^{-\beta\omega/2}}{1 - e^{-\beta\omega}} 忽略与 β\beta 相关的常数因子,这精确给出了量子简谐振子的配分函数: Zeβω/21eβω=n=0eβω(n+1/2)\boxed{ Z \propto \frac{e^{-\beta\omega/2}}{1 - e^{-\beta\omega}} = \sum_{n=0}^{\infty} e^{-\beta\omega(n+1/2)} }

(c) 对于自由标量场,欧几里得拉格朗日量为 LE=12(τϕ)2+12(ϕ)2+12m2ϕ2L_E = \frac{1}{2}(\partial_\tau \phi)^2 + \frac{1}{2}(\nabla \phi)^2 + \frac{1}{2}m^2 \phi^2。泛函积分形式上给出: Z=Dϕexp(0βdτd3xLE)=[det(2+m2)]1/2Z = \int \mathcal{D}\phi \exp\left( - \int_0^\beta d\tau \int d^3x L_E \right) = \left[ \det(-\partial^2 + m^2) \right]^{-1/2} 取对数并利用 lndet=trln\ln \det = \text{tr} \ln,算符 2+m2-\partial^2 + m^2 的本征值为 ωn2+Ep2\omega_n^2 + E_{\mathbf{p}}^2,其中 ωn=2πnβ\omega_n = \frac{2\pi n}{\beta} 为 Matsubara 频率,Ep=p2+m2E_{\mathbf{p}} = \sqrt{\mathbf{p}^2 + m^2}lnZ=12pn=ln(ωn2+Ep2)\ln Z = -\frac{1}{2} \sum_{\mathbf{p}} \sum_{n=-\infty}^{\infty} \ln\left( \omega_n^2 + E_{\mathbf{p}}^2 \right) 为了明确 m2m^2 的依赖关系,对 m2m^2 求导: lnZm2=12pn=1(2πn/β)2+Ep2\frac{\partial \ln Z}{\partial m^2} = -\frac{1}{2} \sum_{\mathbf{p}} \sum_{n=-\infty}^{\infty} \frac{1}{(2\pi n/\beta)^2 + E_{\mathbf{p}}^2} 利用求和公式 n=1n2+a2=πacoth(πa)\sum_{n=-\infty}^{\infty} \frac{1}{n^2 + a^2} = \frac{\pi}{a} \coth(\pi a),令 a=βEp2πa = \frac{\beta E_{\mathbf{p}}}{2\pi},得到: n=1(2πn/β)2+Ep2=β24π2πβEp/2πcoth(βEp2)=β2Epcoth(βEp2)\sum_{n=-\infty}^{\infty} \frac{1}{(2\pi n/\beta)^2 + E_{\mathbf{p}}^2} = \frac{\beta^2}{4\pi^2} \frac{\pi}{\beta E_{\mathbf{p}} / 2\pi} \coth\left(\frac{\beta E_{\mathbf{p}}}{2}\right) = \frac{\beta}{2E_{\mathbf{p}}} \coth\left(\frac{\beta E_{\mathbf{p}}}{2}\right) 因此: lnZm2=pβ4Epcoth(βEp2)\frac{\partial \ln Z}{\partial m^2} = -\sum_{\mathbf{p}} \frac{\beta}{4E_{\mathbf{p}}} \coth\left(\frac{\beta E_{\mathbf{p}}}{2}\right) 另一方面,相对论性自由标量粒子的标准配分函数为 Zstd=peβEp/21eβEpZ_{std} = \prod_{\mathbf{p}} \frac{e^{-\beta E_{\mathbf{p}}/2}}{1 - e^{-\beta E_{\mathbf{p}}}}。对其取对数并对 m2m^2 求导(注意 Epm2=12Ep\frac{\partial E_{\mathbf{p}}}{\partial m^2} = \frac{1}{2E_{\mathbf{p}}}): lnZstdm2=pEp[βEp2ln(1eβEp)]12Ep=pβ4Ep1+eβEp1eβEp=pβ4Epcoth(βEp2)\frac{\partial \ln Z_{std}}{\partial m^2} = \sum_{\mathbf{p}} \frac{\partial}{\partial E_{\mathbf{p}}} \left[ -\frac{\beta E_{\mathbf{p}}}{2} - \ln(1 - e^{-\beta E_{\mathbf{p}}}) \right] \frac{1}{2E_{\mathbf{p}}} = -\sum_{\mathbf{p}} \frac{\beta}{4E_{\mathbf{p}}} \frac{1 + e^{-\beta E_{\mathbf{p}}}}{1 - e^{-\beta E_{\mathbf{p}}}} = -\sum_{\mathbf{p}} \frac{\beta}{4E_{\mathbf{p}}} \coth\left(\frac{\beta E_{\mathbf{p}}}{2}\right) 两者完全一致,证明了: Z=[det(2+m2)]1/2 准确给出了相对论性标量粒子的配分函数\boxed{ Z = \left[ \det(-\partial^2 + m^2) \right]^{-1/2} \text{ 准确给出了相对论性标量粒子的配分函数} }

(d) 对于费米子振子 LE=ψˉψ˙+ωψˉψL_E = \bar{\psi}\dot{\psi} + \omega\bar{\psi}\psi反周期边界条件的合理性:在费米子相空间中,算符的迹可以通过 Grassmann 变量的相干态表示为 trA=dψˉdψeψˉψψAψ\text{tr} A = \int d\bar{\psi} d\psi e^{-\bar{\psi}\psi} \langle -\psi | A | \psi \rangle。由于交换 Grassmann 变量会产生负号,左矢带有负号,这在路径积分中自然导致了反周期边界条件 ψ(τ+β)=ψ(τ)\psi(\tau+\beta) = -\psi(\tau)

引入反周期傅里叶展开,Matsubara 频率为 ωn=(2n+1)πβ\omega_n = \frac{(2n+1)\pi}{\beta}。作用量变为: SE=n=ψˉn(iωn+ω)ψnS_E = \sum_{n=-\infty}^{\infty} \bar{\psi}_n (-i\omega_n + \omega) \psi_n 对 Grassmann 变量的泛函积分给出费米型行列式: Z=n=(iωn+ω)Z = \prod_{n=-\infty}^{\infty} (-i\omega_n + \omega)ω\omega 求导以计算其依赖关系: lnZω=n=1iωn+ω=n=ωωn2+ω2=n=ω(2n+1)2π2β2+ω2\frac{\partial \ln Z}{\partial \omega} = \sum_{n=-\infty}^{\infty} \frac{1}{-i\omega_n + \omega} = \sum_{n=-\infty}^{\infty} \frac{\omega}{\omega_n^2 + \omega^2} = \sum_{n=-\infty}^{\infty} \frac{\omega}{\frac{(2n+1)^2\pi^2}{\beta^2} + \omega^2} 利用求和公式 n=1(2n+1)2+a2=π2atanh(πa2)\sum_{n=-\infty}^{\infty} \frac{1}{(2n+1)^2 + a^2} = \frac{\pi}{2a} \tanh\left(\frac{\pi a}{2}\right),令 a=βωπa = \frac{\beta\omega}{\pi},得到: lnZω=β2tanh(βω2)\frac{\partial \ln Z}{\partial \omega} = \frac{\beta}{2} \tanh\left(\frac{\beta\omega}{2}\right) 积分得到 lnZ=lncosh(βω2)+C\ln Z = \ln \cosh(\frac{\beta\omega}{2}) + C,即: Zcosh(βω2)=eβω/2+eβω/22Z \propto \cosh\left(\frac{\beta\omega}{2}\right) = \frac{e^{\beta\omega/2} + e^{-\beta\omega/2}}{2} 这正是哈密顿量为 Weyl 排序 H=ω2(ψˉψψψˉ)=ω(ψˉψ12)H = \frac{\omega}{2}(\bar{\psi}\psi - \psi\bar{\psi}) = \omega(\bar{\psi}\psi - \frac{1}{2}) 的二能级系统(能级为 ±ω/2\pm \omega/2)的配分函数: Zeβω/2+eβω/2\boxed{ Z \propto e^{\beta\omega/2} + e^{-\beta\omega/2} }

(e) 光子场的配分函数为: Z=DAexp(d4xE14(Fμν)2)Z = \int \mathcal{D}A \exp \left( - \int d^4x_E \frac{1}{4} (F_{\mu\nu})^2 \right) 在 Feynman 规范(ξ=1\xi=1)下,加入规范固定项 12(μAμ)2\frac{1}{2}(\partial_\mu A_\mu)^2。通过分部积分,规范玻色子的拉格朗日量变为 12Aμ(2)Aμ\frac{1}{2} A_\mu (-\partial^2) A_\mu。对 4 个分量的 AμA_\mu 积分给出: ZA=[det(2)]4/2=[det(2)]2Z_A = \left[ \det(-\partial^2) \right]^{-4/2} = \left[ \det(-\partial^2) \right]^{-2} 同时必须引入 Faddeev-Popov 鬼场 c,cˉc, \bar{c},其拉格朗日量为 cˉ(2)c\bar{c}(-\partial^2)c。由于鬼场是为了消除 AμA_\mu 中非物理的纵向和时间分量,它们必须与 AμA_\mu 具有相同的边界条件(即周期性边界条件)。鬼场是 Grassmann 变量,其路径积分给出正的行列式: Zghost=det(2)Z_{ghost} = \det(-\partial^2) 总的配分函数为: Z=ZA×Zghost=[det(2)]2det(2)=[det(2)]1Z = Z_A \times Z_{ghost} = \left[ \det(-\partial^2) \right]^{-2} \det(-\partial^2) = \left[ \det(-\partial^2) \right]^{-1} 根据 (c) 问的结果,单个无质量标量场(m=0m=0)的配分函数为 [det(2)]1/2\left[ \det(-\partial^2) \right]^{-1/2}。因此,光子场的配分函数等价于两个独立的无质量标量场的配分函数之积: Z=(peβp/21eβp)2=peβp(1eβp)2\boxed{ Z = \left( \prod_{\mathbf{p}} \frac{e^{-\beta|\mathbf{p}|/2}}{1 - e^{-\beta|\mathbf{p}|}} \right)^2 = \prod_{\mathbf{p}} \frac{e^{-\beta|\mathbf{p}|}}{(1 - e^{-\beta|\mathbf{p}|})^2} } 这精确给出了光子气体的量子统计结果,且正确地给出了光子的 2 个横向物理偏振态的计数。