Constraining neutrino mass in dark energy dark matter interaction and comparison with 2018 Planck results

In this paper, we investigate the constraints on the total neutrino mass ∑mν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sum m_{\nu }$$\end{document} in a cosmological model in which dark energy and neutrinos are coupled such that the mass of the neutrinos and potentials are function of the scalar field as mν=m0exp(αϕmpl)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\nu }=m_{0}\exp (\frac{\alpha \phi }{m_{pl}})$$\end{document} and V(ϕ)=mpl4exp(-λϕmpl)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V(\phi )=m_{pl}^{4}\exp (\frac{-\lambda \phi }{m_{pl}})$$\end{document} respectively. The observational data used in this work include the type Ia supernovae (SN) observation (Pantheon compilation), CC, CMB and BAO data. We find that the neutrino mass is tightly constrained to ∑mν<0.125\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sum m_{\nu }< 0.125$$\end{document} eV 95% Confidence Level (C.L.) and the effective extra relativistic degrees of freedom to be Neff=2.955-0.12+0.11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_{eff}=2.955^{+0.11}_{-0.12}$$\end{document} 68% C.L in agreement with the Standard Model prediction Neff=3.046\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ N_{eff} = 3.046$$\end{document}, matter-radiation equality, zeq=338925+24\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_{eq}=3389^{+24}_{25}$$\end{document} (68% C.L). These results are in good agreement with the results of Planck 2018 where the limit of the total neutrino mass is ∑mν<0.12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sum m_{\nu }<0.12$$\end{document} eV (95% C.L., TT, TE, EE + lowE + lensing + BAO) , Neff=2.99-0.17+0.17\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_{eff}=2.99^{+0.17}_{-0.17}$$\end{document} (68% C.L., TT, TE, EE + lowE + lensing + BAO) and zeq=338721+21\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_{eq}=3387^{+21}_{21}$$\end{document} (68% C.L TT, TE, EE + lowE + lensing + BAO).


