Time Without Time: A Coherence-Rate Field Theory of Gravitation

Introduction

The operational meaning of time remains one of the least examined foundations of gravitational theory. General Relativity provides a consistent geometric relation between intervals, but the notion of temporal flow itself is external to the formalism: time is a coordinate, not a dynamical quantity. The link between physical clocks, their rates, and the structure of spacetime is formulated geometrically but not derived from a more primitive concept. This raises the question of whether time dilation and redshift can be expressed through a physically measurable rate field, rather than being introduced purely as geometric properties.

In parallel, several modern approaches—including relational and causal-set frameworks, thermodynamic derivations of field equations, and informational bounds on dynamical rates—suggest that temporal structure may have an underlying operational origin. These perspectives do not contradict General Relativity, but they motivate the search for formulations in which the rate of local physical processes plays the primary role. Within such a view, geometry encodes correlations between processes, while the processes themselves—their frequencies, coherence, and synchronization—constitute the more primitive layer.

In this work we develop a coherence-based formulation of gravitational phenomena. We introduce a scalar field \(\omega(x)\) characterizing the local rate at which physical systems accumulate proper time relative to a universal reference scale. This field is defined operationally through ratios of measurable clock frequencies. Rather than prescribing a metric a priori, we derive an effective geometry consistent with the definition of \(\omega(x)\) and show that gravitational redshift, time dilation, and the Newtonian limit emerge from variations in the coherence rate.

We emphasize that the coherence field \(\omega\) does not replace the tensorial degrees of freedom of General Relativity. The framework is not a scalar alternative to GR in the sense of Nordström-type models, but a scalar–rate extension that preserves the metric structure and reduces exactly to GR in the limit \(\omega\to 1\). The scalar degree of freedom acts as an operational layer encoding clock-rate modulation, rather than as a standalone gravitational field.

The goal is not to replace General Relativity but to provide a scalar-rate framework that is observationally equivalent in the weak-field regime and that offers a distinct underlying interpretation: spacetime structure reflects variations in local temporal rates. The theory is formulated through a covariant action for \(\omega(x)\), yielding a Klein–Gordon type equation with matter coupling fixed by Newtonian consistency. We compute the induced metric, derive parametrized post-Newtonian (PPN) parameters, and show agreement with General Relativity at first PPN order. Deviations are controlled by the gradients of \(\omega\) and by a characteristic coherence mass \(m_\omega\).

Finally, we outline a set of experimental signatures accessible to optical-clock comparisons, atom interferometry, pulsar timing arrays, moving-clock campaigns, and gravitational-wave phase measurements. These protocols provide falsifiable predictions and permit constraints on \(\omega(x)\) and its couplings. The approach thus unifies an operationally motivated view of time with a fully testable scalar-field modification of metric gravity.

Operational Definition of the Coherence Rate

The starting point of the framework is the recognition that gravitational redshift, time dilation, and clock-rate variation are fundamentally measurable through ratios of local clock frequencies. Any formulation in which time is treated operationally must begin with such ratios rather than with geometric intervals. This motivates the introduction of a scalar quantity characterizing the rate at which physical systems accumulate proper time relative to a universal reference scale.

Reference rate and upper bound

We postulate the existence of a universal reference temporal rate \(\nu_{\rm ref}\), representing the maximum synchronization rate at which local physical processes can be consistently compared. A natural choice for such a scale is the Planck frequency,

\[\nu_{\rm ref} \equiv \sqrt{\frac{c^5}{\hbar G}} \approx 1.855 \times 10^{43}\,\mathrm{Hz},\]

which provides an invariant bound derived entirely from the fundamental constants. No empirical assumption is required: \(\nu_{\rm ref}\) serves only as a reference against which operational ratios are defined.

Definition of the local coherence rate

Let \(\nu_{\rm loc}(x)\) denote the physical frequency of any sufficiently stable clock located at spacetime point \(x\). We define the coherence rate \(\omega(x)\) as the dimensionless ratio

\[\omega(x) \equiv \frac{\nu_{\rm loc}(x)}{\nu_{\rm ref}} = \frac{d\tau(x)}{dt_{\rm ref}},\]

where \(d\tau(x)\) is the proper-time increment registered by the local clock and \(dt_{\rm ref} \equiv \nu_{\rm ref}^{-1}\) is the reference interval. The definition is operational: \(\omega(x)\) can be obtained from any two clocks via

\[\frac{\omega(x_2)}{\omega(x_1)} = \frac{\nu_{\rm loc}(x_2)}{\nu_{\rm loc}(x_1)}.\]

Thus \(\omega(x)\) captures observable redshift relations without the need to specify a metric at this stage.

Relation to measured time dilation

In general relativity the gravitational redshift between two stationary clocks at positions \(x_1\) and \(x_2\) satisfies

\[\frac{\nu_{\rm loc}(x_2)}{\nu_{\rm loc}(x_1)} = \sqrt{\frac{g_{00}(x_2)}{g_{00}(x_1)}}.\]

Identifying \(\omega(x)\propto \nu_{\rm loc}(x)\), we obtain the correspondence

\[\omega(x) = \sqrt{g_{00}(x)} \qquad \text{(GR limit)}.\]

This shows that \(\omega(x)\) reproduces the observable content of \(g_{00}\) in stationary configurations. Departures from this relation encode possible non-metric or scalar-tensor effects.

Interpretation

The coherence rate \(\omega(x)\) should be thought of not as an additional physical clock but as a field summarizing the relative rates of all possible local clocks. Any theory that uses time operationally must contain such a quantity, whether explicitly or implicitly. Here it is promoted to a dynamical field governing gravitational redshift and temporal variation, with its own action and field equations (to be developed in subsequent sections).

The coherence-rate perspective offers an additional methodological route for constructing temporal and spatial intervals without presupposing a metric or a propagation constant as primitive notions. In conventional formulations, time, distance, and the invariant speed of light are introduced operationally through interdependent measurement procedures. This is fully consistent within General Relativity, yet it leaves open the question of whether these quantities may be derived from a more elementary operational substrate. The coherence-rate framework shows that metric structure can, in principle, emerge from comparisons of local dynamical rates alone. Rather than modifying or replacing General Relativity, this approach suggests a complementary foundation that may be useful for further theoretical exploration, especially in regimes where operational definitions become primary (precision clocks, interferometric networks, and quantum-limited measurements).

Comparison with the Metric Formulation

For clarity, Table 1 summarizes the correspondence between quantities in General Relativity and their operational expressions in the coherence–rate framework.

Correspondence between GR quantities and the coherence–rate formulation.

Quantity	General Relativity	Coherence–Rate Framework
Time dilation	\(\sqrt{g_{00}}\)	\(\omega\)
Newtonian potential	\(\Phi\)	\(c^2(\omega-1)\)
Weak–field metric	\(g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}\)	\(g_{\mu\nu}(\omega)\) via (3.1)
Field equation	\(G_{\mu\nu}=8\pi G T_{\mu\nu}\)	\(\Box \omega - m_\omega^2 \omega = \alpha T\)
Propagating modes	tensor (spin–2)	scalar (spin–0) + induced geometry

Effective Geometry Induced by the Coherence Rate

The operational definition of the coherence rate \(\omega(x)\) specifies how proper time increments are related to a universal reference interval but does not, by itself, impose a geometric structure. In this section we show how an effective spacetime metric follows from the requirement that \(\omega(x)\) encode measurable clock-rate differences and that spatial intervals be defined in a manner consistent with local isotropy and the existence of a standard of rest.

Temporal part of the effective metric

By definition, the proper-time increment recorded by a clock located at \(x\) satisfies

\[d\tau(x) = \omega(x)\, dt,\]

where \(t\) is the coordinate time associated with the chosen reference slicing. This relation implies immediately that the temporal component of any compatible metric must satisfy

