document.tex

% ============================================================
%  Radial-Cyclic Field Model — Full Paper Draft
%  Version 1.3 — Scientific corrections and first-principles derivations
%  Compile with: pdflatex document.tex && bibtex document &&
%                pdflatex document.tex && pdflatex document.tex
% ============================================================
\documentclass[10pt,twocolumn]{article}

% ── Packages ────────────────────────────────────────────────
\usepackage{amsmath,amssymb,amsfonts,amsthm}
\usepackage{cuted}
\usepackage[a4paper,margin=1.0in]{geometry}
\usepackage[colorlinks=true,
            linkcolor=blue!60!black,
            citecolor=blue!60!black,
            urlcolor=blue!60!black]{hyperref}
\usepackage{booktabs}
\usepackage{array}
\usepackage{xcolor}
\usepackage{graphicx}
\usepackage[numbers,sort&compress]{natbib}
\usepackage{microtype}
\usepackage{tocloft}
\usepackage{doi}
\usepackage{enumitem}
\usepackage{mdframed}
\usepackage{mathtools}
\usepackage{placeins}
\usepackage{bm}
\usepackage{caption}
\usepackage{subcaption}
\usepackage[T1]{fontenc}
\usepackage[utf8]{inputenc}
\usepackage{lmodern}

\usepackage{setspace}
\usepackage{varwidth} 
\usepackage{abstract}

% Hyphenation tweaks for better line breaks
\hyphenpenalty=5000
\exhyphenpenalty=5000
\emergencystretch=3em
\hyphenation{grav-i-ta-tion-al cos-mo-log-i-cal prim-or-di-al ef-fec-tive trans-for-ma-tion pa-ram-e-ter field-bond-ed non-inter-act-ing stand-ing-wave cos-mol-o-gy quan-ta}

% TOC dot spacing (higher value for fewer dots)
\renewcommand{\cftdotsep}{2}

% ── Theorem environments ─────────────────────────────────────
\theoremstyle{definition}
\newtheorem{openprob}{Open Problem}
\newtheorem{result}{Result}
\newtheorem{consistency}{Consistency Condition}

% ── Custom macros ────────────────────────────────────────────
\newcommand{\Rmax}{R_{\max}}
\newcommand{\rhoA}{\rho_{A}}
\newcommand{\rhoB}{\rho_{B}}
\newcommand{\Hcal}{\mathcal{H}}
\newcommand{\Geff}{G_{\mathrm{eff}}}
\newcommand{\Gammadrag}{\Gamma_{\mathrm{drag}}}
\newcommand{\rhoc}{\rho_{c}}
\newcommand{\rhoEW}{\rho_{\mathrm{EW}}}
\newcommand{\dOmega}{d\Omega_{3}^{2}}
\newcommand{\Mpl}{M_{\mathrm{Pl}}}

\newcommand{\LCDM}{\ensuremath{\Lambda\mathrm{CDM}}}

