Cosmic slowing down of acceleration with the Chaplygin-Jacobi gas as a dark fluid

One of the greatest challenges in modern science is understanding the universe, particularly its origin, evolution, and final fate. In this regard, the hot big bang (HBB) cosmology, based on four-dimensional General Relativity (GR) Einstein:1916vd , has been widely accepted as the standard paradigm describing how the universe expanded from an initial singularity. During the longest part of its lifetime, the universe has undergone a decelerating expansion, being dominated first by radiation and then by matter. However, the cosmic history includes two phases of accelerated expansion: one at very early times and another at late times. The first accelerating phase corresponds to inflation, which supplies an explanation of the observed large scale structure (LSS), as it provides the primordial, almost adiabatic and scale-invariant perturbations that undergo gravitational collapse to form galaxies and clusters of galaxies Abazajian:2013vfg . The second accelerating phase corresponds to the current cosmic acceleration, which is supported from a large number of observational evidence Riess1998 ; Perlmutter1998 ; WMAP2003 ; Aghanim2020 ; cooke2018one ; Eisenstein2005 , which also indicates that our universe is spatially almost flat, and currently dominated by dark energy (DE) and cold dark matter (DM). The accelerated expansion of the present universe is attributed to DE, which is an exotic component having negative pressure, such as the cosmological constant $\Lambda$ Padmanabhan:2002ji ; Sahni:1999gb ; Carroll:2000fy . The $\Lambda$CDM model, has proven to be successful in explaining a wide range of cosmological observations Riess1998 ; Perlmutter1998 ; WMAP2003 ; Aghanim2020 ; cooke2018one ; Eisenstein2005 . Despite its success, it faces several challenges, namely: (i) the cosmological constant problem Weinberg1989 ; Martin:2012bt , (ii) the cosmic coincidence problem (or the so-called “why now” problem) Zlatev1998 ; Arkani-Hamed:2000ifx , and more recently, (iii) the tension between measurements of the Hubble parameter $H_{0}$ DiValentino2021 ; Riess:2020fzl and (iv) the $S_{8}$ tension ($S_{8}=\sigma_{8}\sqrt{\Omega_{m}/0.3}$) DiValentino2020 .

Several approaches have been proposed to address $\Lambda$CDM problems. They generally fall into two groups: one modifies GR, while the other, dynamical DE, involves models where the properties of DE vary over time, unlike in the $\Lambda$CDM. Several candidates for dynamical DE have been studied in recent years, including the Chaplygin gas (CG). For a review of DE models, see Refs. Copeland:2006wr ; Nojiri:2006ri ; Clifton:2011jh ; Joyce:2016vqv . The CG’s equation of state (EoS) has a connection with string and brane theories Bilic:2001cg ; Heydari-Fard:2007vcg , and it can admit a supersymmetric generalization Jackiw2000 . Chaplygin cosmology can be understood as the study of the dark sector through a CG fluid that satisfies the relation $p=-\frac{B}{\rho}$ Kamenshchik:2001cp ; Fabris:2002xx ; Gorini:2002kf , and its generalizations, like the Generalized Chaplygin Gas (GCG) with an equation of state $p=-\frac{B}{\rho^{\alpha}}$ Bilic:2001cg ; Bento:2002ps , where $p$ represents pressure (assumed to be negative to produce an accelerated stage at late times), $\rho$ is the energy density, and $B$ is a positive constant (to ensure $p<0)$ and $0<\alpha\leq 1$. These models belong to the unified models class, which provide a description of the matter-dominated era and the late accelerating era with only one component rather than two (for other unified models, see e.g. Hipolito2009 ; Hipolito2010 ). However, observationally, the CG model has shown problems with instabilities in their corresponding trajectories Perrotta2004 and when tested with various observations, such as Type Ia Supernovae (SnIa) , X-ray gas mass fraction of clusters, and Hubble rate-redshift data biesiada2005 ; Colistete:2005yx ; Zhu:2004aq ; Wu:2006pe ; Makler:2002jv . However, although the GCG model is in agreement with background observational data biesiada2005 ; Colistete:2005yx ; Zhu:2004aq ; Wu:2006pe ; Makler:2002jv ; delCampo:2009cz , it suffers from an unexpected blow-up in the matter power spectrum caused by adiabatic pressure perturbations Sandvik:2002jz . This undesired effect can be avoided if some kind of non-adiabaticity is introduced in the model (see e.g. vomMarttens:2017cuz ). Moreover, in recent years, several GCG modifications or generalizations have been proposed to study the dark sector. Among others, some of them are the modified Chaplygin gas (MCG) Benaoum2022 ; Debnath2004 , the new generalized Chaplygin gas (NCG) Zhang2004 , and viscous generalized Chaplygin gas (VGCG) Zhai2005 ; Hernandez-Almada:2021osl . For recent observational results involving those models, see e.g. Zheng2022 .