\[g_{00}(x) = \omega^2(x).\]

This identifies \(\omega(x)\) with the redshift factor in stationary configurations and ensures consistency with the operational comparison of clock frequencies.

Constraints on the spatial sector

The spatial part of the metric is not fixed by clock comparisons alone. We therefore impose two physically motivated requirements:

1. Local spatial isotropy. In the rest frame of the local observer, spatial intervals must be rotationally symmetric. This requires \[g_{ij}(x) = -\,\sigma(x)\,\delta_{ij}\] for some positive scalar function \(\sigma(x)\).

2. Recovery of the Newtonian limit. In weak, static fields where \(\omega(x)=1+\epsilon(x)\) with \(|\epsilon|\ll 1\), geodesic motion of slow test bodies must reduce to \[\frac{d^2 \mathbf{x}}{dt^2} = -\nabla \Phi, \qquad \Phi = c^2 \epsilon + \mathcal{O}(\epsilon^2).\]

These conditions uniquely determine \(\sigma(x)\) up to higher-order disformal terms.

Derivation of the spatial scale factor

Consider the spacetime line element

\[ds^2 = \omega^2(x) c^2 dt^2 - \sigma(x)\, d\mathbf{x}^2.\]

For a slowly moving test particle \((|\mathbf{v}| \ll c)\), we have

\[ds \simeq \omega c\,dt \left[1-\frac{\sigma\,\mathbf{v}^2}{2\omega^2c^2}\right].\]

The Lagrangian for spatial motion becomes

\[L = \frac{\sigma}{2\omega} \mathbf{v}^2 - c^2(\omega -1),\]

and the Euler–Lagrange equations give the acceleration

\[\frac{d^2 x^i}{dt^2} = -\frac{c^2}{2}\,\partial_i \left(\frac{\omega^2}{\sigma}\right) + \mathcal{O}(v^2).\]

To recover Newtonian dynamics with potential \(\Phi = c^2(\omega-1)\), we require

\[\frac{\omega^2}{\sigma} = \omega^2 + \mathcal{O}(\epsilon^2),\]

which implies

\[\sigma(x) = \omega^{-2}(x).\]

Thus the spatial metric is uniquely determined as

\[g_{ij}(x) = -\omega^{-2}(x)\,\delta_{ij}.\]

Resulting effective metric

The full effective metric compatible with the operational definition of \(\omega\) is

\[ds^2 = \omega^2(x)c^2 dt^2 - \omega^{-2}(x)\, d\mathbf{x}^2. \tag{3.1} \label{eq:metric}\]

Expanding for \(\omega=1+\epsilon\) yields

\[g_{00} = 1 + 2\epsilon + \mathcal{O}(\epsilon^2), \qquad g_{ij} = -\left(1 - 2\epsilon + \mathcal{O}(\epsilon^2)\right)\delta_{ij},\]

and identifying \(\epsilon = \Phi/c^2\) reproduces the standard weak-field limit of General Relativity. Therefore the metric [eq:metric] is fully consistent with observed gravitational redshift, time dilation, and Newtonian dynamics.

Example: Static Spherically Symmetric Configuration

To illustrate the structure of the induced geometry, consider a weak, static, spherically symmetric configuration. As shown in Sec. 5, the coherence rate outside a point mass satisfies \[\omega(r) = 1 - \frac{GM}{c^2 r} + O\!\left(\frac{1}{r^2}\right).\]

Substituting this into the effective metric (3.1) yields \[ds^2 = \left(1 - \frac{2GM}{c^2 r}\right)c^2 dt^2 - \left(1 + \frac{2GM}{c^2 r}\right) d\mathbf{x}^2 + O\!\left(\frac{1}{r^2}\right),\] which matches the isotropic weak–field Schwarzschild metric at first post–Newtonian order. Thus the operational reconstruction from \(\omega(x)\) reproduces the standard static gravitational field in the appropriate limit.

Comments on disformal extensions

More general metrics of the disformal form

\[g_{\mu\nu}(\omega,\partial\omega) = \omega^2 \eta_{\mu\nu} + \Lambda(\omega)\frac{\partial_\mu\omega\,\partial_\nu\omega}{m_\omega^2}\]

are allowed without violating Lorentz symmetry in the high-coherence limit \(\omega\to 1\), but they contribute only at second order in gradients and are therefore negligible for laboratory and Solar-System scales. The minimal form Eq. [eq:metric] is adopted throughout the remainder of the paper.

Coherence Field Action and Field Equations

Having established the effective metric induced by the operationally defined coherence rate \(\omega(x)\), we now construct a covariant field theory for \(\omega\) itself. The aim is to obtain a dynamical equation that

(i) is consistent with diffeomorphism invariance, (ii) reproduces the Newtonian limit derived in Sec. 5, (iii) reduces to General Relativity in weak fields, and (iv) yields observable and testable deviations controlled by physical parameters.

Principles guiding the construction

The action for \(\omega(x)\) must satisfy the following physical and mathematical requirements:

1. Covariance. The theory must be formulated on an arbitrary spacetime \((M,g_{\mu\nu})\).

2. Locality. Only \(\omega\) and its first derivatives may appear, except in potential terms.

3. Stability. The kinetic term must have the correct sign to avoid ghost propagation.

4. Minimality. No higher-derivative operators or non-local terms are introduced at the fundamental level.

5. Matter coupling. The interaction between \(\omega\) and matter must reduce to the Newtonian Poisson equation in the appropriate limit.

Metric Variations and the Status of the Effective Metric

In the coherence–rate formulation the metric is not introduced as an independent dynamical field. Instead, \(g_{\mu\nu}\) is an effective geometric structure induced by the operational definition of \(\omega(x)\) described in Sec. 3. The action (4.1) is therefore varied only with respect to \(\omega\), while the metric must be treated as a functional \(g_{\mu\nu}[\omega]\).

Variations of the total action thus take the form \[\delta S = \int d^4 x \sqrt{-g} \left[ \frac{\delta \mathcal{L}}{\delta \omega} + \frac{\partial \mathcal{L}}{\partial g_{\mu\nu}} \frac{\partial g_{\mu\nu}}{\partial \omega} \right]\delta\omega,\] which makes explicit that the induced variation of the metric is included through the dependence \(g_{\mu\nu} = g_{\mu\nu}(\omega)\).

This reflects the operational foundation of the framework: geometric quantities are reconstructed from measurable clock–rate relations rather than postulated a priori. The effective metric retains full covariance, and diffeomorphism invariance is preserved because \(g_{\mu\nu}[\omega]\) transforms as a standard tensor field under coordinate changes.

Minimal covariant action

The simplest action satisfying the above criteria is

\[S = \int d^4x \sqrt{-g}\left[ \frac{M_\omega^2}{2}\, g^{\mu\nu}(\nabla_\mu \omega)(\nabla_\nu \omega) - V(\omega) - \alpha\,\omega\, T \right], \tag{4.1} \label{eq:action}\]

where:

• \(M_\omega\) is the coherence-field normalization scale (with dimensions of mass),

• \(V(\omega)\) is a self-interaction potential,

• \(T = T^{\mu}{}_{\mu}\) is the trace of the matter energy–momentum tensor,

• \(\alpha\) is a dimensionless coupling constant to be fixed by Newtonian consistency.

The coupling \(\omega T\) is the lowest-dimension covariant scalar interaction between \(\omega\) and matter. It is analogous to Jordan-frame scalar–tensor theories, but with the crucial difference that \(\omega\) has a direct operational meaning as a physical clock-rate field.

Field equation

Varying [eq:action] with respect to \(\omega\) yields the Euler–Lagrange equation

