Propagation of Tensor Perturbation in Horndeski-like Gravity

Abstract

Scalar–tensor theories have shown promise in many sectors of cosmology. However, recent constraints from the speed of gravitational waves have put severe limits on the breadth of models such classes of theories can realize. In this work, we explore the possibility of a Horndeski Lagrangian that is equipped with two dilaton fields. The evolution of a two-dilaton coupled cosmology is not well known in the literature. We explore the tensor perturbations in order to assess the behavior of the model against the speed of the gravitational wave constraint. Our main result is that this model exhibits a class of cosmological theories that is consistent with this observational constraint.

1. Introduction

Over a century ago, Einstein’s General Theory of Relativity provided a revolution in our understanding of gravity and spacetime, laying the foundation for modern cosmology [1,2,3]. Since its inception, the theory has led to numerous groundbreaking discoveries, including the development of the standard cosmological model, which emerged in the latter half of the 20th century. This model is rooted in two pivotal observations: the recession of galaxies, first identified by Edwin Hubble [4], and the discovery of the Cosmic Microwave Background (CMB) radiation by Penzias and Wilson [5]. These observations provided compelling evidence for the Big Bang paradigm and established a framework for understanding the large-scale structure and evolution of the universe.

Recent advances in observational cosmology, particularly in the study of the CMB and gravitational waves (GWs), have opened new avenues for testing these theoretical frameworks. For instance, the stochastic background of primordial gravitational radiation predicted by cosmology has been shown to be consistent with Planck data [6,7] and may fall within the sensitivity range of current and future GW detectors such as LIGO, Virgo, and LISA [8]. Notably, pre-Big Bang models, as described in [9], occupy a wide region of parameter space that allows for the generation of a tensor perturbation spectrum compatible with these observations.

In addition to Einstein’s theory, alternative gravitational frameworks such as Lovelock gravity, Horndeski gravity, and Gauss-Bonnet gravity have found applications in cosmology [10,11,12,13,14,15,16,17,18]. These models are particularly attractive for inflationary scenarios, as they naturally incorporate scalar fields arising from the gravitational sector. Recent precision measurements from the Dark Energy Survey (DES) and Planck satellite observations [6] have refined our understanding of this expansion, revealing subtle tensions in the Hubble constant ($H_0$) when derived from early versus late universe probes. These discrepancies potentially signal physics beyond the standard $\mathsf{\Lambda }$CDM cosmological model, possibly involving modified dark energy equations of state [19,20] or additional relativistic species in the early universe. In this work, we propose a novel cosmological background within the framework of Horndeski-like gravity, focusing on the direct inflationary production of gravitational radiation. Our approach begins with the amplification of tensor metric fluctuations, leading to the formation of a cosmic background of relic gravitons. This setup provides a promising avenue for exploring the interplay between modified gravity theories and the observational signatures of early-universe physics.

Alternative gravitational frameworks, such as Lovelock gravity, Horndeski gravity, and Gauss-Bonnet gravity [10,11,12,13,14,15,16,17,18], have emerged as compelling extensions to Einstein’s general relativity, offering new insights in cosmology. These theories are particularly appealing in the context of early-universe physics, as they naturally accommodate scalar fields that arise from the gravitational sector, providing a robust theoretical foundation for inflationary models. Unlike standard inflationary scenarios, where scalar fields are often introduced ad hoc, these frameworks integrate them as intrinsic components of the modified gravitational dynamics.

In this work, we explore a novel cosmological background within the framework of Horndeski-like gravity, focusing on its implications for the direct inflationary production of gravitational waves. Specifically, we investigate the amplification of tensor metric fluctuations during inflation, which leads to the generation of a stochastic background of relic gravitons. This mechanism not only enriches our understanding of the interplay between modified gravity and early-universe physics but also offers a testable prediction for future gravitational wave observatories, such as LISA and LIGO [8]. By leveraging the unique features of Horndeski-like gravity, our approach provides a promising avenue for probing the observational signatures of alternative gravitational theories and their role in shaping the dynamics of the primordial universe.

Recently, applications of the Horndeski gravity in cosmology [12,13,14,15,16] has been achieved through the following Lagrangian

(1)$\begin{array}{c}\hfill {\mathcal{L}}_{\mathrm{H}}=(R-2\mathsf{\Lambda })-{\displaystyle \frac{1}{2}}(\mathsf{\alpha }{g}_{\mu \nu }-\gamma G_{\mu \nu }){\nabla }^{\mu }\varphi {\nabla }^{\nu }\varphi ,\end{array}$

where R, $G_{\mu \nu }$ and $\mathsf{\Lambda }$ are the scalar curvature, the Einstein tensor, and the cosmological constant respectively, and where $\varphi =\varphi (r)$ is a scalar field. Also, $\mathsf{\alpha }$ and $\gamma$ are coupling constants. Using the Lagrangian in Ref. (1) equipped with two dilatons, we proposed the scenario of the propagation of tensor perturbations considering the action

(2)$\begin{array}{c}\hfill S=-{\displaystyle \frac{1}{16\pi G}}\int d^{4}x\sqrt{\left|g\right|}{e}^{-\varphi }({\mathcal{L}}_{\mathrm{H}}-{\displaystyle \frac{1}{2}}{\nabla }_{\mu }\chi {\nabla }^{\mu }\chi +V(\varphi ,\chi )).\end{array}$

To reproduce the standard Einstein-Horndeski-like equations, we need to add the boundary term $S_{\mathsf{\Sigma }}$:

(3)$\begin{array}{c}{S}_{\mathsf{\Sigma }}=-{\displaystyle \frac{1}{8\pi G}}\int d^{3}x\sqrt{\left|g\right|}{e}^{-\varphi }({\mathcal{L}}_{\mathrm{H},\mathsf{\Sigma }}+{\mathcal{L}}_{ct}),\hfill \\ {\mathcal{L}}_{\mathrm{H},\mathsf{\Sigma }}=(K-\mathsf{\Sigma })-{\displaystyle \frac{\gamma }{4}}{\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi K^{\mu \nu }\hfill \end{array}$

(4)$\begin{array}{ccc}& & -{\displaystyle \frac{\gamma }{4}}({\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi n^{\mu }{n}^{\nu }-{(\nabla \varphi )}^2)K,\hfill \\ & & {\mathcal{L}}_{ct}=c_0+c_{1}R+c_2{R}^{ij}{R}_{ij}\hfill \end{array}$

(5)$\begin{array}{ccc}& & +c_3{R}^2+b_1{({\mathsf{\partial }}_i\varphi {\mathsf{\partial }}^i\varphi )}^2+\cdots .\hfill \end{array}$

Here $K_{\mu \nu }=g_{\mu }^{\phantom{\mu }\beta }{\nabla }_{\beta }{n}_{\nu }$ is the extrinsic curvature, $K=g^{\mu \nu }{K}_{\mu \nu }$ is the trace of the extrinsic curvature, $n^{\mu }$ is an outward pointing unit normal vector to the boundary, $\mathsf{\Sigma }$ is the boundary tension. ${\mathcal{L}}_{ct}$ is the boundary counterterms.

In this work, we will address the decelerated expansion (i.e., from inflationary to standard evolution), which causes an amplification of the quantum fluctuations of the various background fields (the background fields in our case will be the gravity fields Horndeski) and can produce a large number of various types of radiation [9,10]. The spectral properties of this radiation are strongly correlated with the kinematics of the inflationary phase (or more precisely in our work, with the Horndeski kinetic term of gravity) [9]. Direct (or indirect) observations of such a primordial component of cosmic radiation could thus provide us with important information about inflationary dynamics, testing the predictions of the various cosmological scenarios.

2. The Overarching Theory

In this section, we develop a framework for inflationary kinematics within the context of Horndeski gravity, a generalized scalar-tensor theory that extends the standard cosmological paradigm (see, e.g., [18]). Inflationary dynamics in this framework can exhibit deviations from conventional field-theoretic models, particularly in scenarios inspired by cosmology. These deviations are critical for identifying unique signatures inspired by models and understanding their implications for early-universe physics, including the generation of primordial perturbations and their imprints on the large-scale structure (see, e.g., [21,22]). The effective action for the gravitational sector in our Horndeski-like framework is given by:

(6)$\begin{array}{c}\hfill S_{\mathrm{total}}=S+S_{\mathsf{\Sigma }}+S_m,\end{array}$

where $S_m$ represents the matter action. By varying this action (6), we derive the Einstein-Horndeski equations of motion:

(7)$\begin{array}{c}\hfill G_{\mu \nu }+\mathsf{\Lambda }{g}_{\mu \nu }=2T_{\mu \nu },\end{array}$

where the energy-momentum tensor $T_{\mu \nu }$ includes contributions from scalar fields and matter sources. Specifically, $T_{\mu \nu }$ is expressed as:

(8)$\begin{array}{c}\hfill T_{\mu \nu }=\mathsf{\alpha }{T}_{\mu \nu }^{(1)}+\gamma T_{\mu \nu }^{(2)}+T_{\mu \nu }^{(3)}-g_{\mu \nu }V(\varphi ,\chi ),\end{array}$

with $T_{\mu \nu }^{(1)}$, $T_{\mu \nu }^{(2)}$, and $T_{\mu \nu }^{(3)}$ given by:

$\begin{array}{ccc}& & T_{\mu \nu }^{(1)}={\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi -{\displaystyle \frac{1}{2}}{g}_{\mu \nu }{\nabla }_{\lambda }\varphi {\nabla }^{\lambda }\varphi +\hfill \\ & & {\nabla }_{\mu }\chi {\nabla }_{\nu }\chi -{\displaystyle \frac{1}{2}}{g}_{\mu \nu }{\nabla }_{\lambda }\chi {\nabla }^{\lambda }\chi ,\hfill \\ & & T_{\mu \nu }^{(2)}={\displaystyle \frac{1}{2}}{\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi R-2{\nabla }_{\lambda }\varphi {\nabla }_{(\mu }\varphi R_{\nu )}^{\lambda }-{\nabla }^{\lambda }\varphi {\nabla }^{\rho }\varphi R_{\mu \lambda \nu \rho }\hfill \\ & & -g_{\mu \nu }\left[-{\displaystyle \frac{1}{2}}({\nabla }^{\lambda }{\nabla }^{\rho }\varphi )({\nabla }_{\lambda }{\nabla }_{\rho }\varphi )+{\displaystyle \frac{1}{2}}{(\square \varphi )}^2-({\nabla }_{\lambda }\varphi {\nabla }_{\rho }\varphi )R^{\lambda \rho }\right]\hfill \\ & & -({\nabla }_{\mu }{\nabla }^{\lambda }\varphi )({\nabla }_{\nu }{\nabla }_{\lambda }\varphi )+({\nabla }_{\mu }{\nabla }_{\nu }\varphi )\square \varphi +{\displaystyle \frac{1}{2}}{G}_{\mu \nu }{(\nabla \varphi )}^2,\hfill \\ & & T_{\mu \nu }^{(3)}=e^{\varphi }{T}_{\mu \nu }^{(m)}.\hfill \end{array}$

Describing the kinematic and geometric contributions of the scalar fields $\varphi$ and $\chi$. These terms encapsulate the non-trivial interactions between the scalar fields and the spacetime curvature, as well as the coupling to matter sources. For instance, the scalar charge density $\mathsf{\sigma }$ of the sources is defined via the variation of the matter action with respect to $\varphi$:

(9)$\begin{array}{c}\hfill {\delta }_{\varphi }{S}_m=-{\displaystyle \frac{1}{2}}\int d^{4}x\sqrt{\left|g\right|}\mathsf{\sigma }\delta \varphi .\end{array}$

This formalism allows us to explore potential modifications to the standard cosmological model, particularly in the context of inflationary dynamics and their observational consequences. By connecting the theoretical predictions of Horndeski gravity to observable phenomena, such as the spectral index of primordial perturbations Figure 1, we aim to provide a comprehensive framework for testing these models against current and future cosmological data (e.g., [7]).

In order to study our configuration presented in the diagram Figure 1 at greater depth, we start with the flat Friedmann–Lemaître–Robertson–Walker (FLRW) metric of the form

(10)$d{s}^2=-d{t}^2+a^2(t){\delta }_{ij}d{x}^{i}d{x}^j,$

where $a(t)$ is the scale factor, and where the scalar field depends on cosmic time only. With this in hand, we find that the Friedmann equations using the Einstein-Horndeski equations of motion (7) in the form

(11)$H^2={\displaystyle \frac{4\kappa \mathsf{\Lambda }+\mathsf{\alpha }{\dot{\varphi }}^2(t)+{\dot{\chi }}^2(t)-2\mathsf{\sigma }exp(\varphi )+2V(\varphi ,\chi )}{3(4\kappa -3\gamma {\dot{\varphi }}^2(t))}},$

and

(12)$\begin{array}{ccc}& & (-4\kappa +\gamma {\dot{\varphi }}^2(t))H^2+4\gamma H\dot{\varphi }(t)\ddot{\varphi }(t)+2{\displaystyle \frac{\ddot{a}(t)}{a(t)}}(-4\kappa +\gamma {\dot{\varphi }}^2(t))+\hfill \\ & & (4\kappa \mathsf{\Lambda }-2\mathsf{\sigma }exp(\varphi )+2V(\varphi ,\chi )-{\dot{\chi }}^2(t)-\mathsf{\alpha }{\dot{\varphi }}^2(t))=0.\hfill \end{array}$

The equations of motion of the scalar field given an FLRW background are written as

(13)$\begin{array}{c}\hfill \ddot{\varphi }+3H\dot{\varphi }+{\displaystyle \frac{6\gamma \dot{\varphi }H\dot{H}}{\mathsf{\alpha }+3\gamma H^2}}+{\displaystyle \frac{1}{\mathsf{\alpha }+3\gamma H^2}}{\displaystyle \frac{dV(\varphi ,\chi )}{d\varphi }}=0,\end{array}$

(14)$\begin{array}{c}\ddot{\chi }(u)+3H\dot{\chi }+{\displaystyle \frac{dV(\varphi ,\chi )}{d\chi }}=0,\hfill \end{array}$

where we recovery the usual results of Refs. [16,23] for $V(\varphi ,\chi )=constant$ and $\chi =0$. The motivation to consider these two scalar fields, in particular the field $\chi$ is to find analytical solutions using the first-order formalism approach. With this in hand, we can reduce the equations of motion that are of second order, that is

(15)$\begin{array}{c}\hfill \dot{\varphi }=-W_{\varphi },\\ \hfill \dot{\chi }=-W_{\chi },\\ \hfill H=W.\end{array}$

Employing $H=W$, we can explicitly demonstrate that Hubble’s parameter evolves over time as a function of $W[\varphi (t)]$. This formulation not only provides a deeper understanding of the time dependence of Hubble’s parameter but also offers a more systematic framework for analyzing its behavior. Thus, firstly, combining Equation (11) with Equation (12) for simplicity, we find

(16)$W{W}_{\varphi \varphi }+{\displaystyle \frac{W_{\varphi }^2}{2}}+{\displaystyle \frac{4}{3}}{W}^2-\beta -2{\displaystyle \frac{W_{\chi }^2}{\gamma W_{\varphi }^2}}=0,$

and now combining the equations Equation (11) with Equation (13), we have

(17)$W{W}_{\varphi \varphi }+{\displaystyle \frac{W_{\varphi }^2}{2}}+{\displaystyle \frac{4}{3}}{W}^2-\beta -{\displaystyle \frac{3}{4\gamma }}{\displaystyle \frac{W_{\chi }{W}_{\chi \varphi }}{W{W}_{\varphi }}}=0,$

where $\beta =(4\kappa -\mathsf{\alpha })/\gamma$. However, we can see that for consistency when we compare the Equations (16) and (17), we get the following constraint on the superpotential

(18)${\displaystyle \frac{3}{8}}{\displaystyle \frac{W_{\chi \varphi }}{W}}={\displaystyle \frac{W_{\chi }}{W_{\varphi }}}.$

Following [13], we can choosing the simplest superpotential $W(\varphi ,\chi )=e^{\sqrt{{\displaystyle \frac{8}{3}}}(\varphi +\chi )}$ where this solution is satisfy by the Equation (15) which is consistent with Equation (18) given by

(19)$\begin{array}{ccc}& & \varphi (t)={\displaystyle \frac{3}{8}}\ln (t),\hfill \end{array}$

(20)$\begin{array}{ccc}& & \chi (t)={\displaystyle \frac{3}{8}}\ln (t),\hfill \end{array}$

(21)$\begin{array}{ccc}& & H(t)={\displaystyle \frac{9}{16}}{\displaystyle \frac{1}{t}}.\hfill \end{array}$

(22)$\begin{array}{ccc}& & a(t)=a_0{t}^{9/16},\hfill \end{array}$

which is the general solution of the system. Such solution can be found thanks to the field $\chi$, which provides an analytical system, the other hand, if $\chi =0$ with $\varphi \ne 0$, according to Equation (18) $W(\varphi ,\chi )=$ constant. Of course in the regime $\chi \to 0$, we recover some description of cosmological inflation in Einstein-Horndeski gravity, but performed by numerical computation [14,16,23]; the idea is to provide an analytical description of the Horndeski gravity [17]. Besides, the advantage of our model is that we do not choose the potential, this emerges naturally from first-order formalism with the scalar field $\chi$.

3. Propagation of Tensor Perturbations

We now analyze the tensor perturbations, which encapsulate the spectral properties of radiation in the context of cosmological evolution. Following the methodologies outlined in Refs. [12,13,14,24], we employ the TT (transverse and traceless) gauge, where ${\mathsf{\partial }}^{\mu }{h}_{\mu \nu }=0$ and $h_{\mu }^{\mu }=0$, with the scalar field perturbation $\delta \varphi$ set to zero. Under these conditions, the equations governing tensor fluctuations can be derived. These tensor modes are directly linked to observable phenomena in the recent universe, such as the polarization patterns in the Cosmic Microwave Background (CMB) and the gravitational wave background. By studying these perturbations, we gain insights into the early universe’s dynamics and the imprints left on the large-scale structure of the cosmos.

Tensor perturbations are crucial for understanding the universe’s evolution. They manifest as gravitational waves, which can be detected through their influence on the CMB polarization [25] or through direct observations by gravitational wave detectors like LIGO, Virgo, and future missions such as LISA [8]. These observations provide a window into the physics of the early universe, including inflationary models and potential deviations from standard cosmology. The equations for the tensor fluctuations given by

(23)$\begin{array}{ccc}& & T(t){\ddot{h}}_{\mu \nu }+B(t){\dot{h}}_{\mu \nu }-a^{-2}{\nabla }^2{h}_{\mu \nu }=0,\hfill \end{array}$

(24)$\begin{array}{ccc}& & T(t)={\displaystyle \frac{4\kappa -\gamma {\dot{\varphi }}^2}{4\kappa +\gamma {\dot{\varphi }}^2}},\hfill \end{array}$

(25)$\begin{array}{ccc}& & B(t)={\displaystyle \frac{3H(4\kappa -\gamma {\dot{\varphi }}^2)-2\gamma \dot{\varphi }\ddot{\varphi }-4\kappa \dot{\varphi }}{4\kappa +\gamma {\dot{\varphi }}^2}}.\hfill \end{array}$

The expression in Equation (23) has a similar form as in Ref. [14] except for the coefficients of $T(t)$ and $B(t)$. In the linear regime, the components of $h_{\mu \nu }$ can be decoupled and satisfy an evolution equation in the Einstein-Horndeski frame. We can note that the above equations are obtained from the lowest-order effective action and are valid in the low-energy limit. However, in our case, we only have second-order equations with Horndeski corrections. Another way to derive the expression in Equation (23) is through the ADM formalism [14]. For this, we can start with the following action

(26)$\begin{array}{c}\hfill S_{\mathrm{tensor}}^{(2)}={\displaystyle \frac{1}{16\pi G}}\int dt{d}^{3}x{a}^3{e}^{-\varphi }\left[G_T{\dot{h}}_{ij}^2-{\displaystyle \frac{F_T}{a^2}}{({\mathsf{\partial }}_k{h}_{ij})}^2\right],\end{array}$

where the TT gauge is given by ${\mathsf{\partial }}^i{h}_{ij}=0, h_i^i=0$. With this, we can find the equations for the tensor fluctuations as

(27)$\begin{array}{c}\hfill T(t){\ddot{h}}_{ij}+B(t){\dot{h}}_{ij}-a^{-2}{\nabla }^2{h}_{ij}=0,\end{array}$

with

(28)$\begin{array}{c}\hfill T(t)\equiv {\displaystyle \frac{G_T}{F_T}},B(t)\equiv {\displaystyle \frac{{\dot{G}}_T+3H{G}_T}{F_T}}.\end{array}$

We can find that the explicit form of the coefficients $G_T$ and $F_T$ reads [26]

(29)$\begin{array}{ccc}& & G_T\equiv 2\left[G_4-2X{G}_{4X}-X\left(H\dot{\varphi }{G}_{5X}-G_{5\varphi }\right)\right]\hfill \\ & & =2\left(G_4+X{G}_{5\varphi }\right),\hfill \end{array}$

(30)$\begin{array}{c}{F}_T\equiv 2\left[G_4-X\left(\ddot{\varphi }{G}_{5X}+G_{5\varphi }\right)\right]\hfill \\ =2\left(G_4-X{G}_{5\varphi }\right).\hfill \end{array}$

Here, the coefficients are simplified under our choice of $G_i$’s in the action (2). In this case, we have

(31)$\begin{array}{c}\hfill T(t)={\displaystyle \frac{4\kappa -\eta {\dot{\varphi }}^2}{4\kappa +\eta {\dot{\varphi }}^2}},\hspace{1em}B(t)={\displaystyle \frac{3H(4\kappa -\eta {\dot{\varphi }}^2)-2\eta \dot{\varphi }\ddot{\varphi }}{4\kappa +\eta {\dot{\varphi }}^2}}.\end{array}$

The motivation to write the action in Equation (34) is to find the polarization mode $h_A(t,x_{\mu })$ as a function of conformal time $d\eta =dt/a$. Under these considerations the action in Equation (34) becomes

(32)$\begin{array}{c}\hfill S_{\mathrm{tensor}}^{(2)}=\int d\eta d^{3}x{z}^2(\eta )\left[G_T{h}_{ij}^{\prime 2}-F_T{({\mathsf{\partial }}_k{h}_{ij})}^2\right],\end{array}$

where

(33)$\begin{array}{c}\hfill z(\eta )={\displaystyle \frac{a{e}^{-\varphi /2}}{\sqrt{16\pi G}}},\end{array}$

describes the dynamics of a scalar field variable h with a non-minimal coupling to a time-dependent background field $z(\eta )$, which also called the “pump field”. In our cosmic background in Horndeski-like gravity, this field represents a dilaton and an external and internal scale factor (known as “moduli” field). Now, considering the following transformation $u=zh$, we can rewrite the action in Equation (32) as

(34)$\begin{array}{c}\hfill S_{\mathrm{tensor}}^{(2)}=\int d\eta d^{3}x\left[G_T{u}^{\prime 2}+F_{T}u{\nabla }^{2}u+F_T{\displaystyle \frac{z^{″}}{z}}{u}^2\right],\end{array}$

the variation of the action (32) leads to the canonical Schrodinger-like evolution equation as in the form

(35)$\begin{array}{c}\hfill u^{″}-{\displaystyle \frac{F_T}{G_T}}({\nabla }^2+U(\eta ))u=0;U(\eta )={\displaystyle \frac{z″}{z}},\end{array}$

where the gravitational wave speed is defined in terms of

(36)$\begin{array}{c}\hfill c_{GW}^2\equiv {\displaystyle \frac{F_T}{G_T}}.\end{array}$

Beyond, the effective action in Equation (32)

(37)$\begin{array}{ccc}& & S_{\mathrm{tensor}}^{(2)}=\int d\eta d^{3}x\mathcal{L}\hfill \end{array}$

(38)$\begin{array}{ccc}& & \mathcal{L}={\displaystyle \frac{1}{2}}\left[G_T{u}^{\prime 2}-F_T{({\nabla }_{i}u)}^2+F_T{\displaystyle \frac{z^{″}}{z}}{u}^2\right].\hfill \end{array}$

The field can be quantized starting from the definition of the usual momentum density given by

(39)$\begin{array}{c}\hfill 𝒫={\displaystyle \frac{\mathsf{\partial }\mathcal{L}}{\mathsf{\partial }{u}^{\prime }}}=G_T{u}^{\prime }.\end{array}$

Now, we can impose (equal-time) canonical commutation relations

(40)$\begin{array}{c}\hfill [u(x_i,\eta ),𝒫(x_i^{\prime })]=i{\delta }^3(x-x^{\prime }),\end{array}$

on the $\eta =constant$ hypersurface. In our prescription, the classical variable u can be promoted to a field operator and, thus, it can be expanded over a complete set of solutions of the classical solutions in Equation (35). Through the Fourier modes

(41)$\begin{array}{c}\hfill u(x_i,\eta )=\int {\displaystyle \frac{d^{3}k}{{(2\pi )}^3}}[a_k{u}_k{e}^{i{k}_i{x}^i}+a_k^{†}{u}_k^{*}{e}^{-i{k}_i{x}^i}],\end{array}$

we can write the eigenvalue equation

(42)$\begin{array}{c}\hfill u_k^{″}+c_{GW}^2(k^2+U(\eta ))u_k=0;U(\eta )={\displaystyle \frac{z″}{z}}.\end{array}$

Using the Fourier modes, we can provide the commutation relation as

(43)$\begin{array}{ccc}& & G_T[a_k,a_{k^{\prime }}]=G_T[a_k^{†},a_{k^{\prime }}^{†}]=0\hfill \end{array}$

(44)$\begin{array}{ccc}& & G_T[a_k,a_{k^{\prime }}^{†}]={\delta }^3(k-k^{\prime }),\hfill \end{array}$

where we can apply the usual interpretation with $a_k$ annihilation operator and $a_k^{†}$ creation operator. The normalization condition is given by

(45)$\begin{array}{c}\hfill u_k{u}_k^{\prime *}-u_k^{\prime }{u}_k^{*}=i{G}_T,\end{array}$

with $G_T=4\kappa -\gamma a^{-2}{\varphi }^{\prime 2}=4\kappa -\gamma /(a_0^2{\eta }^2)$. The orthonormality is given by the scalar product

(46)$\begin{array}{c}\hfill ⟨{\overline{u}}_k|{\overline{u}}_{k^{\prime }}⟩\equiv -i\int d^{3}x({\overline{u}}_k{\overline{u}}_{k^{\prime }}^{\prime *}-{\overline{u}}_k^{\prime }{\overline{u}}_{k^{\prime }}^{*})={\delta }^3(k-k^{\prime }).\end{array}$

The normalization condition in Equation (45) can now be applied to fix the initial amplitude of the metric’s quantum fluctuations in a typical inflationary scenario. For the asymptotic regime with $\eta \to -\infty$, we have that the $G_T=4\kappa$, i.e., we recover the usual case with $\kappa =1/4$. Besides, using the equations of $\varphi (t)$ and $a(t)$ rewrite in terms of $\eta$ by $d\eta =dt/a$, we can give the form of z as $z=a_0/\sqrt{16\pi G}$, provide $U(\eta )=0$ and we have a free field equation

(47)$\begin{array}{c}\hfill u_k^{″}+k^2{c}_{GW}^2{u}_k=0.\end{array}$

In this case, the expression in Equation (47) has oscillating solutions. In the inflationary regime, the solution exhibits a positive frequency mode on the initial hypersurface. This behavior arises because the scale factor of the Universe, $a(t)$, grows exponentially during inflation, outpacing the growth of the Hubble radius, $H^{-1}$ Figure 2. Consequently, the physical wavelength of each mode stretches faster than the Hubble radius, leading to a scenario where the mode’s wavelength eventually exceeds the Hubble radius—a process commonly referred to as the mode “leaving the horizon”. This phenomenon is critical for understanding the generation of primordial perturbations, as modes that exit the horizon during inflation become “frozen” and later re-enter during the radiation or matter-dominated eras, leaving observable imprints on the Cosmic Microwave Background (CMB) and the large-scale structure of the Universe. Once inflation ends, the dynamics of these modes are governed by the post-inflationary evolution of the Universe [27,28,29].

The general solution for Equation (47) is given by

(48)$\begin{array}{c}\hfill u_k={\displaystyle \frac{e^{i{k}^2{c}_{GW}^2\eta }}{\sqrt{k{c}_{GW}}}},\end{array}$

the gravitational wave speed in the regime $\eta \to -\infty$, becomes

(49)$\begin{array}{c}\hfill c_{GW}^2\equiv {\displaystyle \frac{F_T}{G_T}}={\displaystyle \frac{4\kappa +\gamma /(a_0^2{\eta }^2)}{4\kappa -\gamma /(a_0^2{\eta }^2)}}\to 1=c_{light}^2.\end{array}$

In the context of early universe cosmology, the normalized positive frequency solution in Equation (48) plays a pivotal role in defining the initial vacuum state of tensor fluctuations. This vacuum state corresponds to the lowest-energy configuration associated with the action in Equation (32). For a time-dependent cosmological background, the effective mass term, given by $z^{\prime \prime }/z=0$, ensures that the notion of positive frequency remains invariant over time. As a result, the vacuum state is uniquely defined and globally consistent across the entire cosmological manifold. This property implies that a state initially devoid of particles (i.e., quanta of tensor fluctuations) remains empty throughout the evolution of the universe. Such a framework is critical for understanding the quantum origins of cosmological perturbations and their imprint on observable phenomena, such as the Cosmic Microwave Background (CMB) anisotropies and primordial gravitational waves (see, e.g., [30,31]). These insights provide a foundation for connecting theoretical predictions with empirical observations Figure 3, thereby advancing our understanding of the universe’s earliest moments.

4. Conclusions and Discussions

The constraints set on the speed of propagation of gravitational waves has turned out to be a stringent constraint on general classes of scalar-tensor theories. In the present work, we explore the possibility of a two-dilaton field scenario. Such scenarios have become interesting in recent years with the appearance of problems such as the cosmic tensions issue. Here, the scalar fields may offer a dual action in which one has an action in the early Universe while the other is free to act more globally over time. In our case, we do not constrain these effects but we do explore the tensor perturbations of this setting with a focus on assessing where the class of theories is constrained in any way by the speed constraint on gravitational wave propagation.

Our Horndeski-like model is defined in Equation (6), which contains a two-dilaton action (2) together with the regular matter action $S_m$ and a boundary component (3) that is important for preserving a healthy general relativistic limit (3). The eventual equations of motion are then derived in Section 2 where the dynamical variables are expressed at background level through the FLRW metric. Indeed, as an example, we provide a particular solution for one simple expression of the superpotential that describes the two-dilaton fields. Ultimately, for this case, we find the simple solution expressed in Equations (19)–(22).

The propagation of tensor perturbations in this framework is governed by modified equations of motion that incorporate the effects of the two dilaton fields. These modifications are consistent with the stringent constraints on the speed of gravitational waves, as established by recent multi-messenger observations. The model’s ability to satisfy these constraints while simultaneously addressing the $H_0$ tension highlights its robustness and viability as an alternative to the standard $\mathsf{\Lambda }$CDM paradigm. The unique features of the Horndeski-like gravity framework, including its intrinsic scalar fields and modified gravitational dynamics, offer a promising avenue for probing the observational signatures of alternative cosmological models. In conclusion, the two-dilaton model within the Horndeski-like gravity framework represents a significant step forward in addressing the $H_0$ discrepancy and advancing our understanding of the early universe. By leveraging the interplay between tensor perturbations and modified gravity, this approach provides a comprehensive framework for testing the predictions of alternative cosmological theories against current and future observational data. Future work should focus on refining the model’s parameters and exploring its implications for other cosmological tensions, such as those related to dark energy and the large-scale structure of the universe.

The cosmic expansion history, characterized by the Hubble parameter evolution discussed in our diagram Figure 3, fundamentally connects with gravitational wave propagation through the underlying spacetime geometry. Tensor perturbations, representing the spin-2 field disturbances of the FLRW metric, provide a complementary probe to the scalar sector measurements that drive the Hubble tension [32]. Our analysis in the transverse-traceless (TT) gauge yields the tensor perturbation propagation equation shown in Equation (27), which reveals how gravitational waves traverse the expanding cosmos. Remarkably, the propagation speed derived in Equation (36) approaches the speed of light in its asymptotic limit as expressed in Equation (49), aligning with the stringent constraints from GW170817 and its electromagnetic counterpart [33]. This consistency strengthens our modified gravity framework while simultaneously addressing the $H_0$ discrepancy illustrated in our analysis. The tension between early and late universe measurements may find resolution in the tensor sector’s behavior, complementing potential explanations from early dark energy models [34] or running vacuum scenarios [35]. Our future investigations will extend to scalar perturbations, which govern structure formation and could impose additional constraints on viable cosmological models, potentially narrowing the parameter space where Hubble tension solutions may exist.

The study of tensor perturbations in modified gravity frameworks, such as Horndeski-like gravity, has become increasingly significant in addressing key cosmological challenges. Tensor modes, as spin-2 perturbations of the metric, provide a unique observational window into the dynamics of the early universe, particularly through their imprints on the Cosmic Microwave Background (CMB) and gravitational wave observatories. These modes are instrumental in constraining viable cosmological models, as they directly probe the interplay between the gravitational sector and the underlying spacetime geometry.

In future work, we will address the integration of viscous fluid cosmology into the Horndeski-like framework represents a promising direction for future research. By bridging the gap between modified gravity theories and effective fluid dynamics, this approach has the potential to address some of the most pressing challenges in modern cosmology, including the nature of dark energy, the origin of cosmic acceleration, and the resolution of the Hubble tension.

Author Contributions

All authors contributed to the study’s conception and design. Material preparation, data collection, and analysis were performed by F.F.S. and J.L.S. The first draft of the manuscript was written by F.F.S. and J.L.S. commented on previous versions. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The datasets generated and/or analyzed during the current study are available in the arxiv repository, https://arxiv.org/abs/2401.03558 (accessed on 8 March 2025).

Conflicts of Interest

The authors declare no conflicts of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1 The diagram has been successfully constructed and visualized. It illustrates the connections between Horndeski gravity, inflationary kinematics, and their observational implications, including primordial perturbations.

Figure 2 The diagram has been constructed and visualizes the relationship between the scale factor, Hubble radius, and mode wavelength during inflation. It shows how the mode’s wavelength grows faster than the Hubble radius, eventually leaving the horizon, and re-entering later.

Figure 3 The diagram illustrating the cosmological timeline, the invariance of the vacuum state, and its connection to observable phenomena like the CMB and gravitational waves.

References

1. Haco, S.; Hawking, S.W.; Perry, M.J.; Strominger, A. Black Hole Entropy and Soft Hair. J. High Energy Phys.; 2018; 12, 98. [DOI: https://dx.doi.org/10.1007/JHEP12(2018)098]

2. Haco, S.; Perry, M.J.; Strominger, A. Kerr-Newman Black Hole Entropy and Soft Hair. arXiv; 2019; arXiv: 1902.02247

3. Haco, S. Large Gauge Transformations and Black Hole Entropy. Ph.D. Thesis; University of Cambridge: Cambridge, UK, 2019; [DOI: https://dx.doi.org/10.17863/CAM.48185]

4. Hubble, E.P. A relation between distance and radial velocity among extra-galactic nebulae. Proc. Natl. Acad. Sci. USA; 1929; 15, pp. 168-173. [DOI: https://dx.doi.org/10.1073/pnas.15.3.168] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16577160]

5. Penzias, A.A.; Wilson, R.W. A Measurement of excess antenna temperature at 4080-Mc/s. Astrophys. J.; 1965; 142, pp. 419-421. [DOI: https://dx.doi.org/10.1086/148307]

6. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S. . Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys.; 2020; 641, A6.Erratum in Astron. Astrophys. 2021, 652, C4 [DOI: https://dx.doi.org/10.1051/0004-6361/201833910]

7. Akrami, Y.; Arroja, F.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S. . Planck 2018 results. X. Constraints on inflation. Astron. Astrophys.; 2020; 641, A10. [DOI: https://dx.doi.org/10.1051/0004-6361/201833887]

8. Auclair, P.; Blanco-Pillado, J.J.; Figueroa, D.G.; Jenkins, A.C.; Lewicki, M.; Sakellariadou, M.; Sanidas, S.; Sousa, L.; Steer, D.A.; Wachter, J.M. . Probing the gravitational wave background from cosmic strings with LISA. J. Cosmol. Astropart. Phys.; 2020; 4, 34. [DOI: https://dx.doi.org/10.1088/1475-7516/2020/04/034]

9. Gasperini, M. On the initial regime of pre-big bang cosmology. J. Cosmol. Astropart. Phys.; 2017; 9, 1. [DOI: https://dx.doi.org/10.1088/1475-7516/2017/09/001]

10. Sberna, L. Early-universe cosmology in Einstein-scalar-Gauss-Bonnet gravity. arXiv; 2017; arXiv: 1708.01150

11. Papallo, G. Causality and the Initial Value Problem in Modified Gravity. Ph.D. Thesis; University of Cambridge: Cambridge, UK, 2018; [DOI: https://dx.doi.org/10.17863/CAM.24726]

12. Santos, F.F.; Sokoliuk, O.; Baransky, A. Holographic Complexity of Braneworld in Horndeski Gravity. Fortschr. Phys.; 2023; 71, 2200141. [DOI: https://dx.doi.org/10.1002/prop.202200141]

13. Santos, F.F.; Pourhassan, B.; Saridakis, E.N. de Sitter Versus Anti-de Sitter in Horndeski-Like Gravity. Fortschr. Phys.; 2024; 72, 2300228. [DOI: https://dx.doi.org/10.1002/prop.202300228]

14. Santos, F.F.; Brito, F.A. Domain walls in Horndeski gravity. arXiv; 2004; arXiv: 2105.00343[DOI: https://dx.doi.org/10.1016/j.physletb.2024.138543]

15. Santos, F.F. Rotating black hole with a probe string in Horndeski Gravity. Eur. Phys. J. Plus; 2020; 135, 810. [DOI: https://dx.doi.org/10.1140/epjp/s13360-020-00805-x]

16. Santos, F.F.; Neves, R.M.P.; Brito, F.A. Modeling dark sector in Horndeski gravity at first-order formalism. Adv. High Energy Phys.; 2019; 2019, 3486805. [DOI: https://dx.doi.org/10.1155/2019/3486805]

17. Harko, T.; Lobo, F.S.N.; Saridakis, E.N.; Tsoukalas, M. Cosmological models in modified gravity theories with extended nonminimal derivative couplings. Phys. Rev. D; 2017; 95, 044019. [DOI: https://dx.doi.org/10.1103/PhysRevD.95.044019]

18. Horndeski, G.W.; Silvestri, A. 50 Years of Horndeski Gravity: Past, Present and Future. Int. J. Theor. Phys.; 2024; 63, 38. [DOI: https://dx.doi.org/10.1007/s10773-024-05558-2]

19. Riess, A.G.; Yuan, W.; Macri, L.M.; Scolnic, D.; Brout, D.; Casertano, S.; Jones, D.O.; Murakami, Y.; Breuval, L.; Brink, T.G. . A Comprehensive Measurement of the Local Value of the Hubble Constant with 1 km s⁻¹ Mpc⁻¹ Uncertainty from the Hubble Space Telescope and the SH0ES Team. Astrophys. J. Lett.; 2022; 934, L7. [DOI: https://dx.doi.org/10.3847/2041-8213/ac5c5b]

20. Abbott, T.M.; Aguena, M.; Alarcon, A.; Allam, S.; Alves, O.; Amon, A.; Andrade-Oliveira, F.; Annis, J.; Avila, S.; Bacon, D. . Dark Energy Survey Year 3 results: Cosmological constraints from galaxy clustering and weak lensing. Phys. Rev. D; 2022; 105, 023520. [DOI: https://dx.doi.org/10.1103/PhysRevD.105.023520]

21. Byrnes, C.T.; Wands, D. Curvature and isocurvature perturbations from two-field inflation in a slow-roll expansion. Phys. Rev. D; 2006; 74, 043529. [DOI: https://dx.doi.org/10.1103/PhysRevD.74.043529]

22. Linde, A.D. Particle Physics and Inflationary Cosmology; Contemporary Concepts in Physics; CRC Press: Boca Raton, FL, USA, 1990; Volume 5, pp. 1-362.

23. Rinaldi, M. Mimicking dark matter in Horndeski gravity. Phys. Dark Univ.; 2017; 16, 14. [DOI: https://dx.doi.org/10.1016/j.dark.2017.02.003]

24. Brito, F.A.; Santos, F.F. Braneworlds in Horndeski gravity. Eur. Phys. J. Plus; 2022; 137, 1051. [DOI: https://dx.doi.org/10.1140/epjp/s13360-022-03270-w]

25. Namikawa, T.; Naruko, A.; Saito, R.; Taruya, A.; Yamauchi, D. Unified approach to secondary effects on the CMB B-mode polarization. J. Cosmol. Astropart. Phys.; 2021; 2021, 29. [DOI: https://dx.doi.org/10.1088/1475-7516/2021/10/029]

26. Kobayashi, T. Horndeski theory and beyond: A review. Rept. Prog. Phys.; 2019; 82, 086901. [DOI: https://dx.doi.org/10.1088/1361-6633/ab2429]

27. Armendariz-Picon, C.; Fontanini, M.; Penco, R.; Trodden, M. Where does Cosmological Perturbation Theory Break Down?. Class. Quantum Gravity; 2009; 26, 185002. [DOI: https://dx.doi.org/10.1088/0264-9381/26/18/185002]

28. Flanagan, E.E.; Hughes, S.A. The Basics of gravitational wave theory. New J. Phys.; 2005; 7, 204. [DOI: https://dx.doi.org/10.1088/1367-2630/7/1/204]

29. Kubota, K.i.; Arai, S.; Mukohyama, S. Propagation of scalar and tensor gravitational waves in Horndeski theory. Phys. Rev. D; 2023; 107, 064002. [DOI: https://dx.doi.org/10.1103/PhysRevD.107.064002]

30. Hu, W.; Sugiyama, N.; Silk, J. The Physics of microwave background anisotropies. Nature; 1997; 386, pp. 37-43. [DOI: https://dx.doi.org/10.1038/386037a0]

31. Page, D.N. No superluminal expansion of the universe. Class. Quantum Gravity; 2009; 26, 127001. [DOI: https://dx.doi.org/10.1088/0264-9381/26/12/127001]

32. Perivolaropoulos, L.; Skara, F. Challenges for ΛCDM: An update. New Astron. Rev.; 2022; 95, 101659. [DOI: https://dx.doi.org/10.1016/j.newar.2022.101659]

33. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B. . Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A. Astrophys. J. Lett.; 2017; 848, L13. [DOI: https://dx.doi.org/10.3847/2041-8213/aa920c]

34. Smith, T.L.; Poulin, V.; Amin, M.A. Oscillating scalar fields and the Hubble tension: A resolution with novel signatures. Phys. Rev. D; 2020; 101, 063523. [DOI: https://dx.doi.org/10.1103/PhysRevD.101.063523]

35. Givans, J.J.; Hirata, C.M. Redshift-space streaming velocity effects on the Lyman-α forest baryon acoustic oscillation scale. Phys. Rev. D; 2020; 102, 023515. [DOI: https://dx.doi.org/10.1103/PhysRevD.102.023515]