Summary: A smooth quartic action has exactly coherent extrema, an explicit complete mode Hessian and finite-mode orbital stability, with the smooth-contact frequency multiplier.

A specified action on the complete mode band

The cumulative tangent-mode carrier at physical depth dd is

HN=Sym⁡N(Cd+1),DN=(N+dd),N≥1. \mathcal H_N=\operatorname{Sym}^N(\mathbb C^{d+1}), \qquad D_N=\binom{N+d}{d},\qquad N\ge1.

Its generating space consists of dd complexified spatial tangent directions and one homogenizing coefficient. Its projective ray space is an auxiliary mode fibre. It does not add physical spatial dimensions, or replace the actual electron spatial manifold by CP6\mathbb{CP}^6. The self-contact module defines the physical base and mode fibre separately.

At a fixed spatial base point, or after a separately specified spatial profile reduction, retain every coefficient c∈HNc\in\mathcal H_N. For a unit generating vector uu, write ∣u;N⟩=u⊗N|u;N\rangle=u^{\otimes N} and define the normalized coherent transform

sc([u])=DN⟨u;N∣c⟩,∫∣sc∣2 dνFS=∥c∥2=J,(1) s_c([u])=\sqrt{D_N}\langle u;N|c\rangle, \qquad \int |s_c|^2\,d\nu_{\rm FS}=\|c\|^2=J, \tag{1}

where the Fubini--Study measure has total volume one. The phase of this transform depends on the lift of [u][u]; its absolute square is well defined.

Specify the canonical finite-mode action

SN[c]=∫dt[i2(c†c˙−c˙†c)−EN(c)],EN(c)=κ2∫∣sc∣4 dνFS,κ∈R.(2) S_N[c]=\int dt\left[ \frac i2(c^\dagger\dot c-\dot c^\dagger c)-E_N(c)\right], \qquad E_N(c)=\frac{\kappa}{2}\int |s_c|^4\,d\nu_{\rm FS}, \qquad \kappa\in\mathbb R. \tag{2}

The coefficient kinetic metric in (2) is part of this declared action. Equation (1) alone does not prove that a physical spatial Dirac reduction has this metric; its canonical normalization must be derived separately. No omitted spatial kinetic term or spinor component is silently identified with (2).

Exact symmetric product formula

Let T:HN⊗HN→H2NT:\mathcal H_N\otimes\mathcal H_N\to\mathcal H_{2N} be orthogonal symmetrization of all 2N2N tensor indices, with the target identified with its symmetric subspace. It satisfies TT†=ITT^\dagger=I and T†T=Π2NT^\dagger T=\Pi_{2N}, the projector onto that subspace. Put

λN=DN2D2N,Q(c)=T(c⊗c),Ach=T(h⊗c). \lambda_N=\frac{D_N^2}{D_{2N}},\qquad Q(c)=T(c\otimes c),\qquad A_c h=T(h\otimes c).

Applying the coherent resolution of identity at degree 2N2N gives

∫∣sc∣4 dνFS=λN∥Q(c)∥2,ic˙=κλNAc†Q(c).(3) \boxed{\int |s_c|^4\,d\nu_{\rm FS} =\lambda_N\|Q(c)\|^2,\qquad i\dot c=\kappa\lambda_N A_c^\dagger Q(c).} \tag{3}

Indeed, the product of the two degree-NN coherent transforms is (DN/D2N)sQ(c)(D_N/\sqrt{D_{2N}})s_{Q(c)}, and δQ=2Acδc\delta Q=2A_c\delta c. Varying (2) with respect to c‾\overline c produces the second identity, including its coefficient. It is the complete cubic vector field on HN\mathcal H_N, not its evaluation on a frozen ray.

For an orthonormal occupation basis ∣α⟩|\alpha\rangle, ∣α∣=N|\alpha|=N, the map can be computed without a tensor of four mode indices:

T(∣α⟩⊗∣β⟩)=tαβ∣α+β⟩,tαβ=(Nα)(Nβ)(2Nα+β).(4) T(|\alpha\rangle\otimes|\beta\rangle) =t_{\alpha\beta}|\alpha+\beta\rangle, \qquad t_{\alpha\beta} =\sqrt{\frac{\binom N\alpha\binom N\beta} {\binom{2N}{\alpha+\beta}}}. \tag{4}

The multi-index Vandermonde identity makes the sum of tαβ2t_{\alpha\beta}^2 at any fixed α+β\alpha+\beta equal to one. This checks the coisometry normalization independently.

