1. Introduction
For over thirty years, cosmologists have accumulated compelling evidence that cosmic expansion is accelerating. More specifically, this acceleration appears to be dominated by a cosmological constant term, commonly denoted as .
This
term can be understood in three distinct ways: (1) as a fundamental (or effective) modification of Einstein’s classical general relativity (GR), denoted
, (2) as an effective source term from a dark energy (DE) component,
, or (3) as a boundary term,
. These three possible origins can be illustrated by writing the Einstein–Hilbert action for classical GR with the corresponding additional terms:
The first two terms in the first integral represents the classical GR Lagrangian with matter-energy content
as source term [
1]. The third term corresponds to the fundamental cosmological constant
[
2]. The fourth term is the DE or quintessence source term, given by
, for a single effective scalar field
. For a homogeneous perfect fluid with density
and pressure
in its ground state (where
), this reduces to
(see also Equations (B66)–(B68) in [
3]). This illustrates how
and
provide totally degenerate interpretations for
.
The last integral represents the Gibbons–Hawking–York (GHY) boundary term [
4,
5,
6], where
K is the trace of the extrinsic curvature at the boundary
. As shown in [
7], for a finite Friedman–Lemaître–Robertson–Walker (FLRW) metric with total mass–energy
, this term results in:
where
, with
. This provides a fundamentally different origin for
compared to
or
. As noted in [
8], the
boundary term cancels the potential contributions of both
and
, solving the fine tuning and coincidence problems [
9].
Regardless of its origin, the
term is typically interpreted as a repulsive force between galaxies that counteracts gravity, driving the accelerated expansion of the universe. This phenomenon is often cited as one of the most profound challenges in contemporary physics and may offer a crucial observational pathway to understanding quantum gravity (e.g., see [
10] and references therein).
Cosmic acceleration is typically measured using the adimensional coefficient
q, defined as
, where
. If the universe follows an equation of state with
, this leads to
. For regular matter or radiation, where
, we would expect deceleration in the expansion (
) due to gravity. However, measurements from various sources such as galaxy clustering, Type Ia supernovae, and CMB consistently show an expansion that is asymptotically approaching
or
(e.g., see [
11] and references therein for a review of more recent results, including weak gravitational lensing). This aligns well with a cosmological constant
, where
approaches
and
q approaches 1. What does this mean?
The term dark energy (DE) was introduced by [
12] to refer to any component with
. However, there is no fundamental understanding of what DE is or why we measure a term with
. A natural candidate for DE is
, which is equivalent to
and can also be thought of as the ground state of a field (the DE), similar to the inflation but with a much smaller (≃
) energy scale.
can also be a fundamental constant in GR, but this has some other complications ([
13,
14,
15]).
This paper critically examines the conventional concept of cosmic acceleration and proposes an alternative framework for understanding cosmic expansion dynamics. In
Section 2, we establish the notation and derive standard definitions for cosmic expansion in the comoving frame. In
Section 3, we demonstrate the dependence of these standard definitions on the observer’s frame, highlighting the lack of covariance and potential for misinterpretation in the commonly used concept of cosmic acceleration.
Section 4 and
Section 5 introduce an alternative definition for cosmic acceleration, which is grounded in the
event horizon. In
Section 6, we compare both definitions to observational data, demonstrating that our proposed approach offers greater consistency with empirical observations.
Appendix A provides a detailed exposition of the correct method for defining 4D acceleration in relativity based on the geodesic deviation equation. We also elaborate on the idea that
corresponds to a friction (attractive) force that decelerates cosmic events and revisit the Newtonian limit to show that
corresponds to an additional (attractive) Hooke’s term to the inverse square gravitational law, envisioning a “rubber band Universe”.
Finally, we conclude with a summary and discussion, emphasizing the significance of our findings for cosmological theory and observational practice and suggesting avenues for further research and exploration in the field.
3. De Sitter Phase
The FLRW metric with
asymptotically tends to a constant:
, which corresponds to exponential inflation and de Sitter metric, which can also be written as:
This form corresponds to a static 4D hypersphere of radius
. In this rest frame, events can only travel a finite distance
within a static 3D surface of the imaginary 4D hypersphere.
This implies that there exists a frame duality, allowing us to equivalently describe the de Sitter space either as static in proper or physical coordinates or as exponentially expanding in comoving coordinates . In the static frame , events are constrained within a limited region of the hypersphere, while in the expanding frame , distances and coordinates evolve with time following an exponential expansion characterized by the de Sitter horizon .
This frame duality can be understood as a Lorentz boost that results in both length contraction and time dilation. If we define the coordinate
, the radial velocity gives us the Hubble law
, leading to a Lorentz factor given by
where
. In the rest frame
, an observer sees the moving fluid element
contracted by the Lorentz factor
in the radial direction and experiences a time dilation by
: i.e.,
and
. More formally, we need to find a change of variables from comoving coordinates
in the FLRW metric of Equation (
3) to rest-frame de Sitter coordinates
, where
(see [
16]):
which agrees with Equation (
7) in [
17], with:
where
. This form reproduces the static de Sitter metric Equation (
6) when
. It also shows that
t retains its time-like character as we cross inside
. This is to be contrasted with the event horizon of the Schwarzschild metric, which requires a change of variables as we cross inside the horizon.
The spatial part of the light element in Equation (
8) is illustrated in
Figure 2. Geometrically, it corresponds to the metric of a hypersphere of radius
that expands towards a constant radius
, which corresponds to an event horizon (see also
Section 5 below and Appendix B in [
16]). In the above rest (de Sitter) frame, the FLRW background is asymptotically static, indicating no expansion or acceleration, while in the comoving frame, there is cosmic acceleration (
). This observation highlights that the concept of cosmic acceleration commonly used in cosmology critically depends on the chosen frame of reference.
4. Event Acceleration
The interpretation of cosmic acceleration in Equation (
5) is solely based on the definition for acceleration
in Equation (
5). We next show that such a definition corresponds to events without a cause-and-effect connection, and this lead us to the wrong picture of what is happening. We will then introduce a more physical alternative.
Consider the distance between two events corresponding to the light emission of a galaxy at
and the reception somewhere in its future
. The photon travels following an outgoing radial null geodesic
, which, from Equation (
3), implies
. This situation is depicted in
Figure 3. We can define a 3D space-like distance
d based in the comoving separation
:
This is, in fact, the distance that corresponds to the acceleration given by
in Equation (
5) because
and
, where the derivative is with respect to
, the time at emission. Such a distance corresponds to the distance between
and
so that
. These events lack causal connection and are beyond observation. While using
d is not inherently incorrect, it involves extrapolating observed events (like luminosity distance) into non-observable realms. Essentially,
d aligns with a non-local theory of gravity or the Newtonian approximation, where action at a distance occurs with an infinite speed of light.
We can instead use the the distance traveled by the photon:
Note that we use units of
, so this should be read as
. But cosmic acceleration is zero
for such a distance because
.
So the usual definition currently used by cosmologists in Equation (
10) corresponds to events that are space-like, i.e., at a fixed comoving separation or fixed cosmic time
. It only takes into account the change in the distance due to the expansion of the universe. To have a measurement of cosmic acceleration that is closer to actual observations, we need to use the distance between events that are causally connected, i.e., that not only takes into account how much the universe has expanded but also how long it has taken for the two events to be causally connected.
To this end, we should use the proper future light-cone distance
obtained from
in the FLRW of Equation (
3) as
and
(see, e.g., Equation (A1) in [
18]):
Note that the term with the integral is not
, but it corresponds to the coordinate distance
traveled by light between the two events, including the effect of cosmic expansion. Thus, we argue that we should use
R instead of
d in Equation (
10) as a measure of the distance in cosmology to define the cosmic acceleration and expansion rate. The differences between these three distances are illustrated in
Figure 3.
Using
R as a distance is equivalent to a simple change of coordinates in the FLRW metric of Equation (
3) from comoving coordinates
to physical coordinates
:
which is just Minkowski’s metric in spherical coordinates with a radius
.
We then have that
, and we define the expansion rate between null events as:
where the additional term
corresponds to a friction term. There is ambiguity in this definition because
R in Equation (
12) depends also on the time
used to define
R. To break this ambiguity, we arbitrarily fix
R to be the distance to
(which corresponds to a possible future event horizon):
where
. The quantity
corresponds to a friction term because it opposes the expansion
H that generates the term. It just originates from the change of frame. As we will see in the next section, this choice implies that
is zero unless
. So this new invariant way to define cosmic expansion reproduces the standard definition when
. But for
, we have that the event expansion halts
(blue line in
Figure 1) due to the friction term (red line) for
, while the standard Hubble rate definition approaches a constant
(black line). This might seem irrelevant at first look, but the resulting physical interpretation is quite different. In the standard definition of
H, the expansion with
becomes asymptotically exponential (or inflationary expansion), while in our new definition of
, the expansion becomes static (as in the static de Sitter metric).
The event acceleration can then be measured as:
The correct way to define 4D acceleration in relativity is based on the geodesic deviation equation: Equation (
A1). The relation to
q and
will be discussed in
Appendix A.
As before, for
, the friction term
makes little difference between
q and
. For
, the friction term asymptotically cancels the
term in
(i.e., Equation (
5)) so that
is always negative, no matter how large
is (
and
). The net effect of the
term is to bring the expansion of events to a faster stop (
) than in the case with gravity alone. This is illustrated in
Figure 1. The
term produces a faster deceleration (than with gravity alone). This corresponds to an attracting (and not repulsive) force, as explained in more detail in
Appendix A.
5. Event Horizon
What is more relevant to understand the meaning of
is that the additional deceleration brings the expansion to a halt within a finite proper distance between the events, creating an event horizon (EH). The EH is the maximum distance that a photon emitted at time
can travel following the outgoing radial null geodesic:
This is illustrated in
Figure 4, which also demonstrates how inflation and the horizon problem (i.e., the observation that CMB measurements detect super-horizon
frozen perturbations) occurred within
. All the observable Universe (green line) is contained inside
, and we can therefore not measure anything outside. In particular, we can not measure any inhomogeneities or anisotropies outside
R from the inside, even if the density is discontinuous at
R [
16].
For
, we have
, so there is no EH. But for
, we have that
(red line in
Figure 1). We can then see that
corresponds to a causal horizon or boundary term. The analog force behaves like a rubber band between observed galaxies (null events) that prevents them from crossing some maximum stretch (i.e., the EH). We can interpret such a force as a boundary term that just emerges from the finite speed of light (see
Appendix A and Eq.
1).
The FLRW metric with
asymptotically tends to the de Sitter metric in Equation (
6). This form corresponds to a static 4D hypersphere of radius
. So in this (rest) frame, events can only travel a finite distance
within a static 3D surface of the imaginary 4D hypersphere. The region inside
is causally disconnected from the outside. In the context of the FLRW framework, this condition corresponds to
, where
is a radial (space-like) distance. This condition implies that the expansion interpretation is valid only as long as
, indicating that it does not make sense for larger values where we cannot transition from
to
. Essentially, beyond this threshold, the cosmological interpretation of expansion breaks down due to the causal disconnection imposed by the horizon defined at
.
As shown in Equation (
8), this frame duality can be understood as a Lorentz boost. An observer in the rest frame sees the moving fluid element
contracted by the Lorentz factor
. This duality is better understood using our new measures for the expansion rate
and cosmic deceleration
based on the distance between causal events.
6. Comparison to Data
We show next how to estimate the new measure of cosmic acceleration
using direct astrophysical observations. As an example, consider the Supernovae Ia (SNIa) data as given by the “Pantheon Sample” compilation ([
19]) consisting of 1048 SNIas between
. Each SNIa provides a direct estimate of the luminosity distance
at a given measured redshift
z. This corresponds to the comoving look-back distance:
so that
directly gives us
. The second derivative gives us the acceleration:
is given by the model prediction in Equation (
17) (arbitrarily fixed at
in both data and models). We adopt here the approach presented in [
20], which used an empirical fit to the luminosity distance measurements based on a third-order logarithmic polynomial:
where
. The authors of [
20] find a good fit using
and
for the full SNIa “Pantheon Sample”. We use these values of
A and
B and their corresponding errors to estimate
H,
q, and
using the above relations. The results for
H and
q are shown as shaded cyan regions in the left panel of
Figure 5. They are compared to the LCDM predictions (dashed lines) in Equations (
4) and (
5).
There is very good agreement in
for
. At
, the
estimates are also consistent with the
predictions. But the detailed
evolution with redshift in the SNIa data does not seem to follow any of the model predictions, especially for
. The
estimates are too steep compared to the different models’ predictions. If we compare instead the
estimates (see right panel of
Figure 5), we find much better agreement with the model predictions. This seems to validate our
approach, but it is not clear from this comparison alone if this is caused by the fitting function used in Equation (
20).
To test this further, we use measurements of the radial BAO data to estimate
. Such measurements give us a direct estimate of
(as first demonstrated by [
21]), so they have the advantage over SNIa in that we only need to do a first-order derivative to estimate
q or
:
As an illustration, we utilize the
measurements presented in Table 2 of [
22], which provide an expanded dataset that includes results up to 2023. This compilation of
values is shown as red points with
error bars in the left panel of
Figure 5 and is labeled as “<2023”. Additionally, we incorporate four new data points from DESI2024 (LRG1, LRG2, LRG3+ELG1, and ELG2) derived from Table 18 of [
23], using a sound horizon
from CMB measurements [
24] to obtain
, along with one data point from the Ly-alpha forest as described in Equation (7.3) of [
25]. These are displayed as blue points with
error bars in the left panel of
Figure 5 and are labeled as “DESI2024”.
The full compilation includes measurements from galaxy clustering () and the Ly-alpha forest in quasars (). The combination of these distinct redshift ranges provides a robust measurement of at intermediate redshifts (), where discrepancies in the supernova Type Ia data are observed when comparing the traditional q model and the model (as discussed earlier). The radial BAO offers strong constraints on cosmic acceleration that are independent of potential calibration errors in or sampling biases from small-area surveys. This level of precision is something that is not yet achievable with current supernova datasets, but it will be exciting to see how upcoming wider and deeper surveys might address this issue in the near future.
We fit a quadratic polynomial with inverse variance weighting to the radial BAO data:
In units of Km/s/Mpc, we find
,
and
, with strong covariance between the errors (the cross-correlation coefficient between
and
is
). The value of
is in good agreement with the Planck CMB fit [
24] but is in some tension with the SNIa local calibration:
(see [
26]). This corresponds to either a local calibration problem (in SNIa, in radial BAO, or in both) or tension in the
CDM model at different times or distances (see, e.g., [
27]). We ignore this normalization problem here and just focus on the evolution of
to measure cosmic acceleration
q or
(which are fairly independent of
).
In the right panel of
Figure 5, we show (as shaded regions) the measurements for
given by combining Equation (
20) with Equation (
19) and Equation (
22) with Equation (
21). The measurements clearly favor models with large negative cosmic event acceleration
, which supports our interpretation of
as a friction term.
Comparing the left and right panels in
Figure 5, we see that both
q and
are roughly consistent with models with
(or
), in good concordance with
in the upper left panel of
Figure 5.
Even when the underlying model for
q and
is the same, note how the measured
q and
data have different tensions with the model predictions as a function of redshift. In particular, the radial BAO and SNIa datasets show inconsistencies among them for
q around
. This is a known tension (see, e.g., Figure 17 in [
28]). But note how this tension disappears when we use the corresponding estimates for
. Thus, data are more consistent with the
than with the
q description.
One would expect that a perfect realizations of the LCDM model in Equation (
4) would produce consistent results in both
q and
. But deviations from LCDM and systematic effects can produce tensions in data, especially if we use a parametrization like
q that refers to events that we never observe. The
q and
parametrizations of acceleration are more general than the particular LCDM model, and the fact that data prefer
is an important indication. Data live in the light-cone, which corresponds to
rather than
q. At
, the difference between a light-cone and space-like separations is very significant, and any discrepancies in the data or model will be more pronounced in
q modeling.
We conclude that the data show some tensions with LCDM predictions (as indicated by q) but confirm that cosmic expansion is clearly decelerating (as indicated by ) so that events are trapped inside an event horizon ().
7. Discussion and Conclusions
In our exploration, we have demonstrated that the commonly interpreted
term, thought to drive cosmic acceleration (as discussed in
Section 2), actually leads to a quicker cosmic deceleration of events compared to the influence of gravity alone (as explained in
Section 4). This relates to the nature of the event horizon (EH) (see
Section 5) that results from an expansion dominated by
. It suggests that
might not be a new form of dark or vacuum energy
([
13,
14,
15]) or a modification of gravity
but rather a boundary or surface term
in the corresponding action (see Equation (
1)).
The measured
in our cosmic expansion exhibits behavior analogous to the interior of a Schwarzschild black hole (BH), particularly under the assumption of nearly empty space beyond
(see [
16]). This analogy is illustrated in
Figure 2 and
Figure 4, which challenge the conventional interpretation of the FLRW metric. The FLRW model assumes that the background density
remains constant both inside and outside
, despite the absence of causal connections between these regions. How could this be possible?
The absence of causality suggests that the universe could be inhomogeneous outside
while remaining homogeneous within
, as indicated by cosmic expansion. This concept diverges from the traditional horizon problem as described by [
3,
29], in which the standard cosmological model envisions a homogeneous and isotropic universe fragmented into approximately
causally disconnected regions in the past [
30].
Inflationary theory, initially proposed by [
31] and further developed by [
30,
32,
33], addresses the problem of cosmological fragmentation. It posits a period of exponential expansion driven by a state of energy characterized by the ground state of a field that translates into an effective negative pressure (similar to
in Equation (
1)). This foundational phase allows scales to initially exit the Hubble horizon (as depicted in the yellow regions of
Figure 2 and
Figure 4) and to re-enter post-inflation. Despite its success, the origins of this inflationary period remain elusive, posing a significant mystery in theoretical physics. Furthermore, inflationary theory does not address the causality outside
(see
Figure 4), raising critical considerations for our understanding of the Universe’s expansion.
The event horizon
measured with
(i.e., Equation (
17), which is equivalent to the presence of
) also tells us that there is a finite mass
trapped within
. If we assume that the space outside
is relatively empty, such a finite mass
provides the explanation for the observed
and, therefore, for
: i.e.,
as resulting from a boundary term (see Equation (
1)). This black hole Universe (BHU) model provides a new and completely classical explanation for the cosmological constant
within GR. It explains why
is so small but not zero: because
is so large but not infinite. Yet it also raises new fundamental questions: If our local universe has fallen inside its own gravitational radius
, why is our universe now expanding and not collapsing?
The BHU interpretation, where the expansion happens inside
(and is therefore a local solution [
16]), opens the way to a new conjecture for the origin of cosmic expansion. Instead of emerging from a singular Big Bang (a global solution), it could result from the cold collapse of a large and low-density (local) cloud into a black hole. Such a collapse originates from a small initial over-density within a flat background, which can be modeled as a local FLRW closed curvature
solution, where the curvature radius
is the comoving radius of the initial cloud. Note that [
16] focused on the junction conditions for a flat
case (which is the right approximation for our expanding universe or a flat expanding perturbation). The corresponding case with curvature during collapse was solved in §12.5.1 of [
34], with the only difference with respect to the
being the replacement of the adimensional
by
. The total relativistic mass
of the collapsing cloud is also given by
, which relates to both
(as a boundary term:
) and the initial curvature radius
. As the collapse approaches the almost singular ground state, the curvature increases as
, which, together with positive background acceleration (from the degenerate ground state), enables the bounce to happen. After the bounce, the inflationary expansion erases the curvature term, while
will eventually dominate the bouncing expansion. This is illustrated in
Figure 6, which is reproduced from [
9]. Unlike traditional models that lead to a singularity, this model suggests that a Big-Bang-like explosion—termed the Big Bounce—prevents such an outcome. This Big Bounce could be driven by neutron degeneracy pressure, which occurs when densely packed neutron matter reaches a ground state governed by the Pauli exclusion principle. It could also be the result of a similar ground state happening at higher energies, like in standard cosmic inflation.
Mirroring the dynamics predicted by cosmic inflation, a ground state acts like a relativistic fluid with negative pressure in a closed cloud (
). The combination of these two ingredients (positive curvature and positive background acceleration) not only halts the gravitational collapse but also catalyzes a fast rebound (exponential expansion that erases the original curvature), initiating the expansive phase of a flat Big Bang. The expansion drives the system away from the ground state, but the system returns to regular radiation- and, later, matter-domination phases. This expansion is eventually stopped by another quasi-de Sitter phase, this time caused by the finite mass of the system
. Crucially, this quantum mechanism (Pauli exclusion principle) violates the strong energy condition
(but not the weak one
) in classical general relativity (GR) within a closed metric (
) and sidesteps the singularity GR theorems proposed by Hawking and Penrose (e.g., [
35]), thus presenting a novel solution to a pivotal issue in cosmological theory.
This idea is further validated by the observed large-scale cut-off in the scale-invariance spectrum of metric perturbations, as observed in the CMB sky (see [
36]). Such a cut-off is measured to be 66 degrees, which corresponds to the
radius at recombination (see dotted horizontal line in
Figure 4) projected in the CMB sky. Recent research ([
37]) has further revealed that several large-scale persistent CMB temperature anomalies originate from parity asymmetry
. A groundbreaking explanation posits that the microscopic laws of quantum physics adhere to
symmetry in a way that preserves causality and promotes curved spacetimes. Cosmic evolution disrupts the
symmetry, resulting in the observed
asymmetry. This idea was originally applied to inflationary quantum fluctuations but can be equally applied to the BHU Big Bounce picture explained above, as they are both defined by a period of quasi-de Sitter expansion.
Additionally, we can conjecture from this notion that the interior dynamics of any other BH (e.g., stellar, binary, or galactic) could also result from a similar BHU solution: a classical and non-singular FLRW expanding interior (that becomes asymptotically de Sitter, i.e., static in the rest frame). The mass (equivalent to ) boundary term in the BHU can then be interpreted as the actual physical mechanism that prevents anything from escaping the BH interior: i.e., it prevents the inside-out crossing of the BH event horizon , which asymptotically results from .
That the measured
term is fixed by the total mass
of our universe is in good agreement with the physical interpretation presented here that
, in the rest frame, corresponds to a friction (attractive) force that decelerates cosmic events. In
Appendix A, we elaborate on this idea and revisit the Newtonian limit to show that
corresponds to an additional (attractive) Hooke’s term to the inverse square gravitational law.