Cosmological solutions of chameleon scalar field model

We investigate cosmological solutions of the chameleon model with a non-minimal coupling between the matter and the scalar field through a conformal factor with gravitational strength. By considering the spatially flat FLRW metric and the matter density as a non-relativistic perfect fluid, we focus on the matter-dominated phase and the late-time accelerated-phase of the universe. In this regard, we manipulate and scrutinize the related field equations for the density parameters of the matter and the scalar fields with respect to the e-folding. Since the scalar field fluctuations depend on the background and the field equations become highly non-linear, we probe and derive the governing equations in the context of various cases of the relation between the kinetic and potential energies of the chameleon scalar field, or indeed, for some specific cases of the scalar field equation of state parameter. Thereupon, we schematically plot those density parameters for two different values of the chameleon non-minimal coupling parameter, and discuss the results. In the both considered phases, we specify that, when the kinetic energy of the chameleon scalar field is much less than its potential energy (i.e., when the scalar field equation of state parameter is ≃-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\simeq - 1 $$\end{document}), the behavior of the chameleon model is similar to the ΛCDM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda CDM$$\end{document} model. Such compatibility suggests that the chameleon model is phenomenologically viable and can be tested with the observational data.


Introduction
It has been suggested that the universe starts from an extremely rapid accelerated-phase after the Planck era called the inflationary epoch that remedy important problems of the standard cosmological model (the ΛCDM model), see, e.g., Refs.[1]- [7].The cosmic microwave background observations are well consistent with the predictions of such an inflation hypothesis [8].In the literature, a lot of work has been carried to investigate the inflationary scenario, see, e.g., Refs.[9]- [17].After the inflationary epoch, observations predict that the universe has a decelerating era during the early phase of its evolution [18].In fact, the universe undergoes a radiation dominated epoch followed by a matter dominated one.Afterward, the present accelerated-phase of the universe starts which is dubbed as a dark energy era.Indeed, observations -such as the type Ia supernova, the baryon acoustic oscillations, and the gravitational waves -confirm the late-time cosmic acceleration [8,18]- [25].But the driver of such acceleration has not been identified by observations and has no fully accepted theoretical model.However, since ordinary matter cannot accelerate the universe, such dark energy (if it exists) must be exotic matter for which many candidates have been proposed and investigated.In most of these models, a scalar field has been employed as dark sector with a dynamical equation of state.On the other hand, modified gravitational theories (e.g.Refs.[26]- [41] and references therein) have also been proposed mainly to describe the accelerated evolution of the universe without any sector as dark energy.For these two issues see, e.g., Refs.[42]- [60] and references therein.Furthermore, other observations -such as the behavior of the galactic rotation curves and the mass discrepancy in clusters of galaxies -indicate the consideration of an exotic matter called dark matter [61]- [64].A lot of alternative efforts have also been performed on various modifications to the Einstein field equations in order to deal with the question of dark matter [45,47,65]- [71] and references therein.Howsoever, the nature of both of these dark sectors (which constitute about 95% of the universe [8,18] and are well consistent with the observational data [72]) is one of the most important issues in physics.
Although general relativity provides accurate predictions in describing some cosmological phenomena [73], the ΛCDM is the most acceptable model due to its high conformity with the observational data of cosmology [74].However, this model still confronts some challenges [75]- [83] and other alternatives have been proposed for it.The scalar-tensor theories of gravitation, which extend general relativity by introducing a scalar field, have become the most popular alternative to the Einstein gravitational theory [84]- [94].
In this work, we concentrate on the chameleon scalar field theory, wherein unlike the common nonminimally coupled scalar field theories, a chameleon field couples with matter sector rather than the geometry.Such an interaction between the scalar field and matter field causes the effective potential of the scalar field depends on the matter density of the environment.Much work has been done on this theory, see, e.g., Refs.[10,37,59, 95]- [113] and references therein.This change of properties of the mass of scalar field acts as a buffer with respect to observational bounds.Accordingly, in this type of theory, a screening mechanism is usually aimed at modifying general relativity on large scales.Hence, the chameleon scalar field theory can be employed to describe the accelerated expansion of the universe while the effects of non-minimal coupling parameter are hidden in the small scale gravitational experiments, see, e.g., Refs.[95]- [98,101,102].
We have previously investigated the chameleon model during inflation in Refs.[10,107], wherein we showed through the proposed scenario that the effects of inflation and chameleon can be described via a single scalar field during the inflation and late-time.Also, we have probed the role of the chameleon scalar field as dark energy in Ref. [59], where it was claimed that such a model justifies dark energy with stronger confirmation.Now, in this work, we intend to investigate the cosmological solutions of the chameleon model in the hope that these solutions can better justify the observational data.Of course, since in the chameleon model, the properties of the scalar fluctuations depend on the background, its field equations of motion become highly non-linear.However, rather than solving the resulted equations analytically, we restrict the solutions to proceed.Indeed, for various cases of the relation between the kinetic and potential energies of the chameleon scalar field (or actually, for some specific cases of the scalar field equation of state parameter), we solve the non-minimal coupled gravity equations to understand the evolution of the scalar field.Meanwhile, some attempts have been done to find analytic solutions for the chameleon model, wherein these analytic solutions of the chameleon model exist only for highly symmetric source shapes such as spheres, plates, and ellipses [95,114,115]).If the shape of the matter source is irregular, it will be more difficult to find analytical solutions to the equations of motion.Nevertheless, in Ref. [116], a software package called SELCIE has been introduced that provides tools for constructing an arbitrary system of mass distributions and then computing the corresponding solution to the chameleon scalar field equation.In addition, through the dynamical systems technique [117,118], one can also study the full non-linearity of these cosmological models.In this technique, usually normalized dimensionless new variables are introduced, then by finding the fixed/critical points of the system and their stability, the evolution of the system can be pictured qualitatively near these points.The dynamical systems technique has also been performed for the chameleon scalar field, see, e.g., Refs.[100,103,105,112].
The work is organized as follows.In the next two sections, while benefiting from our previous works [10,59,107], we first introduce the chameleon model and then, by considering the spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) metric, we obtain the corresponding field equations of motion.Thereupon, in Sec.III, we manipulate and scrutinize the related field equations with respect to the e-folding for better and more use in the later sections.In Secs.IV and V, the cosmological solutions of the chameleon model are respectively investigated for the matter-dominated phase and the late-time accelerated-phase of the universe for various cases of the relation between the kinetic and potential energies of the chameleon scalar field.Indeed, we probe the results of the chameleon model for some specific cases of the scalar field equation of state parameter.Meanwhile, in Sec.V, we take a brief look at the chameleon scalar field profile.At last, we conclude the work in Sec.VI with the summary of the results.