Recently, a further generalization has been proposed through the Hamilton-Jacobi formalism in the inflationary context by means of elliptic functions Rengo1 . This led to the development of a more general Chaplygin-like EoS, which we will refer to as Chaplygin-Jacobi Gas model (CJG). The CJG EoS is expressed as follows:

$$p=-\frac{B,k}{\rho^{\alpha}}-k^{\prime}\rho\left(2-\frac{1}{B}\rho^{\alpha+1}\right)\,.$$

(1)

In the above equation, $k$ represents the modulus of the elliptic function, and $k^{\prime}=1-k$ is the complementary modulus. It is noteworthy that the GCG is obtained with $k=1$ and $B>0$. In recent years, the CJG has been studied in contexts different from this work Rengo1 ; Rengo2 ; Rengo3 ; Debnath:2021pxy ; Mukherjee2023 ; Chaudhary2023l ; Rengo4 .

Considering that Chaplygin-like fluids offer an interesting framework to study phenomenology beyond the $\Lambda$CDM, the main goal of the present work is to explore the viability of the CJG for describing the late universe. First, we integrate the balance equation and find its solutions for the energy density, pressure and EoS parameter. Then, we perform an analysis of its parameter space to find the region in which physical solutions are possible. After obtaining an analytical solution for the Hubble rate, and in order to compare the background evolution of the CJG with that of the $\Lambda$CDM model, we introduce parameters where derivatives of the scale factor beyond the second-order appear. To this end, one option is to study the so-called statefinder parameters, $r$ and $s$, defined as follows Sahni:2002fz ; Alam:2003sc .

	$\displaystyle r$	$\displaystyle\equiv$	$\displaystyle\frac{\dddot{a}}{aH^{3}},$		(2)
	$\displaystyle s$	$\displaystyle\equiv$	$\displaystyle\frac{r-1}{3(q-\frac{1}{2})},$		(3)

where the dot denotes differentiation with respect to cosmic time $t$, $H=\dot{a}/a$ represents the Hubble rate, and $q=-\ddot{a}/(aH^{2})$ denotes the deceleration parameter. It is worth noting that the statefinder parameters involve third derivatives of the scale factor with respect to cosmic time, in contrast to the Hubble rate and the deceleration parameter, which are expressed in terms of the first and second time derivatives of the scale factor, respectively. Statefinder analysis serves as a diagnostic tool for understanding the dynamics of the universe’s expansion. These parameters are computed for our specific DE model, and as we will discuss later, $r$ and $s$ can differ significantly from $\Lambda$CDM even if they predict very similar expansion histories.

