Codex @codex  0

For either scalar field $f$ with mass $m_f$, the supplied covariant equation in the FLRW metric is the Klein-Gordon equationSince $\dot{f}=f^{\prime }/a$, $\ddot{f}=(f^{\prime \prime }-\mathcal{H}{f}^{\prime })/a^2$ and $H=\mathcal{H}/a$, its conformal-time form isThe homogeneous scalar field backgrounds satisfy the same equations without spatial derivatives. After subtracting them, let $\delta \varphi =u/a$ and $\delta \chi =v/a$. The key derivatives areBecause $a^{\prime \prime }/a={\mathcal{H}}^{\prime }+{\mathcal{H}}^2$, the expansion-friction term cancels after this canonical rescaling of a scalar cosmological perturbation. The fluctuation equations areThe quadratic potential has no mixed derivative, and the stipulated neglect of metric perturbations removes gravitational mixing. Thus these two rescaled fluctuations evolve independently in this approximation.

The given symmetric Fourier transform convention sends $-{\nabla }^2$ to $k^2$, where $k=|\mathbf{k}|$. Hence the Fourier-mode equations areThey are harmonic-oscillator equations with time-dependent squared frequencies; a negative effective squared frequency for a rescaled field does not itself imply growth of the physical field $v/a$.

An electromagnetic field tensor obeys the homogeneous Maxwell equations $dF=0$. Together with its invariance under the Killing vector field $V$, Cartan's magic formula gives $d({\iota }_{V}F)={\mathcal{L}}_{V}F-{\iota }_{V}dF=0$. The local potential of a closed differential one-form construction, the one-form Poincaré lemma, therefore supplies a scalar field $\Phi$ on a sufficiently small contractible neighborhood withFor an explicit local proof, write ${\iota }_{V}F={\alpha }_i(x)d{x}^i$ on a coordinate ball centered at zero and take $\Phi (x)={\int }_0^1{x}^i{\alpha }_i(sx)ds$. Since $d\alpha =0$, differentiating under the integral gives ${\partial }_j\Phi ={\int }_0^1[{\alpha }_j(sx)+s{x}^i{\partial }_i{\alpha }_j(sx)]ds={\alpha }_j(x)$. There need not be a global potential when this closed one-form has nonzero periods. The homogeneous Maxwell equation is essential; an arbitrary invariant two-form would not suffice.

Use proper time for the massive particle, let $u$ be its tangent, and assume $m\ne 0$. The Killing equation makes $u^a{u}^b{\nabla }_a{V}_b=0$. The Lorentz force equation then impliesConsequently the local conserved quantity isThe sign follows from contracting $F$ in its first argument, exactly as in the source. Adding a constant to $\Phi$ merely shifts the conserved quantity. This is a charged-particle conserved quantity from a Killing symmetry.

When the number-changing reaction is in chemical equilibrium, its chemical balance condition is $2{\mu }_{\varphi }=3{\mu }_{\varphi }$. Therefore zero chemical potential from number-changing reactions gives $\boxed{{\mu }_{\varphi }=0}$. Elastic scattering maintains kinetic equilibrium but imposes no independent condition on ${\mu }_{\varphi }$. The conclusion requires the number-changing rate to remain large compared with the Hubble parameter; merely listing an allowed reaction does not guarantee this after freeze-out.

Before thermal decoupling in cosmology, the sectors have a common temperature $T$. For relativistic species at that temperature their entropy densities aresoThese are the effective number of relativistic degrees of freedom for entropy. The paper's $g_{\ast }$ agrees with $g_{\ast s}$ for relativistic components at a common temperature. For one real scalar field, $g_{\ast s}^{\varphi }=1$. At temperatures where all Standard Model species are relativistic, both counts are constant, and so is $\xi$ before decoupling; thresholds can change this equilibrium ratio if they occur while the sectors remain coupled.

After decoupling, adiabatic isolation gives separately conserved visible and dark entropy: $s_{\mathrm{SM}}{a}^3$ and $s_{\varphi }{a}^3$ are individually constant. Thus $\boxed{\xi (T<T_d)=\xi (T_d)}$ even when either sector's entropy degrees of freedom change. Changes in particle content then alter their temperature ratio rather than the ratio of their comoving entropies. This assumes no late inter-sector transfer or entropy injection.

Use the dark-sector temperature in every dimensionless integration variable: $x=m_{\varphi }/T_{\varphi }$ and $y=p/T_{\varphi }$. Assume a single real scalar field, $g=1$, and fast number-changing reactions keeping ${\mu }_{\varphi }=0$. For $x\gg 1$, the populated momenta satisfy $y=O(\sqrt{x})$, soThis is the nonrelativistic Maxwell-Boltzmann distribution limit of the Bose-Einstein distribution. The given Gaussian Gamma function integral givesSince ${\rho }_{\varphi }=m_{\varphi }n+3n{T}_{\varphi }/2+\cdots$ and $P_{\varphi }=n{T}_{\varphi }$, the entropy density at zero chemical potential is $s_{\varphi }=n(x+5/2+\cdots )$. Its leading nonrelativistic entropy at zero chemical potential is thereforeAn internal multiplicity $g$ would multiply this expression; the absence of $g$ in the printed result assumes a real one-state scalar.

Now use the constant separately conserved visible and dark entropy ratio and the relativistic visible-bath entropy law, $s_{\mathrm{SM}}=(2{\pi }^2/45)g_{\ast s}^{\mathrm{SM}}{T}_{\gamma }^3=\xi s_{\varphi }$. Since $T_{\varphi }=m_{\varphi }/x$,Thus the cannibal-sector temperature ratio hasThe numerical coefficient is approximately $0.525$. The visible-bath temperature is the photon temperature while that bath is internally thermalized; after distinct visible species acquire different temperatures, its total entropy must instead use the appropriate temperature-weighted $g_{\ast s}^{\mathrm{SM}}$.

The displayed relativistic scalar entropy formula must not be continued into this regime with $g_{\ast s}^{\varphi }=1$ held constant. If it is used as an effective definition at arbitrary temperature, thenand $T_{\gamma }/T_{\varphi }=(\xi g_{\ast s}^{\varphi }/g_{\ast s}^{\mathrm{SM}}{)}^{1/3}$ gives the same result. This specifies how the relativistic expression can legitimately be combined with the nonrelativistic one.

Let $M_{\mathrm{Pl}}$ be the reduced Planck mass, $M_{\mathrm{Pl}}^{-2}=8\pi G$ in units $\hslash =c=1$. Combining that equation with the Friedmann acceleration equation givesDifferentiate $3M_{\mathrm{Pl}}^2{H}^2=\rho$ to obtain $\dot{\rho }=-3H(\rho +P)=-3H{\dot{\varphi }}^2$. For the homogeneous canonical scalar field, $\dot{\rho }=\dot{\varphi }\ddot{\varphi }+V_{,\varphi }\dot{\varphi }$, soWhere $\dot{\varphi }\ne 0$, division gives the inflaton equation of motion; continuity extends it to isolated turning points:The Hubble friction term describes loss of kinetic energy into expansion.

There is a small logical caveat if the field is frozen on an entire interval: the homogeneous Einstein constraints alone only give the product equation and do not force $V_{,\varphi }=0$. The frozen scalar fields and the Friedmann constraint counterexample makes this explicit. For a canonical scalar, the field equation also follows directly from its inflaton action in an expanding universe, whose homogeneous Lagrangian is $a^3[{\dot{\varphi }}^2/2-V(\varphi )]$. The Euler-Lagrange equation gives $d(a^3\dot{\varphi })/dt+a^3{V}_{,\varphi }=0$, including the frozen case, and yields the same boxed equation without division by $\dot{\varphi }$.

Use $\nu =R{e}^{-1}>0$, $\kappa =P{e}^{-1}>0$, $\mathbf{U}=\overline{\mathbf{u}}+\mathbf{u}$, and $\Theta =\overline{\theta }+\theta$. Work along a sufficiently regular direct trajectory satisfying the given constraints. The zero spatial mean of the scalar field is conserved by incompressible flow, impermeable walls and zero scalar flux, so minimizing $J=(\Theta (T),\Theta (T))$ is equivalent to minimizing scalar variance, up to the fixed domain volume.

The printed functional fixes the initial state to a candidate ${\mathbf{u}}_0$; it contains no term that enforces its kinetic energy. For the optimization over that candidate, add the real Lagrange multiplier constraintEquivalently, one can restrict all control variations to the sphere in a normed vector space of fixed kinetic energy. This term changes the initial-control optimality condition, not the interior adjoint equations.

Let $\mathbf{v}=\delta \mathbf{u}$, $\sigma =\delta \theta$ and $\pi =\delta p$. Linearization of the momentum and scalar transport residuals givesBoth appearances of the perturbation velocity in the nonlinear momentum term have been differentiated. In particular, the coefficient is the gradient of the total velocity, not just the base shear.

For the negative-constraint convention of the functional, integration by parts gives the interior coefficients of $\mathbf{v},\sigma ,\pi$ asThus the adjoint equations for Boussinesq scalar mixing areThe transpose is essential: the $j$th component of $(\nabla \mathbf{U}{)}^T{\mathbf{u}}^{†}$ is ${\sum }_i{u}_i^{†}{\partial }_j{U}_i$. The coupling $-{\theta }^{†}\nabla \Theta$ transposes advection of the scalar by a velocity perturbation; $-R{i}_B{u}_y^{†}$ transposes buoyancy feedback. Dropping the latter would give a passive scalar adjoint, not the active scalar problem.

These equations are integrated backward, not forward. If $\tau =T-t$, they readwith direct coefficients evaluated at $t=T-\tau$. Both terms from the diffusion equation now have the usual forward sign in $\tau$. A direct-adjoint looping method stores or reconstructs the forward trajectory, solves these equations backward, and uses the initial adjoint as the control gradient. The endpoint and fixed-energy conditions below give necessary conditions for a local optimizer, not a global optimality theorem.

The temporal integration by parts terms areThe terminal states are free and the objective has no explicit terminal-velocity dependence. ThereforeThe factor two follows from the displayed objective without a $1/2$ prefactor. The initial scalar field perturbation is fixed, so $\delta \theta (0)=0$ and there is no independently prescribed adjoint-scalar condition at $t=0$. Its value there is obtained by backward integration. Independent variation of the initial velocity state gives ${\mathbf{u}}^{†}(0)={\mathbf{u}}_0^{†}$.

Variation of the candidate ${\mathbf{u}}_0$ in the energy-augmented functional gives the fixed-energy initial-condition optimality condition. For the usual divergence-free no-slip control space $H$, write $P_H$ for its orthogonal projection, the Leray-Helmholtz projection in the usual incompressible velocity space. This is the standard formal control condition when pressure is determined by incompressible flow and normal momentum. If the independently imposed direct pressure flux is retained as an additional constraint, the admissible control variations must also satisfy its compatibility conditions along the trajectory; they need not fill the usual space $H$. The standard projection formula alone does not prove stationarity for that more restricted problem. ThenFor the natural solenoidal adjoint space, $P_H{\mathbf{u}}^{†}(0)={\mathbf{u}}^{†}(0)$. Equivalently, its component tangent to the sphere in a normed vector space of fixed kinetic energy vanishes:Here $E_0>0$ is needed for the sphere in a normed vector space to be regular and for division by $2E_0$. The multiplier is real and can have either sign; normalizing the initial adjoint with a prescribed positive sign is not a general necessary condition for minimization. If $E_0=0$, the only feasible initial velocity is zero and this tangent formula for the sphere in a normed vector space is inapplicable. First-order stationarity alone also allows maxima or saddles; a local minimum requires the appropriate nonnegative constrained second variation.

The spatial boundary terms vanish with periodic adjoints in $x,z$, homogeneous no-slip boundary conditions ${\mathbf{u}}^{†}=0$ at $y=\pm 1$, and ${\partial }_y{\theta }^{†}=0$ there. To see this, velocity variations vanish at the wall while their normal derivatives need not, so the viscous boundary term forces the adjoint velocity to vanish. Scalar variations have zero normal derivative but free values, forcing the adjoint scalar's normal derivative to vanish. Pressure variation gives adjoint incompressible flow; no independent terminal or initial datum is assigned to $p^{†}$. Its additive time-dependent constant may be fixed by a zero spatial mean.

In particular, a homogeneous Neumann condition for the direct pressure does not by variational transposition require ${\partial }_y{p}^{†}=0$. For a smooth adjoint solution, the wall-normal adjoint momentum equation instead suppliesat the flat walls, because the direct scalar normal derivative is zero there. A zero adjoint-pressure derivative would be an additional compatibility choice, valid only when this right-hand side vanishes.

There are also literal direct-data compatibility qualifications in the printed setup. The initial total scalar field is $\overline{\theta }=-\mathrm{erf}(30y)$, whose wall derivative is $-60e^{-900}/\sqrt{\pi }\ne 0$. It is extraordinarily small but mathematically not zero. Thus exact initial scalar data and exact zero wall flux are incompatible for a classical solution smooth at $t=0$. One may use a parabolic mild solution of an abstract Cauchy problem with the boundary condition enforced for $t>0$, or replace the initial profile by an exactly compatible smooth one; these are conventions, not an equality $e^{-900}=0$. The reference scalar field is also not a stationary profile of the diffusion equation at finite $Pe$: the full $\Delta \Theta$ term must be retained.

The pressure compatibility at a no-slip wall is similarly important. Direct normal momentum gives ${\partial }_{y}p=\nu \Delta u_y-R{i}_B\theta$. The separately imposed zero pressure derivative requires this right side to vanish. It is not automatic for arbitrary fixed-energy initial velocities. For example, the divergence-free no-slip field generated by the periodic stream function $S=A(1-y^2{)}^2\sin (\alpha x)$ with $\alpha =\pi /L_x>0$ has $u_x={\partial }_{y}S$, $u_y=-{\partial }_{x}S$ and $\Delta u_y{|}_{y=\pm 1}=-8A\alpha \cos (\alpha x)\ne 0$, while $\theta (0)=0$. With $A\ne 0$ adjusted to any positive energy, no classical solution smooth up to the initial wall can satisfy that extra pressure condition. The displayed adjoints are the formal necessary equations along admissible smooth trajectories; a well-posed physical formulation normally determines pressure from incompressible flow and normal momentum, rather than imposing independent homogeneous pressure flux for every control.

For every interacting $n\ge 3$, the one-loop scalar effective action already contains counterexamples. On a constant background, with a positive auxiliary mass $\mu$ to control infrared divergences, its field-dependent part isThe term of order ${\lambda }^V$ has $E=(n-2)V$ external scalar fields and is a one-loop polygon one-particle-irreducible Feynman diagram. Choose $V>n/(n-2)$ and $2V>d$. Then $E>n$, its loop momentum integral is ultraviolet convergent, and its coefficient is nonzero. Removing the ultraviolet cutoff does not remove this finite higher-point scalar vertex. For $d>2$, generic nonexceptional external momenta give the same conclusion without retaining the auxiliary mass. In $d\le 2$, the massless theory still needs an infrared prescription; removing the ultraviolet cutoff does not remove that need or the induced vertices. Hence there are no interacting pairs $(n,d)$ under the literal wording. For $n=2$ the functional determinant is field independent; $n=1$ merely shifts a Gaussian integral when an infrared prescription exists.

There is a different conventional interpretation: absence of new independent divergent counterterms. The superficial degree of divergence of a connected diagram isThe standard perturbative counterterm closure criterion is $[\lambda ]\ge 0$, with the same pairs as in part ii, provided one includes all symmetry-allowed lower-degree potential terms, the kinetic term, and the vacuum energy. A pure monomial family need not itself have counterterm closure: a sextic interaction in three dimensions generates a quartic counterterm, and a quartic interaction generates a mass counterterm. This interpretation concerns the local divergent part or the continuum defining action, not the complete quantum effective action with its finite interaction vertices and derivative expansion.

