The positive collective kinetic coefficient used in the inverse source model must have an origin in the same field. It does not follow merely from allowing the occupied carrier wave to respond. That wave's exact frequency-dependent response can be computed with all its signed modes retained, and its dominant low-frequency time-derivative coefficient has the opposite sign.

This is a calculation of a forced wave response. It is not negative particle inertia. The complete positive-frequency particle and its sustaining motion still require their joint energy and momentum.

The retained geometric mode and the full response

Use the actual finite projected linked normal operator, including its spatially computed hopping integrals and finite-hopping frequency corrections. The constant-polynomial deformation has the two Hermitian source matrices O1,O2O_1,O_2 derived in the reciprocal source calculation. In its stated scaled coordinate,

H(a)=H0+∑jajOj.(1) H(a)=H_0+\sum_j a_jO_j. \tag{1}

Let H0ψn=λnψnH_0\psi_n=\lambda_n\psi_n, with unit norm and λn>0\lambda_n>0. For every other positive or negative eigenstate define

νl=λl−λn,blj=⟨ψl,Ojψn⟩.(2) \nu_l=\lambda_l-\lambda_n,\qquad b_{lj}=\langle\psi_l,O_j\psi_n\rangle. \tag{2}

The mean gj,n=⟨ψn,Ojψn⟩g_{j,n}=\langle\psi_n,O_j\psi_n\rangle contributes the diagonal phase response. Subtracting that mean from the forcing removes the common phase direction from the shape calculation.

For a real drive with Fourier convention e−iωte^{-i\omega t}, solve both the ordinary and conjugate wave equations. Their combined current response is

δgj(ω)=∑kχjk(ω)δak(ω),χjk(ω)=∑l≠n[blj∗blkω−νl−bljblk∗ω+νl].(3) \boxed{\delta g_j(\omega)=\sum_k\chi_{jk}(\omega)\delta a_k(\omega), \qquad \chi_{jk}(\omega)=\sum_{l\ne n} \left[ \frac{b_{lj}^*b_{lk}}{\omega-\nu_l} -\frac{b_{lj}b_{lk}^*}{\omega+\nu_l} \right].} \tag{3}

An upper-half-plane frequency specifies the retarded response. The poles are the actual rotating wave frequencies. No signed partner has been replaced by a positive-frequency approximation.

Static, Berry and time-derivative coefficients

Below the first nonzero response pole, expansion of (3) gives

χ(ω)=K+iωF−ω2I+O(ω3),(4) \chi(\omega)=K+i\omega F-\omega^2 I+O(\omega^3), \tag{4}

Kjk=−2∑lRe⁡(blj∗blk)νl,Fjk=−2∑lIm⁡(blj∗blk)νl2,Ijk=2∑lRe⁡(blj∗blk)νl3.(5) \boxed{ \begin{aligned} K_{jk}&=-2\sum_l\frac{\operatorname{Re}(b_{lj}^*b_{lk})}{\nu_l},\\ F_{jk}&=-2\sum_l\frac{\operatorname{Im}(b_{lj}^*b_{lk})}{\nu_l^2},\\ I_{jk}&=2\sum_l\frac{\operatorname{Re}(b_{lj}^*b_{lk})}{\nu_l^3}. \end{aligned}} \tag{5}

Here KK is the static frequency Hessian of the retained linear deformation. FF is its antisymmetric Berry curvature in the convention Aj=i⟨ψn,∂jψn⟩A_j=i\langle\psi_n,\partial_j\psi_n\rangle. II is the coefficient of a˙ja˙k/2\dot a_j\dot a_k/2 in the wave-induced local effective Lagrangian. The entire expression carries the wave's common canonical norm when that norm is not one.

These coefficients are not independent phenomenological parameters. They follow from the same signed resolvent. Nor can the derivative expansion replace (3) near a pole. A signed coefficient in that expansion does not establish an instability of the full Hermitian forced-wave problem.

The source matrices are derivatives with respect to the scaled aa. For the physical polynomial coordinate ϵ\epsilon, the conversion a=γκdϵ/Cda=\gamma\kappa_d\epsilon/C_d must also be applied to (5). Both quadratic coefficients acquire the corresponding squared coordinate factor. Their units cannot be inferred from an untransformed numerical table.

The transverse partner fixes the leading sign

In the large-hopping two-endpoint limit, O1O_1 is diagonal on the selected positive state and O2O_2 connects it to its negative partner. That partner has ν=−2λn\nu=-2\lambda_n and unit transverse matrix element. Therefore

χ22(ω)=4λn4λn2−ω2,K22=1λn,I22=−14λn3.(6) \boxed{ \chi_{22}(\omega)=\frac{4\lambda_n}{4\lambda_n^2-\omega^2},\quad K_{22}=\frac1{\lambda_n},\quad I_{22}=-\frac1{4\lambda_n^3}.} \tag{6}

The exact finite carrier retains the small fast-band contributions as well. At γ=10000\gamma=10000, 130-digit calculations give the following dimensionless values in the common coordinate convention:

Maintained control Induced transverse coefficient I22I_{22}
Depth-three down −0.250000026526-0.250000026526
Depth-four up −7.40741078150×10−5-7.40741078150\times10^{-5}
Depth-six electron −3.90775627193×10−14-3.90775627193\times10^{-14}
Depth-six muon −4.42087631510×10−21-4.42087631510\times10^{-21}
Depth-ten tau −9.31151593660×10−25-9.31151593660\times10^{-25}

The depth-five neutrino controls and all remaining maintained quark controls obey the same sign and inverse-cube leading law. Both surface signs are included. Across the full scan, −4λn3I22-4\lambda_n^3I_{22} remains close to one; the largest departure is about 2.2×10−42.2\times10^{-4}. The exact finite frequency is used in that ratio, rather than substituting the ideal simplex count.

The complete calculation also retains the much smaller longitudinal coefficient and nonzero Berry curvature. Real and imaginary source quadratures make the symmetric cross coefficient zero in this convention; their antisymmetric response is retained in (3)–(5).

Consequence for source reconstruction

Eliminating this occupied wave does not generate the positive, band-independent geometric inertia assumed in the collective completion. It supplies (5), with the dominant contribution (6). A positive complete geometric kinetic term would therefore need the sustaining modes' own contribution, or a reconstruction retaining the full response instead of assigning this coefficient an independent positive value.

This conclusion applies to the retained linear geometric deformation and the full finite carrier. The full normal complement, the physical kinetic law of the supporting texture, higher derivatives of the deformed spatial operator and any additional moving-basis terms have not been included in that claim. They must be computed before identifying a complete geometric source or a physical particle's inertia.

Verification

carrier_induced_response.py reconstructs every signed carrier eigenvector from the actual normal-mode hopping integrals. An independent inversion of the full ordinary and conjugate matrices agrees with (3), component by component.

carrier_induced_checks.py also extracts (5) from finite-frequency responses and checks every maintained control at depths 3, 4, 5, 6 and 10 on both surface signs. Its 130-check report passes at 130-digit precision. No amplitude was fitted and no physical relative sector scale was inferred from these dimensionless coefficients.