Chameleon Model with Scalar Field
We consider an action of the chameleon model with a scalar field in four dimensions when there are a minimal coupling between the Einstein gravity and that scalar field, and a non-minimal coupling between matter species with it as In this action, ϕ is a scalar field, V (ϕ) is a self-interacting potential, ψ (i) s are various matter fields, L (i) s are Lagrangians of matter fields, g(i)µν s are matter field metrics that are conformally related to the Einstein frame metric as Here, β (i) s are dimensionless constants that represent different non-minimal coupling parameter between the scalar field and each matter species, however in this work, we only consider a single matter component.Also, the lower case Greek indices run from zero to three, g is the determinant of the metric, R is the Ricci scalar, and the reduced Planck mass is M Pl ≡ (8πG) −1/2 ≈ 10 27 eV in the natural units of h = 1 = c.Typically in the literature, the common power-law chameleon potential is usually used, where n is a positive or negative integer constant 1 and M is some positive constant mass scale.The variation of action (1) with respect to the scalar field gives the field equation where the box symbol 2 ≡ ∇ α ∇ α corresponds to the metric g µν , and T (m) µν is the energy-momentum tensor of matter (which is conserved in the Jordan frame) defined as Also, the variation with respect to the metric g µν associated to the Einstein frame gives where µν is the energy-momentum tensor of the scalar field as In addition, we assume the matter field as a perfect fluid in the Jordan frame with the linear barotropic equation of state p(m) = w ρ(m) .Hence, one obtains the trace of the energy-momentum tensor (5), with the signature +2, as where the relation between the matter density in the Einstein frame with the Jordan frame is The matter density ρ(m) is conserved in the Jordan frame, i.e.
however, it is not conserved in the Einstein frame.Nevertheless, the mathematical quantity independent of the scalar field and with the same equation of state parameter, is a conserved quantity in this frame, i.e. ρ + 3H (1 + w) ρ = 0.
Of course, this ρ(t) is not a physical matter density.Furthermore, substituting definition (11) into Eq.( 4) yields the dynamic of the scalar field governed by an effective potential, i.e. where which depends explicitly on the matter density.Additionally, the mass of the chameleon field, which is sufficiently large to evade local constraints [95,120], is In the following, we derive the relevant field equations in order to investigate the cosmology of such a chameleon model.

