THE GOLDEN BALLROOM/BUNKER

\documentclass[11pt]{article} \usepackage{amsmath, amssymb} \usepackage{graphicx} \usepackage{hyperref} \usepackage{geometry} \geometry{margin=1in} \title{FRCMFD: Metric Emergence, PPN Calculation, and Empirical Constraints} \author{FRCMFD Collaboration} \date{May 19, 2026} \begin{document} \maketitle \begin{abstract} We derive the weak-field limit of a version-locked FRCMFD field equation, construct an emergent effective metric from response gradients, and compare the resulting phenomenology with empirical constraints. The weak-field sector yields a screened Poisson equation and a conformal effective metric, from which the Newtonian potential is obtained as a linear functional of the field perturbation. At linear order, the conformal metric implies a PPN parameter $\gamma = -1$, which is incompatible with Cassini solar system bounds. Including a disformal term in the metric ansatz does not significantly modify $\gamma$ at solar system scales under weak-field assumptions. Independent fits to SPARC galaxy data show no statistically significant evidence for Yukawa-type screening in the baryonic residuals. The framework remains mathematically consistent but observationally constrained; any viable branch requires nonlinear screening beyond the linear analysis presented here. \end{abstract} \tableofcontents \section{Core Field Equation (Version-Locked)} We start from the version-locked FRCMFD substrate field equation with additive matter coupling: \begin{equation} \partial_\mu \partial^\mu \Psi + \mu \Psi + \lambda \Psi^3 = \kappa S, \label{eq:field} \end{equation} where: \begin{itemize} \item $\Psi$ is the substrate excitation field, \item $\mu$ is a restoring coefficient (assumed $\mu < 0$), \item $\lambda$ is a nonlinear saturation coefficient, \item $\kappa$ is a source coupling coefficient, \item $S$ is a matter source density. \end{itemize} \subsection{Background solution (vacuum)} For $S = 0$ and a static, homogeneous configuration, Eq.~\eqref{eq:field} reduces to \begin{equation} \mu \Psi_0 + \lambda \Psi_0^3 = 0. \end{equation} Assuming $\Psi_0 \neq 0$, this gives \begin{equation} \Psi_0^2 = -\frac{\mu}{\lambda}, \label{eq:background} \end{equation} which requires $\mu < 0$ and $\lambda > 0$ for a real $\Psi_0$. \subsection{Linearized perturbation} Write \begin{equation} \Psi = \Psi_0 + \delta\Psi, \qquad |\delta\Psi| \ll |\Psi_0|. \end{equation} Substituting into Eq.~\eqref{eq:field} and expanding to first order in $\delta\Psi$ yields \begin{equation} \partial_\mu \partial^\mu \delta\Psi + (\mu + 3\lambda \Psi_0^2)\,\delta\Psi = \kappa S. \label{eq:linearized} \end{equation} Using Eq.~\eqref{eq:background}, \begin{equation} \mu + 3\lambda \Psi_0^2 = \mu + 3\lambda\left(-\frac{\mu}{\lambda}\right) = -2\mu. \end{equation} Define \begin{equation} m_{\text{eff}}^2 \equiv -2\mu > 0, \end{equation} so that \begin{equation} \partial_\mu \partial^\mu \delta\Psi + m_{\text{eff}}^2 \delta\Psi = \kappa S. \label{eq:klein_gordon} \end{equation} This is a sourced Klein--Gordon equation for the perturbation $\delta\Psi$. \subsection{Static weak-field limit} In the static limit ($\partial_t \delta\Psi = 0$), Eq.~\eqref{eq:klein_gordon} becomes \begin{equation} -\nabla^2 \delta\Psi + m_{\text{eff}}^2 \delta\Psi = \kappa S, \label{eq:screened_poisson} \end{equation} or equivalently \begin{equation} \nabla^2 \delta\Psi - m_{\text{eff}}^2 \delta\Psi = -\kappa S. \label{eq:screened_poisson2} \end{equation} This is a screened Poisson (Yukawa-type) equation for $\delta\Psi$. \section{Emergent metric from response gradients} \subsection{Response gradient tensor} Define the rank-2 tensor constructed from field gradients: \begin{equation} R_{\mu\nu} = \partial_\mu \Psi\,\partial_\nu \Psi. \label{eq:R_tensor} \end{equation} This object is analogous to the kinetic part of a scalar-field stress tensor and serves as a natural building block for an emergent effective metric. \subsection{Disformal metric ansatz} We postulate an effective metric of disformal form: \begin{equation} g_{\mu\nu}^{\text{eff}} = A(\Psi)\,\eta_{\mu\nu} + B(\Psi)\,\partial_\mu \Psi\,\partial_\nu \Psi, \label{eq:disformal} \end{equation} where $\eta_{\mu\nu} = \mathrm{diag}(-1, +1, +1, +1)$ is the Minkowski metric, and $A(\Psi)$, $B(\Psi)$ are scalar functions to be specified. \subsection{Weak-field expansion} Write \begin{equation} \Psi = \Psi_0 + \delta\Psi, \qquad |\delta\Psi| \ll |\Psi_0|. \end{equation} Expanding Eq.~\eqref{eq:disformal} to first order in $\delta\Psi$ gives \begin{equation} g_{\mu\nu}^{\text{eff}} \approx (A_0 + A_1 \delta\Psi)\,\eta_{\mu\nu} + B_0\,\partial_\mu \delta\Psi\,\partial_\nu \delta\Psi, \label{eq:weak_expansion} \end{equation} where \begin{equation} A_0 = A(\Psi_0), \qquad A_1 = \left.\frac{dA}{d\Psi}\right|_{\Psi_0}, \qquad B_0 = B(\Psi_0). \end{equation} The term proportional to $B_0$ is quadratic in $\delta\Psi$ and its gradients, so at strictly linear order in the perturbation amplitude it can be neglected. We then have \begin{equation} g_{\mu\nu}^{\text{eff}} \approx (A_0 + A_1 \delta\Psi)\,\eta_{\mu\nu}. \end{equation} Choosing units such that $A_0 = 1$ (Minkowski background), this simplifies to \begin{equation} g_{\mu\nu}^{\text{eff}} \approx (1 + A_1 \delta\Psi)\,\eta_{\mu\nu}. \label{eq:conformal} \end{equation} Thus, at linear order the effective metric is conformally related to Minkowski. \subsection{Newtonian potential identification} In the weak-field limit, the time-time component of the metric is written as \begin{equation} g_{00}^{\text{eff}} = -(1 + 2\Phi), \label{eq:g00_def} \end{equation} where $\Phi$ is the Newtonian gravitational potential. From Eq.~\eqref{eq:conformal}, \begin{equation} g_{00}^{\text{eff}} \approx -(1 + A_1 \delta\Psi). \end{equation} Equating the two expressions, \begin{equation} 1 + 2\Phi = 1 + A_1 \delta\Psi \quad \Rightarrow \quad \Phi = \frac{A_1}{2}\,\delta\Psi. \label{eq:Phi_propto_deltaPsi} \end{equation} Thus, within this conformal ansatz, $\Phi$ is linearly proportional to $\delta\Psi$; this relation is now a consequence of the metric structure rather than an independent ansatz. \subsection{Matching to the field equation} We require that the propagation of $\delta\Psi$ as encoded in Eq.~\eqref{eq:klein_gordon} be consistent with the effective metric structure. In particular, the d'Alembertian constructed from $g_{\mu\nu}^{\text{eff}}$ should reduce to the flat-space d'Alembertian at leading order. This fixes the normalization of $A_1$ up to a constant factor. We parameterize this as \begin{equation} A_1 = \frac{2\kappa \Psi_0}{v^2}, \label{eq:A1_determined} \end{equation} where $v$ is the propagation speed (often set to $c$ in natural units). Substituting into Eq.~\eqref{eq:Phi_propto_deltaPsi} gives \begin{equation} \Phi = \frac{\kappa \Psi_0}{v^2}\,\delta\Psi. \label{eq:Phi_final} \end{equation} \subsection{Effective Newtonian constant} Comparing the screened Poisson equation for $\delta\Psi$ with the Poisson equation for $\Phi$, \begin{equation} \nabla^2 \Phi = 4\pi G_{\text{eff}} \rho + \text{(screening terms)}, \end{equation} and using Eq.~\eqref{eq:screened_poisson2} together with Eq.~\eqref{eq:Phi_final}, one finds that the effective Newtonian constant can be written as \begin{equation} G_{\text{eff}} = \frac{\kappa \Psi_0}{4\pi v^2}, \label{eq:G_eff} \end{equation} up to the usual caveats about screening and the precise identification of $S$ with a mass density $\rho$. \section{Disformal term and action-level considerations} \subsection{Provisional scalar-field action} A minimal scalar-field Lagrangian density that reproduces Eq.~\eqref{eq:field} in flat spacetime is \begin{equation} \mathcal{L} = \frac{1}{2}\partial_\mu\Psi\,\partial^\mu\Psi - V(\Psi) + \kappa S \Psi, \label{eq:action} \end{equation} with potential \begin{equation} V(\Psi) = \frac{1}{2}\mu \Psi^2 + \frac{1}{4}\lambda \Psi^4. \label{eq:potential} \end{equation} Varying this action with respect to $\Psi$ yields Eq.~\eqref{eq:field} in flat spacetime. \subsection{Disformal metric ansatz revisited} Following disformal constructions in the literature, one may consider \begin{equation} g_{\mu\nu}^{\text{eff}} = A(\Psi)\,\eta_{\mu\nu} + B(\Psi)\,\partial_\mu\Psi\,\partial_\nu\Psi, \end{equation} with $A(\Psi)$ and $B(\Psi)$ chosen such that matter fields couple minimally to $g_{\mu\nu}^{\text{eff}}$ while $\Psi$ retains a canonical kinetic term in an appropriate frame. In this document we treat specific forms of $A(\Psi)$ and $B(\Psi)$ as model choices rather than fully derived consequences. For example, one may consider \begin{equation} A(\Psi) = 1 - \frac{2\kappa \Psi}{v^2}, \qquad B(\Psi) = \frac{2\kappa}{v^2 \Psi_0^2}, \label{eq:A_B_example} \end{equation} as a concrete parametrization consistent with the linear relation in Eq.~\eqref{eq:A1_determined}. These forms should be regarded as illustrative rather than uniquely determined. \section{PPN parameter $\gamma$} \subsection{Standard PPN form} In the parametrized post-Newtonian (PPN) formalism, the weak-field metric for a static, spherically symmetric source is written as \begin{equation} g_{00} = -(1 + 2\Phi), \qquad g_{0i} = 0, \qquad g_{ij} = (1 - 2\gamma \Phi)\,\delta_{ij}, \label{eq:ppn_metric} \end{equation} where $\gamma$ is a PPN parameter. General Relativity predicts $\gamma = 1$. \subsection{Linear-order result from the conformal metric} From the conformal metric in Eq.~\eqref{eq:conformal}, we have \begin{equation} g_{ij}^{\text{eff}} \approx (1 + A_1 \delta\Psi)\,\delta_{ij} = (1 + 2\Phi)\,\delta_{ij}, \label{eq:conformal_gij} \end{equation} using Eq.~\eqref{eq:Phi_propto_deltaPsi}. Comparing with Eq.~\eqref{eq:ppn_metric}, \begin{equation} 1 + 2\Phi = 1 - 2\gamma \Phi \quad \Rightarrow \quad 2\Phi = -2\gamma \Phi. \end{equation} For nonzero $\Phi$, this implies \begin{equation} \gamma = -1. \label{eq:gamma_minus1} \end{equation} Thus, at linear order in the conformal metric, the PPN parameter $\gamma$ is fixed to $-1$. \subsection{Solar system constraint} Cassini tracking experiments constrain $\gamma$ via Shapiro time delay measurements. The bound can be expressed as \begin{equation} |\gamma - 1| < 2.3 \times 10^{-5}. \label{eq:cassini} \end{equation} The linear conformal result $\gamma = -1$ is therefore incompatible with this constraint by many orders of magnitude. \section{Disformal corrections and their impact on $\gamma$} \subsection{Including the disformal term} For a static, spherically symmetric configuration, the spatial part of the metric including the disformal term can be written as \begin{equation} g_{ij}^{\text{eff}} = A(\Psi)\,\delta_{ij} + B(\Psi)\,\partial_i\Psi\,\partial_j\Psi. \end{equation} In radial coordinates, this yields \begin{equation} g_{rr}^{\text{eff}} = A(\Psi) + B(\Psi)\,(\partial_r\Psi)^2, \qquad g_{\theta\theta}^{\text{eff}} = A(\Psi)\,r^2. \end{equation} \subsection{Effective $\gamma$ with disformal term} In a PPN-like identification, one can define an effective $\gamma$ from the radial component as \begin{equation} \gamma_{\text{eff}} \equiv \frac{g_{rr}^{\text{eff}} - 1}{2\Phi}, \end{equation} in the regime where $|\Phi| \ll 1$. Substituting \begin{equation} g_{rr}^{\text{eff}} = A(\Psi) + B(\Psi)\,(\partial_r\Psi)^2, \end{equation} and using $A(\Psi) \approx 1 + 2\Phi$ at linear order, we obtain schematically \begin{equation} \gamma_{\text{eff}} \approx 1 + \frac{B_0\,(\partial_r\Psi)^2}{2\Phi}. \end{equation} \subsection{Magnitude of disformal corrections at solar system scales} At Earth orbit around the Sun, one has $|\Phi| \sim 10^{-8}$ and $|\partial_r\Phi| \sim 10^{-16}\,\mathrm{m}^{-1}$ as order-of-magnitude estimates. Consequently, $(\partial_r\Phi)^2/|\Phi|$ is extremely small. For reasonable choices of $B_0$ and the coupling parameters, the disformal contribution to $\gamma_{\text{eff}}$ is negligible compared to the leading conformal term. In such a regime, the effective $\gamma$ remains close to the conformal value $\gamma = -1$, and the Cassini bound is not satisfied. \subsection{Nonlinear screening as a possible escape} The above analysis assumes that the perturbation $\delta\Psi$ remains small compared to $\Psi_0$ and that the linearized treatment is valid. A possible way to evade the $\gamma = -1$ result is if nonlinear terms in the field equation become important in the solar system, effectively screening the conformal contribution. Schematically, one would require terms such as \begin{equation} \lambda \Psi_0 (\delta\Psi)^2 \sim \mu \delta\Psi, \end{equation} which implies \begin{equation} \delta\Psi \sim \Psi_0. \end{equation} At Earth orbit, $|\Phi| \sim 10^{-8}$, and via Eq.~\eqref{eq:Phi_final} this suggests that $\Psi_0$ would need to be of comparable order for nonlinear screening to be relevant. This parameter regime is not excluded by the present derivation but is not yet demonstrated to yield $\gamma \approx 1$; a full nonlinear analysis would be required. \section{Empirical constraints from SPARC data} \subsection{Model comparison: flat vs Yukawa} Using SPARC baryonic residuals $\Delta\gamma_{\text{resid}}$ and a sample of $N = 80$ galaxies, one can test for distance-dependent trends consistent with Yukawa screening. A flat (distance-independent) model and a Yukawa-type model are compared via their residual sum of squares (RSS): \begin{center} \begin{tabular}{|l|l|} \hline Model & RSS \\ \hline Flat (constant) & 80.0000 \\ Yukawa & 79.5243 \\ \hline \end{tabular} \end{center} The Yukawa model has one additional parameter (the screening scale), so information criteria are used to penalize model complexity. \subsection{Information criteria} The differences in Akaike and Bayesian information criteria are \begin{align} \Delta\mathrm{AIC} &= \mathrm{AIC}_{\text{Flat}} - \mathrm{AIC}_{\text{Yukawa}} = -1.52, \\ \Delta\mathrm{BIC} &= \mathrm{BIC}_{\text{Flat}} - \mathrm{BIC}_{\text{Yukawa}} = -3.90. \end{align} For $N = 80$, the BIC penalty for one extra parameter is $\ln N \approx 4.38$, which dominates the small reduction in RSS. Both criteria favor the flat model; the BIC difference indicates at least moderate evidence against the Yukawa model in this dataset. \subsection{Bootstrap uncertainty on the screening parameter} A bootstrap analysis of the Yukawa screening parameter $\mu$ yields \begin{equation} \mu = 12.47 \pm 13.53\ \mathrm{Mpc}^{-1}. \end{equation} Since the uncertainty exceeds the central value, the confidence interval includes $\mu = 0$. This indicates that the screening parameter is not statistically identifiable from the data at the current sensitivity. \subsection{Mathematically precise statement} Within the SPARC $\Delta\gamma_{\text{resid}}$ sample analyzed here, there is no statistically significant evidence for a Yukawa-type distance dependence. A flat, distance-independent model is preferred by both AIC and BIC, and the Yukawa screening parameter $\mu$ is not constrained away from zero. This constitutes a null result for Yukawa screening in this particular weak-field test. \section{Current status summary} The mathematical and empirical status of the framework can be summarized as follows: \begin{itemize} \item The version-locked field equation and its background solution are well-defined. \item Linearization around the background yields a sourced Klein--Gordon equation with an effective mass $m_{\text{eff}}^2 = -2\mu$. \item The static weak-field limit leads to a screened Poisson equation for $\delta\Psi$. \item A disformal metric ansatz with a conformal leading term yields an effective metric at linear order. \item The Newtonian potential $\Phi$ is derived as a linear functional of $\delta\Psi$ within this ansatz. \item The effective Newtonian constant $G_{\text{eff}}$ can be expressed in terms of $\kappa$, $\Psi_0$, and $v$. \item At linear order, the conformal metric implies a PPN parameter $\gamma = -1$, incompatible with Cassini bounds. \item Disformal corrections are suppressed in weak, slowly varying fields and do not significantly alter $\gamma$ at solar system scales under the assumptions used. \item SPARC data show no evidence for Yukawa screening in the baryonic residuals; a flat model is statistically preferred. \item Nonlinear screening remains a possible but unproven mechanism to reconcile the framework with solar system constraints. \end{itemize} \section{Conclusion} We have presented a mathematically consistent weak-field derivation for a version-locked FRCMFD field equation, constructed an emergent effective metric from response gradients, and evaluated the resulting phenomenology against solar system and galactic-scale constraints. The weak-field sector yields a screened Poisson equation and a conformal effective metric in which the Newtonian potential is proportional to the field perturbation. However, the linear-order conformal metric predicts $\gamma = -1$, which is strongly excluded by Cassini. Disformal corrections are negligible at solar system scales under the weak-field assumptions adopted here. Independent analysis of SPARC baryonic residuals finds no statistically significant evidence for Yukawa screening. Any viable branch of the framework must therefore rely on nonlinear screening or additional structure beyond the linear analysis presented. These requirements are mathematically clear but not yet demonstrated. \section*{One-line summary} The weak-field FRCMFD framework is mathematically consistent but, in its linear conformal form, predicts $\gamma = -1$ and no detectable Yukawa signal in SPARC data; reconciling it with solar system tests requires nonlinear screening beyond the present derivation. \begin{thebibliography}{99} \bibitem{Cassini} Bertotti, B., Iess, L., \& Tortora, P. 2003, Nature, 425, 374. \bibitem{Will2014} Will, C. M. 2014, Living Rev. Relativity, 17, 4. \bibitem{Bekenstein2004} Bekenstein, J. D. 2004, Phys. Rev. D, 70, 083509. \end{thebibliography} \end{document}