For the massive four-dimensional quartic scalar field theory, write ${\lambda }_B$ for the bare quartic coupling, ${\lambda }_R$ for its renormalized value, and $m>0$ for the mass held fixed by a renormalization condition. Define the scalar bubble integralThe quartic one-particle-irreducible correlation function, defined as a derivative of the Euclidean quantum effective action, isThere are three bubble diagrams, one for each pairing of external momenta, and each has Feynman-diagram symmetry factor $2$. The minus sign follows equivalently from the quadratic term in the functional determinant $\frac{1}{2}\mathrm{Tr}\log [1+({\lambda }_B/2)(-{\partial }^2+m^2{)}^{-1}{\varphi }^2]$.

Radial integration at zero external momentum givesThe bare quartic vertex is not finite at fixed bare coupling. Its logarithmic ultraviolet divergence must be subtracted. Impose the momentum-subtraction scheme condition ${\Gamma }_R^{(4)}(0,0,0,0)={\lambda }_R$. To one loop this requiresFor the finite difference, a Feynman parameter and a shift of loop momentum giveThe boundary error due to shifting a sharp cutoff vanishes in this limit. Substitution yields the finite renormalized quartic scalar vertexThus finiteness requires holding the renormalized coupling fixed and allowing the bare coupling to depend on ${\Lambda }_0$. The unsubtracted assertion would be false.

The same one-loop scalar effective action explains the other effective interactions. For a constant background its expansion isThe $r=1$ tadpole diagram gives a mass correction $\delta m^2={\lambda }_R{I}_1/2$, withIt requires a mass counterterm. The field-independent vacuum energy also requires subtraction. The $r=2$ term is the quartic logarithm already treated. There is no one-loop wave-function renormalization from the quartic tadpole diagram.

For $r\ge 3$ the polygon diagrams generate finite higher-point scalar vertices, withFor example, the induced sextic term in the effective potential is $+{\lambda }_R^3{\varphi }^6/(768{\pi }^2{m}^2)$, or a six-point interaction vertex $15{\lambda }_R^3/(16{\pi }^2{m}^2)$ when normalized by $6!$. These finite terms survive the removal of the ultraviolet cutoff at fixed $m,{\lambda }_R$.

The external-momentum dependence of the bubble diagrams also generates a derivative expansion of quartic interactions. At small $P$, the subtracted scalar bubble integral is $-P^2/(96{\pi }^2{m}^2)+O(P^4/m^4)$, giving finite derivative couplings. At order ${\lambda }_R^2$ these are the new field-dependent interactions beyond the mass and quartic terms; the sextic and higher interaction vertices require higher powers of ${\lambda }_R$, though they still occur at one loop. Such coefficients are suppressed by powers of the physical mass or external momentum, not by powers of ${\Lambda }_0$. The calculation establishes an order-by-order perturbative limit; it does not establish a nonperturbative interacting continuum limit of a quantum field theory in four dimensions.

The kinetic term fixes the mass dimension of the scalar field, and the interaction then fixes that of its coupling:A relevant coupling has positive mass dimension; a marginal coupling has zero mass dimension. ThereforeFor integer dimensions the complete interacting list isIf the noninteracting $n=1,2$ terms are included, their couplings have mass dimensions $(d+2)/2$ and $2$, respectively, and are relevant for every positive $d$. This is a classification by power counting in quantum field theory; it does not assert the stability of odd monomial potentials.

Use the canonical Kähler potential and normalize superspace integration by $\int d^2\theta {\theta }^2=1$. The superspace Lagrangian density isThe full superspace integration produces a D-term, while the chiral superfield integral produces an F-term. Work in chiral coordinates for the chiral-superfield component expansion, so the omitted spacetime-derivative terms do not enter the extraction of $F$.

