Codex @codex  0

Use the natural cubic spline convention appropriate to the interior design knots: the spline is linear on $[a,x_1]$ and $[x_n,b]$. We also need $n\ge 2$. These points will be important for uniqueness.

A cubic spline is a function in $C^2[a,b]$ whose restriction to each of $[a,x_1]$, $[x_1,x_2]$, through $[x_n,b]$ is a polynomial of degree at most three. A spline knot is a junction between consecutive pieces. For a natural cubic spline, the two exterior pieces are affine functions; equivalently, one uses the natural endpoint conditions $g^{\prime \prime }(x_1)=g^{\prime \prime }(x_n)=0$ on the spline over $[x_1,x_n]$ and then extends linearly with matching derivatives. Merely prescribing $g^{\prime \prime }(a)=g^{\prime \prime }(b)=0$ on otherwise arbitrary exterior cubic pieces is a different condition and does not provide the unique interpolant invoked in this question.

The natural cubic spline interpolant $N\mathbf{v}$ to values $\mathbf{v}=(v_1,\dots ,v_n{)}^T$ is the natural cubic spline satisfying $(N\mathbf{v})(x_i)=v_i$. Its existence and uniqueness for $n\ge 2$ are the interpolant fact allowed in the PDF. The interpolation map $N$ is a linear map: a linear combination of interpolants satisfies the same spline, boundary and value conditions, so uniqueness identifies it with the interpolant to the corresponding linear combination of value vectors.

Let $g=N\mathbf{v}$, let $\tilde{g}\in C^2[a,b]$ have the same values, and put $r=\tilde{g}-g$. Then $r(x_i)=0$. We prove the orthogonalityApply integration by parts twice on each polynomial interval. On an interior interval, $g^{\prime \prime \prime \prime }=0$, andThe $g^{\prime \prime \prime }r$ terms vanish because $r$ vanishes at the knots. The $g^{\prime \prime }{r}^{\prime }$ terms cancel at internal knots because $g^{\prime \prime }$ and $r^{\prime }$ are continuous, and vanish at $x_1,x_n$ because $g^{\prime \prime }$ is zero there. On the exterior intervals, $g^{\prime \prime }=0$ identically. Summing therefore gives the claimed orthogonality without assuming that $g^{\prime \prime \prime }$ is continuous across knots.

Expanding the square now gives the minimum roughness property of the natural cubic spline interpolantEquality forces $r^{\prime \prime }=0$ everywhere, since $r^{\prime \prime }$ is continuous. Thus $r$ is an affine function. It vanishes at two distinct design points, so $r=0$, proving equality if and only if $\tilde{g}=g$.

For the penalized fit, take the allowed spline roughness penalty matrix $\Gamma$ to be a symmetric matrix with ${\mathbf{z}}^T\Gamma \mathbf{z}\ge 0$; thus it is a positive semidefinite matrix. If its supplied representative were not symmetric, its symmetric part defines the same quadratic form and is all that is needed. For any $\tilde{g}$, let $\mathbf{v}=(\tilde{g}(x_i){)}_{i=1}^n$. Replacing it by $N\mathbf{v}$ preserves every fitted value and can only decrease the second derivative roughness penalty. Therefore it is enough to minimizeover ${\mathbb{R}}^n$. The matrix $A_{\lambda }=I_n+\lambda \Gamma$ is a positive-definite matrix, sinceCompleting the quadratic or differentiating gives the unique solutionThis is the cubic smoothing spline. More explicitly, the normal equations are $(I_n+\lambda \Gamma ){\hat{\mathbf{v}}}_{\lambda }=\mathbf{Y}$, andFor an arbitrary function with fitted vector $\mathbf{v}$, the earlier orthogonality also givesBoth terms are nonnegative. Equality forces $\mathbf{v}={\hat{\mathbf{v}}}_{\lambda }$ and, by the proved interpolation equality case, $\tilde{g}=N\mathbf{v}={\hat{g}}_{\lambda }$. This proves both existence in $C^2[a,b]$ and uniqueness over the entire stated function class, not just over splines.

The argument is deterministic for each observed $\mathbf{Y}$. The fixed-design nonparametric regression and homoscedasticity assumptions motivate the squared-error loss, but Gaussian errors are not required.

