Solar system and laboratory tests of the spectral action scale

Abstract

We derive the observational consequences of the spectral action S=Trf(D2/Λ2)S=\mathrm{Tr}\,f(D^{2}/\Lambda^{2}) for weak-field gravity, using the one-loop form factors computed with the full Standard Model content. The linearized field equations yield a modified Newtonian potential with two Yukawa corrections whose amplitudes (4/3-4/3 and +1/3+1/3) are fixed entirely by the spin decomposition of the graviton propagator, and whose ranges are set by calculable effective masses m2=Λ60/13m_{2}=\Lambda\sqrt{60/13} (spin-2) and m0=Λ/6(ξ1/6)2m_{0}=\Lambda/\sqrt{6(\xi-1/6)^{2}} (spin-0), where Λ\Lambda is the spectral cutoff and ξ\xi is the Higgs non-minimal coupling. The 1/r1/r singularity of the Newtonian potential is regularized: V(r)/VN(r)0V(r)/V_{\mathrm{N}}(r)\to 0 as r0r\to 0, and Newton’s law is recovered for r1/Λr\gg 1/\Lambda. We extract the post-Newtonian parameter γ(r)\gamma(r) and confront it with seven gravitational experiments spanning torsion-balance, time-delay, Casimir, satellite, and lunar-ranging measurements. The strongest bound arises from the Eöt-Wash experiment, yielding Λ>2.565×103eV\Lambda>2.565\times 10^{-3}\,\mathrm{eV} (1/Λ<77μm1/\Lambda<77\,\mu\mathrm{m}) at 95% CL. All solar-system tests are satisfied with exponential margin. The bound is determined by the SM content at conformal coupling ξ=1/6\xi=1/6 and depends weakly on ξ\xi otherwise.

spectral action, post-Newtonian parameters, modified gravity, short-range gravity, Standard Model

I Introduction

The spectral action principle [Chamseddine:1996zu, Connes:2006qj] generates the bosonic action from the spectrum of the Dirac operator: S=Trf(D2/Λ2)S=\mathrm{Tr}\,f(D^{2}/\Lambda^{2}), where Λ\Lambda is a spectral cutoff scale. Beyond the Einstein–Hilbert term, the one-loop heat kernel expansion generates curvature-squared corrections characterized by nonlocal form factors [Barvinsky:1985an, Codello:2012kq]. These form factors, when computed with the full Standard Model (SM) spectrum (4 real scalars, 45/245/2 Dirac fermions, 12 gauge bosons), yield the SM-determined Weyl-squared coefficient αC=13/120\alpha_{C}=13/120 [Alfyorov:2026ff].

The curvature-squared structure modifies the graviton propagator, leading to a modified Newtonian potential at distances r1/Λr\lesssim 1/\Lambda, in direct analogy with the Stelle theory of quadratic gravity [Stelle:1976gc, Stelle:1978]. Unlike Stelle gravity, where the quadratic couplings are free parameters, the spectral action determines them from the particle content.

In this paper we derive the observational consequences: the post-Newtonian parameter γ(r)\gamma(r), the modified Newtonian potential, and the constraints on Λ\Lambda from solar-system, torsion-balance, Casimir, and satellite experiments. Our central result is that the Eöt-Wash torsion-balance experiment [Lee:2020] sets the strongest bound, Λ>2.565×103eV\Lambda>2.565\times 10^{-3}\,\mathrm{eV} at 95% CL, corresponding to a Yukawa range of 36μm36\,\mu\mathrm{m}. All other experiments are either exponentially suppressed or below the sensitivity threshold.

The paper is organized as follows. Section II derives the modified potential from the linearized field equations. Section III extracts the post-Newtonian parameter γ(r)\gamma(r). Section IV confronts the predictions with experiments. Section VI discusses the results and their implications. We conclude in Sec. VII.

We use the Lorentzian signature (,+,+,+)(-,+,+,+), natural units c==1c=\hbar=1, and restore SI units when comparing with experiment.

II Modified Newtonian Potential

II.1 From the spectral action to linearized field equations

The one-loop spectral action at 𝒪(2)\mathcal{O}(\mathcal{R}^{2}) in the Weyl basis is [Alfyorov:2026ff]

