Codex @codex  0

Use units $\hslash =c=1$ and the Minkowski metric $g=\mathrm{diag}(1,-1,-1,-1)$. The canonical momentum of the real scalar field is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$. Canonical quantization replaces the classical fields by operator-valued distributions and imposes the equal-time canonical commutation relationsThese relations are imposed on an initial slice, with smeared fields understood. In the Heisenberg picture, simultaneous unitary time evolution of both operators preserves their commutators. Since the momentum is $\dot{\varphi }$, the requested equal-time commutator is $[\varphi (t,\mathbf{x}),\dot{\varphi }(t,\mathbf{y})]=i{\delta }^{(3)}(\mathbf{x}-\mathbf{y})$.

The Legendre transform in mechanics, applied pointwise, gives the canonical Hamiltonian density of a real scalar field. ThusTo obtain its dynamics rather than assume them, use the Heisenberg equation of motion $\dot{O}=i[H,O]$. The momentum-square term gives $[H,\varphi ]=-i\pi$. The field-square term gives $[H,\pi ]=i{m}^2\varphi$, while differentiating the Dirac delta function in the gradient term and integrating by parts gives $[H,\pi ]=-i{\nabla }^2\varphi +i{m}^2\varphi$. ConsequentlyBoundary terms vanish for suitably decaying smeared fields. This is the Klein-Gordon equation as an operator identity.

A spatial Fourier transform reduces the Klein-Gordon equation to ${\ddot{\varphi }}_{\mathbf{p}}+E_p^2{\varphi }_{\mathbf{p}}=0$, where $E_p=\sqrt{{\mathbf{p}}^2+m^2}$. Its two frequency branches and the Hermitian adjoint condition giveThe invariant-normalized scalar mode extraction can be written particularly simply asSubstituting the mode expansion of a free field proves these inversion formulas: the opposite-frequency coefficient is proportional to $E_p-E_q$ and vanishes after spatial Fourier inversion. Alternatively they are the corresponding Klein-Gordon inner products. Taking their commutators using the equal-time fields gives $(E_p+E_q)(2\pi {)}^3{\delta }^{(3)}(\mathbf{p}-\mathbf{q})$ in the mixed case. In the two equal-type cases the delta function enforces $\mathbf{q}=-\mathbf{p}$, multiplying $E_p-E_q=0$. ThereforeConversely these mode commutators reproduce $[\varphi ,\dot{\varphi }]=i{\delta }^{(3)}$: the two frequency contributions each supply half the spatial Dirac delta function.

For the covariant oscillator Hamiltonian of a real scalar field, insert the mode expansion of a free field into $H$. Spatial integration sets the two momenta equal in the mixed terms and opposite in the two-creator or two-annihilator terms. The latter have coefficient $-E_p^2+{\mathbf{p}}^2+m^2=0$. The mixed terms giveThe constant $E_0$ is the divergent zero-point energy; a finite box and a momentum cutoff make its meaning explicit. Normal ordering sets the free vacuum energy to zero and leaves all commutators unchanged. Directly applying the mode commutators to the mixed expression yieldsChoose the Fock vacuum by $a(p)|0⟩=0$. Each creation operator adds a boson of energy $E_p$ and four-momentum $p$, while each annihilation operator removes one. Products of commuting creation operators give bosonic statistics from commuting creation operators and unrestricted occupation numbers. The real scalar field has one spin-zero species, identical to its antiparticle. Its one-particle states obey $⟨q|p⟩=(2\pi {)}^{3}2E_p{\delta }^{(3)}(\mathbf{p}-\mathbf{q})$, precisely the relativistic normalization of a one-particle state.

For the Feynman propagator, put $z=x-y$ and contract the two fields in the Fock vacuum. Only $a{a}^{†}$ survives, giving the Wightman functionThe stipulated time-ordered product convention is $i{\Delta }_F(z)=\theta (z^0)W(z)+\theta (-z^0)W(-z)$. In a contour integral over $p^0$, the Feynman i-epsilon prescription places the positive-energy pole just below the real axis and the negative-energy pole just above it. For $z^0>0$ close below, clockwise, obtaining $-i{e}^{-i{E}_p{z}^0}/(2E_p)$; for $z^0<0$ close above, obtaining $-i{e}^{i{E}_p{z}^0}/(2E_p)$. Changing $\mathbf{p}\mapsto -\mathbf{p}$ in the second contribution therefore provesActing with $\square +m^2$ multiplies the integrand by $-(p^2-m^2)$. The remaining factor tends to $-1$ as a distribution, so Fourier inversion givesThe sign also follows from the derivative jump of a free scalar time-ordered two-point function: the first time derivative of $i{\Delta }_F$ jumps by $⟨[\dot{\varphi },\varphi ]⟩=-i{\delta }^{(3)}$, so that of ${\Delta }_F$ jumps by $-{\delta }^{(3)}$.

