Study of parametrized dark energy models with a general non-canonical scalar field

In this paper, we consider various dark energy models in the framework of a non-canonical scalar field with a Lagrangian density of the form L(ϕ,X)=f(ϕ)XXMPl4α-1-V(ϕ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{L}(\phi , X)=f(\phi )X{\left( \frac{X}{M^{4}_{Pl}}\right) }^{\alpha -1} - V(\phi )$$\end{document}, which provides the standard canonical scalar field model for α=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1$$\end{document} and f(ϕ)=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\phi )=1$$\end{document}. In this particular non-canonical scalar field model, we carry out the analysis for α=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =2$$\end{document}. We then obtain cosmological solutions for constant as well as variable equation of state parameter (ωϕ(z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }(z)$$\end{document}) for dark energy. We also perform the data analysis for three different functional forms of ωϕ(z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }(z)$$\end{document} by using the combination of SN Ia, BAO, and CMB datasets. We have found that for all the choices of ωϕ(z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }(z)$$\end{document}, the SN Ia +\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$+$$\end{document} CMB/BAO dataset favors the past decelerated and recent accelerated expansion phase of the universe. Furthermore, using the combined dataset, we have observed that the reconstructed results of ωϕ(z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }(z)$$\end{document} and q(z) are almost choice independent and the resulting cosmological scenarios are in good agreement with the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model (within the 1σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\sigma $$\end{document} confidence contour). We have also derived the form of the potentials for each model and the resulting potentials are found to be a quartic potential for constant ωϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }$$\end{document} and a polynomial in ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} for variable ωϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _{\phi }$$\end{document}.


Introduction
One of the biggest challenges in modern cosmology is understanding the nature of the dark energy (DE), which seems to be responsible for the observed accelerated expansion phase of the universe at the present epoch [1,2]. Among the many candidates for DE, the cosmological constant ( ) emerges as the most natural and the simplest possibility. However,cosmology suffers from the so-called "fine tuning" and "cosmic coincidence" problems [3,4]. These theoretical problems motivated cosmologists to think beyond the cosmologa e-mail: abdullaalmamon.rs@visva-bharati.ac.in b e-mail: sudipta.das@visva-bharati.ac.in ical constant and explore other unknown components which may be responsible for the late-time accelerated expansion phase of the universe. The scalar field models have played an important leading role as a candidate of DE due to its dynamical nature and simplicity. Till now, a variety of scalar field DE models have been proposed, such as quintessence (canonical scalar field), k-essence, phantom, tachyon, dilatonic dark energy, and so on (for details, see Ref. [5] and the references therein). But the origin and nature of DE still remains completely unknown, despite many years of research.
It is strongly believed that the universe had a rapid exponential expansion phase during a short era in the very early epoch. This is known as inflation [6,7]; it can give a satisfactory explanation to the problems of the Hot Big Bang cosmology (for example, the horizon, flatness, and monopole problems). Generally, cosmologists realized this inflationary scenario by using a single canonical scalar field called the "inflaton", which has a canonical kinetic energy term (φ 2 2 ) in the Lagrangian density. In the literature, there also exist some inflationary models in which the kinetic energy term is different from the standard canonical scalar field case (instead of the standard formφ 2 2 ). Such models are commonly known as the non-canonical scalar field models of inflation. Such noncanonical scalar fields have been found to have many attractive features compared to the canonical scalar field case, for example, the slow-roll conditions can be achieved more easily as compared to the canonical case. Many interesting possibilities with these models have been recently studied in the literature (see Refs. [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23]). It has been first shown in Refs. [8,9] that the k-essence model (which belongs to an important class of non-canonical scalar field models) is capable of generating inflation in the early epoch. Later, Chiba et al. [10] showed that such models can equally effectively describe a DE scenario. Since the nature of DE is completely unknown, it is quite reasonable to consider a non-canonical scalar field as a candidate for DE component and check for the viability of such models. Within the framework of a non-canonical scalar field, in this work, we shall try to obtain an observationally viable cosmological model to analyze the behavior of the deceleration parameter (q) and the equation of state (EoS) parameter (ω φ ) for describing the expansion history of the universe. The motivation for this work is discussed in detail in Sect. 2. As already mentioned, as the nature of DE is unknown to us, we eventually have no firm idea regarding whether the EoS parameter of DE is a constant quantity or whether it is dynamical in nature. In this connection, the most effective choice is to assume a specific functional form for the dark energy EoS parameter ω φ as a function of the redshift z (for details see Sect. 3.2). To study the non-canonical scalar field DE model in a more general framework, in this paper, we have considered both possibilities. First, we have studied the model for a constant EoS parameter ω φ , which is in the range −1 < ω φ < − 1 3 so as to obtain acceleration. Second, we have considered three different choices for ω φ (z) in order to cover a wide range of the DE evolution. We have then solved the field equations and analyzed the respective cosmological scenarios for all the cases. For all the models, the deceleration parameter q is found to exhibit an evolution from early deceleration to late-time acceleration phase of the universe. This feature is essential for the structure formation of the universe. For all the models, we have also derived the potential V (φ) in terms of the scalar field φ by considering a specific parametrization of f (φ). In order to compare the theoretical models of DE (for dynamical ω φ ) with the observations, we have used the SN Ia, BAO, and CMB dataset to constrain the various model parameters (for details see Appendix A). We have found that the combined dataset favors the CDM model within the 1σ confidence contour. We give the detailed results of this work in Sect. 4.
The present paper is organized in the following way. In the next section, we introduce some basic equations of a general non-canonical scalar field model and also discuss the motivation of this work. We then obtain the general solutions of the field equations for a particular choice of the function f (φ) and for different forms of the EoS parameter ω φ . In Sect. 4, we summarize the results of this work. Finally, the conclusions of this work are presented in Sect. 5. Additionally, for completeness, we perform the combined data analysis in Appendix A and find the observational constraints on ω φ (z) and q(z) using the SN Ia, BAO, and CMB datasets.