\[M_\omega^2 \Box \omega + \frac{dV}{d\omega} = \alpha\, T, \tag{4.2} \label{eq:eom}\]

where \(\Box = g^{\mu\nu}\nabla_\mu\nabla_\nu\) is the covariant d’Alembertian.

In the weak-field regime, expanding \(\omega = 1 + \delta\omega\) and \(V(\omega) = \frac{1}{2}M_\omega^2 m_\omega^2 (\delta\omega)^2 + \cdots\) leads to

\[\Box \delta\omega - m_\omega^2 \delta\omega = \frac{\alpha}{M_\omega^2}\,T, \tag{4.3} \label{eq:KG}\]

a Klein–Gordon equation with source.

Stress-energy tensor

Varying the action with respect to \(g_{\mu\nu}\) gives the stress-energy tensor of the coherence field:

\[T^{\omega}_{\mu\nu} = M_\omega^2\left( \nabla_\mu\omega\nabla_\nu\omega - \frac{1}{2}g_{\mu\nu} (\nabla\omega)^2 \right) + g_{\mu\nu} V(\omega) - \alpha g_{\mu\nu}\,\omega T. \tag{4.4} \label{eq:stress-tensor}\]

Diffeomorphism invariance of the total action requires

\[\nabla_\mu \left(T^{\mu\nu}_{\rm matter} + T^{\mu\nu}_\omega \right) = 0. \tag{4.5}\]

In regions where \(\omega\) varies slowly, \(T^{\mu\nu}_\omega\) is small and ordinary conservation of matter stress-energy is recovered to high accuracy.

Physical interpretation

The coherence field modifies the gravitational sector through two mechanisms:

1. It alters the effective metric via Eq. [eq:metric] of Sec. 3.

2. It contributes additional stress-energy through Eq. [eq:stress-tensor].

As a result, the theory behaves as a scalar–tensor modification of general relativity where the scalar field has a direct operational interpretation as the local clock-rate modulation. The potential \(V(\omega)\) determines large-scale dynamics and the mass \(m_\omega\) controls the range of deviations from GR.

Vacuum and stationary solutions

A natural vacuum is defined by

\[\omega = 1, \qquad V'(1) = 0,\]

corresponding to a uniform redshift-free state. Perturbations around this vacuum obey the linearized Klein–Gordon equation [eq:KG].

Static, spherically symmetric solutions satisfy

\[\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{d\omega}{dr}\right) - m_\omega^2 \omega = \frac{\alpha}{M_\omega^2}\,\rho(r)c^2.\]

In the massless limit \(m_\omega \to 0\) and for a point source, this yields

\[\omega(r) = 1 - \frac{GM}{c^2 r},\]

which reproduces the Newtonian potential via \(\Phi_{\rm eff} = c^2 (\omega - 1)\).

Newtonian Limit and Coupling Determination

To establish consistency with Newtonian gravity, we analyze the static, weak-field regime of the coherence field theory. In this limit the metric and the field \(\omega\) deviate only slightly from their vacuum values, and matter sources are characterized by non-relativistic mass density \(\rho\).

Weak-field expansion

Let \[\omega(x) = 1 + \delta\omega(x), \qquad |\delta\omega| \ll 1,\] and expand the metric [eq:metric] to first order: \[g_{00} = 1 + 2 \delta\omega + \mathcal{O}(\delta\omega^2), \qquad g_{ij} = -\left(1 - 2\delta\omega\right)\delta_{ij} + \mathcal{O}(\delta\omega^2).\]

In the Newtonian approximation, spatial velocities of test particles satisfy \(|\mathbf{v}| \ll c\). The geodesic equation reduces to \[\frac{d^2 x^i}{dt^2} = -\partial_i \Phi_{\rm eff} + \mathcal{O}(v^2/c^2),\] where \(\Phi_{\rm eff}\) is the effective Newtonian potential.

Comparing with the non-relativistic limit of the metric, we identify \[\Phi_{\rm eff}(x) = c^2\delta\omega(x) + \mathcal{O}\bigl((\delta\omega)^2\bigr). \tag{5.1} \label{eq:PhiEff}\]

Thus the Newtonian potential is directly encoded in the first-order deviation of the coherence rate \(\omega\) from its vacuum value.

The mass parameter \(m_\omega\) should be interpreted in the same manner as the scalar mass in conventional scalar–tensor and Horndeski-type theories. It is not postulated as a fixed constant but treated as an empirically constrained parameter entering the effective field theory. Its role as a regulator of coherence-field range and Yukawa-type corrections is standard in EFT approaches to gravity.

Static field equation

In the static, non-relativistic limit with negligible pressure (\(T \simeq \rho c^2\)) the field equation [eq:KG] becomes \[M_\omega^2 \nabla^2 \delta\omega - M_\omega^2 m_\omega^2\,\delta\omega = -\,\alpha\,\rho c^2. \tag{5.2} \label{eq:staticKG}\]

The term proportional to \(m_\omega^2\) produces Yukawa-type corrections. For distances much shorter than \(m_\omega^{-1}\), the mass term can be neglected: \[M_\omega^2 \nabla^2 \delta\omega \simeq -\,\alpha\,\rho c^2. \tag{5.3} \label{eq:5.3}\]

Matching to the Poisson equation

Newtonian gravity requires that the potential \(\Phi\) satisfy \[\nabla^2 \Phi = 4\pi G\rho. \tag{5.4} \label{eq:Poisson}\]

Using Eq. [eq:PhiEff], we have \[\nabla^2 \Phi_{\rm eff} = c^2 \nabla^2 \delta\omega. \tag{5.5} \label{eq:5.5}\]

Substituting Eq. [eq:staticKG] (with \(m_\omega L \ll 1\)) yields \[c^2 \nabla^2 \delta\omega = -\,\frac{\alpha c^2}{M_\omega^2}\,\rho c^2. \tag{5.6} \label{eq:5.6}\]

Matching this with the Poisson equation [eq:Poisson] requires \[4\pi G \rho = -\,\frac{\alpha c^4}{M_\omega^2}\,\rho.\]

This fixes the coupling uniquely as \[\boxed{ \alpha = \frac{4\pi G M_\omega^2}{c^2} }. \tag{5.7} \label{eq:coupling}\]

No fine-tuning is involved: the coupling is determined entirely by the requirement that test-particle motion reproduce the Newtonian limit.

Resulting Newtonian potential

With [eq:coupling], the solution of Eq. [eq:staticKG] for a point mass \(M\) and \(m_\omega \to 0\) is \[\delta\omega(r) = -\frac{GM}{c^2 r}. \tag{5.8} \label{eq:domegaPoint}\]

By Eq. [eq:PhiEff], the Newtonian potential becomes \[\Phi_{\rm eff}(r) = -\,\frac{GM}{r},\] as required.

For finite \(m_\omega\), the solution is Yukawa-suppressed: \[\delta\omega(r) = -\frac{GM}{c^2 r} e^{-m_\omega r}, \tag{5.9} \label{eq:5.9}\] leading to short-range deviations from Newtonian gravity.

Consistency with General Relativity

The identification [eq:PhiEff] implies that the coherence field reproduces \[g_{00}(r) = 1 + \frac{2\Phi(r)}{c^2}\] to leading order. Together with the spatial metric in Sec. 3, this ensures that weak-field gravitational redshift, time dilation, and test-body motion agree with General Relativity at first post-Newtonian order.

Therefore the Newtonian limit provides a strong consistency condition that determines the coupling \(\alpha\) and confirms that the theory is compatible with all classical, weak-field gravitational phenomena.

Parametrized Post-Newtonian Analysis

To test the compatibility of the coherence-field framework with Solar-System observations, we derive the parametrized post-Newtonian (PPN) parameters in the static, weak-field limit. The effective metric obtained in Sec. 3 provides the starting point, and the small deviations \(\delta\omega\) derived in Sec. 5 determine the post-Newtonian potentials.

