Summary: The critical conformal scalar operator on a round even-dimensional sphere has exactly the normalized simplex spectrum, while its spinor realization, spatial geometry and master interaction remain additional requirements.

An exact geometric occurrence of the mass-count polynomial

There is a known natural differential operator whose normalized eigenvalues are exactly

S(n,d)=(n+d−1d).(1) S(n,d)=\binom{n+d-1}{d}. \tag{1}

For even dd, it is the critical conformal scalar operator on a round dd-sphere. The result follows from spatial conformal covariance and can also be verified by direct differentiation. It gives a distinct mathematical origin for the same sequence used in the mass-spectrum construction.

This calculation specifies a round physical spatial geometry SRdS_R^d, with undetermined radius RR. Its directions are ordinary macroscopic spatial directions; there is no compactification or hidden small radius. Spherical harmonic coefficients do not add dimensions.

The operator acts on scalar spatial functions. Its appearance in the master spinor interaction, and the identification of its harmonic degree with an occupied particle mode, must be derived separately.

Conformal covariance fixes the recurrence

Work first at unit radius. Let PP commute with round-sphere rotations, and let its scalar eigenvalue on the complete degree-ℓ\ell harmonic space be pℓp_\ell. Rotation covariance alone leaves those eigenvalues free. Impose the stronger critical conformal law

P(C∗f)=ΩdC∗(Pf),C∗g=Ω2g.(2) P(C^*f)=\Omega^d C^*(Pf),\qquad C^*g=\Omega^2g. \tag{2}

This is the scalar law with input conformal weight zero and output weight −d-d. Its infinitesimal form, for LVg=2ωg\mathcal L_Vg=2\omega g, is

[P,V]f=dωPf.(3) [P,V]f=d\omega Pf. \tag{3}

Choose the height coordinate t=Xd+1t=X_{d+1} on the unit sphere and V=∇tV=\nabla t. Since ∇2t=−tg\nabla^2t=-tg, its conformal factor is ω=−t\omega=-t.

For ℓ≥1\ell\ge1, the complex harmonic Hℓ=(X1+iX2)ℓH_\ell=(X_1+iX_2)^\ell satisfies

VHℓ=−ℓ tHℓ.(4) VH_\ell=-\ell\,tH_\ell. \tag{4}

The ambient polynomial tHℓtH_\ell is itself harmonic and homogeneous of degree ℓ+1\ell+1, because HℓH_\ell is independent of tt. Applying (3) to (4) gives

−ℓ(pℓ+1−pℓ)tHℓ=−d pℓtHℓ. -\ell(p_{\ell+1}-p_\ell)tH_\ell =-d\,p_\ell tH_\ell.

Thus

pℓ+1pℓ=ℓ+dℓ(ℓ≥1, p1≠0).(5) \boxed{\frac{p_{\ell+1}}{p_\ell}=\frac{\ell+d}{\ell}} \qquad(\ell\ge1,\ p_1\ne0). \tag{5}

Applying (3) to the constant function also gives p0=0p_0=0. Iteration yields

pℓp1=∏j=1ℓ−1j+dj=(ℓ+d−1d).(6) \boxed{ \frac{p_\ell}{p_1} =\prod_{j=1}^{\ell-1}\frac{j+d}{j} =\binom{\ell+d-1}{d}. } \tag{6}

The last expression is exactly the cumulative Hockey-Stick count in (1). Here it is an eigenvalue ratio forced by the specified conformal law, rather than a count being substituted for an eigenvalue. The chosen complex harmonic is a convenient way to derive the recurrence; rotation covariance gives the same value on every polarization of its harmonic space.

A local differential operator realizes that law

Let Δ=−∇a∇a\Delta=-\nabla^a\nabla_a be the nonnegative scalar Laplacian, and write d=2hd=2h. The critical GJMS operator has order dd; critical means that its differential order equals the spatial dimension. On the round sphere it is

Pd,R=∏j=1h[Δ+(h−j)(h+j−1)R2].(7) \boxed{ P_{d,R}= \prod_{j=1}^{h} \left[\Delta+\frac{(h-j)(h+j-1)}{R^2}\right]. } \tag{7}

The factors and the conformal transformation law are the standard Branson factorization, derived explicitly through stereographic projection in Graham, Proposition 1. The leading term in (7) is Δd/2\Delta^{d/2}.

On a degree-ℓ\ell scalar harmonic, R2Δ=ℓ(ℓ+d−1)R^2\Delta=\ell(\ell+d-1). Each factor becomes

ℓ(ℓ+d−1)+(h−j)(h+j−1)=(ℓ+h−j)(ℓ+h+j−1).(8) \ell(\ell+d-1)+(h-j)(h+j-1) =(\ell+h-j)(\ell+h+j-1). \tag{8}