Introduction
The accelerating expansion of the universe [1][2][3][4] is one of the must surprising discoveries in cosmology. Also the observations of Cosmic Microwave Background (CMB) anisotropies indicate that the universe is flat and the total energy density is very close to the critical one [5].
According to General Relativity, the dynamic of the universe is dominated by a new (dark) energy form with negative pressure. There are prominent candidates for DE such as the cosmological constant [6,7] in which dark energy takes the form of a cosmological constant and dark matter is taken to be cold, in other words having an equation of state equal to zero. While CDM fits the available data very well, it suffers from a number of issues that motivate the study of altera e-mail: salehi.a@lu.ac.ir (corresponding author) natives. These include the fine-tuning [7] and coincidence [8] problems. In addition, there are certain tensions between early and late-universe observations in CDM. The presentday expansion rate of the universe, H 0 and the growth of structure, quantified by σ 8 , can be calculated using the bestfit CDM parameters to cosmological data, including the CMB. This gives rise to a smaller H 0 and a larger σ 8 than the results of local, late-universe measurements (for a recent discussion see Ref. [9]. A popular class of modifications to CDM is quintessence [10,11] , in which the cosmological constant is set to zero and a scalar field φ is introduced whose dynamical properties produce a negative equation of state giving rise to the observed late-time accelerated expansion of the universe. Normally it is assumed that the scalar field does not interact with dark matter. However there is no reason why this must be the case, and the consequences of relaxing this assumption have been widely studied. See Ref. [12] with references therein and [13][14][15][16][17] for a discussion of recent research on interacting dark energy. The other candidates are phantom (field with negative energy) [18] that explains the cosmic accelerating expansion. Meanwhile, the accelerating expansion of universe can also be obtained through modified gravity [19], brane cosmology and so on . On the other hand, to explain the early and late time acceleration of the universe. It is most often the case that such fields interact with matter; directly due to a matter Lagrangian coupling, indirectly through a coupling to the Ricci scalar or as the result of quantum loop corrections [44][45][46][47][48]. If the scalar field self-interactions are negligible, then the experimental bounds on such a field are very strong; requiring it to either couple to matter much more weakly than gravity does, or to be very heavy [49][50][51][52]. Unfortunately, such a scalar field is usually very light and its coupling to matter should be tuned to extremely to small values in order not to be a conflict with the Equivalence Principal [53]. The discovery of the accelerated expansion of the universe is also the main challenge for particle physics [54]. It requires new physics for the explanation of dark energy. Neutrinos were first shown to have mass in observations of neutrino flavour oscillations [55,56], the presence of which demands that at least two of the neutrino states are massive [57]. While the attempts of the laboratory experiments of particle physics to measure the absolute masses of neutrinos, have always been facing great challenges [58][59][60][61][62][63][64][65][66][67][68], the cosmological observations are more prone to be capable of measuring the absolute masses of neutrinos [65][66][67], since massive neutrinos can leave rich signatures on the cosmic microwave background (CMB) anisotropies and the large-scale structure (LSS) formation at different epochs of the cosmic evolution [68]. Recently some studies have attempted to constrain the total neutrino mass m ν and as well as the effective number of relativistic degrees of freedom (N e f f ) using cosmological observations  Also the cosmological consequences of interacting dark energy and dark matter have been widely studied , however, we are interested in consideration that the role of neutrino is prominent. In this work, we implement a cosmological model which proposed by [167] and extend by [168] to constrain total neutrino mass with observations. In this model dark energy and neutrinos are coupled such that the mass of the neutrinos is a function of the scalar field m ν = m 0 exp( αφ m pl ). The scalar field plays the role of dark energy and drives the late time accelerated expansion of the universe. The motivation of such consideration has been investigated by [168][169][170][171][172][173][174][175][176][177].
Here we consider a generalized model of [167] which allows both dark matter and neutrino interact with dark energy with different interacting couplings β and α. We constrain on ( m ν , N e f f , h, m , b h 2 , c h 2 , α, β, λ, ω) using observational data data. The structure of the paper is as follows. Section 2 introduces the cosmological model explored here while Sect. 3 describes the methodology and the measurements exploited in our data analyses. Section 4 presents our results and we conclude the article in Sect. 5

The model
The expansion rate of the Universe is given by Hubble parameter H =ȧ a , where a is scale factor and an overdot denote cosmic time derivative. We assume a spatially-flat Friedmann-Robertson-Walker Universe filled with baryons (b), radiation (r ), dark energy (ϕ), dark matter (dm) and neutrinos (ν). The baryons and radiation are regarded as noninteraction fields. The Friedmann equations which follow from Einstein field equations are as Where, ρ = ρ b + ρ r + ρ dm + ρ υ + ρ ϕ and p = p b + p r + p dm + p υ + p ϕ . From the conservation of the energy momentum tensor it follows the evolution equation for the total energy density: The baryons are treated like dust ( p b = 0) and the barotropic equation of state for the radiation field is ( p r = 1 3 ρ r ). Once both of them have no interaction with other components, the evolution equations for their energy densities are : respectively. Since the dark energy is modeled as a scalar field φ its energy density and pressure are given by where V (ϕ) denotes the potential of the scalar field. In this paper, we consider the interactions between dark matter and dark energy as the following evolution equatioṅ where β stands for the coupling constant between dark matter and dark energy. The neutrinos are understood to be massless and relativistic particles in the past, but with the coupling to the dark energy, they have acquired mass and became nonrelativistic, having an oscillating mass behavior at low redshifts. In this paper we follow the idea that proposed by [167] and extended by [168]. In the cosmological context, neutrinos cannot be described as fluid. Instead, one must solve the distribution function f (x i , p i , τ ) in phase space (where τ is the conformal time). Considering the case that neutrinos are collisionless, the distribution function f does not depend explicitly on time. Solving the Boltzmann equation, one can then calculate the energy density stored in neutrinos ( f 0 is the background neutrino distribution function): The evolution equation for its energy density according to [168] where α denotes coupling constant which can be related neutrino massm ν with relation α = d ln m ν dϕ . Furthermore, from the resulting equations it is possible to obtain evolution equations for the scalar field as For most of the Universe's history, the neutrinos are highly relativistic and (ρ ν − 3 p ν ) ≈ 0 such that the scalar field and the neutrinos are effectively uncoupled, here only the coupling parameter β is important . After the neutrinos become non-relativistic hence the coupling parameter α also becomes important. We consider an exponential potential where λ is a dimensionless parameter that determines the slope of the potential. The motivation for choosing these functions have been investigated in [169]. Also we define ω = P ν ρ ν . In order to simplify the field equations, we introduce following new variables, Hence, the equations of the autonomous dynamical system can be derived as, Where, N = ln a. In term of the new dynamical variable, we also have, In term of new variable the Friedmann equation (1) puts a constraint on new variables as We demonstrate that for the flat Friedmann-Robertson-Walker model the dynamics can be reduced to the form of the six dimensional autonomous dynamical system where by exerting the constraint (14) it reduces to five dimensional dynamical system. The parameters α, β, λ, ω are the free parameters of the model. While the main advantage of the dynamical system methods is that without knowledge of an exact solution it is possible to investigate the properties of the solutions as well as their stability, But this method is not necessarily to check the stability of the system. Rather, even if our goal is to solve equations numerically, this method is a useful method. Because in the solution of the original equations (Friedman and field equations), we are faced with differential equations of order 2 and higher, which are relatively more complex where not only the initial conditions but also the first and second order initial conditions must be determined for any dynamical variable. For example, because of advent of theφ in field equation, we need ϕ(0),φ(0) and ϕ(0) and due toḢ we need (a(0), H (0) andḢ (0) for numerical solutions, however, when the equations are introduced in terms of a set of the first order equations, only the initial condition for the new variable need to be determined x 1 (0)..x 5 (0), this makes numerical solution of the equations easier. On the other hand the variables in the dynamical system usually are dimensional lees and in most cases, or at list in the case of Friedman equation we have Pre information about the range of initial conditions. For example, we get Hence in numerical solution, it is not necessary to cover a large area of m (0) but we focus on a small area , approximately between 0.2 and 0.4. It makes the analysis easier. The importance of this issue becomes more reveal in observational cosmology where the initial condition play important role in the evolution of the universe However, the question that may arise is that why we don't implement this method for any set of high order differential equations. The answer is that although any second-order differential equation is equivalent to two first-order equations, however the suitable choice of new variables and converting equations to the first order may not be easy except in special cases, but if we could do that then we are faced with a set of first order equations which are more easer to solve.

Observational data, analysis and results
In what follows, first, we briefly describe the observational data sets used to constrain the parameters of the models under consideration.
• Pantheon: The use of type Ia supernovae (SNe) as standard candles has been of critical importance to cosmology, leading to the discovery of cosmic acceleration [54,178]. In this paper, we use the new "Pantheon" sample of Scolnic et al. [179], which adds 276 supernovae from the Pan-STARRS1 Medium Deep Survey at 0.03 < z < 0.65 and various low-redshift and HST samples to give a total of 1048 supernovae spanning the redshift range 0.01 < z < 2.3 The luminosity distance d L can be calculated by In order to incorporate the Eq. (15) with the dynamical system equations of (12), it can be rewriten in terms of the following differential equations   (12) must be coupled and solved simultaneously. Hence the distance modulus also can be obtained as μ th (z) = 5 log(x d ) + 42.38. we can compute the χ 2 -statistics for each case. Therefore, we proceed to define the following quantities: where N is the number of data points, σ i is the uncertainty associated with each measurement. • Cosmic Microwave Background(CMB): The observations of temperature anisotropies in the CMB provide a valuable independent test for the reality of dark energy at the recombination epoch z 1090. The photons were coupled to baryons and electrons before that red shift and decoupled right after. Due to the fact that in the Boltzmann and Einstein equations all the components of the universe are coupled, in order to extract information from the full spectrum, demanding numerical simulations are needed. A convenient and efficient way to summarize information from the CMB data, without using the full spectrum, is by employing the so called CMB shift parameters or distance priors. The CMB shift parameter R, given by [180,181] where E(z) = H (z) H 0 andz rec is the redshift of recombination z rec = 1090 [182]. The parameter R ties up the angular diameter distance to the last scattering surface, the comoving size of the sound horizon at z = 1091.3 and the angular scale of the first acoustic peak in CMB power spectrum of temperature fluctuations [180,181]. The updated value of R from WMAP5 is R obs = 1.710 ± 0.019 [183]. The χ 2 C M B for the CMB data is where the corresponding 1σ errors is σ R = 0.019.

• Baryon acoustic oscillations BAO data
For BAO data, from the measurement of the BAO peak in the distribution of SDSS luminous red galaxies, we define parameter A as [184] where z b = 0.35. The SDSS BAO measurement [184] gives A obs = 0.469(n s /0.98).0.35 ± 0.017, where the scalar spectral index is taken to be n s = 0.965as measured by Planck 2018 [185]. The parameter A is nearly model-independent and imposes the robust constraint as complement to SNIa data. The χ 2 for the BAO data is where the corresponding 1σ errors is σ A = 0.017. Table  1 shows the best fitted model parameters and initial conditions in both power law and exponential cases. • CC We use the cosmic chronometers (CC) data set comprising of 36 measurements spanning the redshift range z ≤ 2.36, recently compiled in [186] In order to put constraints on the parameters of the model we must note that the model has five independent variables (x 1 , x 2 , x 3 , x 4 , x 5 ) which according to 11 are equivalent to ( b , ν , r , dm , 1/2 φ ) as well as four free parameters of the model (α, β, λ, ω). Hence in order to solve the equation numerically the five initial conditions (x 1 (0), .., x 5 (0)) and value of the parameters must be known. In observational measurements one or more parameter are added to the free parameters. For example in numerical analysis using Pantheon data the two new variables x d = d L and x h = H and for CC, CMB and BAO data the variable x h = H are added to the free parameters. The other parameters are expressed in terms of the main parameters and can be constrained indirectly. For example m ν can be related to the main parameters (h, ν ) as where h is the reduced Hubble constant (the Hubble constant H 0 = 100h km/s/Mpc. Hence if the parameters (h, ν ) are constrained then the parameter m ν is constrained automatically.
The relativistic energy density in the early universe include the contributions from photons and neutrinos, and possibly other extra relativistic degrees of freedom, called dark radiation. The effective number of relativistic species, including neutrinos and any other dark radiation, is defined by a parameter, N e f f , for which the standard value is 3.046 corresponding to the case with threegeneration neutrinos and no extra dark radiation [187]. If the value of N e f f is beyond 3.046, it indicates that there is some dark radiation other than three-generation active neutrinos. The behaviour of dark radiation is exactly equivalent to massless neutrinos. Thus, the total radiation energy density in the Universe is given by ρ r = ρ γ 1 + 7 8 4 11 4 3 N e f f (26) where ρ γ is the energy density of photons. We parametrize the relativistic degrees of freedom using the effective number of neutrino species, N e f f . This quantity can be written in terms of the matter density, m h 2 , and the redshift of matter-radiation equality z eq as [182] N e f f = 3.04 + 7.44 m h 2 0.1308

Results
Throughout this section we will present the results obtained within the two different IDE scenario

IDE+ m ν
The results for the cosmological parameters within this interacting dark energy model are shown in Table 1. Figure 1 also show the parametric space at 68 %CL and 95%CL for some selected parameters for the different observational data sets. In order to compare our results with those obtained by Planck 2018 [185], we have listed some of the Planck 2018 results [185] in Table 2. For this case the free parameters are From the analyses of the Pantheon data alone, as shown in Table 1, we find that This result is very close to the result of [121], the case IDE+ m ν , using (CMB + Pantheon + CC data) with m ν < 0.255 eV at 95% CL. Both the model is similar to ours and the data used includes Pantheon data. The result is also closed to the result of [122], the case interacting vacuum scenario (IVS)+ m ν , with m ν < 0.277 eV at 95% CL) and the result of Planck 2018 [185] the case (TT, TE, EE + lowE + lensing) with m ν < 0.241 eV at 95% CL. Using CC data, we find that which is close to the result of [121], the case IDE+ m ν using (CMB+Pantheon+CC data) with m ν < 0.159 eV at 95 % CL and the result of Planck 2018 [185] the case (TT, lowE + BAO) with m ν < 0.16 eV at 95% CL.
Using CMB+BAO, we find that This result is very close to the result of [121], the case IDE+ m ν , using (CMB data) with m ν < 0.313 eV 95% CL and comparable with and comparable with results of [185] the case TT, TE, EE + lowE[CamSpec] with m ν < 0.38 eV 95% CL Fig. 1 The constraints at the 68% and 95% CL two-dimensional contours for selected cosmological parameters and parameters of the model parameters in IDE+ m ν scenario for the Pantheon, CC, CMB + BAO and Pantheon + CC + CMB + B AO dataset For combination of full data, Pantheon + CC + CMB + BAO, we find The result is very close to results of [185] TT, TE, EE + lowE + lensing + BAO with m ν < 0.12 eV at 95% CL and case TT, TE, EE + lowE + BAO with m ν < 0.13 eV at 95% CL. It is also comparable with that obtained by [121] with m ν < 0.156 eV at 95% CL using same model IDE+ m ν and same dada (Pantheon + CC + CMB + BAO) We also put constraint on coupling parameters (α, β, λ, ω). The lower panel of Fig. 1 shows (68.3%, 95.%) confidence levels for the parameters (α, β) and (λ, ω) for Pantheon, CC and CMB + BAO and combination of the data Pantheon + CC + CMB + BAO. The results also have been listed in Table 1 m pl ), the positive value of λ indicates that potential V (ϕ) is a monotonically decreasing function of ϕ. and has a negative gradient. It is also important to note that the best fitted values of λ for both individual and combination data values are very close to that deter-   mined by upper bounds on early dark energy ( λ ≥ 10) [161]. In fact, when the potential energy approaches a constant , where the evolution of the cosmon field stops close to a value φ t which is characteristic for the transition between the two different cosmological epochs, it acts similar to a cosmological constant and causes the accelerated expansion and for λφ t m pl 276 the cosmological constant has a value compatible with observation. This amount gives λ ≥ 10 (for more discussion see [161]) For most of the Universe's history, the neutrinos are highly relativistic and (ρ ν − 3 p ν ) ≈ 0 such that the scalar field and the neutrinos are effectively uncoupled, here only the coupling parameter β is important. After the neutrinos become non-relativistic hence the coupling parameter α also becomes important. Constraining on parameter α, we have obtained α = 15.5 +9.50 −9.50 , α = 7.01 +10.05 −10.05 , α = 20.02 +10.20 −10.20 and λ = 11.1 +3.60 −3.60 at 68% CL for Pantheon, CC, CMB + BAO and Pantheon + CC + CMB + BAO respectively. As investigated by [188], the following conditions must be met to give rise to growing neutrino quintessence: • V (ϕ) must have a negative gradient in order to cause the value of the scalar field to increase with time. This gradient must be sufficiently steep that ϕ reaches large enough values in the late Universe to act as dark energy. • |α| must be sufficiently large when the neutrinos become non-relativistic that β(ρ ν − 3 p ν ) is able to act as a strong enough restoring force to stop the evolution of ϕ in Eq. (10).
The best fitted of (α, λ) satisfy the above condition. Despite the small value of ν the neutrinos are important for the evolution of the cosmon due to their large coupling α For both individual and combination of the dataset, we find that the large value for coupling parameter β. We find |β| > 19| at (95% CL. This indicates there is a strong interaction between dark matter and dark energy.

IDE+ m ν + N e f f
The results for the cosmological parameters within this interacting dark energy model are shown in Table 3 Fig. 2 The constraints at the 68% and 95% CL two-dimensional contours for selected cosmological parameters and parameters of the model in IDE+ m ν + N ef f scenario for the Pantheon, CC, CMB + BAO and Pantheon + CC + CMB + BAO dataset however the bounds on some of the parameters have been changed slightly.
From the analyses of the Pantheon data alone, as shown in Table 3, we find that Which in comparison of previous model, IDE+ m ν , with m ν < 0.253 eV have been changed significantly, however it is close to the result of [189] where using same model, interacting scenario IDE1p + m ν + N e f f and using Planck 2018 data, have been obtained m ν < 0.438 eV at 95% CL The interesting result of our analysis is that for this model, the most stringent upper limit we have on this parameter is obtained for the data set combination C M B + B AO + Pantheon + CC in which m ν < 0.125 eV 95%C L (33) This result is in good agreement with the results of Planck 2018 [185], where the limit of the total neutrino mass is m ν < 0.12 eV at (95% C.L. using TT, TE, EE + lowE + lensing + BAO data) and close to the result of [121], where they also find the most stringent upper limit on this parameter for the same model, IDE+ m ν and the same data (CMB + Pantheon + CC data) with m ν < 0.15 eV at 95% CL. For combination data we also find which is very close to the result of Planck 2018 [185] with N e f f = 2.96 +0.34 −0.33 at 68%CL, the case TT,TE,EE,LowE +lensing+BAO and close to result of [121] with N e f f = 3.02 +0. 34 −0.33 using the same data and same model. For this case, the best fitted values for λ, have obtained as λ = 9.7 +0.
The results are close to those obtained in previous case, IDE+ m ν . Hence same as the previous case, the positive values of λ indicate that potential V (ϕ) is a monotonically decreasing function of ϕ. and has a negative gradient as well as these value are consistence with that predicted by observation.

Conclusion
In this paper, we have explored possible extensions of the Interacting Dark Energy, where the dark energy and the dark matter fluids interact with each other.
We have considered a cosmological model in which dark energy and neutrinos are coupled such that the mass of the neutrinos and potentials are function of the scalar field as m ν = m 0 exp( αφ m pl ) and V (φ) = m 4 pl exp( λφ m pl ) respectively. While the theoretical aspects of the neutrino-dark energy interaction have been studied in previous studies [167][168][169]188], here we have extended the model such that not only neutrino but also dark matter interact with dark energy with different coupling terms Q de−dm = β m pφ ρ dm and Q de−ν = α m p (ρ ν − 3 p ν )φ and focused on observational aspect of the model. We have exploited the most recent publicly available cosmological observations, which include the Supernovae Type Ia Pantheon data and measurements of the Hubble parameter from Cosmic Chronometers. CMB data, Baryon Acoustic Oscillations data (BAO) to put constrain on parameters on m ν , N e f f , z eq , m , c , b , h , α, β, λ, ω We find that while the results of individuals and combination of datasets are closed, the most stringent upper limit we have on this parameter is obtained for the data set combination C M B + B AO + Pantheon + CC in which m ν < 0.125 eV 95%C L This result is in good agreement with the results of Planck 2018 [185], where the limit of the total neutrino mass is m ν < 0.12 eV at (95% C.L. using TT,TE,EE+lowE+lensing+BAO data) and close to the result of [121], where they also find the most stringent upper limit on this parameter ..for the same model, IDE+ m ν and same data (CMB + Pantheon + CC data) with m ν < 0.15 eV at 95% CL.
This study also investigated one of the main problem in growing neutrinos and cosmological selection. As point out by Amendola and Wetterich [190], the most crucial observational issues can be understood by understanding on constant parameters (λ, α, β) it will be a challenge to measure them or to falsify the growing matter scenario. For neutrino growing matter a determination of λ and α would fix the neutrino mass, allowing for an independent test of this hypothesis by comparing with laboratory experiments. The values obtained for these parameters in this study are close to those require for growing matter mechanism.
For both model IDE + m ν and IDE+ m ν + N e f f and for all dataset, Both individual and combination, we find the mean value of λ as λ 10. This value is very close to that determined by upper bounds on early dark energy (λ ≥ 10) [161]. As also point out by V (ϕ) must have a negative gradient in order to cause the value of the scalar field to increase with time. This gradient must be sufficiently steep that ϕ reaches large enough values in the late Universe to act as dark energy. Since we get V (φ) = m 4 pl exp( −λφ m pl ), the best fitted values of λ which is positive satisfies this condition.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The paper has no external data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .