Summary: The written scalar Dirac self-interaction gives an explicit real-linear response with ordinary and conjugate mode couplings determined by the higher spatial field.

The existing nonlinear term

Use the scalar self-interaction present in the projected Dirac action. In units ℏ=1\hbar=1, its Hamiltonian form is

i∂tΨ=H0Ψ+N[Ψ],N[Ψ]=κ(Ψ†βΨ)βΨ.(1) i\partial_t\Psi=H_0\Psi+\mathcal N[\Psi], \qquad \mathcal N[\Psi]=\kappa(\Psi^\dagger\beta\Psi)\beta\Psi. \tag{1}

Here β=β†\beta=\beta^\dagger, β2=I\beta^2=I, κ\kappa is real, and H0H_0 is the specified self-adjoint spatial Dirac operator including whichever fixed background terms have actually been adopted. Its domain and the spatial measure must make the following pairings valid. The energy is

E[Ψ]=⟨Ψ,H0Ψ⟩+κ2∫(Ψ†βΨ)2dμ.(2) \mathcal E[\Psi] =\langle\Psi,H_0\Psi\rangle +\frac\kappa2\int(\Psi^\dagger\beta\Psi)^2d\mu. \tag{2}

This is one field on the nested macroscopic spatial manifold. Spinor components describe amplitudes over that space, not independent spatial dimensions. Equation (1) is a specified contribution to V2 dynamics, not the completed master action or a derivation of occupied particle geometries. No mediator particle is introduced.

Nonlinear Dirac equations with this scalar structure and linearization about solitary waves are treated in Boussaïd and Comech. The sign convention here is fixed directly by (2).

Differentiate the complete source

At a spatial field ϕ\phi, set

s=ϕ†βϕ∈R,v=βϕ. s=\phi^\dagger\beta\phi\in\mathbb R, \qquad v=\beta\phi.

For a variation η\eta, the scalar changes by δs=v†η+η†v\delta s=v^\dagger\eta+\eta^\dagger v. Therefore

DN[ϕ]η=κ[sβη+v(v†η+η†v)].(3) D\mathcal N[\phi]\eta =\kappa\left[s\beta\eta +v(v^\dagger\eta+\eta^\dagger v)\right]. \tag{3}

Equivalently,

DN[ϕ]η=ANη+BNη‾,AN=κ(sβ+vv†),BN=κvvT.(4) \boxed{ D\mathcal N[\phi]\eta=A_N\eta+B_N\overline\eta, \qquad A_N=\kappa(s\beta+vv^\dagger), \qquad B_N=\kappa vv^T.} \tag{4}

The sβs\beta term is the variation of the final field factor. The two rank-one terms come from changing the source scalar itself. Retaining only sβηs\beta\eta would freeze the higher particle's source while it is being perturbed. The complete derivative is real-linear; generally DN[iη]≠iDN[η]D\mathcal N[i\eta]\ne iD\mathcal N[\eta].

The matrices obey AN=AN†A_N=A_N^\dagger and BN=BNTB_N=B_N^T. Spatial variation of s(x)s(x) and v(x)v(x) retains the target's local shape and polarization. Replacing them by spatial averages is an additional approximation.

Linearized dynamics about a stationary profile

Suppose Ψ0=e−iμtϕ\Psi_0=e^{-i\mu t}\phi is an actual solution of (1), so

(H0+κsβ)ϕ=μϕ.(5) (H_0+\kappa s\beta)\phi=\mu\phi. \tag{5}

In the same rotating frame put Ψ=e−iμt(ϕ+η)\Psi=e^{-i\mu t}(\phi+\eta). The linearized equation is

i∂tη=Aη+Bη‾,A=H0−μ+AN,B=BN.(6) i\partial_t\eta=A\eta+B\overline\eta, \qquad A=H_0-\mu+A_N, \qquad B=B_N. \tag{6}

Its doubled representation is

i∂t(ηη‾)=L(ηη‾),L=(AB−B‾−A‾).(7) i\partial_t \begin{pmatrix}\eta\\\overline\eta\end{pmatrix} =\mathcal L \begin{pmatrix}\eta\\\overline\eta\end{pmatrix}, \qquad \mathcal L= \begin{pmatrix}A&B\\-\overline B&-\overline A\end{pmatrix}. \tag{7}

The second row is the conjugate of the first; it does not add another physical field or another set of spatial dimensions. With Σ3=diag⁡(I,−I)\Sigma_3=\operatorname{diag}(I,-I),

H=Σ3L=(ABB‾A‾)=H†.(8) \mathcal H=\Sigma_3\mathcal L =\begin{pmatrix}A&B\\\overline B&\overline A\end{pmatrix} =\mathcal H^\dagger. \tag{8}

This is the quadratic-energy structure. L\mathcal L itself need not be Hermitian or have a real spectrum. Stability cannot be inferred merely from the Hermiticity of H0H_0. The global phase tangent supplies the consistency identity

L(iϕ−iϕ‾)=0,(9) \mathcal L \begin{pmatrix}i\phi\\-i\overline\phi\end{pmatrix}=0, \tag{9}

