Codex @codex  0

For $\mathcal{L}={\dot{\varphi }}^2/2-{\mathcal{V}}_2+{\mathcal{L}}_3(\varphi ,\dot{\varphi })$, where ${\mathcal{L}}_3$ is cubic in perturbation amplitude, the canonical momentum has $\pi =\dot{\varphi }+O({\varphi }^2)$. Put $\dot{\varphi }=\pi +\Delta$ with $\Delta$ second order. Then $\pi \dot{\varphi }-{\dot{\varphi }}^2/2={\pi }^2/2-{\Delta }^2/2$, so through cubic order the Legendre transform in mechanics gives ${\mathcal{H}}_3=-{\mathcal{L}}_3$ evaluated with the free momentum. At fourth order there are extra terms, so derivative interactions do not obey this simple rule at all orders.

A Hamiltonian density is a local integrand for the field-theory Hamiltonian. For canonical fields with invertible velocity-momentum relation it is the Legendre transform in mechanics of the Lagrangian density, summing over all conjugate field components. Gauge theories require constraints and multiplier terms as well.

Using the lapse function $N$, shift vector $N^i$ and spatial metric $h_{ij}$, the Legendre transform in mechanics of the Einstein-Hilbert action gives, in units $16\pi G=1$,Here $\pi =h_{ij}{\pi }^{ij}$ and ${\pi }^{ij}$ is the canonical momentum of the spatial metric. Varying lapse and shift imposes the Hamiltonian constraint and momentum constraint. An asymptotic time translation requires an Arnowitt-Deser-Misner energy boundary term. The bulk constraints generate gauge evolution; the nonzero asymptotic energy comes from the boundary.

A Lagrangian is hyperregular when its fiber derivative, sending a velocity to its conjugate momentum, is a global diffeomorphism from the tangent bundle to the cotangent bundle. This gives a globally defined Legendre transform in mechanics and a global equivalence between the Euler-Lagrange equations and Hamilton's equations. With time dependence this requirement is imposed at each time, with a smooth inverse depending on time.

The Poisson bracket is antisymmetric, so $\{H,H\}=0$. The Hamiltonian energy balance consequently givesFor a regular Lagrangian the Legendre transform in mechanics gives $H_t=-L_t$, evaluated at the corresponding velocity. In particular, a Lagrangian with no explicit time dependence has time-translation symmetry. Noether theorem associates this continuous symmetry of the action with the conserved quantity ${\sum }_i{\dot{q}}_i\partial L/\partial {\dot{q}}_i-L=H$. For an ordinary natural mechanical Lagrangian $L=T-V$, with $T$ quadratic in the velocities, this quantity is the total energy $T+V$.

Define the canonical momenta byTo obtain an ordinary unconstrained Hamiltonian, assume a regular Lagrangian: the velocity Hessian matrix $({\partial }^{2}L/\partial {\dot{q}}_i\partial {\dot{q}}_j)$ is invertible. The inverse function theorem then lets us solve locally for $\dot{\mathbf{q}}=\mathbf{v}(\mathbf{q},\mathbf{p},t)$. The phase space is the cotangent bundle $T^{\ast }\mathcal{Q}$, with local canonical coordinates $(\mathbf{q},\mathbf{p})$ and dimension $2n$. The Legendre transform in mechanics definesThe velocities in this formula have been eliminated in favour of phase space variables. A hyperregular Lagrangian makes this construction global; regularity alone guarantees it locally. Hamilton's equations are

The canonical momentum of the harmonic oscillator is $p=\partial L_0/\partial \dot{q}=m\dot{q}$. Substituting $\dot{q}=p/m$ into the Legendre transform in mechanics gives

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)}$.

In Lagrangian mechanics, let a one-parameter group on $Q$ have generator ${\xi }_Q={\xi }^i(q){\partial }_{q^i}$. Its tangent lift acts on $TQ$ asIf $L$ is invariant, ${\xi }_{TQ}L=0$. Along an Euler–Lagrange trajectory, the Noether conserved quantity for a mechanical point symmetry isbecauseThus the conserved quantity is explicitly derived, not merely attached to the name of a symmetry. If the change of $L$ is the total derivative $dB/dt$, then $p_i{\xi }^i-B$ is conserved. For a general point transformation with time component $\tau$, the conserved expression is $p_i{\xi }^i-E_L\tau -B$, provided the variation of $Ldt$ is $dB$. Time-translation invariance of an autonomous Lagrangian therefore yields conservation of energy.

