Kinematic Reconstruction of $Λ(t)$CDM Models

In this work, we have \textbf{analysed} two kinematic parametrizations for $Λ(t)$CDM models, namely, the linear expansions $Λ(z)=Λ_0+Λ_1z$ and $Q(z)=Q_0+Q_1z$, where $Q$ is the interaction term. In the case of the $Q(z)$ parametrization, we have also tested the particular case of a constant interaction term, $Q(z)=Q_0$. In order to constrain the free parameters of these models, we have used Cosmic Chronometers (CC), SNe Ia data (Pantheon+&SH0ES) and BAO data. As a general result, we have found weak constraints over the free parameters of the analysed models. In the case of $Λ(z)$, we have found for the $Λ$ variation parameter, $Ω_{\Lambda1}\equiv\frac{Λ_1}{3H_0^2}=0.02\pm0.14$. In the case of the $Q(z)$ parametrization, we have worked with the dimensionless interaction term $\gQ(z)\equiv\frac{8πGQ(z)}{3H_0^3}$, from which we have found $\gQ_0 = 2.2 \pm 2.7$ and $\gQ_1 = -6.2 \pm 7.6$. In the particular case of a constant interaction term, we have found $\gQ_0 = 0.18 \pm 0.7$. All these constraints are at 68% c.l. The constraints we have obtained are compatible with the standard $Λ$CDM model, although still providing a large margin for $Λ$ variation.

💡 Research Summary

In this paper the authors investigate phenomenological extensions of the standard ΛCDM cosmology in which the vacuum energy density Λ is allowed to vary with cosmic time and to exchange energy with the pressure‑less matter component. They adopt a purely kinematic reconstruction strategy: rather than deriving Λ(t) or the interaction term Q(t) from a specific underlying theory, they assume the simplest possible redshift dependence, namely a linear expansion in the observed redshift z. Two separate parametrizations are studied.

The first one treats the vacuum energy density as Λ(z)=Λ0+Λ1 z. Using the definition Q=−ṽΛ/(8πG) together with the continuity equation for matter, the authors derive analytic expressions for the matter density ρm(z) and the Hubble rate H(z). After normalising by the present‑day critical density they obtain a compact Friedmann equation that depends only on three dimensionless parameters: the present matter density Ωm, the linear Λ‑variation parameter ΩΛ1≡Λ1/(3H0²) and the usual flatness condition ΩΛ0=1−Ωm. In this formulation ΩΛ1 directly quantifies any deviation from a constant Λ.

The second parametrization focuses on the interaction term itself, writing Q(z)=Q0+Q1 z. Substituting this ansatz into the continuity equations and differentiating the Friedmann equation leads to a second‑order non‑linear differential equation for H(z). Because no closed‑form solution exists, the authors recast the problem as a first‑order dynamical system for the dimensionless Hubble function E(z)=H(z)/H0 and its derivative u(z)=dE/dz, with initial conditions E(0)=1 and u(0)=3Ωm/2 derived from the present‑day Friedmann equation. Numerical integration of this system yields H(z) for any chosen (Q0,Q1) pair.

To constrain the three free parameters in each model the authors combine three complementary data sets: (i) Cosmic Chronometers (CC), which provide direct measurements of H(z) from differential ages of passively evolving galaxies; (ii) the Pantheon+ & SH0ES compilation of Type Ia supernovae, supplying distance moduli over a wide redshift range; and (iii) Baryon Acoustic Oscillation (BAO) measurements, which give standard‑ruler distances and help break degeneracies between Ωm and the interaction parameters. They perform a Markov Chain Monte Carlo (MCMC) analysis, assuming flat priors and a spatially flat universe (k=0), and quote 68 % confidence intervals.

The resulting constraints are remarkably weak, reflecting the limited sensitivity of the current data to the chosen linear deviations. For the Λ‑variation model they obtain ΩΛ1=0.02±0.14, fully compatible with a constant Λ. For the interaction model they find Q0=2.2±2.7 and Q1=−6.2±7.6 (in the dimensionless form Q̃=8πGQ/(3H0³)). When the interaction is forced to be constant (Q1=0) the constraint tightens slightly to Q0=0.18±0.7. All these values are statistically indistinguishable from the ΛCDM limit Q0=Q1=0.

The authors discuss the implications of these findings. First, the large uncertainties indicate that the combination of CC, SNe Ia, and BAO data—while powerful for standard cosmology—does not yet have the discriminatory power to detect modest time‑dependence in Λ or a weak matter–vacuum coupling. Second, the linear‑in‑z ansatz is deliberately minimal; more complex functional forms (higher‑order polynomials, logarithmic terms, or model‑motivated expressions from scalar‑field theories) remain viable and could be explored once data precision improves. Third, the analysis assumes exact spatial flatness and neglects radiation, which is justified at low redshift but could bias results if high‑z CC measurements become available.

Methodologically, the paper’s strength lies in its clear separation of the two physical effects (Λ evolution versus interaction) and its transparent kinematic framework that avoids theoretical bias. However, the treatment of systematic uncertainties—particularly those inherent to the CC method (stellar population synthesis models, age‑dating techniques) and to supernova standardisation—could be elaborated further, as they may inflate the apparent parameter errors.

In conclusion, this work provides a proof‑of‑concept that Λ(t)CDM models can be confronted with data using a simple, model‑independent kinematic reconstruction. The current constraints are compatible with the standard ΛCDM paradigm but leave ample room for future, more precise observations (e.g., DESI, Euclid, JWST‑based high‑z chronometers) to tighten the bounds on ΩΛ1, Q0 and Q1, potentially revealing subtle dynamics in the dark sector.