The first factors cover ℓ,…,ℓ+h−1\ell,\ldots,\ell+h-1 and the second cover ℓ+h,…,ℓ+d−1\ell+h,\ldots,\ell+d-1. Therefore

Pd,RYℓ,a=ℓ(ℓ+1)⋯(ℓ+d−1)RdYℓ,a=d!RdS(ℓ,d)Yℓ,a.(9) \boxed{ P_{d,R}Y_{\ell,a} =\frac{\ell(\ell+1)\cdots(\ell+d-1)}{R^d}Y_{\ell,a} =\frac{d!}{R^d}S(\ell,d)Y_{\ell,a}. } \tag{9}

The label aa runs over the complete scalar harmonic multiplet. Its multiplicity is (ℓ+dd)−(ℓ+d−2d)\binom{\ell+d}{d}-\binom{\ell+d-2}{d}, with the second term zero when ℓ<2\ell<2. That multiplicity is separate from the eigenvalue (9). The constant harmonic has eigenvalue zero.

For the two charged-lepton depths, the operators are particularly explicit:

P6,R=Δ(Δ+4R2)(Δ+6R2),P10,R=Δ(Δ+8R2)(Δ+14R2)(Δ+18R2)(Δ+20R2).(10) \begin{aligned} P_{6,R} &=\Delta\left(\Delta+\frac4{R^2}\right) \left(\Delta+\frac6{R^2}\right),\\ P_{10,R} &=\Delta\left(\Delta+\frac8{R^2}\right) \left(\Delta+\frac{14}{R^2}\right) \left(\Delta+\frac{18}{R^2}\right) \left(\Delta+\frac{20}{R^2}\right). \end{aligned} \tag{10}

Define the dimensionless operator P^d=RdPd,R/d!\widehat P_d=R^dP_{d,R}/d!, which normalizes the first positive harmonic eigenvalue to one. At the stated integer indices:

Spatial dimension dd Scalar harmonic degree ℓ\ell Eigenvalue of P^d\widehat P_d
66 1313 18 56418\,564
66 3535 3 838 3803\,838\,380
1010 2323 64 512 24064\,512\,240

The degrees in the table have not been selected by (2) or (7). Identifying ℓ=n\ell=n with a particle's occupied mode is an additional physical step. The harmonic degree ℓ\ell is also different from the oscillator cutoff N=n−1N=n-1. No change of particle dimensionalities or unrestricted mode search is being used.

At the level of scalar spatial operators, a trial family Md=μdP^dM_d=\mu_d\widehat P_d has precisely the proposed eigenvalue sequence. It is one differential operator for all harmonic degrees, with a common scale in each dimension. Its cumulative harmonic compressions obey the ordinary consistency law. This is not the auxiliary rank-one contact family.

Normalizing by d!d! is an explicit spectral convention. It does not derive the dimensional coefficient μd\mu_d, nor the charged-lepton relation μ10=μ6\mu_{10}=\mu_6. A master reduction must fix those physical normalizations across sectors.

Verify it by differentiating spatial functions

In stereographic coordinates, put s=∣x∣2s=|x|^2 and Ω=2/(1+s)\Omega=2/(1+s). At unit radius, the critical conformal law gives

Pdf=Ω−d(−ΔRd)d/2f.(11) P_d f=\Omega^{-d}(-\Delta_{\mathbb R^d})^{d/2}f. \tag{11}

For a zonal harmonic,

f(s)=Cℓ(d−1)/2(1−s1+s),ΔRdf=4sf′′+2df′,(12) f(s)=C_\ell^{(d-1)/2}\left(\frac{1-s}{1+s}\right),\qquad \Delta_{\mathbb R^d}f=4s f''+2d f', \tag{12}

where CℓαC_\ell^\alpha is a Gegenbauer polynomial. Exact rational differentiation of (11)–(12) gives (9), including the assigned degrees.

A second family retains a complex angular phase. With the homogeneous Euclidean harmonic Hℓ=(x1+ix2)ℓH_\ell=(x_1+ix_2)^\ell,