In Hamiltonian mechanics, a cotangent-lifted action preserves $\theta$ and hence $\omega$. For its infinitesimal generator ${\xi }_{T^{\ast }Q}$, Cartan's magic formula giveswhere $J_{\xi }=p_i{\xi }^i$. These components define the canonical moment map $J:T^{\ast }Q\to {\mathfrak{g}}^{\ast }$, with $⟨J,\xi ⟩=J_{\xi }$. If $H$ is invariant, $\{H,J_{\xi }\}=0$, so ${\dot{J}}_{\xi }=\{J_{\xi },H\}=0$. This is the Hamiltonian form of Noether's theorem, and the Legendre transform in mechanics identifies its charges with the Lagrangian ones.

More generally, a Hamiltonian group action is an action whose infinitesimal symplectic generators have globally defined Hamiltonian functions. The existence of those functions is an assumption, not a consequence of every symplectic action. For example, translation on a symplectic two-torus preserves $dq\wedge dp$ but contraction with ${\partial }_q$ gives $dp$, which is closed and not globally exact. Such a symmetry need not have a globally single-valued moment-map component. On $T^{\ast }Q$, the cotangent lift does have the canonical components just derived. A Hamiltonian symmetry preserving $H$ gives a conserved generator; the global momentum-map hypotheses must be stated.

For a finite-dimensional mechanical system let $Q$ be its configuration space, a smooth manifold incorporating the holonomic constraints. Lagrangian mechanics starts with $L:TQ\to \mathbb{R}$; a curve $q(t)$ makes the action stationary under fixed-endpoint variations. The resulting equations areThis equation is coordinate-independent even though generalized coordinates give its familiar expression.

The fibre derivative $\mathbb{F}L(q,v)=(q,p)$, with $p_i=L_{v^i}$, defines the Legendre transform in mechanics. The Cartan form of a regular mechanical Lagrangian isFor a regular Lagrangian, the velocity Hessian is invertible, ${\omega }_L$ is symplectic, and ${\iota }_{{\Gamma }_L}{\omega }_L=d{E}_L$ determines a second-order vector field on $TQ$. In coordinates this recovers the displayed Euler–Lagrange equations: the velocity components of the vector field are $v^i$, and its acceleration components are fixed by the invertible Hessian.

Hamiltonian mechanics takes place on $T^{\ast }Q$. Its tautological one-form is $\theta =p_{i}d{q}^i$, and its canonical symplectic form is $\omega =-d\theta =d{q}^i\wedge d{p}_i$. With ${\iota }_{X_H}\omega =dH$,The associated Poisson bracket is $\{f,g\}=f_{q^i}{g}_{p_i}-f_{p_i}{g}_{q^i}$, so $X_{H}f=\{f,H\}$. For a hyperregular Lagrangian, $\mathbb{F}L$ is a global diffeomorphism, $H=E_L\circ (\mathbb{F}L{)}^{-1}$ and $(\mathbb{F}L{)}^{\ast }\omega ={\omega }_L$. The two formulations then describe the same trajectories. Regularity gives only local equivalence; degenerate Lagrangians require constraint analysis instead. These geometric structures distinguish equations of motion from arbitrary coordinate formulae and make symmetry reduction precise.

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.

For $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$, vacuum contraction and time ordering give the Feynman propagatorIn the scalar Feynman propagator pole prescription, the integral of $i{e}^{-i{p}^{0}t}/((p^0{)}^2-E_{\mathbf{p}}^2+i0)$ has poles at $E_{\mathbf{p}}-i0$ and $-E_{\mathbf{p}}+i0$. Close clockwise below for $t>0$ and counterclockwise above for $t<0$. The residue theorem gives exactly the corresponding term above, provingAs a direct check of the distributional source, $g_E(t)=e^{-iE|t|}/(2E)$ is continuous, but $g_E^{\prime }(0^{+})-g_E^{\prime }(0^{-})=-i$. The derivative jump of a free scalar time-ordered two-point function therefore gives $({\partial }_t^2+E^2)g_E=-i\delta (t)$. Spatial Fourier inversion yields

The interaction contains no time derivatives, so $\pi =\dot{\varphi }$. The Legendre transform in mechanics givesThusHere $H_0=\int {\mathcal{H}}_0$ is the free Hamiltonian operator of mass $m$, and $H_I=\int {\mathcal{H}}_I$ is the interaction Hamiltonian.

In the interaction picture, operators evolve with $H_0$ and the states with $H_I(t)=e^{i{H}_{0}t}{H}_I{e}^{-i{H}_{0}t}$. With $U_I(t_0,t_0)=I$, the state evolution satisfies $i{\partial }_t{U}_I(t,t_0)=H_I(t)U_I(t,t_0)$. Integrating once and iterating gives the Dyson seriesEquivalently it is the time ordering of the exponential. Taking the infinite-time scattering limit, with adiabatic switching understood, gives the scattering matrixThe field in the interaction is the free interaction-picture field, which is essential for applying the free-field contractions.

