Summary: The full response of a two-frequency spatial Dirac field reduces to a frequency-indexed moment equation that retains scalar-center motion, normal propagation and the conditions for source-sideband cancellation.

Follow the complete periodic background

The two-frequency spatial field has an unchanging scalar binding source and a changing number-current pattern. Its linearized field equation nevertheless has periodic coefficients. The response can be derived without assigning its bound frequency to an external forcing frequency or discarding one of its two components.

All coordinates describe ordinary macroscopic spatial directions of the one nested manifold. The bound coordinates are x∈Rdx\in\mathbb R^d and yy is one existing containing direction. The field-component indices, source moments and temporal harmonics below are mathematical descriptions of that same field; they add no physical spatial dimensions or particles. The present normal momentum is along the direction used to construct the opposite-frequency partner.

Let Hb=−iαi∂i+β(m0+asource∣x∣2)H_b=-i\alpha_i\partial_i+\beta(m_0+a_{\rm source}|x|^2) be the self-consistent bound operator, with {Hb,αy}=0\{H_b,\alpha_y\}=0. Choose normalized radial fields with possibly different rest polarizations,

Hbϕ1=Mϕ1,Hbϕ2=Mϕ2,Ψ0(t)=aϕ1e−iMt+bαyϕ2eiMt.(1) H_b\phi_1=M\phi_1,\qquad H_b\phi_2=M\phi_2,\qquad \Psi_0(t)=a\phi_1e^{-iMt}+b\alpha_y\phi_2e^{iMt}. \tag{1}

Their complete radial frame has the pointwise scalar identities ϕj†βϕj=σ0\phi_j^\dagger\beta\phi_j=\sigma_0 and ϕ1†βαyϕ2=0\phi_1^\dagger\beta\alpha_y\phi_2=0. Thus the binding density is sσ0s\sigma_0, where s=∣a∣2−∣b∣2>0s=|a|^2-|b|^2>0. Both profiles must be solved at that scalar weight and the same interaction coefficient. Set ∣a∣2+∣b∣2=1|a|^2+|b|^2=1 for the unit canonical norm used in the examples.

Retain shape and scalar-center variation

Use Hermitian spatial source operators

OA=β(xTAx),Oi=βxi,(2) O_A=\beta(x^TAx),\qquad O_i=\beta x_i, \tag{2}

where AA runs over a Frobenius-orthonormal basis of all symmetric spatial tensors, including breathing. Define CσC_\sigma as the column operator whose columns are OAψσO_A\psi_\sigma, with ψ+=aϕ1\psi_+=a\phi_1 and ψ−=bαyϕ2\psi_-=b\alpha_y\phi_2. Then

C(t)=∑σ=±1Cσe−iσMt.(3) C(t)=\sum_{\sigma=\pm1}C_\sigma e^{-i\sigma Mt}. \tag{3}

The real symmetric moment curvature is constant on this background. Its shape block is g0Ig_0I, and its dipole block is −2asourceI/Qσ-2a_{\rm source}I/Q_\sigma, with Qσ=∫Ψ0†βΨ0 dx≠0Q_\sigma=\int\Psi_0^\dagger\beta\Psi_0\,dx\ne0. Translational response derives that second block by varying the scalar center. An unchanged background center does not justify omitting its perturbations.

For the complexified ordinary and conjugate perturbation Z=(η,η∗)TZ=(\eta,\eta^*)^T, the full equation at normal momentum qq is

i∂tZ=[B0(q)+S(t)κT(t)]Z+S(t)f(t),(4) \boxed{ i\partial_tZ=\left[\mathcal B_0(q)+S(t)\kappa T(t)\right]Z +S(t)f(t),} \tag{4}

with

Hq=Hb+qαy,B0(q)=(Hq00−H−q∗),S(t)=(C(t)−C(t)∗),T(t)=(C(t)†,C(t)T),h=T(t)Z.(5) \begin{gathered} H_q=H_b+q\alpha_y,\qquad \mathcal B_0(q)=\begin{pmatrix}H_q&0\\0&-H_{-q}^*\end{pmatrix},\\ S(t)=\begin{pmatrix}C(t)\\-C(t)^*\end{pmatrix},\qquad T(t)=(C(t)^\dagger,C(t)^T),\qquad h=T(t)Z. \end{gathered} \tag{5}