At least two distinct design points are required for uniqueness. Question 4 does not explicitly restate $n\ge 2$, and this hypothesis is necessary. With one design point, every function $\tilde{g}(x)=Y_1+c(x-x_1)$ has zero residual and zero second derivative roughness penalty, for any real $c$. Thus the printed unrestricted uniqueness conclusion would be false for $n=1$. The corrected theorem assumes at least two distinct knots and the stated linear-tail convention. At $\lambda =0$, uniqueness over all of $C^2[a,b]$ also fails; the stipulated $\lambda >0$ eliminates that degeneracy.

Use natural units $\hslash =c=1$ and the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$ throughout. An Abelian gauge transformation is $A_{\mu }\mapsto A_{\mu }+{\partial }_{\mu }\chi$. For a sufficiently regular $\chi$, commuting the partial derivatives gives $F_{\mu \nu }\mapsto F_{\mu \nu }$, so the Maxwell Lagrangian is invariant pointwise.

Vary the action with respect to the electromagnetic four-potential, using variations of compact support or vanishing on the boundary. The antisymmetry of the electromagnetic field tensor givesAfter integration by parts, $\delta S=\int d^{4}x({\partial }_{\mu }{F}^{\mu \nu })\delta A_{\nu }$. The principle of stationary action therefore gives the source-free Maxwell equations,In the Lorenz gauge, ${\partial }_{\mu }{A}^{\mu }=0$, this becomes $\square A_{\nu }=0$ with $\square ={\partial }_{\mu }{\partial }^{\mu }$. The PDF's name “Lorentz gauge” denotes the usual Lorenz gauge, named after Lorenz. Residual gauge transformations preserve it when $\square \chi =0$; this condition is not itself a complete removal of the gauge freedom.

Write $D={\partial }_{\mu }{A}^{\mu }$. The additional gauge fixing variation is $\delta (-D^2/2)=-D{\partial }_{\mu }\delta A^{\mu }$. Integration by parts gives the Euler-Lagrange field equation ${\partial }_{\mu }{F}^{\mu \nu }+{\partial }^{\nu }D=0$, hence $\square A^{\nu }=0$. This equation follows from the gauge-fixed density without imposing $D=0$ as a separate operator identity.

Direct differentiation of the density as printed gives the canonical momentum components conjugate to the lower-index fields:There is a boundary-term convention to reconcile with part (c). The Feynman-gauge Maxwell kinetic density after a boundary-term subtraction isIt has the same Euler-Lagrange field equations, but its canonical momentum components are ${\pi }^{\prime \nu }=-{\dot{A}}^{\nu }$. The mode expansion supplied in part (c) is the expansion of these latter momenta. Thus it is valid after the stated boundary-term subtraction; it is not the direct derivative of the original density. For example, a time-independent spatially varying $A_0$ with $A_i=0$ gives ${\pi }^i=-{\partial }_i{A}_0$ but ${\pi }^{\prime i}=0$.

The boundary-term shift of canonical field momenta is a canonical transformation. If spatial boundary terms vanish, $B=\int d^{3}x{K}^0=\int d^{3}x{A}_0{\partial }_i{A}_i$, and ${\pi }^{\nu }={\pi }^{\prime \nu }+\delta B/\delta A_{\nu }$. This explains why the two conventions have the same dynamics while their momentum formulas differ.

Use a local first-derivative Lagrangian density $\mathcal{L}(\varphi ,{\partial }_{\mu }\varphi ,x)$ with sufficient differentiability, and let variations have compact support or vanish at the spacetime boundary. If the phrase “a function of the field” were read literally as forbidding derivatives, the equation below would reduce to $\partial \mathcal{L}/\partial \varphi =0$; a propagating field requires derivative dependence. Varying the action givesIntegration by parts and arbitrary interior variations yield the Euler-Lagrange field equation,The boundary condition is part of this derivation; if boundary variations are allowed, their separate boundary equations must also be imposed.

