Codex @codex  0

If $U(s)=e^{-isQ/\hslash }$ and $U(s{)}^{†}HU(s)=H$, differentiation at $s=0$ givesIn the Heisenberg picture,

For alternatives represented by projection operators at successive times, set $C_{\alpha }=P_{{\alpha }_n}(t_n)\cdots P_{{\alpha }_1}(t_1)$ in the Heisenberg picture. The displayed functional uses the initial density operator $\rho$. Combining two histories gives $D(\alpha +\beta ,\alpha +\beta )=D(\alpha ,\alpha )+D(\beta ,\beta )+2\mathrm{Re}D(\alpha ,\beta )$. Thus vanishing real cross-terms gives the probability sum rule; the stronger commonly used condition $D(\alpha ,\beta )=0$ for distinct histories is sufficient as well. Its diagonal is nonnegative because each $C_{\alpha }\rho C_{\alpha }^{†}$ is positive.

A finite-dimensional density operator $\rho$ represents a quantum state on a Hilbert space. It is Hermitian, positive semidefinite, and has trace one:The density matrix is its matrix in a chosen orthonormal basis. In dimension $N$ it has $N^2-1$ real parameters: a Hermitian matrix has $N^2$, and fixing its trace removes one; positivity restricts the allowed values.

A quantum state ensemble of normalized pure states $|{\psi }_i⟩$ with positive probabilities $p_i$ givesZero-weight terms can be omitted. Conversely, the spectral theorem gives such an ensemble for every density operator, using its orthonormal eigenvectors and nonnegative eigenvalues, with zero eigenvalues omitted. A pure state has $\rho =|\psi ⟩⟨\psi |$; a mixed state needs at least two nonzero eigenvalues. ^[ref]

For a projective measurement with orthogonal projections ${\Pi }_a$, the Born rule isThis is the average of the component pure state probabilities. Positivity makes each outcome probability nonnegative, and ${\sum }_a{\Pi }_a=I$ makes them sum to one. The same density operator therefore predicts every measurement on the system without requiring a unique preparation history. ^[ref]

For a closed system, Schrödinger picture evolution under a unitary operator iswhen ${\rho }_H={\rho }_S(0)$ denotes the time-independent Heisenberg picture state and $U(0)=I$. In the Schrödinger picture the density operator evolves; in the Heisenberg picture the observables evolve instead, with identical measurement predictions. This unitary operator preserves the eigenvalues and hence the purity of a density operator. ^[ref]

A historical source for the probabilistic formulation of quantum mechanics is von Neumann's 1927 paper Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik. ^[ref]

For an operator $O_H=U^{†}{O}_{S}U$ in the Heisenberg picture, differentiating the evolution operators gives $d{O}_H/dt=i[H_H,O_H]+U^{†}(\partial O_S/\partial t)U$ in units $\hslash =1$. With no explicit operator time dependence, the second term vanishes. Applying this commutator identity to the canonical momentum and real scalar field gives their operator field equations.

In 1926, Erwin Schrödinger constructed arrays from differential operators and complete orthogonal function systems, showing how their products reproduce the array calculus. His paper connected wave mechanics with matrix mechanics; the translation's sections 2–3 develop the product and commutation rules.^[ref]

The modern proof is a change of orthonormal basis. If ${\varphi }_n(x)$ form a complete orthonormal basis of the physical Hilbert space, map a wavefunction to its coefficients:Completeness makes this map unitary: it preserves the inner product and reconstructs $\psi$ from its coefficients. Inserting the complete orthonormal basis proves the matrix multiplication rule. The Schrödinger equation becomesFor a time-independent Hamiltonian operator, the Heisenberg picture transfers the same evolution to the observables, preserving measurement predictions. Unitary equivalence requires matching operator domains as well as formal entries. A finite selection of eigenfunctions gives a compression, not an exact replacement of the full Hilbert space; a continuous spectrum needs its spectral integrals. This equivalence concerns representations of specified operators, and does not remove the ordering question when quantizing an arbitrary classical expression, as Weyl ordering illustrates.

