Barrow's VSL Cosmology

• Barrow's VSL theory is a cosmological framework that models the speed of light as a time-dependent parameter, fundamentally altering universal dynamics.
• It employs a power-law relation, c(t) = c₀ aⁿ, where the exponent n impacts key observables like cosmic distances and redshift drift.
• Empirical constraints show n is nearly zero, indicating the current universe exhibits minimal deviations from a constant speed of light.

Barrow's Varying Speed of Light (VSL) theory postulates a time-dependent speed of light as a solution to foundational problems in cosmology, such as the horizon and flatness problems, and as an alternative to the standard cosmological constant-driven late-time acceleration paradigm. The canonical Barrow model posits a power-law dependence of the vacuum speed of light on the cosmic scale factor, expressed as $c(t) = c_0 a^n(t)$ (or equivalently $c(z) = c_0(1+z)^{-n}$), with the exponent $n$ as a fundamental parameter. This ansatz modifies key cosmological observables and introduces distinctive dynamics in the evolution of the universe, affecting both background and perturbative quantities.

1. Theoretical Framework of Barrow's VSL Model

Barrow's VSL theory operates within a spatially flat Friedmann-Robertson-Walker (FRW) metric where the speed of light is promoted to a cosmic scalar field or parameterized function of time: $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$ The modified Einstein equations take the form

$G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$

leading to altered Friedmann and acceleration equations: $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$ The novel term $-2H\dot{c}/c$ is absent in standard cosmology and encodes the impact of a time-varying $c$.

Energy-momentum conservation is also revised: $$\dot{\rho} + 3H\left(\rho + \frac{p}{c^2}\right) = \frac{3H^2}{4\pi G} \frac{\dot{c}}{c},$$ where changes in $c$ act as a source term for cosmic "matter creation" (Qi et al., 2014).

The power-law ansatz is central: $c(z) = c_0(1+z)^{-n}$0 with $c(z) = c_0(1+z)^{-n}$1 controlling the direction and rate of secular evolution of $c(z) = c_0(1+z)^{-n}$2. In generalized formulations, $c(z) = c_0(1+z)^{-n}$3 may be written as $c(z) = c_0(1+z)^{-n}$4 with $c(z) = c_0(1+z)^{-n}$5 (Lee, 2024), or via alternative parametric forms such as the CPL-style and "magnetically-triggered" transitions (Salzano et al., 2016).

2. Modifications to Cosmological Observables and Kinematics

A variable $c(z) = c_0(1+z)^{-n}$6 modifies the computation of cosmological distances. For luminosity distance,

$c(z) = c_0(1+z)^{-n}$7

where the dimensionless Hubble parameter $c(z) = c_0(1+z)^{-n}$8 is changed by both the time dependence of $c(z) = c_0(1+z)^{-n}$9 and altered continuity equations: $n$0 with pressureless matter $n$1 and dark energy $n$2 given by the above integral.

For redshift-drift observables (Balcerzak et al., 2013), the rate is

$n$3

which introduces an $n$4-dependent correction: for $n$5, dust components acquire negative pressure and the cosmological constant becomes phantom-like; $n$6 boosts CDM-like behavior.

The luminosity–distance–redshift relation is further modified in the power-law framework: $n$7 along the special locus $n$8, which empirically emerges in supernova fits (Nguyen, 9 Jan 2026).

3. Observational Constraints and Empirical Performance

Comprehensive likelihood analyses utilizing supernova Ia (Union 2.1, Pantheon), BAO, OHD, and CMB shift parameters yield stringent bounds on the allowed variation of $n$9 (Qi et al., 2014, Nguyen, 2020, Nguyen, 9 Jan 2026). With power-law models (Barrow ansatz), the best-fit exponent is extremely small: $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$0 indicating near-perfect constancy of $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$1 over the observable universe. Reconstruction of $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$2 with this bound shows that for redshift $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$3, the variation is negligible ($$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$4), and at the CMB recombination epoch ($$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$5), the deviation is only $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$6 (Qi et al., 2014).

Stochastic approaches using BAO and cosmic chronometer data (covering $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$7) further emphasize the statistical rejection of the classical Barrow-VSL model; the reduced chi-square and AIC/BIC metrics favor a strictly constant speed of light (Zhang et al., 2024).

However, alternative parametrizations and the inclusion of galaxy-scale effects (local expanding systems, "yardstick" correction $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$8) show high-likelihood degeneracies along $$ds^2 = -c^2(t)dt^2 + a^2(t)\left[dr^2 + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)\right].$$9, yielding empirical fits comparable in quality to standard $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$0CDM (Nguyen, 2020, Nguyen, 9 Jan 2026). This degeneracy implies a universal synchrony between $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$1 and $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$2: $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$3 which has profound kinematic consequences absent in $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$4CDM (Nguyen, 9 Jan 2026).

4. Physical and Cosmological Implications

Barrow's VSL models have far-reaching implications for classical cosmology:

• Late-time acceleration without $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$5: Along the empirical synchrony $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$6, high-$G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$7 supernovae appear dimmer due to modified kinematics, reproducing the apparent acceleration without invoking dark energy (Nguyen, 9 Jan 2026, Nguyen, 2020).
• Horizon and flatness problems: Early-universe epochs with large (positive) $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$8 can prevent the formation of particle/event horizons, removing the necessity for inflation; comoving integrals diverge globally under the synchrony law (Nguyen, 9 Jan 2026).
• Resolution of Hubble tension: Allowing for monotonic variation in the local gravitational scale for bound objects induces a shift in the effective $G_{\mu\nu} = \frac{8\pi G}{c^4(t)} T_{\mu\nu},$9 determined at different redshifts, consistent with the observation that high-$$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$0 estimates are $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$1 below local values (Nguyen, 2020).
• Generalized Copernican Principle: The condition $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$2 leads to cosmological self-invariance in time, with the Riemann tensor and Ricci scalar independent of the particular epoch (Nguyen, 9 Jan 2026).
• Novel conformally flat metric: By appropriate rescaling under the synchrony law, the metric becomes manifestly conformal to Minkowski space, eliminating built-in cosmological horizons (Nguyen, 9 Jan 2026).

5. Comparison with Minimal and Extended VSL Models

Modern refinements, such as the minimally extended VSL (meVSL) model (Lee, 2024), retain the essential Barrow parametrization ($$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$3), but require compensatory variation of gravitational coupling ($$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$4) to preserve the Einstein-Hilbert action constant ($$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$5). This avoids explicit breaking of local Lorentz invariance and maintains the conservation law ($$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$6), with all dimensional constants co-varying to ensure operational consistency at each epoch.

The meVSL scenario produces algebraic corrections to Friedmann, continuity, and all observables, directly testable by cosmic chronometers, distance duality, and SNeIa time-dilation. Current constraints allow $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$7–0.3, consistent with constancy but not excluding an O(10\%) cosmic drift (Lee, 2024). In contrast, the original Barrow–Magueijo models often entail strong Lorentz-violation, bimetric structures, and explicit energy non-conservation, which are typically unconstrained by late-time data (Salzano et al., 2016).

6. Observational Diagnostics and Distinctive Predictions

Barrow's VSL modifies observable diagnostics such as:

• Geometrical diagnostic $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$8: In VSLDE, $$3H^2 = 8\pi G \rho, \quad 2\frac{\ddot{a}}{a} + H^2 - 2H\frac{\dot{c}}{c} = -\frac{8\pi G}{c^2} p.$$9 remains nearly constant and indistinguishable from $-2H\dot{c}/c$0CDM for small $-2H\dot{c}/c$1 (Qi et al., 2014).
• Angular-diameter distance maxima: The location $-2H\dot{c}/c$2 of $-2H\dot{c}/c$3 maximum shifts in VSL models, permitting precision tests against $-2H\dot{c}/c$4CDM in planned BAO and cosmic chronometer surveys (Salzano et al., 2016).
• Phantom dark energy signatures: For $-2H\dot{c}/c$5, effective equations of state (e.g., $-2H\dot{c}/c$6) become phantom ($-2H\dot{c}/c$7), altering the evolution of cosmic components and potentially impacting CMB and nucleosynthesis observables (Balcerzak et al., 2013).
• Redshift drift: The predicted drift rates are below the current detection threshold unless $-2H\dot{c}/c$8, rendering future extremely high-precision observations necessary for empirical discrimination (Balcerzak et al., 2013).

7. Current Status and Prospects for Future Constraints

Empirical studies using SNIa, BAO, CMB, and chronometer datasets deliver a robust constraint $-2H\dot{c}/c$9, indicating negligible cosmological variation of the speed of light to within $c$0 at all redshifts probed (Qi et al., 2014, Zhang et al., 2024). Bayesian evidence comparisons reveal that, despite the theoretical flexibility, constant-$c$1 models provide statistically superior fits across most scenarios, with only specific, extended VSL models attaining substantial evidence in favor over $c$2CDM (Salzano et al., 2016).

Future redshift-drift surveys, 21cm intensity mapping, gravitational-wave standard sirens, and high-precision BAO will further tighten constraints, with the potential to probe $c$3 at $c$4 level or better. Observational signatures such as the synchronized law $c$5 and the absence of cosmic horizons, if confirmed, would necessitate reformulations of gravitational dynamics fundamentally distinct from $c$6CDM (Nguyen, 9 Jan 2026).

---

In summary, Barrow's varying speed of light cosmology establishes a rigorous alternative to standard-model cosmology through the introduction of the simple ansatz $c$7, fundamentally altering the kinematic and dynamic structure of the universe. While the present observational epoch strongly favors a constant $c$8, ongoing empirical studies and future precision measurements remain pivotal in testing its cosmological validity and probing the conceptual boundaries of gravitational theory.