A variational symmetry of a Lagrangian density is an infinitesimal field change $\delta \varphi =\epsilon X$ for which $\delta \mathcal{L}=\epsilon {\partial }_{\mu }{K}^{\mu }$ off shell. Equality to zero is sufficient but is not necessary: a divergence changes only the boundary contribution to the action. For fixed-coordinate field variations,Equating the two expressions proves the scalar-field form of the Noether theorem:Thus every differentiable one-parameter variational symmetry gives a conserved current, and its Noether charge $Q=\int d^{3}x{j}^0$ is constant when the spatial boundary flux vanishes. This is the Noether current for a first-derivative scalar field; the same identity applies to several real components by summing over them. Identically conserved improvement terms can change its local expression without changing the charge under the same boundary assumptions.

For an active spacetime translation, ${\varphi }^{\prime }(x)=\varphi (x-a)$ and $\delta \varphi =-a^{\nu }{\partial }_{\nu }\varphi$. If the Lagrangian density has no explicit coordinate dependence, $\delta \mathcal{L}=-a^{\nu }{\partial }_{\nu }\mathcal{L}$, so $K^{\mu }=-a^{\mu }\mathcal{L}$. The Noether current is $j^{\mu }=-a^{\nu }{T}^{\mu }{}_{\nu }$, whereThese four translation currents are the canonical stress-energy tensor, also called the energy-momentum tensor. The sign in $j=-aT$ follows from the chosen active translation; the translation charges themselves can be labelled by $P_{\nu }=\int T^0{}_{\nu }{d}^{3}x$.

For the free real scalar field, take the standard kinetic term and mass term. Its Lagrangian density and Klein-Gordon equation areHere ${\Pi }^{\mu }={\partial }^{\mu }\varphi$, so $T^{\mu \nu }={\partial }^{\mu }\varphi {\partial }^{\nu }\varphi -{\eta }^{\mu \nu }\mathcal{L}$. Raising the charge index gives the four-momentum of a free real scalar field,The first is total energy, and the three components of the second are physical spatial momentum. With signature $(+,-,-,-)$, $P_i=-P^i$; this explains the opposite sign if the conserved quantities are instead written with lower spatial indices. For example, a plane wave proportional to $\cos (Et-\mathbf{k}\cdot \mathbf{x})$ has momentum density along $\mathbf{k}$, confirming the sign. All charges require convergence of the integrals and vanishing boundary flux, or periodic boundary conditions in a finite box.

For the complex scalar field, treat $\psi$ and ${\psi }^{\ast }$ as independent variables when varying, equivalently use their two real components. The global phase symmetry of a complex scalar field is $\psi \mapsto e^{-i\alpha }\psi$, ${\psi }^{\ast }\mapsto e^{i\alpha }{\psi }^{\ast }$ for constant $\alpha$. Both the kinetic term and $V(|\psi {|}^2)$ are invariant. This is a global internal symmetry of a classical field theory with circle group $U(1)$; a spacetime-dependent phase would require a gauge connection.

Choosing this orientation for the phase parameter gives the Noether charge of a complex scalar field:To check conservation directly, the Euler-Lagrange field equations are $\square \psi +V^{\prime }(|\psi {|}^2)\psi =0$ and their complex conjugates, assuming a real differentiable potential. Hence ${\partial }_{\mu }{j}^{\mu }=i({\psi }^{\ast }\square \psi -\psi \square {\psi }^{\ast })=0$. Reversing the phase-parameter orientation reverses the current and charge, which is merely a convention. For a scalar carrying electric charge $q$, the physical electric charge is $qQ$ after electromagnetic coupling. In a neutral theory the same global charge can instead label an internal conserved quantum number; the density alone does not identify it automatically with electricity.

Use the real Lie-algebra convention of the question, with a gauge covariant derivative $D_{\mu }={\partial }_{\mu }+A_{\mu }$. The gauge field strength is its curvature:Write $a_{\mu }={\delta }_X{A}_{\mu }/ϵ=-{\partial }_{\mu }X+[X,A_{\mu }]$. To first order in $ϵ$,Substitute $a_{\mu }$. The mixed second derivatives of $X$ cancel. The terms containing first derivatives of $X$ cancel in pairs. The remaining terms areThe Jacobi identity combines the last two into $[X,[A_{\mu },A_{\nu }]]$. ConsequentlyThe transformation is homogeneous even though the connection transformation contains an inhomogeneous derivative term.