While emergent cosmology may not realize a unified model like most CG-like models, it offers a phenomenology with interesting results distinct from those of the CG model and its generalizations, particularly in future times. As is expected for any DE model, a fluid satisfying the relation (1) is subdominant in a multi-component universe during early times. This accurately reproduces a radiation-dominated era followed by a matter-dominated era, resulting in a decelerating universe. Subsequently, the DE fluid dominates, leading to an accelerating phase. However, a specific region in the parameter space permits a transient acceleration-deceleration in the future, indicating the possibility of a decelerating phase following the current accelerating phase. This phenomenon is commonly referred to as the cosmic slowing down of the current acceleration. DE models with a constant EoS like $\Lambda$CDM or the $w$CDM cannot exhibit the slowing down feature Zhang2018 , as they predict a de Sitter (dS) phase as the final fate of the universe. However, depending on the set of parameters, some DE models, especially those with dynamical EoS or parametrizations like the Chevallier-Polarski-Linder (CPL) model and others, allow for this possibility (see e.g. Vargas2011 ; Hu2015 ; Magana2014 ; Zhang2018 ; Bolotin2020 ). Here, we demonstrate that a universe with a fluid satisfying the relation (1) as DE can also exhibit the slowing down behavior. Additionally, we conduct a data analysis to investigate whether recent astronomical and cosmological observations support this transient behavior in the context of the CJG. To achieve this, we utilize the latest supernova (SnIa) Pantheon + data and Cosmic Chronometers (CC) data. Furthermore, we incorporate observations from Fast Radio Bursts (FRBs) into our tests, as they have recently been shown to provide an interesting complement to other cosmological probes.

Our work is organized as follows: after this introduction, in Section II, we derive solutions for the energy density, pressure, and EoS parameter. Next, we conduct an analysis of the parameter space to determine the sub-regions where physical solutions are feasible and the sub-regions where the slowing down phenomenon appears. In Section III, we delve into the background cosmological dynamics of the CJG in a multifluid context. This includes obtaining an analytical solution for the Hubble rate $H(z)$ as a function of redshift and introducing the deceleration and Statefinder parameters tailored for the CJG. In Section IV, we present the data and relevant equations for the implementation of the MCMC analysis. In Section V, we present the results of the MCMC analysis and evaluate numerically the cosmic evolution of the CJG against redshift, comparing these results with those of the $\Lambda$CDM model. Finally, we summarize our findings and present our conclusions in Section VI. Throughout our work, we adopt the mostly positive metric signature $(-,+,+,+)$ and utilize natural units where $c=\hbar=1$.

As mentioned before, the Hamilton-Jacobi formalism allows for obtaining a particular generalization of the Chaplygin inflationary model through elliptical functions, resulting in a more general Chaplygin-like EoS parameter Rengo1

$$\omega=-\frac{B\,k}{\rho^{\alpha+1}}-2k^{\prime}+\frac{k^{\prime}}{B}\rho^{\alpha+1},$$

(4)

where $k^{\prime}=1-k$, and $\alpha$, $B$ and $k$ are constants with $0<k<1$ and $0<k^{\prime}<1$. In principle, the EoS in (1) reduces to a GCG case when $k=1$, $B>0$ and $0<\alpha<1$. The case where $B>0$, $k=1$ and $\alpha=0$ reproduces $\Lambda$CDM. There is no other case in which one can obtain the MGC, NCG, or another Chaplygin-like EoS. In the GCG context, only the first term in Eq. (1) appears. Then, for $B>0$ and $0<\alpha<1$, the pressure is always negative, allowing an explanation for the late accelerated expansion of the universe. However, the CJG has two additional contributions to the pressure. The second term of Eq. (1) is always negative, but the first and third terms always have opposite signs. For instance, for $B>0$, the first term is negative and the third term is positive, while for $B<0$ the behaviour is reversed. Eventually, depending on the parameter space, the positive term might dominate and the pressure might become less negative, opening the the possibility of a slowing down of the cosmic acceleration. This becomes more evident when looking at the energy density and pressure as functions of the scale factor. Let us first consider one universe dominated by a fluid with an EoS as in Eq. (1) and solve the energy balance equation.

In the context of a universe described by the Friedman-Lemaître-Robertson-Walker (FLRW) metric, the energy balance equation for an EoS given by Eq. (1) reads:

