Skip to content
Shadow Theory

Appendix B 4 October 2026

Preparation estimates and transfer to the enlarged Hamiltonian

Reading position 33 of 37

B Preparation estimates and transfer to the enlarged Hamiltonian

All momenta in this appendix have the constant Galilean carrier removed. A superscript gg additionally removes the bounded preparation phase θ\theta. Thus an ordinary clock momentum here means pS−Mvp_S-Mv in the laboratory frame. We give the finite coefficient construction as well as its numerical enclosures; the latter are consequences of the construction, not extra hypotheses.

B.1 Coefficient conventions and the Gaussian residual

Write T=0.104T=0.104, T0=0.03T_0=0.03, Li=10−4L_i=10^{-4}, Lf=10−2L_f=10^{-2}, M=10M=10, v=1v=1, K=100K=100, and λ=252T0/Li2\lambda=2^{52}T_0/L_i^2, so that m=ℏλm=\hbar\lambda. Lengths and times in this appendix are in metres and seconds. Let bpb_p denote the preparation dilation, to distinguish it from the writer displacement. Set a=bp′/bpa=b_p'/b_p, f∗=0.201f_*=0.201, F∗=0.080802F_*=0.080802, g1=0.0005g_1=0.0005, g2=0.00025g_2=0.00025 and γ∗=141\gamma_*=141. Here F(r)=2∫0rfF(r)=2\int_0^r f is the bounded preparation profile; the subscript on F∗F_* in this subsection is unrelated to the writer force. Here b0(s)=1+99B9(s/0.1)b_0(s)=1+99B_9(s/0.1) has its constant extensions, and bpb_p is its convolution with the fixed 10−510^{-5}-scale mollifier. The explicit clamped ninth-degree polynomial and its fixed positive mollifier give the derivative bounds

B1=99(315/128)/0.1,B2=99(2520/64)/0.12,B3=99(7560/16+5040/64)/0.13,B4=99(1360800)/0.14.\begin{align}B_1&=99(315/128)/0.1,& B_2&=99(2520/64)/0.1^2,\tag{14}\\ B_3&=99(7560/16+5040/64)/0.1^3,& B_4&=99(1360800)/0.1^4. \tag{15}\end{align}

The sharper relative bounds bp′≤140bpb_p'\le140b_p and ∣bp′′∣≤24000bp|b_p''|\le24000b_p follow from polynomial positivity before convolution. One finite verification is to express each of 140b0−b0′140b_0-b_0', 24000b0−b0′′24000b_0-b_0'', 24000b0+b0′′24000b_0+b_0'' in Bernstein form on dyadic subintervals of [0,1][0,1]; bisect a subinterval if any coefficient is negative. All leaves are nonnegative by depth five. For reproducibility, if P(x)=∑j=0npjxjP(x)=\sum_{j=0}^n p_jx^j, its Bernstein coefficient of index kk on [u,u+w][u,u+w] is

∑j=0k(kj)(nj)∑l=jnpl(lj)ul−jwj. \sum_{j=0}^k\frac{\binom{k}{j}}{\binom{n}{j}} \sum_{l=j}^n p_l\binom{l}{j}u^{l-j}w^j.

Positive convolution preserves these inequalities. Put

a0=140,a1=24000+a02=43600,a2=B3+3a0(24000)+2a03=70141750,a3=B4+4a0B3+3(24000)2+12a02(24000)+6a04=1387431060000.\begin{align}a_0&=140,\tag{16}\\ a_1&=24000+a_0^2=43600,\tag{17}\\ a_2&=B_3+3a_0(24000)+2a_0^3=70141750,\tag{18}\\ a_3&=B_4+4a_0B_3+3(24000)^2+12a_0^2(24000)+6a_0^4 =1387431060000. \tag{19}\end{align}

These bound ∣a(j)∣|a^{(j)}|. Let B5B_5 be the largest absolute Bernstein coefficient on [0,1][0,1] of d5b0/ds5d^5b_0/ds^5 and set

a4=B5+5a0B4+10(24000)B3+20a02B3+30a0(24000)2+60a03(24000)+24a05. a_4=B_5+5a_0B_4+10(24000)B_3+20a_0^2B_3 +30a_0(24000)^2+60a_0^3(24000)+24a_0^5.

The normalized clock envelope has derivative norms bounded by

p0=1,p1=q16/7 π/(2σ),p2=q512/35 [π/(2σ)]2,p3=q1024/7 [π/(2σ)]3,p4=q(2π/σ)4,\begin{aligned}p_0&=1,&p_1&=q\sqrt{16/7}\,\pi/(2\sigma),& p_2&=q\sqrt{512/35}\,[\pi/(2\sigma)]^2,\\ p_3&=q\sqrt{1024/7}\,[\pi/(2\sigma)]^3,& p_4&=q(2\pi/\sigma)^4, \end{aligned}

where σ=0.001\sigma=0.001, q=1.000001q=1.000001. The convolution scale is 10−1010^{-10} and q(1−10−10p1/q)>1q(1-10^{-10}p_1/q)>1, which proves the normalization allowance. These are weak-derivative estimates for the compact envelope, not a band-limit assertion.

Let A=2rexp⁡[−r2/(2L2)]/(π1/4L3/2)A=2r\exp[-r^2/(2L^2)]/(\pi^{1/4}L^{3/2}), L=Libp(S)L=L_ib_p(S), and G=ϕ(S−Sc−vt)AG=\phi(S-S_c-vt)A. Set z=(r/L)2z=(r/L)^2, Mj=(2j+1)!!/2jM_j=(2j+1)!!/2^j and νj=M2j\nu_j=\sqrt{M_{2j}}. Thus ∥zjA∥=νj\|z^jA\|=\nu_j. For a polynomial P=∑cjzjP=\sum c_jz^j define N(P)=∑∣cj∣νj\mathcal N(P)=\sum|c_j|\nu_j. This coefficient norm is deliberately subadditive; signs of independently bounded derivatives cannot cancel. Use

H1=z−32,H2=z2−5z+94,H3=z3−212z2+794z−278. H_1=z-\tfrac32,\quad H_2=z^2-5z+\tfrac94,\quad H_3=z^3-\tfrac{21}2z^2+\tfrac{79}4z-\tfrac{27}8.

Put k=B1k=B_1, k1=B2+B12k_1=B_2+B_1^2, k2=B3+3B1B2+2B13k_2=B_3+3B_1B_2+2B_1^3, and

C=λLi2(100B2+B12)/2,Ck=λLi2(B1B2+B13)/2,C2=λLi2(100B3+3B1B2+2B13)/2,C3=λLi2(100B4+4B1B3+3B22+12B12B2+6B14)/2.\begin{aligned}C&=\lambda L_i^2(100B_2+B_1^2)/2,\\ C_k&=\lambda L_i^2(B_1B_2+B_1^3)/2,\\ C_2&=\lambda L_i^2(100B_3+3B_1B_2+2B_1^3)/2,\\ C_3&=\lambda L_i^2(100B_4+4B_1B_3+3B_2^2 +12B_1^2B_2+6B_1^4)/2. \end{aligned}

The three nonnegative coefficient bounds

r0=p2+3kp1+32k1+94k2,r1=2kp1+k1+5k2+2Cp1+3Ck+C2,r2=k2+2Ck\begin{aligned}r_0&=p_2+3kp_1+\tfrac32 k_1+\tfrac94k^2,\\ r_1&=2kp_1+k_1+5k^2+2Cp_1+3C_k+C_2,\\ r_2&=k^2+2C_k \end{aligned}

give the residual norm coefficient R0=∑j=02rjνjR_0=\sum_{j=0}^2r_j\nu_j. For its radial derivative set

Rr=∑j=02rj[2jM2j−1+M2j+1], R_r=\sum_{j=0}^2r_j\bigl[2j\sqrt{M_{2j-1}}+\sqrt{M_{2j+1}}\bigr],

where the first summand is zero at j=0j=0. This follows equally by the unitary radial reduction of the three-dimensional Gaussian gradient; its integration-by-parts identity includes the reduced amplitude's factor rr. For the clock derivative set

U1=kN(H1),U2=k1N(H1)+k2N(H2),U3=k2N(H1)+3kk1N(H2)+k3N(H3),RS=p3+3p2U1+3p1U2+U3+2p2Cν1+4p1kCN(zH1)+3p1C2ν1+2C[k1N(zH1)+k2N(zH2)]+3kC2N(zH1)+C3ν1.\begin{aligned}U_1&=k\mathcal N(H_1),\\ U_2&=k_1\mathcal N(H_1)+k^2\mathcal N(H_2),\\ U_3&=k_2\mathcal N(H_1)+3kk_1\mathcal N(H_2)+k^3\mathcal N(H_3),\\ R_S&=p_3+3p_2U_1+3p_1U_2+U_3 +2p_2C\nu_1+4p_1kC\mathcal N(zH_1)+3p_1C_2\nu_1\\ &\quad+2C[k_1\mathcal N(zH_1)+k^2\mathcal N(zH_2)] +3kC_2\mathcal N(zH_1)+C_3\nu_1. \end{aligned}