Coherent states are precisely the sharp extrema

Orthogonal projection cannot increase a norm. Consequently

∫∣sc∣4 dνFS≤λN∥c∥4.(5) \boxed{\int|s_c|^4\,d\nu_{\rm FS} \le\lambda_N\|c\|^4.} \tag{5}

For nonzero cc, equality holds exactly when cc is a coherent tensor, c=Jeiϕz⊗Nc=\sqrt J e^{i\phi}z^{\otimes N} with ∥z∥=1\|z\|=1. Here is a direct proof of the equality condition. Normalize cc and suppose c⊗cc\otimes c lies in the fully symmetric 2N2N-index subspace. Exchanging one index between the two copies then has expectation one. That expectation is tr⁡(ρ12)\operatorname{tr}(\rho_1^2), where ρ1\rho_1 is the one-index contraction of ∣c⟩⟨c∣|c\rangle\langle c|. Since ρ1\rho_1 is positive with trace one, it must have rank one. A tensor with a pure one-index contraction factors across that index; symmetry forces all remaining indices to carry the same vector. The converse follows immediately from z⊗N⊗z⊗N=z⊗2Nz^{\otimes N}\otimes z^{\otimes N}=z^{\otimes2N}. For N=1N=1 every vector is already coherent.

This is the fourth-moment instance of coherent-state concentration. The broader convex-function inequality is due to Lieb and Solovej, with equality and quantitative results developed by Nicola, Riccardi and Tilli. The elementary projector proof above supplies the particular result needed for the action without an entropy argument.

At fixed JJ, the coherent orbit is therefore the complete set of global energy maxima for κ>0\kappa>0, and global minima for κ<0\kappa<0. This is variational selection in the specified finite coefficient space; it does not establish the preferred sign of a physical scalar interaction.

Exact coherent branch and nonlinear frequency

For c0=Jz⊗Nc_0=\sqrt J z^{\otimes N}, (3) gives

c(t)=e−iΩNtc0,ΩN=κJλN=κJDN2D2N.(6) \boxed{c(t)=e^{-i\Omega_N t}c_0,\qquad \Omega_N=\kappa J\lambda_N =\kappa J\frac{D_N^2}{D_{2N}}.} \tag{6}

Every coefficient equation is satisfied. The ray [z][z] is stationary and only its common phase evolves. Coherence therefore supplies an exact stationary branch of the smooth action, without replacing its spectrum by that of a point insertion.

The norm JJ is conserved by the phase symmetry. All mode-frame rotation moments are also conserved. In tensor notation,

ρ1=1N(⟨c,aj†aic⟩)ij,tr⁡ρ1=J,ρ˙1=0.(7) \rho_1=\frac1N\bigl(\langle c,a_j^\dagger a_i c\rangle\bigr)_{ij}, \qquad \operatorname{tr}\rho_1=J, \qquad \dot\rho_1=0. \tag{7}

The operators aj†aia_j^\dagger a_i act on the symmetric occupation indices. Equation (7) follows from the U(d+1)U(d+1) invariance of (2). This mode-fibre symmetry is an assumption of the isotropic completion, not an additional set of physical spatial rotations. In particular, conservative evolution does not generally capture an arbitrary pattern into a coherent minimum: its energy and the one-index contraction are fixed.

Complete perturbation equation

Choose z=e0z=e_0 by mode-frame covariance and use a frame rotating with ΩN\Omega_N. For c=e−iΩNt(c0+η)c=e^{-i\Omega_Nt}(c_0+\eta), define the normalized Toeplitz eigenvalue ratio

rj=N(N−1)⋯(N−j+1)2N(2N−1)⋯(2N−j+1),r0=1,0≤j≤N.(8) r_j=\frac{N(N-1)\cdots(N-j+1)} {2N(2N-1)\cdots(2N-j+1)},\qquad r_0=1, \qquad 0\le j\le N. \tag{8}

Here jj counts occupations outside coordinate zero. Its complex multiplicity is (j+d−1d−1)\binom{j+d-1}{d-1}. The full linearized equation is

iη˙0=κJλN(η0+η‾0),iη˙j=κJλN(2rj−1)ηj,1≤j≤N.(9) \boxed{ \begin{aligned} i\dot\eta_0&=\kappa J\lambda_N(\eta_0+\overline\eta_0),\\ i\dot\eta_j&=\kappa J\lambda_N(2r_j-1)\eta_j, \qquad 1\le j\le N. \end{aligned}} \tag{9}