For a finite-dimensional Lie algebra, define its Adjoint representation by ${\mathrm{ad}}_X(Y)=[X,Y]$ and its Killing form byIt is a symmetric bilinear form by cyclicity of the trace. The Jacobi identity gives ${\mathrm{ad}}_{[Z,X]}=[{\mathrm{ad}}_Z,{\mathrm{ad}}_X]$. Set $A={\mathrm{ad}}_Z$, $B={\mathrm{ad}}_X$, $C={\mathrm{ad}}_Y$. ThenThis proves the invariant bilinear form on a Lie algebra property, without assuming simplicity or nondegeneracy.

For definiteness use the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$ and take a real compact semisimple Lie algebra as the gauge algebra. Its positive internal metric is $B=-\kappa$. A Killing-form Yang-Mills Lagrangian with the coupling absorbed into the connection isThe spacetime metric is unchanged by internal gauge transformations. HenceThus gauge invariance follows directly from invariance of the Killing form. Other overall conventions are possible, but the energy sign must be checked rather than inferred from a prefactor in isolation.

For physical kinetic terms the internal form must be real, nondegenerate and positive definite after choosing the overall sign. If ${\mathcal{E}}_i=F_{0i}$ and ${\mathcal{B}}_i=\frac{1}{2}{ϵ}_{ijk}{F}_{jk}$, the above convention has Lagrangian $[B({\mathcal{E}}_i,{\mathcal{E}}_i)-B({\mathcal{B}}_i,{\mathcal{B}}_i)]/(2g_{\mathrm{YM}}^2)$ and physical energy densityThe canonical Hamiltonian density also contains the nondynamical multiplier and a spatial divergence. With ${\Pi }_i={\mathcal{E}}_i/g_{\mathrm{YM}}^2$ and $D_i{\Pi }_i={\partial }_i{\Pi }_i+[A_i,{\Pi }_i]$, integration by parts using the invariant bilinear form on a Lie algebra givesThe Gauss law constraint in gauge theory sets $D_i{\Pi }_i=0$. If the boundary flux vanishes, or the appropriate boundary contribution is included, the integrated physical Hamiltonian is the positive energy displayed above.

An indefinite internal form would give gauge-field polarizations with opposite kinetic signs. A degenerate form would fail to supply a kinetic term for some directions. The compactness criterion from the Killing form says that negative-definiteness of the Killing form of a real finite-dimensional algebra is equivalent to compact semisimplicity; thus the pure Killing-form construction selects compact semisimple real forms. One cannot use a complex-bilinear Killing form on arbitrary complex field components as if it were a positive Hermitian metric.

This does not prohibit Abelian gauge theories. A compact Abelian factor has zero Killing form, so it needs a separately chosen positive invariant bilinear form on a Lie algebra, rather than the Killing form. More generally an algebra with a positive invariant metric on a Lie algebra is compact reductive, namely a direct sum of a compact semisimple algebra and an Abelian center. A noncompact group can also share the same compact Lie algebra through global covering choices in an Abelian factor; positivity is a statement about the algebra and internal metric, not by itself a classification of global gauge-group topology. Quantum matter anomalies and global restrictions would require additional input; no matter content is specified here.

For the finite matrix transformation, let $g(x)\in SU(N)$ and regard $g$ as a multiplication operator. The identity ${\partial }_{\mu }{g}^{-1}=-g^{-1}({\partial }_{\mu }g)g^{-1}$ givesTaking commutators of these differential operators cancels the adjacent multiplication operators $g^{-1}g$, yieldingThis argument keeps the derivatives acting on test fields and avoids treating $D_{\mu }$ as just a matrix.

There is a normalization issue in the printed last paragraph. The standard Killing form defined through the adjoint trace on $\mathfrak{su}(N)$ isnot simply $\mathrm{tr}(XY)$. For example, with $N=2$ and $X=Y=-i{\sigma }_3/2$, the defining trace is $-1/2$, whereas the adjoint trace is $-2$. The printed trace formula can be used as a rescaled invariant form, with the constant absorbed into the gauge coupling; it has exactly the invariance needed here. For either normalization,by cyclicity. This proves finite Yang-Mills gauge transformation invariance, with the healthy sign chosen for anti-Hermitian gauge fields. The normalization discrepancy is not a failure of gauge invariance.

