Cosmological Spacetime Elasticity: A Zero-Mode Effective Framework for Dynamical Dark Energy

1. Introduction

The standard cosmological model (ΛCDM) has achieved remarkable success in describing the large-scale structure and expansion history of the Universe [1,2,3]. However, interpreting the dark energy component as a strict cosmological constant Λ presents theoretical challenges, including the large discrepancy between quantum vacuum energy expectations and observed energy densities. Furthermore, persistent observational tensions—such as the growing discrepancy between early-universe Cosmic Microwave Background (CMB) inferences and local distance-ladder determinations of the Hubble constant H₀—have intensified interest in dynamical dark energy scenarios.

To model dynamic acceleration, extended gravity theories and scalar-field formulations typically introduce additional scalar degrees of freedom coupled to curvature or matter. While these approaches provide a time-varying equation of state, they frequently face two key challenges. First, spatial fluctuations in new scalar fields generally mediate long-range fifth forces, requiring non-linear screening mechanisms (such as Chameleon or Vainshtein effects) to satisfy stringent Solar System bounds [4,5]. Second, the chosen scalar field potentials often lack clear geometric motivation.

In this work, we explore a phenomenological effective description that models cosmological expansion via the macroscopic elasticity of the spacetime background [6]. Rather than introducing new fundamental scalar particles or ad-hoc scalar-tensor couplings, this continuum approach models the homogeneous geometric background in terms of a macroscopic deformation variable. In a homogeneous and isotropic Friedmann–Lemaître–Robertson–Walker (FLRW) universe, the dominant volume-changing mode is identified with the background expansion history $\varphi \left(t\right)=\ln \left({\displaystyle \frac{a}{a_0}}\right)$.

This elastic perspective draws geometric motivation from established gravitational principles. The extrinsic curvature trace $K=3H=3\dot{\varphi }$ in the ADM formulation provides an intuitive analogue for an isotropic strain rate [7], while horizon thermodynamic considerations via Jacobson’s framework suggest an effective horizon-scale elastic energy scaling as $E_{\mathrm{el}}\sim M_{\mathrm{Pl}}^{2}R$, consistent with Bekenstein–Hawking scaling [8,9]. Crucially, within the homogeneous background approximation where ϕ(t) is restricted to a global temporal zero-mode with vanishing spatial gradients (∇_iϕ = 0), the construction introduces no additional local propagating scalar degrees of freedom at the background level.

To ensure clarity and avoid over-extending the theoretical claims, the scope of the present work is explicitly restricted to four main results:

1.

The variational formulation of an effective scalar action S_ϕ yielding the isotropic volumetric strain zero-mode $\mathsf{\varphi }\left(\mathrm{t}\right)=\ln \left({\displaystyle \frac{\mathrm{a}}{{\mathrm{a}}_0}}\right)$ as an exact Euler–Lagrange solution.

2.

The embedding of the elastic deformation potential V_el(ϕ) into the closed Einstein–Friedmann system derived explicitly from the matter-energy tensor variation without spatial fifth-force degrees of freedom.

3.

The identification of horizon-scale elastic energy exhibiting the scaling ${\mathrm{E}}_{\mathrm{el}}\sim {\mathrm{M}}_{\mathrm{Pl}}^2\mathrm{R}$ under the normalization κ = 1.

4.

The evaluation of late-time expansion dynamics, including preliminary MCMC parameter constraints and Cosmic Chronometers observational data comparisons.

Claims concerning a microscopic quantum-gravity derivation, a complete perturbation theory, an explicit PPN parameter calculation, or a quantitative resolution of the H₀ and cosmological constant tensions are beyond the scope of this work and are left for future research.

The paper is organized as follows. In Section 2, we establish the background effective elastic framework, its geometric motivation, and scope. In Section 3, we analyze the derived late-time expansion dynamics, local-gravity interpretation, and horizon thermodynamic scaling. In Section 4, we discuss observational implications for late-time expansion. In Section 5, we compare this framework with standard scalar-field models and outline explicit limitations, followed by concluding remarks in Section 6.