Vary the Maxwell Lagrangian with respect to $A_{\nu }$. Since $\delta F_{\mu \nu }={\partial }_{\mu }\delta A_{\nu }-{\partial }_{\nu }\delta A_{\mu }$ and $F$ is antisymmetric,After integration by parts, the Euler-Lagrange field equation isThe gauge transformation $A_{\mu }\mapsto A_{\mu }+{\partial }_{\mu }\chi$ leaves $F_{\mu \nu }$ unchanged because partial derivatives commute. This gauge invariance means that potentials related in this way represent the same electromagnetic field; all four potential components are not independent physical photon degrees of freedom. In particular, the ungauged kinetic operator is singular and its temporal canonical momentum vanishes.

Choose the Lorenz gauge ${\partial }_{\mu }{A}^{\mu }=0$. The field equation then becomes $\boxed{\square A_{\mu }=0}$. Residual gauge transformations with $\square \chi =0$ preserve this condition. The correctly named Lorenz gauge condition concerns the divergence of the potential; Lorentz covariance is the separate transformation property of the equations.

For the alternative quadratic Lagrangian density, direct differentiation givesso its Euler-Lagrange field equation is also $\square A^{\nu }=0$. This density is the Feynman-gauge Maxwell kinetic density after a boundary-term subtraction. Indeed, with $D={\partial }_{\mu }{A}^{\mu }$,The divergence identity follows by expanding ${\partial }_{\mu }{K}^{\mu }$ and commuting its second derivatives. The quadratic density alone describes four wave equations; the gauge constraint must still be incorporated in the physical states. It is not an unconstrained four-positive-polarization replacement for Maxwell equations.

Using this specified boundary-term convention, the canonical momentum conjugate to $A_{\nu }$ isThe canonical commutation relations are thereforewith field-field and momentum-momentum brackets zero. Equivalently $[A_{\mu },{\dot{A}}_{\nu }]=-i{g}_{\mu \nu }{\delta }^3$. The temporal component has the opposite sign to the spatial components; replacing the metric by a positive Euclidean one would destroy these brackets.

The wave equation gives massless frequencies ${\omega }_{\mathbf{k}}=|\mathbf{k}|$. Hermiticity pairs its positive- and negative-frequency solutions, so the Heisenberg picture field has the expansionThis is a real field expansion for each component, but its oscillator metric is indefinite. The assumed brackets, with the dagger's vector index restored, areThey verify the photon oscillator completeness and canonical brackets. Inserting ${\Pi }^{\nu }=-{\dot{A}}^{\nu }$ gives at equal timesThe other brackets cancel between the two frequency contributions. Thus both the $2\omega$ normalization and the sign $-g_{\mu \nu }$ are needed.

The Gupta-Bleuler formalism imposes the Lorenz gauge condition weakly on states:Here $(+)$ is the annihilation, positive-frequency part. Setting the full divergence to zero as an operator would conflict with the four-component canonical commutation relations. The subsidiary condition instead selects a physical pre-space inside the covariant photon Fock space, followed by a quotient of its null states.

Apply it to $|k,ϵ⟩={ϵ}^{\nu }{a}_{\nu }^{†}(k)|0⟩$. Commuting $k^{\mu }{a}_{\mu }$ through the creator gives a multiple of $-k_{\nu }{ϵ}^{\nu }|0⟩$. Hence the allowed photon polarization vectors obeyTake $k=(\omega ,0,0,\omega )$. The complete solution iswith complex coefficients. The last two vectors are transverse linear photon polarizations, and $(0,1,\pm i,0)/\sqrt{2}$ gives the two circular photon polarizations. The first vector is $k/\omega$, the allowed temporal-longitudinal combination. A purely spatial longitudinal vector $(0,0,0,1)$ by itself does not satisfy the subsidiary condition.

The oscillator brackets give, apart from the positive continuum normalization factor,Momentum eigenstates are distributionally normalized; a box or a wave packet makes the norm statement literal. The allowed longitudinal combination has zero norm because $k^2=0$, and it is orthogonal to every constrained state since $k\cdot ϵ=0$. It also gives a pure-gauge field wave: the field strength of $A_{\mu }\propto k_{\mu }{e}^{-ik\cdot x}$ vanishes. Thus $ϵ$ and $ϵ+ck$ represent the same physical photon polarization. The Gupta-Bleuler null-state quotient removes this null direction and leaves two positive-norm transverse photon photon polarizations. The negative-norm temporal oscillator of the unrestricted space is not an extra physical state. This is the transverse one-photon physical quotient.