Γ(1)=116π2d4xg[αCF^1(z)CμνρσCμνρσ+αR(ξ)F^2(z,ξ)R2],\Gamma^{(1)}=\frac{1}{16\pi^{2}}\int\mathrm{d}^{4}x\,\sqrt{-g}\,\Bigl{[}\alpha% _{C}\,\hat{F}_{1}(z)\,C_{\mu\nu\rho\sigma}C^{\mu\nu\rho\sigma}\\ +\alpha_{R}(\xi)\,\hat{F}_{2}(z,\xi)\,R^{2}\Bigr{]}, (1)

where z/Λ2z\equiv\Box/\Lambda^{2}, the normalized form factors satisfy F^i(0)=1\hat{F}_{i}(0)=1, and the coefficients are [Alfyorov:2026ff]

αC=13120,αR(ξ)=2(ξ16)2.\alpha_{C}=\frac{13}{120}\,,\qquad\alpha_{R}(\xi)=2\!\left(\xi-\frac{1}{6}% \right)^{\!2}\,. (2)

Linearizing the total action S=SEH+Γ(1)S=S_{\mathrm{EH}}+\Gamma^{(1)} around flat space gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu} and decomposing via Barnes–Rivers spin projectors yields [Stelle:1976gc, vanDam:1970vg] the dressed propagator denominators

ΠTT(z)\displaystyle\Pi_{\mathrm{TT}}(z) =1+1360zF^1(z),\displaystyle=1+\frac{13}{60}\,z\,\hat{F}_{1}(z)\,, (3)
Πs(z,ξ)\displaystyle\Pi_{\mathrm{s}}(z,\xi) =1+6(ξ16)2zF^2(z,ξ),\displaystyle=1+6\!\left(\xi-\frac{1}{6}\right)^{\!2}z\,\hat{F}_{2}(z,\xi)\,, (4)

for the spin-2 (transverse-traceless) and spin-0 (scalar trace) sectors, respectively. Both satisfy Π(0)=1\Pi(0)=1, recovering GR in the infrared. At conformal coupling ξ=1/6\xi=1/6, the scalar sector decouples identically: Πs1\Pi_{\mathrm{s}}\equiv 1.

II.2 Yukawa potentials

For a static point source T00=Mδ(3)(x)T_{00}=M\,\delta^{(3)}(\vec{x}), the linearized field equations in Fourier space yield the temporal and spatial metric potentials

g00=(1+2Φ),gij=(1+2Ψ)δij.g_{00}=-(1+2\Phi)\,,\qquad g_{ij}=(1+2\Psi)\delta_{ij}\,. (5)

The potentials factorize through Newton kernels

KΦ(z)=43ΠTT(z)13Πs(z,ξ),K_{\Phi}(z)=\frac{4}{3\,\Pi_{\mathrm{TT}}(z)}-\frac{1}{3\,\Pi_{\mathrm{s}}(z,% \xi)}\,, (6)
KΨ(z)=23ΠTT(z)+13Πs(z,ξ),K_{\Psi}(z)=\frac{2}{3\,\Pi_{\mathrm{TT}}(z)}+\frac{1}{3\,\Pi_{\mathrm{s}}(z,% \xi)}\,, (7)

with KΦ(0)=KΨ(0)=1K_{\Phi}(0)=K_{\Psi}(0)=1 (Newtonian limit). The coefficients 4/34/3, 2/32/3, and ±1/3\pm 1/3 are fixed by the spin decomposition and are universal for all quadratic-gravity theories [Stelle:1976gc, Edholm:2016].

In the local Yukawa approximation (F^i1\hat{F}_{i}\approx 1, valid for k2Λ2k^{2}\ll\Lambda^{2}), the propagator denominators reduce to ΠTT1+(13/60)z\Pi_{\mathrm{TT}}\approx 1+(13/60)z and Πs1+6(ξ1/6)2z\Pi_{\mathrm{s}}\approx 1+6(\xi-1/6)^{2}z, yielding the partial-fraction decomposition

KΦloc(z)=43m22k2+m2213m02k2+m02,K_{\Phi}^{\mathrm{loc}}(z)=\frac{4}{3}\,\frac{m_{2}^{2}}{k^{2}+m_{2}^{2}}-% \frac{1}{3}\,\frac{m_{0}^{2}}{k^{2}+m_{0}^{2}}\,, (8)

where the effective masses are