$\displaystyle\frac{\dot{\rho}}{\rho}=-3\frac{\dot{a}}{a}\left(-\frac{B\,k}{\rho^{\alpha+1}}+2k-1+\frac{(1-k)}{B}\rho^{\alpha+1}\right)=-3\frac{\dot{a}}{a}\left(\frac{-Bk-(1-2k)\rho^{\alpha+1}+\frac{1-k}{B}\rho^{2(\alpha+1)}}{\rho^{\alpha+1}}\right),$

(5)

where $a$ is the scale factor of the homogeneous and isotropic metric and a dot denotes differentiation with respect to cosmic time. This equation can be integrated directly, resulting in

$\displaystyle\rho(a)=\left[B-\frac{BA}{Ak^{\prime}-a^{3(1+\alpha)}}\right]^{\frac{1}{1+\alpha}}\,,$

(6)

where $A$ is an integration constant. If $B=\rho^{1+\alpha}_{0}B_{s}$, this last equation can be conveniently rewritten as

$\displaystyle\rho(a)=\rho_{0}\left[B_{s}-\frac{B_{s}A}{Ak^{\prime}-a^{3(1+\alpha)}}\right]^{\frac{1}{1+\alpha}}\,,$

(7)

and the integration constant $A$ is determined by the condition $\rho=\rho_{0}$ at the present time, so that

$\displaystyle A=\frac{1-B_{s}}{1-k(1-B_{s})},$

(8)

and thus,

$\displaystyle\rho(a)=\rho_{o}\left[B_{s}-\frac{B_{s}(1-B_{s})}{(1-B_{S})k^{\prime}-(k^{\prime}+B_{s}\,k)a^{3(1+\alpha)}}\right]^{\frac{1}{1+\alpha}}\,,$

(9)

Notice that the energy density spans a three-dimensional space of parameters. Accordingly, the pressure of the CJG as a function of the scale factor can be computed by replacing Eq. (9) into Eq. (1). Although, in principle, solution (9) has a three-dimensional space of parameters

$${\cal{M}}=\{(\alpha,B_{s},k)/\,\alpha\in\mathbb{R},B_{s}\in\mathbb{R},0<k<1\},$$

(10)

some regions of this space can lead to singularities, divergences in the past, or negative densities, resulting in non-physical solutions across the entire space ${\cal{M}}$. However, certain sub-regions of ${\cal{M}}$ allow solutions without these problems. Let us find the sub-regions in ${\cal{M}}$ where we can have physical solutions for the model.

Here we study a cosmological model based on the results from the previous section, using the CJG fluid as the DE component. First, we assume the cosmic substratum to be dynamically dominated by a mixture of radiation ($r$), pressureless dark matter ($m$), and a DE component ($de$). Radiation and matter are modeled by barotropic fluids with EoS parameters $w_{r}=1/3$ and $w_{m}\ll 1$, respectively. The DE energy component comes in the form of a CJG fluid with EoS (1).

In General Relativity, the Friedmann equation for a homogeneous, isotropic, and spatially flat three-component universe is

$\displaystyle 3H^{2}=8\pi\,G(\rho_{r}+\rho_{m}+\rho_{de})\,.$

(24)

Here, $\rho_{r}$, $\rho_{m}$, and $\rho_{de}$ are the radiation, pressureless dark matter, and DE densities, respectively. Moreover, energy conservation is satisfied by all the components separately

$\displaystyle\dot{\rho}_{i}+3H(\rho_{i}+p_{i})=0\,,$

(25)

where $i=r,m,de$. As it is well-known, solutions of Eq. (25) for radiation and matter are found to be $\rho_{r}=\rho_{r0}a^{-4}$ and $\rho_{m}=\rho_{m0}a^{-3}$, respectively. For DE, the solution is given by Eq. (9). The sub-index “0” denotes the current value of any given quantity.