Scrutinizing Field Equations
In this section, we intend to manipulate and scrutinize the cosmological equations of the chameleon model while considering the spatially flat FLRW metric in the Einstein frame, namely where a(t) is the scale factor as a function of the cosmic time t.In this respect, first, the field equation (13) gives the wave equation where we have also considered the scalar field just as a function of the cosmic time, dot denotes the derivative with respect to this time, and H(t) ≡ ȧ/a is the Hubble parameter.Furthermore by metric (16), Eq. ( 6) gives the corresponding Friedmann and Raychaudhuri equations as and Also, relation (7) gives the energy and the pressure densities as Accordingly, Eqs. ( 18) and ( 19) can be rewritten as where ρ (tot) ≡ ρ (ϕ) + ρ (m) and p (tot) ≡ p (ϕ) + p (m) are the total energy and pressure densities.The field equations ( 21) and ( 22) indicate that ρ (tot) is conserved in the Einstein frame, i.e.
However, in general, the scalar and matter fields are not separately conserved in this frame, and an interacting term stands among those, i.e. and wherein the interacting term depends on the scalar field and the background matter density of the environment.However, for an ultra-relativistic matter, with w = 1/3, such an interacting term is identically zero.That is, for the radiation dominated epoch of the universe, the scalar and matter fields are separately conserved.
In the following analysis, we consider the matter density as a non-relativistic perfect fluid,2 i.e. dust matter with w = 0. Accordingly, Eqs. ( 24) and (25) where w (ϕ) = p (ϕ) /ρ (ϕ) .Moreover, we present the equations via the dimensionless density parameters defined as where ρ Pl is the critical density of the universe at the present-time.Hence, Eq. ( 21) can be rewritten as which in turn imposes constraint Ω on the initial values at the present-time.
On the other hand, by employing the e-folding variable we have Also from relations ( 20) and ( 28), we obtain It has been shown [120] that the monotonic increase of the matter coupling factor e βϕ M Pl , when corresponds to φ > 0, leads to a minimum for V eff (ϕ).Moreover, It has been stated [121] that φ > 0 is satisfactory from the astrophysical point of view [122].Hence, we confine our investigation 3 to the positive sign in relation (33) to proceed, and the generalization to the negative sign can be performed with similar considerations.Thus, Eqs. ( 26) and (27) read and where the prime denotes the derivative with respect to the e-folding.In addition, in the case of w = 0, the defined total equation of state parameter for this model is Also, using the deceleration parameter, i.e. q ≡ −äa/ ȧ2 = −ä/(aH 2 ), while substituting Eqs. ( 21) and ( 22), with employing the definition of the total equation of state parameter, into it, leads to Obviously if q < 0 (which is equal to w (tot) < −1/3), it will describe an accelerated-phase in evolution of the universe.
In the continuation of the work, we intend to investigate the cosmological solutions of the chameleon model for two important phases of the cosmos.
On the other hand, in the matter-dominated phase, Eq. ( 21) reduces to and hence, Eqs. ( 34) and (35) read and In the above, one obviously has ρ (m) > 0 and H > 0, wherein we have assumed that4 ρ (ϕ) > 0, i.e. φ2 /2 >| V (ϕ) | for negative potentials.Also, we have assumed that which, with relations (20) and the assumption ρ (ϕ) > 0, it is already satisfied.At this stage, it is more instructive to define two new parameters y ≡ √ Ω (m) and x ≡ √ Ω (ϕ) .Thus, we can rewrite Eqs. ( 42) and (43) as the following coupled first derivative equations for x and y, namely and However, rather than solving these coupled equations analytically, we restrict the solutions to proceed.Indeed, in the following subsections, we investigate the solutions of these coupled differential equations for various cases of the relation between the kinetic and potential energies of the chameleon scalar field.Actually, we probe the results of the chameleon model for four specific cases of the scalar field equation of state parameter.