m2=Λ60132.148Λ,m_{2}=\Lambda\sqrt{\frac{60}{13}}\approx 2.148\,\Lambda\,, (9)
m0(ξ)=Λ6(ξ1/6)2(ξ1/6).m_{0}(\xi)=\frac{\Lambda}{\sqrt{6(\xi-1/6)^{2}}}\qquad(\xi\neq 1/6)\,. (10)

The inverse Fourier transform, via the Gradshteyn–Ryzhik identity

2π0sin(kr)km2k2+m2dk=1emr,\frac{2}{\pi}\int_{0}^{\infty}\frac{\sin(kr)}{k}\,\frac{m^{2}}{k^{2}+m^{2}}\,% \mathrm{d}k=1-e^{-mr}\,, (11)

gives the central result:

V(r)VN(r)=143em2r+13em0r.\boxed{\frac{V(r)}{V_{\mathrm{N}}(r)}=1-\frac{4}{3}\,e^{-m_{2}r}+\frac{1}{3}\,% e^{-m_{0}r}\,.} (12)

The corresponding spatial potential (from KΨK_{\Psi}) is

Ψ(r)ΨN(r)=123em2r13em0r.\frac{\Psi(r)}{\Psi_{\mathrm{N}}(r)}=1-\frac{2}{3}\,e^{-m_{2}r}-\frac{1}{3}\,e% ^{-m_{0}r}\,. (13)

Note the different coefficients: the spin-2 contribution enters with 2/32/3 (not 4/34/3) and the spin-0 with 1/3-1/3 (not +1/3+1/3), reflecting the different projections of the graviton propagator onto g00g_{00} and gijg_{ij}.

At conformal coupling, the scalar term is absent:

V(r)VN(r)|ξ=1/6=143em2r.\frac{V(r)}{V_{\mathrm{N}}(r)}\bigg{|}_{\xi=1/6}=1-\frac{4}{3}\,e^{-m_{2}r}\,. (14)

The Yukawa amplitudes 4/3-4/3 (spin-2) and +1/3+1/3 (spin-0) are fixed by the spin structure of the graviton propagator and are independent of Λ\Lambda, ξ\xi, or any coupling constant. This is a key distinction from generic Yukawa-modified gravity scenarios, where the coupling strengths are free parameters.

II.3 Properties of the potential

  1. 1.

    Finite at the origin. The potential ratio (12) satisfies V(0)/VN(0)=14/3+1/3=0V(0)/V_{\mathrm{N}}(0)=1-4/3+1/3=0: the Newtonian 1/r1/r singularity is cancelled and the physical potential V(r)GM(4m2/3m0/3)V(r)\to-GM(4m_{2}/3-m_{0}/3) is finite as r0r\to 0.

  2. 2.

    Newtonian recovery. Both Yukawa terms decay exponentially, so V(r)VN(r)V(r)\to V_{\mathrm{N}}(r) for r1/Λr\gg 1/\Lambda.

  3. 3.

    Mass ratio. At minimal coupling ξ=0\xi=0: m2/m0=10/130.877m_{2}/m_{0}=\sqrt{10/13}\approx 0.877, fixed by the SM spectrum.

  4. 4.

    Scalar decoupling. At ξ=1/6\xi=1/6, the scalar mass diverges and the potential reduces to (14). The Newtonian singularity reappears: V(r)/VN(r)1/3V(r)/V_{\mathrm{N}}(r)\to-1/3 as r0r\to 0.

The potential ratio is plotted in Fig. 1 for several values of ξ\xi.

Refer to caption
Figure 1: Modified Newtonian potential V(r)/VN(r)V(r)/V_{\mathrm{N}}(r) as a function of rΛr\Lambda for several values of the non-minimal coupling ξ\xi. At conformal coupling ξ=1/6\xi=1/6 (dashed), the scalar mode decouples and the potential dips to 1/3-1/3 before recovering Newton’s law. For all ξ1/6\xi\neq 1/6, the potential is finite at the origin (V(0)=0V(0)=0) and approaches VNV_{\mathrm{N}} exponentially for r1/Λr\gg 1/\Lambda.

III Post-Newtonian Parameter γ\gamma

III.1 Definition and derivation

The Eddington–Robertson–Schiff parameter γ\gamma measures the ratio of the spatial and temporal metric potentials:

γ(r)Ψ(r)Φ(r).\gamma(r)\equiv\frac{\Psi(r)}{\Phi(r)}\,. (15)

In GR, γ=1\gamma=1 identically. Since ΦN=ΨN\Phi_{\mathrm{N}}=\Psi_{\mathrm{N}} in the Newtonian limit, the ratio of Eqs. (13) and (12) gives