To obtain it directly, write VV for multiplication by ∣sc0∣2|s_{c_0}|^2 followed by projection. Its eigenvalue in sector jj is JλNrjJ\lambda_Nr_j, as follows from the exact FS monomial integral. The variation of ∣s∣2s|s|^2s has both terms

δ(∣s∣2s)=2∣s∣2δs+s2δs‾. \delta(|s|^2s)=2|s|^2\delta s+s^2\delta\overline s.

The second term has a nonzero projected matrix entry only when the two occupation multi-indices sum to 2Ne02Ne_0, so only the radial entry survives. Subtracting the rotating-frame frequency gives (9). Dropping this conjugate term would give an incorrect radial equation.

Since r1=1/2r_1=1/2, the dd complex tangent directions to the coherent ray are zero modes. For j≥2j\ge2, 0<rj<1/20<r_j<1/2, and all remaining frequencies are real and have the same sign opposite to κ\kappa. When N≥2N\ge2, the smallest nonzero absolute frequency is

∣Ωtransverse,min∣=∣κ∣JλNN2N−1.(10) |\Omega_{\rm transverse,min}| =|\kappa|J\lambda_N\frac{N}{2N-1}. \tag{10}

The radial equation contains the phase zero mode and the phase drift caused by changing the conserved norm. On a fixed-norm space modulo phase this does not supply an exponentially growing perturbation.

Finite-mode nonlinear orbital stability

At fixed nonzero JJ, the state sphere is compact. For N≥2N\ge2 the nonnegative deficit

D(c)=λNJ2−∫∣sc∣4 dνFS \mathcal D(c)=\lambda_NJ^2-\int|s_c|^4\,d\nu_{\rm FS}

vanishes exactly on the compact coherent orbit. It has a strictly positive minimum on every closed subset a positive distance from that orbit. Norm and energy conservation therefore imply nonlinear orbital stability of the coherent set under (3), for either nonzero sign of κ\kappa. The polynomial vector field exists globally because the norm remains bounded. This argument supplies finite-mode nonlinear stability, rather than inferring it solely from real linearized frequencies.

The conclusion concerns the whole coherent set, including phase and ray rotations. It is not attraction to one prepared ray. Physical spatial variation, omitted bands, Dirac components, nuclear sources and additional interaction terms require their own stability analysis.

Consequence for the mass construction

For unit coefficient norm and a band-independent κ\kappa, the derived frequency multiplier is DN2/D2ND_N^2/D_{2N}. The occupied point-contact multiplier is DND_N. Their ratio is

λNDN=DND2N=∏r=1dN+r2N+r, \frac{\lambda_N}{D_N} =\frac{D_N}{D_{2N}} =\prod_{r=1}^d\frac{N+r}{2N+r},

which depends on NN and tends to 2−d2^{-d} for large NN at fixed dd. At the specified electron and muon bands (12,34)(12,34) in depth six, (6) gives the smooth-action ratio 137.029686760779137.029686760779, while the point count gives 206.764705882353206.764705882353. These are the distinct prescriptions compared in the contact-kernel audit.

Thus coherent extrema and finite-mode stability can be established for a specific action while its mass-sequence multiplier remains different. Converting (6) into a physical rest mass additionally requires the spatial Dirac equation, its normalized profile and the identification of its pole.

The homogeneous coherent action uses a normalized square root of the quartic current integral to obtain the occupied coefficient DND_N. Its full variation has a different response, with a cancellation of the radial conjugate term. It supplies exact nonconstant spatial solutions for the specified independent bundle, while the changed functional and its normalization require a master origin. This does not alter the quartic action's result (6).

Reciprocal extension

The complete pair action retains both full mode families. A generic pair leaves the coherent subset, and conservation relates its shape deformation to changing density contact. The aligned stationary set has a separate complete stability calculation.

The current projection identity now expresses this same smooth quartic as an exact weighted sum of all reduced tensor-current purities. It identifies the coefficients supplied by the projection, including signed full-spinor currents and different grades. This gives a direct operator comparison with the shell-weighted self-mass action, rather than inferring that action from the coherent count alone.

coherent_contact_variation_checks.py checks independent FS quadrature, the full nonlinear gradient and both Hessian blocks, symmetric-product normalization, the coherent stationary equation, all perturbation sectors and complete finite-mode evolution. Exact rational formulas also cover large bands without allocating their matrices. The supporting algebra is in coherent_mode_algebra.py.