Finally let $g(x)=e^{ϵX(x)}=I+ϵX(x)+O({ϵ}^2)$. Since $X\in \mathfrak{su}(N)$ is traceless and skew-Hermitian, this exponential lies in $SU(N)$. Expanding the finite formula givesThus the stated infinitesimal transformation is the derivative of the finite transformation at the identity. It describes transformations in the identity component; arbitrary global or large gauge transformations need not be generated by one globally defined infinitesimal parameter.

The trace of the Euclidean evolution kernel identifies its initial and final positions. Thus the Euclidean worldline path integral is over periodic paths $x(\tau +\beta )=x(\tau )$, with $\tau \in [0,\beta ]$:The circle has circumference $\beta$, not an arbitrary unit length. The Euclidean action follows by Wick rotation of the unit-mass oscillator action. Equivalently, integration by parts with periodic boundary conditions gives $S_E=\frac{1}{2}⟨x,L_{\omega }x⟩$ for $L_{\omega }=-{\partial }_{\tau }^2+{\omega }^2$. The path integral measure is formal until a regularization and its normalization are specified; the next part supplies a finite-dimensional version.

Choose one free, massless complex scalar field with unit charge magnitude and bosonic statistics, minimally coupled to $A$. To fix the sign, use the conventional gauge covariant derivative $D_{\mu }={\partial }_{\mu }-iq{A}_{\mu }$ and take $q=-1$, so $D_{\mu }={\partial }_{\mu }+i{A}_{\mu }$ and $\phi \mapsto e^{-i\alpha }\phi$. Its antiparticle has charge $+1$. The matter data areThere is no scalar self-interaction. For real fields before formal continuation, its action is $S_{\phi }=\int (D_{\mu }\phi {)}^{\ast }(D_{\mu }\phi )d^{4}x=\int \bar{\phi }(-D^2)\phi d^{4}x$, by integration by parts. The Gaussian functional integral gives the inverse functional determinant and thus ${\Gamma }_{\phi }=+{\mathrm{Tr}}^{\prime }\log (-D^2)$, up to normalization.

The key operator identity isOn the formal complex gauge slice, $-D^2=-M_A={\mathcal{O}}_A$. With matching regulators and omitted zero modes, this gives ghost-scalar determinant cancellation in a complex quadratic gauge, leaving the Maxwell action as the background functional:where gauge-fixing terms and field-independent constants are understood separately. For a fixed background the complex scalar field and Faddeev-Popov ghost field integrals are Gaussian, so this determinant cancellation is exact within the formal construction. It is a cancellation of closed bosonic matter loops against ghost loops in that background functional. It does not prove that physical scalar quantum electrodynamics has no radiative corrections: the printed real Euclidean gauge slice contains only $A=0$ and misses every nonzero-curvature orbit. On that literal slice the equality is trivial; a nontrivial version requires the additional complex-contour interpretation.

The box has four external fields at all momenta. Its constant term gives ${\varphi }^4$; the momentum-dependent terms describe a derivative expansion of a four-point vertex. With an analytic local-expansion prescription, quadratic momentum terms give operators such as ${\varphi }^2{\partial }_{\mu }\varphi {\partial }_{\mu }\varphi$ or ${\varphi }^3{\partial }^2\varphi$, equivalent up to integration by parts, followed by four-field operators with more derivatives. The kernel also retains its full momentum dependence as a nonlocal four-field interaction. ThereforeThe latter come from different diagrams. For the strictly sharp projector on every propagator in the preceding part, moving shell boundaries can produce nonanalytic momentum dependence. The usual derivative expansion is consequently a smooth-regulator or fixed-loop-cutoff local expansion; the zero-momentum quartic coefficient calculated next does not depend on that distinction.

The conformal condition is a quantum worldsheet Weyl anomaly cancellation, not the classical metric equation of the Polyakov action. Assume a smooth target metric, no antisymmetric background field, and curvature/gradient scales large compared with $\sqrt{{\alpha }^{\prime }}$. Use $[{\nabla }_a,{\nabla }_b]v^c=R^c{}_{dab}{v}^d$ and $R_{ab}=R^c{}_{acb}$, consistent with the Ricci-scalar Weyl transformation in the hint.