Basic framework
Usually, the scalar field models are characterized by a general action which has the following functional form: where R is the Ricci scalar, and L(φ, X ) is the Lagrangian density, which is an arbitrary function of the scalar field φ and its kinetic term X . The kinetic term X is defined as X = 1 2 ∂ μ φ∂ μ φ, which is a function of time only. The last term, S m , represents the action of the background matter. Throughout this paper we shall work in natural units, such that 8π G = c = 1.
The expressions for the energy density (ρ φ ) and pressure ( p φ ) associated with the scalar field are given by where X = 1 2φ 2 . In general, the Lagrangian for a scalar field model can be represented as (Melchiorri et al. [24]) where f (φ) and F(X ) are arbitrary functions of φ and X , respectively. V (φ) is the potential for the scalar field φ. Let us consider a homogeneous, isotropic, and spatially flat FRW universe which is characterized by the following line element: where a(t) is the scale factor of the universe. With the FRW geometry, the equations of motion take the form where H =ȧ a denotes the Hubble parameter, an overdot indicates differentiation with respect to the time coordinate t, and ρ m represents the energy density of the matter component of the universe, F X ≡ ∂ F ∂ X and F X X ≡ ∂ 2 F ∂ X 2 . It deserves mentioning that Eq. (4) includes all the popular single scalar field models. It reduces to a canonical scalar field model when f (φ) = constant = 1 and F(X ) = X . Again, it describes a pure k-essence model when V (φ) = 0 and a phantom scalar field model when f (φ) = 1 and F(X ) = −X . It is interesting to note that Eq. (4) reduces to the general non-canonical scalar field model [L(φ, X ) = F(X ) − V (φ)] when f (φ) = 1. This type of non-canonical scalar field models was proposed by Fang et al. [21]. They studied several aspects of this type of scalar fields for different forms of F(X ). Recently, these types of non-canonical scalar field models have gathered attention due to their simplicity. Unnikrishnan et al. [11] have showed that for noncanonical scalar field models, the slow-roll conditions can be more easily satisfied compared to the canonical inflationary theory. They have shown that such models (with quadratic and quartic potentials) are more consistent with the current observational constraints relative to the canonical inflation. They have also shown that such non-canonical models can drop the tensor-to-scalar ratio rather than their canonical counterparts. In fact, a lot of work has been done in the framework of the non-canonical inflationary scenario in the early epoch [8][9][10][11][12][13][14][15][16][17][18][19][20][21]. Furthermore, Franche et al. [25] showed that the non-canonical scalar fields are the most universal case with a general Lagrangian density satisfying certain conditions. These interesting properties of non-canonical scalar field models motivated us to study the cosmological aspects of such fields in a more general framework in the context of dark energy. In the literature, a large number of functional forms of L(φ, X ) have been proposed so far; see for example [11,[26][27][28]. In our earlier work [22,23], we have considered a Lagrangian density of the following form: which can be obtained from the general form of the Lagrangian density [11,27,28] for α = 2 and M Pl = The above equation describes a purely canonical scalar field Lagrangian density In this present work, we try to extend our previous work in [22,23] by considering a general non-canonical scalar field model which has the following Lagrangian density: Here, following [22,23] also, we consider α = 2 and M Pl = In this case, the energy density and pressure of the scalar field are given by It is evident from Eqs. (13) and (14) that the usual definition of ρ φ (= 1 2φ 2 + V (φ)) and p φ (= 1 2φ 2 − V (φ)) for a standard canonical scalar field model gets modified due to the Lagrangian density (12). Also, the equations of motion for this Lagrangian (Eqs. (6)-(9)) turn out to be The solution for ρ m from Eq. (18) is obtained: where ρ m0 is the matter density at the present time, z = a 0 a −1 is the redshift, and the present value of the scale factor a 0 is normalized to unity.
One of the important quantities in cosmology is the dark energy EoS parameter ω φ = p φ ρ φ , which, in our case, is given by From Eq. (17), one can then obtain the expression for the energy density of the scalar field as where ρ φ0 is an integration constant. Now, the Friedmann equation can be written in the following integrated form: where H 0 is the Hubble parameter at the present epoch, parameters of the matter and scalar field (or dark energy), respectively, at the present epoch. For this model, from Eqs. (13) and (14), the potential can be expressed (as a function of redshift z) as and In order to solve the field equations analytically we will proceed as follows. Out of Eqs. (15)-(18), only three are independent equations in view of the Bianchi identity (with five unknown quantities: H , ρ m , f (φ), V (φ), and φ). So naturally one has to assume two relationships among the different variables to solve the system of equations.
Following the above argument, in this paper, we assume that the quantity f has the functional form where f 0 is an arbitrary constant. It deserves mentioning here that the above parametrization of f (φ) helps us to close the system of equations. With this input, Eq. (24) can be written in the following integral form: where φ 0 is an arbitrary constant of integration and . Another important observable quantity, the deceleration parameter q(z), can also be expressed in terms of H (z) as which describes the evolution of our universe.
We shall now concentrate on the dark energy EoS parameter ω φ (z). If a function of ω φ (z) is given, then we can find the evolution of ρ φ (z) from Eq. (21). As a result, we can also find the evolutions of H (z), q(z), V (z), and φ(z). Inverting φ(z) into z(φ) and using Eq. (23), one can then obtain the potential V (φ) in terms of φ. As already mentioned, we have considered a specific parametrization of f (φ) and still we need another assumption to match the number of unknown parameters with the number of independent equations. With this freedom, we choose different functional forms for ω φ (z), the equation of state parameter. In the next section, we try to obtain the functional forms of various cosmological parameters for different choices of ω φ (z) and study their cosmological implications.

Theoretical models
In this section, we shall consider two phenomenological DE models for obtaining the current acceleration of the universe in the framework of a general non-canonical scalar field theory.
3.1 Model I: Accelerating universe driven by a constant EoS parameter for dark energy (−1 < ω φ < − 1 3 ) In this model, we shall investigate the properties of an accelerated expanding universe driven by a non-canonical scalar field dark energy with a constant EoS parameter. Recent observations suggest that the dark energy EoS ω φ is very close to −1 and the approximate bound on ω φ is −1.1 ≤ ω φ ≤ −0.9 [29,30]. Keeping this limit in mind, we choose a constant ω φ in the limit −1 < ω φ < − 1 3 , which ensures that the model does not deviate much from a CDM model.
In this case, the energy density of DE can be obtained (from Eq. (21)): where ρ φ0 is the integration constant which represents the dark energy density at the present time. Then the corresponding Friedmann equation becomes Inserting the Hubble parameter H (z) (as given in the above equation) into Eq. (27), we obtain the deceleration parameter as where and κ 2 = 1− m0 m0 . The plot of q(z) against z is shown in Fig. 1 for different values of ω φ (within the range −1 < ω φ < − 1 3 ) and m0 = 0.3. Figure 1 shows that q(z) crosses its transition point from its positive value regime to the negative value regime in the recent past, which is consistent with the independent measurements reported by several authors (see Ref. [31] and the references therein).
With the help of Eqs. (26) and (28), one can solve for the scalar field in the flat FRW universe as With the help of Eqs. (23) and (31), we find the form of the potential in terms of φ as where V 0 = . Thus, the constant ω φ model leads to a quartic potential. For φ 0 = 0, the above potential is similar to the potential used by Linde [7] in the context where f α = It is evident from the above equation that the function f (φ) can be arranged as a series expansion in powers of (φ − φ 0 ).

Model II: Accelerating universe driven by time-dependent EoS parameter for dark energy
We shall now focus on the second possibility i.e., the EoS parameter ω φ is dynamical in nature. For this purpose, in this subsection, we have considered three different choices of ω φ to study the model in a more general way. If ω φ is dynamical in nature, then one way to study models going beyond the cosmological constant is by using a particular functional form for the dark energy EoS parameter ω φ (z). However, a large number of functional forms for ω φ (z) have appeared in the literature [32][33][34][35][36][37][38][39][40][41]. Usually, the parametrized form of ω φ (z) is written as [5] where the ω n 's are arbitrary constants and the x n (z)'s are functions of the redshift z. The numerical values of the ω n 's can be found by fitting to the observational data. Following the first order expansions in Eq. (34), several authors con-sidered many functional forms for x n (z) to investigate the evolution of ω φ (z). For example: (i) ω φ (z) = ω 0 = constant (as we have discussed in model I) for x 0 (z) = 1 and x n = 0 (n ≥ 1). (ii) ω φ (z) = ω 0 + ω 1 z i.e., linear redshift parametrization [32,33], for x n (z) = z n with n ≤ 1. (iii) ω φ (z) = ω 0 +ω 1 log(1+z) i.e., logarithmic parametrization [36], for x n (z) = [log(1 + z)] n with n ≤ 1. (iv) ω φ (z) = ω 0 +ω 1 z (1+z) i.e., CPL parametrization [34,35], for x n = z 1+z n with n = 1.
There are many more. It is worth mentioning that the parametrizations (ii) and (iii) diverge at high redshifts, whereas the fourth one blows up in the future, when z → −1.
It should be noted that the assumed parametrization would lead to possible biases in the study of the evolution of the DE but in the absence of any information regarding the true nature of DE, these parametrizations provide some insight regarding the possible nature of DE component and are worth studying. In this paper, we shall consider two different divergence-free functional forms of ω φ (z) which do not diverge in the future (z → −1). In addition, we shall also consider the linear redshift parametrization of ω φ (z) for the statistical model comparisons with the divergence-free parametrizations, at low redshifts. In order to explore the evolution of DE, we shall also try to reconstruct the deceleration parameter q(z) using Eq. (27) for these different choices in Sect. 4.

• Assumption I
Here, we consider the linear redshift parametrization of the EoS parameter ω φ [32,33], which has the following functional form: where ω 0 represents the present value of ω φ and the second term measures the variation of ω φ with respect to z. Inserting for ω φ (z) from Eq. (35) into Eq. (21), we obtain Now Eq. (22) can be written as In the context of DE (as it is a late-time phenomenon), the above choice of ω φ (z) has been widely used due to its simplicity and we find for the present parametrization of f given in Eq. (25) that the potential turns out to be a polynomial in φ.

• Assumption II
Next, we propose where ω 2 and ω 3 are arbitrary constants to be fixed by observations. It is easy to see that the EoS parameter reduces to for z → +∞ (early epoch), for z → −1 (far future).
• Assumption III The next choice adopted in this paper was suggested by Alam et al. [39,40], which is the functional form This choice is exact and gives the cosmological constant where ρ φ0 is the present value of ρ φ .
In this case, the Hubble parameter is expressed as We have then solved Eqs. (23) and (26) and However, in general for Model II (which includes ansatz I, II, and III) with the particular choice of Eq. (25), one can write the potential V as a polynomial in φ in the following manner: where n > 0, the V i 's are constants, and the values of these parameters are different for different choices of ω φ (z). Interestingly, we have found that it is a generalization of other well-known potentials (see [5] and the references therein), for example, a constant potential or a power-law potential. We have also found that the parametrization (25) leads to the quantity f (φ) as an exponential, polynomial, and power law in φ for choices I, II, and III, respectively. In the following section, we shall use these choices to discuss the possibility of constraining ω φ (z) and q(z) from observations.

Results
Following the statistical analysis (see Appendix A), in this section, we present the fitting results for different choices of the EoS parameter for DE. Figure 5 shows the 1σ and 2σ confidence contours for each choice (I, II, and III) using the SN Ia + BAO/CMB dataset. The best-fit values of the model parameters and ω φ (z = 0) for these different choices are given in Table 1.
Using those best-fit values, we have reconstructed the deceleration parameter q(z) for each model and the results are plotted in Fig. 6.
It is evident from Fig. 6 that q(z) shows a smooth transition from a decelerated (q > 0) to an accelerated (q < 0) phase of expansion of the universe at the transition redshift z t = 0.38 (for ansatz I), 0.36 (for ansatz II), and 0.43 (for ansatz III) for the best-fit models. These results are in good agreement with those obtained by several authors based on various other considerations [42][43][44].
Furthermore, we also show the reconstructed evolution history of the EoS parameter in Fig. 7 for each choice of ω φ (z).
We have also plotted the total EoS parameter, which is defined as ω tot (z) = p φ ρ φ +ρ m , as a function of z for these choices (see the inset diagram of Fig. 7). From Table 1, we have found that the current values of ω φ (z) for the best-fit DE models are very close to −1, i.e., the models do not deviate very far from the CDM model (ω = −1) at the present epoch. However, as indicated in Table 1, the present parametrized model favors a phantom model (ω φ < −1) in the 2σ limit and thus requires further attention.

Conclusions
In this work, we have studied various non-canonical scalar field DE models in a spatially flat, homogeneous, and isotropic FRW space-time. In this framework, we have obtained the general solutions of the field equations for different choices of the EoS parameter. For completeness, we have also investigated how the joint analysis of the SN Ia + BAO/CMB datasets constrains the redshift evolutions of q(z) and ω φ (z) for different choices of ω φ (z) (as given in Model II). In Fig. 5, we have also shown the 1σ and 2σ contour plots of the pairs (ω 0 , ω 1 ) (upper panel), (ω 2 , ω 3 ) (middle panel), and (A 1 , A 2 ) (lower panel) for the ansatzes I, II, and III, respectively. In this analysis, we have also calculated the best-fit values of the free parameters (as shown by large dots in Fig. 5) and it has been found that the chosen values of these parameters (which were chosen for solving the parametric relations in Appendix B) are well fitted within the 1σ confidence contour (as shown by small dots in Fig. 5).
We have shown that the deceleration parameter q undergoes a smooth transition from its deceleration phase (q > 0, at high z) to an acceleration phase (q < 0, at low z) for all of the considered parametrized models. However, as mentioned in the previous section, the value of z t , where the signature flip of q (from the decelerating to an accelerat-ing expansion phase) takes place has been calculated and the results obtained are consistent with the present day cosmological observations. From the SN Ia + BAO/CMB analysis, we have also found q(z = 0) = −0.56, −0.64, and −0.60 for ansatzes I, II, and III, respectively, which also agree very well with the recent observational results.
From Table 1, we have observed that the EoS parameter ω φ (z = 0) ≈ −1, but slightly less than −1 for all three choices (as discussed in Sect. 4). As we have seen ω φ (z = 0) ≈ −1, our models do not deviate very far from the CDM model (see also Fig. 7), which is currently known as the standard model for modern cosmology. In order to gain more physical insight into these time evolutions of the EoS parameter, we have also plotted the reconstructed total EoS parameter ω tot (z) in Fig. 7 (see the inset diagram of Fig.  7). For each choice, this figure shows that ω tot (z) attains the required value of − 1 3 around z = 0.62 (within 1σ confidence level) and remains always greater than −1 up to the present epoch. These scenarios also agree very well with the observational data.
However, the models presented here are restricted because the form of f (φ) chosen was ad hoc (as given in Eq. (25)) and did not follow from any principle. In this regard, we have mentioned earlier that we make this choice in order to close the system of equations. With this choice of f (φ), we have derived the form of the potential V (φ) in terms of φ for different models. We have found that Model I leads to a quartic potential, whereas Model II leads to a polynomial potential for each choice of ω φ (z). We have seen that, with a suitable choice of V i 's for the potential (as given in Eq. (51)), it is possible to reproduce the other well-known potentials in the context of DE. However, many possibilities are opened up to accommodate a physically viable potential for other parametrizations of f (φ) or f (H ). Finally, we would like to emphasize that all the considered models provide a deceleration for high redshift and an acceleration for low redshift as required for the structure formation of the universe. However, these results are completely independent of any choice of f (φ). With the increase of more good quality observational data at low, intermediate, and high redshifts, the constraints on z t (or q(z)) and ω φ (z) are expected to get improved in the near future.
for assumptions I (i = 0, j = 1) and II (i = 2, j = 3). For the sake of simplicity, we have reduced the three dimensional parameter space ( m0 , ω i , ω j ) into the two dimensional plane (ω i , ω j ) by fixing m0 to some constant value. On the other hand, h(z) contains four free parameters ( m0 , A 0 , A 1 , and A 2 ) for assumption III. In this case, we have also reduced the four dimensional parameter space ( m0 , A 0 , A 1 , A 2 ) into the two dimensional plane (A 1 , A 2 ) by fixing m0 and A 0 to some constant values. Now, we can deal with only two free parameters for each ansatz and will perform a χ 2 analysis of the SN Ia + BAO/CMB dataset. The values of the model parameters at which χ 2 m (the minimum value of χ 2 function) is obtained are the best-fit values of these parameters for the joint analysis of the observational datasets from SN Ia, BAO, and CMB measurements.
To reconstruct V (φ), we proceed as follows. One can easily find that it is not possible to express V (φ) in terms of φ explicitly, because φ(z) poses a very complicated form for each choice of ω φ (z). These equations only give a parametric representation of V (φ), which cannot be solved analytically. Therefore, one can plot the potential V (φ) against φ for some arbitrary values of the model parameters. After this, one can obtain the form of V (φ) by using a fitting function to fit the corresponding plot. Following this procedure, we have plotted V (φ) as a function of φ for these choices (I, II, and III), which are shown in Figs. 2, 3, and 4. Similarly, using the parametric relations [ f (z), φ(z)] for each choice of ω φ (z), we also obtain the form of f (φ) by a numerical method for some given values of the model parameters.