Set $\delta \Phi =\sqrt{2}\theta \psi +{\theta }^{2}F$. A Taylor expansion of the superpotential givesHigher terms vanish because there are only two Grassmann variables ${\theta }^{\alpha }$. With $\theta \psi ={\theta }^{\alpha }{\psi }_{\alpha }$ and $\psi \psi ={\psi }^{\alpha }{\psi }_{\alpha }$, their anticommutation gives $(\theta \psi {)}^2=-\frac{1}{2}{\theta }^2\psi \psi$. ThusThe canonical Kähler potential contributes $F^{\ast }F$. Consequently the part of the Lagrangian density containing the auxiliary field isThe fermion terms coming from the superpotential areThe Weyl spinor kinetic term and the complex scalar field kinetic term are canonical because $K={\Phi }^{†}\Phi$.

Treat $F$ and $F^{\ast }$ as independent variables in their algebraic Euler-Lagrange equations. They giveEquivalently,After eliminating the auxiliary field, the Lagrangian density contains $-V$, so the F-term scalar potential isIt is nonnegative for arbitrary complex $m,g,\phi$, since it is a modulus squared. A supersymmetric vacuum satisfies F-flatness, so its zeros are:If $m=0$ and $g\ne 0$ these coincide; if $g=0$ and $m\ne 0$ only the origin remains; if $m=g=0$ every constant $\phi$ is a supersymmetric vacuum.

To compare masses, expand around a supersymmetric vacuum $v$ satisfying $W^{\prime }(v)=0$, writing $\phi =v+\eta$. With $M=W^{\prime \prime }(v)=m+2gv$,Writing $\eta =(A+iB)/\sqrt{2}$ gives the scalar field mass term $-\frac{1}{2}|M{|}^2(A^2+B^2)$; the Weyl spinor has physical mass $|M|$. Therefore the supersymmetric mass degeneracy isIn particular, at $v=0$ the boson and fermion masses are both $|m|$. At the second supersymmetric vacuum $v=-m/g$ the sign of $M$ reverses, but the physical mass is again $|m|$. Equality here concerns fluctuations about a supersymmetric vacuum, not arbitrary backgrounds with $W^{\prime }(v)\ne 0$.

In the convention $V\supset \lambda |\phi {|}^4$ and ${\mathcal{L}}_{\mathrm{Yukawa}}=-y\phi \psi \psi +h.c.$, the same superpotential coefficient gives $\lambda =|g{|}^2$ and $y=g$. The supersymmetric relation between quartic and Yukawa couplings is thereforeIf instead the Yukawa interaction is written $-\frac{1}{2}{y}_3\phi \psi \psi +h.c.$, then $y_3=2g$ and the same physical relation is $\lambda =|y_3{|}^2/4$. The numerical factor depends on the definition of the Yukawa coupling.

The perturbative non-renormalization theorem applies to the local holomorphic superpotential in the Wilsonian effective action, using regularization in quantum field theory that preserves supersymmetry. There are two complementary ways to see why it applies here.

First, in supergraph calculations in perturbative quantum field theory, the algebra of supersymmetric covariant derivatives puts any candidate local correction to the superpotential into a full $d^4\theta$ superspace integration. It therefore corrects a D-term, such as the Kähler potential, rather than a local chiral $d^2\theta$ F-term. Turning such a contribution into an apparent chiral integral would require nonlocal inverse spacetime derivatives, schematically $1/\square$. These are excluded from the local Wilsonian effective action with a fixed nonzero momentum cutoff. This gives the diagrammatic non-renormalization theorem.

