% Save as 085T_Landmark_VIID.tex
% Compile with pdflatex twice. No external bibliography is required.
%
% Source declaration: formulas from 085T_Landmark_VIIC_v1.tex.
% Polish pass: Blocks 1--4 OK. Math frozen. Closed answer: Yes.
% Individual audit marks appear beside the four checked calculations.
% Build status: compilation was not executed in this environment.

\documentclass[11pt,a4paper]{article}

\usepackage[T1]{fontenc}
\usepackage[utf8]{inputenc}
\usepackage{lmodern}
\usepackage{microtype}
\usepackage[margin=25mm]{geometry}
\usepackage{amsmath,amssymb,amsthm,mathtools}
\usepackage{booktabs,array,tabularx,enumitem}
\usepackage[hidelinks]{hyperref}

\setlength{\emergencystretch}{3em}
\setlist[itemize]{leftmargin=*,itemsep=3pt,topsep=4pt}
\setlist[enumerate]{leftmargin=*,itemsep=3pt,topsep=4pt}
\numberwithin{equation}{section}

\newtheorem{theorem}{Theorem}[section]
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}

\newcommand{\R}{\mathbb R}
\newcommand{\C}{\mathbb C}
\newcommand{\D}{\mathcal D}
\newcommand{\Ezero}{\mathsf E_0}
\newcommand{\U}{\mathsf U}
\newcommand{\eps}{\varepsilon}
\newcommand{\dd}{\,\mathrm d}
\newcommand{\supp}{\operatorname{supp}}
\newcommand{\ev}{\operatorname{ev}_0}
\newcommand{\Rel}{\mathrm{Rel}}
\newcommand{\la}{\langle}
\newcommand{\ra}{\rangle}
\newcommand{\cC}{\mathfrak C}
\newcommand{\Ps}{\mathcal H}
\newcommand{\Om}{\mathrm d\Omega}
\newcommand{\Os}{\operatorname{Os}}

\title{\textbf{Finite part of the regulated Cauchy-reflection triangle}\\[0.4em]
\large A strong numerical distributional limit for the VII-C family}
\author{}
\date{}

\usepackage{tikz}
\begin{document}
\maketitle

\begin{abstract}
For the insertion-kernel regularization fixed in VII-C, we prove
existence of the strong distributional finite part of its numerical
reflection-triangle coefficient. The proof applies to arbitrary
compact joint tests. An explicit normalization ledger identifies the
leading reflected contribution with the prescribed light-ray residue,
including its density and every Fourier factor. A joint spatial
Fourier argument gives a regulator-independent ultraviolet
decomposition: the radial coefficient differs from that residue by
\(O(r^{-1})\) in a fixed finite-order distribution seminorm, and the
remaining regulator-placement error is
\(O(\eps(1+|\log\eps|))\).
The estimates retain the full phase and are uniform through the
intersection of reflected collinearity with the first--third
observation diagonal. For two admissible normalized radial
mollifiers, the complete conversion is \(c_{ab}r\), with an explicit
real coefficient, as a corollary of the ultraviolet decomposition.
The result is a numerical subtraction prescription, not an algebra
construction.
\end{abstract}

\begin{figure}[ht]
\centering
\begin{tikzpicture}[x=1cm,y=1cm]
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (a) at (0,0) {Regulated triangle};
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (b) at (5,0) {Finite part at fixed scale};
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (c) at (10,0) {Regulator removed};
\node[draw,dashed,rounded corners,align=center,text width=9.6cm,font=\small,inner sep=3pt] (s) at (5,-2.35) {The finite part is not the limit with the regulator removed.};
\draw[->] (a) -- (b);
\draw[->] (b) -- (c);
\draw[->] (b) -- (s);
\end{tikzpicture}
\caption{A finite part is not the removed regulator. The three boxes name objects this manuscript keeps separate. The dashed box is the stated stop. Not a proof, and not an admission of a quantum coefficient.}
\label{fig:scope-VII-D}
\end{figure}

\begin{figure}[ht]
\centering
\begin{tikzpicture}[x=1cm,y=1cm]
\node[draw,rounded corners,align=center,text width=4.3cm,font=\small,inner sep=3pt] (a) at (0,0) {Regulated triangle};
\node[draw,rounded corners,align=center,text width=4.3cm,font=\small,inner sep=3pt] (b) at (6.4,0) {Regulator removed};
\node[draw,dashed,rounded corners,align=center,text width=9.6cm,font=\small,inner sep=3pt] (s) at (3.2,-2.35) {The finite part is not the limit with the regulator removed.};
\draw[->] (a) -- (s);
\draw[->] (b) -- (s);
\end{tikzpicture}
\caption{A finite part is not the removed regulator. The boxes name objects this manuscript keeps separate. The dashed box is the stated stop. Not a proof, and not an admission of a quantum coefficient.}
\label{fig:split-VII-D}
\end{figure}

\begin{figure}[ht]
\centering
\begin{tikzpicture}[x=1cm,y=1cm]
\node[draw,rounded corners,align=center,text width=7.2cm,font=\small,inner sep=3pt] (a) at (0,0) {Regulated triangle};
\node[draw,rounded corners,align=center,text width=7.2cm,font=\small,inner sep=3pt] (b) at (0,-1.7) {Finite part at fixed scale};
\node[draw,rounded corners,align=center,text width=7.2cm,font=\small,inner sep=3pt] (c) at (0,-3.4) {Regulator removed};
\node[draw,dashed,rounded corners,align=center,text width=7.2cm,font=\small,inner sep=3pt] (s) at (0,-5.2) {The finite part is not the limit with the regulator removed.};
\draw[->] (a) -- (b);
\draw[->] (b) -- (c);
\draw[->] (c) -- (s);
\end{tikzpicture}
\caption{Reading downward is not a transfer. The boxes name objects this manuscript keeps separate. The dashed box is the stated stop. Not a proof, and not an admission of a quantum coefficient.}
\label{fig:chain-VII-D}
\end{figure}

\begin{figure}[ht]
\centering
\begin{tikzpicture}[x=1cm,y=1cm]
\node[draw,rounded corners,align=center,text width=4.2cm,font=\small,inner sep=3pt] (a) at (0,0) {Regulated triangle};
\node[align=center,font=\small] (m) at (5.2,0) {not the same object};
\node[draw,rounded corners,align=center,text width=4.2cm,font=\small,inner sep=3pt] (b) at (10.4,0) {Regulator removed};
\node[draw,dashed,rounded corners,align=center,text width=9.6cm,font=\small,inner sep=3pt] (s) at (5.2,-2.35) {Finite part at fixed scale. The finite part is not the limit with the regulator removed.};
\draw (a) -- (m);
\draw (m) -- (b);
\draw[->] (m) -- (s);
\end{tikzpicture}
\caption{Regulated triangle is not Regulator removed. The boxes name objects this manuscript keeps separate. The dashed box is the stated stop. Not a proof, and not an admission of a quantum coefficient.}
\label{fig:pair-VII-D}
\end{figure}

\begin{figure}[ht]
\centering
\begin{tikzpicture}[x=1cm,y=1cm]
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (a) at (0,0) {Named: Regulated triangle};
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (b) at (5,0) {Used here: Finite part at fixed scale};
\node[draw,rounded corners,align=center,text width=3.15cm,font=\small,inner sep=3pt] (c) at (10,0) {Not claimed: Regulator removed};
\node[draw,dashed,rounded corners,align=center,text width=9.6cm,font=\small,inner sep=3pt] (s) at (5,-2.5) {The finite part is not the limit with the regulator removed.};
\draw[->] (a) -- (b);
\draw[->] (b) -- (c);
\draw[->] (b) -- (s);
\end{tikzpicture}
\caption{Three statuses stay separate. The boxes name objects this manuscript keeps separate. The dashed box is the stated stop. Not a proof, and not an admission of a quantum coefficient.}
\label{fig:gate-VII-D}
\end{figure}

Figure~\ref{fig:split-VII-D} records one view of the same separation. Figure~\ref{fig:chain-VII-D} records one view of the same separation. Figure~\ref{fig:pair-VII-D} records one view of the same separation. Figure~\ref{fig:gate-VII-D} records one view of the same separation. The drawings add no hypothesis.


Figure~\ref{fig:scope-VII-D} records that separation.
The drawing adds no hypothesis.



\section*{Scope and fixed data}