Metric expansion

Let \[\omega = 1 + \epsilon, \qquad |\epsilon| \ll 1,\] with \[\epsilon = \frac{\Phi_{\rm eff}}{c^2}.\]

From Eq. (3.1), the metric components are

\[g_{00} = \omega^2 = 1 + 2\epsilon + 2\epsilon^2 + \mathcal{O}(\epsilon^3), \tag{6.1} \label{eq:6.1}\]

\[g_{ij} = -\omega^{-2}\delta_{ij} = -\left(1 - 2\epsilon + 3\epsilon^2 + \mathcal{O}(\epsilon^3)\right)\delta_{ij}. \tag{6.2} \label{eq:6.2}\]

These expressions must be matched to the standard PPN metric in isotropic gauge:

\[g_{00} = 1 - \frac{2U}{c^2} + \frac{2\beta U^2}{c^4} + \cdots, \tag{6.3} \label{eq:6.3}\]

\[g_{ij} = -\left(1 + \frac{2\gamma U}{c^2} + \cdots\right)\delta_{ij}, \tag{6.4} \label{eq:6.4}\]

where \(U\) is the Newtonian potential.

Identification of the Newtonian potential

From Sec. 5, we have the identification

\[\Phi_{\rm eff} = -\, U, \qquad \epsilon = -\frac{U}{c^2}. \tag{6.5} \label{eq:6.5}\]

Substituting [eq:6.5] into [eq:6.1]–[eq:6.2] gives:

\[g_{00} = 1 - \frac{2U}{c^2} + \frac{2U^2}{c^4} + \mathcal{O}\!\left(\frac{U^3}{c^6}\right), \tag{6.6} \label{eq:6.6}\]

\[g_{ij} = -\left( 1 + \frac{2U}{c^2} + \mathcal{O}\!\left(\frac{U^2}{c^4}\right) \right)\delta_{ij}. \tag{6.7} \label{eq:6.7}\]

This matches the PPN form.

Determination of \(\gamma\)

Comparing [eq:6.7] with [eq:6.4], we obtain

\[\gamma = 1 + \mathcal{O}\!\left((\nabla\omega)^2\right). \tag{6.8} \label{eq:6.8}\]

All first-order PPN constraints are satisfied. Deviations arise only from disformal and higher-gradient terms.

Physical meaning of \(\gamma\).

In the standard PPN framework, the parameter \(\gamma\) quantifies the amount of spatial curvature produced per unit Newtonian potential. Its observational role is most clearly seen in the deflection of light and in the Shapiro time delay. Although we do not derive these expressions explicitly here, the identification of \(g_{ij}\) in Eqs. (6.2) and (6.7) is sufficient to guarantee that the leading-order light–propagation effects agree with General Relativity whenever \(\gamma = 1\).

Determination of \(\beta\)

Comparing [eq:6.6] with [eq:6.3], we obtain

\[\beta = 1 + \mathcal{O}\!\left( m_\omega^{-2}\nabla^2 \omega,\; (\nabla\omega)^2 \right). \tag{6.9} \label{eq:6.9}\]

Thus \(\beta\) also agrees with General Relativity at leading order.

Consistency with Solar-System tests

Observations from the Cassini tracking experiment impose

\[|\gamma - 1| \lesssim 2.3 \times 10^{-5}. \tag{6.10} \label{eq:6.10}\]

Binary pulsar tests constrain

\[|\beta - 1| \lesssim 10^{-4}. \tag{6.11} \label{eq:6.11}\]

Using [eq:6.8]-[eq:6.9], these translate into bounds on the coherence-field gradients:

\[|\nabla\omega| \ll 10^{-15}, \qquad m_\omega^{-1} \ll 10^{11}\,\mathrm{m}. \tag{6.12} \label{eq:6.12}\]

Therefore all Solar-System and binary-pulsar observations are satisfied provided that coherence gradients are sufficiently small on astronomical scales and that the coherence mass \(m_\omega\) is not ultralight.

Summary

The coherence field reproduces General Relativity at first post-Newtonian order. Deviations appear only at second order and are suppressed either by coherence gradients or by the range \(m_\omega^{-1}\) of the scalar interaction. This makes the framework compatible with current observational tests while allowing distinct signatures in laboratory, interferometric, and high-precision timing experiments.

Linearized Spectrum and Propagating Modes

To understand the dynamical degrees of freedom and stability properties of the coherence field \(\omega(x)\), we analyze perturbations around the vacuum solution. This provides the dispersion relation, group velocity, and causal structure of the theory, and determines the conditions under which the coherence field represents a healthy propagating mode.

Linearization around the vacuum