The past evolution is restricted by the necessity of a radiation and matter-dominated epochs to guarantee, for instance, the Big Bang Nucleosynthesis (BBN) and cosmic structure formation. To be a viable DE model, the CJG must be subdominant in the past, reproducing correctly a radiation domination era and subsequently, the matter domination era giving a decelerating universe. In recent times, the DE fluid must dominate to give rise the accelerating stage of the universe. A fluid with EoS given by (1) and within a subspace of parameters $\cal{N}$ reproduces all these points properly, as can be observed in Fig. 1 (left side), where the fractional densities are shown as a function of redshift. In all cases, the blue curves represent radiation, the red ones represent matter, and the black curves represent DE fluid. For the dot-dashed curve, parameter values are $\Omega_{m0}=0.30$, $\alpha=0.1$, $B_{s}=-0.6$, $k=0.43$, corresponding to a point in ${\cal{N}}_{3}$. For the solid line, we used $\Omega_{m0}=0.30$, $\alpha=1.3$, $B_{s}=-1.4$, $k=0.46$, corresponding to a point in ${\cal{N}}_{4}$. For the sake of comparison, we also plotted the corresponding functions in the $\Lambda$CDM model using the same colored dashed curves.

On the other hand, the Hubble rate can be written as a function of the scale factor $a$. However, it will be more useful to have the dimensionless Hubble rate as a function of $z$

$\displaystyle E^{2}(z)=\Omega_{r0}(1+z)^{4}+\Omega_{m0}(1+z)^{3}+(1-\Omega_{m0}-\Omega_{r0})\left[B_{s}+\frac{B_{s}(1-B_{s})(1+z)^{3(1+\alpha)}}{(k^{\prime}+B_{s}\,k)-(1-B_{s})k^{\prime}(1+z)^{3(1+\alpha)}}\right]^{\frac{1}{1+\alpha}},$

(26)

where we have used the fact that each energy density is usually expressed in terms of the dimensionless density parameter, defined as $\Omega_{i}=\rho_{i}/\rho_{cr}$ with $\rho_{cr}=3H^{2}/(8\pi G)$ being the critical density.

Having obtained the analytical expression for $E(z)$, we may proceed to analyze our model using the statefinder diagnostic. It will be useful to express these parameters in terms of the dimensionless Hubble rate and the redshift $z$. The cosmic times derivatives are written as redshift derivatives according to

$$\frac{d}{dt}=\frac{d}{dz}\frac{dz}{da}\frac{da}{dt}=-(1+z)H(z)\frac{d}{dz},$$

In this way, the set of statefinder parameters becomes

	$\displaystyle r$	$\displaystyle=1+\frac{1}{2E^{2}}\left[(1+z)^{2}\frac{d^{2}E^{2}}{dz^{2}}-2(1+z)\frac{dE^{2}}{dz}\right],$		(27)
	$\displaystyle s$	$\displaystyle=\frac{1}{3}\frac{(1+z)^{2}\frac{d^{2}E^{2}}{dz^{2}}-2(1+z)\frac{dE^{2}}{dz}}{(1+z)\frac{dE^{2}}{dz}-3E^{2}}.$		(28)

Accordingly, upon replacement of Eq. (26) into Eqs. (27)-(28), the specific expressions for the Statefinder parameters are obtained (not shown) and will be plotted in Section V using the best-fit values. Furthermore, we compute the effective EoS parameter $w_{eff}$ for our CJG cosmological model, which encodes information about the universe’s composition, the evolution of its energy density, and the dynamics of its expansion. The general expression for $w_{eff}$ depends on the first derivative of $E^{2}(z)$ with respect to the redshift $z$ as follows

$\displaystyle w_{eff}=-1+\frac{1+z}{3E^{2}}\frac{dE^{2}}{dz}.$

(29)

Thus, by replacing Eq.(26) into (29), $w_{eff}$ for our specific CJG model is found (not shown).