f(x)=2ℓHℓ(x)(1+s)−ℓ,ΔRd(Hℓh(s))=Hℓ[4sh′′+(2d+4ℓ)h′].(13) f(x)=2^\ell H_\ell(x)(1+s)^{-\ell},\qquad \Delta_{\mathbb R^d}(H_\ell h(s)) =H_\ell\left[4s h''+(2d+4\ell)h'\right]. \tag{13}

It gives the same eigenvalue independently of the zonal calculation. Thus the spectral identity is supported by both the conformal recurrence and the actual spatial differential action.

Scalar covariance is not automatically spinor covariance

The standard untwisted spinor Dirac spectrum on a round sphere is ±(k+d/2)/R\pm(k+d/2)/R, with k=0,1,…k=0,1,\ldots, as discussed in round-sphere Dirac counts. Its conformal differential powers have a different factorization: at unit radius they include

D∏j=1m(D2−j2).(16) \mathscr D\prod_{j=1}^{m}(\mathscr D^2-j^2). \tag{16}

The spinorial intertwining family has multipliers proportional to

Γ(d/2+k+1/2−λ)Γ(d/2+k+1/2+λ).(17) \frac{\Gamma(d/2+k+1/2-\lambda)} {\Gamma(d/2+k+1/2+\lambda)}. \tag{17}

For even dimension the intertwiner maps between the appropriate chiral bundles; it is not the scalar operator acting on independent component functions. These conventions and (16) are given in Clerc and Ørsted, Proposition 7.1 and Theorem 7.1.

At the parameter λ=−d/2\lambda=-d/2, the magnitude of (17) becomes

(k+12)(k+32)⋯(k+d−12).(18) \left(k+\frac12\right) \left(k+\frac32\right)\cdots \left(k+d-\frac12\right). \tag{18}

Its consecutive ratio is (k+d+1/2)/(k+1/2)(k+d+1/2)/(k+1/2), rather than (5) with integer ℓ\ell. An overall normalization does not remove this spectral shift. Consequently the successful scalar construction cannot simply be declared to be the standard untwisted spinor mass operator. A different physical twist or composite scalar source requires its own derivation.

An exact transfer through a complete Killing-spinor frame does produce a local spinor operator with the scalar spectrum. Its opposite Dirac signs carry adjacent scalar degrees, and its real scalar mass insertion changes the coefficient. This is a distinct operator construction, with its action and physical scope checked separately.

Acting on the scalar current changes the nonlinear equation

Another possible insertion would let PdP_d act on the real scalar current J=ψ†βψJ=\psi^\dagger\beta\psi. On a closed spatial manifold, consider the explicitly different interaction

EP(ψ)=η2∫JPdJ dμ.(19) \mathcal E_P(\psi)=\frac{\eta}{2}\int J P_d J\,d\mu. \tag{19}

Self-adjointness gives

δEP=η∫(PdJ) δJ dμ,δEPδψ†=η(PdJ)βψ.(20) \delta\mathcal E_P=\eta\int(P_dJ)\,\delta J\,d\mu,\qquad \frac{\delta\mathcal E_P}{\delta\psi^\dagger} =\eta(P_dJ)\beta\psi. \tag{20}

The variation is multiplication by PdJP_dJ; it is not PdψP_d\psi. Moreover,

Pd1=0,∫PdJ dμ=0.(21) P_d1=0,\qquad \int P_dJ\,d\mu=0. \tag{21}

Equation (19) contains derivatives of the current. It is an additional interaction, beyond the written algebraic scalar term tested by the local derivative calculation. Any claim that it arises by eliminating modes must identify those modes of the same master field and retain the induced kinetic normalization.

The complete response on a zero-current Killing background provides one exact spatial test of this proposed interaction. Its simplex coefficient enters collective frequencies; any nonzero derivative-current coefficient gives unstable high harmonics on that specified branch.

The exact result is therefore a geometric origin for the simplex scalar spectrum, with conformal covariance fixing the integer recurrence. The physical mass-operator problem still requires the actual spatial geometry, a scalar or spinor realization inherited from the master field, a selected harmonic mode, and a dynamical realization whose frequency matches the known mass-frequency and whose inertial response agrees with it. Converting mass to frequency units is already available; it is not a separate unresolved derivation.

The surface-bending construction gives a different geometric occurrence of the same sequence: it is the Hessian of an extrinsic curvature energy on a surface of dimension d−2d-2 in ordinary dd-dimensional space. Its complete normal pair and nonlinear curvature variation provide a separate dynamical test from the scalar operator here.

Verification

conformal_simplex_spectrum_checks.py checks the exact factorizations, the complete conformal commutator on zonal harmonics, the normalized simplex sequence, and independent stereographic differentiation of both zonal and complex angular fields. Exact rational arithmetic covers the assigned degrees without using a fitted spectral polynomial.

Those operator identities do not select a physical geometry or prove that the master action contains this conformal interaction.

A Distinct Flat-Ball Realization

The flat-ball count clock is the restricted subcritical order-d operator on an auxiliary sphere of dimension d+2, with a nonzero constant clock and its own polynomial multiplicities. The differential bridge passes through the critical operator here and derives its exact zonal kinetic normalization. The auxiliary sphere is a mathematical representation, not another physical sector manifold.