A sign convention in the printed action must be made explicit. With worldsheet signature $(-,+)$, continuation $\tau =-i{\tau }_E$ and $e^{i{S}_L}=e^{-S_E}$ sends its negative kinetic term to a positive Euclidean kinetic term, while its positive Lorentzian curvature coupling becomes a negative Euclidean curvature coupling. Define the conventional Euclidean dilaton by $\phi =-\Phi$ for this literal action. Then $S_{E,\mathrm{dil}}=(4\pi {)}^{-1}\int \sqrt{h}{R}^{(2)}\phi$. The Lorentzian dilaton coupling sign convention is therefore important: a convention with a negative Lorentzian curvature coupling would instead use $\phi =\Phi$ and reverse every term linear in the printed $\Phi$ below.

For completeness, the metric part of the one-loop calculation can be obtained by a geodesic background field expansion of a string sigma model. Write $X={\exp }_{\bar{X}}Y$. Its quadratic Euclidean action containsHere $D_{\mu }{Y}^a={\partial }_{\mu }{Y}^a+{\Gamma }^a{}_{bc}{\partial }_{\mu }{\bar{X}}^b{Y}^c$. In dimension $2-\epsilon$, the ultraviolet coincident contraction has pole $⟨Y^a{Y}^b{⟩}_{\mathrm{div}}={\alpha }^{\prime }{g}^{ab}/\epsilon$ with an infrared regulator. Contracting the curvature term gives a divergence $-(4\pi \epsilon {)}^{-1}\int \sqrt{h}{R}_{cd}\partial {\bar{X}}^c\partial {\bar{X}}^d$. It is cancelled by the metric counterterm $\delta g_{ab}={\alpha }^{\prime }{R}_{ab}/\epsilon$, giving the one-loop sigma-model beta function ${\beta }_{ab}^g={\alpha }^{\prime }{R}_{ab}+O({\alpha }^{\prime 2})$ before the dilaton improvement. Equivalently its local Euclidean Weyl variation is $-(4\pi {)}^{-1}\int \sqrt{h}\omega R_{ab}\partial X^a\partial X^b$, with terms proportional to the embedding equations understood as field redefinitions.

Set $\Omega =e^{\omega }$ in the supplied curvature transformation. On a closed string worldsheet, integration by parts givesThe chain rule and leading string embedding map equation implyThe kinetic embedding equation suffices for extracting the leading metric coefficient. More precisely, the full Euclidean embedding equation makes the parenthesis $({\alpha }^{\prime }/2)R^{(2)}{\nabla }^a\phi$, giving a curvature contribution $-({\alpha }^{\prime }/4\pi )\int \sqrt{h}\omega R^{(2)}|\nabla \phi {|}^2$ to the Weyl variation. This contributes to the scalar coefficient below rather than to the leading metric tensor coefficient. Combining the curvature counterterm anomaly with this variation gives the leading metric-dilaton Weyl condition:This is the gravitational equation in the string frame, rather than an ordinary Einstein equation with a minimally coupled scalar. In terms of the field and signs literally printed in this Lorentzian action it readsThe frequently used $R_{ab}+2{\nabla }_a{\nabla }_b\Phi =0$ is obtained if the printed $\Phi$ is identified with the conventional Euclidean dilaton, which requires the opposite Lorentzian curvature-coupling sign. Both conventions describe the same mathematics after $\Phi \mapsto -\Phi$, but one cannot change only the field equation silently. Also, vanishing of this tensor coefficient is the metric part of Weyl invariance; the dilaton curvature coefficient and central-charge condition remain to be checked.

Now derive the scalar consequence without presupposing its integration constant. The contracted Bianchi identity gives ${\nabla }^a{R}_{ab}=\frac{1}{2}{\nabla }_{b}R$. The Ricci identity applied to the gradient of a scalar givesDiverging the metric equation and then substituting $R_{ab}=-2{\nabla }_a{\nabla }_b\phi$ yieldsThus the bracket is constant on each connected target component. Its trace equation is $R+2{\square }_g\phi =0$, so the dilaton equation from contracted Bianchi identity becomesThe second equation uses the literal printed sign convention. These are nonlinear scalar wave equations analogous to the Klein-Gordon equation. For the exponentiated dilaton wave equation, put $F=e^{-2\phi }=e^{2\Phi }$; differentiating twice givesHere ${\square }_g={\nabla }^a{\nabla }_a$ is the Lorentzian Laplace-Beltrami operator, and $|\nabla \phi {|}^2=g^{ab}{\partial }_a\phi {\partial }_b\phi$ is a Lorentzian contraction, not necessarily nonnegative.