γ(r)=123em2r13em0r143em2r+13em0r.\boxed{\gamma(r)=\frac{1-\frac{2}{3}\,e^{-m_{2}r}-\frac{1}{3}\,e^{-m_{0}r}}{1-% \frac{4}{3}\,e^{-m_{2}r}+\frac{1}{3}\,e^{-m_{0}r}}\,.} (16)

III.2 Limiting behavior

For r1/m2r\gg 1/m_{2} (which includes all solar-system distances for any viable Λ\Lambda), both exponentials are negligible and γ1\gamma\to 1 with corrections

γ(r)123em2r(r1/m2).\gamma(r)-1\approx\frac{2}{3}\,e^{-m_{2}r}\qquad(r\gg 1/m_{2})\,. (17)

The leading correction is positive, spin-2 dominated, and exponentially suppressed.

At r=0r=0 with ξ1/6\xi\neq 1/6:

γ(0)=2m2+m04m2m0.\gamma(0)=\frac{2m_{2}+m_{0}}{4m_{2}-m_{0}}\,. (18)

At minimal coupling (ξ=0\xi=0), γ(0)1.098\gamma(0)\approx 1.098; at conformal coupling (ξ=1/6\xi=1/6), γ(0)=1\gamma(0)=-1, reflecting the absent scalar cancellation.

III.3 Effective PPN table

For Λ\Lambda satisfying the experimental bounds derived below, the departure from GR is negligible at all solar-system distances. The effective PPN parameters are listed in Table 1. The preferred-frame parameters α1=α2=α3=0\alpha_{1}=\alpha_{2}=\alpha_{3}=0 and the conservation-law parameters ζi=0\zeta_{i}=0 follow from the diffeomorphism invariance of the spectral action [Will:2014, Chamseddine:1996zu].

Table 1: Effective PPN parameters of the spectral action, evaluated at solar-system distances r1/Λr\gg 1/\Lambda. The parameter β\beta requires the nonlinear 𝒪(h2)\mathcal{O}(h^{2}) field equations and is not derived here.
Parameter GR value SCT value Source
γ\gamma 11 1+𝒪(em2r)1+\mathcal{O}(e^{-m_{2}r}) This work
β\beta 11 Not yet derived
ξPPN\xi_{\mathrm{PPN}} 0 0 Diffeo. inv.
α1,2,3\alpha_{1,2,3} 0 0 Diffeo. inv.
ζ1,2,3,4\zeta_{1,2,3,4} 0 0 Conserv. law

IV Experimental Bounds on Λ\Lambda

The spectral scale Λ\Lambda sets the Yukawa range λi=1/mi\lambda_{i}=1/m_{i}. Experiments that probe deviations from the inverse-square law or from GR constrain Λ\Lambda from below. We consider five classes of experiments, ordered by the strength of the resulting bound.

IV.1 Torsion-balance experiments (Eöt-Wash)

The most precise tests of the gravitational inverse-square law at sub-millimeter distances are torsion-balance experiments by the Eöt-Wash group [Lee:2020, Kapner:2007, Adelberger:2009]. The most recent result [Lee:2020] excludes Yukawa deviations with |α|1|\alpha|\geq 1 at ranges λ38.6μm\lambda\geq 38.6\,\mu\mathrm{m} (95% CL).

The SCT spin-2 Yukawa has coupling |α1|=4/3>1|\alpha_{1}|=4/3>1, so the exclusion curve is crossed at a shorter range than the |α|=1|\alpha|=1 boundary. From the published exclusion contour at |α|=4/3|\alpha|=4/3:

λ1=1m2<35.8μm,\lambda_{1}=\frac{1}{m_{2}}<35.8\,\mu\mathrm{m}\,, (19)

giving

