Geometric and dynamical assumptions

Choose a fixed local product description of nested spatial submanifolds. Coordinates xμx^\mu comprise the shared time and smaller spatial sector; yay^a comprise the additional ordinary spatial directions. They all describe the same physical manifold. The coordinate grouping is relative to the smaller sector and gives the containing particle no inherent division among its spatial dimensions.

Assume a block metric

GAB=(gμν(x)00hab(x,y)),∣G∣=∣g∣h,h=det⁡(hab)>0. G_{AB}= \begin{pmatrix}g_{\mu\nu}(x)&0\\0&h_{ab}(x,y)\end{pmatrix}, \qquad \sqrt{|G|}=\sqrt{|g|}\sqrt h,\qquad h=\det(h_{ab})>0.

The integration domain FF in yy is fixed. The base metric is independent of yy. These hypotheses specify the projection being calculated; a general embedding, moving integration domain, or off-diagonal metric requires its corresponding geometric terms.

Let TAB=TBAT^{AB}=T^{BA} be the full source satisfying ∇A(G)TAB=0\nabla_A^{(G)}T^A{}_B=0. One sufficient origin is a diffeomorphism-invariant matter action with all its participating fields on shell. In the convention

δSm=12∫∣G∣ TABδGAB dx dy, \delta S_{\rm m}=\frac12\int\sqrt{|G|}\,T^{AB}\delta G_{AB}\,dx\,dy,

the metric change under a compactly supported vector field ξ\xi is δGAB=2∇(AξB)\delta G_{AB}=2\nabla_{(A}\xi_{B)}. On-shell invariance and integration by parts give 0=−∫∣G∣(∇ATAB)ξB0=-\int\sqrt{|G|}(\nabla_AT^{AB})\xi_B, hence conservation. This is the standard action identity derived in Carroll's GR notes, §5, equations 5.35–5.38. All relevant interaction contributions to TT must be retained.

Exact integrated balance

Define the tangential stress integrated over the full spatial fibre:

Tˉμν(x)=∫Fh Tμν(x,y) dpy. \bar T^\mu{}_\nu(x)=\int_F\sqrt h\,T^\mu{}_\nu(x,y)\,d^py.

Expanding covariant conservation for a symmetric tensor gives the coordinate identity

∂A ⁣(∣G∣TAν)=12∣G∣TAB∂νGAB. \partial_A\!\left(\sqrt{|G|}T^A{}_\nu\right) =\frac12\sqrt{|G|}T^{AB}\partial_\nu G_{AB}.

To see the factor 1/21/2, the connection contribution with the lower index is ΓAνCTAC\Gamma^C_{A\nu}T^A{}_C. After lowering the contracted index, symmetry of TABT^{AB} cancels the two oppositely signed metric derivatives in the Christoffel symbol, leaving TAB∂νGAB/2T^{AB}\partial_\nu G_{AB}/2. The upper-index connection trace supplies the volume derivative on the left.

Integrate the identity over FF and divide by ∣g∣\sqrt{|g|}. The term containing ∂νgαβ\partial_\nu g_{\alpha\beta} combines with the base volume derivative into ∇μ(g)Tˉμν\nabla^{(g)}_\mu\bar T^\mu{}_\nu. Consequently,

∇μ(g)Tˉμν=12∫Fh Tab∂νhab dpy−Fν,Fν=∫∂FnaTaν dSh. \boxed{ \nabla^{(g)}_\mu\bar T^\mu{}_\nu =\frac12\int_F\sqrt h\,T^{ab}\partial_\nu h_{ab}\,d^py -\mathcal F_\nu, \qquad \mathcal F_\nu=\int_{\partial F}n_aT^a{}_\nu\,dS_h. }

The normal and surface measure are those of the fibre metric at fixed xx. For a noncompact fibre, Fν\mathcal F_\nu means the limiting flux through an exhaustion by bounded domains. Require convergence of the displayed integrals and permission to differentiate them, for example by locally uniform integrable bounds on the integrands and their base derivatives. The same derivation can be formulated weakly against compactly supported base test fields.

If hh is independent of xx and the boundary flux vanishes, then ∇μ(g)Tˉμν=0\nabla^{(g)}_\mu\bar T^\mu{}_\nu=0. Finite polynomial combinations of oscillator eigenmodes give useful decaying test sources when the remaining coefficients are controlled. The actual stress tensor determines which decay estimates are needed.