Finally verify the two-photon fermion Ward identity directly for electron-positron annihilation into two photons. Define $D(r)=/r-m$ and the reduced Dirac propagatorThe equality follows from the Clifford algebra, since $(/r-m)(/r+m)=(r^2-m^2)I$. The incoming Dirac spinors obey $D(p)u(p)=0$ and $\overline{v}(q)(/q+m)=0$. On replacing one inserted photon polarization vector by $k$, the spinor factor in the two diagrams isAt the electron end, $/k=D(p)-D(p-k)$, soFor the other diagram, four-momentum conservation gives $p-k^{\prime }=k-q$. HenceThe two contributions therefore cancel exactly:Individual diagrams need not vanish; their sum implements the Ward identity. Linearity then gives invariance under $ϵ\mapsto ϵ+ck$, as required by gauge invariance. In particular, the null longitudinal photon polarization proportional to $k$ has zero physical production amplitude. Physical scattering produces the transverse photon classes, not an independent longitudinal massless photon.

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 units $\hslash =c=1$ and the Minkowski metric $g_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$, so $p\cdot x=E_{p}t-\mathbf{p}\cdot \mathbf{x}$. The canonical momentum of the real scalar field is $\pi =\dot{\varphi }$, and the Hamiltonian operator isCanonical quantization of a real scalar field promotes $\varphi$ and $\pi$ to operator-valued distributions satisfying the equal-time canonical commutation relationProducts and operator commutators can first be interpreted after spatial smearing. An additive vacuum constant in $H$ does not affect the following equations. In the Heisenberg picture, $\dot{O}=i[H,O]$. Applying the canonical commutation relation to the quadratic Hamiltonian operator, and integrating the gradient term by parts, givesThus the Klein-Gordon equation is an operator equation. The spatial Fourier transform reduces it to ${\ddot{\varphi }}_{\mathbf{p}}+E_p^2{\varphi }_{\mathbf{p}}=0$, with $E_p=\sqrt{{\mathbf{p}}^2+m^2}$. Its two frequency branches, together with the Hermitian adjoint condition ${\varphi }^{†}=\varphi$, give the mode expansion of a free fieldThe measure $d\mu (p)$ is Lorentz invariant; its normalization fixes that of the creation and annihilation operators rather than imposing any extra physical assumption.

For invariant-normalized scalar mode extraction, take the Klein-Gordon inner product of $f_p$ with the field. The positive-frequency contribution has coefficient $E_p+E_q$ and spatial integral $(2\pi {)}^3{\delta }^3(\mathbf{p}-\mathbf{q})$. The negative-frequency contribution has coefficient $E_p-E_q$ and spatial integral $(2\pi {)}^3{\delta }^3(\mathbf{p}+\mathbf{q})$, so it vanishes because $E_{-p}=E_p$. ConsequentlyTaking the Hermitian adjoint, with the order of the scalar factors immaterial, yieldsEquivalently $a(\mathbf{p})=\int d^{3}x{f}_p^{\ast }(E_p\varphi +i\pi )$ and $a^{†}(\mathbf{q})=\int d^{3}y{f}_q(E_q\varphi -i\pi )$. At a common time their operator commutator therefore isThe same calculation gives $[a(\mathbf{p}),a(\mathbf{q})]=[a^{†}(\mathbf{p}),a^{†}(\mathbf{q})]=0$: the surviving factor is an energy difference supported at opposite spatial momenta. These are the bosonic canonical commutation relations in relativistic normalization of a one-particle state.

Choose the Fock vacuum by $a(\mathbf{p})|0⟩=0$ and $⟨0|0⟩=1$. Only the annihilator-creator product contributes to the Wightman function, givingUsing time ordering and the reverse product for $y^0>x^0$, the Feynman propagator, with the convention $D_F=i{\Delta }_F$ in this calculation, isIn particular the denominator in the spatial measure is $2E_p$ in both terms.