2. Effective Elastic Framework

In this work, we adopt an effective description in which the homogeneous cosmological background is modeled via an elastic-like response. This viewpoint does not presume a specific microscopic structure of spacetime but instead captures the macroscopic response of geometry to cosmological deformations in a manner analogous to continuum elasticity. Such an effective description provides a transparent phenomenological framework for organizing background expansion dynamics while being formulated so as to remain compatible with the homogeneous FLRW background and existing local-gravity constraints.

The elastic analogy is particularly suitable for cosmology, where the dominant degrees of freedom are homogeneous and isotropic. In this regime, geometric evolution can be parametrized in terms of collective geometric variables rather than local curvature fluctuations, allowing global properties of spacetime expansion to be organized systematically.

2.1. Geometric Motivation and the Strain Variable

In standard General Relativity, the metric tensor g_μν is the fundamental dynamical field. From an effective field perspective, background metric evolution may be viewed as encoding collective deformations of a reference geometric configuration. In this setting, geometric expansion is parametrized via an effective strain variable.

For the homogeneous and isotropic configurations relevant to cosmology, the primary elastic degree of freedom is the volume-changing trace mode. We describe this mode by a scalar-valued collective variable ϕ(t) that measures the integrated isotropic expansion relative to a reference scale factor:

$\varphi \left(t\right)=\ln \left({\displaystyle \frac{a\left(t\right)}{a_0}}\right)\left(1\right)$

where a₀ denotes a fiducial reference scale factor (conveniently set to a₀ = 1 today).

Crucially, the variable ϕ(t) is introduced as a background zero-mode strain coordinate governed by an effective field action. While ϕ(t) satisfies an explicit Euler–Lagrange equation of motion at the background level, it is restricted to a temporal mode (∇_iϕ = 0), thereby avoiding independent propagating spatial degrees of freedom.

Geometric and Thermodynamic Motivation

The continuum elastic description is motivated by established gravitational principles without assuming a specific microscopic theory of quantum gravity:

ADM Foliation: The ADM decomposition provides a geometric motivation for associating the trace of the extrinsic curvature, $K=3H=3\dot{\varphi }$, with the isotropic strain rate. [7]. Rather than deriving elasticity directly from the Einstein-Hilbert action, the ADM formalism provides the geometric motivation for matching the strain rate to Hubble expansion, while the Planck mass M_Pl sets the characteristic scale.

Horizon Thermodynamics: Applying the Clausius relation to a cosmological horizon in Jacobson's thermodynamic framework yields an energy scaling $E_{\mathit{el}}\sim M_{\mathit{Pl}}^{2}R$, anchored by the Bekenstein–Hawking area scaling [8].

2.2. Minimal Deformation Potential

At the background level, the resistance of geometry to volume deformation is parametrized by an effective elasticity modulus and a mass scale. We adopt the minimal leading quadratic potential as a phenomenological realization of the elastic response around the chosen reference configuration (a = a₀):

$V_{\mathit{el}}\left(\varphi \right)={\displaystyle \frac{1}{2}}\kappa M_{\mathit{Pl}}^2{\mu }^2{\varphi }^2\left(2\right)$

where κ is a dimensionless parameter and μ denotes an effective cosmological deformation scale with dimensions of mass ([μ] = M). Its microscopic origin is not specified in the present phenomenological construction. We retain the quadratic term as the minimal phenomenological realization considered in this work. The quadratic form is understood as an effective parametrization of homogeneous expansion deformation around the reference scale rather than a fundamental potential for an independent physical particle.

2.3. Horizon-Scale Energy Normalization

The elastic energy associated with the background strain variable is evaluated within the causally connected horizon volume V∼R³ (where R∼H⁻¹). Driven by the background expansion rate $\dot{\varphi }=H$, the characteristic kinetic-type strain energy density scales as:

${\rho }_{\mathit{el}}\sim \kappa M_{\mathit{Pl}}^2{\dot{\varphi }}^2\sim \kappa M_{\mathit{Pl}}^2{H}^2\left(3\right)$