Let \(M=\R^{1+3}\), \(M_+=\{t>0\}\), and \(X=M_+^3\).
We ask whether
\[
 u_\mu^\rho
 =
 \lim_{\eps\downarrow0}
 \left(
 T_\eps^\rho
 -
 \log\!\left(\frac1{\mu\eps}\right)r
 \right),
 \qquad \mu>0,
\]
exists strongly in \(\D'(X)\), with the counterterm class fixed to
\[
                              \cC=\R r.
\]
The family \(T_\eps^\rho\) is precisely the numerical VII-C family:
both slots of each insertion are regularized with the same fixed
mollifier \(\rho\), while the single middle observable is unsmeared.

The formulas on the following convention sheet are from
\texttt{085T\_Landmark\_VIIC\_v1.tex}~\cite{VIIC}.
These displayed formulas define the family studied in this paper.
The proof requires no companion methods manuscript and no
unspecified normalization.

We reserve \(\mu\) for the subtraction scale and write
\(\vartheta\) for the mollifier profile denoted by \(\mu\) in VII-C.
The symbol \(g\) denotes the middle test; \(s=y^0\) is a time
coordinate.

\clearpage

\section*{Source convention sheet}
\begingroup
\footnotesize
\setlength{\abovedisplayskip}{4pt}
\setlength{\belowdisplayskip}{4pt}
\setlength{\abovedisplayshortskip}{3pt}
\setlength{\belowdisplayshortskip}{3pt}

\paragraph{Insertion and placement
(\texttt{eq:B}, \texttt{eq:C}, \texttt{eq:Beps} in VII-C).}
\begin{equation}
\begin{gathered}
 \eta=\operatorname{diag}(-1,1,1,1),\qquad
 \square=-\partial_t^2+\Delta,\qquad
 (Ch)^\mu=\partial_\nu h^{\mu\nu}
       -\tfrac12\partial^\mu(\eta^{\alpha\beta}h_{\alpha\beta}),
 \qquad L=\Ezero C,\\
 \square\Ezero j=j,\qquad
 (\Ezero j)|_{t=0}=0,\qquad
 \partial_t(\Ezero j)|_{t=0}=0,\\
 B_f=-\int f(x)(Lh)^\mu(x)\partial_\mu\phi(x)\dd^4x
       =\la b_f,\phi\otimes h\ra,\qquad W_g=\int g h_{00}^2,\\
 \rho\in C_c^\infty(\R^4),\quad
 \rho\geq0,\quad \rho\text{ real, even and Euclidean-radial},
 \quad \int\rho=1,\\
 S_\eps=\rho_\eps*,\qquad
 \rho_\eps(x)=\eps^{-4}\rho(x/\eps),\qquad
 B_{f,\eps}=B_f[S_\eps h,S_\eps\phi],\qquad
 b_{f,\eps}=(\rho_\eps\otimes\rho_\eps)*b_f.
\end{gathered}
\tag{C.1}\label{sc:placement}
\end{equation}

\paragraph{Even covariance
(\texttt{eq:Hhh}, \texttt{eq:Hscalar}, \texttt{eq:positive}).}
\begin{equation}
\begin{gathered}
 H_{hh\,\alpha\beta\rho\sigma}
 =P_{\alpha\beta\rho\sigma}\Delta_0^+,\qquad
 P_{\alpha\beta\rho\sigma}
 =\eta_{\alpha(\rho}\eta_{\sigma)\beta}
       -\tfrac12\eta_{\alpha\beta}\eta_{\rho\sigma},\\
 H_{\phi\phi}=\Delta_m^+,\qquad H_{h\phi}=H_{\phi h}=0,\qquad
 \Delta_m^+(x-y)=
 \int\frac{\dd^3q}{(2\pi)^3\,2\omega_m(q)}
 e^{-i\omega_m(q)(t-s)+iq\cdot(\mathbf x-\mathbf y)},\\
 \omega_m(q)=\sqrt{|q|^2+m^2},\qquad
 C_x^\mu H_{hh\,\rho\sigma}
 =\delta^\mu_{(\rho}\partial_{\sigma)}^x\Delta_0^+,\\
 C_x^\mu C_y^\nu H_{hh}
 =-\tfrac12\eta^{\mu\nu}\square_x\Delta_0^+=0,
 \qquad H_{ww}=0.
\end{gathered}
\tag{C.2}\label{sc:covariance}
\end{equation}

\paragraph{Zero-Cauchy mode and mixed covariance
(\texttt{eq:K}).}
\begin{equation}
\begin{gathered}
 \Ezero(e^{-iPt})
 =-\frac{it}{2P}e^{-iPt}+\frac{i\sin(Pt)}{2P^2},
 \qquad P=|p|,\\
 K(x,y)=
 \int\frac{\dd^3p}{(2\pi)^3}e^{ip\cdot(\mathbf x-\mathbf y)}
 \left[
 \left(-\frac{t}{4P}+\frac{i}{8P^2}\right)e^{-iP(t-s)}
 -\frac{i}{8P^2}e^{iP(t+s)}
 \right].
\end{gathered}
\tag{C.3}\label{sc:K}
\end{equation}
The factors \(P^{-1}\) and \(P^{-2}\) are locally integrable in
three momentum dimensions. No graviton positivity is assumed.

\paragraph{Profiles and spectral identity
(\texttt{eq:F}, \texttt{eq:spectral}, \texttt{lem:spectral}).}
\begin{equation}
\begin{gathered}
 \vartheta(u)=\widehat\rho(-u,u\mathbf n),\qquad
 M_\eps(r)=
 \vartheta\!\left(\eps\sqrt{r^2+m^2/2}\right),\\
 K_\eps=K\text{ with the multiplier }\vartheta(\eps P),\qquad
 F_{f,\eps}(y,q)=
 \int f(x)K_\eps(x,y)e^{-i\omega_m(q)t+iq\cdot\mathbf x}\dd^4x,\\
 \ev[\hbar^3](B_{f,\eps}\star_H W_g\star_H B_{f',\eps})
 =\frac1{(2\pi)^3}\int g(y)\dd^4y\int
 \omega_m(q)M_\eps(|q|)^2
 F_{f,\eps}(y,q)\overline{F_{f',\eps}(y,q)}\dd^3q
 \quad(f,f',g\text{ real}).
\end{gathered}
\tag{C.4}\label{sc:spectral}
\end{equation}
There is one graph type. Its factor \(2\) counts the assignments of
the two middle legs. The middle observable is not regularized.

\paragraph{Locked residue
(\texttt{eq:ray}, \texttt{eq:residue}, \texttt{thm:log}).}
\begin{equation}
\begin{gathered}
 A_f(y,\mathbf n)
 =\int_\R f(t,\mathbf y+(t+s)\mathbf n)\dd t,\qquad
 c_0=\frac1{64(2\pi)^3},\\
 \mathcal C(f,g,f')
 =c_0\int g(y)\dd^4y\int_{S^2}
 A_f(y,\mathbf n)A_{f'}(y,\mathbf n)\Om(\mathbf n),\\
 \mathcal G_\eps(f,g,f')
 =\mathcal C(f,g,f')\log(\eps^{-1})+O(1).
\end{gathered}
\tag{C.5}\label{sc:residue}
\end{equation}
Here \(\Om\) is ordinary surface area measure:
\(\int_{S^2}\Om=4\pi\).

\paragraph{Off-shell identity
(\texttt{eq:reconstruction}).}
\begin{equation}
 h=\U(\gamma_0h,\gamma_1h)+\Ezero j,\qquad j=\square h.
 \tag{C.6}\label{sc:off}
\end{equation}
The source is retained. Equation \eqref{sc:covariance} does not set
a classical field to zero.

\endgroup
\clearpage

\section{The exact distribution and intersection coordinates}
\label{sec:exact}

Write
\[
 x=(t,\mathbf x),\qquad y=(s,\mathbf y),\qquad
 x'=(t',\mathbf x').
\]
The scalar mass \(m\geq0\) and each admissible mollifier are fixed
independently of \(\eps\).

Our Fourier convention is
\[
 \widehat v(k)=\int e^{-ik\cdot z}v(z)\dd z,
 \qquad
 v(z)=\int\frac{\dd k}{(2\pi)^{\dim z}}
                 e^{ik\cdot z}\widehat v(k).
\]
The dot product in this convention is the coordinate pairing, not
the Lorentzian metric pairing. Radiality of \(\rho\) is Euclidean.

For a compact \(K\subset X\), let
\[
 p_{K,N}(\Phi)
 =\max_{|\gamma|\leq N}\sup_X|\partial^\gamma\Phi|,
 \qquad \Phi\in\D_K(X).
\]
Tests are extended by zero to the ambient Euclidean space when
taking Fourier transforms. Their supports are separated from the
boundary of \(X\).

All estimates below can be taken with the single test order \(64\).
Its size is inessential; fixing it makes independence from the
cutoff and the radial parameter explicit.

Strong convergence in \(\D'(X)\) means uniform convergence on
bounded subsets of \(\D(X)\). A bounded subset of \(\D(X)\) is
contained and bounded in some fixed-support test space. Thus
\begin{equation}
 |\la U_\eps-U,\Phi\ra|
 \leq C_K a(\eps)p_{K,64}(\Phi),
 \qquad a(\eps)\longrightarrow0,
 \label{eq:strongtemplate}
\end{equation}
on every compact implies strong convergence.

\subsection{Positive-cutoff definition and regulator placement}

Define
\begin{equation}
 D_\eps(z)
 =
 \int_{\R^3}\frac{\dd^3q}{(2\pi)^3}
 \frac{\omega_m(q)}2 M_\eps(|q|)^2
 e^{-i\omega_m(q)z^0+iq\cdot\mathbf z}.
 \label{eq:Deps}
\end{equation}

\begin{proposition}[Exact positive-cutoff distribution and placement]
\label{prop:positive}
The field regularization and the mixed-kernel convolution in
\eqref{sc:placement} define the same insertion. Their numerical
coefficient is
\begin{equation}
 T_\eps^\rho(x,y,x')
 =
 2K_\eps(x,y)\overline{K_\eps(x',y)}D_\eps(x-x').
 \label{eq:Texact}
\end{equation}
For every \(\eps>0\), this is a smooth function on \(X\), and
\[
 |\la T_\eps^\rho,\Phi\ra|
 \leq C_{K,\eps,\rho,m}p_{K,0}(\Phi),
 \qquad \Phi\in\D_K(X).
\]
Its pairing on real product tests is exactly
\eqref{sc:spectral}.
\end{proposition}

\begin{proof}
\paragraph{Field and kernel convolution.}
For \(\check\rho(z)=\rho(-z)\), transposition of convolution gives
\[
 \la b_f,S_\eps\phi\otimes S_\eps h\ra
 =
 \la(\check\rho_\eps\otimes\check\rho_\eps)*b_f,
          \phi\otimes h\ra.
\]
Since \(\rho\) is even, this is the kernel convolution displayed
in \eqref{sc:placement}. The identity holds on arbitrary smooth
configurations: the original mixed kernel has compact support,
and the convolution kernels are compactly supported.

\paragraph{Plane-wave multiplier identities.}
The coefficients of \(C\) are constant, so
\[
                         C S_\eps h=S_\eps C h
\]
for every smooth \(h\). For any constant tensor \(v\),
\[
 S_\eps\bigl(v e^{-iPz^0+ip\cdot\mathbf z}\bigr)
 =
 \widehat\rho\bigl(\eps(-P,p)\bigr)
        v e^{-iPz^0+ip\cdot\mathbf z}
 =
 \vartheta(\eps P)
        v e^{-iPz^0+ip\cdot\mathbf z}.
\]
This applies in particular to the de Donder polarizations in
\eqref{sc:covariance}. Applying the fixed operator \(\Ezero C\)
to this equality retains the scalar multiplier
\(\vartheta(\eps P)\). This proves the displayed formula for
\(K_\eps\), without commuting \(S_\eps\) through \(\Ezero\).

For the scalar leg, Euclidean radiality gives
\[
 \widehat\rho\bigl(\eps(-\omega_m(q),q)\bigr)
 =
 \vartheta\!\left(\eps\sqrt{|q|^2+m^2/2}\right)
 =M_\eps(|q|).
\]
The two scalar insertion slots therefore supply \(M_\eps^2\).
The derivative factors are
\[
 (-i\omega_m(q))(+i\omega_m(q))=\omega_m(q)^2;
\]
multiplying the scalar spectral density \(1/(2\omega_m(q))\)
gives \eqref{eq:Deps}. The two insertion signs also multiply to
\(+1\). The middle observable is untouched.

\paragraph{Smoothness and the product-test identity.}
For fixed positive \(\eps\), every coordinate derivative of the
momentum integrands in \(K_\eps\) and \(D_\eps\) is integrable,
locally uniformly in the observations. At zero momentum the only
graviton denominators are \(P^{-1}\) and \(P^{-2}\), locally
integrable in three dimensions. The scalar density is locally
integrable also for \(m=0\). At infinity the profiles provide
arbitrary polynomial decay.

Thus \eqref{eq:Texact} is a product of smooth functions. Its
compact-test bound follows by a supremum on a compact
neighborhood of \(K\). Integrating against \(f\otimes g\otimes f'\)
and using \eqref{eq:Deps} yields \eqref{sc:spectral}.
\end{proof}

\paragraph{Complex tests.}
The distribution \eqref{eq:Texact} is complex-linear in its joint
test. On complex product tests, the spectral factors are
\[
 F_{f,\eps}(y,q)\,
 \overline{F_{\overline{f'},\eps}(y,q)}.
\]
The formula with \(\overline{F_{f',\eps}}\) alone applies to real
\(f'\). Kernel conjugation in \eqref{eq:Texact} does not conjugate
the test.

\subsection{Intersection coordinates}

On \(\mathbf x\neq\mathbf y\), define
\begin{equation}
\begin{gathered}
 a=t+s,\qquad b=t'+s,\qquad
 \mathbf v=\mathbf x-\mathbf y,\qquad
 \mathbf n=\frac{\mathbf v}{|\mathbf v|},\qquad
 \sigma=|\mathbf v|-a,\\
 \mathbf w=\mathbf x'-\mathbf y-\frac ba(\mathbf x-\mathbf y),
 \qquad \mathbf z=\mathbf y+s\mathbf n,\qquad d=t-t'.
\end{gathered}
\label{eq:chart}
\end{equation}
The chart domain has \(a>0\) and \(a+\sigma>0\).
The inverse spatial formulas are
\[
 \mathbf y=\mathbf z-s\mathbf n,\qquad
 \mathbf x=\mathbf y+(a+\sigma)\mathbf n,\qquad
 \mathbf x'=\mathbf y+\frac ba(a+\sigma)\mathbf n+\mathbf w.
\]
Consequently,
\begin{equation}
 \Rel=\{\sigma=0,\mathbf w=0\},\qquad
 \Delta_{13}=\{d=0,\mathbf w=0\},\qquad
 \Rel\cap\Delta_{13}=\{\sigma=0,\mathbf w=0,d=0\}.
 \label{eq:chartrelations}
\end{equation}

For fixed times, the linear transformation to
\((\mathbf y,\mathbf v,\mathbf w)\) has determinant of absolute
value one. Replacing \(\mathbf y\) by
\(\mathbf z=\mathbf y+s\mathbf n\) is triangular with unit
determinant. Polar coordinates in \(\mathbf v\) therefore give
\begin{equation}
 \dd X
 =(a+\sigma)^2
 \dd^3\mathbf z\,\Om(\mathbf n)\,
 \dd t\,\dd s\,\dd t'\,\dd\sigma\,\dd^3\mathbf w .
 \label{eq:density}
\end{equation}
The full transformation, including the times, has the same
Jacobian because its time block is the identity.
Replacing \(t'\) by \(d=t-t'\) has absolute Jacobian one.
There is no geometric \(1/d\).

\section{The locked residue}
\label{sec:residue}

Let
\[
 Z=\R^3\times S^2\times(0,\infty)^3
\]
with coordinates \((\mathbf z,\mathbf n,t,s,t')\), and set
\begin{align}
 \Psi(\mathbf z,\mathbf n,t,s,t')
 &=
 \bigl((t,\mathbf z+t\mathbf n),
       (s,\mathbf z-s\mathbf n),
       (t',\mathbf z+t'\mathbf n)\bigr),
 \label{eq:Psi}\\
 \dd\nu
 &=\dd^3\mathbf z\,\Om(\mathbf n)\,\dd t\,\dd s\,\dd t'.
 \label{eq:nu}
\end{align}

\begin{definition}[Prescribed VII-C residue]
\label{def:residue}
The subtraction distribution is
\begin{equation}
 r=c_0\Psi_*\nu,\qquad c_0=\frac1{64(2\pi)^3}.
 \label{eq:rdef}
\end{equation}
Equivalently, for an arbitrary joint test,
\begin{equation}
\begin{split}
 \la r,\Phi\ra
 =c_0\int &\dd s\,\dd^3\mathbf y\,\dd t\,\dd t'
 \int_{S^2}\\
 &\Phi\bigl(t,\mathbf y+(t+s)\mathbf n;\
           s,\mathbf y;\
           t',\mathbf y+(t'+s)\mathbf n\bigr)
 \Om(\mathbf n).
\end{split}
\label{eq:rjoint}
\end{equation}
\end{definition}

The map \(\Psi\) is a proper embedding over compact subsets of
\(X\). Its inverse on its image is given by the observation times
and
\[
 \mathbf n=\frac{\mathbf x-\mathbf y}{t+s},
 \qquad \mathbf z=\mathbf y+s\mathbf n.
\]
On the preimage of a compact in \(X\), the time coordinates remain
bounded and positive, the spatial coordinate \(\mathbf z\) is
bounded, and \(\mathbf n\in S^2\). The preimage is closed in this
bounded parameter region and hence compact.

It follows directly that \(r\) is a positive order-zero
distribution, with support exactly \(\Rel=\Psi(Z)\).
Its product-test pairing is \eqref{sc:residue}.
The separated-support nonzero witness is already established in
VII-C, \texttt{lem:separated}.

Relative to the ambient density \eqref{eq:density},
\begin{equation}
                  r=c_0a^{-2}\delta(\sigma)\delta^{(3)}(\mathbf w).
 \label{eq:rlocal}
\end{equation}
Indeed, multiplying this expression by \((a+\sigma)^2\) and
integrating in \(\sigma,\mathbf w\) gives exactly
\eqref{eq:rjoint}. There is no \(\delta(d)\) in this measure.

If the same embedding and parameter density are written as
\(r_\Sigma=c_\Sigma\Psi_*\nu\), their normalization comparison is
\[
                      r_\Sigma=64(2\pi)^3c_\Sigma\,r.
\]
This paper always subtracts \eqref{eq:rdef}.
The fact that the actual leading graph coefficient equals this
prescribed distribution is proved in
Proposition~\ref{prop:ledger}; it is not built into the definition
of the radial family.

\section{Joint Fourier estimates and convergent sectors}
\label{sec:Fourier}

\subsection{The joint spatial transform}

For fixed \((t,s,t',\mathbf y)\), define
\begin{equation}
\begin{split}
 \Ps_\Phi(t,s,t',\mathbf y;\alpha,\beta)
 =\int_{\R^6}
 &\Phi(t,\mathbf x;s,\mathbf y;t',\mathbf x')\\
 &{}\times
 e^{i\alpha\cdot(\mathbf x-\mathbf y)
    -i\beta\cdot(\mathbf x'-\mathbf y)}
 \dd^3\mathbf x\,\dd^3\mathbf x'.
\end{split}
\label{eq:jointFT}
\end{equation}
This preserves all correlations of the joint test.

We shall use the fixed momentum-decay exponent \(N=12\).
Applying
\((1-\Delta_{\mathbf x})^6(1-\Delta_{\mathbf x'})^6\)
under the integral gives
\begin{equation}
 |\partial_t^j\partial_{t'}^{j'}\Ps_\Phi|
 \leq C_Kp_{K,64}(\Phi)
 (1+|\alpha|)^{-12}(1+|\beta|)^{-12},
 \qquad j+j'\leq8.
 \label{eq:FTbound}
\end{equation}
Only \(24+j+j'\leq32\) derivatives are needed for this estimate;
order \(64\) also accommodates the fixed localizations used below.
The variables \((t,s,t',\mathbf y)\) range over a fixed compact
projection of \(K\).

\begin{lemma}[Shifted momentum bounds]
\label{lem:convolution}
Let \(r\geq1\), \(\boldsymbol\nu\in S^2\),
\(\ell\in\{1,2\}\), \(k\in\{0,1,2\}\), and \(N>k+3\).
Then
\begin{align}
 \int_{\R^3}
 \frac{(1+|\alpha|)^{k-N}}
      {|\alpha-r\boldsymbol\nu|^\ell}\dd^3\alpha
 &\leq C r^{-\ell},
 \label{eq:shiftbound}\\
 \int_{|\alpha|>r/2}
 \frac{(1+|\alpha|)^{k-N}}
      {|\alpha-r\boldsymbol\nu|^\ell}\dd^3\alpha
 &\leq C r^{3-\ell+k-N}.
 \label{eq:shifttail}
\end{align}
The constants are uniform in \(r,\boldsymbol\nu\).
For \(0\leq r\leq1\), the first integral is uniformly bounded.
In particular, all these assertions apply with \(N=12\).
\end{lemma}

\begin{proof}
On \(|\alpha|\leq r/2\), the denominator is at least
\((r/2)^\ell\), and the remaining weight is integrable.

On \(|\alpha|>r/2\), first consider
\(|\alpha-r\boldsymbol\nu|\leq r/2\).
Here \(|\alpha|\) is comparable to \(r\), and
\[
 \int_{|p|\leq r/2}|p|^{-\ell}\dd^3p
 =C_\ell r^{3-\ell}.
\]
This part is bounded by \(Cr^{3-\ell+k-N}\).
In the remaining part the denominator gives \(Cr^{-\ell}\),
while
\[
 \int_{|\alpha|>r/2}(1+|\alpha|)^{k-N}\dd^3\alpha
 \leq Cr^{3+k-N}.
\]
This proves the tail bound. Since \(N>k+3\), the tail also
satisfies the full-integral bound.

For \(0\leq r\leq1\), split into a fixed-size neighborhood of
\(\alpha=r\boldsymbol\nu\) and its complement.
Local integrability of \(|p|^{-\ell}\) treats the first region,
and the fixed polynomial weight treats the second.
\end{proof}

\paragraph{Certificate.}
These bounds include all momentum ratios and all directions.
They use no time separation, regulator decay, or angular
stationary-phase hypothesis.

\subsection{Exact branch integrals}

Set
\[
 a_D(t,P)=-\frac{t}{4P}+\frac{i}{8P^2},\qquad
 a_R(t,P)=-\frac{i}{8P^2},\qquad
 \epsilon_D=-1,\quad \epsilon_R=1.
\]
Throughout the momentum calculation we use the single substitution
\begin{equation}
 p=\alpha-q,\qquad p'=\beta-q,\qquad
 P=|\alpha-q|,\qquad P'=|\beta-q|.
 \label{eq:substitution}
\end{equation}
For \(A,B\in\{D,R\}\), define
\begin{equation}
\begin{split}
 Q_\eps^{AB}(q;\Phi)
 ={}&\int \dd t\,\dd s\,\dd t'\,\dd^3\mathbf y
 \int\frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6}
 \vartheta(\eps P)\vartheta(\eps P')\\
 &{}\times a_A(t,P)\overline{a_B(t',P')}
 e^{i[(\epsilon_A P-\omega_m(q))t
       -(\epsilon_B P'-\omega_m(q))t'
       +(P-P')s]}
 \Ps_\Phi(t,s,t',\mathbf y;\alpha,\beta).
\end{split}
\label{eq:Qbranches}
\end{equation}
For \(\eps=0\), the two profile factors are replaced by one.
For each fixed \(q\), the integral is absolutely convergent after
taking the joint spatial transform.

The exact positive-cutoff identity is
\begin{equation}
 \la T_\eps^\rho,\Phi\ra
 =
 \frac1{(2\pi)^3}\int_{\R^3}
 \omega_m(q)M_\eps(|q|)^2
 \sum_{A,B\in\{D,R\}}Q_\eps^{AB}(q;\Phi)\dd^3q .
 \label{eq:jointbranches}
\end{equation}
The spatial signs follow from
\[
 p\cdot(\mathbf x-\mathbf y)
 -p'\cdot(\mathbf x'-\mathbf y)
 +q\cdot(\mathbf x-\mathbf x')
 =
 \alpha\cdot(\mathbf x-\mathbf y)
 -\beta\cdot(\mathbf x'-\mathbf y).
\]

Define \(S_\eps^\rho\) by taking all four branches in
\eqref{eq:jointbranches} on \(|q|\leq1\), and only
\(DD,DR,RD\) on \(|q|>1\).
Define \(S\) by the same prescription with all multipliers
replaced by one. The high-\(|q|\) integrals defining \(S\) are
understood after performing the compact time integrations in
\eqref{eq:Qbranches}.

\begin{lemma}[Low-frequency and direct sectors]
\label{lem:S}
The distribution \(S\) exists, is independent of \(\rho\), and
satisfies
\begin{equation}
 |\la S_\eps^\rho-S,\Phi\ra|
 \leq C_{K,\rho,m}\eps\,p_{K,64}(\Phi),
 \qquad 0<\eps\leq\tfrac12.
 \label{eq:Srate}
\end{equation}
\end{lemma}

\begin{proof}
The profiles are bounded and globally Lipschitz, with value one
at zero. Telescoping the product of multipliers gives
\begin{equation}
 \left|
 M_\eps(|q|)^2
 \vartheta(\eps P)\vartheta(\eps P')-1
 \right|
 \leq C_{\rho,m}\eps(1+|q|+P+P').
 \label{eq:multdifference}
\end{equation}

For \(|q|\leq1\), \eqref{eq:FTbound} and the bounded-\(r\) part
of Lemma~\ref{lem:convolution} give an integrable bound.
The extra powers in \eqref{eq:multdifference} are harmless:
\(P\leq|\alpha|+1\) and \(P'\leq|\beta|+1\).
Thus this sector has a cutoff-free limit with \(O(\eps)\) error.

For \(r=|q|\geq1\), suppose the left branch is direct.
Its \(t\)-frequency is \(-(P+\omega_m(q))\).
Integrating eight times by parts in \(t\) gains
\[
                         (P+\omega_m(q))^{-8}\leq r^{-8}.
\]
The amplitude is at most linear in \(t\); its differentiated
coefficients are still bounded by
\(C_K(P^{-1}+P^{-2})\).
The other amplitude is bounded by
\(C_K(P'^{-1}+P'^{-2})\).
If only the right branch is direct, use \(t'\) instead.

The shifted bounds show that the two momentum integrations
contribute at most \(Cr^{-2}\). Since
\(r^2\omega_m(r)\leq C_mr^3\), the unregulated radial majorant
is \(Cr^{-7}\).
For the regulated difference,
\[
 P\leq r+|\alpha|,\qquad P'\leq r+|\beta|,
\]
and \eqref{eq:multdifference} costs at most \(C\eps r\)
after the weighted integrations.
The resulting radial majorant is \(C\eps r^{-6}\).
Both are integrable.

All integrations by parts are performed against the full compact
joint test. No restriction to an observation diagonal is used.
\end{proof}

\paragraph{Certificate.}
This lemma holds on every compact in \(X\), including ordinary
diagonals and the reflected intersection. It covers bounded
scalar momentum and every sector with a direct branch. The test
order is \(64\), independent of \(\eps\); the error is
\(O_{K,\rho,m}(\eps)\), without a lower bound on \(|t-t'|\).

\section{The reflected radial family: normalization and Lemma A}
\label{sec:radial}

For \(q=r\boldsymbol\nu\), \(r\geq1\), set
\[
 P=|\alpha-r\boldsymbol\nu|,\qquad
 P'=|\beta-r\boldsymbol\nu|,\qquad
 \omega(r)=\sqrt{r^2+m^2},
\]
in accordance with \eqref{eq:substitution}.

The radial families are defined intrinsically from the exact
branch integrals:
\begin{align}
 \la A(r),\Phi\ra
 &=
 \frac{r^3\omega(r)}{(2\pi)^3}
 \int_{S^2}Q_0^{RR}(r\boldsymbol\nu;\Phi)
                 \Om(\boldsymbol\nu),
 \label{eq:Aintrinsic}\\
 \la A_\eps^\rho(r),\Phi\ra
 &=
 \frac{r^3\omega(r)}{(2\pi)^3}
 \int_{S^2}Q_\eps^{RR}(r\boldsymbol\nu;\Phi)
                 \Om(\boldsymbol\nu).
 \label{eq:Aepsintrinsic}
\end{align}
No residue constant has been inserted in these definitions.
They give the exact decomposition
\begin{equation}
 T_\eps^\rho
 =
 S_\eps^\rho+
 \int_1^\infty M_\eps(r)^2A_\eps^\rho(r)\frac{\dd r}{r}.
 \label{eq:rawsplit}
\end{equation}

The use of \(r\) as the radial integration variable is distinct
from its use for the fixed residue distribution.
In expressions such as \(A(r)-r\), the second term is the
distribution \eqref{eq:rdef}.

\subsection{Normalization ledger}

\begin{proposition}[Normalization ledger]
\label{prop:ledger}
The leading reflected--reflected radial coefficient derived from
\eqref{eq:Aintrinsic} equals the prescribed distribution
\eqref{eq:rdef} on all joint tests. In particular, its constant is
\(c_0=1/(64(2\pi)^3)\), with the parameter density in
\eqref{eq:nu}.
\end{proposition}

\begin{proof}
The calculation is recorded in six identities.

\paragraph{L1: exact pointwise kernel.}
Proposition~\ref{prop:positive} gives
\begin{equation}
 T_\eps^\rho
 =
 2K_\eps(x,y)\overline{K_\eps(x',y)}D_\eps(x-x').
 \tag{L1}\label{ledger:1}
\end{equation}

\paragraph{L2: reflected amplitudes and conjugation.}
The reflected terms of the two mixed covariances have amplitudes
\begin{equation}
 a_R(t,P)=-\frac{i}{8P^2},\qquad
 \overline{a_R(t',P')}=\frac{i}{8P'^2},\qquad
 a_R(t,P)\overline{a_R(t',P')}=\frac1{64P^2P'^2}.
 \tag{L2}\label{ledger:2}
\end{equation}
The corresponding time phases are, separately,
\(e^{iP(t+s)}\) and \(e^{-iP'(t'+s)}\).
The right spatial phase is
\(e^{-ip'\cdot(\mathbf x'-\mathbf y)}\).
These are conjugations of the kernel, not of \(\Phi\).

\paragraph{L3: scalar density.}
The scalar factor is
\begin{equation}
 D_\eps(z)
 =
 \int\frac{\dd^3q}{(2\pi)^3}
 \frac{\omega_m(q)}2M_\eps(|q|)^2
 e^{-i\omega_m(q)z^0+iq\cdot\mathbf z}.
 \tag{L3}\label{ledger:3}
\end{equation}

\paragraph{L4: combined numerical prefactor.}
Consequently the reflected--reflected integrand has prefactor
\begin{equation}
 2\cdot\frac1{64P^2P'^2}\cdot\frac{\omega_m(q)}2
 =
 \frac{\omega_m(q)}{64P^2P'^2},
 \qquad\text{with measure }\frac{\dd^3q}{(2\pi)^3}.
 \tag{L4}\label{ledger:4}
\end{equation}
The assignment factor in L1 has now been exhausted. No further
factor of two occurs in \(A(r)\).

\paragraph{L5: polar measure and radial normalization.}
With ordinary sphere area measure,
\begin{equation}
 \dd^3q=r^2\dd r\,\Om(\boldsymbol\nu),\qquad
 \frac{r^2\omega(r)}{P^2P'^2}\dd r
 =
 \frac{r^3\omega(r)}{P^2P'^2}\frac{\dd r}{r}.
 \tag{L5}\label{ledger:5}
\end{equation}
This explains the extra \(r^3\) in
\eqref{eq:Aintrinsic}, rather than adding a normalization by hand.

\paragraph{L6: translation and Fourier inversion.}
For fixed \(q\), \eqref{eq:substitution} is a translation and has
Jacobian one. The two mixed-covariance measures become exactly
\(\dd^3\alpha\,\dd^3\beta/(2\pi)^6\).
The Fourier inversion identity is
\begin{equation}
\begin{gathered}
 \left|
 \det\frac{\partial(p,p')}{\partial(\alpha,\beta)}
 \right|=1,\qquad
 \frac{\dd^3p\,\dd^3p'}{(2\pi)^6}
 =
 \frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6},\\
 \int\frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6}
 e^{-ia\boldsymbol\nu\cdot\alpha+ib\boldsymbol\nu\cdot\beta}
 \Ps_\Phi(t,s,t',\mathbf y;\alpha,\beta)\\
 =
 \Phi(t,\mathbf y+a\boldsymbol\nu;\
      s,\mathbf y;\
      t',\mathbf y+b\boldsymbol\nu).
\end{gathered}
\tag{L6}\label{ledger:6}
\end{equation}
In particular, the two factors \((2\pi)^{-3}\) supplied by the
mixed covariances are the normalization of the two Fourier
inversions. They cancel the factors \((2\pi)^3\) from the
unnormalized Fourier delta identities; they are not canceled a
second time against another Fourier measure.

For the leading coefficient, replace the amplitude factors
\(P^{-2},P'^{-2}\) by \(r^{-2},r^{-2}\) and \(\omega(r)\) by \(r\).
In the phase, however, retain the first-order expansions
\[
 P=r-\boldsymbol\nu\cdot\alpha,\qquad
 P'=r-\boldsymbol\nu\cdot\beta,\qquad \omega(r)=r.
\]
Substituting \(P=P'=r\) into the phase itself would be incorrect.
The retained time phase is exactly the inversion phase in L6.

The only remaining numerical factor is therefore
\[
                         \frac1{(2\pi)^3}\frac1{64}=c_0.
\]
After L6, the remaining measure is
\(\dd t\,\dd s\,\dd t'\,\dd^3\mathbf y\,\Om(\boldsymbol\nu)\).
Changing \(\mathbf y\) to
\(\mathbf z=\mathbf y+s\boldsymbol\nu\) has unit Jacobian, so this
is precisely the pushforward density \eqref{eq:nu}. We have proved
\begin{equation}
 \boxed{
 \begin{gathered}
 A_{\mathrm{lead}}=r\quad\text{in }\D'(X),
 \qquad c_0=\frac1{64(2\pi)^3},\\
 \text{leading RR contribution}
   =\la r,\Phi\ra\,\frac{\dd r}{r}.
 \end{gathered}
 }
 \label{eq:ledgerconclusion}
\end{equation}
% Polish Block 1: OK.
% 2*(1/64)*(omega/2)=omega/64; the scalar (2*pi)^(-3)
% remains after the two spatial Fourier inversions.
% The translations have Jacobian one and give the positive ray signs.

The cutoff factors
\(M_\eps(r)^2\vartheta(\eps P)\vartheta(\eps P')\)
have not been used in this leading identification.
They are treated by the replacement estimate below.
Direct-branch amplitudes do not enter this calculation.
\end{proof}

As a consequence of the ledger, the intrinsic definition
\eqref{eq:Aintrinsic} can now be written
\begin{equation}
\begin{split}
 \la A(r),\Phi\ra
 ={}&c_0r^3\omega(r)
 \int_{S^2}\Om(\boldsymbol\nu)
 \int\dd t\,\dd s\,\dd t'\,\dd^3\mathbf y\\
 &{}\times
 \int\frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6P^2P'^2}
 e^{i[(P-\omega(r))t-(P'-\omega(r))t'+(P-P')s]}
 \Ps_\Phi(t,s,t',\mathbf y;\alpha,\beta).
\end{split}
\label{eq:Adef}
\end{equation}
The formula for \(A_\eps^\rho(r)\) has the two additional factors
\(\vartheta(\eps P)\vartheta(\eps P')\).

\subsection{Lemma A: exact phase and complete error budget}

\begin{lemma}[Joint-test radial asymptotic: Lemma A]
\label{lem:A}
For every compact \(K\subset X\),
\begin{equation}
 |\la A(r)-r,\Phi\ra|
 \leq C_{K,m}r^{-1}p_{K,64}(\Phi),
 \qquad r\geq1.
 \label{eq:Ar}
\end{equation}
In particular,
\[
 |\la A(r),\Phi\ra|\leq C_{K,m}p_{K,64}(\Phi).
\]
The estimate holds without a positive lower bound on
\(|t-t'|\).
\end{lemma}

\begin{proof}
We retain the exact phase, form an exact difference, and only then
estimate its terms.

\paragraph{1. Exact phase.}
Set
\[
 \delta_r^{\boldsymbol\nu}(\alpha)
 =
 |\alpha-r\boldsymbol\nu|-r+\boldsymbol\nu\cdot\alpha,
 \qquad
 \varphi_0=-a\boldsymbol\nu\cdot\alpha
                  +b\boldsymbol\nu\cdot\beta.
\]
The temporal phase in \eqref{eq:Adef} is exactly
\begin{equation}
 \varphi_0
 +a\delta_r^{\boldsymbol\nu}(\alpha)
 -b\delta_r^{\boldsymbol\nu}(\beta)
 +(r-\omega(r))d.
 \label{eq:timephase}
\end{equation}

Before the spatial integration defining \(\Ps_\Phi\), the exact
full phase in the incidence chart is
\begin{equation}
\begin{split}
 \Theta_r={}&
 (a\alpha-b\beta)\cdot(\mathbf n-\boldsymbol\nu)
 +\sigma\mathbf n\cdot\left(\alpha-\frac ba\beta\right)
 -\beta\cdot\mathbf w\\
 &{}+a\delta_r^{\boldsymbol\nu}(\alpha)
       -b\delta_r^{\boldsymbol\nu}(\beta)
       +(r-\omega(r))d.
\end{split}
\label{eq:intersectionphase}
\end{equation}
Thus both angular variables and all of
\((\sigma,\mathbf w,d)\) are retained.
The phase \(\varphi_0\) is the leading \emph{temporal} phase in
Cartesian Fourier space; the spatial phase is already inside
\(\Ps_\Phi\). One does not delete the first line of
\eqref{eq:intersectionphase} in an incidence-coordinate integral.

\paragraph{2. Exact central--tail decomposition.}
Let
\[
 \mathcal C_r=\{|\alpha|\leq r/2,\ |\beta|\leq r/2\},
 \qquad
 \lambda_r=\frac{\omega(r)}r,\qquad
 J_r=\frac{r^4}{P^2P'^2},
\]
and abbreviate, within this proof,
\[
 u_\delta=a\delta_r^{\boldsymbol\nu}(\alpha)
                -b\delta_r^{\boldsymbol\nu}(\beta),
 \qquad
 u_\omega=(r-\omega(r))d.
\]
All integrals over \(\mathcal C_r\) or its complement below include
the remaining variables, with measure
\[
 \Om(\boldsymbol\nu)\,
 \dd t\,\dd s\,\dd t'\,\dd^3\mathbf y\,
 \frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6}.
\]
Write \(\Ps_\Phi\) at its full arguments from
\eqref{eq:jointFT}. Define
\begin{align}
 I_{\mathrm{cent}}^{\mathrm{lead}}
 &=c_0\int_{\mathcal C_r}e^{i\varphi_0}\Ps_\Phi,
 \label{eq:Icentral}\\
 E_P
 &=c_0\int_{\mathcal C_r}
 e^{i\varphi_0}\Ps_\Phi\,
 \lambda_r(J_r-1)e^{i(u_\delta+u_\omega)},
 \label{eq:EPdef}\\
 E_\delta
 &=c_0\int_{\mathcal C_r}
 e^{i\varphi_0}\Ps_\Phi\,
 \lambda_r(e^{iu_\delta}-1)e^{iu_\omega},
 \label{eq:EdeltaDef}\\
 E_\omega
 &=c_0\int_{\mathcal C_r}
 e^{i\varphi_0}\Ps_\Phi\,
 \bigl(\lambda_re^{iu_\omega}-1\bigr),
 \label{eq:EomegaDef}\\
 T_{\mathrm{tail}}
 &=c_0\int_{\mathcal C_r^c}
 e^{i\varphi_0}\Ps_\Phi\,
 \lambda_rJ_re^{i(u_\delta+u_\omega)},\qquad
 R_{\mathrm{tail}}^{\mathrm{lead}}
 =c_0\int_{\mathcal C_r^c}e^{i\varphi_0}\Ps_\Phi.
 \label{eq:taildefs}
\end{align}
These definitions precede all absolute-value estimates.
The algebraic identity behind them is
\begin{equation}
\begin{split}
 \lambda_rJ_re^{i(u_\delta+u_\omega)}-1
 ={}&
 \lambda_r(J_r-1)e^{i(u_\delta+u_\omega)}
 +\lambda_r(e^{iu_\delta}-1)e^{iu_\omega}\\
 &+\bigl(\lambda_re^{iu_\omega}-1\bigr).
\end{split}
\label{eq:exacterroridentity}
\end{equation}
Consequently,
\begin{align}
 \la A(r),\Phi\ra
 &=
 I_{\mathrm{cent}}^{\mathrm{lead}}
 +E_P+E_\delta+E_\omega+T_{\mathrm{tail}},
 \label{eq:Acentralexact}\\
 I_{\mathrm{cent}}^{\mathrm{lead}}
 &=
 \la r,\Phi\ra-R_{\mathrm{tail}}^{\mathrm{lead}}.
 \label{eq:Icentralexact}
\end{align}
The second equality is the full Fourier inversion proved in the
ledger, minus its leading tail. In particular, the central leading
term is not counted again after subtracting \(r\).
With
\[
 E_{\mathrm{tail}}
 =T_{\mathrm{tail}}-R_{\mathrm{tail}}^{\mathrm{lead}},
\]
the exact difference is
\begin{equation}
 \la A(r)-r,\Phi\ra
 =E_P+E_\delta+E_\omega+E_{\mathrm{tail}}.
 \label{eq:budgetidentity}
\end{equation}
% Polish Block 2: OK.
% Expanding the three central errors gives lambda*J*exp(i(u_delta+u_omega))-1.
% E_omega has no J; subtracting the leading tail gives exactly the stated sum.

\paragraph{3. Amplitude error \(E_P\).}
On \(\mathcal C_r\),
\[
 \frac r2\leq P,P'\leq\frac{3r}2,\qquad
 \left|\frac{r^2}{P^2}-1\right|
 =\frac{|r-P|(r+P)}{P^2}
 \leq C\frac{|\alpha|}{r}.
\]
The corresponding bound holds for \(P'\). Since each ratio
\(r^2/P^2\) is bounded by \(4\),
\begin{equation}
 |J_r-1|
 \leq C\frac{|\alpha|+|\beta|}{r},
 \qquad J_r\leq16.
 \label{eq:amplitudeerror}
\end{equation}
Also \(\lambda_r\leq\sqrt{1+m^2}\) for \(r\geq1\).
Using \eqref{eq:FTbound} in \eqref{eq:EPdef} gives
\begin{equation}
\begin{split}
 |E_P|
 &\leq
 \frac{C_{K,m}}r p_{K,64}(\Phi)
 \int_{\R^6}
 \frac{|\alpha|+|\beta|}
 {(1+|\alpha|)^{12}(1+|\beta|)^{12}}
 \dd^3\alpha\,\dd^3\beta\\
 &\leq \frac{C'_{K,m}}r p_{K,64}(\Phi).
\end{split}
\label{eq:EPbound}
\end{equation}

\paragraph{4. Quadratic phase error \(E_\delta\).}
For \(|\alpha|\leq r/2\), the exact identity
\begin{equation}
 \delta_r^{\boldsymbol\nu}(\alpha)
 =
 \frac{|\alpha|^2-(\boldsymbol\nu\cdot\alpha)^2}
 {|\alpha-r\boldsymbol\nu|+r-\boldsymbol\nu\cdot\alpha}
 \quad\text{implies}\quad
 |\delta_r^{\boldsymbol\nu}(\alpha)|
 \leq\frac{|\alpha|^2}{r}.
 \label{eq:deltaexact}
\end{equation}
The denominator is at least \(r\) in this region.
On a compact observation set, \(a,b\) have fixed upper bounds.
Therefore
\[
 |u_\delta|
 \leq \frac{C_K}{r}(|\alpha|^2+|\beta|^2).
\]
Applying \(|e^{iu}-1|\leq|u|\) to the exact difference in
\eqref{eq:EdeltaDef}, and then \eqref{eq:FTbound}, yields
\begin{equation}
\begin{split}
 |E_\delta|
 &\leq \frac{C_{K,m}}r p_{K,64}(\Phi)
 \int_{\R^6}
 \frac{|\alpha|^2+|\beta|^2}
 {(1+|\alpha|)^{12}(1+|\beta|)^{12}}
 \dd^3\alpha\,\dd^3\beta\\
 &\leq \frac{C'_{K,m}}r p_{K,64}(\Phi).
\end{split}
\label{eq:EdeltaBound}
\end{equation}
No expansion of \(\delta_r^{\boldsymbol\nu}\) is made in the tail.

\paragraph{5. Mass and time-difference error \(E_\omega\).}
The two relevant bounds are
\begin{equation}
 \omega(r)-r=\frac{m^2}{\omega(r)+r}
 \leq\frac{m^2}{2r},\qquad
 \lambda_r-1=\frac{m^2}{r(\omega(r)+r)}
 \leq\frac{m^2}{2r^2}.
 \label{eq:massbounds}
\end{equation}
The exact factor in \eqref{eq:EomegaDef} splits as
\[
 \lambda_re^{iu_\omega}-1
 =
 (\lambda_r-1)e^{iu_\omega}+(e^{iu_\omega}-1).
\]
In particular,
\begin{equation}
 |u_\omega|
 =|(r-\omega(r))d|
 \leq \frac{m^2|d|}{2r}.
 \label{eq:dmass}
\end{equation}
Only an upper bound on \(|d|\) over the compact observation
projection is used. At \(d=0\), this phase error vanishes exactly.
Equations \eqref{eq:FTbound} and
\eqref{eq:massbounds}--\eqref{eq:dmass} give
\begin{equation}
 |E_\omega|
 \leq C_{K,m}(r^{-1}+r^{-2})p_{K,64}(\Phi).
 \label{eq:EomegaBound}
\end{equation}
This is the required estimate through \(\Delta_{13}\); it does not
use division by \(d\).

\paragraph{6. Exact tail \(E_{\mathrm{tail}}\).}
Let \(N=12\). The union bound for
\(\mathcal C_r^c\), together with
\eqref{eq:shiftbound}--\eqref{eq:shifttail} for \(\ell=2,k=0\),
gives the complete majorant
\begin{align}
 |T_{\mathrm{tail}}|
 &\leq C_{K,m}p_{K,64}(\Phi)\,
 \underbrace{r^4}_{\text{prefactor}}
 \left(
 \underbrace{r^{1-N}}_{\alpha\text{-tail}}
 \underbrace{r^{-2}}_{\text{full }\beta}
 +
 \underbrace{r^{-2}}_{\text{full }\alpha}
 \underbrace{r^{1-N}}_{\beta\text{-tail}}
 \right)\notag\\
 &\leq C'_{K,m}r^{3-N}p_{K,64}(\Phi)
 =C'_{K,m}r^{-9}p_{K,64}(\Phi).
 \label{eq:actualtail}
\end{align}
Here the prefactor bound is \(r^3\omega(r)\leq C_mr^4\).
The shifted-tail exponent is \(3-2+0-N=1-N\), whereas
the full shifted integral is bounded by \(Cr^{-2}\).

This estimate includes \(P\approx0\) and \(P'\approx0\), which
occur in the shifted tails. Those singularities are treated by
local integrability, not by a phase expansion.

For the leading tail there are no shifted denominators. Directly,
\[
 \int_{|\alpha|>r/2}(1+|\alpha|)^{-N}\dd^3\alpha
 \leq Cr^{3-N},
\]
and similarly for \(\beta\). Hence
\begin{equation}
 |R_{\mathrm{tail}}^{\mathrm{lead}}|
 +|E_{\mathrm{tail}}|
 \leq C_{K,m}r^{3-N}p_{K,64}(\Phi)
 =C_{K,m}r^{-9}p_{K,64}(\Phi).
 \label{eq:tailbudget}
\end{equation}
% Polish Block 3: OK.
% Actual tail: r^4 * r^(1-N) * r^(-2) = r^(3-N) = r^(-9).
% Leading tail: no shifted denominator, also O(r^(3-N)).
% N=12 and test order 64 are unchanged.

\paragraph{7. Coordinate change: \(E_{\mathrm{chart}}=0\).}
The estimate has been proved in Cartesian Fourier space.
On a fixed incidence-coordinate patch, rewriting the same
distribution with \eqref{eq:chart} introduces exactly the density
\eqref{eq:density} and the phase \eqref{eq:intersectionphase}.
It introduces no additional error:
\[
                              E_{\mathrm{chart}}=0.
\]
For a direct justification, first insert compact smooth cutoffs
in \(\alpha,\beta\). The remaining momentum singularities
\(P^{-2},P'^{-2}\) are locally integrable, so the ordinary
coordinate change is legitimate. Remove those auxiliary cutoffs
after the spatial integration defining \(\Ps_\Phi\);
\eqref{eq:FTbound} and Lemma~\ref{lem:convolution} justify that
limit. This identifies the coordinate expression with the original
distribution, without restricting it to a diagonal.
The Jacobian has no \(1/d\).

\paragraph{8. Sum of the errors.}
The complete budget is
\[
\begin{array}{c|c}
 \text{term} & \text{bound divided by }p_{K,64}(\Phi)\\
 \hline
 E_P & C_{K,m}r^{-1}\\
 E_\delta & C_{K,m}r^{-1}\\
 E_\omega & C_{K,m}(r^{-1}+r^{-2})\\
 E_{\mathrm{tail}} & C_{K,m}r^{-9}\\
 E_{\mathrm{chart}} & 0
\end{array}
\]
Summing in \eqref{eq:budgetidentity} proves \eqref{eq:Ar}.
The uniform bound on \(A(r)\) follows because \(r\) is an order-zero
distribution.
\end{proof}

\paragraph{Certificate.}
Lemma A holds on every compact of \(X\), with test order \(64\),
independently of the regulator. The proof covers the comparable
central sector and all noncomparable shifted tails, uniformly in
the momentum direction. It uses only upper bounds on \(a,b,|d|\).
Thus it applies unchanged on neighborhoods of
\(\Rel\cap\Delta_{13}\).

The estimate is an equality and bound on all joint tests.
In particular, an additional order-\(r^0\) contribution proportional
to \(\delta(d)\), a normal derivative of the residue, or any other
nonzero distribution could not be hidden in its error:
it would contradict \eqref{eq:Ar} on a test detecting that
distribution.

\section{Replacement and the complete ultraviolet split}
\label{sec:UV}

\subsection{Replacement of the two graviton multipliers}

\begin{lemma}[Insertion-regulator replacement]
\label{lem:replacement}
For every compact \(K\subset X\),
\begin{equation}
 \left|
 \la A_\eps^\rho(r)-\vartheta(\eps r)^2A(r),\Phi\ra
 \right|
 \leq C_{K,\rho,m}\eps\,p_{K,64}(\Phi),
 \qquad r\geq1,\quad 0<\eps\leq1.
 \label{eq:Aregerror}
\end{equation}
Only boundedness and a global Lipschitz bound for \(\vartheta\)
are needed for this estimate.
\end{lemma}

\begin{proof}
\paragraph{1. Exact multiplier difference.}
Before taking absolute values, use the identity
\begin{equation}
\begin{split}
 &\vartheta(\eps P)\vartheta(\eps P')-\vartheta(\eps r)^2\\
 &\quad=
 \vartheta(\eps P')
       \bigl(\vartheta(\eps P)-\vartheta(\eps r)\bigr)
 +\vartheta(\eps r)
       \bigl(\vartheta(\eps P')-\vartheta(\eps r)\bigr).
\end{split}
\label{eq:multtelescoping}
\end{equation}
Each difference has the exact integral remainder
\begin{equation}
 \vartheta(\eps P)-\vartheta(\eps r)
 =
 \eps(P-r)\int_0^1
 \vartheta'\bigl(\eps(r+\tau(P-r))\bigr)\dd\tau,
 \label{eq:meanvalueexact}
\end{equation}
and likewise with \(P'\).

Using the phase notation from the proof of Lemma A, inserting the
left-hand side of \eqref{eq:multtelescoping} into \eqref{eq:Adef}
gives exactly
\begin{equation}
\begin{split}
 &\la A_\eps^\rho(r)-\vartheta(\eps r)^2A(r),\Phi\ra\\
 &\quad=
 c_0r^3\omega(r)
 \int
 \frac{\vartheta(\eps P)\vartheta(\eps P')
             -\vartheta(\eps r)^2}{P^2P'^2}
 e^{i[\varphi_0+u_\delta+u_\omega]}\Ps_\Phi\\
 &\qquad\qquad{}\times
 \Om(\boldsymbol\nu)\,
 \dd t\,\dd s\,\dd t'\,\dd^3\mathbf y\,
 \frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6}.
\end{split}
\label{eq:replacementexact}
\end{equation}
Here the phase is the same exact phase used in Lemma A.
No approximation of it is made.

Since
\[
 |P-r|\leq|\alpha|,\qquad |P'-r|\leq|\beta|,
\]
equations \eqref{eq:multtelescoping}--\eqref{eq:meanvalueexact}
imply
\begin{equation}
 \left|
 \vartheta(\eps P)\vartheta(\eps P')-\vartheta(\eps r)^2
 \right|
 \leq
 \|\vartheta\|_\infty\|\vartheta'\|_\infty
 \eps(|\alpha|+|\beta|).
 \label{eq:multLipschitz}
\end{equation}
This uses no decay when \(P\) or \(P'\) is small.

\paragraph{2. Central region.}
On the region \(\mathcal C_r\) of Lemma A,
\(r^3\omega(r)/(P^2P'^2)\leq C_m\).
Taking absolute values in the already formed difference
\eqref{eq:replacementexact} gives
\begin{equation}
\begin{split}
 \left|\text{central part of \eqref{eq:replacementexact}}\right|
 &\leq C_{K,\rho,m}\eps\,p_{K,64}(\Phi)
 \int_{\R^6}
 \frac{|\alpha|+|\beta|}
 {(1+|\alpha|)^{12}(1+|\beta|)^{12}}
 \dd^3\alpha\,\dd^3\beta\\
 &\leq C'_{K,\rho,m}\eps\,p_{K,64}(\Phi).
\end{split}
\label{eq:replacementcentral}
\end{equation}

\paragraph{3. Shifted tails, including \(P\approx0\).}
Use the same complement \(\mathcal C_r^c\) as in Lemma A,
and retain \(N=12\).
For the \(\alpha\)-tail, apply \eqref{eq:shifttail} to the
\(\alpha\)-integral and \eqref{eq:shiftbound} to the full
\(\beta\)-integral. The extra momentum factor is absorbed using
\(|\alpha|\leq1+|\alpha|\) or \(|\beta|\leq1+|\beta|\).
The following table records the parameters of those two estimates.
Every row includes the prefactor
\(r^3\omega(r)\leq C_mr^4\); the common factor
\(C_{K,\rho,m}p_{K,64}(\Phi)\) is suppressed.

\begin{center}
\small
\renewcommand{\arraystretch}{1.25}
\begin{tabularx}{\linewidth}{
@{}c
>{\centering\arraybackslash}X
>{\centering\arraybackslash}X
>{\centering\arraybackslash}X@{}}
\toprule
Extra factor &
\(\alpha\)-tail: \((k,\ell)\) and bound &
Full \(\beta\): \((k,\ell)\) and bound &
Radial majorant\\
\midrule
\(\eps|\alpha|\) &
\((1,2):\ Cr^{2-N}\) &
\((0,2):\ Cr^{-2}\) &
\(C\eps r^{4-N}\)\\
\(\eps|\beta|\) &
\((0,2):\ Cr^{1-N}\) &
\((1,2):\ Cr^{-2}\) &
\(C\eps r^{3-N}\)\\
\bottomrule
\end{tabularx}
\end{center}

Indeed, the tail powers are
\(3-2+1-N=2-N\) and \(3-2+0-N=1-N\).
The full-integral estimate with \(k=1,\ell=2\) still gives
\(Cr^{-2}\), since \(N=12>1+3\); no additional power of \(r\)
is lost there.
The \(\beta\)-tail has the same two bounds after interchanging
\(\alpha\) and \(\beta\). The union bound therefore yields
\begin{equation}
\begin{split}
 \left|\text{tail part of \eqref{eq:replacementexact}}\right|
 &\leq
 C_{K,\rho,m}\eps(r^{4-N}+r^{3-N})p_{K,64}(\Phi)\\
 &=C_{K,\rho,m}\eps(r^{-8}+r^{-9})p_{K,64}(\Phi)\\
 &\leq C'_{K,\rho,m}\eps\,p_{K,64}(\Phi).
\end{split}
\label{eq:replacementtail}
\end{equation}
% Polish Block 4: OK.
% Weighted alpha-tail: eps*r^4*r^(2-N)*r^(-2)=eps*r^(4-N).
% Cross-weighted tail: eps*r^4*r^(1-N)*r^(-2)=eps*r^(3-N).
% The full k=1, ell=2 bound remains O(r^(-2)).
% At N=12 both tails are O(eps*r^(-8)) or better; Lipschitz only.

The shifted-tail bounds include the locally integrable singularities
at \(P=0\) and \(P'=0\).
Combining \eqref{eq:replacementcentral} and
\eqref{eq:replacementtail} proves the claim.
\end{proof}

\paragraph{Certificate.}
The replacement estimate holds on all of \(K\), including
\(\Rel\cap\Delta_{13}\), with order \(64\), uniformly in \(r\).
Its central error is \(O(\eps)\), and its tail is
\(O(\eps(r^{-8}+r^{-9}))\).
No small-argument decay of a graviton multiplier is used.
The estimate does not replace any multiplier in the direct-sector
distribution \(S_\eps^\rho\), and it does not yet replace the scalar
multiplier \(M_\eps\).

\subsection{The ultraviolet decomposition}

\begin{proposition}[Regulator-independent ultraviolet decomposition]
\label{prop:UV}
The exact family has the decomposition
\begin{equation}
 \boxed{
 T_\eps^\rho
 =
 S+\int_1^\infty\vartheta(\eps r)^4A(r)\frac{\dd r}{r}
   +E_\eps^\rho ,
 }
 \label{eq:UV}
\end{equation}
where \(S,A(r)\) are independent of \(\rho\), and
\begin{equation}
 |\la E_\eps^\rho,\Phi\ra|
 \leq C_{K,\rho,m}\eps(1+|\log\eps|)p_{K,64}(\Phi),
 \qquad 0<\eps\leq\tfrac12.
 \label{eq:Eerror}
\end{equation}
\end{proposition}

\begin{proof}
The ledger identifies the leading coefficient with \(r\), and
Lemma A gives an integrable radial remainder \(A(r)-r\).
For the finite-cutoff identity, subtract the radial term in
\eqref{eq:UV} from \eqref{eq:rawsplit}.
The exact error is the sum of \(S_\eps^\rho-S\) and
\begin{align}
 &\int_1^\infty
 M_\eps(r)^2
 \bigl(A_\eps^\rho(r)-\vartheta(\eps r)^2A(r)\bigr)
 \frac{\dd r}{r},
 \label{eq:errorfirst}\\
 &\int_1^\infty
 \vartheta(\eps r)^2
 \bigl(M_\eps(r)^2-\vartheta(\eps r)^2\bigr)
 A(r)\frac{\dd r}{r}.
 \label{eq:errorsecond}
\end{align}

The first difference is \(O(\eps)\) by Lemma~\ref{lem:S}.
For \eqref{eq:errorfirst}, Lemma~\ref{lem:replacement} and
\[
 |M_\eps(r)|\leq C_{\rho,N}(1+\eps r)^{-N}
\]
give the bound
\[
 C_{K,\rho,m}\eps p_{K,64}(\Phi)
 \int_1^\infty(1+\eps r)^{-N}\frac{\dd r}{r}
 \leq
 C'_{K,\rho,m}\eps(1+|\log\eps|)p_{K,64}(\Phi).
\]
It is the scalar multiplier that supplies this large-\(r\) decay.

The scalar replacement is a separate identity and estimate:
\begin{equation}
\begin{split}
 \sqrt{r^2+m^2/2}-r
 &=
 \frac{m^2/2}{\sqrt{r^2+m^2/2}+r}
 \leq\frac{m^2}{4r},\\
 |M_\eps(r)-\vartheta(\eps r)|
 &\leq \frac{m^2}{4}\|\vartheta'\|_\infty\eps r^{-1}.
\end{split}
\label{eq:scalarreplacement}
\end{equation}
The boundedness of the profiles and of \(A(r)\) therefore makes
\eqref{eq:errorsecond} \(O(\eps)\), since
\(\int_1^\infty r^{-2}\dd r<\infty\).
All estimates use the same compact-test seminorm.
\end{proof}

\section{Intersection and sector audit}
\label{sec:audit}

\subsection{The full intersection pairing}

For a test supported in an incidence-coordinate patch, the exact
pairing corresponding to \eqref{eq:Adef} is
\begin{equation}
\begin{split}
 \la A(r),\Phi\ra
 ={}&c_0r^3\omega(r)\,
 \Os\!\int
 \frac{(a+\sigma)^2\widetilde\Phi\,e^{i\Theta_r}}{P^2P'^2}\\
 &{}\times
 \dd^3\mathbf z\,\Om(\mathbf n)\,
 \dd t\,\dd s\,\dd d\,\dd\sigma\,\dd^3\mathbf w\,
 \Om(\boldsymbol\nu)\,
 \frac{\dd^3\alpha\,\dd^3\beta}{(2\pi)^6},
\end{split}
\label{eq:intersectionpairing}
\end{equation}
where \(t'=t-d>0\) and \(\Theta_r\) is
\eqref{eq:intersectionphase}.
The oscillatory notation means the Cartesian-Fourier distribution
\eqref{eq:Adef} rewritten by the cutoff-and-coordinate procedure
in Lemma A. It is not a new product of singular distributions.

The leading Fourier inversion imposes
\[
 \mathbf x-\mathbf y=a\boldsymbol\nu,\qquad
 \mathbf x'-\mathbf y=b\boldsymbol\nu.
\]
Only after that inversion may one conclude
\[
 \boldsymbol\nu=\mathbf n,\qquad \sigma=0,\qquad\mathbf w=0.
\]
The positive radius \(a=t+s\), not \(a-b=d\), determines the
first angular constraint. Neither this identification nor its
density requires division by \(d\).

\subsection{Normal jets as a consequence, not a second proof}

On a fixed patch, write \(v=(\sigma,\mathbf w)\), and let \(u\)
denote all tangential variables, including \(d\).
A fixed local cutoff and the integral Taylor formula give
\[
 \widetilde\Phi(u,v)
 =
 \chi(v)\widetilde\Phi(u,0)+\sum_{j=1}^4v_j\Phi_j(u,v),
\]
with smooth compactly supported \(\Phi_j\) and finite-order
seminorm bounds on a slightly larger patch.
The residue kills the normal-vanishing terms.
Lemma A already controls those terms and the trace term uniformly
through \(d=0\).
This observation does not restrict \(A(r)\) to the relation and
is not used as an alternative proof of the radial estimate.

\subsection{Observation coverage}

The estimates are global in the Cartesian observation variables,
so they cover the following regions.

\begin{enumerate}[label=\textup{(\roman*)}]
\item
Off \(\Rel\), the residue vanishes on supported tests.
Equations \eqref{eq:UV} and \eqref{eq:Ar} then give convergence
without a logarithmic subtraction.

\item
Near \(\Rel\setminus\Delta_{13}\), the same estimates apply
without using a lower bound on \(|d|\).

\item
Near \(\Rel\cap\Delta_{13}\), the exact phase and density are
\eqref{eq:intersectionphase} and \eqref{eq:density}.
The mass-phase bound \eqref{eq:dmass} is uniform through \(d=0\).

\item
The remaining ordinary diagonals and the region
\(\mathbf x=\mathbf y\), where the incidence chart is unavailable,
are covered by the Cartesian Fourier argument.
Moreover,
\(\Rel\cap\Delta_{12}=\Rel\cap\Delta_{23}=\varnothing\) in \(X\),
because the corresponding reflected spatial distances are
\(t+s>0\) and \(t'+s>0\).
\end{enumerate}

These are assertions on joint tests supported in open observation
neighborhoods, not tests confined to positive-codimension strata.
A finite smooth partition on any compact gives the corresponding
regional estimates with the same test order.

\subsection{Integration-sector coverage}

\begin{center}
\small
\renewcommand{\arraystretch}{1.14}
\begin{tabularx}{\linewidth}{
@{}>{\raggedright\arraybackslash}p{0.32\linewidth}
>{\raggedright\arraybackslash}X@{}}
\toprule
Sector & Controlling calculation\\
\midrule
All frequencies large and comparable &
The normalization ledger and the central errors
\(E_P,E_\delta,E_\omega\).\\

Only a proper subset large &
Bounded \(q\) belongs to Lemma~\ref{lem:S}.
For large \(q\), a small graviton momentum forces a shifted tail,
covered by \eqref{eq:shifttail}.\\

Ratios tending to zero or infinity &
The shifted-tail bounds include both limits and
\(P=0,P'=0\), with their locally integrable denominators.\\

A direct branch present &
Time integration by parts in a frequency
\(P+\omega(r)\) or \(P'+\omega(r)\), independently of \(d\).\\

Angular degeneracy &
Uniform estimates on \(S^2\), not angular stationary phase.
The two directions in \(\Theta_r\) are kept distinct.\\

A graviton cutoff argument tends to zero &
The exact telescoping and mean-value identities
\eqref{eq:multtelescoping}--\eqref{eq:meanvalueexact};
only boundedness and Lipschitz continuity are used.\\

Bounded frequencies &
The common distribution \(S\), with \(O(\eps)\) convergence.\\
\bottomrule
\end{tabularx}
\end{center}

The partitions at \(|q|=1\) and
\(|\alpha|,|\beta|=r/2\) are independent of \(\eps\).
No derivative of a cutoff-dependent partition is omitted.

\paragraph{Scope audit.}
This numerical subtraction leaves the VII-C coefficient-continuity
obstruction under its stated normal-topology hypotheses unchanged
\cite[\texttt{thm:coefficient}]{VIIC}.

\section{Strong finite part and dyadic certificates}
\label{sec:finitepart}

Define the real profile constant
\begin{equation}
 \kappa_\rho
 =
 \int_0^1\frac{\vartheta(u)^4-1}{u}\dd u
 +\int_1^\infty\frac{\vartheta(u)^4}{u}\dd u.
 \label{eq:kappa}
\end{equation}
The profile is real, even, Schwartz, and equals one at zero.
Hence \(\vartheta(u)^4-1=O(u^2)\), and both integrals converge.

Define the common distribution
\begin{equation}
 U_m
 =
 S+\int_1^\infty\bigl(A(r)-r\bigr)\frac{\dd r}{r}.
 \label{eq:Ucommon}
\end{equation}
This integral is defined by its joint-test pairings.
They are measurable in \(r\), and Lemma A gives the integrable
bound
\[
 \left|\la A(r)-r,\Phi\ra\right|\frac1r
 \leq C_{K,m}r^{-2}p_{K,64}(\Phi).
\]
Thus \eqref{eq:Ucommon} defines a continuous distribution.
Its finite radial truncations converge strongly, with tail at
most \(C_{K,m}R^{-1}p_{K,64}(\Phi)\) above radius \(R\).

\begin{theorem}[Strong finite part with quantitative remainder]
\label{thm:finitepart}
For every admissible fixed mollifier \(\rho\) and every \(\mu>0\),
\begin{equation}
 \boxed{
 u_\mu^\rho=U_m+(\kappa_\rho+\log\mu)\,r
 }
 \label{eq:finitepart}
\end{equation}
is the strong distributional limit
\[
 T_\eps^\rho-\log\!\left(\frac1{\mu\eps}\right)r
 \longrightarrow u_\mu^\rho
 \qquad\text{in }\D'(X).
\]
More precisely,
\begin{equation}
\begin{split}
 &\left|
 \left\langle
 T_\eps^\rho-\log\!\left(\frac1{\mu\eps}\right)r-u_\mu^\rho,
 \Phi
 \right\rangle
 \right|\\
 &\qquad\leq
 C_{K,\rho,m}\eps(1+|\log\eps|)p_{K,64}(\Phi),
 \qquad 0<\eps\leq\tfrac12.
\end{split}
\label{eq:finalrate}
\end{equation}
The estimate includes compacts meeting
\(\Rel\cap\Delta_{13}\).
\end{theorem}

\begin{proof}
The ultraviolet decomposition gives
\begin{align}
 T_\eps^\rho
 ={}&S+
 \int_1^\infty\vartheta(\eps r)^4
               \bigl(A(r)-r\bigr)\frac{\dd r}{r}\notag\\
 &{}+
 \left(\int_1^\infty\vartheta(\eps r)^4\frac{\dd r}{r}\right)r
 +E_\eps^\rho.
 \label{eq:subtracteddecomp}
\end{align}
Changing variables in the scalar radial integral gives
\begin{equation}
\begin{split}
 \int_1^\infty\vartheta(\eps r)^4\frac{\dd r}{r}
 &=
 \log(\eps^{-1})+\kappa_\rho
 -\int_0^\eps\frac{\vartheta(u)^4-1}{u}\dd u\\
 &=
 \log(\eps^{-1})+\kappa_\rho+O_\rho(\eps^2).
\end{split}
\label{eq:scalarfinite}
\end{equation}
Furthermore,
\[
 |\vartheta(\eps r)^4-1|
 \leq C_\rho\min\{(\eps r)^2,1\}.
\]
Lemma A therefore bounds the difference between the first radial
integral in \eqref{eq:subtracteddecomp} and the integral in
\eqref{eq:Ucommon} by
\[
 C_{K,\rho,m}p_{K,64}(\Phi)
 \int_1^\infty
 \frac{\min\{(\eps r)^2,1\}}{r^2}\dd r
 \leq C'_{K,\rho,m}\eps p_{K,64}(\Phi).
\]
Together with \eqref{eq:Eerror} and \eqref{eq:scalarfinite}, this
proves \eqref{eq:finitepart}--\eqref{eq:finalrate}.
The compact-test estimate implies uniform convergence on every
bounded subset of \(\D(X)\).
\end{proof}

\subsection{Shell and within-shell bounds}

Let \(\eps_n=2^{-n}\), \(n\geq1\), and define
\begin{align}
 E_n^\rho
 &=T_{\eps_{n+1}}^\rho-T_{\eps_n}^\rho-(\log2)r,
 \label{eq:dyadicE}\\
 J_{n,\eta}^\rho
 &=T_{2^{-(n+\eta)}}^\rho-T_{2^{-n}}^\rho
       -\eta(\log2)r,\qquad 0\leq\eta\leq1.
 \label{eq:dyadicJ}
\end{align}
The positive-cutoff integrals are continuous in the cutoff,
locally uniformly on bounded sets of tests, so \(J_{n,\eta}^\rho\)
is continuous in \(\eta\).

\begin{corollary}[Summable shell certificate]
\label{cor:shell}
For every compact \(K\subset X\),
\begin{align}
 |\la E_n^\rho,\Phi\ra|
 &\leq a_{n,K}p_{K,64}(\Phi),\\
 \sup_{0\leq\eta\leq1}|\la J_{n,\eta}^\rho,\Phi\ra|
 &\leq b_{n,K}p_{K,64}(\Phi),
\end{align}
where one may take
\begin{equation}
 a_{n,K}=b_{n,K}=C_{K,\rho,m}(n+1)2^{-n}.
 \label{eq:shellbounds}
\end{equation}
In particular, \(\sum_n a_{n,K}<\infty\) and \(b_{n,K}\to0\).
\end{corollary}

\begin{proof}
Apply \eqref{eq:finalrate} with \(\mu=1\) to the two terms in
each difference. For the within-shell estimate,
\(2^{-(n+\eta)}\in[2^{-n-1},2^{-n}]\), so the bound is uniform
in \(\eta\).
\end{proof}

\paragraph{Regional certificate.}
Localizing by any fixed smooth partition from
Section~\ref{sec:audit} preserves the order \(64\).
On each of its four observation regions, the shell and
within-shell bounds are therefore
\(C_{K,\rho,m}(n+1)2^{-n}\).
The constants may depend on the localization but not on \(n\).
For intersection neighborhoods they require no condition
\(|t-t'|\geq\delta>0\).

\section{Complete conversion and the real counterterm class}
\label{sec:conversion}

Let \(\rho_a,\rho_b\) be two admissible fixed mollifiers.
Each is used on both slots of both insertions, as prescribed on
the convention sheet.

\begin{theorem}[Complete regulator conversion]
\label{thm:conversion}
The difference has the strong limit
\begin{equation}
 \boxed{
 \begin{gathered}
 L_{ab}
 =\lim_{\eps\downarrow0}(T_\eps^a-T_\eps^b)
 =c_{ab}r,\\
 c_{ab}
 =\int_0^\infty
       \frac{\vartheta_a(u)^4-\vartheta_b(u)^4}{u}\dd u
 \in\R.
 \end{gathered}
 }
 \label{eq:conversion}
\end{equation}
For every common scale \(\mu>0\),
\[
                         u_\mu^a-u_\mu^b=L_{ab}.
\]
This is the complete conversion distribution: no further
normal derivative, nonconstant weight on \(\Rel\), or additional
term supported on \(\Rel\cap\Delta_{13}\) occurs.
\end{theorem}

\begin{proof}
The distribution \(U_m\) in \eqref{eq:finitepart} is independent
of the profile. Hence
\[
 L_{ab}=(\kappa_{\rho_a}-\kappa_{\rho_b})r.
\]
Combining the two integrals in \eqref{eq:kappa} gives the
convergent integral in \eqref{eq:conversion}.
It is real because both profiles are real.

Applying \eqref{eq:finalrate} to both regulators proves strong
convergence of their complete difference on all joint tests.
This leaves no distribution on a smaller observation stratum
undetected.
\end{proof}

\subsection{Perpendicular and reality checks}

Choose a real compact test \(\psi\) with \(\la r,\psi\ra=1\).
Such a test exists since \(r\) is a nonzero positive distribution.
Then
\[
 \la L_{ab},\psi\ra=c_{ab}\in\R,\qquad
 L_{ab}^{\perp}
 =L_{ab}-\la L_{ab},\psi\ra r=0.
\]
Thus both conditions for the prescribed real class
\(\cC=\R r\) are checked.
The weighted-pushforward representation is completely identified:
\[
                         L_{ab}=c_{ab}c_0\Psi_*\nu.
\]
Its weight is constant, as a consequence of the full conversion
limit rather than only the leading logarithmic coefficient.

\subsection{Adjoint symmetry}

Let \(\tau(x,y,x')=(x',y,x)\), and define
\[
 \la U^\dagger,\Phi\ra
 =
 \overline{\la U,\overline{\Phi\circ\tau}\ra}.
\]
The real multiplier in \eqref{eq:Deps} implies
\[
                         \overline{D_\eps(z)}=D_\eps(-z).
\]
Together with \eqref{eq:Texact}, this gives
\[
                         (T_\eps^\rho)^\dagger=T_\eps^\rho.
\]
The residue is real and invariant under \(\tau\), by interchanging
\(t,t'\) in \eqref{eq:rjoint}. Strong passage to the limit gives
\[
 r^\dagger=r,\qquad
 (u_\mu^\rho)^\dagger=u_\mu^\rho,\qquad
 L_{ab}^\dagger=L_{ab}.
\]
These symmetries have been proved directly; they are not used as
an assumption replacing the reality calculation.

\subsection{Reference profile, dilation, and cocycle}

Fix one admissible reference mollifier \(\rho_*\). Then
\[
 u_\mu^\rho
 =
 u_\mu^{\rho_*}
 +\left(
   \int_0^\infty
   \frac{\vartheta_\rho(u)^4-\vartheta_*(u)^4}{u}\dd u
  \right)r.
\]
In particular,
\[
 L_{ab}+L_{bc}=L_{ac},\qquad
 L_{aa}=0,\qquad L_{ba}=-L_{ab}.
\]

If
\[
 \rho_b(z)=\lambda^{-4}\rho_a(z/\lambda),\qquad \lambda>0,
\]
then the original insertion regularizations give the exact
finite-cutoff identity
\[
 T_\eps^b=T_{\lambda\eps}^a,\qquad
 \vartheta_b(u)=\vartheta_a(\lambda u).
\]
Changing variables in \eqref{eq:kappa} yields
\[
 \kappa_{\rho_b}=\kappa_{\rho_a}-\log\lambda,\qquad
 u_\mu^a-u_\mu^b=(\log\lambda)r.
\]
The subtraction-scale dependence is likewise
\[
 u_{\mu_2}^\rho-u_{\mu_1}^\rho
 =\log(\mu_2/\mu_1)\,r.
\]

\section{Closed numerical statement}
\label{sec:closed}

\begin{theorem}[Finite part for the prescribed VII-C family]
\label{thm:main}
Fix the convention-sheet data, a mass \(m\geq0\), and the class of
real, nonnegative, even, normalized, Euclidean-radial compact
mollifiers. Use one fixed mollifier on both slots of each
insertion, and keep the middle observable unsmeared.

Then:
\begin{enumerate}[label=\textup{(\roman*)}]
\item
For each \(\eps>0\), the numerical family is the smooth
distribution \eqref{eq:Texact}, with the prescribed product-test
pairings.

\item
The leading reflected coefficient, including its density, is
exactly
\[
                              r=\frac1{64(2\pi)^3}\Psi_*\nu.
\]
It is obtained from the exact branch integral by the normalization
ledger, not imposed on that integral.

\item
For every \(\mu>0\), the prescribed subtracted family has a strong
limit in \(\D'(M_+^3)\), with the fixed-order estimate
\eqref{eq:finalrate}. This includes the full reflected-diagonal
intersection.

\item
The regulator-independent ultraviolet decomposition is
\eqref{eq:UV}. Its radial and regulator-placement errors satisfy
\[
\begin{aligned}
 \la A(r)-r,\Phi\ra
 &=O_{K,m}(r^{-1})p_{K,64}(\Phi),\\
 \la E_\eps^\rho,\Phi\ra
 &=O_{K,\rho,m}\!\bigl(\eps(1+|\log\eps|)\bigr)p_{K,64}(\Phi).
\end{aligned}
\]

\item
For any two admissible mollifiers, the complete conversion is
\eqref{eq:conversion}. In particular,
\(L_{ab}^{\perp}=0\) and \(c_{ab}\in\R\).
The fixed class \(\cC=\R r\) suffices.
\end{enumerate}
These conclusions concern only the numerical prescription defined
in this paper.
\end{theorem}

\begin{proof}
Proposition~\ref{prop:positive} proves (i).
Proposition~\ref{prop:ledger} proves (ii) as an equality on arbitrary
joint tests.
Lemma~\ref{lem:S}, Lemma A, and
Lemma~\ref{lem:replacement} give the ultraviolet decomposition
and its uniform error estimates.
Theorem~\ref{thm:finitepart} proves (iii), while
Section~\ref{sec:audit} records the observation and momentum
coverage. Theorem~\ref{thm:conversion} proves (v).
\end{proof}

\begin{thebibliography}{9}

\bibitem{VIIC}
\emph{Landmark VII-C:
A mixed-Wick obstruction to a normal-topology boundary extension
for sharp Cauchy harmonic insertions}.
Manuscript source
\texttt{085T\_Landmark\_VIIC\_v1.tex}.
The source labels used here include
\texttt{eq:B}, \texttt{eq:C}, \texttt{eq:Hhh},
\texttt{eq:Hscalar}, \texttt{eq:positive}, \texttt{eq:K},
\texttt{eq:Beps}, \texttt{eq:reconstruction},
\texttt{lem:spectral}, \texttt{thm:log},
\texttt{eq:residue}, \texttt{lem:separated},
and \texttt{thm:coefficient}.

\end{thebibliography}

\end{document}