which follows immediately from (5). A moving orbit or a background with several temporal frequencies instead requires time-dependent coefficients; a stationary spectral analysis cannot silently replace that motion.

Spatial coupling between mode families

Let fa(x)f_a(x) be photon-sector test modes and uλ(x)u_\lambda(x) target variation modes, normalized by the action's kinetic metric. The ordinary and conjugate matrix elements of the interaction derivative are

Kaλ=∫fa†ANuλ dμ,Caλ=∫fa†BNuλ‾ dμ.(10) K_{a\lambda}=\int f_a^\dagger A_Nu_\lambda\,d\mu, \qquad C_{a\lambda}=\int f_a^\dagger B_N\overline{u_\lambda}\,d\mu. \tag{10}

Thus a target amplitude vector qq produces the linearized source Kq+Cq‾Kq+C\overline q. The spatial integrals include polarization, relative phase, and all of the higher-dimensional profile present in (4). Any cross-family kinetic matrix element of H0H_0 must also be retained.

Interchanging target and photon mode labels gives Hermitian reciprocity for KK and transpose symmetry for CC:

Kλa=Kaλ‾,Cλa=Caλ.(11) K_{\lambda a}=\overline{K_{a\lambda}}, \qquad C_{\lambda a}=C_{a\lambda}. \tag{11}

For the second equality, use BNT=BNB_N^T=B_N and transpose the scalar pairing. In the doubled mode space these are precisely the symmetries of the real quadratic energy. They do not require treating the conjugate variation as an independent particle.

A zero of one matrix element excludes that particular linear response. It is not a zero of the complete interaction unless all relevant ordinary, conjugate, and kinetic blocks vanish. The criterion must be applied to the actual nonlinear source, not to the parity of a bare field alone; angular kernel parity gives an explicit example of this distinction.

Which temporal frequencies the source contains

For a finite prescribed field h(x,t)=∑jcjϕj(x)e−iΩjth(x,t)=\sum_j c_j\phi_j(x)e^{-i\Omega_jt}, equation (1) gives

N[h]=κ∑j,ℓ,mc‾jcℓcm(ϕj†βϕℓ)βϕme−i(Ωℓ+Ωm−Ωj)t.(12) \mathcal N[h] =\kappa\sum_{j,\ell,m} \overline c_jc_\ell c_m (\phi_j^\dagger\beta\phi_\ell)\beta\phi_m e^{-i(\Omega_\ell+\Omega_m-\Omega_j)t}. \tag{12}

Projection onto an output mode faf_a uses the full cubic spatial overlap

Ta;jℓm=κ∫fa†βϕm(ϕj†βϕℓ) dμ.(13) T_{a;j\ell m} =\kappa\int f_a^\dagger\beta\phi_m (\phi_j^\dagger\beta\phi_\ell)\,d\mu. \tag{13}

The frequency combination and spatial compatibility follow from the same written term. Frequency matching alone does not make (13) nonzero. Several triples reaching the same output frequency interfere coherently.

Equations (10)–(13) assume a linear output-mode inclusion in the master field. If a photon is instead realized as a nonlinear collective observable of that field, derive its tangent map and induced equations. For example a bilinear observable contains frequency differences, whereas the cubic field source in (12) contains the displayed three-term combinations. Their different frequency rules cannot be interchanged without deriving that output map.

Scope of the concrete calculation

Complete scalar response determines the full pointwise rank, singular values and inverse of (3). At nonzero scalar density the complete source derivative is invertible; at scalar nullity it retains one real response direction. Those facts concern the full source before a spatial channel projection and before adding H0−μH_0-\mu.

Transition-density coupling derives the corresponding nonlocal source derivative and full quartic mode tensor for the projected bilocal density interaction, with both patterns retained.

The operator DN[ϕ]D\mathcal N[\phi] is now explicit for this term. The remaining physical inputs are the actual higher-sector solution, the photon embedding or collective map, all additional action terms, and the boundary conditions selecting incoming and outgoing modes. The matrix formulas apply at any depth supporting the specified Dirac structure and do not themselves single out depth five.

Local scalar contact normalization adds the spatial virial and an exact normalized massive profile for the same local source convention. It makes explicit why a scalar vertex coefficient does not by itself determine a stationary particle frequency.

Verification

Local polarized pulse coupling retains the exact cross-current of two traveling components of this same field. Its complete local equations give source zeros and both reciprocal weak-scattering maps without first postulating a bilocal spatial kernel.

Causal scalar response solves this complete derivative about a specified uniform spatial background, including every scalar angular mode, both signed propagation channels and the conserved zero mode. The background is a control, not an identified particle.

polarization_coupling_checks.py checks the displayed identities with independent spatial quadrature, energy variations, or full finite-system evolution as appropriate. These checks verify the declared mathematics and benchmarks, not empirical couplings.

Source frequencies in spatial coupling tracks the distinct temporal frequencies in the bilinear density, complete cubic field source and perturbations relative to the stationary carrier.