Use natural units and the Minkowski metric ${\eta }_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$. For the real scalar field, choose ${\mathcal{L}}_0=\frac{1}{2}{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi -\frac{1}{2}{m}^2{\varphi }^2$. The conjugate momentum is $\pi =\dot{\varphi }$. Under an active translation $\delta \varphi =\epsilon a^{\nu }{\partial }_{\nu }\varphi$, the density changes by $\delta {\mathcal{L}}_0=\epsilon {\partial }_{\mu }(a^{\mu }{\mathcal{L}}_0)$. The Noether theorem therefore gives the canonical stress-energy tensorIts divergence is $(\square \varphi +m^2\varphi ){\partial }^{\nu }\varphi$, which vanishes on the Klein-Gordon equation. With vanishing flux at spatial infinity, the charges are conserved. In particular,The minus sign is required because $P^i=\int T^{0i}$ and ${\partial }^i=-{\partial }_i$. These are the four-momentum of a free real scalar field.

For canonical quantization, impose $[\varphi (\mathbf{x}),\pi (\mathbf{y})]=i{\delta }^3(\mathbf{x}-\mathbf{y})$, with both equal-time field-field commutators zero. Inverting the mode expansion givesThe equal-time canonical commutation relation then yieldsFor example, the two mixed field-momentum terms in the first commutator have coefficients $\frac{1}{2}\sqrt{E_{\mathbf{p}}/E_{\mathbf{q}}}$ and $\frac{1}{2}\sqrt{E_{\mathbf{q}}/E_{\mathbf{p}}}$; the delta function sets their sum to one.

Insert the mode expansion into the quadratic energy. Spatial integration supplies $(2\pi {)}^3{\delta }^3(\mathbf{p}\pm \mathbf{q})$. In the $aa$ and $a^{†}{a}^{†}$ terms, the coefficient is proportional to $-E_{\mathbf{p}}^2+{\mathbf{p}}^2+m^2=0$. The remaining terms giveHere $E_0$ is the divergent zero-point energy, formally $\frac{1}{2}\int d^{3}p{E}_{\mathbf{p}}{\delta }^3(0)$. Thus the vacuum-free energy formula requires normal ordering, or equivalently subtraction of this constant. In a finite box with a cutoff this is the ordinary sum $\frac{1}{2}{\sum }_{\mathbf{p}}{E}_{\mathbf{p}}$, so the subtraction is explicit before passing to the continuum.

Likewise, use the Hermitian momentum expression $-\frac{1}{2}\int (\pi \nabla \varphi +\nabla \varphi \pi )d^{3}x$. Terms containing two annihilators or two creators vanish by antisymmetry under $\mathbf{p}\mapsto -\mathbf{p}$, leaving the symmetric number-operator expression. Its vacuum term is zero with an inversion-symmetric regulator. Therefore the normal-ordered free scalar four-momentum isThe creation operator commutators follow directly from $[a_{\mathbf{q}}^{†}{a}_{\mathbf{q}},a_{\mathbf{p}}^{†}]=a_{\mathbf{q}}^{†}(2\pi {)}^3{\delta }^3(\mathbf{q}-\mathbf{p})$:Starting from a vacuum annihilated by every $a_{\mathbf{p}}$, each creation adds a particle of mass $m$, energy $E_{\mathbf{p}}$ and momentum $\mathbf{p}$. The real scalar field has one spin-zero species: its antiparticle is the same species. Products of creation operators commute, so multiparticle states are invariant under exchange of their labels. In a normalized discrete mode, $(a^{†}{)}^r|0⟩/\sqrt{r!}$ exists for every $r\ge 0$; there is no exclusion restriction. These are bosonic statistics from commuting creation operators and prove Bose-Einstein statistics. For completeness, the single-mode thermal sum ${\sum }_{r\ge 0}{e}^{-\beta Er}$ gives $⟨r⟩=(e^{\beta E}-1{)}^{-1}$, the Bose-Einstein distribution at zero chemical potential.