These are obtained by differentiating the complete gauged residual

Rgℏ=ℏ2M[GSS+2iθSGS+iθSSG]+iaϕ[(f−r)Ar+(f′−1)A/2]+(Q−Qc)G/ℏ.\frac{R^g}{\hbar}=\frac{\hbar}{2M} [G_{SS}+2i\theta_SG_S+i\theta_{SS}G] +ia\phi[(f-r)A_r+(f'-1)A/2]+(Q-Q_c)G/\hbar. (20)

In particular, no derivative of a restored phase is added to this gauged residual.

Here are all coefficients needed to include its tails and the weak trap. Write ωi=(λLi2)−1\omega_i=(\lambda L_i^2)^{-1}, q2=mωi2/2q_2=m\omega_i^2/2, q0=3ℏωi/2q_0=3\hbar\omega_i/2, q2S=4a0q2q_{2S}=4a_0q_2, q0S=2a0q0q_{0S}=2a_0q_0, rc=0.201r_c=0.201, c1=5000c_1=5000, and ϵ20=1012/7100\epsilon_{20}=10^{12}/7^{100}. Set

Q0=q2rc2+q0,Qr=2q2rc+c1Q0,QS=q2Src2+q0S+c1Q0,Qt=q2Lf2+q0,QtS=q2SLf2+q0S.\begin{aligned}Q_0&=q_2r_c^2+q_0,&Q_r&=2q_2r_c+c_1Q_0,\\ Q_S&=q_{2S}r_c^2+q_{0S}+c_1Q_0,& Q_t&=q_2L_f^2+q_0,& Q_{tS}&=q_{2S}L_f^2+q_{0S}. \end{aligned}

The complete residual rates are

d=ℏR0/(2M)+(λωi2Lf2/2+3ωi/2+3a0)ϵ20,er=ℏ2Rr/(2MLi)+[c1Qt+2q2Lf+Qt/Li+ℏa0(3/0.2+7/Li+2000)]ϵ20,eS=ℏ2RS/(2M)+[QtS+Qt(p1+5a0/2)+ℏ{3(a1+a0p1)+19a02/2}]ϵ20.\begin{align}d&=\hbar R_0/(2M) +(\lambda\omega_i^2L_f^2/2+3\omega_i/2+3a_0)\epsilon_{20},\tag{21}\\ e_r&=\hbar^2R_r/(2ML_i) +[c_1Q_t+2q_2L_f+Q_t/L_i +\hbar a_0(3/0.2+7/L_i+2000)]\epsilon_{20},\tag{22}\\ e_S&=\hbar^2R_S/(2M) +[Q_{tS}+Q_t(p_1+5a_0/2) +\hbar\{3(a_1+a_0p_1)+19a_0^2/2\}]\epsilon_{20}. \tag{23}\end{align}

The Gaussian tail recurrence ∫20∞yne−y2dy=20n−1e−400/2+(n−1)∫20∞yn−2e−y2dy/2\int_{20}^\infty y^n e^{-y^2}dy =20^{n-1}e^{-400}/2+(n-1)\int_{20}^\infty y^{n-2}e^{-y^2}dy/2 proves the common tail norm used here; e2>7e^2>7 makes its stated bound rational.

B.2 First moments and preparation current

In the gauge the velocity coefficients are ur=afu_r=af and uS=v+δu_S=v+\delta, δ=ma′F/(2M)\delta=ma'F/(2M). Put

dr=ma1f∗/M,dS=ma2F∗/(2M),urS=a1f∗,zr=Qr+ℏ2(a0 4000+ma2f∗/M),zS=QS+ℏ2(a1+ma3F∗/(2M)).\begin{aligned}d_r&=ma_1f_*/M,&d_S&=ma_2F_*/(2M),& u_{rS}&=a_1f_*,\\ z_r&=Q_r+\tfrac\hbar2(a_0\,4000+ma_2f_*/M),\\ z_S&=Q_S+\tfrac\hbar2(a_1+ma_3F_*/(2M)). \end{aligned}

If W0=ℏ235/T0W_0=\hbar2^{35}/T_0 and W1=ℏ(252/256+5⋅235)/(T0Li)W_1=\hbar(2^{52}/256+5\cdot2^{35})/(T_0L_i), the remote activation forces are bounded by Fr=W1+QrF_r=W_1+Q_r, FS=25000(W0+Q0)+QSF_S=25000(W_0+Q_0)+Q_S. Let Pj=∥pj(Ψg−G)∥P_j=\|p_j(\Psi^g-G)\| and D=∥Ψg−G∥D=\|\Psi^g-G\|. Differentiating the symmetrized drift gives

D(t)≤dt,Pr′≤a0Pr+drPS+Frn+er+zrdt,PS′≤urSPr+dSPS+FSn+eS+zSdt,n′≤g0+PS/(Mg1),g0=ℏ[p1+a03/2]/(Mg1).\begin{align}D(t)&\le dt,\tag{24}\\ P_r'&\le a_0P_r+d_rP_S+F_rn+e_r+z_rdt,\tag{25}\\ P_S'&\le u_{rS}P_r+d_SP_S+F_Sn+e_S+z_Sdt,\tag{26}\\ n'&\le g_0+P_S/(Mg_1),\tag{27}\\ g_0&=\hbar[p_1+a_0\sqrt{3/2}]/(Mg_1). \tag{28}\end{align}

The collar defining nn translates at v+∥δ∥∞v+\|\delta\|_\infty on the right and v−∥δ∥∞v-\|\delta\|_\infty on the left; its convective contribution is nonpositive. Both collars vanish on the initial stock. Their final full-one edges are respectively below 0.10350000012620.1035000001262 and above 0.10049999987380.1004999998738, so that they cover the remote activation and the left phase support respectively. Thus the force term FjnF_jn is valid throughout preparation, including the exact wave's tails.

For the sharper clock estimate take α=10−12\alpha=10^{-12} and

s=er+eS/K+αg0,c=zr+zS/K,E∗=(11/4)15. s=e_r+e_S/K+\alpha g_0,\qquad c=z_r+z_S/K,\qquad E_*=(11/4)^{15}.

The three scaled row sums a0+Kdr+Fr/αa_0+Kd_r+F_r/\alpha, urS/K+dS+FS/(Kα)u_{rS}/K+d_S+F_S/(K\alpha) and αK/(Mg1)\alpha K/(Mg_1) are less than 141141. A positive supersolution gives

B(t)=e141t(st+cdt2/2),Pr(t)≤B(t),PS(t)≤KB(t),BT=E∗(sT+cdT2/2),BI=E∗(sT2/2+cdT3/6),n1T=g0T+KBI/(Mg1).\begin{align}B(t)&=e^{141t}(st+cdt^2/2),& P_r(t)&\le B(t),&P_S(t)&\le KB(t),\tag{29}\\ B_T&=E_*(sT+cdT^2/2),& B_I&=E_*(sT^2/2+cdT^3/6),\tag{30}\\ n_{1T}&=g_0T+KB_I/(Mg_1). \tag{31}\end{align}

Using e141T<E∗e^{141T}<E_* is legitimate since 141T<15141T<15 and e<11/4e<11/4. These terminal majorants are increasing and apply at every earlier time.

For the smaller radial estimate and the preparation current retain also the mass-weighted system. Put ℓ=10−6\ell=10^{-6},

e0=er/m+eS/M,c0=zr/m+zS/M,s0=e0+ℓg0,ZT=E∗(s0T+c0dT2/2),ZI=E∗(s0T2/2+c0dT3/6),nw=g0T+ZI/(Mg1).\begin{aligned}e_0&=e_r/\sqrt m+e_S/\sqrt M,& c_0&=z_r/\sqrt m+z_S/\sqrt M,\\ s_0&=e_0+\ell g_0,& Z_T&=E_*(s_0T+c_0dT^2/2),\\ Z_I&=E_*(s_0T^2/2+c_0dT^3/6),& n_w&=g_0T+Z_I/(\sqrt M g_1). \end{aligned}

Indeed the weighted momentum plus ℓn\ell n has growth at most 141141: the drift row bound is a0+m/Ma1f∗+dSa_0+\sqrt{m/M}a_1f_*+d_S, increased by ℓ/(Mg1)\ell/(\sqrt M g_1), and (Fr/m+FS/M)/ℓ<141(F_r/\sqrt m+F_S/\sqrt M)/\ell<141. Consequently the ordinary radial error at the handoff is bounded by

Pr,o=m ZT+ma0f∗nw,PS,o=KBT+ma1F∗n1T/2,Do=dT.P_{r,o}=\sqrt m\,Z_T+ma_0f_*n_w, \qquad P_{S,o}=KB_T+ma_1F_*n_{1T}/2, \qquad D_o=dT. (32)

The phase terms here are localized to the reflected collar; no phase is set to zero on an exact wave tail.

For the Gaussian rank let kr=6/(7Li)k_r=6/(7L_i), kS=2a0k_S=2a_0, cr=kr3/2/Lic_r=k_r\sqrt{3/2}/L_i, cS=kS(p1+a03/2)c_S=k_S(p_1+a_0\sqrt{3/2}). The exact material-rank identity and the bilinear current expansion give

Io=(cr/λ+ℏcS/M)dT2/2+(kr/m+kS/M)ZI,β=mLi2(B1B2+a0B12)T/M+Ta0 32000/7200.\begin{align}I_o&=(c_r/\lambda+\hbar c_S/M)dT^2/2 +(k_r/\sqrt m+k_S/\sqrt M)Z_I,\tag{33}\\ \beta&=mL_i^2(B_1B_2+a_0B_1^2)T/M +Ta_0\,32000/7^{200}. \tag{34}\end{align}

The first term of β\beta bounds KSδK_S\delta pointwise; the second bounds the nonquadratic radial tail. The remaining current has absolute value before integration over any live coordinate. This proves Io<1.974548968⋅10−19I_o<1.974548968\cdot10^{-19} and β<2.638660⋅10−13\beta<2.638660\cdot10^{-13}.

B.3 Second spatial moments and the second collar

Here is an explicit construction of every second-order residual coefficient. Continue HjH_j by H0=1H_0=1, Hj+1=(z−3/2)Hj−2zHj′H_{j+1}=(z-3/2)H_j-2zH_j', and take absolute coefficients before adding independently bounded terms. The positive Gaussian clock polynomials are

A0=1,A1=a0∣H1∣,A2=a1∣H1∣+a02∣H2∣,A3=a2∣H1∣+3a0a1∣H2∣+a03∣H3∣,A4=a3∣H1∣+(4a0a2+3a12)∣H2∣+6a02a1∣H3∣+a04∣H4∣.\begin{aligned}A_0&=1,& A_1&=a_0|H_1|,\\ A_2&=a_1|H_1|+a_0^2|H_2|,& A_3&=a_2|H_1|+3a_0a_1|H_2|+a_0^3|H_3|,\\ A_4&=a_3|H_1|+(4a_0a_2+3a_1^2)|H_2| +6a_0^2a_1|H_3|+a_0^4|H_4|. \end{aligned}

Here ∣P∣|P| means coefficientwise absolute value. To evaluate radial norms, form Pj,r(y)=(∂y−y)ry2j+1P_{j,r}(y)=(\partial_y-y)^r y^{2j+1} and set

Rr(j;s)=Li−r∑e∣[ye]Pj,r∣Me−1+s,M−1=2. \mathcal R_r(j;s)=L_i^{-r} \sum_e |[y^e]P_{j,r}|\sqrt{M_{e-1+s}},\qquad M_{-1}=2.

Only indices at least −1-1 occur. For a factor F/L2F/L^2 and r≤2r\le2 replace this by

RrF(j)=∑l=0r(rl)clLi−lRr−l(j;2−l),(c0,c1,c2)=(1,2,2). \mathcal R_r^F(j)=\sum_{l=0}^r\binom rl c_lL_i^{-l} \mathcal R_{r-l}(j;2-l), \qquad(c_0,c_1,c_2)=(1,2,2).

This is the product rule using F≤r2F\le r^2, ∣F′∣≤2r|F'|\le2r, ∣F′′∣≤2|F''|\le2. Define

Jn,r=∑j=0n(nj)pn−j∑l[zl]Aj Rr(l;0), J_{n,r}=\sum_{j=0}^n\binom nj p_{n-j} \sum_l[z^l]A_j\,\mathcal R_r(l;0),

and Jn,rFJ_{n,r}^F by replacing Rr\mathcal R_r with RrF\mathcal R_r^F. With Cj∗=λLf2aj/2C_j^*=\lambda L_f^2a_j/2 for 1≤j≤41\le j\le4, the sources are

err=ℏ32M(J2,2+2C1∗J1,2F+C2∗J0,2F)+10−100,erS=ℏ32M(J3,1+2C1∗J2,1F+3C2∗J1,1F+C3∗J0,1F)+10−100,eSS=ℏ32M(J4,0+2C1∗J3,0F+5C2∗J2,0F+4C3∗J1,0F+C4∗J0,0F)+10−100.\begin{aligned}e_{rr}&=\frac{\hbar^3}{2M} (J_{2,2}+2C_1^*J_{1,2}^F+C_2^*J_{0,2}^F)+10^{-100},\\ e_{rS}&=\frac{\hbar^3}{2M} (J_{3,1}+2C_1^*J_{2,1}^F+3C_2^*J_{1,1}^F+C_3^*J_{0,1}^F) +10^{-100},\\ e_{SS}&=\frac{\hbar^3}{2M} (J_{4,0}+2C_1^*J_{3,0}^F+5C_2^*J_{2,0}^F +4C_3^*J_{1,0}^F+C_4^*J_{0,0}^F)+10^{-100}. \end{aligned}

For example the last formula differentiates GSS+2iθSGS+iθSSGG_{SS}+2i\theta_SG_S+i\theta_{SS}G twice and retains coefficients 2,5,4,12,5,4,1. The nonquadratic and trap tails after two derivatives have Gaussian degree at most seven and total frequency coefficient below 103010^{30}; their momentum contribution is less than ℏ21030ϵ20<10−100\hbar^2 10^{30}\epsilon_{20}<10^{-100} in each row. This can also be checked by the general product-rule construction below.

For clarity we spell out the second-force coefficients. Set c2=1.28⋅108c_2=1.28\cdot10^8, qt1=q2Src2+q0Sq_{t1}=q_{2S}r_c^2+q_{0S}, qt2=(4a1+16a02)q2rc2+(2a1+4a02)q0q_{t2}=(4a_1+16a_0^2)q_2r_c^2+(2a_1+4a_0^2)q_0. Then

Qrr=2q2+4c1q2rc+c2Q0,QrS=c1Qr+2q2Src+c1qt1,QSS=c2Q0+2c1qt1+qt2.\begin{aligned}Q_{rr}&=2q_2+4c_1q_2r_c+c_2Q_0,\\ Q_{rS}&=c_1Q_r+2q_{2S}r_c+c_1q_{t1},\\ Q_{SS}&=c_2Q_0+2c_1q_{t1}+q_{t2}. \end{aligned}

Put γ1=25000\gamma_1=25000, γ2=(2520/64)/(0.0001)2\gamma_2=(2520/64)/(0.0001)^2, and

W2=ℏT0Li2[252(1+(15/16)217)+5⋅252/128+105⋅235],Frr=W2+Qrr,FrS=γ1(W1+Qr)+QrS,FSS=γ2(W0+Q0)+2γ1QS+QSS.\begin{aligned}W_2&=\frac{\hbar}{T_0L_i^2} [2^{52}(1+(15/16)2^{17})+5\cdot2^{52}/128+10^5\cdot2^{35}],\\ F_{rr}&=W_2+Q_{rr},\\ F_{rS}&=\gamma_1(W_1+Q_r)+Q_{rS},\\ F_{SS}&=\gamma_2(W_0+Q_0)+2\gamma_1Q_S+Q_{SS}. \end{aligned}

Let vjrv_{jr} and vjSv_{jS} denote the component derivatives of (af,ma′F/(2M))(af,ma'F/(2M)). The second-derivative ceilings used below are

rrrSSSura0 4000a1a2f∗δma1/Mma2f∗/Mma3F∗/(2M)div⁡ua0 107+ma2/Ma1 4000+ma3f∗/Ma2+ma4F∗/(2M) \begin{array}{c|ccc} &rr&rS&SS\\\hline u_r&a_0\,4000&a_1&a_2f_*\\ \delta&ma_1/M&ma_2f_*/M&ma_3F_*/(2M)\\ \operatorname{div}u&a_0\,10^7+ma_2/M& a_1\,4000+ma_3f_*/M&a_2+ma_4F_*/(2M) \end{array}

Write drdiv=a0 4000+ma2f∗/Md^{\rm div}_r=a_0\,4000+ma_2f_*/M and dSdiv=a1+ma3F∗/(2M)d^{\rm div}_S=a_1+ma_3F_*/(2M) for the first divergence ceilings. If Ujk,Vjk,DjkU_{jk},V_{jk},D_{jk} denote the three rows of this table, all scalar coefficients of the Hessian comparison are

e2=err+erS/K+eSS/K2,cB=ℏ(Urr+KVrr+drdiv)+2Qr+{ℏ[UrS+KVrS+(Kdrdiv+dSdiv)/2]+KQr+QS}/K+{ℏ[USS+KVSS+KdSdiv]+2KQS}/K2,cD=ℏ22(Drr+DrS/K+DSS/K2)+ℏ(Qrr+QrS/K+QSS/K2),cn=ℏ(Frr+FrS/K+FSS/K2)+2Frℏ/(g2K)+4FSℏ/(g2K2),cA=2Fr+2FS/K.\begin{aligned}e_2&=e_{rr}+e_{rS}/K+e_{SS}/K^2,\\ c_B&=\hbar(U_{rr}+KV_{rr}+d^{\rm div}_r)+2Q_r\\ &\quad+\{\hbar[U_{rS}+KV_{rS}+(Kd^{\rm div}_r+d^{\rm div}_S)/2] +KQ_r+Q_S\}/K\\ &\quad+\{\hbar[U_{SS}+KV_{SS}+Kd^{\rm div}_S]+2KQ_S\}/K^2,\\ c_D&=\tfrac{\hbar^2}{2}(D_{rr}+D_{rS}/K+D_{SS}/K^2) +\hbar(Q_{rr}+Q_{rS}/K+Q_{SS}/K^2),\\ c_n&=\hbar(F_{rr}+F_{rS}/K+F_{SS}/K^2) +2F_r\hbar/(g_2K)+4F_S\hbar/(g_2K^2),\\ c_A&=2F_r+2F_S/K. \end{aligned}

These follow directly by commuting two momenta through {uj,pj}/2+Qc+VA\{u_j,p_j\}/2+Q_c+V_A. For example the mixed commutator contributes ℏ∑l∥ul,jk∥Pl\hbar\sum_l\|u_{l,jk}\|P_l and ℏ(∥∂jdiv⁡u∥Pk+∥∂kdiv⁡u∥Pj)/2\hbar(\|\partial_j\operatorname{div}u\|P_k+ \|\partial_k\operatorname{div}u\|P_j)/2; both occur above. The three homogeneous scaled row sums are below 282282.

With Z2=max⁡(∥pr2Eg∥,∥prpSEg∥/K,∥pS2Eg∥/K2)Z_2=\max(\|p_r^2E^g\|,\|p_rp_SE^g\|/K, \|p_S^2E^g\|/K^2), nested clock cutoffs give

Z2′≤282Z2+cAn1Z2+e2+cBB+cDdt+cnn1. Z_2'\le282Z_2+c_A\sqrt{n_1}\sqrt{Z_2} +e_2+c_BB+c_Ddt+c_nn_1.

Indeed a cutoff supported where the first collar is one satisfies ∥χprEg∥≤n1∥pr2Eg∥\|\chi p_rE^g\|\le\sqrt{n_1\|p_r^2E^g\|} and ∥χpSEg∥≤n1∥pS2Eg∥+2ℏn1/g2\|\chi p_SE^g\|\le\sqrt{n_1\|p_S^2E^g\|}+2\hbar n_1/g_2. These inequalities are integration by parts, not finite propagation. Let bn=K/(Mg1)b_n=K/(Mg_1), γ=141\gamma=141, and define

Jn=g0π2γ3/2+bns/2(γ/2)2+3πbncd/64(γ/2)5/2,Jf=e22γ+cB(s/γ2+cd/γ3)+cDd(2γ)2+cn{g0/(2γ)2+bn(s/γ3+cd/γ4)},Z2T=e2γT(cAJn/2+Jf)2,HSS=K2Z2T.\begin{aligned}J_n&=\frac{\sqrt{g_0\pi}}{2\gamma^{3/2}} +\frac{\sqrt{b_ns/2}}{(\gamma/2)^2} +\frac{3\sqrt{\pi b_ncd/6}}{4(\gamma/2)^{5/2}},\\ J_f&=\frac{e_2}{2\gamma} +c_B(s/\gamma^2+cd/\gamma^3) +\frac{c_Dd}{(2\gamma)^2}\\ &\quad+c_n\{g_0/(2\gamma)^2+b_n(s/\gamma^3+cd/\gamma^4)\},\\ Z_{2T}&=e^{2\gamma T}(c_AJ_n/2+\sqrt{J_f})^2, \qquad H_{SS}=K^2Z_{2T}. \end{aligned}

For the last display, extend positive integrals to infinity after using n1(t)≤g0t+bneγt(st2/2+cdt3/6)n_1(t)\le g_0t+b_ne^{\gamma t}(st^2/2+cdt^3/6). The comparison follows by setting y=e−γtZ2y=e^{-\gamma t}\sqrt{Z_2}: the sum of the separate positive supersolutions for y′≤A+B/(2y)y'\le A+B/(2y) is a supersolution, with zero handled by regularization. The factor eγTe^{\gamma T} is enclosed by its 90-term Taylor sum and the remaining positive geometric tail with ratio γT/91\gamma T/91; its upper enclosure is squared when evaluating e2γTe^{2\gamma T}.

The second collar amplitude and the physical phase restoration are

n2T=TMg2[n1THSS+2ℏDo/g1],LS=n2THSS+2ℏn1T/g2,Lr=n2THSS/K2,HSSo=HSS+2q1LS+(q12+ℏq2∗)n2T,q1=ma1F∗/2,q2∗=ma2F∗/2.\begin{align}n_{2T}&=\frac{T}{Mg_2} [\sqrt{n_{1T}H_{SS}}+2\hbar D_o/g_1],\tag{35}\\ L_S&=\sqrt{n_{2T}H_{SS}}+2\hbar n_{1T}/g_2, &L_r&=\sqrt{n_{2T}H_{SS}/K^2},\tag{36}\\ H_{SS}^{o}&=H_{SS}+2q_1L_S+(q_1^2+\hbar q_2^*)n_{2T}, &q_1&=ma_1F_*/2,\quad q_2^*=ma_2F_*/2. \tag{37}\end{align}

Both reflected second collars fit before the relevant remote support: their final edges are below 0.1037500001270.103750000127 and above 0.1002499998730.100249999873. These same increasing majorants hold for all preceding times. In particular the ordinary derivatives on the writer support are bounded by LS,LrL_S,L_r, because both the helper and the phase vanish there.

B.4 Fourth spatial closure and its complete coefficient ledger

Only ordinary spatial derivatives are used here. There is no assumption that the unmollified radial reference lies in the domain of a fourth power of its Hamiltonian. Use the commuting scaled momenta Dr=prD_r=p_r, DS=pS/KD_S=p_S/K, and Zj=max⁡a+b=j∥DraDSbEg∥Z_j=\max_{a+b=j}\|D_r^aD_S^bE^g\|.

The higher coefficient construction is finite. Define qj=bp(j)/bpq_j=b_p^{(j)}/b_p. Starting with P0(q)=q1P_0(q)=q_1, form

Pj+1=∑l=16(ql+1−q1ql)∂qlPj,0≤j<6. P_{j+1}=\sum_{l=1}^6(q_{l+1}-q_1q_l)\partial_{q_l}P_j, \qquad 0\le j<6.

Evaluate the absolute coefficients at

(qˉ1,…,qˉ7)=(140,24000,B3,B4,B5,4B5/10−5,300B5/10−10) (\bar q_1,\ldots,\bar q_7) =(140,24000,B_3,B_4,B_5,4B_5/10^{-5},300B_5/10^{-10})

to obtain aˉj≥∣a(j)∣\bar a_j\ge|a^{(j)}|. For j≤4j\le4 these agree with the bounds above. The last two entries use convolution of b0(5)b_0^{(5)} with the first and second derivative of the mollifier. Similarly set p5=4p4/10−10p_5=4p_4/10^{-10}, p6=300p4/10−20p_6=300p_4/10^{-20}. The normalized unit bump has derivative polynomials generated by

U0(x,u)=1,Uj+1=∂xUj+2xu2∂uUj−2xu2Uj,u=(1−x2)−1.U_0(x,u)=1,\qquad U_{j+1}=\partial_xU_j+2xu^2\partial_uU_j-2xu^2U_j, \qquad u=(1-x^2)^{-1}. (38)

Its normalization is greater than 1/41/4; hence its jjth derivative supremum is bounded by 4∑l,k∣[xluk]Uj∣ kk/(8/3)k4\sum_{l,k}|[x^lu^k]U_j|\,k^k/(8/3)^k, with the k=0k=0 factor defined as one. This proves the quoted mollifier bounds (the first derivative also uses unimodality), as well as all tent jets below.

Let the positive partial Bell coefficients be

B0,0=1,Bn,k=∑j=1n−k+1(n−1j−1)aˉj−1Bn−j,k−1. \mathcal B_{0,0}=1,\qquad \mathcal B_{n,k}=\sum_{j=1}^{n-k+1}\binom{n-1}{j-1} \bar a_{j-1}\mathcal B_{n-j,k-1}.

All unspecified entries are zero. For n=0n=0 put A0=1\mathcal A_0=1; otherwise An=∑k=1nBn,k∣Hk∣\mathcal A_n=\sum_{k=1}^n\mathcal B_{n,k}|H_k|. For ε=0,1\varepsilon=0,1, define the computable Gaussian norm majorant

Jn,r(ε)=∑j=0n(nj)pn−j∑l[zl]Aj Li−r∑e∣[ye](∂y−y)ry2l+1+2ε∣Me−1. \mathcal J_{n,r}^{(\varepsilon)}= \sum_{j=0}^n\binom nj p_{n-j} \sum_l[z^l]\mathcal A_j\,L_i^{-r} \sum_e\left|[y^e](\partial_y-y)^r y^{2l+1+2\varepsilon}\right| \sqrt{M_{e-1}}.

For the higher-core calculation one may use π≤22/7\pi\le22/7 in pjp_j and ℏ≤hˉ=1.055⋅10−34\hbar\le\bar h=1.055\cdot10^{-34}. The fourth physical source coefficient is bounded by

s4=hˉ52M∑n=04K−n[Jn+2,4−n(0)+∑j=1n+2cn,jλLf2aˉj2Jn+2−j,4−n(1)],cn,j=2(nj−1)+(nj−2),\begin{align}s_4&=\frac{\bar h^5}{2M}\sum_{n=0}^4K^{-n}\Bigl[ \mathcal J_{n+2,4-n}^{(0)} +\sum_{j=1}^{n+2}c_{n,j}\frac{\lambda L_f^2\bar a_j}{2} \mathcal J_{n+2-j,4-n}^{(1)}\Bigr],\tag{39}\\ c_{n,j}&=2\binom n{j-1}+\binom n{j-2}, \tag{40}\end{align}

where a binomial coefficient outside its range is zero. Direct rational arithmetic and outward square roots give s4<7.141⋅10−127s_4<7.141\cdot10^{-127}. The physical phase-square cancellation in the gauge is essential to this small coefficient.

We next give a complete way to check the nonquadratic tail bound, including higher derivatives of F−r2F-r^2. Let bjb_j denote the bump derivative suprema just constructed and use

(f0,f1,f2,f3,f4,f5)=(0.201,1,4000,107,2b2/0.0013,2b3/0.0014),F0=0.080802,Fj=2fj−1(1≤j≤6). \begin{gathered} (f_0,f_1,f_2,f_3,f_4,f_5) =(0.201,1,4000,10^7,2b_2/0.001^3,2b_3/0.001^4),\\ F_0=0.080802,\qquad F_j=2f_{j-1}\quad(1\le j\le6). \end{gathered}

The disjoint jumps of the unsmoothed tent give these bounds. The scaled velocity derivative sum at order jj is explicitly

Vj=max⁡0≤k≤jaˉkfj−kKk+max⁡0≤k≤jKhˉλaˉk+1Fj−k2MKk.V_j=\max_{0\le k\le j}\frac{\bar a_k f_{j-k}}{K^k} +\max_{0\le k\le j}\frac{K\bar h\lambda\bar a_{k+1}F_{j-k}} {2MK^k}. (41)

It gives V2<560001V_2<560001, V3<1.400001⋅109V_3<1.400001\cdot10^9, V4<3.669751⋅1013V_4<3.669751\cdot10^{13}, V5<2.786850⋅1018V_5<2.786850\cdot10^{18}.

For an entirely algebraic tail check replace every Gaussian norm in Jn,r(0)\mathcal J_{n,r}^{(0)} by the absolute sum of its polynomial coefficients; call the result Cn,r\mathcal C_{n,r}. Keep powers of yy as formal factors. The derivatives of F−r2F-r^2 have coefficient ceilings ΔF=(2,2,4,2f2,2f3)\Delta F=(2,2,4,2f_2,2f_3), and those of f−rf-r have Δf=(2,2,f2,f3,f4,f5)\Delta f=(2,2,f_2,f_3,f_4,f_5). The unit logistic cutoff has order-jj derivative bounded by

tj=8j!∑k=1j(j−1k−1)2k(j+k)j+k(8/3)j+k(1≤j≤4). t_j=8j!\sum_{k=1}^j\binom{j-1}{k-1}2^k \frac{(j+k)^{j+k}}{(8/3)^{j+k}}\quad(1\le j\le4).

To prove it on x≤1/2x\le1/2, put u=1/xu=1/x: the exponent is at most −u+2-u+2, its jjth derivative at most 2j!uj+12j!u^{j+1}, and apply the Bell product rule and max⁡une−u≤nn/(8/3)n\max u^ne^{-u}\le n^n/(8/3)^n. Symmetry handles the other half. Set t0=1t_0=1 and τj=tj/0.001j\tau_j=t_j/0.001^j. For p=2,4p=2,4 define Jp(0)=1J_p(0)=1, Jp(n)=∑k=1nBn,kpkJ_p(n)=\sum_{k=1}^n\mathcal B_{n,k}p^k. With A2=4⋅10−41A_2=4\cdot10^{-41} and A0=1.2⋅10−48A_0=1.2\cdot10^{-48}, put

q(n,r)=A2J4(n)(Lf2,2Lf,2)r+1r=0A0J2(n),q^(n,r)=q(n,r)+∑c=0n∑e=0r(nc)(re)τcτeq(n−c,r−e).\begin{aligned}q(n,r)&=A_2J_4(n)(L_f^2,2L_f,2)_r +\one_{r=0}A_0J_2(n),\\ \widehat q(n,r)&=q(n,r)+ \sum_{c=0}^n\sum_{e=0}^r\binom nc\binom re \tau_c\tau_e q(n-c,r-e). \end{aligned}

The three-entry sequence is zero outside r=0,1,2r=0,1,2. A bound for the sum of all fourth-derivative tail frequency coefficients is the following finite positive sum, with r=4−nr=4-n:

T=∑n=04∑a=0n∑b=0r(na)(rb)[1.425⋅10−122Maˉa+1ΔFbCn−a+1,r−b+1.425⋅10−124Maˉa+2ΔFbCn−a,r−b+aˉaΔfbCn−a,r−b+1+12aˉaΔfb+1Cn−a,r−b+1034q^(a,b)Cn−a,r−b]<1070.\begin{align}\mathcal T={}&\sum_{n=0}^4\sum_{a=0}^n\sum_{b=0}^{r} \binom na\binom rb\Bigl[ \frac{1.425\cdot10^{-12}}{2M}\bar a_{a+1}\Delta F_b \mathcal C_{n-a+1,r-b}\tag{42}\\ &+\frac{1.425\cdot10^{-12}}{4M}\bar a_{a+2}\Delta F_b \mathcal C_{n-a,r-b} +\bar a_a\Delta f_b\mathcal C_{n-a,r-b+1}\tag{43}\\ &+\tfrac12\bar a_a\Delta f_{b+1}\mathcal C_{n-a,r-b} +10^{34}\widehat q(a,b)\mathcal C_{n-a,r-b}\Bigr]<10^{70}. \tag{44}\end{align}

This is simply the product rule applied respectively to the two phase terms, transport tail, and Q−QcQ-Q_c in (20). Taking absolute coefficients before forming C\mathcal C prevents cancellation between independently bounded clock jets. As a rough independent size check, pj≤109jp_j\le10^{9j}, aˉj<103109j\bar a_j<10^3 10^{9j}, and Gaussian mixed coefficients through order six are below 1012109j10^{12}10^{9j}. The transport term has at most five derivative slots and fewer than 10510^5 product summands, giving 103+3+12+45+5=106810^{3+3+12+45+5}=10^{68}. In the phase term ℏθ\hbar\theta replaces the large phase coefficient by the physical mass: six derivative slots give at most 10−12+3+3+12+54+5=106510^{-12+3+3+12+54+5}=10^{65}. The pure GSSG_{SS} term is smaller. The explicit positive sum above also checks the weak potential term. No nonzero tail is discarded.

The degree is at most twelve. Integration by parts gives

∥1y≥20y12A∥2≤2⋅20251−25/800e−400<1040/7200. \|\one_{y\ge20}y^{12}A\|^2 \le\frac{2\cdot20^{25}}{1-25/800}e^{-400} <10^{40}/7^{200}.

Thus the additional fourth physical source is stail=hˉ41090/7100<10−130s_{\rm tail}=\bar h^4 10^{90}/7^{100}<10^{-130}. The same finite product rule for the quiet potential gives, for total order j≤4j\le4,

Qj≤10−40109j.Q_j\le10^{-40}10^{9j}. (45)

For an explicit verification use the preceding Jp(n)J_p(n) and cutoffs, replace (Lf2,2Lf,2)(L_f^2,2L_f,2) by (1,1,2)(1,1,2) on the cutoff support, and evaluate ∑a=0n∑b=0r(na)(rb)τaτb[A2J4(n−a)(1,1,2)r−b+1r=bA0J2(n−a)]\sum_{a=0}^n\sum_{b=0}^r\binom na\binom rb \tau_a\tau_b[A_2J_4(n-a)(1,1,2)_{r-b} +\one_{r=b}A_0J_2(n-a)] for n+r=jn+r=j.

For the remote potential define Wj=hˉ wj/(T0Lij)W_j=\bar h\,w_j/(T_0L_i^j), where

w0=235,w1=252/256+5⋅235,w2=252(1+(15/16)217)+5⋅252/128+105⋅235,w3=2112,w4=2160.\begin{aligned}w_0&=2^{35},\qquad w_1=2^{52}/256+5\cdot2^{35},\\ w_2&=2^{52}(1+(15/16)2^{17})+5\cdot2^{52}/128+10^5\cdot2^{35},\\ w_3&=2^{112},\qquad w_4=2^{160}. \end{aligned}

Set

(γ0,…,γ4)=(1,25000,2520/6410−8,7560/16+5040/6410−12,136080010−16),Rj=max⁡0≤k≤jK−k[γkWj−k+∑l=0k(kl)γlQj−l].\begin{aligned}(\gamma_0,\ldots,\gamma_4)&=\left(1,25000, \frac{2520/64}{10^{-8}}, \frac{7560/16+5040/64}{10^{-12}}, \frac{1360800}{10^{-16}}\right),\\ R_j&=\max_{0\le k\le j}K^{-k} \left[\gamma_kW_{j-k}+ \sum_{l=0}^k\binom kl\gamma_lQ_{j-l}\right]. \end{aligned}

These bounds follow by differentiating the displayed capped radial potential and activation polynomial; the last two deliberately loose bounds require only its bounded fourth weak derivative. All remote supports lie at S≥0.104S\ge0.104.

Add two smooth analytic cutoffs of width g=10−4g=10^{-4}, with final transition intervals (0.103751,0.103851)(0.103751,0.103851) and (0.103851,0.103951)(0.103851,0.103951). They fit between the second old collar and activation. Their first two derivative bounds are 5/g5/g and 128/g2128/g^2. Put ac=hˉ/(Kg)a_c=\bar h/(Kg). With nn the second-collar amplitude, integration by parts and interpolation give local derivative bounds

N1≤nZ2+10acn,N2≤n(Z4+18acZ2),N3≤N2Z4+10acN2.\begin{align}N_1&\le\sqrt{nZ_2}+10a_cn,\tag{46}\\ N_2&\le\sqrt n(\sqrt{Z_4}+18a_c\sqrt{Z_2}),\tag{47}\\ N_3&\le\sqrt{N_2Z_4}+10a_cN_2. \tag{48}\end{align}

For example the squared second-derivative norm is at most nZ4+20acnZ2Z4+306ac2nZ2nZ_4+20a_cn\sqrt{Z_2Z_4}+306a_c^2nZ_2; its square root is bounded as displayed since 306<182306<18^2. The global interpolation Z3≤Z2Z4Z_3\le\sqrt{Z_2Z_4} moves one commuting derivative between the two factors in its squared norm. Odd extension removes the radial boundary term.

The useful time envelopes, with Z1T=BTZ_{1T}=B_T, are

Z1(t)≤Z1Te−γ(T−t),Z2(t)≤Z2Te−2γ(T−t),n(t)≤Ant3/2eγt+Bnt2e3γt/2+Cnt2,An=K2Z2Tg0/(Mg2eγT),Bn=K2Z2Tn1T/(Mg2Te3γT/2),Cn=2hˉDT/(Mg2g1T),DT=3.006684⋅10−11.\begin{aligned}Z_1(t)&\le Z_{1T}e^{-\gamma(T-t)},& Z_2(t)&\le Z_{2T}e^{-2\gamma(T-t)},\\ n(t)&\le A_nt^{3/2}e^{\gamma t}+B_nt^2e^{3\gamma t/2}+C_nt^2,\\ A_n&=\sqrt{K^2Z_{2T}g_0}/(Mg_2e^{\gamma T}),\\ B_n&=\sqrt{K^2Z_{2T}n_{1T}}/(Mg_2Te^{3\gamma T/2}),\\ C_n&=2\bar h D_T/(Mg_2g_1T),\qquad D_T=3.006684\cdot10^{-11}. \end{aligned}

Here DTD_T is used only in the fourth-order majorant, not to replace the sharper handoff norm in the final probability accounting. Writing Z=Z4Z=Z_4, four commutations yield

Z′≤4γZ+α(t)Z3/4+β(t)Z1/2+f(t),α=4R1n1/4,β=4R1[18ac n1/4Z21/4+10acn1/2]+6hˉR2n1/2+(8hˉV2+4Q1)Z2,f=(720R1ac2+108hˉR2ac)nZ2+4hˉ2R3(nZ2+10acn)+hˉ3R4n+(7hˉ2V3+6hˉQ2)Z2+(3hˉ3V4+4hˉ2Q3)Z1+(hˉ4V5/2+hˉ3Q4)dt+s4+stail.\begin{align}Z'&\le4\gamma Z+\alpha(t)Z^{3/4}+\beta(t)Z^{1/2}+f(t),\tag{49}\\ \alpha&=4R_1n^{1/4},\tag{50}\\ \beta&=4R_1[\sqrt{18a_c}\,n^{1/4}Z_2^{1/4}+10a_cn^{1/2}] +6\bar hR_2n^{1/2}+(8\bar hV_2+4Q_1)\sqrt{Z_2},\tag{51}\\ f&=(720R_1a_c^2+108\bar hR_2a_c)\sqrt{nZ_2} +4\bar h^2R_3(\sqrt{nZ_2}+10a_cn)+\bar h^3R_4n\tag{52}\\ &\quad+(7\bar h^2V_3+6\bar hQ_2)Z_2 +(3\bar h^3V_4+4\bar h^2Q_3)Z_1 +(\bar h^4V_5/2+\bar h^3Q_4)dt+s_4+s_{\rm tail}. \tag{53}\end{align}

The factors 8,7,3,1/28,7,3,1/2 in the drift terms are the sum of the first-order drift and its divergence product-rule coefficients. Putting y=e−γtZ1/4y=e^{-\gamma t}Z^{1/4} and adding the separate positive supersolutions gives

Z4(T)1/4≤eγT[14∫0Tαe−γtdt+12∫0Tβe−2γtdt+(∫0Tfe−4γtdt)1/4].Z_4(T)^{1/4}\le e^{\gamma T} \left[\tfrac14\int_0^T\alpha e^{-\gamma t}dt +\sqrt{\tfrac12\int_0^T\beta e^{-2\gamma t}dt} +\left(\int_0^Tf e^{-4\gamma t}dt\right)^{1/4}\right]. (54)

This includes the nonlinear third-derivative coupling; it is not an exponential multiplying a terminal forcing value.

For completeness all integrals in this expression can be evaluated by one finite rule. For 0<p≤10<p\le1 use subadditivity on the three terms of n(t)n(t); combine with Z2(t)qZ_2(t)^q and the integrating factor, and extend each positive integral to infinity. If its time power is u≥0u\ge0 and its decay is ρ>0\rho>0, with j=⌊u⌋j=\lfloor u\rfloor, use

∫0∞tue−ρtdt≤j!(j+1)u−j/ρu+1. \int_0^\infty t^ue^{-\rho t}dt \le j!(j+1)^{u-j}/\rho^{u+1}.

This is Hölder interpolation of adjacent integer moments. All decay rates in (54) are positive. Integer roots are enclosed rationally; no numerical gamma function is needed. For a final far cutoff of width 0.010.01 before S=0.14145S=0.14145, put

P4=K4Z4(T),LSS=K2n2T[Z4(T)+18hˉZ2T/(K⋅0.01)].P_4=K^4Z_4(T),\qquad L_{SS}=K^2\sqrt{n_{2T}} [\sqrt{Z_4(T)}+18\bar h\sqrt{Z_{2T}}/(K\cdot0.01)]. (55)

The unshortened positive constructions define the values used in later interfaces; their shortened displays are P4<1.075567⋅10−76P_4<1.075567\cdot10^{-76} and LSS<2.534075⋅10−57L_{SS}<2.534075\cdot10^{-57} uniformly through preparation.

B.5 Original-stock transfer to the writer

Let Ψo\Psi_o be the exact preparation without the writer, and Ψn\Psi_n the enlarged exact wave, with their common initial state Ψn(0)=Ψo(0)g0\Psi_n(0)=\Psi_o(0)g_0. Here g0(Z)=π−1/4e−Z2/2g_0(Z)=\pi^{-1/4}e^{-Z^2/2}. Write Δ=Ψn−Ψog0\Delta=\Psi_n-\Psi_og_0. The ground oscillator energy is subtracted and the additional potential is J=(ℏ/T0)Π1[Cw(s)−Fw(s)Z]J=(\hbar/T_0)\Pi_1[C_w(s)-F_w(s)Z]. Thus

iℏ∂tΔ=HnΔ+JΨog0,Δ(0)=0. i\hbar\partial_t\Delta=H_n\Delta+J\Psi_og_0,\qquad\Delta(0)=0.

No comparison flow or initial law is substituted in this identity. The writer is supported on S≥0.14145S\ge0.14145. Throughout preparation, the old helper and its bounded phase both vanish on this support. The source amplitudes and momentum norms are therefore bounded by n2T,LS,Lr,LSSn_{2T},L_S,L_r,L_{SS} established above.

The following constants are exact rational majorants. Put μ=0.01\mu=0.01, w=0.003w=0.003 and

(b0,b1,b2,b3,b4)=24(1,315128w,252064w2,7560/16+5040/64w3,1360800w4). (b_0,b_1,b_2,b_3,b_4)=24\left(1,\frac{315}{128w}, \frac{2520}{64w^2},\frac{7560/16+5040/64}{w^3}, \frac{1360800}{w^4}\right).

In this subsection bjb_j denotes a writer derivative bound, not the preparation dilation or bump derivative. Define

Fj=bj/μ+μbj+2(0≤j≤2),C0=b02/(2μ)+μb0b2/2,C1=b0b1/μ+μ(b1b2+b0b3)/2,C2=(b12+b0b2)/μ+μ(b22+2b1b3+b0b4)/2,c0=b0+μb1,c1=b1+μb2.\begin{aligned}F_j&=b_j/\mu+\mu b_{j+2}\quad(0\le j\le2),\\ C_0&=b_0^2/(2\mu)+\mu b_0b_2/2,\\ C_1&=b_0b_1/\mu+\mu(b_1b_2+b_0b_3)/2,\\ C_2&=(b_1^2+b_0b_2)/\mu +\mu(b_2^2+2b_1b_3+b_0b_4)/2,\\ c_0&=b_0+\mu b_1,& c_1&=b_1+\mu b_2. \end{aligned}

These bound Fw,CwF_w,C_w and their dimensionless derivatives. Hereafter Fold=8⋅10−6+10−17F_{\rm old}=8\cdot10^{-6}+10^{-17} and Fold,r=10−6F_{{\rm old},r}=10^{-6}. These force ceilings follow by differentiating the displayed potential: the preparation clock force is bounded by

ma2(0.080802)/2+ma0a1(0.201)2+m2a1a2(0.080802)2/(4M)<8⋅10−6, ma_2(0.080802)/2+ma_0a_1(0.201)^2 +m^2a_1a_2(0.080802)^2/(4M)<8\cdot10^{-6},

and its radial force by ma1(0.201)+ma02(0.201)+m2a12(0.080802)(0.201)/(2M)ma_1(0.201)+ma_0^2(0.201)+m^2a_1^2(0.080802)(0.201)/(2M). Adding the quiet and activation derivatives gives the stated ceilings; their clock contribution is below 10−1710^{-17}. A second differentiation gives the conservative bound Fold,SS≤1F_{{\rm old},SS}\le1: its preparation part is

ma3(0.080802)/2+m(a12+a0a2)(0.201)2+m2(a22+a1a3)(0.080802)2/(4M), ma_3(0.080802)/2+m(a_1^2+a_0a_2)(0.201)^2 +m^2(a_2^2+a_1a_3)(0.080802)^2/(4M),

with QSS+γ2(W0+Q0)+2γ1qt1Q_{SS}+\gamma_2(W_0+Q_0)+2\gamma_1q_{t1} added.

For concise exact formulas set no=n2Tn_o=n_{2T} and

dΔ=(C0+F0/2)no/T0,qΔ=(C0+3/2F0)no/T0,rS=(C0+F0/2)LS/T0+ℏ(C1+F1/2)no/T02,rr=(C0+F0/2)Lr/T0.\begin{align}d_\Delta&=(C_0+F_0/\sqrt2)n_o/T_0,& q_\Delta&=(C_0+\sqrt{3/2}F_0)n_o/T_0,\tag{56}\\ r_S&=(C_0+F_0/\sqrt2)L_S/T_0 +\hbar(C_1+F_1/\sqrt2)n_o/T_0^2,\tag{57}\\ r_r&=(C_0+F_0/\sqrt2)L_r/T_0. \tag{58}\end{align}

The exact oscillator identities ∥Zg0∥=∥pZg0∥=1/2\|Zg_0\|=\|p_Zg_0\|=1/\sqrt2 and ∥Z2g0∥2=∥pZZg0∥2=3/4\|Z^2g_0\|^2=\|p_ZZg_0\|^2=3/4 prove these coefficients. Unitary Duhamel and oscillator rotation in the direct sum of its two quadratures then give increasing bounds

DΔ(t)=dΔt,QΔ(t)=qΔt+F0dΔt2/(2T0),PΔ,S(t)=rSt+FolddΔt2/2+ℏT02[F1(qΔt2/2+F0dΔt3/(6T0))+C1dΔt2/2],PΔ,r(t)=rrt+Fold,rdΔt2/2.\begin{align}D_\Delta(t)&=d_\Delta t,\tag{59}\\ Q_\Delta(t)&=q_\Delta t+F_0d_\Delta t^2/(2T_0),\tag{60}\\ P_{\Delta,S}(t)&=r_St+F_{\rm old}d_\Delta t^2/2\tag{61}\\ &\quad+\frac\hbar{T_0^2} [F_1(q_\Delta t^2/2+F_0d_\Delta t^3/(6T_0)) +C_1d_\Delta t^2/2],\tag{62}\\ P_{\Delta,r}(t)&=r_rt+F_{{\rm old},r}d_\Delta t^2/2. \tag{63}\end{align}

Here QΔQ_\Delta bounds (∥ZΔ∥2+∥pZΔ∥2)1/2(\|Z\Delta\|^2+\|p_Z\Delta\|^2)^{1/2}, not an unbounded coordinate times a bare norm estimate. In the clock inequality the term F1QΔF_1Q_\Delta is retained explicitly.

Let B=∂Z+ZB=\partial_Z+Z. Its exact equations have sources −Π1FwΨn-\Pi_1F_w\Psi_n and −2Π1FwBΨn-2\Pi_1F_wB\Psi_n for BΨnB\Psi_n and B2ΨnB^2\Psi_n, respectively. Because BB annihilates the old product wave, integration gives

B1(t)=F0(not+dΔt2/2)/T0,B2(t)=F02(not2+dΔt3/3)/T02,Ah=B1(T)+c0[no+DΔ(T)],A2h=B2(T)+2c0B1(T)+c02[no+DΔ(T)].\begin{align}B_1(t)&=F_0(n_ot+d_\Delta t^2/2)/T_0,\tag{64}\\ B_2(t)&=F_0^2(n_ot^2+d_\Delta t^3/3)/T_0^2,\tag{65}\\ A_h&=B_1(T)+c_0[n_o+D_\Delta(T)],\tag{66}\\ A_{2h}&=B_2(T)+2c_0B_1(T)+c_0^2[n_o+D_\Delta(T)]. \tag{67}\end{align}

The last quantity bounds the norm of the squared moving annihilator; it is not the square of its norm. These are operator identities on the full oscillator Gaussian domain, with no truncation of ZZ. For f=BΨnf=B\Psi_n, oscillator algebra gives ∥Zf∥2+∥pZf∥2=∥Bf∥2+∥f∥2\|Zf\|^2+\|p_Zf\|^2=\|Bf\|^2+\|f\|^2. Consequently, with integrals over [0,T][0,T],

PB=Fold∫B1+ℏT02[F1∫(B2+B1)+C1∫B1]+F0T0[LST+TPΔ,S(T)]+ℏF1T02(noT+dΔT2/2),Uh=[PB+c0(LS+PΔ,S(T))+ℏc1(no+DΔ(T))/T0]/ℏ.\begin{aligned}P_B={}&F_{\rm old}\int B_1 +\frac\hbar{T_0^2}\left[F_1\int(B_2+B_1)+C_1\int B_1\right]\\ &+\frac{F_0}{T_0}[L_ST+TP_{\Delta,S}(T)] +\frac{\hbar F_1}{T_0^2}(n_oT+d_\Delta T^2/2),\\ U_h={}&[P_B+c_0(L_S+P_{\Delta,S}(T)) +\hbar c_1(n_o+D_\Delta(T))/T_0]/\hbar. \end{aligned}

This proves ∥∂S(AΨn)(T)∥≤Uh\|\partial_S(A\Psi_n)(T)\|\le U_h. All displayed integrals are polynomials; for example ∫B1=F0(noT2/2+dΔT3/6)/T0\int B_1=F_0(n_oT^2/2+d_\Delta T^3/6)/T_0 and ∫B2=F02(noT3/3+dΔT4/12)/T02\int B_2=F_0^2(n_oT^3/3+d_\Delta T^4/12)/T_0^2.

The exact symbolic entrance values used subsequently are

Dh=3.006683⋅10−11+DΔ(T),nh=n1T+DΔ(T),Ph=ℏp1+PS,o+PΔ,S(T),Pr,h=Pr,o+PΔ,r(T)+ϵq,0≤ϵq<10−100,Rh=1+QΔ(T)<2.\begin{align}D_h&=3.006683\cdot10^{-11}+D_\Delta(T),\tag{68}\\ n_h&=n_{1T}+D_\Delta(T),\tag{69}\\ P_h&=\hbar p_1+P_{S,o}+P_{\Delta,S}(T),\tag{70}\\ P_{r,h}&=P_{r,o}+P_{\Delta,r}(T)+\epsilon_q, \qquad 0\le\epsilon_q<10^{-100},\tag{71}\\ R_h&=1+Q_\Delta(T)<2. \tag{72}\end{align}

Here ϵq\epsilon_q is the quiet-reference comparison error in the radial momentum, as bounded in the activation appendix. In particular the larger norm 3.006684⋅10−113.006684\cdot10^{-11} used in the fourth-moment estimate is not substituted for DhD_h.

B.6 Localized repair of the transfer Hessian

A coarse second-momentum bound is useful only to generate source-free left collars. We give it explicitly to avoid circular reuse of the improved result. Let P=PΔ,S(T)P=P_{\Delta,S}(T), D=DΔ(T)D=D_\Delta(T), Q=QΔ(T)Q=Q_\Delta(T), and

rSQ=(C0+3/2F0)LS/T0+ℏ(C1+3/2F1)no/T02,QS=F0PT/T0+Fold∫QΔ+ℏT02[C1∫QΔ+F1∫(B2+2B1+2DΔ)]+rSQT,rSS(L)=(C0+F0/2)L/T0+2ℏ(C1+F1/2)LS/T02+ℏ2(C2+F2/2)no/T03,WΔ=2ℏT(C1P+F1QS)/T02+ℏ2T(C2D+F2Q)/T03,Hc=rSS(HSS)T+2FoldPT+ℏDT+WΔ.\begin{aligned}r_{SQ}&=(C_0+\sqrt{3/2}F_0)L_S/T_0 +\hbar(C_1+\sqrt{3/2}F_1)n_o/T_0^2,\\ Q_S&=F_0PT/T_0+F_{\rm old}\int Q_\Delta +\frac\hbar{T_0^2}\left[C_1\int Q_\Delta +F_1\int(B_2+2B_1+2D_\Delta)\right]+r_{SQ}T,\\ r_{SS}(L)&=(C_0+F_0/\sqrt2)L/T_0 +2\hbar(C_1+F_1/\sqrt2)L_S/T_0^2 +\hbar^2(C_2+F_2/\sqrt2)n_o/T_0^3,\\ W_\Delta&=2\hbar T(C_1P+F_1Q_S)/T_0^2 +\hbar^2T(C_2D+F_2Q)/T_0^3,\\ H_c&=r_{SS}(H_{SS})T+2F_{\rm old}PT+\hbar DT+W_\Delta. \end{aligned}

The quadrature estimate ∥Z2Δ∥\|Z^2\Delta\| and its companion are controlled by B2+2B1+2DΔB_2+2B_1+2D_\Delta, so that QSQ_S bounds the mixed weighted momentum without replacing ZZ by a bounded operator. Differentiating the inhomogeneous equation twice proves ∥pS2Δ∥≤Hc\|p_S^2\Delta\|\le H_c; the coarse value is below ℏ2(2.486519⋅1027)\hbar^2(2.486519\cdot10^{27}).

Choose decreasing left cutoffs with transitions [0.126,0.136][0.126,0.136] and [0.116,0.126][0.116,0.126], width g=0.01g=0.01. The writer forcing vanishes on both supports and every old clock force is contained in their full-one regions, because it vanishes for S≥0.106S\ge0.106. Positive drift favours exit from each left region. Continuity and one integration by parts give

nΔ1=TP/(Mg),nΔ2=TMg[nΔ1Hc+2ℏD/g].n_{\Delta1}=TP/(Mg),\qquad n_{\Delta2}=\frac{T}{Mg} [\sqrt{n_{\Delta1}H_c}+2\hbar D/g]. (73)

For the new Hessian H(t)=∥pS2Δ(t)∥H(t)=\|p_S^2\Delta(t)\| the old-force momentum is instead localized by ∥1oldforcepSΔ∥≤nΔ2H(t)+2ℏnΔ1/g\|\one_{\rm oldforce}p_S\Delta\| \le\sqrt{n_{\Delta2}H(t)}+2\hbar n_{\Delta1}/g. Use the proved far old-wave Hessian LSSL_{SS} in the source. The differential inequality has the form H′≤2αH+fH\prime\le2\alpha\sqrt H+f, where α=FoldnΔ2\alpha=F_{\rm old}\sqrt{n_{\Delta2}}. The function [αt+∫0tf]2[\alpha t+\sqrt{\int_0^t f}]^2 is a supersolution, because its derivative is at least 2αH+f2\alpha\sqrt H+f when evaluated at that function. Bounding the nonnegative source integral gives

HΔ=[FoldTnΔ2+rSS(LSS)T+4FoldℏnΔ1T/g+ℏnΔ2T+WΔ]2.H_\Delta=\left[F_{\rm old}T\sqrt{n_{\Delta2}}+ \sqrt{r_{SS}(L_{SS})T+4F_{\rm old}\hbar n_{\Delta1}T/g +\hbar n_{\Delta2}T+W_\Delta}\right]^2. (74)

It proves ∥pS2Δ(T)∥≤HΔ\|p_S^2\Delta(T)\|\le H_\Delta and HΔ/ℏ2<1.828510⋅1019H_\Delta/\hbar^2<1.828510\cdot10^{19}. The coarse Hessian only generated the amplitude collar; it was not relabelled as the improved source or inserted into this nonlinear old-force term.

B.7 Handoff enclosures and the preparation event

The constructions above give the following outward decimal enclosures. Every endpoint in the table is a rational number. The unrounded expressions, not shortened displays, define the quantities when used in the final event sum.

QuantityStrict upper bound
DoD_o3.00668265254473406785⋅10−113.00668265254473406785\cdot10^{-11}
PS,oP_{S,o}6.98319079275307899984⋅10−326.98319079275307899984\cdot10^{-32}
n1Tn_{1T}6.24924697825783043607⋅10−306.24924697825783043607\cdot10^{-30}
HSSH_{SS}3.23800849281954087759⋅10−493.23800849281954087759\cdot10^{-49}
HSSoH_{SS}^{o}3.24499145042629721024⋅10−493.24499145042629721024\cdot10^{-49}
n2Tn_{2T}5.97036914776425754723⋅10−385.97036914776425754723\cdot10^{-38}
P4P_41.07556686771573853653⋅10−761.07556686771573853653\cdot10^{-76}
LSSL_{SS}2.53407404062459981525⋅10−572.53407404062459981525\cdot10^{-57}
DΔ(T)D_\Delta(T)2.767839⋅10−302.767839\cdot10^{-30}
PΔ,S(T)P_{\Delta,S}(T)7.597256⋅10−367.597256\cdot10^{-36}
PΔ,r(T)P_{\Delta,r}(T)2.08385966319223143790⋅10−372.08385966319223143790\cdot10^{-37}
AhA_h5.04959232864420059458⋅10−245.04959232864420059458\cdot10^{-24}
A2hA_{2h}1.22824296520291159974⋅10−171.22824296520291159974\cdot10^{-17}
UhU_h276122.969821611569892276122.969821611569892
PhP_h3.20281666862106903351⋅10−313.20281666862106903351\cdot10^{-31}
HΔ/ℏ2H_\Delta/\hbar^21.82850947100693947781⋅10191.82850947100693947781\cdot10^{19}

B.8 Preparation history under the enlarged exact flow

The preparation rank K(S,r)K(S,r) is independent of ZZ; hence its material derivative under the enlarged flow has no pointer-current term. Its explicit drift β\beta in (34) is unchanged. The error in its current must nevertheless be recomputed for that flow. In the bounded gauge the additional momentum ceilings are

drg=PΔ,r(T)+ma0(0.201)DΔ(T),dSg=PΔ,S(T)+ma1(0.080802)DΔ(T)/2. d_r^g=P_{\Delta,r}(T)+ma_0(0.201)D_\Delta(T),\qquad d_S^g=P_{\Delta,S}(T)+ma_1(0.080802)D_\Delta(T)/2.

The exact normalized wave stays in the derivative-error factor of the bilinear expansion, retaining the quadratic error terms. Consequently

ΔI=T[(ℏcr/m+ℏcS/M)DΔ(T)+(kr/m)drg+(kS/M)dSg]<1.304454833⋅10−22.\Delta I=T\left[(\hbar c_r/m+\hbar c_S/M)D_\Delta(T) +(k_r/m)d_r^g+(k_S/M)d_S^g\right] <1.304454833\cdot10^{-22}. (75)

Under the single original joint domination constant C=270000/67499C=270000/67499, NN fixed radial cuts therefore have preparation failure probability bounded by

Ep(N)=2CNβ+2C2N(Io+ΔI)+C[Do+DΔ(T)]2.E_p(N)=2CN\beta+2C\sqrt{2N(I_o+\Delta I)} +C[D_o+D_\Delta(T)]^2. (76)

This follows by taking initial rank collars of width β+ϵ\beta+\epsilon, using Markov only on the unknown absolute current, and minimizing in ϵ\epsilon. The last term covers the terminal clock suffix where the helper vanishes. It is an estimate on the enlarged flow with its original initial law; closeness of two waves alone has not been used to assert history agreement. The initial conditional-CDF interface costs, separately,

Ecdf(N)=2CN[Do+DΔ(T)].E_{\rm cdf}(N)=2CN[D_o+D_\Delta(T)]. (77)

For the final cut family, including the exterior ray, use N=256001N=256001 in both formulas. Evaluating the finite expressions gives

Ep(256001)<0.000003084949592860634224432805,Ecdf(256001)<0.0000615780135255954031515710. \begin{aligned} E_p(256001)&<0.000003084949592860634224432805,\\ E_{\rm cdf}(256001)&<0.0000615780135255954031515710. \end{aligned}

The pointer remains part of the one jointly dominated auxiliary stock; no independent pointer cap, sector sample, equilibrium reset, or new law at handoff enters these estimates.