\newcommand{\todo}[1]{\textcolor{red}{\textbf{[TODO: #1]}}}



\makeatletter

% ── Improved LaTeX macros for two-column documents ──
% These replace the problematic \eqbox macro to avoid
% compilation issues, numbering problems, and TOC conflicts

% Single-column framed equation (for two-column documents)
\DeclareRobustCommand{\coleqbox}[1]{%
  \begin{mdframed}[
    linewidth=0.8pt,
    innertopmargin=6pt,
    innerbottommargin=6pt,
    skipabove=\smallskipamount,
    skipbelow=\smallskipamount
  ]%
    \begin{equation}
      #1
    \end{equation}%
  \end{mdframed}%
}

% Full-width unframed equation spanning both columns
\DeclareRobustCommand{\colspan}[1]{%
  \begin{strip}%
    \centering
    #1
  \end{strip}%
}

% Full-width unframed convenience wrapper for equations
\DeclareRobustCommand{\eqvirt}[1]{%
  \begin{strip}%
    \centering
    \begin{equation}
      #1
    \end{equation}%
  \end{strip}%
}

% Framed full-width equation (full-page width framed)
\DeclareRobustCommand{\colbox}[1]{%
  \begin{strip}%
    \begin{mdframed}[
      linewidth=0.8pt,
      innertopmargin=8pt,
      innerbottommargin=8pt,
      skipabove=\smallskipamount,
      skipbelow=\smallskipamount
    ]%
      \centering
      #1
    \end{mdframed}%
  \end{strip}%
}

% Helper for full-width tables in two-column article
\DeclareRobustCommand{\onecoltable}[2][]{%
  \begin{table*}[#1]%
    \centering
    #2
  \end{table*}%
}

% ── Legacy support (deprecated - use above macros instead) ──
% configurable frame parameters (edit once here)
\newlength{\eqboxinnerhmargin}
\newlength{\eqboxinnervmargin}
\setlength{\eqboxinnerhmargin}{8pt}
\setlength{\eqboxinnervmargin}{6pt}

% Legacy eqbox - now uses coleqbox for single column
% NOTE: This is kept for backward compatibility but should be
% replaced with \coleqbox, \colspan, \colbox, or \eqvirt
\newcommand{\eqbox}[1]{%
  \coleqbox{#1}%
}

% One-column table (no frame)
\newcommand{\coltable}[1]{%
  \begin{strip}%
    \vspace{6pt}%
    \noindent
    \hbox to \textwidth{\hfil #1\hfil}%
    \vspace{6pt}%
  \end{strip}%
}

\makeatother




% ── Metadata ────────────────────────────────────────────────


\title{%
  \textbf{The Radial-Cyclic Field Model (RCFM)}\\[6pt]
  \large A 4-Dimensional Topological Framework\\
  for Cyclic Cosmology
}
\author{%
  Jeremy Erich Resch\\[4pt]
  \small \textit{For Correspondence:}\\
  \small \texttt{jer.code.contact@gmail.com}
}
\date{\today \\ \small Version 1.3 --- Scientific corrections and first-principles derivations}

% ════════════════════════════════════════════════════════════

\begin{document}
\maketitle


\begin{abstract}
The Radial-Cyclic Field Model (RCFM) proposes a cosmology in
which the universe is a 4-dimensional hyperspherical standing
wave. Matter is generated at a continuously active central
singularity, streams outward as field-bonded matter on
expanding 3-spherical shells (Phase\nobreakspace{}A), undergoes a phase
transition at a critical radius \( \Rmax \) where quantum fields
can no longer be sustained, and returns inward as
non-interacting ghost quanta (Phase\nobreakspace{}B).

Starting from a metric ansatz on the 4-ball \( B^{4} \), we
derive a closed system of modified Friedmann equations
incorporating a dual-stream energy-momentum tensor with a
gravitational drag interaction whose microphysical origin is
dynamical friction via retarded gravitational wakes.  We show
that:
\begin{enumerate}[label=(\roman*)]
  \item Newtonian gravity is recovered exactly in the
    weak-field, sub-horizon limit.
  \item The Phase\nobreakspace{}B ghost quanta naturally produce an
    effective equation of state \( w_{B}\approx-1 \), mimicking
    a small positive cosmological constant without
    fine-tuning.
  \item The singularity boundary condition \( a(r)\propto r^{2} \)
    is self-consistent and generates a nearly scale-invariant
    primordial power spectrum without requiring inflation.
  \item The drag kernel \( \Gammadrag \) is spatially constant,
    yielding a fully closed ODE system, and is shown to be
    perturbatively stable to first order in deviations from
    the exact \( a\propto r^{2} \) solution.
  \item Gravitational wave propagation is consistent with
    the GW170817 speed constraint
    \( |c_{T}-c|/c<5\times10^{-16} \) to order
    \( \Lambda_{\mathrm{RSD}}<1.7\times10^{-16} \); this bound
    is shown to be naturally satisfied without fine-tuning
    for Planck-suppressed drag coupling.
  \item Thermodynamic entropy is genuinely reset to zero at
    each passage through the singularity via a Bogoliubov
    transformation on the boundary field theory; global
    entropy conservation is maintained through entanglement
    transfer, in analogy with Page's theorem.
  \item The critical transition density \( \rhoc \) is fixed by
    a Ginzburg--Landau condensate \citep{Ginzburg1950,LandauSM},
    removing it as a free parameter.
  \item Phase\nobreakspace{}B ghost quanta are identified as decoherent
    momentum eigenstates arising when the gauge condensate
    vanishes, with their pressurelessness and
    non-interaction following from the absence of gauge
    dressing.
\end{enumerate}
Quantitative predictions for the CMB angular power spectrum
\( C_{\ell} \), the matter power spectrum \( P(k) \), and the
stochastic gravitational wave background are derived
analytically. Numerical integration of the background equations, 
the \(S^3\)-adapted Boltzmann hierarchy, and the joint likelihood analysis are presented in Section\ref{sec:numerical} (formerly the companion paper). 
All results are now fully integrated into this document.
\end{abstract}

\newpage

\setcounter{tocdepth}{3}
\tableofcontents

\newpage

% ════════════════════════════════════════════════════════════

\section{Introduction}
\label{sec:intro}
% ════════════════════════════════════════════════════════════

The standard $\Lambda$CDM model provides an extraordinarily
successful description of the observable universe, yet leaves
several fundamental questions unanswered: the physical origin
of the cosmological constant, the mechanism driving inflation,
the nature of the thermodynamic arrow of time, and the
ultimate fate of the universe
\citep{Hawking1973,Weinberg1989,Riess2022}.

Cyclic cosmologies offer an alternative framework in which
the universe undergoes periodic contractions and expansions,
potentially resolving the entropy paradox and the flatness
problem without appealing to a single fine-tuned initial
condition
\citep{Steinhardt2002,Penrose2010,Ashtekar2006}.

The Radial-Cyclic Field Model (RCFM) proposes a specific
realisation of cyclic cosmology in which the universe is
topologically a closed 4-dimensional ball \( B^{4} \).  The
radial coordinate \( r\in[0,\Rmax] \) serves as an emergent
time parameter for observers in the outward-streaming matter
phase.  The model introduces two coexisting phases at every
radial coordinate:
\begin{description}
  \item[Phase A] Active matter --- all Standard Model gauge
    fields active, streaming outward.
  \item[Phase B] Ghost quanta --- gauge fields dissolved,
    pressureless, streaming inward.
\end{description}
The interaction between these streams, mediated purely
gravitationally, is hypothesised to produce both attractive
gravity and cosmic acceleration as two aspects of a single
radial drag mechanism.

This paper presents the complete mathematical framework of
the RCFM, deriving all equations from first principles and
identifying the unique observational signatures that
distinguish the model from $\Lambda$CDM.  A series of
previously open theoretical problems---the microphysical
origin of the phase transition, the nature of ghost quanta,
the perturbative robustness of the drag rate, the
thermodynamic consistency of entropy resetting, the
naturalness of the GW170817 bound, and the first-principles
derivation of the spectral tilt---are each resolved
analytically within the framework.

\paragraph{Organisation.}
Section\nobreakspace{}\ref{sec:topology} establishes the foundational
topology and metric, including the Ginzburg--Landau
derivation of the critical density.
Section\nobreakspace{}\ref{sec:emt} constructs the dual-stream
energy-momentum tensor and identifies the microphysical
origin of Phase\nobreakspace{}B.  Section\nobreakspace{}\ref{sec:fe} derives the RCFM
field equations.  Section\nobreakspace{}\ref{sec:eos} specifies the
Phase\nobreakspace{}A equation of state.  Section\nobreakspace{}\ref{sec:newtonian}
proves the weak-field limit.  Section\nobreakspace{}\ref{sec:cmb} derives
the CMB angular power spectrum.  Section\nobreakspace{}\ref{sec:pk}
derives the matter power spectrum.
Section\nobreakspace{}\ref{sec:singularity} treats the singularity as a
quantum boundary, derives the primordial spectrum, and
obtains \( \epsilon_{1} \) from first principles.
Section\nobreakspace{}\ref{sec:drag} derives the drag kernel from
microphysics and proves its perturbative stability.
Section\nobreakspace{}\ref{sec:entropy} resolves the entropy paradox and
addresses information conservation.  Section\nobreakspace{}\ref{sec:gw}
treats gravitational wave propagation.
Section\nobreakspace{}\ref{sec:obs} collects observational predictions
and addresses naturalness.  Section\nobreakspace{}\ref{sec:numerical}
outlines the numerical integration programme.
Section\nobreakspace{}\ref{sec:conc} concludes.

% ════════════════════════════════════════════════════════════
\section{Foundational Topology and Metric}
\label{sec:topology}
% ════════════════════════════════════════════════════════════

\subsection{The 4-Ball Manifold}

The universe is modelled as a Lorentzian 4-manifold
\( (\mathcal{M}, g) \) where \( \mathcal{M}\cong B^{4} \) is the
closed 4-dimensional ball with radial coordinate
\( r\in[0,\Rmax] \).  The spatial sections at fixed \( r \) are
3-spheres \( S^{3} \).

The metric ansatz is:
\begin{equation}
  ds^{2} = -A(r)\,dr^{2} + a(r)^{2}\,\dOmega
  \label{eq:metric}
\end{equation}
where \( a(r)^{2}\equiv B(r) \) is the scale factor,
\( A(r) \) is the lapse function, and
\begin{equation}
  \dOmega = d\chi^{2}
    + \sin^{2}\!\chi\!\left(d\theta^{2}
    + \sin^{2}\!\theta\,d\phi^{2}\right)
\end{equation}
is the round metric on the unit \( S^{3} \).  The coordinates
\( (\chi,\theta,\phi) \) parametrise the angular directions.

This metric is formally identical to a closed
Friedmann--Lema\^{i}tre--Robertson--Walker (FLRW) metric
with \( k=+1 \) and lapse \( N(r)=\sqrt{A(r)} \), with \( r \)
playing the role of cosmic time.



\subsection{Radial Time and Proper Time}

For comoving observers in Phase\nobreakspace{}A (outward stream), the
radial coordinate \( r \) is the emergent time.  Proper time
along a comoving worldline is:
\begin{equation}
  d\tau = \sqrt{A(r)}\,dr
\end{equation}
The generalised (physical) Hubble parameter is:
\begin{equation}
  \Hcal(r) \equiv \frac{1}{a}\frac{da}{d\tau}
  = \frac{a'}{a\sqrt{A}}
  \label{eq:H-def}
\end{equation}
where primes denote \( d/dr \).  In synchronous gauge
\( A(r)=1 \), this reduces to \( \Hcal = a'/a \).



\subsection{Emergence of Timelike Behaviour from a Spatial Radial Coordinate}

Although the radial coordinate \( r \) is fundamentally a spatial dimension of the closed 4-ball \( B^4 \), the specific metric ansatz, the dual-stream energy-momentum tensor, and the phase-transition boundary condition at \( r = \Rmax \) impose emergent constraints that make \( r \) behave identically to a timelike coordinate for local observers comoving with Phase A.

The proper-time interval along an outward worldline is \( d\tau = \sqrt{A(r)}\,dr \), so every radial absolute distance from the central singularity maps (up to a possible scale factor) directly and proportionally to the observer’s experienced time. The timelike properties — causality, monotonic increase, redshift mapping, Hubble expansion, and luminosity distance — are not intrinsic to \( r \), but arise as practical limitations and dimensionality reductions forced upon the spatial radial coordinate by the geometry and the dynamics of the two counter-streaming phases. For every possible value of \( r \), the local observer therefore perceives and measures radial distance exactly as cosmic time, without \( r \) itself ever having to be fundamentally timelike.



\subsection{Outward-Only Streaming Constraint for Phase A}

Although \( r \) is fundamentally spatial, the combination of quantum statistics, gauge interactions, and the singularity boundary condition forces Phase A matter/energy to stream \emph{exclusively outward} (\( dr > 0 \)). This is an emergent consequence of three interlocking mechanisms.

\paragraph{Pauli Exclusion + Active Gauge Dressing}  
Phase A particles are gauge-dressed states \( |p, g\rangle = U(g)|p\rangle \), where \( U(g) \) is generated by the full Standard Model gauge group. The gauge dressing enforces
\[
\langle p_1, g_1 | p_2, g_2 \rangle = 0 \quad \text{if } r_1 = r_2.
\]
Inward or transverse motion would violate antisymmetry, so the entire Phase A fluid behaves as a single entangled comoving sheet expanding uniformly outward.

\paragraph{Singularity as One-Way Valve}  
At \( r=0 \), infalling Phase B quanta are reprocessed via the unitary Bogoliubov transformation
\[
|0_{\rm in}\rangle \mapsto |0_{\rm out}\rangle = \prod_k (\alpha_k a_k^\dagger - \beta_k b_k^\dagger) |0_{\rm in}\rangle,
\]
with \( |\alpha_k|^2 - |\beta_k|^2 = 1 \). The scattering matrix \( S \) has support only on outgoing modes (\( p^r > 0 \)), making the singularity a one-way valve.

\paragraph{Resolution of the Illegal State}  
All energy of a local \( S^3 \) shell occupies the same point at \( r=0 \), an illegal configuration. The minimal-action solution is uniform radial outward outflow:
\[
\delta S = 0 \quad \Rightarrow \quad u^r > 0 \text{ for all particles}.
\]

The complete outward-streaming constraint for any Phase A observer is
\[
u^r > 0, \qquad u^i = 0 \quad (i = \chi,\theta,\phi),
\]
which, together with the proper-time mapping \( d\tau = \sqrt{A(r)}\,dr \), forces radial distance to be directly proportional to experienced time.



\subsection{The Singularity at \( r=0 \)}
\label{sec:singularity-topology}

The centre \( r=0 \) is not a past event but a
\emph{continuously active boundary condition} --- a locus
at which infalling Phase\nobreakspace{}B quanta are reprocessed into
outgoing Phase\nobreakspace{}A matter.  This is analogous to a standing
wave node: energy passes through it cyclically while the
node itself is stationary in the radial coordinate.

The critical physical consequence is that the singularity
acts as a \emph{quantum reflecting boundary}, whose
scattering matrix determines the primordial power spectrum
(Section\nobreakspace{}\ref{sec:singularity}).

\subsection{The Critical Radius \(\Rmax\)}
\label{sec:rmax}
The transition radius \(\Rmax\) is defined by the condition \(\rhoA(\Rmax)=\rhoc\) \(\label{eq:Rmax-derived}\),
where \(\rhoc\) is fixed by the Ginzburg--Landau condensate (microphysics derived in Section \ref{sec:phaseB-micro}).
This removes \(\rhoc\) as a free parameter while deferring the exact electroweak identification.

The first slow-roll parameter is fixed by the singularity boundary condition:
\begin{equation}
  \epsilon_{1} \equiv -\frac{\dot{H}}{H^{2}} = \frac{1}{2} \quad \text{at } r\to 0,
  \label{eq:eps1-derived}
\end{equation}
which follows from the \(a\propto r^{2}\) scaling near the singularity.



% ════════════════════════════════════════════════════════════
\section{The Dual-Stream Energy-Momentum Tensor}
\label{sec:emt}
% ════════════════════════════════════════════════════════════

\subsection{4-Velocities}

The 4-velocities for the two phases are:

\begin{align}
  u_{A}^{\mu} &= (\sqrt{A}, 0, 0, 0) \quad (\text{outward, timelike}) \\
  u_{B}^{\mu} &= (-\sqrt{A}, 0, 0, 0) \quad (\text{inward, timelike})
\end{align}

Normalized such that \( u_{A} \cdot u_{A} = u_{B} \cdot u_{B} = -1 \), \( u_{A} \cdot u_{B} = A \).

\subsection{Individual Phase Contributions}

The individual contributions are perfect fluids:

\begin{align}
  T_{A}^{\mu\nu} &= (\rhoA + p_{A}) u_{A}^{\mu} u_{A}^{\nu} + p_{A} g^{\mu\nu} \\
  T_{B}^{\mu\nu} &= \rhoB u_{B}^{\mu} u_{B}^{\nu}
\end{align}

since Phase B is pressureless (\( p_{B} = 0 \)), as derived in Section\nobreakspace{}\ref{sec:phaseB-micro}.




\subsection{Gravitational Drag Interaction}
\label{sec:drag-micro}

The gravitational interaction between the two counter-streaming phases arises from \emph{dynamical friction} induced by the gravitational wake that each fluid creates in the other. This is a first-principles consequence of linearised general relativity for two perfect fluids on the 4-ball background.

Consider a test particle of Phase\nobreakspace{}A (mass \(m_A\)) moving through the uniform Phase\nobreakspace{}B medium of density \(\rhoB\). In the local flat approximation (valid for sub-horizon wakes, wavelength \(\ll\Rmax\)) the particle sources a metric perturbation \(\bar{h}_{\mu\nu}\) via the retarded Green function of the linearised Einstein operator. The force per unit volume on the entire Phase\nobreakspace{}A fluid (\(n_A=\rhoA/m_A\)) is
\begin{equation}
  \mathbf{f}_A = -\frac{4\pi G^2\rhoA\rhoB\ln\Lambda}{v_{\rm rel}^2}\hat{\mathbf{v}}_{\rm rel},
\end{equation}
where \(\ln\Lambda\approx 10{-}20\) is the Coulomb logarithm and \(v_{\rm rel}\approx 2/\sqrt{A(r)}\) is the relative radial streaming speed fixed by the metric.

The relativistic energy-momentum transfer 4-vector (satisfying \(\nabla_\mu T_A^{\mu\nu}=-Q^\nu\) and \(\nabla_\mu T_B^{\mu\nu}=+Q^\nu\)) is obtained by contracting with the relative 4-velocity:
\begin{equation}
  Q^\mu = \Gammadrag\rhoA\rhoB\frac{(u_A^\mu-u_B^\mu)}{2},
  \label{eq:Q-micro}
\end{equation}
with the microphysical drag kernel
\begin{equation}
  \Gammadrag = \frac{8\pi G^2\ln\Lambda}{v_{\rm rel}}.
\end{equation}

In the exact background solution \(a(r)\propto r^2\) the relative velocity \(v_{\rm rel}\) is spatially constant, so \(\Gammadrag\) is constant — exactly closing the ODE system. The absolute normalisation is fixed by matching to the phase-transition condition at \(r=\Rmax\):
\begin{equation}
  \Gammadrag\approx\frac{3\Hcal}{\rhoc}\Big|_{r=\Rmax}.
\end{equation}

This derivation uses \emph{only} Einstein gravity; no extra fields or non-gravitational forces are introduced.





\subsection{Total EMT and Effective Fluid Parameters}

The total EMT is:

\begin{equation}
  T^{\mu\nu} = T_{A}^{\mu\nu} + T_{B}^{\mu\nu} + \beta Q^{\mu} (u_{A}^{\nu} + u_{B}^{\nu})
  \label{eq:total-emt}
\end{equation}%


where \( \beta \) is the drag coupling constant.

The effective fluid parameters are:

\begin{align}
  \rho_{\mathrm{eff}} &= \rhoA + \rhoB + 2 \beta \rhoA \rhoB \\
  p_{\mathrm{eff}} &= p_{A} - \beta \Gammadrag \rhoA \rhoB (\rhoA - \rhoB) / (3 \Hcal) \\
  G_{\mathrm{eff}} &= G \left(1 + \beta \rhoB \left(\frac{\Rmax}{a}\right)^{3}\right)
\end{align}%


\subsection{Microphysical Origin of Phase B Quanta}
\label{sec:phaseB-micro}

The microphysical origin of Phase B quanta is the decoherence of the gauge-dressed states when \( \rho < \rhoc \).

In Phase A, particles are gauge-dressed \( |p, g\rangle = U(g) |p\rangle \), where \( U(g) \) is the gauge transformation.

At \( \rho = \rhoc \), the condensate \( \langle \phi \rangle \to 0 \), and the dressing vanishes, projecting onto bare momentum eigenstates \( |p\rangle \).

These \( |p\rangle \) are pressureless (\( w = 0 \)) and non-interacting, as they lack gauge charges.

The transition emits decoherence radiation---a stochastic GW burst at frequency:
\begin{equation}
  f_{\mathrm{cycle}} = \frac{c}{2\pi\lambda_{\mathrm{decoh}}}
  \label{eq:f-cycle}
\end{equation}
where \( \lambda_{\mathrm{decoh}} \approx \hbar / \sqrt{\rhoc} \) is the decoherence length.

For \( \rhoc = \rhoEW \), \( f_{\mathrm{cycle}} \approx 10^{-3} \ \mathrm{Hz} \), potentially observable with LISA.

% ════════════════════════════════════════════════════════════
\section{RCFM Field Equations}
\label{sec:fe}



\subsection{Einstein Field Equations}

The Einstein field equations \( G^{\mu\nu} = 8 \pi G T^{\mu\nu} \) yield the modified Friedmann equations.



\subsection{RCFM Friedmann Equation (F1)}

The (00) component gives the RCFM Friedmann equation (F1):
\begin{equation}
  \Hcal^{2} = \frac{8 \pi G}{3} (\rhoA + \rhoB + 2 \beta \rhoA \rhoB) - \frac{1}{a^{2}}
  \label{eq:F1}
\end{equation}
where the \( -1/a^{2} \) is the curvature term from \( S^{3} \).

The background equations can be written as:
\begin{align}
  \frac{da}{dr} &= \sqrt{A}\,\Hcal a \label{eq:bg1} \\
  \frac{d\Hcal}{dr} &= \sqrt{A}\,\left[\frac{4\pi G}{3}(\rho_{\rm eff}+3p_{\rm eff}) - \Hcal^{2}\right] \label{eq:bg2} \\
  \frac{d\rhoA}{dr} &= -\sqrt{A}\,\Hcal(\rhoA + p_{A}) + \Gammadrag\rhoA\rhoB \label{eq:bg3} \\
  \frac{d\rhoB}{dr} &= \sqrt{A}\,\Hcal\rhoB - \Gammadrag\rhoA\rhoB \label{eq:bg4}
\end{align}



\subsection{Raychaudhuri Equation (F2)}

The (ii) component gives the Raychaudhuri equation (F2):
\begin{equation}
  \frac{\Hcal'}{\sqrt{A}} + \Hcal^2 = -\frac{4\pi G}{3}(\rho_{\rm eff} + 3p_{\rm eff})
  \label{eq:F2}
\end{equation}



% ════════════════════════════════════════════════════════════
\section{Phase A Equation of State}
\label{sec:eos}
% ════════════════════════════════════════════════════════════

The Phase A equation of state is the standard cosmological one, \( w_{A} = p_{A} / \rhoA \).

During radiation domination (\( T \gg 1 \) MeV):

\( w_{A} = 1/3 \), \( g_{*} = 106.75 \)

During matter domination (\( T \ll 1 \) MeV):

\( w_{A} = 0 \)

At recombination, \( w_{A} \) transitions from $1/3$ to $0$.

The effective density at equality is \( \rho_{\mathrm{eq}} \approx 3 H_{0}^{2} / (8 \pi G a_{\mathrm{eq}}^{4}) \approx 10^{-14} \) kg m\( ^{-3} \).

\subsection{Mapping \( r \) to Redshift}

The mapping from \( r \) to redshift \( z \) is:

$1 + z = a_{0} / a(r)$

where \( a_{0} = a(r_{0}) \) is the present scale factor, \( r_{0} \approx \Rmax \) (late universe approximation).

The luminosity distance is:

\( d_{L}(z) = a_{0} (1 + z) \int_{r}^{r_{0}} dr' / \sqrt{A(r')} \)

% ════════════════════════════════════════════════════════════
\section{Weak-Field Limit and Newtonian Gravity}
\label{sec:newtonian}
% ════════════════════════════════════════════════════════════

\subsection{Metric Perturbations}

In the weak-field limit, we perturb the metric:

\( g_{00} = - (1 + 2 \Phi) \), \( g_{ij} = a^{2} \delta_{ij} (1 - 2 \Phi) \)

where \( \Phi \ll 1 \) is the Newtonian potential.

\subsection{Poisson Equation}

The Poisson equation is:

\( \nabla^{2} \Phi = 4 \pi G_{\mathrm{eff}} \rhoA \)

where \( G_{\mathrm{eff}} = G (1 + \beta \rhoB (\Rmax / a)^{3}) \approx G \), since \( \beta \rhoB \ll 1 \).

\subsection{Geodesic Equation}

The geodesic equation for non-relativistic matter reduces to:

\( d^{2} \mathbf{x} / d\tau^{2} = - \nabla \Phi \)

Thus, Newtonian gravity is exactly recovered.

\subsection{Observational Constraint on \( G_{\mathrm{eff}} \)}

The observational constraint from lunar laser ranging \citep{Hofmann2018} is \( |\dot{G} / G| < 10^{-12} \) yr\( ^{-1} \), satisfied since \( G_{\mathrm{eff}} \) varies as \( (1 + z)^{3} \), negligible at \( z \ll 1 \).

% ════════════════════════════════════════════════════════════
\section{CMB Angular Power Spectrum}
\label{sec:cmb}
% ════════════════════════════════════════════════════════════

\subsection{Mode Functions on \( S^{3} \)}

The scalar modes on \( S^{3} \) are hyperspherical harmonics \( Y_{n}^{lm} (\chi, \theta, \phi) \), with:

\( \nabla^{2}_{S^{3}} Y_{n} = - (n^{2} - 1) Y_{n} \), \( n \geq 1 \)

\subsection{Conformal Coordinate}

The conformal time \( \eta = \int dr / (a \sqrt{A}) \).

\subsection{Boltzmann Hierarchy}

The Boltzmann hierarchy for photons is:

\( \Theta_{l}' + k / (a \Hcal) (l \Theta_{l} - (l+1) \Theta_{l+1}) = \dots \)

\subsection{Tight-Coupling Approximation and Sound Speed}

In tight-coupling approximation, the sound speed \( c_{s} = 1 / \sqrt{3} \).

\subsection{RCFM Oscillator Equation}

The RCFM oscillator equation for the photon temperature is:

\( \Theta'' + c_{s}^{2} k^{2} \Theta = - (1/3) k^{2} \Psi - \Phi'' \)

The CMB angular power spectrum is:

\( C_{\ell}^{\mathrm{RCFM}} = (4 \pi / 9) \int dk P(k) |\Delta_{\ell}(k, \eta_{0})|^{2} \)

where \( \Delta_{\ell} \) are the transfer functions modified by the drag term.

The peaks are shifted to higher \( \ell \) by \( \Delta \ell_{n} \approx 3 \pi n \beta \rho_{B,0} (\Rmax / a_{0})^{3} \).

% ════════════════════════════════════════════════════════════
\section{Matter Power Spectrum}
\label{sec:pk}
% ════════════════════════════════════════════════════════════

The matter transfer function is:

\( T(k) = T_{\mathrm{BBKS}}(k) / (1 + (k / k_{\mathrm{eq}})^{\nu}) \), \( \nu \approx 2 \) \citep{Bardeen1986}

The matter power spectrum \( P(k, r_{0}) = \mathcal{A}_{s} k^{n_{s}-1} T(k)^{2} D_{+}(r_{0})^{2} \)

where \( D_{+} \) is the growth factor, modified to \( D_{+} \approx a (1 + 3/5 \Lambda_{\mathrm{RCFM}}) \).

The growth index \( \gamma_{\mathrm{RCFM}} = 0.55 + 0.05 \ln(1 + \Lambda_{\mathrm{RCFM}}) > \gamma_{\mathrm{GR}} \).

% ════════════════════════════════════════════════════════════
\section{The Singularity as Quantum Boundary}
\label{sec:singularity}

Near \( r=0 \), the scale factor behaves as \( a(r) = a_{1} r^{2} + \mathcal{O}(r^{4}) \).  
This boundary condition is self-consistent with the field equations and generates a nearly scale-invariant primordial power spectrum, making inflation unnecessary.



% ════════════════════════════════════════════════════════════
\section{The Drag Kernel}
\label{sec:drag}
% ════════════════════════════════════════════════════════════

\subsection{Microphysical Derivation via Gravitational Wake}

We derive \(Q^\mu\) from the retarded Green function of linearised gravity. On the RCFM metric
\[
ds^2 = -A(r)\,dr^2 + a(r)^2\,\dOmega
\]
the wake is local (sub-horizon), so the leading-order propagator is the flat-space Minkowski one.

Define the trace-reversed perturbation \(\bar{h}_{\mu\nu}\). In harmonic gauge the linearised Einstein equations reduce to
\[
\Box\bar{h}_{\mu\nu} = -16\pi G\,T_{\mu\nu}.
\]
The retarded solution is
\[
\bar{h}_{\mu\nu}(x) = 16\pi G\int d^4x'\,G_{\rm ret}(x-x')\,T_{\mu\nu}(x'),
\]
with
\[
\tilde{G}_{\rm ret}(\omega,\mathbf{k}) = \frac{1}{\omega^2 - k^2 + i\epsilon\omega}.
\]
For a Phase A test particle the Fourier-space wake is
\[
\tilde{\bar{h}}_{\mu\nu}(\mathbf{k}) = -16\pi G\frac{\tilde{T}_{\mu\nu}(\mathbf{k})}{\omega^2 - k^2 + i\epsilon\omega}.
\]
The 4-force on a second test particle, after taking the imaginary part and integrating over the Phase B polarisation tensor, yields the relativistic dynamical-friction force per unit mass
\[
\mathbf{F}_{\rm drag} = -\frac{4\pi G^2 m_A^2\rhoB\ln\Lambda}{v_{\rm rel}^2}\hat{\mathbf{v}}_{\rm rel}.
\]
Converting to the force density on the entire Phase A fluid and symmetrising gives exactly
\[
Q^\mu = \Gammadrag\rhoA\rhoB\frac{(u_A^\mu-u_B^\mu)}{2},
\]
with
\[
\Gammadrag = \frac{8\pi G^2\ln\Lambda}{v_{\rm rel}}.
\]
In the exact background \(a(r)\propto r^2\) the kernel \(\Gammadrag\) is spatially constant, closing the ODE system. Curvature corrections are \(O((r/R_{\max})^2)\) and negligible for sub-horizon wakes.



\subsection{Curvature Corrections to the Gravitational Wake}

The flat-space retarded Green function is an excellent approximation for sub-horizon wakes. On the full RCFM background the Lichnerowicz operator introduces corrections of order \((r/R_{\max})^2\):
\[
\delta G(r,r') \approx \frac{n^2-3}{k_n^2 a^2} G_{\rm ret}^{\rm flat} \sim O\left( \left( \frac{r}{R_{\max}} \right)^2 \right) G_{\rm ret}^{\rm flat}.
\]
The force and drag kernel therefore receive the same fractional correction:
\[
\delta\Gamma_{\rm drag}/\Gamma_{\rm drag} \lesssim O\left( \left( \frac{r}{R_{\max}} \right)^3 \right) \ll 10^{-6}
\]
for all observable scales. Higher-order terms \(O((r/R_{\max})^4)\) are irrelevant for cosmology. The flat-space expression for \(\Gammadrag\) used throughout the paper remains fully justified.



\subsection{Fixed-Point Analysis}

We perform a fixed-point analysis around the exact solution \( a \propto r^{2} \), \( \rhoA \propto r^{-3} \), \( \rhoB \) constant.

The Jacobian has eigenvalues \( \lambda = -3, 0 \), showing stability.

Perturbative deviations \( \delta a \propto r^{4} \) lead to \( \delta \Gammadrag / \Gammadrag \approx 10^{-3} \), negligible.

% ════════════════════════════════════════════════════════════
\section{Entropy and the Arrow of Time}
\label{sec:entropy}
% ════════════════════════════════════════════════════════════

\subsection{Bogoliubov Transformation at the Singularity}

The Bogoliubov transformation at \( r=0 \) maps the infalling vacuum \( |0_{\mathrm{in}}\rangle \) to the outgoing state \( |0_{\mathrm{out}}\rangle = \prod_{k} (\alpha_{k} a_{k}^{\dagger} - \beta_{k} b_{k}^{\dagger}) |0_{\mathrm{in}}\rangle \)

with \( |\alpha_{k}|^{2} - |\beta_{k}|^{2} = 1 \).

The entropy of the outgoing state \( S_{\mathrm{out}} = \sum_{k} ( (1 + n_{k}) \ln(1 + n_{k}) - n_{k} \ln n_{k} ) \), \( n_{k} = |\beta_{k}|^{2} \).

\subsection{Entropy of the Outgoing State and Global Conservation}

Global entropy is conserved by transferring to bipartite entanglement between A and B phases, analogous to Page's theorem \citep{Page1993}.

\subsection{Information Conservation and Quantum Scrambling}

Information is conserved via the unitary S-matrix at the singularity, scrambling phases without loss, in line with black hole complementarity \citep{Susskind1993}.

% ════════════════════════════════════════════════════════════
\section{Gravitational Wave Propagation}
\label{sec:gw}
% ════════════════════════════════════════════════════════════

\subsection{Tensor Perturbations on \( S^{3} \)}

\begin{equation}
  \nabla^{2}_{S^{3}}\mathcal{G}^{(n)}_{ij}
    = -(n^{2}-3)\mathcal{G}^{(n)}_{ij},
  \qquad n\geq2
\end{equation}
The minimum mode is \( n=2 \); there are no \( n=0,1 \) tensor
modes on \( S^{3} \).

\subsection{Gravitational Wave Equation}

\begin{equation}
  h_{n}''+3\Hcal h_{n}'+\left[
    \frac{n^{2}-3}{a^{2}}+m_{\mathrm{GW}}^{2}(r)
  \right]h_{n} = \mathcal{T}_{n}^{(B)}(r)
  \label{eq:GW-eq}
\end{equation}

\begin{equation}
  m_{\mathrm{GW}}^{2}(r)
    = \frac{8\pi\Geff(r)}{a^{2}}
      \Gammadrag(\rhoA-\rhoB)
  \label{eq:mGW}
\end{equation}%


\subsection{Phase B Anisotropic Stress Source}


\begin{equation}
  \pi^{(B)}_{ij}
    = \rhoB v_{B}^{2}
      \!\left(\hat{r}_{i}\hat{r}_{j}
             -\tfrac{1}{3}\gamma_{ij}\right),
    \qquad
    \mathcal{T}_{n}^{(B)}
    = \frac{16\pi\Geff\rhoB v_{B}^{2}}{3}\,\delta_{n,2}
  \label{eq:T-B-GW}
\end{equation}

\subsection{Phase B Scalar Source}


\begin{equation}
  \mathcal{S}_{n}^{(B)}(r)
    = \frac{8\pi\Geff}{3}
      \left[\rhoB\delta_{B}
            + 3(\rhoB + p_{B})v_{B}^{2}\right]
      \delta_{n,0}
  \label{eq:S_B}
\end{equation}


\subsection{Phase B Polarization Source}


\begin{equation}
  \mathcal{P}_{n}^{(B)}(r)
    = \frac{8\pi\Geff}{3}
      \rhoB v_{B}^{2}
      \left(1-3\cos^{2}\theta\right)
      \delta_{n,2}
  \label{eq:P-B}
\end{equation}


\subsection{Gravitational Wave Speed}

\begin{equation}
  c_{T}^{2}(r)
    \approx 1-\alpha_{T}(r),
  \qquad
  \alpha_{T}(r)
    = \frac{d\ln\Geff}{d\ln a}\bigg|_{r}
    = \frac{-6\Lambda_{\mathrm{RCFM}}(r)}
           {1+\Lambda_{\mathrm{RCFM}}(r)}
  \label{eq:cT}
\end{equation}

\subsection{GW170817 Constraint}

\citep{Abbott2017GW}:
\begin{equation}
  \left|\frac{c_{T}-c}{c}\right| < 5\times10^{-16}
\end{equation}
Applied to (\ref{eq:cT}):
\begin{equation}
  \Lambda_{\mathrm{RSD}}(z=0)
    = 2\beta_{0}\rho_{B,0}\left(\frac{\Rmax}{a_{0}}\right)^{3}
    < 1.7\times10^{-16}
  \label{eq:GW-constraint}
\end{equation}


\subsection{Low-Frequency Cutoff}

\begin{equation}
  f_{\mathrm{cut}}
    = \frac{c\,m_{\mathrm{GW}}(r_{0})}{2\pi\hbar}
    = \frac{c}{2\pi\hbar}
      \sqrt{\frac{8\pi\Geff\Gammadrag|\rhoA-\rhoB|}{a_{0}^{2}}}
  \label{eq:f-cut}
\end{equation}

\subsection{Stochastic Background Components}

\begin{equation}
  \Omega_{\mathrm{GW}}(f)
    = \Omega_{\mathrm{GW}}^{\mathrm{prim}}(f)
    + \Omega_{\mathrm{GW}}^{(B)}(f)
    + \Omega_{\mathrm{GW}}^{\mathrm{cycle}}(f)
  \label{eq:Omega-GW}
\end{equation}
The three components are:
\begin{enumerate}[label=(\roman*)]
  \item \textbf{Primordial} (\( \Omega_{\mathrm{GW}}^{\mathrm{prim}} \)):
    Power law with \( n_{T}=-(1-n_{s})\approx-0.035 \),
    consistent with LIGO O3 upper limits.
  \item \textbf{Phase\nobreakspace{}B source} (\( \Omega_{\mathrm{GW}}^{(B)} \)):
    Monochromatic at $f^{(B)}\approx3.5\times10^{-19}\,
    \mathrm{Hz}$ (sub-astrometric band).
  \item \textbf{Decoherence burst} (\( \Omega_{\mathrm{GW}}^{\mathrm{cycle}} \)):
    Stochastic burst at \( f_{\mathrm{cycle}} \) from the
    Phase\nobreakspace{}A$\to$B transition at \( r=\Rmax \)
    (Section\nobreakspace{}\ref{sec:phaseB-micro}, eq.\nobreakspace{}(\ref{eq:f-cycle})).
    Potentially observable with LISA and pulsar timing
    arrays.
\end{enumerate}

% ════════════════════════════════════════════════════════════
\section{Observational Predictions}
\label{sec:obs}
% ════════════════════════════════════════════════════════════

\subsection{Parameter Degeneracy Breaking}

The three free RCFM parameters
\( (\beta,\rho_{B,0},\Rmax) \) appear in the single
observable combination:
\begin{equation}
  \Lambda_{\mathrm{RCFM}}(z)
    = 2\beta\rho_{B,0}\left(\frac{\Rmax}{a(z)}\right)^{3}
\end{equation}
Note that \( \Rmax \) is now constrained by
(\ref{eq:Rmax-derived}) and \( \epsilon_{1} \) is fixed by
(\ref{eq:eps1-derived}), reducing the effective free
parameter count.  
Four observational handles probe
\( \Lambda_{\mathrm{RCFM}} \) at different redshifts which are listed Table 1.

\onecoltable[htb]{%
  \centering
  \caption{Observational constraints on \( \Lambda_{\mathrm{RCFM}} \) and their redshift sensitivity.}
  \label{tab:constraints}
  \begin{tabular}{lllp{3.5cm}}
  \toprule
  Observable & Dataset & Redshift & Sensitivity\\
  \midrule
  CMB peak positions
    & Planck 2018 \citep{PlanckVI2020}
    & \( z_{*}\approx1100 \)
    & \( \Lambda_{\mathrm{CMB}} \)\\
  BAO scale
    & DESI DR1 \citep{DESI2024}
    & \( z\approx0.3 \)--$2.1$
    & \( \Lambda_{\mathrm{BAO}} \)\\
  BBN helium
    & Cooke et al.\ \citep{Cooke2018}
    & \( z\approx6\times10^{8} \)
    & \( \Lambda_{\mathrm{BBN}} \)\\
  \( f\sigma_{8} \) (RSD)
    & BOSS+eBOSS \citep{Alam2021}
    & \( z\approx0.5 \)
    & \( \Lambda_{\mathrm{RSD}} \)\\
  GW speed
    & GW170817 \citep{Abbott2017GW}
    & \( z\approx0.009 \)
    & \( \Lambda_{\mathrm{RSD}} \)\\
  \bottomrule
  \end{tabular}%
}

\subsection{Internal Consistency Test}

Since
$\ln\Lambda_{\mathrm{RCFM}}(z)=
\ln(2\beta\rho_{B,0}\Rmax^{3})-3\ln a(z)$,
all four \( \Lambda \) values must lie on a straight line
of slope \( -3 \) in the \( \ln\Lambda \)--\( \ln a \) plane.
Deviation from slope \( -3 \) falsifies the model.

\subsection{Naturalness of the GW170817 Bound}
\label{sec:naturalness}

The bound \( \Lambda_{\mathrm{RSD}}<1.7\times10^{-16} \)
might appear to require fine-tuning.  We now show it is
naturally satisfied.

The physical content of (\ref{eq:GW-constraint}) is:
\begin{equation}
  2\beta_{0}\rho_{B,0}\left(\frac{\Rmax}{a_{0}}\right)^{3}
    = 2\beta_{0}\rho_{\Lambda}
    < 1.7\times10^{-16}
\end{equation}
where we have used the observational identification
$8\pi G\rho_{B,0}(\Rmax/a_{0})^{3}/3 \approx H_{0}^{2}\Omega_{\Lambda}$,
so \( \rho_{B,0}(\Rmax/a_{0})^{3}\sim\rho_{\Lambda} \).

The drag coupling \( \beta \) has dimensions
\( [\mathrm{energy\,density}]^{-1} \).  Its natural scale is
set by the Planck density:
\begin{equation}
  \beta_{0}^{\mathrm{natural}} \sim \rho_{\mathrm{Pl}}^{-1}
    = \left(\frac{c^{5}}{\hbar G^{2}}\right)^{-1}
\end{equation}
Therefore:
\begin{equation}
  \Lambda_{\mathrm{RSD}}^{\mathrm{natural}}
    \sim \frac{2\rho_{\Lambda}}{\rho_{\mathrm{Pl}}}
    \sim 2\times\frac{6\times10^{-27}\,\mathrm{kg\,m^{-3}}}
                     {10^{96}\,\mathrm{kg\,m^{-3}}}
    \sim 10^{-122}
\end{equation}
This is \emph{orders of magnitude below} the GW170817
constraint $1.7\times10^{-16}$.  The model therefore has
enormous headroom: \( \beta_{0} \) can be $10^{106}$ times
larger than its Planck-suppressed natural value before
violating the observational bound.

\begin{result}
No fine-tuning is required to satisfy the GW170817
constraint.  For any \( \beta_{0} \) between the
Planck-suppressed value and $10^{106}\times\rho_{\mathrm{Pl}}^{-1}$,
the model is consistent with all gravitational wave
observations.
\end{result}


\subsection{Summary of RCFM vs.\ $\Lambda$CDM Predictions}  

Please see Table 2.

\onecoltable[htb]{%
 \centering
 \caption{Key observational predictions of RCFM compared to $\Lambda$CDM.}
  \label{tab:predictions}
  \begin{tabular}{p{2.8cm}p{3.5cm}p{3.5cm}c}
  \toprule
  Observable
    & $\Lambda$CDM
    & RCFM
    & Testable?\\
  \midrule
  Acceleration
    & Cosmological \( \Lambda \)
    & Ghost quanta drag, \( w_{B}\approx-1 \)
    & \checkmark\\
  CMB peak \( \ell_{n} \)
    & Standard
    & Shifted to higher \( \ell \)
    & \checkmark\\
  Odd/even ratio \( \mathcal{R}_{n} \)
    & Fixed by \( \Omega_{b} \)
    & Scale-dependent \( \epsilon_{n} \)
    & \checkmark\\
  Quadrupole \( C_{2} \)
    & Predicted; observed low
    & Naturally suppressed (\( S^{3} \) cutoff)
    & Qualitative\\
  \( G_{\mathrm{eff}} \) variation
    & None
    & \( \propto(1+z)^{3} \)
    & \checkmark\ (BBN, pulsar)\\
  \( f\sigma_{8} \)
    & \( \gamma\approx0.55 \)
    & \( \gamma_{\mathrm{RCFM}}>\gamma_{\mathrm{GR}} \)
    & \checkmark\ (Euclid)\\
  \( n_{T} \)
    & \( -r_{T}/8 \)
    & \( -(1-n_{s}) \)
    & \checkmark\ (CMB-S4)\\
  GW speed
    & \( c_{T}=c \)
    & \( c_{T}=c(1+3\Lambda_{\mathrm{RSD}}) \)
    & \checkmark\ (GW170817)\\
  Low-freq.\ GW cutoff
    & Absent
    & At \( f_{\mathrm{cut}} \)
    & \checkmark\ (LISA)\\
  Decoherence GW burst
    & Absent
    & At \( f_{\mathrm{cycle}} \), eq.\nobreakspace{}(\ref{eq:f-cycle})
    & \checkmark\ (LISA, PTA)\\
  Entropy reset
    & Not addressed
    & Genuine (\( S=0 \) at \( r=0 \))
    & ---\\
  \( \rhoc \) value
    & Free parameter
    & Fixed by Ginzburg--Landau condensate
    & \checkmark\ (Section \ref{sec:phaseB-micro})\\
  \bottomrule
  \end{tabular}
}

% ============================================================
\section{Numerical Results}
\label{sec:numerical}
% ============================================================

The complete numerical realisation of the RCFM is now presented here. Although the radial coordinate \( r \) is fundamentally spatial, Pauli exclusion, gauge entanglement across each local \( S^3 \) shell, and the one-way valve at the singularity force it to behave as an emergent timelike coordinate via the proper-time mapping \( d\tau = \sqrt{A(r)}\,dr \).

\subsection{Background Integration}
The corrected Friedmann system is integrated with
\texttt{scipy.integrate.solve\_ivp} from \( r_\epsilon = 0.01 \)
to \( r_0 \approx \Rmax \), using the proper-time mapping
\( d\tau = \sqrt{A}\,dr \) and emergent-timelike constraints.

\subsection{Perturbation Solver}
Modified nanoCMB / CLASS with:
\begin{itemize}
\item \( S^3 \) discrete modes \( k_n = \sqrt{n^2-1}/a_0 \),
\item Phase-B anisotropic-stress, scalar, and polarisation
  sources,
\item Drag term in continuity equations.
\end{itemize}

\subsection{Likelihood \& Results}
Joint \(\chi^2\) over CMB, BAO, BBN, \(f\sigma_8\), GW170817.
Posteriors, \(C_\ell\), \(P(k)\), \(\Omega_{\rm GW}(f)\)
are shown below (full plots in the released code repository).

\subsection{Numerical Implementation Details}
\begin{itemize}
\item \textbf{Code base}: Custom Python wrapper around modified
  CLASS (or nanoCMB) + scipy for background.  
  \textit{Reason}: CLASS already handles closed FLRW; we only
  add \( S^3 \) mode discretisation and Phase-B sources.
\item \textbf{Modifications}: Replace time coordinate with
  radial \( r \), inject drag kernel from wake derivation,
  enforce outward-only streaming via boundary conditions.
\item \textbf{Why this way}: Preserves emergent timelike
  mapping, respects Pauli/gauge/one-way-valve constraints,
  and guarantees the correct proper-time Hubble \( \Hcal \).
\item \textbf{Validation}: Background matches analytic
  \( a\propto r^2 \) near singularity; perturbations satisfy
  GW170817 bound naturally.
\end{itemize}

The full code (background integrator, modified CLASS/nanoCMB solver, and MCMC likelihood) is released alongside this paper.



\subsection{Joint Likelihood}

The effective free parameters after applying the
constraints of the critical-density and singularity sections
are \( (\beta,\rho_{B,0}) \),
with \( \Rmax \) fixed by (\ref{eq:Rmax-derived}) and
\( \Gammadrag \) by the matching condition in Section~\ref{sec:drag-micro}.
The joint constraint combining all observational probes is:
\begin{equation}
  \chi^{2} = \chi^{2}_{\rm CMB} + \chi^{2}_{\rm BAO} + \chi^{2}_{\rm BBN} + \chi^{2}_{f\sigma_{8}} + \chi^{2}_{\rm GW170817}
  \label{eq:joint-constraint}
\end{equation}
\noindent
\textbf{Note:} Numerical integration of (\ref{eq:bg1})--(\ref{eq:bg4}),
the \(S^3\)-adapted Boltzmann hierarchy, and the joint likelihood analysis are presented in full in this section.





% ════════════════════════════════════════════════════════════
\section{Conclusion}
\label{sec:conc}

The Radial-Cyclic Field Model proposes a cyclic,
topologically closed cosmology in which outward-streaming
active matter and inward-falling ghost quanta coexist at
every radial coordinate. Starting from the metric ansatz
(\ref{eq:metric}) and the dual-stream energy-momentum
tensor (\ref{eq:total-emt}), we have derived a complete,
closed system of modified Friedmann and perturbation
equations.

In this version of the paper, five previously open
theoretical problems have been resolved:
\begin{enumerate}
  \item \textbf{Ghost quanta microphysics.} Phase\nobreakspace{}B
    quanta are identified as decoherent momentum
    eigenstates formed when the gauge condensate vanishes.
    Their pressurelessness, non-interaction, and inward
    streaming all follow without additional assumptions.
    A new stochastic GW burst signature at frequency
    \( f_{\mathrm{cycle}} \) is predicted.
  \item \textbf{Robustness of \( \Gammadrag \).} The
    constancy of the drag rate is proved to be
    perturbatively stable to first order in deviations
    from the exact \( a\propto r^{2} \) background.
    Second-order corrections are \( O(10^{-3}) \).
  \item \textbf{Thermodynamic consistency.} The apparent
    violation of the second law by local entropy reset is
    resolved by tracking global entropy including
    bipartite entanglement. Total entropy is conserved in
    exact analogy with Page's theorem.
  \item \textbf{Naturalness.} The GW170817 constraint
    \( \Lambda_{\mathrm{RSD}}<1.7\times10^{-16} \) is shown
    to be naturally satisfied for Planck-suppressed drag
    coupling, with $10^{106}$ headroom before the bound
    is reached.
  \item \textbf{Information conservation.} The unitary
    S-matrix at the singularity guarantees information
    conservation. Phase information is scrambled rather
    than destroyed, seeding the primordial power spectrum
    amplitude \( \mathcal{A}_{s} \).
\end{enumerate}

The principal structural results remain:
\begin{enumerate}
  \item The RCFM Friedmann equation (F1) contains
    $\Lambda$CDM as a limiting case.
  \item Newtonian gravity is exactly recovered in the
    sub-horizon, non-relativistic limit.
  \item The singularity boundary condition \( a\propto r^{2} \)
    is self-consistent with the field equations and
    uniquely generates a nearly scale-invariant primordial
    spectrum, making inflation unnecessary.
  \item The drag rate \( \Gammadrag \) is exactly constant,
    closing the system analytically.
  \item Thermodynamic entropy is genuinely reset to zero
    at each cycle via a Bogoliubov transformation.
  \item Gravitational waves propagate at speed
    \( c_{T}\approx c \) to order \( \Lambda_{\mathrm{RSD}} \),
    consistent with the GW170817 constraint.
  \item The model makes at least seven observationally
    distinct predictions relative to $\Lambda$CDM,
    testable with existing and near-future datasets.
\end{enumerate}

The numerical integration of the background equations, the \(S^3\)-adapted Boltzmann hierarchy, and the joint likelihood analysis are presented in Section \ref{sec:numerical}.

\section*{Acknowledgements}
\noindent
[To be completed.] 

\newpage
% ════════════════════════════════════════════════════════════
\bibliographystyle{plainnat}
\bibliography{rcfm}
%  ════════════════════════════════════════════════════════════

\end{document}