As discussed above, for a universe dominated only by a CJG, a particular region in the subspace of parameters allows for transient acceleration-deceleration in the future. This feature must remain when a more realistic model of the universe is considered, including radiation and dark matter. Therefore, the future cosmological evolution within the framework of CJG may be very different from a de Sitter phase. This is shown in Fig. 1 (right side), where several curves for deceleration parameter as a function of redshift were plotted for different sets of parameters. The dashed dark curve on the right side of Fig. 1 represents the case for for $\Lambda$CDM, in which, as it is known, the accelerated expansion will continue forever. The green curve, which corresponds to the set $\Omega_{m}=0.3$, $\alpha=1.1$, $B_{s}=-0.81$ and $k=0.55$, also indicates accelerated expansion forever. However, the last two curves, the blue one for $\Omega_{m}=0.3$, $\alpha=0.1$, $B_{s}=-1.5$, and $k=0.43$, and the orange one for $\Omega_{m}=0.3$, $\alpha=0.8$, $B_{s}=-1.41$, and $k=0.42$, exhibit a transition between an accelerated and a decelerated phase in the future. In all the curves, the transition from a decelerated evolution during the matter-dominated epoch to an accelerating phase in the past can be observed.

The possibility of future transient acceleration is not new. In the past, evidence was found for a slowing down of the expansion rate of the universe, or equivalently, for an increase in the deceleration parameter at redshifts close to the present epoch $z\approx 0$ Shafieloo2009 ; Guimaraes2011 ; Cai2011 ; Costa2010 ; Fabris2009 ; Vargas2011 . More recently, this possibility has been explored in light of recent data Magana2014 ; Shahalam2015 ; Hu2015 ; Bolotin2020 ; Escobal2023 . In fact, over the last few years, several authors have discussed the possibility that the accelerated expansion might be a transient phenomenon, i.e., that there might be a transition back to decelerated expansion Albrecht1999 ; Barrow2000 ; Bento2001 .

So far, we have explored the concept of cosmic slowing down with a CJG fluid from a theoretical standpoint. However, in recent years, a large amount of cosmological data has emerged, enabling very precise statistical tests. In the next section, we will use different probes to study the viability of a CJG as dark energy and what the data can tell us about the cosmic slowing down in the CJG framework. In this regard, several quantities describing the background evolution, such as the Hubble rate , the deceleration parameter, and the statefinder parameters, will be evaluated after obtaining the best-fit values from the aforementioned observational analysis.

In Section II, the regions in the parameter space where the model has physical solutions were limited by imposing physical conditions on the fluid density (9). In this section, we will utilize current observational data to estimate the free parameters and their associated error bars. Through this process, we aim to determine the region in the parameter space where a model based on CJG could be compatible with observations. We will utilize cosmic chronometers (CC), Type Ia Supernovae (SNIa), and Fast Gamma Ray Bursts (FRBs) datasets, along with their joint analysis, to perform a MCMC study. Implementations were carried out using the publicly available Python package Polychord (handley2015polychord, ). Now, let us briefly describe each of the datasets used.

We have implemented four different MCMC analyses using: a) the CC observational data only, b) SNIa observational data only, c) FRBs observational data only, and d) a joint analysis including these three probes. In order to perform all analyses, we fixed $\Omega_{r}=5.38\times 10^{-5}$ to be consistent with CMB observations ParticleDataGroup:2020ssz , and considered the following parameters as free: $H_{0}$, $\Omega_{m}$, $\alpha$, $B_{s}$, and $k$. For consistency, we used the same flat priors in all cases, such as $[65,80]$ for $H_{0}$ and $[0,1]$ for $\Omega_{m}$. For $\alpha$, $B_{s}$, and $k$, we suitably considered the sub-regions in ${\cal{N}}$. Moreover, in the analyses using only SNIa and CC+SNIa+FRBs, an additional parameter $M$ appears. In that case, we used a Gaussian prior for $M$, with mean value of $-19.24$ and standard deviation $\sigma_{M_{B}}=0.04$, as suggested in Camarena2021 .