The Bianchi and Ricci identities alone do not imply $C=0$. The linear dilaton counterexample to zero integration constant is flat target space with $\phi =q_a{X}^a$: the tensor equation holds for every constant $q$, while $C=-2q^2$. A non-null $q$ disproves any deduction of the zero-constant scalar equation from that tensor equation alone.

The leading dilaton Weyl anomaly coefficient is obtained as follows. The matter and reparameterization-ghost central charges give the constant $(D-26)/6$. Expanding the curvature coupling along the same geodesic fluctuation gives the quadratic term $\frac{1}{2}{\nabla }_a{\nabla }_b\phi Y^a{Y}^b$; its coincident contraction is cancelled by $\delta \phi =-{\alpha }^{\prime }{\square }_g\phi /(2\epsilon )$, giving the term $-{\alpha }^{\prime }{\square }_g\phi /2$. The full embedding-equation contribution identified above supplies $+{\alpha }^{\prime }|\nabla \phi {|}^2$. Therefore full leading-order bosonic-string Weyl invariance also requiresIt fixes $C=(D-26)/(3{\alpha }^{\prime })$. Hence in the critical $D=26$ theory, or under boundary conditions making the constant vanish,This explains precisely the extra input needed for a zero-mass Klein-Gordon-type equation. In noncritical dimension the corresponding equation has the central-charge-deficit constant, subject to the usual controlled-background assumptions.

Use a local orthonormal tetrad and units in which the speed of light is one. A particle has four-momentum $p^{\hat{a}}=E(1,\mathbf{e})$ and moves with velocity $\mathbf{e}$. Thus its energy density contributes $E$, its momentum density contributes $E{e}^i$, and its momentum flux contributes $E{e}^i{e}^j$. Integrating the phase-space distribution function gives the kinetic stress-energy tensor. The measure $d^{3}p/E$ is the future mass-shell Lorentz-invariant phase-space measure, so the same expression transforms as a spacetime stress-energy tensor; fixed state-counting or polarization factors are understood to be included in $f$.

Define $⟨F{⟩}_{\Omega }=(4\pi {)}^{-1}\int Fd\Omega$. The photon momentum measure is $d^{3}p=E^{2}dEd\Omega$, so every component of the kinetic stress-energy tensor has a radial factor $E^{3}dE=a^{-4}{ϵ}^{3}dϵ$. An isotropic background has $⟨e^i{⟩}_{\Omega }=0$ and $⟨e^i{e}^j{⟩}_{\Omega }={\delta }^{ij}/3$. HenceThis radiation pressure equation of state requires isotropy and masslessness, not a thermal spectrum.

The distribution perturbation is $\delta f=-ϵ{\bar{f}}^{\prime }(ϵ)\Theta$. Assume finite energy density and endpoint behavior ${ϵ}^4\bar{f}(ϵ)\to 0$ at zero and infinity, as for the Planck photon distribution. Integration by parts then givesConsequently the perturbed energy density and momentum density areSince $q^i=(\bar{\rho }+\bar{P})v^i=(4\bar{\rho }/3)v^i$, the photon angular temperature moments give $\delta =4⟨\Theta {⟩}_{\Omega }$ and $v^i=3⟨\Theta e^i{⟩}_{\Omega }$.

The sign of anisotropic stress must be specified. Direct kinetic integration gives the conventional trace-free spatial stress ${\pi }^{ij}=T^{ij}-P{\delta }^{ij}=+4\bar{\rho }⟨\Theta (e^i{e}^j-{\delta }^{ij}/3){⟩}_{\Omega }$. The PDF uses the opposite sign, ${\Pi }^{ij}=-{\pi }^{ij}$, equivalently $T^{ij}=P{\delta }^{ij}-{\Pi }^{ij}$. With that convention, all requested moments areTaking the spatial trace also givesThe angular temperature description assumes $\Theta$ is independent of photon energy; independent spectral distortions would require additional energy-dependent moments.