The star means entrywise conjugation; the conjugate block retains the opposite normal momentum. These are laboratory field frequencies, with no stationary carrier subtracted from HqH_q. The additional source coefficients are ff, while hh contains the induced moments. A physical real perturbation includes the conjugate Fourier component at (−ω,−q)(-\omega,-q).

Equation (4) follows by the complete real-linear differentiation in Dirac source linearization. Its energy Hessian obeys ΣB=(ΣB)†\Sigma\mathcal B=(\Sigma\mathcal B)^\dagger, where Σ=diag⁡(I,−I)\Sigma=\operatorname{diag}(I,-I). This action identity does not make the energy positive. The field coefficients contain frequencies 0,±2M0,\pm2M and have period π/M\pi/M.

A monochromatic source produces three bare moment frequencies

First remove the feedback term κh\kappa h and impose fei(qy−ωt)f e^{i(qy-\omega t)}. Multiplication by (3) drives ordinary fields at ω+σM\omega+\sigma M and conjugate fields at ω−σM\omega-\sigma M. Reading those fields through T(t)T(t) gives moment outputs at ω+ΔM\omega+\Delta M, for Δ=0,±2\Delta=0,\pm2:

ΠΔ(ω,q)=∑σ−ρ=ΔCρ†(ω+σM−Hq)−1Cσ−∑ρ−σ=ΔCρT(ω−σM+H−q∗)−1Cσ∗.(6) \boxed{ \begin{aligned} \Pi_\Delta(\omega,q) ={}&\sum_{\sigma-\rho=\Delta} C_\rho^\dagger(\omega+\sigma M-H_q)^{-1}C_\sigma\\ &-\sum_{\rho-\sigma=\Delta} C_\rho^T(\omega-\sigma M+H_{-q}^*)^{-1}C_\sigma^*. \end{aligned}} \tag{6}

The first sum has ordinary propagation; the second has the conjugate response with its required sign. For example:

moment output ordinary contribution conjugate contribution
ω+2M\omega+2M C−†(ω+M−Hq)−1C+C_-^\dagger(\omega+M-H_q)^{-1}C_+ −C+T(ω+M+H−q∗)−1C−∗-C_+^T(\omega+M+H_{-q}^*)^{-1}C_-^*
ω−2M\omega-2M C+†(ω−M−Hq)−1C−C_+^\dagger(\omega-M-H_q)^{-1}C_- −C−T(ω−M+H−q∗)−1C+∗-C_-^T(\omega-M+H_{-q}^*)^{-1}C_+^*

Normal propagation can be evaluated with the full squared bound operator:

(ζ−Hq)−1=(ζ+Hb+qαy)(ζ2−Hb2−q2)−1,(ζ+H−q∗)−1=(ζ−Hb∗+qαy∗)(ζ2−(Hb∗)2−q2)−1.(7) \begin{aligned} (\zeta-H_q)^{-1} &=(\zeta+H_b+q\alpha_y)(\zeta^2-H_b^2-q^2)^{-1},\\ (\zeta+H_{-q}^*)^{-1} &=(\zeta-H_b^*+q\alpha_y^*) (\zeta^2-(H_b^*)^2-q^2)^{-1}. \end{aligned} \tag{7}

The domain and pole qualifications are those of the full normal response. No radial or angular field components have been eliminated by simply neglecting them.

Repeated feedback couples the frequency sidebands

Let hjh_j and fjf_j denote moment amplitudes at ω+2jM\omega+2jM. Define an operator on their frequency-indexed sequence by

Pjk(ω,q)={Π2(j−k)(ω+2kM,q),∣j−k∣≤1,0,∣j−k∣>1.(8) \mathbb P_{jk}(\omega,q)= \begin{cases} \Pi_{2(j-k)}(\omega+2kM,q),&|j-k|\le1,\\ 0,&|j-k|>1. \end{cases} \tag{8}

The complete source equation is

h=P(f+Kh),h=(I−PK)−1Pf,Kjk=δjkκ.(9) \boxed{h=\mathbb P(f+\mathbb K h),\qquad h=(I-\mathbb P\mathbb K)^{-1}\mathbb P f,\qquad \mathbb K_{jk}=\delta_{jk}\kappa.} \tag{9}

Its repeated feedback can populate more frequencies than the three in (6). A finite harmonic cutoff is an approximation to this sequence; convergence of a physical response or pole must be checked as that cutoff increases. At singular frequencies use the complete field operator or the appropriate limiting resolvent.

There is also an exact finite-cutoff identity useful for verification. Keep field frequencies ω+(2n+1)M\omega+(2n+1)M, −N≤n<N-N\le n<N, and moment frequencies ω+2jM\omega+2jM, −N≤j≤N-N\le j\le N. Let DD be the block diagonal frozen field resolvent denominator, and construct SF,TFS_F,T_F by matching the harmonics in (3)–(5). Then

PN=TFD−1SF,det⁡(D−SFKNTF)=det⁡D det⁡(I−PNKN).(10) \boxed{\mathbb P_N=T_FD^{-1}S_F,\qquad \det(D-S_F\mathbb K_NT_F) =\det D\,\det(I-\mathbb P_N\mathbb K_N).} \tag{10}

The frozen determinant factor retains source-free field modes and pole cancellations. A truncated moment equation and a differently truncated field equation need not agree at their frequency boundaries. Matching both cutoffs gives (10) exactly.

A periodic field can have a simpler scalar readout at rest

Periodic field coefficients do not guarantee a nonzero sideband in every observable. At q=0q=0, let PEP_E be the complete spectral projector of HbH_b at eigenvalue EE in a finite spatial regulator. Since αyPEαy=P−E\alpha_yP_E\alpha_y=P_{-E}, the positive sideband in (6) vanishes identically away from the poles if and only if every residue satisfies

C−†PEC+=(C−†P−EC+)T.(11) \boxed{C_-^\dagger P_EC_+ =\left(C_-^\dagger P_{-E}C_+\right)^T.} \tag{11}

The negative sideband has the corresponding conjugate relation. This criterion follows by expanding both resolvents in spectral projectors. It checks an identity at every pole rather than inferring cancellation from a few sampled frequencies. A continuum statement would require the corresponding equality of spectral measures and valid operator domains.

For the centered Cartesian radial examples at depths three, five and six, all complete projectors satisfy (11) for the scalar second moments and dipoles. The direct source response has no ω±2M\omega\pm2M scalar output at rest, even though its field perturbation contains both carriers and its full field operator is periodic. This is a verified property of the declared regulators; it is not an assumption about an arbitrary source, an arbitrary nonlocal probe or a physical photon current.

For identical rest polarizations, write CC for the unit-profile source columns and χ1\chi_1 for their bare single-frequency response. Since C+=aCC_+=aC and C−=−bαyCC_-=-b\alpha_yC, the diagonal harmonic obeys the exact normal-partner identity

Π0(ω,q)=∣a∣2χ1(ω,q)−∣b∣2χ1(ω,q)T.(12) \boxed{\Pi_0(\omega,q) =|a|^2\chi_1(\omega,q)-|b|^2\chi_1(\omega,q)^T.} \tag{12}

This identity follows directly from (7), including the opposite normal momentum. Different occupied rest polarizations require their respective single-frequency responses. Equation (12) concerns the diagonal harmonic; at nonzero qq it does not remove the other blocks of (9).

In particular, if a shape channel has χ1=A(ω)I+B(ω)T\chi_1=A(\omega)I+B(\omega)T, with TT=−TT^T=-T, then

Π0=sA(ω)I+nB(ω)T,n=∣a∣2+∣b∣2.(13) \boxed{\Pi_0=sA(\omega)I+nB(\omega)T, \qquad n=|a|^2+|b|^2.} \tag{13}

The signed scalar balance and canonical norm multiply different parts of the response. Replacing both by a single small scalar source would lose its polarization-dependent dynamics.

Check growth in the full periodic field

When (11) removes the rest scalar sidebands, there is an autonomous source realization with four full spatial blocks:

L4=diag⁡(Hb−M,Hb+M,−Hb∗+M,−Hb∗−M),T4=(C+†,C−†,C+T,C−T),S4=Σ4T4†,Σ4=diag⁡(I,I,−I,−I),A4=L4+S4κT4.(14) \begin{gathered} L_4=\operatorname{diag}(H_b-M,H_b+M,-H_b^*+M,-H_b^*-M),\\ T_4=(C_+^\dagger,C_-^\dagger,C_+^T,C_-^T),\qquad S_4=\Sigma_4T_4^\dagger,\qquad \Sigma_4=\operatorname{diag}(I,I,-I,-I),\\ \boxed{\mathcal A_4=L_4+S_4\kappa T_4.} \end{gathered} \tag{14}

Its resolvent gives exactly Π0\Pi_0. For a mode (u+,u−,v+,v−)(u_+,u_-,v_+,v_-) of frequency ω\omega, reconstruct the physical ordinary and conjugate fields as

Z(t)=e−iωt(e−iMtu++eiMtu−eiMtv++e−iMtv−).(15) \boxed{Z(t)=e^{-i\omega t} \begin{pmatrix} e^{-iMt}u_++e^{iMt}u_-\\ e^{iMt}v_++e^{-iMt}v_- \end{pmatrix}.} \tag{15}

For a nonreal source-carrying pole, (11) makes its unwanted scalar sidebands vanish, so (15) solves (4) at q=0q=0. Source-free modes of L4L_4 have real frequencies and cannot create exponential growth. Conversely, the closed rest moment equation separates the harmonics of a growing Floquet solution into poles of this same source resolvent. This is a finite-dimensional rest-response reduction under (11), not a claim of a global symmetry for arbitrary fields or propagating modes.

The period is P=π/MP=\pi/M. The vector multiplying e−iωte^{-i\omega t} in (15) changes sign after PP, so its full evolution multiplier is −e−iωP-e^{-i\omega P}. Independent integration of (4) over a full period checks this multiplier and the growth rate Im⁡ω\operatorname{Im}\omega. The evolution preserves the indefinite canonical form even when a mode grows; conservation of that form is not positive-energy stability.

At unit number norm, scalar balance s=0.3s=0.3, one common g0=6.5g_0=6.5 and zero bare mass, the total-degree-one Cartesian regulators give:

depth bound frequency MM largest source-pole growth rate full time-evolution growth rate
33 1.7573941.757394 1.0677341.067734 1.0677341.067734
55 2.7921142.792114 0.9887120.988712 0.9887120.988712
66 3.3016313.301631 0.9670620.967062 0.9670620.967062

These values verify the reduction in the declared spatial regulators. They are not converged particle frequencies. In these examples the periodic field is unstable, so a small residue of one slow response pole cannot be interpreted as a stable weakly coupled particle. Slow shape poles and their signed response derives the slow frequency and residue together and verifies their canonical sign.

Propagation and the second object's response remain essential

In the same spatial calculations, normal momentum makes the scalar shape sidebands nonzero. Their first small-qq term is linear in qq; the positive and negative normal momenta give opposite sideband amplitudes for the tested preparation. Thus the rest cancellation does not authorize using a stationary scalar susceptibility for the whole propagating periodic field.

A reciprocal nested interaction uses the curvature of the common energy, including its off-diagonal shape restriction and scalar-center terms. The readjusted nested pair supplies the stationary source requirement. If two periodic backgrounds have different bound frequencies, their common response generally needs both frequency indices; a single period cannot be assumed for an irrational frequency ratio. Both objects' spatial polarizations and their feedback must be retained.

The source-frequency analysis explains why an actual photon output map is still necessary. A scalar moment cancellation is a statement about that readout, while a photon could sample a different spatial current with different surviving sidebands. Spatial current response derives those output columns from the same field and verifies the normal continuity identity before and after source feedback.

Verification

periodic_dirac_response_checks.py constructs complete centered Cartesian field sources with one common interaction coefficient. Complex finite variations verify (4), including all center, breathing and shape terms. An independent temporal Fourier transform reconstructs the full periodic field operator.

Complete spectral projectors check (11). Dense normal Dirac resolvents independently check (6)–(7). Explicit time-dependent forced fields retain both spatial polarizations; their temporal Fourier transforms recover the predicted moment sidebands. Matched full-field and moment harmonic cutoffs check (9)–(10), retaining ordinary and conjugate fluctuations. Full-period evolution independently verifies the reconstructed growing modes, their Floquet multipliers and preservation of the canonical form. Reports use timestamped filenames and checks remain active under Python optimization. These calculations establish the stated algebra and regulated response; they do not identify physical particle masses or interaction rates.