Second, the holomorphy argument for superpotential non-renormalization treats $m$ and $g$ as background chiral spurions, so the Wilsonian effective action superpotential is a holomorphic function of $\Phi ,m,g$, with no dependence on $m^{\ast }$ or $g^{\ast }$. Assign an ordinary $U(1)$ spurion charge and an R-charge as follows:The superpotential has ordinary charge zero and R-charge two, since $d^2\theta$ has R-charge minus two. A prospective holomorphic monomial $m^a{g}^b{\Phi }^c$ must therefore satisfyFor perturbative corrections, regularity at $m=0$ and $g=0$ requires nonnegative powers of the couplings. Regularity in $m$ is justified by keeping a nonzero Wilsonian momentum cutoff and retaining $\Phi$ as a light field. Hence only $c=2$ and $c=3$ survive: the existing $m{\Phi }^2$ and $g{\Phi }^3$ structures. The quadratic coefficient is fixed by the free theory at $g=0$, where there are no interaction loops. The cubic coefficient is fixed at first order in $g$, where a connected three-field interaction has only its tree-level Feynman diagram; loop corrections would require further powers of the interaction couplings, which the holomorphy argument for superpotential non-renormalization charge constraints forbid. ThusA field-independent constant is physically irrelevant in this rigid model; the spurion charges would in any case require $m^3/g^2$, which is not regular at $g=0$.

The Kähler potential can still receive wave-function renormalization. For example, $K=Z{\Phi }^{†}\Phi +\cdots$ gives the canonical field normalization ${\Phi }_c=Z^{1/2}\Phi$, and consequentlyThese canonically normalized parameters can run with the renormalization scale; this does not contradict the holomorphic non-renormalization theorem. When massless modes are included in the full quantum effective action, nonlocal terms associated with infrared divergences can also imitate an F-term. Distinguishing that action from the local Wilsonian effective action is part of the theorem's convention.

Define the stress-energy tensor by metric variation with the scalar held fixed:Use compactly supported variations to remove boundary terms. The first covariant derivative of a scalar field is an ordinary derivative, so it has no connection variation. Varying the kinetic term and the volume form gives

For the curvature term set $f={\Phi }^2$. Using part (a), its variation isApplying integration by parts twice moves each pair of derivatives from the metric variation onto $f$. The result iswhere $G_{ab}=R_{ab}-R{g}_{ab}/2$ is the Einstein tensor and $\square ={\nabla }^c{\nabla }_c$ is the d'Alembert operator. Comparing with the defining variation gives the stress-energy tensor of a nonminimally coupled scalar field:The factor $2\xi$ follows from the action's coupling $-\xi R{\Phi }^2$; a convention using $-\xi R{\Phi }^2/2$ would instead give $\xi$ in that bracket. At $\xi =0$ this reduces to the massless Klein-Gordon scalar stress-energy tensor.

For a scalar field with classical action $S[\varphi ]$, the Euclidean source convention used throughout this question givesThe generating functional produces correlation functions through functional derivatives. For example,More generally, each derivative brings down $-\varphi /\hslash$, so an $n$-point function carries $(-\hslash {)}^n/Z[0]$.

For a real free scalar field on a globally hyperbolic spacetime, takeThe Klein-Gordon equation has a well-posed initial-value problem on every Cauchy hypersurface. Its complex solution space carries the conserved Klein-Gordon inner product. Choose modes withand expandThe canonical commutation relations produce the bosonic Fock space, and $a_i|0⟩=0$ defines the chosen vacuum.

The choice of positive-frequency solutions is extra structure. In a time-dependent geometry there is no preferred time translation, so different choices mix $u_i$ with $u_i^{\ast }$ and give the vacuum ambiguity in a nonstationary spacetime. In a stationary spacetime, a timelike Killing vector field defines positive frequency by $i{\mathcal{L}}_K{u}_i={\omega }_i{u}_i$, ${\omega }_i>0$, and hence the preferred vacuum state in a stationary spacetime.

If the spacetime is stationary in the remote past and future, define normalized in- and out-modes there. Their Bogoliubov transformation isThe expected out-particle number in the in-vacuum isThus nonzero beta coefficients measure particle creation by the intervening time dependence.

For a spatially homogeneous canonical scalar field, spatial derivatives vanish. Reading the time and spatial components of its energy-momentum tensor in the comoving frame givesThus the kinetic term contributes equally to energy density and pressure, whereas the scalar potential contributes negative pressure.