Λ>2.565×103eV(95% CL, Eöt-Wash).\boxed{\Lambda>2.565\times 10^{-3}\,\mathrm{eV}\quad\text{(95\% CL, E\"{o}t-% Wash)}\,.} (20)

This is the primary laboratory bound on the SCT spectral scale. The corresponding physical length is 1/Λ=76.8μm1/\Lambda=76.8\,\mu\mathrm{m}, and the effective masses are m25.5×103eVm_{2}\approx 5.5\times 10^{-3}\,\mathrm{eV} and m06.3×103eVm_{0}\approx 6.3\times 10^{-3}\,\mathrm{eV} (at ξ=0\xi=0). The scalar Yukawa (|α2|=1/3|\alpha_{2}|=1/3) produces a weaker constraint because |α2|<1|\alpha_{2}|<1 lies below the current exclusion curve at all explored ranges.

IV.2 Cassini Shapiro time delay

The Cassini spacecraft measurement of the Shapiro time delay [Bertotti:2003] constrains |γ1|<2.3×105|\gamma-1|<2.3\times 10^{-5} at r1AUr\simeq 1\,\mathrm{AU}. Using (17):

23em2rAU<2.3×105.\frac{2}{3}\,e^{-m_{2}r_{\mathrm{AU}}}<2.3\times 10^{-5}\,. (21)

With rAU1.496×1011mr_{\mathrm{AU}}\simeq 1.496\times 10^{11}\,\mathrm{m}:

Λ6.3×1018eV(Cassini).\Lambda\gtrsim 6.3\times 10^{-18}\,\mathrm{eV}\quad\text{(Cassini)}\,. (22)

This is fourteen orders of magnitude weaker than the Eöt-Wash bound, reflecting the exponential suppression at astronomical distances.

IV.3 MESSENGER radioscience

The MESSENGER radioscience analysis [Verma:2014] provides |γ1|<2.5×105|\gamma-1|<2.5\times 10^{-5}, yielding a comparable bound:

Λ6.3×1018eV(MESSENGER).\Lambda\gtrsim 6.3\times 10^{-18}\,\mathrm{eV}\quad\text{(MESSENGER)}\,. (23)

IV.4 Casimir experiments

Casimir force measurements at sub-micrometer distances [Decca:2007, Chen:2016] require Yukawa couplings |α|3000|\alpha|\gtrsim 3000 to be detectable, far exceeding the SCT values |α1|=4/3|\alpha_{1}|=4/3 and |α2|=1/3|\alpha_{2}|=1/3. At r=0.1μmr=0.1\,\mu\mathrm{m} with Λ=2.565×103eV\Lambda=2.565\times 10^{-3}\,\mathrm{eV}, the potential deviation is only 0.3%\sim 0.3\%, well below the experimental sensitivity. Casimir experiments provide no constraint on Λ\Lambda.

IV.5 Satellite experiments and Lunar Laser Ranging

Gravity Probe B [Everitt:2011], MICROSCOPE [Touboul:2022], Lunar Laser Ranging [Williams:2004], and atom interferometry experiments probe gravity at centimeter-to-orbital distances, where the SCT Yukawa corrections are suppressed by factors of em2re102e^{-m_{2}r}\sim e^{-10^{2}} (atom interferometry at 1 cm) to e1013e^{-10^{13}} (Earth–Moon). These experiments provide zero effective constraint on Λ\Lambda.

LLR also probes the Nordtvedt parameter ηN=4βγ3\eta_{N}=4\beta-\gamma-3 [Williams:2004], but since β\beta has not been derived in the spectral action framework (it requires nonlinear field equations), this bound cannot yet be evaluated. If β=1\beta=1 at leading order (as in GR), the Nordtvedt parameter vanishes automatically. (MICROSCOPE [Touboul:2022] tests the weak equivalence principle on test masses, not the Nordtvedt strong-EP parameter.)

IV.6 Summary of bounds

Table 2: Lower bounds on the spectral scale Λ\Lambda from gravitational experiments, derived within the local Yukawa approximation. The Eöt-Wash bound is the only experimentally meaningful constraint.
Experiment Λmin\Lambda_{\mathrm{min}} (eV) Constraint
Eöt-Wash [Lee:2020] 2.565×1032.565\times 10^{-3} |α|<4/3|\alpha|<4/3 at 36μ36\,\mum
Cassini [Bertotti:2003] 6.3×10186.3\times 10^{-18} |γ1|<2.3×105|\gamma-1|<2.3\times 10^{-5}
MESSENGER [Verma:2014] 6.3×10186.3\times 10^{-18} |γ1|<2.5×105|\gamma-1|<2.5\times 10^{-5}
Casimir [Decca:2007] No constraint |α||αlim||\alpha|\ll|\alpha_{\mathrm{lim}}|
GP-B [Everitt:2011] No constraint em2r0e^{-m_{2}r}\approx 0
MICROSCOPE [Touboul:2022] Requires β\beta Not yet derived
LLR [Williams:2004] Requires β\beta Not yet derived
Refer to caption
Figure 2: Unified exclusion plot in the (λ,|α|)(\lambda,|\alpha|) plane, combining torsion-balance, Casimir, and satellite constraints. The horizontal lines mark the SCT Yukawa couplings: |α1|=4/3|\alpha_{1}|=4/3 (spin-2, solid) and |α2|=1/3|\alpha_{2}|=1/3 (spin-0, dashed). The Eöt-Wash exclusion curve intersects |α1|=4/3|\alpha_{1}|=4/3 at λ36μm\lambda\approx 36\,\mu\mathrm{m}, setting Λ>2.565×103eV\Lambda>2.565\times 10^{-3}\,\mathrm{eV}.

V Dependence on the non-minimal coupling

The Eöt-Wash bound depends on ξ\xi through m0(ξ)m_{0}(\xi): at conformal coupling ξ=1/6\xi=1/6, the scalar Yukawa decouples and only the spin-2 channel contributes (the bound is then set entirely by m2m_{2}). Away from conformal coupling, both Yukawa terms contribute but the spin-2 channel dominates because |α1|=4/3>|α2|=1/3|\alpha_{1}|=4/3>|\alpha_{2}|=1/3.

The minimum spectral scale as a function of ξ\xi is plotted in Fig. 3. The bound varies by less than 8% across the full range 0ξ10\leq\xi\leq 1, reflecting the dominance of the spin-2 constraint. At ξ=0\xi=0: Λmin=2.565×103eV\Lambda_{\mathrm{min}}=2.565\times 10^{-3}\,\mathrm{eV}; at ξ=1/6\xi=1/6: Λmin=2.38×103eV\Lambda_{\mathrm{min}}=2.38\times 10^{-3}\,\mathrm{eV}.

Refer to caption
Figure 3: Minimum spectral scale Λmin\Lambda_{\mathrm{min}} as a function of the Higgs non-minimal coupling ξ\xi, from the Eöt-Wash constraint. The bound is nearly ξ\xi-independent because the spin-2 Yukawa channel dominates. At ξ=1/6\xi=1/6 (conformal coupling, vertical dashed), the scalar mode decouples and only the spin-2 constraint remains.
Refer to caption
Figure 4: Potential deviation |V(r)/VN(r)1||V(r)/V_{\mathrm{N}}(r)-1| as a function of distance, evaluated at the boundary value Λ=2.565×103eV\Lambda=2.565\times 10^{-3}\,\mathrm{eV}. The deviation is 𝒪(1)\mathcal{O}(1) at r36μmr\lesssim 36\,\mu\mathrm{m} and drops below 10310^{-3} by r0.1mmr\sim 0.1\,\mathrm{mm}. The experimental sensitivity thresholds for Eöt-Wash and Casimir experiments are indicated.

VI Discussion

VI.1 Validity of the Yukawa approximation

The local Yukawa approximation replaces the full propagator denominator ΠTT(z)=1+(13/60)zF^1(z)\Pi_{\mathrm{TT}}(z)=1+(13/60)\,z\,\hat{F}_{1}(z) by its small-zz limit 1+(13/60)z1+(13/60)\,z. The exact dressed propagator has a zero at z02.41z_{0}\approx 2.41 [Alfyorov:2026ff], yielding an exact spin-2 mass m2exact=Λz01.55Λm_{2}^{\mathrm{exact}}=\Lambda\sqrt{z_{0}}\approx 1.55\,\Lambda, smaller than the Yukawa-approximation value m2loc2.15Λm_{2}^{\mathrm{loc}}\approx 2.15\,\Lambda by a factor of 1.381.38. The exact Yukawa range is therefore longer than the approximation predicts, making the stated Eöt-Wash bound on Λ\Lambda conservative: the full nonlocal analysis would yield a tighter constraint.

For solar-system experiments (r1/Λr\gg 1/\Lambda), the Yukawa corrections are exponentially suppressed regardless of the approximation, so the local and exact results are indistinguishable at current precision.

VI.2 Comparison with Stelle gravity

In Stelle’s local quadratic gravity [Stelle:1976gc, Stelle:1978], the potential has the same functional form (12) but with arbitrary masses m2m_{2} and m0m_{0} determined by two free coupling constants. The spectral action eliminates this freedom: the couplings c2=13/60c_{2}=13/60 and 3c1+c2=6(ξ1/6)23c_{1}+c_{2}=6(\xi-1/6)^{2} are computed from the SM particle content, yielding definite effective masses (9) and (10) as functions of Λ\Lambda alone (plus ξ\xi).

VI.3 Ghost interpretation

The spin-2 mass m2m_{2} corresponds to a ghost mode: the coefficient of em2re^{-m_{2}r} in (12) is 4/3<0-4/3<0, indicating that the massive spin-2 graviton couples with wrong-sign residue. This is the well-known Ostrogradsky ghost of higher-derivative gravity [Stelle:1976gc, Woodard:2015zca]. In the spectral action framework, the dressed propagator 1/ΠTT(z)1/\Pi_{\mathrm{TT}}(z) has a pole at z02.41z_{0}\approx 2.41 with residue R20.49R_{2}\approx-0.49, approximately 50% suppressed relative to the Stelle value R2Stelle=1R_{2}^{\mathrm{Stelle}}=-1. The physical interpretation of this ghost (Lee–Wick field, fakeon, or dark matter candidate) is deferred to a dedicated analysis.

VI.4 Parameter-free nature

The modified potential (12) is determined by two inputs: the spectral scale Λ\Lambda and the Higgs non-minimal coupling ξ\xi. The Yukawa amplitudes (4/3-4/3, +1/3+1/3) are fixed by spin decomposition. The mass ratios m2/Λm_{2}/\Lambda and m0/Λm_{0}/\Lambda are fixed by the SM content. This structure means that a single measurement of the Yukawa range determines Λ\Lambda uniquely, and all other predictions follow.

VI.5 Implications for the spectral scale

The Eöt-Wash bound Λ>2.565×103eV\Lambda>2.565\times 10^{-3}\,\mathrm{eV} corresponds to a length scale 1/Λ<77μm1/\Lambda<77\,\mu\mathrm{m}, far above the Planck length P1035m\ell_{P}\sim 10^{-35}\,\mathrm{m}. The spectral scale is therefore not directly constrained to lie near the Planck mass by current experiments; it could be much higher. Future improvements in short-range gravity experiments, targeting λ10μm\lambda\lesssim 10\,\mu\mathrm{m}, would push the bound to Λ102eV\Lambda\gtrsim 10^{-2}\,\mathrm{eV}.

VII Conclusions

We have derived the post-Newtonian parameter and laboratory constraints on the spectral action with Standard Model content. The modified Newtonian potential (12) has its Yukawa amplitudes, with effective masses set by the particle content and the spectral scale Λ\Lambda.

The central experimental result is:

Λ>2.565×103eV(95% CL, Eöt-Wash),\Lambda>2.565\times 10^{-3}\,\mathrm{eV}\quad\text{(95\% CL, E\"{o}t-Wash)}\,,

corresponding to 1/Λ<77μm1/\Lambda<77\,\mu\mathrm{m} (Yukawa range λ1<36μm\lambda_{1}<36\,\mu\mathrm{m}). All solar-system tests are satisfied with exponential margin. Casimir, atom interferometry, and satellite experiments are below their sensitivity thresholds for the SCT coupling strengths.

The spectral action passes all currently available gravitational tests. The most promising avenue for future improvement is the extension of torsion-balance experiments to shorter ranges, which would directly constrain Λ\Lambda in the meV regime.

Acknowledgements.
Numerical results were obtained using mpmath and scipy and verified to 100-digit precision.

AI Disclosure

In preparing this manuscript, large language model tools (specifically Anthropic Claude Opus 4.7 and OpenAI ChatGPT GPT-5.5 Pro) were used for English-language editorial assistance, bibliographic cross-checking, LaTeX markup verification, and the writing of computational code and numerical-verification scripts that were subsequently audited by the author. All mathematical derivations, the choice of observables, the comparison against Cassini, lunar-laser-ranging, and short-range torsion-balance constraints, the figures, the interpretations, and the conclusions are the author’s own. Every numerical value reported in the manuscript was computed by software written by the author (with AI assistance as declared) and was cross-checked against the canonical-validator modules of the upstream Spectral Causal Theory program at 100-digit precision. Every citation was individually verified against arXiv, Crossref, or the journal page of record. AI tools are not listed as authors and bear no authorial responsibility for this manuscript.

Data Availability

The computational tools and experimental data references used in this work are available from the author upon reasonable request.

References