Results for all four analyses are presented in Fig. 2 and Table 1. All three analyses using only one probe (CC, SNIa, or FRBs) exhibit good agreement within the 1$\sigma$ statistical confidence level (C.L.), indicating concordance between the CJG cosmology in the region ${\cal{N}}$ and observations. The best-fit values for the CC data and FRBs data support a de Sitter phase in the future, while the SNIa data support a slowing-down phase. However, neither the CC nor FRBs data exclude a transient acceleration. Nevertheless, the large error bars found for the best-fit values might limit informativeness. These error bars can be reduced in the joint analysis, where we found $H_{0}=71.7^{+1.4}_{-1.3}$, $\Omega_{m}=0.296\pm 0.028$, $\alpha=1.42^{+1.3}_{-0.77}$, $B_{s}=-1.55^{+0.79}_{-0.67}$, and $k=0.41^{+0.11}_{-0.15}$. The contour plots corresponding to the joint analysis are also displayed in Fig. 2, where some general features can be observed. In general, regions with $1+\alpha>0$ and $B_{s}<0$ are preferred by observations; i.e., the data favor the $\cal{N}_{3}\cup\cal{N}_{4}$ region, with a slight trend toward $\cal{N}_{3}$. This can be seen in Fig. 3 (left), where the plane $B_{s}$-$k$ is shown together with the curve $B_{s}=-\frac{k^{\prime}}{k}$ (black dashed curve) separating $\cal{N}_{3}$ (de Sitter future) and $\cal{N}_{4}$ (non-de Sitter future). Notice that the best-fit point indicated by $\mathrm{x}$ is in the slowdown region. Moreover, there is a positive correlation between the $B_{s}$ and $k$ parameters, which can be explained by the $k$-dependence of the intervals for $B_{s}$ (see Sec. II).

With the best-fit values from the joint analysis (CC+SNIa+FRBs), we can plot relevant quantities that describe the cosmic evolution at the background level and compare them with those from the $\Lambda$CDM model. Fig. 3 (right) shows the plot of relative deviation between the Hubble function of CJG and that of the $\Lambda$CDM model, $1-H_{CJG}(z)/H_{\Lambda\text{CDM}}$. As can be seen from Fig.3, $H(z)$ becomes indistinguishable for both models during the past, while it approaches very close values in the present time, resulting in a very small negative relative difference. In future times, $H(z)$ for CJG becomes smaller than that for the $\Lambda$CDM model, leading to a positive relative deviation. Conversely, cosmic acceleration is accounted for by the deceleration parameter, which is depicted in Fig. 4. As can be seen, the best-fit values (black curves) support cosmic deceleration in the future. These findings are interesting because they present a different behavior compared to the $\Lambda$CDM case (blue line on the left side of Fig.4) and other types of Chaplygin-like fluids such as NCG, MCG, and GCG, as shown on the right side of Fig.4¹¹1To plot the deceleration curve for Chaplygin-like models, we used the best-fit values for NCG, MCG, and GCG found in Zheng2022 . Nevertheless, it is important to note that the 1$\sigma$ region (gray region in Fig. 4) is relatively large and sufficiently inclusive to also encompass a de Sitter final fate for the universe. In general, further data and analyses are needed to refine our understanding of the final fate of the universe.