Let \[\omega(x) = 1 + \delta\omega(x), \qquad |\delta\omega|\ll 1,\] and expand the action [eq:action] to quadratic order. Assuming \(V'(1)=0\) and \(V''(1) = M_\omega^2 m_\omega^2\), the quadratic action becomes

\[S^{(2)} = \frac{M_\omega^2}{2} \int d^4x \sqrt{-g} \left[ (\nabla \delta\omega)^2 - m_\omega^2 (\delta\omega)^2 \right]. \tag{7.1} \label{eq:quadAction}\]

Thus \(\delta\omega\) behaves as a canonical scalar field with mass \(m_\omega\).

Linearized equation of motion

In Minkowski spacetime and in the absence of matter sources, the field equation reduces to

\[(\Box - m_\omega^2)\,\delta\omega = 0. \tag{7.2} \label{eq:KGfree}\]

This is the standard Klein–Gordon equation. Therefore the theory contains a single, healthy scalar degree of freedom.

Plane-wave solutions

Solutions of Eq. [eq:KGfree] take the form

\[\delta\omega(t,\mathbf{x}) = A\, e^{-i(\Omega t - \mathbf{k}\cdot \mathbf{x})},\]

giving the dispersion relation

\[\Omega^2 = c^2 k^2 + c^2 m_\omega^2. \tag{7.3} \label{eq:dispersion}\]

This is precisely the relativistic dispersion for a massive scalar. No additional modes or instabilities arise.

Group velocity and causality

The group velocity is

\[v_g = \frac{\partial\Omega}{\partial k} = c\,\frac{k}{\sqrt{k^2 + m_\omega^2}}. \tag{7.4} \label{eq:vg}\]

Thus

\[0 \le v_g < c,\]

and the coherence mode is strictly subluminal, with causal propagation.

In the massless limit \(m_\omega \to 0\), \(v_g\to c\) and the theory reduces smoothly to a relativistic massless scalar.

Stability

The quadratic Hamiltonian associated with [eq:quadAction] is

\[\mathcal{H} = \frac{M_\omega^2}{2} \left[ (\partial_t \delta\omega)^2 + c^2 (\nabla\delta\omega)^2 + m_\omega^2 c^2 (\delta\omega)^2 \right], \tag{7.5} \label{eq:H}\]

which is manifestly positive definite provided

\[M_\omega^2 > 0, \qquad m_\omega^2 \ge 0. \tag{7.6}\]

Therefore the theory has:

• no ghost instabilities (\(M_\omega^2>0\)),

• no tachyonic instabilities (\(m_\omega^2 \ge 0\)),

• no gradient instabilities (correct sign of spatial term).

Role of disformal corrections

If disformal terms of the form

\[\frac{\partial_\mu\omega\,\partial_\nu\omega}{m_\omega^2}\]

are included in the metric (Sec. 3), they contribute only at order \(\mathcal{O}((\nabla\omega)^2)\) to the quadratic action. In the linearized regime \(\delta\omega \ll 1\), they therefore do not modify:

• the dispersion relation,

• the group velocity,

• the stability criteria.

This guarantees that the minimal spectrum [eq:dispersion] is robust.

Summary of dynamical properties

The coherence field propagates as a healthy scalar mode with:

\[\Omega_k^2 = c^2 (k^2 + m_\omega^2),\]

\[0 \le v_g < c,\]

and no ghosts, tachyons, or anomalies. The mass scale \(m_\omega\) controls the range of deviations from General Relativity. The linearized regime is fully stable, and the dynamics are equivalent to those of a standard massive scalar coupled to matter via the trace \(T\).

Experimental Signatures and Constraints

The coherence-rate field \(\omega(x)\) affects all physical processes whose dynamics depend on the accumulation of proper time. Modern precision experiments can therefore probe small deviations from the General-Relativistic redshift and time-dilation relations. In this section we describe the main experimental channels and derive the corresponding constraints on the parameters of the theory.

Clock comparisons and redshift tests

For two stationary clocks at positions \(x_1\) and \(x_2\), the coherence framework predicts

\[\frac{\nu(x_2)}{\nu(x_1)} = \frac{\omega(x_2)}{\omega(x_1)} = 1 + \frac{\Phi_{\rm eff}(x_2)-\Phi_{\rm eff}(x_1)}{c^2} + \mathcal{O}\!\left((\nabla\omega)^2\right). \tag{8.1}\]

Let \(\Delta h\) be the height difference and \(g\) the local gravitational acceleration. Then

\[\Delta \omega \equiv \omega(x_2)-\omega(x_1) \simeq \frac{g\,\Delta h}{c^2} + \delta\omega_{\rm anom}, \tag{8.2}\]

where \(\delta\omega_{\rm anom}\) encodes deviations due to gradients of \(\omega\) beyond GR.

State-of-the-art optical lattice clocks achieve fractional stability

\[\sigma_y \lesssim 10^{-18} \quad \text{per}\ 10^4\ {\rm s}. \tag{8.3}\]

Thus the bound on anomalous coherence gradients is

\[|\nabla\omega|_{\rm anom} \lesssim 10^{-16}\ {\rm m}^{-1}. \tag{8.4} \label{eq:8.4}\]

This constraint matches and strengthens the PPN bound from Sec. 6.

Atom interferometry

Two-space–time-path phase accumulation for a Mach–Zehnder-type atom interferometer under the coherence framework is

\[\Delta\phi = k_{\rm eff} g T^2 + \int_{\Gamma_1}\!\omega\,dt - \int_{\Gamma_2}\!\omega\,dt. \tag{8.5}\]

The first term is the standard GR acceleration phase; the remaining terms are path-dependent coherence contributions.

For interrogation times \(T\sim 0.1\ {\rm s}\), typical phases are

\[\Delta\phi_{\rm GR} \sim 10^6\ {\rm rad}.\]

State-of-the-art sensitivity reaches

\[\delta\phi \sim 10^{-3}\ {\rm rad}, \tag{8.6}\]

yielding the constraint

\[|\nabla\omega|_{\rm anom} \lesssim 10^{-16}\ {\rm m}^{-1} \tag{8.7}\]

for terrestrial-scale baselines.

Moving clocks: aircraft, satellites, and LEO tests

For clocks transported with velocity \(v\), the predicted rate is

\[\omega(v,x) = \omega(x)\left(1 - \frac{v^2}{2c^2} + \delta_{\rm NL}(v,x)\right), \tag{8.8}\]

where \(\delta_{\rm NL}\) is a potential non-linear velocity–coherence correction.

Existing comparison campaigns (airborne and LEO) achieve

\[\left|\frac{\Delta\nu}{\nu}\right| \lesssim 10^{-17}, \tag{8.9}\]

implying

\[|\delta_{\rm NL}(v,x)| \lesssim 10^{-17}. \tag{8.10}\]

This restricts second-order velocity couplings of \(\omega\) to be negligible at orbital speeds \(v\simeq 7\,{\rm km/s}\).

Gravitational-wave phase tracking

If \(\omega\) varies along the photon path of an interferometric detector, the accumulated phase shift is

\[\Delta\phi_\omega = \int_0^L \delta\omega(x)\,k\,dx. \tag{8.11}\]

LIGO/Virgo/KAGRA phase residuals satisfy

\[|\Delta\phi_{\rm res}| \lesssim 10^{-2}\ {\rm rad} \quad \text{in the}\ 100\!-\!1000\ {\rm Hz}\ \text{band}. \tag{8.12}\]

This constrains path-integrated deviations:

\[\int_0^L |\nabla\omega|\,dx \lesssim 10^{-2}. \tag{8.13}\]

For \(L\simeq 4\,{\rm km}\), this gives

\[|\nabla\omega| \lesssim 10^{-6}\ {\rm m}^{-1} \quad (\text{interferometer-scale}). \tag{8.14}\]

Although weaker than clock tests, GW data probe different regimes.

Pulsar timing arrays

For millisecond pulsars,

\[r_\omega = \int \delta\omega(x(t))\,dt \tag{8.15}\]

contributes to timing residuals. PTA timing precision is

\[\sigma_r \lesssim 50\ {\rm ns}. \tag{8.16}\]

Propagation through a region of size \(L\) gives the constraint

\[|\delta\omega| \lesssim 10^{-14} \quad (\text{for Galactic baselines}). \tag{8.17}\]

This constrains ultra-long-range modes (small \(m_\omega\)).

Quantum coherence and potential dependence

For a quantum two-level system with coherence time \(T_2\), the predicted dependence is

\[T_2^{-1}(\Phi) = T_{2,0}^{-1} + \alpha_\omega \frac{\Delta\Phi}{c^2}, \tag{8.18} \label{eq:8.18}\]

where \(\alpha_\omega\) is a material-dependent coefficient induced by \(\omega\)–matter coupling.

Experiments with trapped ions and neutral atoms can reach

\[\delta(T_2^{-1}) \sim 10^{-3}, \tag{8.19} \label{eq:8.19}\]

giving constraints consistent with [eq:8.4].

Combined constraints

Table 1 summarizes the bounds from all channels.

Implications for coherence-field parameters

Using the results of Sec. 5 and Sec. 6, the combined constraints imply:

\[|\nabla\omega| \lesssim 10^{-16}\ {\rm m}^{-1},\]

\[m_\omega^{-1} \ll 10^{11}\ {\rm m} \quad\text{(to satisfy PPN limits)},\]

\[\alpha = \frac{4\pi G M_\omega^2}{c^2} \quad\text{(from Newtonian matching)}.\]

These results show that the coherence field is compatible with all existing weak-field tests and that its deviations are observable in next-generation precision experiments.

While several effects appear at the \(10^{-17}\)–\(10^{-18}\) level, such sensitivities are already being approached by state-of-the-art optical clocks, clock networks, and atom-interferometric platforms. The framework is therefore directly testable within the experimental roadmap of the coming decade, rather than relying on speculative future technology.

Relation to Existing Approaches and Predictive Outlook

The coherence-field framework developed in this work is consistent with General Relativity in the weak-field regime but differs in its underlying operational interpretation. In this section we outline its relation to several established approaches in gravitational physics and quantum foundations, emphasizing structural similarities, points of departure, and potential observational discriminants.

Relation to scalar–tensor and metric theories

At the level of field equations, the coherence rate \(\omega(x)\) plays the role of a scalar field coupled to matter through the trace \(T\) and inducing a conformal modification of the metric. This places the theory within the broader class of scalar–tensor frameworks. However, two distinctions are noteworthy:

1. The field \(\omega\) has a direct operational definition in terms of clock-rate ratios, rather than being introduced as an independent gravitational degree of freedom.

2. The effective metric arises from the requirement of consistency with this operational definition, rather than being posited a priori.

These differences do not alter phenomenology at the weak-field level but provide a distinct interpretational underpinning.

Relation to teleparallel formulations

Teleparallel gravity attributes gravitational effects to torsion rather than curvature. In such approaches, the metric retains its standard form while non-metric variables encode dynamical information. The coherence-field theory differs in that

• \(\omega(x)\) modifies the metric directly through conformal factors,

• torsion is not employed, and

• all deviations appear through scalar modulation rather than vector or tensor structures.

Nevertheless, both frameworks share the idea that gravitational effects can be reinterpreted in terms of modifications to locally measurable dynamical rates.

Relation to causal-set and discrete frameworks

Causal-set theory interprets spacetime as a partially ordered set of elementary events. In such a view, temporal intervals emerge from counting relations. The coherence-field framework is compatible with this perspective in the sense that:

• temporal structure is derived from a fundamental rate,

• \(\omega(x)\) encodes local deviations from a reference rate, and

• redshift relations emerge from comparisons of these rates.

Unlike causal-set models, the coherence field produces a smooth effective metric in continuum regimes, without requiring discreteness at observable scales.

Relation to thermodynamic and entropic gravity

Thermodynamic derivations of Einstein’s equations (e.g., Jacobson’s Einstein-equation-of-state approach) suggest that gravitational dynamics emerge from equilibrium conditions on local causal horizons. The present theory resonates with this perspective through the idea that temporal rates—and hence clock phases—encode energetic gradients. However:

• we employ a dynamical scalar field rather than thermodynamic relations,

• the effective metric is derived from \(\omega\), not from entropy variations,

• the theory remains local and Lagrangian-based.

Thus the frameworks are conceptually related but technically distinct.

Relation to information-theoretic bounds

Information-theoretic limits such as the Margolus–Levitin bound and Bekenstein entropy bounds associate energy with rates of physical state evolution. The coherence-field interpretation reflects this idea inasmuch as energy corresponds to a deviation of local dynamical rates from a reference maximum. However the coherence theory does not assume or require any fundamental limits from quantum information; the connection is purely structural.

Distinctive signatures and predictive differences

Although the theory reproduces General Relativity at first post-Newtonian order, its scalar degree of freedom leads to the following potentially observable differences:

1. Yukawa-type corrections to the Newtonian potential for finite \(m_\omega\), detectable at sub-millimeter or geophysical scales.

2. Nonlinear coherence gradients, which modify redshift measurements at the level of \(10^{-17}\)–\(10^{-18}\) in next-generation optical clock arrays.

3. Path-dependent phase accumulation in atom interferometry beyond the standard \(k_{\rm eff} g T^2\) term.

4. Frequency-dependent modifications to gravitational-wave phase evolution through effective variations in \(\omega\) along photon or graviton trajectories.

5. Stochastic coherence fluctuations, which may contribute to timing-array residuals or noise floors in high-stability oscillators.

These signatures are complementary to Solar-System PPN tests and provide avenues for falsification or confirmation of the framework.

Strong–Field Considerations

Although the present analysis focuses on weak–field and post–Newtonian regimes, it is instructive to outline the expected behavior of the coherence field in strong–gravity environments. The effective metric (3.1) may exhibit redshift horizons where \(\omega \rightarrow 0\), corresponding to vanishing local coherence rate. Such horizons are geometric rather than singular: curvature scalars remain finite provided \(\omega\) and its first derivatives remain finite.

Because the field equations for \(\omega\) are second–order and well–posed, strong–field configurations are expected to behave similarly to scalar–tensor theories: regularity is maintained unless \(\nabla \omega\) diverges. A complete treatment of compact–object and horizon–scale solutions lies beyond the scope of this work, but the absence of higher derivatives in the action ensures that no strong–field pathologies analogous to gradient instabilities or ghost modes arise.

Operational origin of time dilation

In this formulation, gravitational time dilation is not a geometric primitive but an operational statement about the rate at which local physical systems evolve. A clock at position \(x\) accumulates phase at rate \(\nu_{\rm loc}(x)\), and the coherence field encodes the ratio of this rate to the universal reference scale, \(\omega = \nu_{\rm loc}/\nu_{\rm ref}\). A lower value of \(\omega\) means that all local dynamical processes evolve more slowly relative to the reference rate.

Thus “time dilation’’ is a statement about the suppression of local dynamical rates, with the metric (3.1) emerging as a secondary encoding of these observable relations.

Light Deflection

Because the effective metric (3.1) is conformally related to the isotropic weak–field Schwarzschild metric, null geodesics are invariant under the corresponding conformal rescaling. Therefore the bending angle of light in the coherence–rate framework is identical to the General–Relativistic value at first post–Newtonian order:

\[\hat{\alpha} = \frac{4GM}{c^2 b} + \mathcal{O}\!\left(\frac{GM}{c^2 b}\right)^2 .\]

Thus gravitational lensing in the weak–field regime is reproduced exactly, with possible deviations appearing only through higher–order gradient corrections of \(\omega\).

Quantum Considerations

In this work the coherence field is treated as a classical scalar degree of freedom. A full quantum treatment would require specifying the ultraviolet behavior of the potential \(V(\omega)\) and the normalization scale \(M_\omega\).

Small quantum fluctuations \(\delta\omega\) around the vacuum \(\omega=1\) are suppressed by \(M_\omega\) and propagate according to the Klein–Gordon spectrum derived in Sec. 7. No conflict with causal structure arises because the propagation speed remains subluminal.

Possible quantum signatures include stochastic fluctuations of local clock rates and small corrections to interferometric phases. A detailed quantization of \(\omega\) lies outside the present scope but constitutes a natural extension of the framework.

Outlook

Future developments may include the construction of:

• exact spherically symmetric and cosmological solutions,

• quantum extensions of the coherence field and its fluctuations,

• coupling of \(\omega\) to additional matter sectors,

• potential relations to discrete or information-theoretic models of spacetime.

The operational definition of \(\omega\) ensures that all such developments remain tied to measurable quantities, preserving the empirical accessibility of the framework.

Future Work Roadmap

Several directions for further analysis naturally emerge from the coherence–rate perspective:

• Detailed strong–field solutions for compact objects and possible horizon structures where \(\omega \rightarrow 0\).

• Precision redshift tests using next–generation optical–clock networks capable of resolving gradients of order \(10^{-18}\).

• Coherence–induced corrections to atom–interferometric phase shifts beyond the \(k_{\mathrm{eff}} g T^2\) term.

• Gravitational–wave phase evolution in backgrounds with slowly varying \(\omega(x)\).

• Constraints on stochastic fluctuations of \(\omega\) from pulsar timing arrays and quantum–coherence platforms.

These investigations may help determine whether the operational definition of time provided by the coherence field offers useful complementary insights for gravitational phenomenology.

Limitations of the Framework

The coherence–rate formulation is designed as an operationally motivated extension of weak–field gravitational phenomenology. Several limitations should be noted.

First, the effective metric \(g_{\mu\nu}[\omega]\) does not include independent tensorial degrees of freedom; thus the theory does not aim to reproduce the full nonlinear dynamics of General Relativity outside the weak–field regime.

Second, quantum fluctuations of \(\omega\) are not addressed in this work. The scalar field is treated classically, and a coherent quantum treatment would require specifying the ultraviolet completion of the potential \(V(\omega)\).

Third, the framework is currently formulated for low–curvature and post–Newtonian settings. Applications to black holes, neutron stars, and cosmological backgrounds require additional analysis of the full nonlinear solutions for \(\omega(x)\).

These limitations do not affect the weak–field predictions developed in the main text but delineate the scope of applicability of the coherence–rate approach.

Black–body spectrum

Because both energy levels and emitted frequencies scale with the local coherence rate \(\omega(x)\) in the same way, the Planck black–body spectrum remains locally unchanged. Gravitational redshift of thermal radiation is therefore treated exactly as in General Relativity: the spectrum is preserved and only its frequency scale is shifted by the ratio \(\omega(x_2)/\omega(x_1)\) between emitter and observer.

Conclusion

We have developed a coherence-based formulation of gravitational phenomena in which the local accumulation of proper time is described by a scalar rate field \(\omega(x)\). The field is defined operationally through measurable clock-frequency ratios, providing a direct link between dynamical rates and spacetime structure. From this definition, an effective metric consistent with gravitational redshift and time-dilation experiments is uniquely derived.

A covariant action for the coherence field was constructed, yielding a Klein–Gordon type equation with matter coupling fixed unambiguously by the Newtonian limit. The resulting framework behaves as a scalar–tensor theory that reproduces General Relativity at first post-Newtonian order. Deviations are controlled by the spatial gradients of \(\omega\) and the coherence mass \(m_\omega\), with Yukawa-type corrections at distances comparable to \(m_\omega^{-1}\).

We analyzed the linearized spectrum and demonstrated that the coherence field propagates as a healthy scalar mode with subluminal group velocity and no ghost or tachyonic instabilities. The theory is therefore dynamically consistent and stable.

A range of precision experiments—including optical-clock comparisons, atom interferometry, moving clocks, pulsar timing arrays, gravitational-wave phase tracking, and quantum-coherence measurements—places bounds on the magnitude of coherence gradients and on the parameters of the potential. These tests constrain \[|\nabla\omega| \lesssim 10^{-16}\,{\rm m}^{-1}, \qquad m_\omega^{-1} \ll 10^{11}\,\mathrm{m},\] ensuring compatibility with all current weak-field observations.

The framework provides a consistent and testable scalar-rate reformulation of gravitational phenomena grounded in operationally defined temporal structure. Its phenomenology is fully falsifiable: improved clock networks, long-baseline interferometers, and timing arrays will either confirm the predicted coherence-field effects or restrict the parameter space further. Future work may include exact solutions, cosmological extensions, and quantum generalizations of the coherence field.

Overall, the coherence-field approach offers a unified and experimentally accessible link between local dynamical rates and gravitational phenomenology, complementing the geometric interpretation of General Relativity while remaining observationally consistent with it.

Appendix A. Dimensions and Units

Throughout this work we employ SI units with explicit factors of \(c\), \(\hbar\), and \(G\). The coherence rate \(\omega(x)\) is defined as the ratio of two clock frequencies and is therefore dimensionless:

\[[\omega] = 1.\]

The reference frequency is defined as \[\nu_{\rm ref} = \sqrt{\frac{c^5}{\hbar G}}, \qquad [\nu_{\rm ref}] = {\rm s}^{-1}.\]

A.1 Field normalization and coupling

The action for the coherence field is \[S = \int d^4x\,\sqrt{-g}\left[ \frac{M_\omega^2}{2} g^{\mu\nu}(\nabla_\mu\omega)(\nabla_\nu\omega) - V(\omega) - \alpha\,\omega\, T \right].\]

The dimensions of the quantities appearing in the action are:

\[[M_\omega] = {\rm kg}, \qquad [(\nabla\omega)^2] = {\rm m}^{-2}, \qquad [V] = {\rm J\,m^{-3}},\]

and since \(T = T^\mu{}_\mu\) is a trace of the stress-energy tensor,

\[[T] = {\rm J\,m^{-3}}.\]

Because \(\omega\) is dimensionless, the coupling \(\alpha\) must also be dimensionless:

\[[\alpha] = 1.\]

A.2 Potential and mass parameter

The potential is expanded around the vacuum value \(\omega=1\) as

\[V(\omega) = \frac{1}{2} M_\omega^2 m_\omega^2 (\omega-1)^2 + \cdots.\]

The parameter \(m_\omega\) has units of inverse length:

\[[m_\omega] = {\rm m}^{-1}.\]

This identifies \(m_\omega^{-1}\) as the characteristic range of coherence-field modifications. Higher-order coefficients in \(V(\omega)\) scale as powers of \(M_\omega^2 m_\omega^2\).

A.3 Effective Newtonian potential

To leading order, the effective Newtonian potential is defined as

\[\Phi_{\rm eff} = c^2 (\omega - 1), \qquad [\Phi_{\rm eff}] = {\rm m^2\,s^{-2}},\]

consistent with standard gravitational units.

The matching condition derived in Sec. 5 fixes the coupling,

\[\alpha = \frac{4\pi G M_\omega^2}{c^2},\]

which is dimensionless because

\[[G] = {\rm m^3\,kg^{-1}\,s^{-2}},\qquad [M_\omega^2] = {\rm kg^2}.\]

Thus the dimensional factors cancel.

A.4 Summary of key dimensions

\[\begin{aligned} [\omega] &= 1,\\[4pt] [M_\omega] &= {\rm kg},\\[4pt] [m_\omega] &= {\rm m^{-1}},\\[4pt] [V] &= {\rm J\,m^{-3}},\\[4pt] [T] &= {\rm J\,m^{-3}},\\[4pt] [\alpha] &= 1,\\[4pt] [\Phi_{\rm eff}] &= {\rm m^2\,s^{-2}}. \end{aligned}\]

All quantities in the coherence-field action and its derived equations are therefore dimensionally consistent.

Appendix B. Linearized Spectrum and Stability

This appendix provides a detailed derivation of the linearized spectrum, dispersion relation, Hamiltonian structure, and stability properties of the coherence field introduced in Sec. 4. We expand the full action to quadratic order around the vacuum solution and analyze the resulting dynamical modes.

B.1 Expansion of the action

Let the coherence field be decomposed as \[\omega(x) = 1 + \delta\omega(x), \qquad |\delta\omega|\ll 1.\]

Assuming that the potential satisfies \[V(1) = 0, \qquad V'(1) = 0, \qquad V''(1) = M_\omega^2 m_\omega^2,\] the action [eq:action] becomes, to quadratic order, \[S^{(2)} = \frac{M_\omega^2}{2} \int d^4x \sqrt{-g}\, \left[ g^{\mu\nu}\partial_\mu\delta\omega\,\partial_\nu\delta\omega - m_\omega^2 (\delta\omega)^2 \right]. \tag{B.1} \label{eq:quad_action_app}\]

Thus \(\delta\omega\) is a canonical scalar field with mass \(m_\omega\).

B.2 Linearized field equation

Varying [eq:quad_action_app] yields \[(\Box - m_\omega^2)\,\delta\omega = 0. \tag{B.2} \label{eq:KGfree_app}\]

In a locally inertial frame, \(\Box = \eta^{\mu\nu}\partial_\mu\partial_\nu\), and the equation reduces to the standard Klein–Gordon equation.

B.3 Plane-wave modes

Looking for solutions of the form \[\delta\omega(t,\mathbf{x}) = A e^{-i(\Omega t - \mathbf{k}\cdot\mathbf{x})},\] we obtain the dispersion relation \[\Omega^2 = c^2 k^2 + c^2 m_\omega^2. \tag{B.3} \label{eq:dispersion_app}\]

This corresponds to a single relativistic massive scalar mode, with no extra degrees of freedom.

B.4 Group velocity and causal propagation

The group velocity is \[v_g = \frac{\partial\Omega}{\partial k} = c\,\frac{k}{\sqrt{k^2 + m_\omega^2}}. \tag{B.4} \label{eq:group}\]

Thus, \[0 \le v_g < c.\]

Propagation is therefore strictly causal, with luminal speed recovered as \(m_\omega\to 0\).

B.5 Hamiltonian structure

The canonical Hamiltonian density derived from [eq:quad_action_app] is \[\mathcal{H} = \frac{M_\omega^2}{2} \left[ (\partial_t\delta\omega)^2 + c^2 (\nabla\delta\omega)^2 + c^2 m_\omega^2 (\delta\omega)^2 \right]. \tag{B.5} \label{eq:Hamiltonian_app}\]

This is positive definite provided \[M_\omega^2 > 0, \qquad m_\omega^2 \ge 0. \tag{B.6}\]

Hence:

• No ghost instabilities (positive kinetic energy).

• No tachyonic instabilities (real mass parameter).

• No gradient instabilities (correct sign for spatial derivatives).

The coherence field is dynamically stable in all regimes where the quadratic expansion is valid.

B.6 Role of background curvature

In a curved background \((M,g_{\mu\nu})\), the linearized equation becomes \[\left( \Box_g - m_\omega^2 \right)\delta\omega = 0. \tag{B.7}\]

Curvature corrections enter through \(\Box_g\) but do not alter the sign of the kinetic, gradient, or mass terms. Therefore, local stability properties are unchanged.

B.7 Influence of matter sources

In the presence of matter, the full equation is \[(\Box - m_\omega^2)\,\delta\omega = \frac{\alpha}{M_\omega^2} T, \tag{B.8} \label{eq:KG_app_matter}\] as described in Sec. 5. The source term does not modify the propagating spectrum unless \(T\) varies rapidly on coherence-field length scales. In quasi-static or slowly varying matter configurations, the perturbations remain governed by [eq:dispersion_app].

B.8 Summary

The linearized coherence field satisfies a Klein–Gordon equation with mass \(m_\omega\) and has a standard relativistic dispersion relation. The Hamiltonian is positive definite for \(M_\omega^2>0\) and \(m_\omega^2\ge 0\), ensuring dynamical stability. Propagation is subluminal, and no additional propagating degrees of freedom arise. The linearized limit is therefore consistent and robust across all physically relevant regimes.

Appendix C. Cosmological Background Solutions

We analyze homogeneous and isotropic cosmological solutions of the coherence-field equations derived in Sec. 4. The goal is to determine the behavior of \(\omega\) in an FLRW background and identify the conditions under which the coherence field contributes effectively to cosmic expansion.

The purpose of this appendix is not to propose an alternative cosmological model to \(\Lambda\)CDM, but to show that the coherence field admits consistent homogeneous solutions and remains compatible with the standard expansion history. A detailed cosmological model lies beyond the scope of the present work; crucially, for any \(m_\omega>0\) the field dynamically relaxes to \(\omega\to 1\), reproducing the GR limit at late times.

C.1 Background metric and symmetry reduction

Consider a spatially flat FLRW metric, \[ds^2 = c^2 dt^2 - a^2(t)\, d\mathbf{x}^2, \tag{C.1} \label{eq:C.1}\] where \(a(t)\) is the scale factor. Spatial homogeneity and isotropy require that the coherence rate depend only on cosmic time, \[\omega = \omega(t). \tag{C.2} \label{eq:C.2}\]

The matter sector is described by a perfect fluid with energy density \(\rho(t)\) and pressure \(p(t)\), giving stress-energy trace \[T = \rho c^2 - 3p. \tag{C.3}\]

C.2 Reduced action and field equation

Substituting [eq:C.1]–[eq:C.2] into the action [eq:action] and using \(\sqrt{-g} = a^3(t)\) gives the reduced Lagrangian \[L_{\rm eff} = \frac{M_\omega^2 a^3}{2} \left[ \dot{\omega}^2 - m_\omega^2 c^2 (\omega-1)^2 \right] - a^3 \alpha\,\omega T. \tag{C.4}\]

Variation yields the background equation of motion, \[M_\omega^2\left(\ddot{\omega} + 3H\dot{\omega}\right) + M_\omega^2 m_\omega^2 c^2(\omega-1) = \alpha\, T, \tag{C.5} \label{eq:omega_FLRW}\] where \[H \equiv \dot{a}/a. \tag{C.6}\]

C.3 Interpretation

Equation [eq:omega_FLRW] is a damped, driven oscillator with Hubble friction, mass term, and matter-source term. The source-to-friction ratio determines whether \(\omega\) tracks the matter density or evolves independently.

C.4 Static and quasi-static solutions

In the regime where \(\dot{\omega} \approx 0\) and \(\ddot{\omega}\approx 0\), Eq. [eq:omega_FLRW] reduces to \[m_\omega^2 c^2 (\omega-1) = \frac{\alpha}{M_\omega^2} T. \tag{C.7}\]

For pressureless matter (\(p\approx 0\)), \[\omega-1 \simeq \frac{\alpha}{M_\omega^2 m_\omega^2 c^2} \rho c^2 = \frac{4\pi G\rho}{m_\omega^2 c^2}, \tag{C.8} \label{eq:C.8}\] using the coupling \(\alpha = 4\pi G M_\omega^2/c^2\) from Sec. 5.

Thus:

- **If \(m_\omega\) is large**, deviations decay as \(1/m_\omega^2\) and \[\omega \to 1,\] meaning the coherence field *does not participate* in cosmological dynamics.

- **If \(m_\omega\) is small**, \(\omega\) acquires a slow, density-dependent shift.

C.5 Dynamical solutions

Neglecting the source term during radiation domination (\(T \approx 0\)), the equation reduces to \[\ddot{\omega} + 3H\dot{\omega} + m_\omega^2 c^2(\omega-1) = 0. \tag{C.9}\]

Solutions:

1. **Massive coherence field (\(m_\omega \gg H/c\))** \[\omega(t) - 1 \propto a^{-3/2}(t) \cos(m_\omega c t + \phi_0), \tag{C.10}\] exponentially damped by Hubble friction.

2. **Light coherence field (\(m_\omega \ll H/c\))** \[\omega(t) = \omega_0 + \omega_1\, t^{-1} + \cdots, \tag{C.11}\] drifting slowly toward a constant determined by initial conditions.

C.6 Effective energy density and equation of state

The energy density and pressure of \(\omega\) follow from \[\rho_\omega c^2 = \frac{M_\omega^2}{2}\dot{\omega}^2 + \frac{M_\omega^2 m_\omega^2 c^2}{2}(\omega-1)^2 + V_{\rm higher}(\omega), \tag{C.12}\]

\[p_\omega = \frac{M_\omega^2}{2}\dot{\omega}^2 - \frac{M_\omega^2 m_\omega^2 c^2}{2}(\omega-1)^2 - V_{\rm higher}(\omega). \tag{C.13}\]

Thus:

- **For kinetic domination**: \(p_\omega \simeq \rho_\omega c^2\) (stiff fluid). - **For potential domination** (massive field near minimum): \(p_\omega \simeq -\rho_\omega c^2\). - **For slow roll of light fields**: intermediate equation of state.

These results mirror standard scalar-field cosmology.

C.7 Late-time attractor

If \(m_\omega\) is nonzero, cosmic expansion drives \(\omega\) to the minimum of its potential: \[\omega(t) \to 1 \qquad (t \to \infty). \tag{C.14}\]

Thus the coherence field naturally restores the GR limit at late times.

C.8 Summary

Homogeneous cosmological dynamics of \(\omega(t)\) reduce to a damped, driven scalar-field equation. The behavior depends on the mass scale \(m_\omega\):

• Massive coherence field (\(m_\omega \gg H/c\)): rapid decay to \(\omega=1\), negligible cosmological impact.

• Light coherence field (\(m_\omega \ll H/c\)): slow evolution, potentially contributing to effective dark-energy-like behavior.

• Matter-dominated epoch: \(\omega-1 \propto \rho/m_\omega^2\), suppressed unless \(m_\omega\) is ultralight.

In all cases, the system remains dynamically stable and approaches the GR limit whenever the potential has a quadratic minimum.

References