A varying fibre metric

The volume factor belongs inside Tˉ\bar T. The metric term on the right measures exchange with a spatial metric that varies along the base. For a concrete check, take

ds2=−dt2+dx2+e2axdy2,−L≤y≤L,TAB=PδBA, ds^2=-dt^2+dx^2+e^{2ax}dy^2,\qquad -L\leq y\leq L, \qquad T^A{}_B=P\delta^A_B,

with constants a,L>0a,L>0 and PP. This tensor is conserved by metric compatibility. Its tangential boundary flux is zero, while

Tˉxx=2LPeax,∂xTˉxx=2LaPeax. \bar T^x{}_x=2LP e^{ax},\qquad \partial_x\bar T^x{}_x=2LaP e^{ax}.

Since Tyy=Pe−2axT^{yy}=Pe^{-2ax} and ∂xhyy=2ae2ax\partial_xh_{yy}=2ae^{2ax}, the metric term is also 2LaPeax2LaPe^{ax}. Thus a conserved full tensor can have a nonzero integrated tangential divergence supplied exactly by the varying spatial geometry. This example checks the balance law without requiring an infinite fibre volume or a particle assignment.

Sampling a finite portion

For a smooth scalar weight w(x,y)w(x,y), define Tˉwμν=∫Fh wTμν dy\bar T_w^\mu{}_\nu=\int_F\sqrt h\,w T^\mu{}_\nu\,dy. Multiplying the coordinate identity by ww and applying the product rule gives

∇μ(g)Tˉwμν=12∫Fh wTab∂νhab dy−∫∂FwnaTaν dSh+∫Fh[(∂μw)Tμν+(∂aw)Taν]dy. \begin{aligned} \nabla^{(g)}_\mu\bar T_w^\mu{}_\nu ={}&\frac12\int_F\sqrt h\,wT^{ab}\partial_\nu h_{ab}\,dy -\int_{\partial F}w n_aT^a{}_\nu\,dS_h\\ &+\int_F\sqrt h\left[ (\partial_\mu w)T^\mu{}_\nu+(\partial_a w)T^a{}_\nu \right]dy. \end{aligned}

Even a time-independent aperture can exchange tangential momentum through its spatial gradient. For an explicit flat example use the symmetric spatial tensor

Txx=f(x)q′′(y),Tyx=Txy=−f′(x)q′(y),Tyy=f′′(x)q(y). T^{xx}=f(x)q''(y),\qquad T^{yx}=T^{xy}=-f'(x)q'(y),\qquad T^{yy}=f''(x)q(y).

Both divergences vanish by commuting derivatives. With q=e−y2q=e^{-y^2} and w=e−βy2w=e^{-\beta y^2}, β>0\beta>0, the full boundary flux is zero but

∂xTˉwxx=−2βπ(1+β)3/2f′(x)=∫Rw′(y)Tyx(x,y) dy. \partial_x\bar T_w^x{}_x =-\frac{2\beta\sqrt\pi}{(1+\beta)^{3/2}}f'(x) =\int_{\mathbb R}w'(y)T^{yx}(x,y)\,dy.

The sampling term is therefore essential even for a smooth, rapidly decaying source.

Factorized normalization

For an xx-independent fibre metric and the specified factorization Tμν=∣χ(y)∣2tμν(x)T^\mu{}_\nu=|\chi(y)|^2 t^\mu{}_\nu(x) with ∫h∣χ∣2dy=1\int\sqrt h|\chi|^2dy=1, direct integration gives Tˉμν=tμν\bar T^\mu{}_\nu=t^\mu{}_\nu. This is a normalization identity. Kinetic and interaction terms must actually permit that factorization. For a finite coherent mode expansion, a bilinear expression integrates with the Gram matrix ∫h χi∗χjdy\int\sqrt h\,\chi_i^*\chi_jdy; orthonormality removes those cross terms only when the operator acting on the modes permits it.

These formulas provide conservation and normalization tests for an effective source. A metric field equation and its coupling coefficient still come from the specified gravitational action.

spatial_profile_checks.py checks the Christoffel identity, the varying-metric example, finite-boundary flux, weighted sampling, and mode normalization.