The form-dual estimates above and the following estimates solve different input problems. A residual with a spatial form norm can be propagated using actual forces without paying an irrelevant constant carrier energy. A residual with only a Hilbert norm can instead use temporal regularity and a common strong domain. Neither hypothesis implies the other for singular interactions.
Theorem 8.1 (Forced positive form with complete first-order work)
Use (2.1) and K=−∑iλiDi2+U≥I, with scalar U≥1. Assume a common form domain, the true unitary propagator, and a form-stable approximation justifying the commutators below. Suppose K′≤aKK as forms and, for form vectors v,
Their absolute values are at most 2LQ2 and 2Q(G0e+G1Q). Also i[F,U]=∑iai∂iU, whose quadratic form is bounded by Q(u0e+u1Q). Thus, retaining the forcing,
(Q2)′≤2αQ2+2βeQ+2Q∥K1/2r∥.
Regularize Q at zero and use Duhamel's norm bound for e. The integrating factor proves (8.1).
For form data, set Rϵ=(I+ϵK)−1 and Kϵ=K(I+ϵK)−1. The ordered identities
Kϵ′=RϵK′Rϵ,[V+F,Kϵ]=Rϵ[V+F,K]Rϵ
give the same inequality with Qϵ=∥Kϵ1/2E∥. Indeed ∥RϵE∥≤e, ∥K1/2RϵE∥≤Qϵ, and the forcing is bounded by Qϵ∥K1/2r∥. Integrate before taking the monotone spectral limit ϵ↓0. The assumed conforming approximation supplies the propagator and commutator passage for unbounded coefficients. Resolvent algebra alone does not prove those domain premises.
□
When ai=0, this reduces to the spatial-force estimate with α=aK/2+g1 and β=g0. A covariantly constant carrier M(t) has ∇iM=0 and contributes no spatial force. Its phase is still in the true evolution. A laboratory-constant matrix can instead have [Ci,M]=0 and must then be priced. For a scalar dilation a=cx in one flat coordinate, the kinetic contribution is −2c∥∂xE∥2, confirming the commutator sign. Spatially varying mobility or a discontinuous coefficient needs additional terms or a different theorem.
Proposition 8.2 (Transported quadratic form and physical coflow)
Let H0(t) be a real scalar Weyl-ordered quadratic rate Hamiltonian in z=(q,p), p=−i∂q, including its affine terms. Let S(t) be its classical symplectic propagator and m(t) its affine classical solution. For any positive matrix W0, set
W(t)=S(t)W0S(t)T,K(t)=1+41(z−m)TW(t)−1(z−m).
Then K≥1 and, as forms on the transported Schwartz core,
Kt+i[H0,K]=0.(8.2)
For a pure Gaussian with phase gradient mp+R(q−mq), R=RT, and positive imaginary phase matrix I, its covariance has blocks
Writing D0=∂−i(mp+R(q−mq)) and Λ=21I−1, one has K≥D0†ΛD0. The Gaussian ground value of K is 1+d/2.
Proof
If the classical quadratic generator is A=JG0, then W′=AW+WAT and (W−1)′=−ATW−1−W−1A. The affine equation for m cancels the linear terms. Commutators of Weyl quadratics equal their Poisson expressions, since higher Moyal derivatives vanish; these identities prove (8.2). A positive quadratic matrix is a sum of squares of real linear canonical observables, proving positivity. Factoring the displayed covariance as
(1R01)(C00I/2)(10R1)
and inverting proves the form formula. Each centered Gaussian coordinate contributes 1/2 to its quadratic ground value.
□
This is established invariant mathematics [12, 13]; its role here is a positive error metric through an inverted interval. It is dimensionless and does not replace the physical mass in the current. For the same laboratory vector potential A0, define b0,i=(ℏmp,i+ℏ(R(q−mq))i−eiA0,i)/Mi. For any true wave in that constitution,
ji[Ψ]−ρΨb0,i=MiℏIm(Ψ†D0,iΨ).
Consequently its difference from the analogous expression for P is at most
Miℏ(e∥D0,iP∥+N(Λ−1)iiQ)(8.3)
in integrated absolute value. The weighted coordinate Cauchy inequality proves the last factor. A changed potential ΔA adds −eiΨ†ΔAiΨ/Mi to the relative current.
For H=H0+V with Hermitian multiplication V and form forcing r, the exact cancellation (8.2) gives Q′≤g1Q+g0e+∥K1/2r∥ if ∥Λ1/2∇Vv∥≤g0∥v∥+g1∥K1/2v∥. The proof is the preceding commutator calculation, with no separate Kt loss. Resolvent regularization preserves Kϵ,t+i[H0,Kϵ]=0. A canonical drift rewritten in D0 also creates the multiplication term a⋅(mp+R(q−mq)); its force and literal current remain. An internal-branch-dependent covariance would produce further commutators and is outside this scalar invariant statement.
Theorem 8.3 (First-domain causal estimate with a form-controlled time derivative)
Let H(t) be self-adjoint on a fixed domain D, with t↦H(t)∈L(D,H) continuously differentiable for an equivalent fixed graph norm. Assume its unitary common-domain propagator is locally bounded in that graph norm. Let Ds be a fixed finite Hermitian right-clock matrix, and write L(t)X=H(t)X−XDs. For an absolutely continuous Hilbert–Schmidt residual r, with r′∈L1, consider
iδ′=Lδ−r,δ(0)=0.
Suppose K=H+b(t)≥I and, on D,
∥(H′−ω′)v∥≤c1(t)∥K1/2v∥+c0(t)∥v∥
for a declared real differentiable scalar ω and integrable nonnegative c0,c1. For each T define
When these quantities are finite, the causal solution belongs to the common strong domain and satisfies
[0,T]sup∥K1/2δ∥≤ETAT+CT+21ETGT.(8.4)
No spatial form norm of r and no application of H twice to δ is assumed.
Proof
We first justify the domain, rather than use a formal H2 calculation. Set B(t)=(L(t)−i)−1. It maps the Hilbert–Schmidt space boundedly into its common graph domain and B′=−BL′B. If U(t,s) is the two-sided propagator, integration by parts in the graph norm gives
Indeed ∂sU(t,s)=iU(t,s)L(s) and iLB=i−B, so differentiating UBr recovers iUr minus the displayed integrand. This is precisely the Duhamel error. The terms are graph-integrable; absolutely continuous approximation in time proves the identity for the stated data.
Put η=(L−ω)δ−r. Testing against common-domain adjoint solutions, or using (8.5), yields the mild identity
iη′=Lη−i(r′+iωr)+i(H′−ω′)δ,η(0)=−r(0).
It does not assert η∈D. Let ℓ=∥(L−ω)δ∥, e=∥δ∥ and Q=∥K1/2δ∥. Unitarity and the commuting right multiplication give e(t)≤E(t)=∫0t∥r∥ and ∥δDs∥≤ED(t)=∫0t∥rDs∥. The auxiliary equation and Kδ=(L−ω)δ+(b+ω)δ+δDs therefore imply
ℓ(t)≤AT+∫0tc1Q,Q(t)2≤ETℓ(t)+CT.
For ET>0, set M(t)=AT+∫0tc1Q. Then M′≤c1ETM+CT, and differentiation of the square root, with a positive regularizer if necessary, gives ETM+CT≤ETAT+CT+ET∫0tc1/2. This proves (8.4). If ET=0, the causal error is identically zero.
□
The common-domain evolution premise is a standard nonautonomous evolution requirement; the regularity equivalence in [9] is relevant to checking it. The estimate retains the actual right-clock defect ED,T and the covariant temporal combination r′+iωr. A large internal operator is not a removable scalar phase.
For example, with Πj=pj−ej(A0,j+aj(t)) in physical units and rate generator Hkin=∑jΠj2/(2Mjℏ),
A linear force Fj(t)xj is also allowed when fixed positive traps provide a position form: Young's inequality supplies a shift b≥1+∥Vadd−∥/ℏ+∑jFj2/(MjΩj2ℏ) and domination of MjΩj2xj2/(4ℏ). Its additional contribution to c12 is ∑j4∣Fj′∣2/(MjΩj2ℏ). A changing quadratic trap or an unbounded source-coordinate field needs its actual stronger graph bound.
Proposition 8.4 (Coulomb residual and temporal coordinates)
In three relative spatial dimensions, let ϕ be a smooth Gaussian nonzero at a collision. Then ∣x∣−1ϕ∈Lloc2 but ∣x∣−1ϕ∈/Hloc1. Nevertheless, if V(x)=g/∣x∣ is fixed in laboratory coordinates, P2(t,x) is a differentiable quadratic comparison potential and ϕt is a polynomial times a smooth finite Gaussian family, then
f=(V−P2)ϕ,ft=−P2,tϕ+(V−P2)ϕt
are locally in L2 and globally in L2 for the indicated Gaussian polynomial families. On a finite parameter path with nondegenerate covariance they are continuous in L2. Under a homogeneous point change x=L(t)X, with invertible L,
∣∂t(g/∣L(t)X∣)∣≤∥L′L−1∥∣g∣/∣L(t)X∣.
Proof
The radial square integrals near zero have orders ∫01dr for 1/r and ∫01r−2dr for its leading gradient. The nonzero Gaussian value prevents cancellation of the latter leading term. Polynomial factors and Gaussian tails give the global L2 statements; local domination by an integrable r−2 square envelope and uniform Gaussian tail bounds prove continuity. Finally differentiation gives −gx⋅(L′L−1x)/∣x∣3, proving the last inequality.
□
An exponentially small nonzero Gaussian value does not remove the divergent spatial derivative. A translation that moves the collision point instead produces a 1/r2 temporal multiplier. The homogeneous repair retains the physical centroid, linear force, boost and scalar action; recentering them away may restore that moving singularity. These observations select a valid response method, without claiming a small apparatus error.
For completeness a strong-graph variant applies when a changing trap has a graph-relative rather than form-relative time derivative.
Proposition 8.5 (One graph by an ordered inverse lift)
Let H(t) satisfy the common-domain assumptions above, and let A=H+c≥aI for constants a>0,c≥0. Assume γ(t)=∥H′(t)A(t)−1∥∈L1. For iE′=HE+f, E(0)∈D and f∈W1,1([0,T];H), set R(t)=∫0tγ. Then
The auxiliary Y=AE+f obeys iY′=HY+cf+iH′E+if′; the terms containing Hf cancel. Solve its Volterra equation using E=A−1(Y−f) and the bounded perturbation H′A−1. The norm inequality for Y and Gronwall give (8.6). This is also a rigorous construction: (A−1)′=−A−1H′A−1 verifies the original equation weakly for the constructed E, and uniqueness of the Hilbert mild solution identifies it. No separate expression Hf is needed.
□
For a fixed laboratory Coulomb singularity and moving oscillator centres, the time derivative of the Coulomb kernel is zero. The derivative of the trap is quadratic and linear in the laboratory coordinates, so it can be graph-relative under the earned oscillator domain. Proposition 8.5, rather than an unjustified spatial derivative of f, is then the relevant interface. A detector moving with those centres still pays every term in (2.6).