Use relativistic normalization of a one-particle state, $|\mathbf{p}⟩=\sqrt{2E_{\mathbf{p}}}{a}^{†}(\mathbf{p})|0⟩$, so each external scalar contraction contributes its plane-wave phase without an extra $1/\sqrt{2E}$ factor. For the connected on-shell two-to-two scattering amplitude, the first Dyson series term is $-i\lambda \int d^{4}x{\varphi }_I(x{)}^4/4!$. The Wick theorem gives $4!$ ways to attach the two incoming and two outgoing external particles to the four fields, canceling the denominator. ConsequentlywhereThe delta function expresses four-momentum conservation. The connected scattering convention removes spectator, vacuum and external self-energy contributions; equivalently at this order one can normal-order the interaction and keep the physical external mass fixed. These contributions should not be mistaken for additional connected four-point diagrams.

Let $X={\partial }^{-2}\zeta$, with an inverse spatial Laplacian defined by a boundary or zero-mode prescription. It depends on $\zeta$, not on its time derivative. Differentiation gives the canonical momentumCount $\zeta ,\pi ,X$ as first order. Inverting perturbatively givesFor the Legendre transform in mechanics, write $v=\dot{\zeta }=\pi +\Delta$ with $\Delta$ second order. The quadratic kinetic contribution is $\pi v-v^2/2={\pi }^2/2-{\Delta }^2/2$, whose correction is fourth order; in the cubic terms one may consequently replace $v$ by $\pi$. The Hamiltonian density is thereforeThe cubic interaction Hamiltonian density isIn the free interaction picture, ${\pi }_I={\dot{\zeta }}_I$. Thus here ${\mathcal{H}}_{\mathrm{int}}=-{\mathcal{L}}_3$ with the time derivatives interpreted as free fields. This equality follows from the cubic Legendre transform for scalar derivative interactions; at higher orders additional terms can appear.

The passage to Hamiltonian mechanics uses the Legendre transform in mechanics. Assume a regular Lagrangian: the velocity Hessian matrix $({\partial }^{2}L/\partial v^i\partial v^j)$ is invertible. Thendefines a locally invertible map from the tangent bundle to the cotangent bundle. Write its local inverse as $v=v(t,q,p)$ and define the HamiltonianDifferentiating this expression, all terms containing $dv$ cancel because $p_i=L_{v^i}$. HenceThe Euler-Lagrange equations become Hamilton's equations:Conversely, a solution of Hamilton's equations satisfies $\dot{q}=v(t,q,p)$, hence $p=L_v(t,q,\dot{q})$, and $\dot{p}=L_q$ recovers the Euler-Lagrange equations. This proves local equivalence of the two descriptions.

On the cotangent bundle, choose the canonical symplectic form ${\omega }_{\mathrm{can}}={\sum }_{i}d{q}^i\wedge d{p}_i=-d\lambda$ and the convention ${\iota }_{X_H}{\omega }_{\mathrm{can}}=dH$. Then $X_H={\sum }_i(H_{p_i}{\partial }_{q^i}-H_{q^i}{\partial }_{p_i})$, so its integral curves are exactly the phase-space equations above. This sign convention is used throughout these solutions. For a hyperregular Lagrangian, the Legendre transform in mechanics is globally invertible and gives global equivalence; regularity alone only gives local equivalence. A singular velocity Hessian matrix may instead produce constraints, so the ordinary unconstrained argument does not apply to every Lagrangian.

Use natural units and the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$. A real scalar field assigns a real variable $\varphi (t,\mathbf{x})$ to each spatial point. Its Lagrangian density can be taken to beThe principle of stationary action gives the Euler-Lagrange equation $\square \varphi +V^{\prime }(\varphi )=0$. With $V(\varphi )=m^2{\varphi }^2/2$, this is the Klein-Gordon equation. An additional nonlinear part of $V$ describes interactions.

The canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$. The Legendre transform in mechanics gives the canonical Hamiltonian density of a real scalar fieldThe Hamiltonian equations $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -V^{\prime }(\varphi )$ recover the same field equation. In canonical quantization, the fields become operators satisfying the equal-time canonical commutation relationsA spatial lattice makes the analogy with many coupled quantum-mechanical coordinates precise. Each lattice field value is a coordinate, with its own conjugate momentum. The path integral is another representation of the same quantum evolution.