Case φ2 /2 ≫| V (ϕ) |
In this case, one obviously has w (ϕ) ≃ +1, hence Eqs. ( 45) and ( 46) read 2x Now, by applying the upper bound value of the β = 3.7×10 2 (which this value of the non-minimal coupling parameter in the chameleon model has been shown [123] to be consistent with the experimental constrain), the solutions of these equations are and where C 1 and C 2 = x(0) are constants of integrations.Using solutions ( 49) and ( 50), we have plotted Ω (m) = y 2 and Ω (ϕ) = x 2 in Figure 1, wherein (without loss of generality) the present values of the parameters have been used to specify the constants. 5This figure indicates that the behavior of the matter density and the chameleon scalar field density are both damped oscillations, i.e. the value of both densities decreases over time.Note that, the negative values of the e-folding in this figure (and also in the subsequent figures) are because we have assumed the value of the scale factor to be equal to one at the present-time, i.e. a(t 0 ) = 1.We should also remind that the zero values of Ω (m) and Ω (ϕ) must be disregarded in all the figures.On the other hand, by considering the non-minimal coupling parameter as β = 1, Eqs. ( 47) and (48) lead to solutions and where C 3 and C 4 = x(0) are constants of integrations.Also, using solutions ( 51) and ( 52), we have plotted Ω (m) and Ω (ϕ) in Figure 2.This figure indicates that at the beginning of this phase of the universe the matter density dominates, and then at the end of its dominance era, the chameleon scalar field density is dominant.
The solutions to these equations with β = 3.7 × 10 2 are and x (N ) = e − 3 2 N C 5 sin 185 where C 5 = −y(0) and C 6 = x(0) are constants of integrations.Using solutions ( 55) and ( 56), we have plotted Ω (m) and Ω (ϕ) in Figure 3.This figure indicates that the behavior of the matter density and the chameleon scalar field density are both also damped oscillations and their values decrease over time.However, using the non-minimal coupling parameter β = 1, Eqs. ( 53) and ( 54) lead to where C 7 = −y(0) and C 8 = x(0) are constants of integrations.Again, we have plotted Ω (m) and Ω (ϕ) from solutions ( 57) and ( 58) in Figure 4.This figure also indicates that at the beginning of this phase the matter density dominates, and then at the end of its dominance era, the chameleon scalar field density is dominant.

Case φ2
In this case, the chameleon potential must again be positive, and hence relations (39) and (40) will hold.Although, one gets w (ϕ) = −1/3, that is why we have highlighted this case, which also brings almost different results.Accordingly, Eqs. ( 45) and ( 46) read The solutions to these equations with β = 3.
where C 9 and C 10 = x(0) are constants of integrations.This time, using solutions ( 61) and ( 62), we have plotted Ω (m) and Ω (ϕ) in Figure 5.This figure indicates that the behavior of the matter density and the chameleon scalar field density are both also damped oscillations and their values decrease over time.
On the other hand, by considering the non-minimal coupling parameter as β = 1, Eqs. ( 59) and ( 60) lead to solutions and where C 11 and C 12 = x(0) are constants of integrations.Again using solutions ( 63) and ( 64), we have plotted Ω (m) and Ω (ϕ) in Figure 6.This figure also shows that at the beginning of this phase, the matter density dominates as assumed, but at the end of its dominance era, the chameleon scalar field density is dominant.
These equations do not depend on the β parameter, and we have the solutions and where y (0) = Ω (m) 0 and x (0) = Ω (ϕ) 0 are the initial conditions at the present-time. 7Once again, using solutions ( 67) and ( 68), we have plotted Ω (m) and Ω (ϕ) in Figure 7.This figure also indicates that at the beginning of this phase, the matter density dominates as assumed, but at the end of its dominance era, the constant density of the chameleon scalar field is dominant.

Late-Time Accelerated-Phase
In this section, we investigate the evaluation of the late-time accelerated-phase of the universe via the chameleon scalar field.In this regard, in this accelerated-phase of the universe, by reducing the dust matter density over time, we can assume ρ (ϕ) ≫ ρ (m) .Hence, the equation of state parameter (36) reduces to which gives the deceleration parameter (37) as q ≃ 1 + 3w (ϕ) /2.Then, in order to achieve q ≤ 0 in this phase of the universe, one must use w (ϕ) ≤ −1/3, which leads to the constraint φ2 ≤ V (ϕ) with positive potentials in the chameleon model.On the other hand, in the accelerated-phase of the universe, Eq. ( 21) yields Substituting Eq. ( 70) into Eqs.( 34) and (35) gives and Once again, in the following subsections, we probe these solutions for various cases of the relation between the kinetic and potential energies of the scalar field that satisfy the φ2 /2 ≤ V (ϕ) /2 constraint with positive potentials in the chameleon model.Meanwhile, while using the common power-law chameleon potential (3), we take a brief look at the chameleon scalar field profile in this phase.

Case φ2
In this case, we have w (ϕ) = −1/3, and in turn q = 0, which is the transition point between the decelerated-epoch and the accelerated-phase in the evaluation of the universe.Hence, Eqs. ( 71) and (72) and The solutions to these equations are and We have plotted Ω (m) and Ω (ϕ) in Figure 8.This figure shows that, in this phase of the universe, although the matter density increases slightly, its amount is very negligible all along the path, and the chameleon scalar field density, while decreasing, is dominant at all times.On the other hand, for β = 1, Eq. (71) gives Figure 8: For solutions ( 75) and ( 76) with β = 3.7 × 10 2 , the figure schematically (i.e., scale-free) shows the parameters Ω (m) and Ω (ϕ) versus N as solid and dashed lines, respectively.Besides, the present values of the parameters, i.e.Ω (m) 0 ≃ 0.3 and Ω (ϕ) 0 ≃ 0.7, have been used as the initial conditions, and a(t 0 ) = 1.and, while using the assumption Ω (m) ≪ Ω (ϕ) in the late-time accelerated-phase of the universe, Eq. ( 72) yields The solutions to Eqs. ( 77) and (78) are and We have plotted Ω (m) and Ω (ϕ) in Figure 9.This figure indicates that the non-relativistic matter density and the chameleon scalar field density both decrease with increasing the scale factor.Although the decrease in the matter density during this epoch of the evolution of the universe is less than that of the chameleon, the chameleon scalar field density dominates at all times.
Continuing the analysis, let us take a brief look at the chameleon scalar field profile while using the common power-law chameleon potential (3).In the case of this subsection, such a potential leads to Therefore, the chameleon scalar field is where A(t = 0) is an integration constant.Furthermore, the constraint V (ϕ) > 0 indicates that, for the case ϕ > 0, the parameter n can be both odd and even, but for case ϕ < 0, the parameter n can only be even.
5.2 Case φ2 /2 ≪ V (ϕ) In this case, we have w (ϕ) ≃ −1, and in turn q < 0, hence Eqs. ( 71) and (72) read and These equations do not depend on the β parameter, and we have the solutions and We have plotted solutions ( 85) and ( 86) in Figure 10.This figure illustrates that first the matter density dominates, then with the increase of the scale factor, the chameleon scalar field density, while always remaining constant, dominates the late-time accelerated-phase of the universe.Also, in the case of subsection 5.2, the common power-law chameleon potential (3) leads to Hence, the chameleon scalar field is where B(t = 0) is an integration constant.Now, to make our investigations more instructive, let us compare our results with the corresponding one in the ΛCDM model.For the ΛCDM model, we have with the field equations where the energy-momentum tensor of the cosmological constant is defined as This tensor describes a vacuum state with a constant energy density ρ (Λ) and a constant isotropic pressure density p (Λ) as Assuming the FLRW metric ( 16) with the matter density as a dust matter perfect fluid, one obtains the Friedmann equations for the ΛCDM model as Furthermore, the energy-momentum conservation in this model leads to Using the dimensionless density parameters, these are that lead to solution (85) and solution like (86), respectively.Therefore, in the case φ2 /2 ≪ V (ϕ) with positive potentials in the chameleon model, i.e. subsections 4.4 and 5.2, the behavior obtained for the presented chameleon model is similar to the ΛCDM model.