To verify the four-dimensional Fourier transform representation, let $\tau =x^0-y^0$ and integrate the energy variable first:The Feynman i-epsilon prescription puts the positive-energy pole at $E_p-i0$ and the negative-energy pole at $-E_p+i0$. For $\tau >0$, the contour integral closes in the lower half-plane with clockwise orientation. The residue theorem multiplies the residue $i{e}^{-i{E}_p\tau }/(2E_p)$ by $-i$, giving $e^{-i{E}_p\tau }/(2E_p)$. For $\tau <0$, it closes counterclockwise above: multiplying the residue $-i{e}^{i{E}_p\tau }/(2E_p)$ by $i$ gives $e^{i{E}_p\tau }/(2E_p)$. Thus $I(\tau )=e^{-i{E}_p|\tau |}/(2E_p)$. After inserting the spatial exponential, reversing $\mathbf{p}$ in the negative-time term gives exactly the time ordering expression above. HenceThis equality is understood as a distribution limit; the same prescription gives $(\square +m^2)D_F(x-y)=-i{\delta }^4(x-y)$.

All three fields have mass $m$. In the interaction picture the quantum fields evolve with the free Hamiltonian operator, so their equations areThe interactions evolve the states instead. For $\chi \in \{\varphi ,\psi ,\Phi \}$ the appropriate mode expansion of a free field isIndependent species have commuting creation and annihilation operators. The full set of canonical commutation relations isThese follow from the separate equal-time canonical commutation relations for the three real scalar fields. For comparison, in the Heisenberg picture the full Euler-Lagrange field equations contain interactions: $(\square +m^2){\varphi }_H=-g{\varphi }_H{\Phi }_H$, $(\square +m^2){\psi }_H=-g{\psi }_H{\Phi }_H$, and $(\square +m^2){\Phi }_H=-\frac{g}{2}({\varphi }_H^2+{\psi }_H^2)$. These are distinct from the free equations used in the interaction picture.

There are no derivative couplings, so the interaction Hamiltonian is $H_I(t)=-\int d^{3}x{\mathcal{L}}_I(x)$. The interaction picture evolution operator obeysIterating its integral equation gives the Dyson seriesThe ordered integration region is equivalent to the unrestricted region with time ordering and the factor $1/n!$. Taking the scattering limit $t_0\to -\infty$, $t\to \infty$, with the usual adiabatic prescription for perturbation theory, therefore yields the S-matrixAll fields inside this expression are interaction picture fields.

Use relativistic normalization of a one-particle state and take $|i⟩=a_{\varphi }^{†}(\mathbf{p})a_{\psi }^{†}(\mathbf{q})|0⟩$, $|f⟩=a_{\varphi }^{†}({\mathbf{p}}^{\prime })a_{\psi }^{†}({\mathbf{q}}^{\prime })|0⟩$. Define the connected scattering amplitude by $⟨f|S|i{⟩}_{\mathrm{conn}}=i(2\pi {)}^4{\delta }^4(p^{\prime }+q^{\prime }-p-q)\mathcal{M}$. The identity contribution describes no scattering; the term linear in $g$ cannot connect these four external particles. At quadratic order, the connected term in the Dyson series isBy Wick theorem, a connected Wick contraction must attach the two external $\varphi$ legs to ${\varphi }_I^2$ at one vertex, the two external $\psi$ legs to ${\psi }_I^2$ at the other, and contract the remaining two ${\Phi }_I$ fields. The two assignments of species to $x,y$ cancel the Dyson factor $1/2!$. At either vertex, the $2!$ assignments of its identical fields to the incoming and outgoing particle cancel the interaction factor $1/2$. Equivalently, the coefficient $-g^2/8$ is multiplied by $2\times 2\times 2$, giving $-g^2$. Therefore, writing $D_F(x-y)=⟨0|T{\Phi }_I(x){\Phi }_I(y)|0⟩$,Insert the Feynman propagator $D_F(x-y)=\int \frac{d^{4}r}{(2\pi {)}^4}\frac{i{e}^{-ir\cdot (x-y)}}{r^2-m^2+i0}$. The two vertex integrals impose $r=p^{\prime }-p$ and $p^{\prime }+q^{\prime }=p+q$. With the Mandelstam invariant $t=(p-p^{\prime }{)}^2$, this givesThis is distinguishable scalar scattering by a shared mediator with mediator mass $m$. In physical elastic kinematics $t\le 0$, so the denominator has no pole for $m>0$. There is only the displayed $t$-channel tree-level Feynman diagram: an $s$- or $u$-channel would require an absent mixed $\varphi \psi \Phi$ vertex. Each vertex contributes $-ig$ and the exchanged $\Phi$ contributes $i/(t-m^2+i0)$, agreeing with the direct Wick contraction calculation.

Let $U(t,t_0)$ be the unitary Schrödinger propagator. The Heisenberg state is fixed at its initial value, while $A_H=U^{†}{A}_{S}U$. Since ${\psi }_S(t)=U{\psi }_H$, matrix elements and probabilities agree exactly in the two pictures. Differentiating, for an operator with no explicit Schrödinger-time dependence, givesFor constant Hamiltonian operator, $U=e^{-iH(t-t_0)/\hslash }$ and $H_H=H$.

A spatial rotation is represented by $U_R(\theta )=\exp [-i\theta n\cdot J/\hslash ]$, where $J$ is the total angular momentum operator. For a particle $J=L+S$. The commutators $[J_i,x_j]=i\hslash {ϵ}_{ijk}{x}_k$ and $[J_i,S_j]=i\hslash {ϵ}_{ijk}{S}_k$ giveThese are the correct infinitesimal vector rotations, verifying the generator and sign.

For the specified Hamiltonian operator, $p^2$ and $x^2$ commute with every $J_i$. Moreoverby antisymmetry. Thus $[J_i,H]=0$ and every total-angular-momentum component is constant in the Heisenberg picture. The separate torques areThey cancel but need not vanish individually, so orbital and spin angular momentum are generally not separately conserved.

Use natural units $\hslash =c=1$ and the Minkowski metric $g=\mathrm{diag}(1,-1,-1,-1)$. The canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Legendre transform in mechanics givesFor canonical quantization of a real scalar field, impose the equal-time canonical commutation relationsThe Heisenberg picture equation $\dot{O}=i[H,O]$ then gives $\dot{\varphi }=\pi$. For the other equation, the momentum term commutes with $\pi$; commuting the gradient term past $\pi$ and integrating the derivative of the Dirac delta function by parts gives $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$. ConsequentlyThe Fourier transform reduces this Klein-Gordon equation to oscillators of frequency $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$. Hermiticity of the real scalar field pairs their positive- and negative-frequency coefficients, so its general mode expansion can be normalized asTo verify its normalization rather than guess the oscillator algebra, use scalar field oscillator inversion at $t=0$:Its adjoint has the conjugate phase and the opposite sign in front of $i\pi$. The two mixed field-momentum commutators giveThe Fourier representation of the Dirac delta function sets the energies equal in the last step. Similarly, $[a(\mathbf{p}),a(\mathbf{q})]$ has a factor $E_{\mathbf{q}}-E_{\mathbf{p}}$ multiplying ${\delta }^3(\mathbf{p}+\mathbf{q})$ and vanishes; the two-creator commutator vanishes as well.

Substituting the expansion in $H$ and integrating the spatial phases, the two-annihilator and two-creator coefficients cancel by $E_{\mathbf{p}}^2={\mathbf{p}}^2+m^2$. The free real scalar Hamiltonian in oscillator variables isIt differs from $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}a/(2\pi {)}^3$ by its constant vacuum energy; normal ordering removes that constant. Using the oscillator commutators, for either Hamiltonian,Choose the Fock vacuum annihilated by every $a(\mathbf{p})$. Each creation operator adds a particle of energy $E_{\mathbf{p}}$ and momentum $\mathbf{p}$; its annihilation operator removes one. Their mutual commutativity gives bosonic Fock space and symmetric many-particle states. Thus the scalar mass-shell frequency has a particle interpretation, rather than merely describing a classical oscillation.