1. Introduction
The concordance model of cosmology,
Cold Dark Matter (ΛCDM), is thought to work well at explaining several observations of the Universe, including the observed late cosmic acceleration [
1,
2], the structure of the cosmic microwave background (CMB), the baryon acoustic oscillations (BAOs), among others [
3]. In recent years, there have been several issues at hand [
4], e.g., CMB anisotropy anomalies [
5], cosmic dipoles [
4], and the most notable issue: the Hubble constant tension, reaching a statistical discrepancy of ∼5
[
6,
7,
8,
9] between early and late measurements.
Direct measurements of the current Hubble constant
, through late-universe observations, have yielded a value of
km s
−1 Mpc
−1 [
8], while early-universe observations of the CMB give a lower result of
km s
−1 Mpc
−1 [
6], differing by around 8–10%. This discrepancy has positioned itself as an important crisis for the standard ΛCDM model. Consequently, several theories have been proposed to solve, alleviate or give insights into the Hubble tension issue. In the same line of thought, some proposals focus on the nature of dark energy, since if the cosmological constant Λ played a bigger role than previously thought in the early stages of the expansion of the Universe, it could alter the calculation for the early-universe Hubble constant [
10], as a consequence, this would bring it closer to the locally measured
value. On one hand, changing the equation-of-state of dark energy by making it time-dependent, has also been proposed as a possible explanation [
11]. Other theories include adding new interactions, such as interacting dark energy [
12,
13], interacting dark matter [
14] and interacting holographic dark energy [
15]; including sterile neutrinos [
16], axions [
17], or other theorized particles to increase the relativistic degrees of freedom; unifying dark energy and dark matter as a single fluid [
18]; introducing cosmological piecewise functions [
19]; using a statistical lookback time approach [
20], considering the black hole shadows as standard rulers calibrators [
21], among many others [
22,
23,
24]. See, for example, [
25], for a summary and comparison of several proposals that involve interesting solutions to the
tension.
On the other hand, alternative gravity theories also attempt to solve the Hubble tension and other problems within standard cosmology. These theories change General Relativity’s (GR) Einstein-Hilbert action to introduce different dynamics in the resulting equations of motion, which can explain different phases of expansion of the universe and could explain the Hubble tension [
26,
27]. Allowing higher order equations of motion, there are theories such as
theories, which introduce a function of the Ricci scalar,
, to act as the action. In the context of Teleparallel Gravity [
28], it can be considered
theories, which instead use functions of the torsion scalar. Another way to extend GR is by allowing the Lagrangian density to depend not only on the metric but also on another degree of freedom such as a scalar field
, these are the so-called scalar-tensor theories. If we maintain the requirement of working in four dimensions, and having second-order equations of motion—to avoid Ostrogradsky instabilities [
26]—the most general of these theories is Horndeski theory [
29,
30].
In this work, we focus on a special case of Horndeski’s theory first proposed in [
29] called the
Fab Four theory. As a first approach, in [
29] a
self-tuning filter was applied on the original Horndeski action that partially de-couples the curvature from the cosmological constant and reduces the action to only four base Lagrangians, each one a function dependent on combinations of the Riemann tensor and an arbitrary function of the scalar field
. The self-tuning filter consisted in requiring the following [
29]:
The theory admits a patch of Minkowski vacuum for any value of the cosmological constant Λ, and for instantaneous changes in its value;
The theory allows non-trivial cosmologies by introducing additional dynamical terms to the standard one.
According to these two requirements, in [
31], an initial outlook was provided on the types of cosmology and its dynamics derived from the Fab-Four theory. Although it is a fairly recent proposal, the Fab-Four theory has been explored extensively and in different contexts. Specifically, there have been several studies involving the
term in isolation, since it is the one that mainly drives the self-tuning mechanism [
32] and it contains a derivative coupling to the scalar field, e.g., in [
33] the
term (see
Section 2) is studied in detail, specifically its generalization as the
term in Horndeski’s theory which adds a dependency on the field kinetic term
, finding that the resulting
Fab 5 theory can be ghost-free, or stable, or self-tuning, but not all these characteristics at the same time.
Fab 5 cosmology is studied in more detail in [
34,
35,
36]. After the extension of the original Horndeski Lagrangian [
37], adding two more terms in Horndeski’s action, it was shown in [
38] that these new terms can be included in a generalized Fab-Four theory where the arbitrary functions of the field
and
become arbitrary functions of the field and its kinetic term
and
. Additionally, it was shown that these new terms also present self-tuning for screening the cosmological constant. In [
39], it was argued that the measurement of the speed of gravitational waves (GWs) strongly constrained
John-like terms in Horndeski’s action, along with several other scalar-tensor theories and dark energy (DE) models. However, in [
32], it is shown that working in an Effective Field Theory (EFT) framework, the resulting effective action is not constrained as such, and so it still admits
John-like terms, therefore allowing a self-tuning mechanism in the EFT regime and the possibility to be consistent with the observational data [
39,
40]. In order to explore GW constraints, see [
41]. Furthermore, as any alternative theory of gravity could, the Fab Four proposal could be excluded if we demand that the speed of gravitational waves
is exactly the same as the speed of light
c along the evolution of the Universe. Therefore, if we consider small deviations from
c, then several models can become consistent with the observable data (and their constrictions). New self-tuning models known as
well-tempered have been developed and deal with the GW constraints by construction [
42], and also prevent the scalar field from screening not only the cosmological constant but also all other matter sources in the system. A generalized class of
well-tempered models and their self-tuning mechanisms has been discussed in [
43]. In a quantum gravity context, the study of the quantization of the Fab-Four
John term results in a Hamiltonian with fractional powers of the momenta, leading to the use of a conformal fractional derivative approach to get a wide range of solutions for the scale factor and the scalar field [
44].
In this work, we follow closely the derivations performed for the Fab-Four cosmology [
31] to obtain the corresponding equations of motion for a Friedmann-Lemaître-Robertson-Walker (FLRW) background and to apply a dynamical systems approach to solving the dynamical equations. Under this scheme and considering all four components of Fab-Four theory, we also extend this analysis by using a cosmographic approach, which will be useful to obtain kinematic quantities, such as
and
. Finally, with these quantities at hand, we will be capable of setting general requirements on the free parameters to relax the Hubble tension and give a likelihood landscape for future works and analyses of this theory. The core idea behind this process is to facilitate the search for a stable cosmological solution through dynamical systems. The advantage of this kind of solution is the possibility to explore those that fulfill the cosmological constraint equation of interest for the analysis scheme. Once this constraint is satisfied, the analytical solution at hand can be tested with observations to establish a cutoff for the cosmological parameters of interest, namely
and
. Although there may exist more analytical solutions for the set of equations, our cosmographic approach up to second order allows us to narrow down the options and select the ones that address the late-time cosmic acceleration.
This paper is divided as follows: In
Section 2, we present a brief introduction of the results discussed in [
31], where the full equations of motion of the system are given for the Fab-Four theory in a FLRW background. Also, we describe the dynamical system approach for this theory to find the dynamical parameters and the corresponding fixed points. We rewrite the equations in the zero-curvature fixed point of interest and solve them to find viable cosmological scenarios. In
Section 3, we construct the cosmographic series for the Fab-Four theory and find the expressions for the Hubble constant
, and the deceleration parameter
for each cosmological scenario described. Also, we will restrict our analysis to a couple of cases due to the solutions found in the dynamical system approach. In
Section 4, we discuss how our solutions from the Fab-Four theory are able to relax the
tension by following a set of requirements on the free parameters of the theory. Finally, we present our conclusions and considerations for future work in
Section 5.
2. Fab-Four Cosmology as a Dynamical System
The Fab-Four model comes from a restriction set upon Horndeski’s theory, which simplifies it greatly so that the Lagrangian density is now just the sum of four base Lagrangians [
29]:
where each Lagrangian is defined as a combination of derivatives of the scalar field
, combinations of the metric tensor, and an arbitrary function of the field
, with
, namely:
where
is the Einstein tensor,
is the double dual of the Riemann tensor,
R is the Ricci scalar and
is the gauss-Bonnet combination [
31]. Given this notation and considering an FLRW metric whose line element is given by
where
is the scale factor, which denotes the relative size of the co-moving, constant-time hypersurfaces of space-time,
is the element of solid angle and
k is the spatial curvature constant. Here,
gives a flat space, while
(
) corresponds to positive (negative) curvature. We can express the Fab-Four Lagrangians as follows [
31]:
Here,
denotes the derivative of
with respect to the field,
, and
, with
the 4-velocity. Following [
31], we find that in the presence of standard matter and a cosmological constant Λ, the Hamiltonian density is given by the following:
where, for convenience, we can express each Hamiltonian in conformal time by considering
, where
denotes differentiation with respect to
N, such that
:
The equations for the scalar field
can be obtained through varying
which gives the following expressions [
31]:
with
. Notice that this function goes to zero as the latter equations approach a self-tuning solution corresponding to a Ricci flat scenario and the full system of equations that arises is coupled with second-order differential equations, which are extremely difficult to solve analytically. Therefore, a viable path to reach a set of possible solutions is through a dynamical systems approach [
31]. In order to implement this procedure, we consider the following dynamical variables:
notice that if
const. we obtain that
where a matter-like expansion would arise from
, radiation-like expansion from
, and an inflationary expansion coming from
. To build the dynamical system according to the latter variables, we define the following expressions [
31]:
where
is an arbitrary constant that provides freedom in re-defining the scalar field while maintaining equations of motion that do not explicitly depend on
H. This freedom will later be fixed to a particular value in
Section 3.
The evolution equations for these variables are directly obtained by differentiating Equations (21)–(25) with respect to
N. Finally, we obtain:
where
and
denotes the exponents of
H in Equation (23). Knowing how these variables evolve, we can easily find their fixed points. Considering constant
allows for potentials in a power law form. For the
and the
, these fixed points are those where
In this fixed point, combining Equations (30) and (31) yields:
which allows us to directly find the exponents for all four potentials, and they become
Furthermore, the fixed point for
corresponds to
, i.e., fixed points solutions exist for
and
. On this line of thought, we are interested in the case where
, which refers to vanishing spatial curvature, as it is consistent with observations [
6]. In this selected fixed point, by also using Equations (30) and (31), the full scalar field equation Equation (20) can be written as follows:
Similarly, the Hamiltonian constraint Equation (11) can be written as
Following the arguments in [
31], because the self-tuning filter demands that the solutions hold for any value of the cosmological constant, the parameters in the potentials (
,
) should be all independent of
. This assumption constrains the possible values of the terms present in Equations (37) and (38) while satisfying the following:
We can distinguish several values for h that fulfill these requirements, which are associated with different behaviors of a, due to . Two cases of interest in this article are arbitrary for and radiation for .
Under the arbitrary case, we focus on three cases: Firstly, the case of
corresponds to an evolution where
, indicating an accelerated expansion. Secondly, the case of
exhibits an evolution with
, resembling matter-like expansion. Lastly, we consider the radiation case, which is applicable for studying the standard cosmological model. This case encompasses an inflationary period followed by a radiation-dominated epoch, then a matter-like expansion facilitating structure formation, and finally, a late-time inflationary expansion that aligns with our current observations [
3]. The applicability of such cases is due to the fact that we are looking for a viable cosmological scheme that can reproduce the current cosmic dynamics. Of course, in this case, we are considering solely a particular scalar field dynamic minimally coupled to gravity; so, the possibility that this scalar field can imitate dark matter and the effects of dark energy could set up a
toy contemporary universe model.
3. Fab Four Cosmography
In what follows, we proceed with the standard cosmographic methodology [
45,
46,
47] to derive the kinematic parameters, which will be evaluated at current cosmic times. The advantage of considering this approach is the geometry construction behind it. By assuming the structure of the Cosmological Principle, we can factorize the scale factor around the present time. However, we should take into account that we must relax the issues of such an approach: (i) the convergence of the series, which can be attended by using observables in order to restrain the order of it. (ii) the cosmic systematics that could be propagated on the cosmographic parameters. In such a case, we require a control of the numerical fit. Notice that for this analysis, we must select a fixed stable point and two possible cosmological scenarios.
Staying in the fixed point
, we can start by constructing cosmographic series based on the Taylor series expansion context of the scale factor
a, defined as the following:
where
q is defined as
and
is the current deceleration parameter. To obtain
, we first choose an appropriate scenario as described in
Section 2, substitute it into the Hamiltonian constraint Equation (11), and solve the system for the Hubble parameter
H. Furthermore, with
H we can find
q using Equation (43). Finally, we can then evaluate both expressions at the current time to obtain
and
, respectively.
Specifically, we chose to build the cosmographical series because it is a purely kinematic expression where other parameters such as matter density or curvature do not play any role. Therefore, we can focus solely on the kinematic terms without specifying other parameters of the model.
At this point, we consider the radiation-dominant behaviour case (for the early-universe calculation of
), and the arbitrary case (for the late-universe calculation of
). On one hand, for the arbitrary case, we have only two potential functions different from zero (as can be seen in Equations (39)–(41) by letting
h be arbitrary):
where the subindex
arb denotes arbitrary. Therefore, the only Hamiltonian parts considered would be
and
; thus, we can write Equation (11) as follows:
From the latter, we can solve for
H:
which, evaluating at current times, gives us the following:
where we have discarded the negative signs of the fourth roots to ensure a real, positive value of
. Substituting Equation (48) in Equation (43), and evaluating at the current time gives us the following:
On the other hand, for the radiation scenario, the potentials to be considered are as follows:
which come from considering
in Equations (39)–(41); with the corresponding Hamiltonian parts
,
and
. The full Hamiltonian constraint for this case is:
where solving for
H we obtain:
Here, we have defined the quantity:
which will appear frequently in our results. From Equation (55) and using Equation (56) we can obtain values for
and
for the radiation scenario:
Now, in the calculated values for
and
in Equations (49), (50) and (57)–(59) there are several free parameters, namely:
,
,
,
,
h and the scalar field and its derivative at current time
and
. We can choose an arbitrary value of
to fix the freedom in re-defining the scalar field. Following [
31], we can select
such that
. Alternatively, we can choose values of
to ensure we can obtain power law potentials of the form
, which gives an inflationary behaviour. Only the cases where
and
, give potentials of this kind, while
is preferred, as it also fulfills
when
(inflation).
3.1. Case for α = 1
For this case, we consider such
value in Equations (49) and (50), which yields:
and
Now, choosing a value for
h determines how the Hubble parameter depends on
a according to Equation (21), and thus specifies the type of expansion. We can study the cases for
and
, which correspond to matter-dominated expansion and inflation, respectively. Substituting these values of
h in Equations (59) and (60) gives:
which correspond to
, and
which corresponds to the inflationary case
.
Now, the case for radiation-like expansion is given by substituting
in
as in Equations (57)–(59):
3.2. Case for α = 2
Another way to obtain potentials of the form
is with
. Substituting this value on Equations (49), (50), and (57)–(59), and following a similar procedure as in the latter case, we obtain the following:
and for the matter case,
The inflation case gives us the following:
and for the radiation case,
For both values of
, we notice similarities in all forms of
, which are all proportional to
. Additionally, we notice similarities in the expressions for
, where all of them can be written in the same functional form:
where
A,
B and
C are, in general, functions of the parameters
,
,
,
and the constants
,
and
. Considering standard GR with a cosmological constant and a scalar field as the source in an FLRW background, the deceleration parameter has the same functional form.
4. Fab-Four on the H0 Tension
According to the latter results, we can now discuss the possible requirements to relax the value of in the cosmological scenarios discussed.
With the expressions obtained for
and
we can now calculate the predicted discrepancy between both values. The most recent measurements [
6,
8] suggest a discrepancy of around ∼10% between the values for the Hubble constant in the early universe and those in the late universe. As such, we expect that
For
, we can use Equations (63)–(65) to get:
On the other hand, for
, using Equations (71)–(73), we obtain the following
Now, before proceeding, several requirements arise in the process of obtaining expressions for
,
, and
to ensure that the solutions are real. These are summarized in
Table 1.
Taking into account the restrictions given by
Table 1, it is possible to get values for
, by considering the requirements on the five free parameters that appear in the expressions:
and
. We also have another two parameters,
and
, which only appear in the expressions for
.
However, the behaviour of
comes into play. The Fab-Four model was introduced to handle the cosmological constant fine-tuning problem [
29,
31]. As such, instead of the cosmological constant being constant, it is assumed that it is a piecewise constant function, with several phase transitions that change its value. Therefore,
is discontinuous at these phase transitions. This means that the
r.h.s of Equation (11) is discontinuous, but the
l.h.s of this equation is composed of the Hamiltonian parts, which depend only on the functions
a,
, and
. Since
a is continuous, and
as well—by construction—it follows that
must be discontinuous to maintain equality at the phase transitions. Consequently, both
and
contain delta-functions localized at the transition times. These same arguments are used to initially restrict Horndeski’s full theory of gravity [
29].
At these phase-transitions, the self-tuning nature of the Fab-Four theory is active to still allow a standard cosmology to develop [
29]. However, by construction, the self-tuning filter can only work when the curvature is non-zero. Since the latest observations assert a negligible curvature [
6], we assume that we are not presently in a phase transition. For this work, we are only taking the present-day values of these parameters (
,
), which we will assume to be zero since we are not currently in a phase-transition, where there would be a delta-function.
Therefore, substituting
and
, from Equations (64)–(66) and (72)–(74), we obtain:
for each
case, respectively. Now, given a value of
through the latest observations [
8] (
), we can either force the deceleration parameters (evaluated at the current time) of both the early and late universe to be the same, or we can leave them to be different.
In this work, we assume for simplicity, and we leave the other case for future work. Given this restriction, with the expressions for Equations (79) and (81), we will get a fixed value for (depending on the choice of ). Then, using the expressions for Equations (80) and (82), we can constrain the value of one more parameter, say, for example, .
Finally, we are only left with two free parameters:
and either
or
. These are enough to get
, and reducing the number of free parameters is important to enhance the likelihood of this model in Bayesian statistic comparisons [
48].
5. Conclusions
In this work, we consider the possibility of tackling the Hubble tension using the Fab-Four theory. This scalar-tensor theory of gravity applies a self-tuning filter on the more general Horndeski theory, to solve the cosmological constant fine-tuning problem. This leaves a much more manageable theory, consisting only of four base Lagrangians, given by Equations (7)–(10). We follow the procedure discussed in [
31] to obtain a simpler form of the theory, using a dynamical system approach to find the fixed points of the theory, and we concentrate on the
fixed-point, which denotes zero curvature (as supported by observations [
6]). In this fixed point, the equation of motion of the scalar field Equation (37) and the Hamiltonian restriction Equation (38) of the theory are such that we have certain allowed types of cosmological scenarios, among which are the inflationary expansion, radiation-like expansion, and a matter-like phase. This shows the utility of the Fab-Four theory, as
turning on-and-off certain potentials lead naturally to different types of cosmological behaviors, and the full dynamics of the theory allow for a smooth, continuous shift between these types of cosmological evolution.
Furthermore, we construct the Fab-Four cosmography in the zero-curvature fixed point. In particular, we derive analytical expressions for the current Hubble constant
, and the deceleration parameter
, for the matter, inflation, and radiation scenarios, as well as for different choices of the arbitrary constant
. To do this, we take into account the forms of the potentials, substitute them into the Hamiltonian constraint, which results in quartic equations of the Hubble parameter, and then solve for
H. These equations yield four solutions:
(i) in the arbitrary case, we discard two complex roots and a negative one;
(ii) in the radiation case, we can only discard two negative roots, leaving both signs of the expression for
. Examining
Table 1, we also notice that the requirements do not restrict the choice of signs in the expressions for
and
, and both signs can give
, which is why we have left them written explicitly in all expressions. Afterward, we differentiate the results for
H according to Equation (43) to obtain expressions for the deceleration parameter, which all ultimately take the functional form of Equation (75). Using these expressions evaluated at current time, we can compare the Hubble constant in the late universe with the Hubble constant in the early universe (radiation-like expansion).
Finally, we show that it is possible to get the desired discrepancy between the early-universe and the late-universe Hubble constants, of around ∼10%, showing that the Hubble constant could, in reality, be dynamical. Therefore, we have succeeded in our objective to show that it is indeed possible to get a disagreement between the observed Hubble constant if measured through local means, and the one obtained from observations of the CMB, due to the different dynamics put forward by the Fab-Four theory.
The final results contain five arbitrary parameters which are constrained by the requirements in
Table 1. However, using the expressions for the deceleration parameter, we can reduce the number of free parameters in these results to only two constants. In future work, we will propose performing an analysis with observational data to obtain confidence limits on the free variables of the theory. For this, extracting other observables from the Fab Four theory, e.g., expressions for the luminosity distance, and even perturbation theory results such as the angular power spectrum from the CMB would be convenient. Afterward, a possible course of action would be to make a Bayesian statistics comparison of this model’s prediction against the standard ΛCDM to test its plausibility, as well as to further constrain the possible values of the free parameters to match observations [
48]. Furthermore, the possibility to explore beyond Horndeski theories [
26] could be a path in order to tackle the cosmological tensions at early cosmic times.