Conclusions
By considering the spatially flat FLRW line element as the background geometry, we have investigated the cosmological solutions of the chameleon model.In this model, the scalar field non-minimally couples with the matter field, and its interaction with the ambient matter goes through a conformal factor that leads to a dependence of the chameleon mass on the matter density.Accordingly, the equations of motion of the chameleon scalar field become highly non-linear, hence rather than solving the resulted equations analytically, we have restricted the solutions to proceed.
After deriving the corresponding Friedmann and Raychaudhuri equations, we have manipulated and scrutinized the related field equations for the dimensionless density parameters of the matter field and the scalar field with respect to the e-folding while treating the matter density as a nonrelativistic perfect fluid and all variables simply as a function of the cosmic time.Then, we have focused and investigated the cosmological solutions of the chameleon model in the matter-dominated phase and the late-time accelerated-phase of the universe in the context of various cases of the relation between the kinetic and potential energies of the chameleon scalar field, or indeed, for some specific cases of the scalar field equation of state parameter.Thereupon, we have schematically plotted those density parameters versus the e-folding for two different values on the chameleon non-minimal coupling parameter.Meanwhile, we have shown that the assumption of the matter density much more than the chameleon density corresponds to the matter-dominated phase with a decelerated evolution for positive chameleon potential energies.However conversely, the assumption of the chameleon density much more than the matter density can correspond to an accelerated-phase when the scalar field equation of state parameter is less than or equal to −1/3 (in this case, the equality corresponds to the deceleration parameter being zero, i.e. the transition point).This situation leads to the case that the kinetic energy of the chameleon scalar field must be less than or equal to half of its potential energy for positive chameleon potentials.
More clearly, for the described chameleon model in the matter-dominated phase, with the nonminimal coupling parameter β = 1, we have indicated that, when the kinetic energy of the chameleon scalar field is much more than the absolute value of, equal to, and equal to half of its potential energy8 (which respectively correspond to the scalar field equation of state parameter of approximately +1, 0, and −1/3), the matter density dominates at the beginning, but then, at the end of its dominance era, the chameleon scalar field density is dominant.Conversely, if the kinetic energy of the chameleon scalar field is much less than its potential energy for positive chameleon potential energies (i.e., when the scalar field equation of state parameter is approximately −1), the final result will still be the same as above, but for any non-minimal coupling parameter.That is, the result of this case does not depend on the non-minimal coupling parameter.Also, in this phase with the strong non-minimal coupling parameter β = 3.7 × 10 2 , the behavior of the matter density and the chameleon scalar field density, in the first three cases mentioned above, are both damped oscillations, i.e. the value of both densities decreases with increasing scale factor over time.
Furthermore, in the late-time accelerated-phase, first, the plausible assumption of the chameleon density much more than the matter density leads to the total equation of state parameter being almost identical to the scalar field equation of state parameter, the one described at the end of the second paragraph.Then, in the case that the kinetic energy of the chameleon scalar field is equal to half of its potential energy, with the strong non-minimal coupling parameter β = 3.7 × 10 2 , we have indicated that, although the matter density increases slightly, its amount is very negligible all along the path and the chameleon scalar field density, while decreasing, is dominant at all times.In this case with the non-minimal coupling parameter β = 1, the matter density and the chameleon scalar field density both decrease, and although the decrease of the matter density is less than that of the chameleon, the chameleon scalar field density dominates at all times.On the other hand, in the case that the kinetic energy of the chameleon scalar field is much less than its potential energy (i.e., when the scalar field equation of state parameter is approximately −1), the behavior of the matter density and the chameleon scalar field density do not depend on the non-minimal coupling parameter.In this case, we have shown that first the matter density dominates, then with the increase of the scale factor, the chameleon scalar field density, while always remaining constant, dominates the late-time phase.Also, we have taken a brief look at the chameleon scalar field profile while using the common power-law chameleon potential.
Finally, in both the matter-dominated phase and the late-time accelerated-phase, we have specified that, when the kinetic energy of the chameleon scalar field is much less than its potential energy (i.e., when the scalar field equation of state parameter is approximately −1), the behavior obtained for the presented chameleon model is similar to the ΛCDM model.Such compatibility suggests that the chameleon model is phenomenologically viable and can be tested with the observational data.