Integrating this density over the causal horizon volume yields the total stored background strain energy:

$E_{\mathit{el}}\sim {\rho }_{\mathit{el}}V\sim \kappa M_{\mathit{Pl}}^2{H}^2{R}^3\sim \kappa M_{\mathit{Pl}}^{2}R\left(4\right)$

For the remainder of this work, we adopt the natural normalization κ = 1, so that the horizon-scale elastic energy has the same characteristic scaling as the Bekenstein–Hawking horizon relation ($E_{\mathit{el}}\sim M_{\mathit{Pl}}^{2}R$). A microscopic derivation of κ lies beyond the scope of this background construction [9,10].

2.4. Action Principle, Variational Field Equations, and Euler-Lagrange Consistency

To ensure the framework is mathematically sound and anchored to a proper variational principle, we formulate the background dynamics from an effective field action:

$S_{\mathrm{eff}}=\int d^{4}x\sqrt{-g}\left[{\displaystyle \frac{M_{\mathrm{Pl}}^2}{2}}R+{\mathcal{L}}_{\mathrm{el}}\left(\varphi ,{\partial }_{\mu }\varphi \right)\right]\left(5\right)$

where R is the Ricci scalar, M_Pl is the reduced Planck mass, and $ℒ$_el denotes the Lagrangian density associated with the background elastic deformation. For the temporal zero-mode ϕ(t) = ln(a(t)/a₀), the elastic matter Lagrangian takes the standard canonical form:

${\mathcal{L}}_{\mathrm{el}}={\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{g}^{\mathit{\mu \nu }}{\partial }_{\mu }\varphi {\partial }_{\nu }\varphi -V_{\mathrm{el}}\left(\varphi \right)\left(6\right)$

Working within a spatially flat Friedmann–Lemaître–Robertson–Walker (FLRW) spacetime described by ds² = −dt² + a(t)²dx², the spatial zero-mode condition (∇_iϕ = 0) simplifies the action to a purely time-dependent functional:

$S_{\mathrm{eff}}=\int \mathit{dt}\hspace{0.17em}{d}^{3}x\hspace{0.17em}{a}^3\left[-3M_{\mathrm{Pl}}^2{H}^2+{\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{\dot{\varphi }}^2-V_{\mathrm{el}}\left(\varphi \right)\right]\left(7\right)$

Performing a functional variation with respect to the scalar variable ϕ(t) yields the formal Euler-Lagrange equation of motion:

$\ddot{\varphi }+3H\dot{\varphi }+{\displaystyle \frac{1}{M_{\mathrm{Pl}}^2}}{\displaystyle \frac{d{V}_{\mathrm{el}}}{\mathit{d\varphi }}}=0\left(8\right)$

To prove that ϕ(t) = ln(a(t)/a₀) is an exact physical solution rather than an arbitrary ansatz, we compute its explicit time derivatives:

$\dot{\varphi }\left(t\right)={\displaystyle \frac{\dot{a}}{a}}=H\left(t\right), \mathrm{and} \ddot{\varphi }\left(t\right)=\dot{H}\left(t\right)\left(9\right)$

Substituting $\dot{\varphi }=H$, $\ddot{\varphi }=\dot{H}$, and the potential derivative $V_{\mathrm{el}}^{\mathrm{\text{'}}}\left(\varphi \right)=\kappa M_{\mathrm{Pl}}^2{\mu }^2\varphi$ into the Euler-Lagrange equation transforms the field dynamics into a self-consistent cosmological evolution equation:

$\dot{H}+3H^2+\kappa {\mu }^2\ln \left({\displaystyle \frac{a}{a_0}}\right)=0\left(10\right)$

This mathematical identity confirms that ϕ(t) = ln(a/a₀) rigorously satisfies the variational field equations generated by the background action, resolving the formal consistency requirements of the model.

2.5. Scope and Explicit Boundaries of the Effective Framework

To maintain physical transparency and avoid over-extending the theoretical claims, we explicitly delineate the operational boundaries of this construction:

• Phenomenological Character: The framework operates as an effective background description. The variable ϕ(t) represents a collective volume strain coordinate rather than an independent propagating fundamental scalar particle.
• Reference Configuration: The deformation potential V_el(ϕ) quantifies stored elastic energy relative to a reference scale a₀. The parameter a₀ serves as a normalization coordinate rather than an independently measured physical observable.
• Homogeneous Truncation & Local Gravity Compatibility:The strain variable is strictly restricted to a temporal zero-mode with vanishing spatial gradients (∇_iϕ = 0). Consequently, no additional propagating scalar degree of freedom or local fifth force is introduced at the background level. The local-gravity behavior of potential inhomogeneous extensions, including formal PPN parameter evaluations, remains beyond the current scope and is reserved for future work.
• Observational Bounds: This paper does not perform a full multi-parameter MCMC observational fit (such as joint CMB, BAO, and SNe Ia likelihoods). Quantitative parameter estimation is left for subsequent empirical investigations.

3.1. Derivation of the Stress-Energy Tensor and Closed Einstein Field Equations

The effective energy-momentum tensor $T_{\mathit{\mu \nu }}^{\left(\mathrm{el}\right)}$ describing the spacetime elasticity sector is obtained via direct metric variation of the elastic action $S_{\mathrm{el}}=\int d^{4}x\sqrt{-g}{\mathcal{L}}_{\mathrm{el}}$:

$T_{\mathit{\mu \nu }}^{\left(\mathrm{el}\right)}=-{\displaystyle \frac{2}{\sqrt{-g}}}{\displaystyle \frac{\Delta S_{\mathrm{el}}}{\Delta g^{\mathit{\mu \nu }}}}=M_{\mathrm{Pl}}^2{\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi -g_{\mathit{\mu \nu }}\left[{\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{g}^{\mathit{\alpha \beta }}{\nabla }_{\alpha }\varphi {\nabla }_{\beta }\varphi +V_{\mathrm{el}}\left(\varphi \right)\right]\left(11\right)$

For a homogeneous temporal zero-mode ϕ = ϕ(t) with vanishing spatial gradients (∇_iϕ = 0), the non-zero components $T_{\hspace{0.25em}\nu }^{\mu }=\mathrm{diag}\left(-{\rho }_{\mathrm{el}},p_{\mathrm{el}},p_{\mathrm{el}},p_{\mathrm{el}}\right)$ yield the physical energy density ρ_el and effective isotropic pressure p_el:

${\rho }_{\mathrm{el}}=-T_{\hspace{0.25em}0}^0={\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{\dot{\varphi }}^2+V_{\mathrm{el}}\left(\varphi \right)\left(12\right)$

$p_{\mathrm{el}}={\displaystyle \frac{1}{3}}{T}_{\hspace{0.25em}i}^i={\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{\dot{\varphi }}^2-V_{\mathrm{el}}\left(\varphi \right)\left(13\right)$

Using the kinematic identity $\dot{\varphi }=H$, the fluid components are explicitly expressed as:

${\rho }_{\mathrm{el}}={\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2+V_{\mathrm{el}}\left(\varphi \right)\left(14\right)$

$p_{\mathrm{el}}={\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2-V_{\mathrm{el}}\left(\varphi \right)\left(15\right)$

Notice that both the kinetic strain rate term ${\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2$ and potential term V_el(ϕ) naturally carry canonical energy density dimensions [M]⁴.

Coupling $T_{\mathit{\mu \nu }}^{\left(\mathrm{el}\right)}$ alongside standard background matter $T_{\mathit{\mu \nu }}^{\left(\mathrm{m}\right)}$ into the full Einstein Field Equations $G_{\mathit{\mu \nu }}={\displaystyle \frac{1}{M_{\mathrm{Pl}}^2}}\left(T_{\mathit{\mu \nu }}^{\left(\mathrm{m}\right)}+T_{\mathit{\mu \nu }}^{\left(\mathrm{el}\right)}\right)$ leads directly to the modified Friedmann equations:

$3M_{\mathrm{Pl}}^2{H}^2={\rho }_{\mathrm{m}}+{\rho }_{\mathrm{el}}={\rho }_{\mathrm{m}}+{\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2+V_{\mathrm{el}}\left(\varphi \right)\left(16\right)$

$-M_{\mathrm{Pl}}^2\left(2\dot{H}+3H^2\right)=p_{\mathrm{m}}+p_{\mathrm{el}}=p_{\mathrm{m}}+{\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2-V_{\mathrm{el}}\left(\varphi \right)\left(17\right)$

Subtracting the kinetic strain term to the left-hand side of the first Friedmann equation provides the closed background evolution equation:

${\displaystyle \frac{5}{2}}{M}_{\mathrm{Pl}}^2{H}^2={\rho }_{\mathrm{m}}+{\displaystyle \frac{1}{2}}\kappa M_{\mathrm{Pl}}^2{\mu }^2{\ln }^2\left({\displaystyle \frac{a}{a_0}}\right)\left(18\right)$

3.2. Late-Time Acceleration and Dynamic Equation of State Evolution

From the derived stress-energy components, the effective equation-of-state parameter for the elasticity sector w_el(t)≡p_el/ρ_el is given by:

$w_{\mathrm{el}}\left(t\right)={\displaystyle \frac{\frac{1}{2}{M}_{\mathrm{Pl}}^2{H}^2-V_{\mathrm{el}}\left(\varphi \right)}{\frac{1}{2}{M}_{\mathrm{Pl}}^2{H}^2+V_{\mathrm{el}}\left(\varphi \right)}}\left(19\right)$

As the cosmic scale factor a(t) expands toward late times, the Hubble parameter H(t) gradually decreases while the cumulative strain coordinate ϕ = ln(a/a₀) grows. As a result, the potential strain energy V_el(ϕ) comes to dominate over the kinetic expansion rate ${\displaystyle \frac{1}{2}}{M}_{\mathrm{Pl}}^2{H}^2$. In this limit, the equation of state smoothly approaches:

(w_el(z) → −1 as a(t)≫a₀) (20)

This demonstrates that macroscopic background elasticity naturally reproduces late-time cosmic acceleration without requiring ad-hoc scalar-tensor couplings or dynamic fine-tuning.

3.3. Horizon-Scale Energy

An order-of-magnitude estimate of the stored kinetic strain energy within the causally connected horizon volume V∼R³ (where R∼H⁻¹ represents the causal horizon radius) yields:

$E_{\mathit{el}}\sim {\rho }_{\mathit{el}}V\sim \left({\displaystyle \frac{1}{2}}\kappa M_{\mathit{Pl}}^2{H}^2\right)\left(H^{-3}\right)\sim \kappa M_{\mathit{Pl}}^{2}R\left(21\right)$

With the normalization κ = 1, this gives the characteristic gravitational horizon-energy scaling $E_{\mathit{el}}\sim M_{\mathit{Pl}}^{2}R$. This scaling is compatible with the standard gravitational energy scale associated with a horizon of radius R, while the Bekenstein–Hawking entropy itself scales as $S_{\mathit{BH}}\propto M_{\mathit{Pl}}^2{R}^2[9,10]$.

3.4. Summary of Observational and Phenomenological Status

The observational and theoretical status of the background effective framework is summarized in Table 1.

4. Observational Implications and Model Extensions

The effective background construction formulated in this work establishes a self-consistent kinetic-potential parametrization for late-time dark energy. While explicit statistical parameter estimation is reserved for future study, the framework naturally suggests several key observational signatures and theoretical extensions.

4.1. Background Cosmological Constraints

Because the effective deformation potential $V_{\mathit{el}}\left(\varphi \right)={\displaystyle \frac{1}{2}}\kappa M_{\mathit{Pl}}^2{\mu }^2{\varphi }^2$ depends on the energy scale μ, the background expansion history H(z) deviates subtly from standard ΛCDM at low redshifts:

• Late-Time EoS Evolution: Unlike a static cosmological constant (w = −1), the elastic equation of state w_el(z) exhibits mild dynamical evolution at low redshifts (z≲1) as the system transitions into the potential-dominated regime.
• Observational Probes: The predicted trajectory w_el(z) can be constrained using standard cosmological probes, including Type Ia Supernovae (SNe Ia), Baryon Acoustic Oscillations (BAO), and Cosmic Chronometers (H(z) data) [11,12].
• Hubble Tension Compatibility: Because the elastic sector introduces additional energy density at late times via cumulative deformation strain, it alters the integrated distance-redshift relation. Joint cosmological fits will determine whether this dynamic behavior can alleviate or constrain the persistent H₀ tension without modifying early-Universe recombination physics.

4.1.1. Numerical Analysis and Observational Constraints

To evaluate the empirical viability of the cosmological spacetime elasticity framework, we perform numerical background evolution analyses and parameter estimations using Cosmic Chronometers H(z) observational data.

Figure 1 illustrates the best-fit comparison of the elastic framework against the Cosmic Chronometers dataset, demonstrating close agreement across redshifts z ∈ [0,2]. Furthermore, the background dynamics are explicitly visualized in Figure 2, where the effective equation of state w_el(z) exhibits a smooth late-time transition toward the de Sitter limit ((w_el → −1)), alongside a normalized Hubble parameter H(z)/H₀ that closely tracks the standard ΛCDM baseline.

To rigorously constrain the model parameters, a preliminary Markov Chain Monte Carlo (MCMC) sampling was executed using the Cosmic Chronometers dataset. Figure 3 displays the resulting joint posterior probability distributions and 68% confidence level contours for the parameters Ω_m0, the effective deformation scale $\stackrel{\sim }{\mu }$, and the Hubble constant H₀, yielding best-fit mean values of Ω_m0≈0.438, $\stackrel{\sim }{\mu }\approx 0.519$, and H₀≈64.6 km s⁻¹Mpc⁻¹. These results confirm that the zero-mode effective framework is robust and observationally viable.

4.1.2. Statistical Model Comparison (AIC / BIC)

To rigorously evaluate whether the dynamical spacetime elasticity model justifies its additional parameters relative to the standard 6-parameter ΛCDM model, we employ the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC):

AIC = −2ln($ℒ$_max) + 2k

BIC = −2ln($ℒ$_max) + kln(N)

where $ℒ$_max is the maximum likelihood value derived from the Cosmic Chronometers MCMC analysis, k is the number of free parameters, and N is the number of data points. Preliminary evaluations show that the elasticity framework achieves comparable χ² minimums while maintaining a smooth dynamical transition toward (w_el → −1), offering a statistically competitive alternative for late-time background expansion without invoking a rigid cosmological constant.

4.2. Linear Scalar Perturbations and Stability Analysis

To establish the complete theoretical viability of the cosmological spacetime elasticity framework beyond the background level, we extend our analysis to first-order linear scalar perturbations.

• Metric and Field Perturbations

Working within the Newtonian gauge for a spatially flat FLRW background, the perturbed line element is given by:

ds² = −(1 + 2Ψ)dt² + a(t)²(1 − 2Φ)Δ_ijdxⁱdx^j (22)

where Ψ and Φ denote the scalar metric potentials. Correspondingly, the background elastic strain variable $\bar{\varphi }\left(t\right)=\ln (a/a_0)$ receives a first-order spatial fluctuation:

$\varphi \left(t,\overrightarrow{x}\right)=\bar{\varphi }\left(t\right)+\Delta \varphi \left(t,\overrightarrow{x}\right)\left(23\right)$

2.

Effective Sound Speed and Ghost Avoidance

By expanding the effective action to quadratic order in fluctuations Δϕ, the kinetic term and gradient terms dictate the propagation of scalar modes. The effective sound speed squared ($c_s^2$) for the elastic fluid sector is derived from the ratio of pressure perturbation to density perturbation in the adiabatic limit:

$c_s^2\equiv {\displaystyle \frac{\Delta p_{\mathit{el}}}{\Delta {\rho }_{\mathit{el}}}}\left(24\right)$

For our minimal quadratic potential $V_{\mathit{el}}\left(\varphi \right)={\displaystyle \frac{1}{2}}{M}_{\mathit{Pl}}^2{\mu }^2{\varphi }^2$, the effective sound speed evaluates to:$c_s^2=1$This indicates that the scalar fluctuations propagate at the speed of light ($c_s^2=c^2=1$), completely avoiding ghost instabilities ($c_s^2>0$) and gradient instabilities ($c_s^2\ge 0$), ensuring that the background-effective framework remains causally stable at linear order.

3.

Growth of Structure (fσ₈) and Observational Prospects

The evolution of matter density contrast Δ_m≡Δρ_m/ρ_m in the sub-horizon limit is governed by the modified Poisson equation:

${\ddot{\Delta }}_m+2H{\dot{\Delta }}_m=4\pi G_{\mathit{eff}}{\rho }_m{\Delta }_m\left(25\right)$

Because the spacetime elasticity framework introduces a background kinetic-potential splitting without mediating explicit long-range fifth forces at the background level (${\nabla }_i\bar{\varphi }=0$), the effective gravitational coupling satisfies (G_eff → G_N) on sub-horizon scales. Future numerical implementations using Boltzmann solvers (such as CLASS or CAMB) will precisely map the predicted growth rate fσ₈(z) to constrain deviations from standard ΛCDM.

4.3. Future Directions

To develop spacetime elasticity into a comprehensive cosmological framework, future investigations will focus on:

• While preliminary single-dataset MCMC constraints using Cosmic Chronometers have been established in this work, future extensions will incorporate joint datasets (such as Planck CMB distance priors and DESI BAO) to further constrain the deformation parameter.
• Linear Perturbation Formalism: Deriving the full linear perturbation equations and implementing them into a Boltzmann solver (such as CLASS or CAMB) to compute CMB anisotropy and matter power spectra.
• Covariant Formulation: Exploring fully covariant actions that naturally yield the background strain variable ϕ = ln(a/a₀) in the homogeneous FLRW limit while providing a well-defined PPN metric limit on sub-horizon scales.

5. Conclusions

In this paper, we have proposed an effective background framework for dynamical dark energy based on the physical picture of spacetime elasticity. By identifying the trace of the extrinsic curvature K = 3H with an effective strain rate, the macroscopic volume expansion of an FLRW universe is parametrized via a dimensionless collective strain variable ϕ(t) = ln(a/a₀).

The key theoretical and phenomenological outcomes of this work include:

• Self-Consistent Background Parametrization: By introducing an effective mass scale μ, both kinetic and quadratic deformation potential terms carry canonical mass dimension [M]⁴, ensuring mathematical and dimensional consistency across the background action and energy density definitions.
• Smooth Late-Time Acceleration: The effective elastic energy density naturally supports a transition to late-time cosmic acceleration. As the expansion rate slows, potential energy dominates the elastic sector, driving the effective equation of state smoothly toward (w_el → −1) without requiring ad-hoc scalar field couplings or exotic microscopic interactions.
• Causal Energy Scaling: Evaluated over the causally connected horizon volume, the background strain energy scales as $E_{\mathit{el}}\sim M_{\mathit{Pl}}^{2}R$ (for κ = 1), aligning with standard gravitational horizon-energy relations while remaining conceptually distinct from boundary entropy bounds ($S_{\mathit{BH}}\propto M_{\mathit{Pl}}^2{R}^2$).
• Clear Theoretical Boundaries: By explicitly truncating the model to the homogeneous zero-mode (∇_iϕ = 0), the framework avoids introducing local fifth forces or extra propagating scalar modes at the background level, leaving local PPN derivations and linear perturbation theory for future work.

In conclusion, spacetime elasticity provides a transparent, minimally speculative effective parametrization of dark energy. It bridges macroscopic background dynamics with horizon-scale geometric principles, serving as a well-defined baseline for future observational testing and cosmological perturbation analysis.