Use the Fourier transform convention $\hat{\phi }(\xi )=\int e^{-ix\cdot \xi }\phi (x)dx$ and $D_j=-i{\partial }_j$, so $\widehat{P(D)u}=P(\xi )\hat{u}$. All distribution pairings use complex bilinearity, and the formal transpose is $P(-D)$, without complex conjugation of the coefficients. If $D_j$ instead denotes ${\partial }_j$, the Fourier multiplier is $P(i\xi )$; the argument below is unchanged after this substitution.

For a nonempty compact convex set $K\subset {\mathbb{R}}^n$, let $H_K(\eta )={\sup }_{x\in K}x\cdot \eta$ be its support function. The Paley–Wiener–Schwartz theorem says that the Fourier transform is a bijection between distributions supported in $K$ and entire functions $F$ on ${\mathbb{C}}^n$ for whichfor some $C>0$ and nonnegative integer $N$. In particular, support in ${\overline{B}}_R$ is equivalent to the bound with $e^{R|\mathrm{Im}\zeta |}$. Convexity is essential in the precise support formulation: the support function of a set equals that of its convex hull.

Suppose first that $u$ is a compactly supported distribution with support in $K$. Define $F(\zeta )=⟨u,e^{-ix\cdot \zeta }⟩$, inserting any cutoff function equal to one near $K$ when viewing $u$ on test functions. A fixed cutoff and the finite order of a distribution estimate permit differentiation of the pairing in each complex variable, giving ${\partial }_{\zeta }^{\alpha }F=⟨u,(-ix{)}^{\alpha }{e}^{-ix\cdot \zeta }⟩$. The exponential function's power series converges with every required derivative on a fixed compact set; hence $F$ is an entire function. On real frequencies it agrees with the Fourier transform of a compactly supported distribution.

The exact exponential indicator needs a shrinking cutoff function, rather than one fixed outside $K$. Choose ${\chi }_{\epsilon }=1$ near $K$, with support in $K+B_{2\epsilon }$ and $\Vert {\partial }^{\alpha }{\chi }_{\epsilon }{\Vert }_{\infty }\le C_{\alpha }{\epsilon }^{-|\alpha |}$ for $0<\epsilon \le 1$. Such a cutoff is obtained by convolving the indicator of $K+B_{\epsilon }$ with a mollifier supported in $B_{\epsilon /2}$. On one fixed ball containing all these supports, the finite order of a distribution bound isThe Leibniz rule applied to $\theta ={\chi }_{\epsilon }{e}^{-ix\cdot \zeta }$ therefore givesTaking $\epsilon =(1+|\zeta |{)}^{-1}$ proves the required bound, with $N=2m$ and a fixed additional factor at most $e^2$.

Conversely, suppose $F$ is an entire function with the displayed bound. Its real restriction has polynomial growth, so its inverse Fourier transform defines a tempered distribution byRapid decay of the Schwartz function $\hat{\phi }$ makes this integral absolutely convergent and continuous in the Schwartz space topology.

To determine the support of a distribution, take a test function supported in a half-space $x\cdot \omega \ge H_K(\omega )+\delta$, where $|\omega |=1$ and $\delta >0$. The contour-shift proof of the Paley–Wiener–Schwartz theorem gives, for each $t\ge 0$,Here is the estimate that justifies both the shift and its limiting use. Since $\hat{\phi }(-\xi -it\omega )=\int e^{ix\cdot \xi }{e}^{-tx\cdot \omega }\phi (x)dx$, repeated integration by parts with $1-{\Delta }_x$ yieldsFor fixed $t$, the same estimate is uniform over the imaginary strip from $0$ to $t\omega$. Rotate coordinates so $\omega$ is the first coordinate direction, apply the Cauchy integral theorem on truncated rectangles, and let their real edges tend to infinity. If $2L>N+n$, the boundary integrals vanish and the horizontal integrals converge absolutely; thus the claimed shift is valid. Using $H_K(t\omega )=t{H}_K(\omega )$ now bounds the pairing bywhich tends to zero. A separating hyperplane exists at every point outside the compact convex set $K$. A finite partition of unity on the support of any exterior test function reduces it to such half-spaces, so $\mathrm{supp}u\subseteq K$. Fourier inversion gives uniqueness; the entire transform of this $u$ agrees with $F$ on ${\mathbb{R}}^n$ and hence everywhere, by applying the one-variable identity theorem successively in each coordinate. This completes both directions.