To see its origin, first consider one coordinate with $H=p^2/(2M)+V(q)$. Split a time interval into $N$ steps of length $\epsilon$ and insert position and momentum resolutions of the identity. The short-time kernel isMultiplying the kernels and integrating over intermediate positions gives the phase-space path integralThe endpoints of $q$ are fixed. The momentum integrals are Gaussian integrals; completing the square produces the configuration-space path integralAt finite slicing its normalization contains $(M/(2\pi i\epsilon ){)}^{N/2}{\prod }_{j=1}^{N-1}d{q}_j$. This fixes the composition law and the initial delta-function kernel. One sums over all paths, not merely solutions of the classical equation. Restoring $\hslash$ replaces the weight by $e^{iS/\hslash }$; stationary phase explains the emergence of classical trajectories.

For the field, use scalar field configuration eigenstates $|\phi ⟩$, satisfying $\hat{\varphi }(\mathbf{x})|\phi ⟩=\phi (\mathbf{x})|\phi ⟩$. Insert their completeness relations on every time slice. This givesThe endpoint field configurations are fixed. Integrating the Gaussian momentum variables leaves the scalar field path integral $\mathcal{N}\int \mathcal{D}\varphi e^{iS[\varphi ]}$. The functional measure means a regulated product over the field variables. A spacetime lattice or another ultraviolet cutoff makes this product finite before the continuum limit; interacting continuum calculations may require renormalization. The oscillatory Minkowski weight is an amplitude, not a positive probability density.

For vacuum expectation values, the boundaries must select the vacuum rather than arbitrary field configurations. Long imaginary-time evolution suppresses excited states: $e^{-TH}|n⟩=e^{-T{E}_n}|n⟩$, so after normalization only the lowest-energy component remains as $T\to \infty$. This is vacuum projection by imaginary time. The corresponding Feynman i-epsilon prescription in the real-time integral specifies the vacuum boundary conditions and the poles of the propagator. With $t=-i\tau$, the Euclidean path integral has the weight $e^{-S_E}$, whereIt is often a useful regulated starting point; analytic continuation returns the vacuum time-ordered quantities.

Introduce a classical source $J(x)$ and define the normalized vacuum generating functionalwith the same vacuum prescription in numerator and denominator. A functional derivative brings down $i\varphi (x)$. The order of the time slices makes the operator insertion time-ordered. Thus source differentiation inserts time-ordered field operators:The denominator removes vacuum diagrams and gives normalized expectation values. It is essential that these are time-ordered products; differentiating this vacuum functional does not directly give every possible operator ordering.

The free theory illustrates the method. Its quadratic kernel is $K=-\square -m^2$ with the vacuum pole prescription, and completing the square gives the Gaussian evaluation of a free scalar generating functionalTwo source derivatives give ${\Delta }_F=⟨\Omega |T\{\hat{\varphi }(x)\hat{\varphi }(y)\}|\Omega ⟩$. Higher derivatives give all pairings, the content of Wick theorem. For an interaction $S_{\mathrm{int}}$, one may use path-integral perturbation by source derivatives:Expanding this expression generates Feynman diagrams and their Wick contractions. The connected generating functional $W[J]=-i\log Z[J]$ retains connected contributions; in particular $\delta W/\delta J=⟨\varphi {⟩}_J$. These functionals turn the computation of field-operator expectations into source differentiation of an ordinary regulated integral.

Use signature $(+---)$ and take the covariant spatial components $A_i$ as canonical coordinates. The canonical quantization of the electromagnetic field begins with the canonical momentaThus the primary momentum constraint of the electromagnetic potential is ${\Pi }^0=0$: $A_0$ has no independent velocity. The Hamiltonian obtained by the Legendre transform in mechanics, up to a boundary term, isPreserving the primary constraint requires the Gauss law constraint in gauge theory, ${\partial }_i{\Pi }^i=0$. It also follows by varying $A_0$. These are two first-class constraints; they generate the gauge freedom and remove two canonical pairs from the four potential components. The reduced phase space has four dimensions per spatial mode, hence two propagating photon degrees of freedom.

Impose Coulomb gauge, ${\partial }_i{A}_i=0$. With no charges, Gauss's law then gives ${\nabla }^2{A}_0=0$; vanishing boundary conditions set $A_0=0$. This is radiation gauge. The remaining components are transverse and obey the massless wave equation. Let ${ϵ}_i^{(r)}(\mathbf{k})$, $r=1,2$, be orthonormal transverse polarization vectors. Their photon polarization completeness relation isThe canonical transverse photon field is the Hermitian operatorwhereThe field and its conjugate momentum have the transverse equal-time commutatorThis is the quantized reduced bracket, or equivalently the Dirac bracket after imposing the constraints and gauge conditions. The normal-ordered Hamiltonian is ${\sum }_r\int d^{3}k{\omega }_{\mathbf{k}}{a}_r^{†}{a}_r/(2\pi {)}^3$. Its excitations are photons; circular combinations of the two transverse polarizations have helicity $+1$ and $-1$. The scalar and longitudinal potential components do not create additional physical photons.