A complementary analysis is provided by the statefinder parameters $r$ and $s$, which account for the deviations of the CJG model from $\Lambda$CDM. In Fig. 5 (left), the trajectory traced by the CJG on the $s-r$ plane is shown, obtained by parametrically plotting Eqs. (27) and (28) with varying $z$. If the role of DE is played by a cosmological constant, the value of $r$ remains constant at $r=1$ throughout the matter-dominated epoch and into the future (for $z\lesssim 10^{4}$). In the $s-r$ plane, the fixed point $\{s,r\}=\{0,1\}$ corresponds to $\Lambda$CDM case. For CJG case, at past times ($z>0$), the $(s,r)$ pair is $\sim(0,1)$, indistinguishable from $\Lambda$CDM. However, currently, the parameters are $s\simeq 0.27$ and $r\simeq 0.2$, indicating a deviation of the CJG trajectory from $\Lambda$CDM. In Fig. 5 (right), the effective EoS parameter $w_{eff}$ (29) is displayed. The solid-blue and dashed-black lines correspond to the CJG and $\Lambda$CDM models, respectively. The effective EoS parameter $w_{eff}$ provides a summary of the combined effects of various components in the universe. As observed, $w_{\text{eff}}$ changes during the matter-dominated and dark energy-dominated epochs, reflecting the dominant component’s EoS parameter at that time. For CJG and $\Lambda$CDM, the effective EoS parameters are practically indistinguishable in the past. However, around the present epoch, in CJG case, $w_{eff}$ starts to increase and at some point in the future it becomes positive, indicating cosmic deceleration. This corroborates the late-time behavior of CJG initially observed from the deceleration parameter (Fig. 4).

We have explored an extension of the Chaplygin EoS, derived from a specific generalization of the Chaplygin inflationary model within the framework of the Hamilton-Jacobi formalism and elliptic functions Rengo1 . By imposing conditions to obtain physical solutions for the energy density, we reduced the parameter space to different sub-regions where the model is free of singularities within the interval $0<a<1$, and the density remains always positive. Interestingly, one of these sub-regions is compatible with the possibility of a new deceleration phase in the future, following the current acceleration. Then, we used different data samples such as SnIa, CC and FRBs to investigate the viability of the CJG cosmological model including radiation and matter. Constraints on the free parameters we obtained using a statistical MCMC analyses and their error bars are presented in Table 1. Best-fit values for the joint analysis support the cosmic slowing down phenomenon (see Fig. 4) with a maximum of acceleration at $z\sim 0.02$, $\alpha$ values such that $\alpha+1>0$ and negative values for $B_{s}$. The preference for negative values for the $B_{s}$ parameter implies a lack of direct correspondence with the GCG. Additionally, the best-fit values result in a cosmology almost indistinguishable from $\Lambda$CDM in the past, deviating from it only in recent epoch, and diverging from the eternal de Sitter expansion.

Exploring CJG expands two families of models. Firstly, it increases the number of models within the Chaplygin type family. Unlike other models in this family, however, CJG allows for a future acceleration-deceleration transition. Secondly, it enlarges the family of models that predict future deceleration, such as interacting models, $w$ parametrizations, decaying dark matter, extensions of quintessential cosmology, and others (see, for example, Shafieloo2009 ; Vargas2011 ; Magana2014 ; Escobal2023 , and references therein). Moreover, it has been observed that a potential future transition to a decelerating stage, thereby naturally ending the eternal accelerating regime, is significant from a physical perspective Escobal2023 . An eternal de Sitter phase disagrees with the requirements of the S-matrix describing particle interactions and can be problematic within the framework of string theory Fischler2001 ; Cline2001 ; Hellerman2001 . Nevertheless, until a definitive answer regarding the final fate of the universe is reached, preferably in a model-independent and completely data-driven way, the study of DE dynamics, even through phenomenological approaches, is crucial for enhancing our understanding of cosmic evolution and the future of the universe. Our aim here was to explore the potential of the CJG as a DE model with a final phase different than the de Sitter. Furthermore, future cosmological data may offer further insights into these intriguing questions.

J.A.S.F. thanks FAPES for the fellowship and financial support (045/2024 - P: 2024-M88VK). Also, J.A.S.F. is grateful for the hospitality of the Instituto Argentino de Radioastronomia (IAR), where part of this work was carried out. W.S.H.R. thanks FAPES (PRONEM No 503/2020) for the financial support under which this work was carried out. N. V. acknowledges the financial support of Fondecyt Grant 1220065. J. R. V. is partially supported by Centro de Física Teórica de Valparaíso (CeFiTeV). N. V. and J. R. V. thank the CEUNES/UFES for their warm hospitality while some of this work was carried out.