Viable Intermediate Inflation in the Mimetic DBI Model

We study the intermediate inflation in the mimetic Dirac-Born-Infeld model. By considering the scale factor as $a=a_{0}\exp(bt^{\beta})$, we show that in some ranges of the intermediate parameters $b$ and $\beta$, the model is free of the ghost and gradient instabilities. We study the scalar spectral index, tensor spectral index, and the tensor-to-scalar ratio in this model and compare the results with Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data at $68\%$ and $95\%$ CL. In this regard, we find some constraints on the intermediate parameters that lead to the observationally viable values of the perturbation parameters. We also seek the non-gaussian features of the primordial perturbations in the equilateral configuration. By performing the numerical analysis on the nonlinearity parameter in this configuration, we show that the amplitude of the non-gaussianity in the intermediate mimetic DBI model is predicted to be in the range $-16.7<f^{equil}<-12.5$. We show that, with $0<b\leq 10$ and $0.345<\beta<0.387$, we have an instabilities-free intermediate mimetic DBI model that gives the observationally viable perturbation and non-gaussianity parameters.


Introduction
A new approach to General Relativity has been proposed by Chamseddine and Mukhanov in 2003, called mimetic gravity. The important property of their approach is that the conformal symmetry is respected as an internal degree of freedom [1]. Chamseddine and Mukhanov in their interesting proposal, have written the physical metric in terms of an auxiliary metric and a scalar field. In fact, the metric is given by where the free non-dynamical scalar field φ encodes the conformal mode of the gravity. Also, the definition (1) shows that if we perform a Weyl transformation on the auxiliary metric (g µν ), the physical metric (g µν ) remains invariant. Using the equation (1) gives the following constraint on the scalar field [1] g µν φ ,µ φ ,ν = −1 .
In the action of the mimetic gravity introduced in Ref. [1], there is a contribution of the matter fields coupled to g µν that leads to an extra term in the Einstein's field equations. In the sense that the dependence of this extra term to the scale factor is a −3 , it is considered as a source of dark matter. The mimetic gravity scenario has been explored in another mathematical approach, by adding the Lagrange multipliers in the action of the theory [2,3]. In Ref. [4], some ghost-free models of the mimetic gravity have been discussed. The authors of Ref. [5] have considered a potential for the mimetic field in a Lagrange multiplier approach, leading to some interesting results. In fact, they have shown that if we take the appropriate potential terms, it is possible to consider the mimetic field as inflaton, quintessence, or phantom fields. The mimetic gravity has attracted a lot of attention and the authors have extended it to the braneworld scenario [6], non-minimal coupling model [7,8], f (G) gravity [9], Horndeski gravity [10,11], f (R) theories [12,14,13], unimodular f (R) gravity [15] and Galileon gravity [16]. In studying the perturbations in mimetic gravity models, it is necessary to have the ghost and gradient instabilities-free models. In Ref. [17] it is shown that the direct coupling between the curvature of the space-time and the higher derivatives of the mimetic field can help to overcome such instabilities in some ranges of the parameters space. See the papers in Refs. [18,19,20,21,22,23,24,25,26,27,28,29,30] for more works on the (in)stability issue.
On the other hand, to solve some problems of the standard model of cosmology, the inflation paradigm has been introduced, where a single canonical scalar field (inflaton) with a flat potential has been considered. The flat potential causes the slow-roll of the inflaton and enough exponential expansion of the early universe. In this simple model, the dominant modes of the primordial perturbations are predicted to be scale-invariant, adiabatic and gaussian [31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46]. However, there is a lot of attention to the extended inflation models predicting the non-Gaussian distributed perturbations [46,47,48,49,50,51,52,53,54,55].
One interesting inflation model is the DBI (Dirac-Born-Infeld) inflation [56,57,58,59,60,61]. In the string based DBI model, where the radial position of a D3 brane characterizes the scalar field [56,58], it is possible to have large non-Gaussianity that many authors are interested to [58,59,60,61,62,63,64,65,66,67]. Another string-based scalar field is the tachyon field [68,69] which is an interesting scalar field in the inflation models [61,70,71]. We have included the DBI and tachyon field in the mimetic models in our previous works. In Ref. [72], we have assumed that the scalar field of the DBI model is a mimetic scalar field. By considering the Lagrange multiplier approach and assuming the power-law inflation, we have studied the inflation and perturbations in this model and shown that this model in some ranges of the model's parameter space is free of the ghost and gradient instabilities. We have also shown that, in those ranges of the model's parameter space, the values of the perturbations parameters are observationally viable. In fact, with the Mimetic DBI (MDBI) model, there is no need to consider such complicated higher-order terms in the action of mimetic gravity. In Ref. [73], we have considered the tachyon field in the mimetic gravity setup and within the Lagrange multiplier approach. We have adopted both powerlaw and intermediate scale factors and studied the inflation in the tachyon mimetic model. We have shown that, in both cases, the mimetic model is free of instabilities in some ranges of the model's parameters. Also, those ranges of the parameters that lead to instabilities free tachyon mimetic model, give observationally viable perturbation parameters.
In the continuation of the previous works, in this paper, we study other aspects of the MDBI model. Now, we consider the MDBI model with an intermediate scale factor [74,75,76] that gives another type of potential. In this regard, we show that the intermediate MDBI model also is an instabilities-free mimetic model that is consistent with observational data. Therefore, with this new work on the MDBI model, we prove that to have an instabilities-free mimetic model there is no need to restrict ourselves to one type of potential (or scale factor). The MDBI model is a simple mimetic model that gives a viable cosmological model with both power-law and intermediate scale factor (corresponding to two types of potential). This makes the MDBI model more favorable. Also, in this paper, we study the non-gaussian features of the primordial perturbations that is one of the important issues in the inflationary models. The prediction of the intermediate MDBI model for the sound speed of the perturbations and also the amplitude of the perturbations in equilateral configuration is the issue that hasn't been studied in the MDBI model previously.
This paper is organized as follows: In section 2, we review the MDBI model in the Lagrange multiplier approach. In this section, we present the perturbation parameters (the scalar spectral index, the tensor spectral index, and the tensor-to-scalar ratio) that are functions of the Lagrange multiplier. The amplitude of the non-gaussianity in the equilateral configuration is also presented in this section, which is related to the Lagrange multiplier via the sound speed of the model. In section 3, we introduce the intermediate MDBI model where the scale factor is given by a = a 0 exp(bt β ). By this scale factor, we find the Hubble parameter in terms of the intermediate parameters β and b. Then, by obtaining the potential and the Lagrange multiplier in terms of the Hubble parameter, we find the slow-roll parameters in the intermediate MDBI model. Also in this section, we study the model numerically and show that the intermediate MDBI model in some ranges of its parameter space is free of the ghost and gradient instabilities. In section 4, by performing numerical analysis, we compare the results with Planck2018 observational data. In this regard, we show that the intermediate MDBI model for some values of β and b, which lead to the instabilities-free model, gives observationally viable values of the perturbation parameters. We also present some predictions of the model on the non-gaussian feature of the primordial perturbations. In section 5, we present a summary of our work.

Mimetic DBI Model
The action of the DBI mimetic gravity, in the presence of the Lagrange multiplier and a potential term, is given by where R is the Ricci scalar, V (φ) presents the potential of the scalar field, F −1 (φ) is the inverse brane tension (note that, the D3-brane passes a compact manifold which geometry of its throat is related to the F −1 (φ)). Also, κ is the gravitational constant, defined as κ 2 = 8πG c 4 . The parameter α is a coupling parameter which is corresponding to the string theory parameter. Although we adopt a DBI-like Lagrangian in our work, the scalar field φ is not a DBI field. It is a mimetic field obeying the constraint (2). This means that its dimension is [T ] = [M ] −1 , demonstrating the dimension of time (note that we work in natural units where the speed of light is c = 1 and also = 1). Since the action should be dimensionless, to have dimensional consistency, there should be a coupling parameter α with dimension [M ] 4 = [T ] −4 in the Lagrangian term, as shown in equation (3). This parameter is corresponding to the string theory parameter presented in [56]. Also, λ is a Lagrange multiplier by which we enter the mimetic constraint (2) in the action.
If we vary action (3) with respect to the metric, we find the following Einstein's field equations (G µν = κ 2 T µν ) By using the flat FRW metric as the background, the field equations (4) give the following Friedmann equations in this model The equation of motion of the mimetic field φ in the MDBI model is obtained by varying the action Also, by varying action (3) with respect to λ we reach the constraint (2). In this paper, we are going to consider the MDBI setup as an inflation model. The slow-roll parameters in the inflation model are given by the following definitions where c s is the sound speed of the perturbations, defined as c 2 s = P ,X ρ ,X (the subscript ", X" shows derivative with respect to X = − 1 2 ∂ ν φ ∂ ν φ). In our model, the square of sound speed is given by which should satisfy the constraint 0 < c 2 s ≤ c 2 (note that c is the local speed of light). In the case of c ≡ 1, the constraint becomes as 0 < c 2 s ≤ 1 [77,78]. To seek the observational viability of the MDBI model, it is important to compare the values of the perturbation parameters with Planck2018 data [79,80]. The Planck collaboration has considered that in the case of the statistical isotropy, the two-point correlations of the CMB anisotropies are described by the angular power spectra [81,82,83,84,85]. In fact, the following expressions for the contributions from the scalar and tensor perturbations in the CMB angular power spectra have been used [86] C ab,s where a, b = T, E, B. Also, l is the multipole moment number, ∆ s l,A and ∆ T l,A are the transfer functions 2 , and A j (k) (j = d, T ) is the primordial power spectrum, identified by the physics of the primordial universe [86]. In one procedure to compare the inflationary parameters with data, the Planck collaboration has expanded the scalar and tensor power spectra in a model-independent form as [80,86] where, A j is introduced as the amplitude of the scalar (for j = s) or tensor (for j = T ) perturbations. Also, dn j d ln k is the the running of the scalar (for j = s) or tensor (for j = T ) spectral index and d 2 ns is the running of the running of the scalar spectral index. The ratio between the tensor and scalar amplitudes is an important perturbation parameter, named tensor-to-scalar ratio.
To perform numerical analysis and compare the results with Planck2018 data, we should obtain the perturbation parameters in the MDBI model. The scale dependence of the scalar spectral index, at the time of sound horizon exit of the physical scales (c s k = aH), is identified by Calculating the perturbation parameters at pivot scale k = k * causes that the running term doesn't appear in the definition (16). For the MDBI model, the amplitude of the scalar spectral index is given by [72] where Positive values of W s makes the model free of the ghost instability. Now, we find the scalar spectral index as follows which is expressed in terms of the slow-roll parameters. Equation (14) gives the following tensor spectral index The parameter A T , the amplitude of the tensor perturbations, is given by leading to We can find the tensor-to-scalar ratio from equations (15) and (22) as follows or Equation (24) is an important equation, named the consistency relation, which in the simple single filed inflation with a canonical scalar field simplifies to r = 16 .
To get more information about the viability of an inflation model, it is useful to study the nongaussian feature of the primordial perturbations. Although the two-point correlation characterizes the gaussian perturbations, to get the additional statistical information related to the non-gaussian distribution, we should consider three and higher-order correlations. In the interaction picture, the 3-point correlation for the spatial curvature perturbation Ψ is given by [46,50,89] where and A s is given by (17). The parameter G Ψ in equation (26) is defined as where the shapes of the non-gaussianity are given by and Equation (26) shows that the three-point correlator depends on the three momenta k 1 , k 2 and k 3 . Note that these momenta should satisfy the translation and rotational invariance. By defining the following dimensionless parameter, called "nonlinearity parameter", we can measure the amplitude of the non-gaussianity. The nonlinearity parameter depends on the shape of the non-gaussianity. Different values of the momenta lead to different shapes and there is a maximal signal for each shape in a special configuration of the three momenta.
However, as has been demonstrated in Refs. [57,90,91,92], for the k-inflation and higher-order derivative models, the maximal peak of the signal occurs at the equilateral configuration, where k 1 = k 2 = k 3 . In this regard, here also, we obtain the parameter G Ψ in this configuration and at the leading-order as [57,58,59] G equil leading to the following nonlinearity parameter This equation is very important to explore the non-gaussian feature of the primordial perturbations in an inflation model and compare with observational data.
In the following, we study the intermediate inflation in the MDBI model and perform some numerical analysis on this model.

Intermediate MDBI Model
One of the interesting scenarios in inflation models is intermediate inflation. In the intermediate inflation, the scale factor evolves faster than the power-law inflation (a = t p ) but slower than the standard de Sitter inflation (a = exp(Ht)) [74,75,76]. In fact, in the intermediate inflation the evolution of the scale factor is given by where b is a constant and 0 < β < 1. This scale factor leads to the following Hubble parameter To obtain the main perturbation parameters in the intermediate inflation, we should follow [93,94] and find the potential in terms of the Hubble parameter and its derivatives. From now on, we assume F −1 (φ) = V (φ). Note that, in Ref. [58], it has been demonstrated that for AdS 5 throat, the warp factor F −1 (φ) is equal to ς φ 4 . Also, in the case with AdS 5 × X geometry, the potential of a DBI field is quartic. Another interesting case is clarified in Refs. [95,96]. In those papers, the authors have shown that with F ∼ e mφ and V ∼ e −mφ (with m to be a constant), we can get the Lagrangian of the DBI model. Following Refs. [95,96], we consider F −1 (φ) = V (φ) in our calculations. In fact, it is always possible to have viable DBI inflation with these choices of functions. However, in this paper, we construct the potential (and therefore, F) for intermediate inflation from background equations. In this way, there is no need to choose an arbitrary function for V and F. Now, we introduce a new scalar field ϕ. This field is identified by the number of e-folds N and parameterizes the scalar field φ as φ = φ(ϕ). From these points and by using equation (7), we obtain the potential in the intermediate inflation as follows Note that, in the above and forthcoming equations, a prime refers to the derivative of the parameter with respect to N , and also we have H = H(N ). From equations (6) and (37) we get where The slow-roll parameters in terms of the Hubble parameter and its derivatives are given by the following expressions where H ≡ H(N ) and N = 2HH + 3H 2 with In the intermediate inflation, the slow-roll parameter takes the following form   The slow-roll parameters η and s are very long and complicated, so we avoid writing these parameters here. By using the above equations, we can express the perturbation parameters in terms of the model's parameters. In Ref. [72], we have shown that the MDBI model with powerlaw scale factor is free of ghost and gradient instabilities. Now, in this paper, we show numerically that the intermediate MDBI model also, in some ranges of the model's parameter space, is free of the instabilities. In this regard, by using equations (37) and (38) and by assuming F = V −1 , we can obtain the sound speed (equation (10)) in the intermediate MDBI model. Any range of the parameter space in which we have c 2 s > 0, leads to the gradient instability-free intermediate MDBI model. Note that, another constraint on the sound speed is c s ≤ c, where c is the value of the local speed of light. This constraint is required from causality. If we consider the scale factor (35) and perform some numerical analysis, we find that the range of the parameter space lead to 0 < c 2 s ≤ 1 (with c ≡ 1). The result is shown in figure 1. Note that, in this figure and forthcoming figures we have adopted κ = 1 and α = 1. The pink region in the left panel of figure 1 shows the ranges of the parameters β and b leading to the gradient instability-free intermediate MDBI model. Now, we study W s to see if the intermediate MDBI model is free of ghost instability. To this end, we use equations (18), (37), (38) and (39). By performing the numerical analysis, we find the result shown in the right panel of figure 1. In summary, the intermediate MDBI model in some ranges of its parameter space is free of gradient and ghost instabilities, which is a good result.
To ensure there are no instabilities during the whole inflationary evolution, we plot parameters c 2 s and W s versus the e-folds number N , for some sample values of the model's parameter. The results are shown in figure 2. This figure has been plotted with b = 10 and 0.340 < β < 0.400 (we see in the next section that, these adopted values of b and β are observationally viable). As this figure shows, both conditions 0 < c 2 s ≤ 1 and W s > 0 are satisfied during inflation. Now, we study the behavior of the potential (37) versus the scalar field. To this end, we consider  (1)) and from equations (36) and (37) we obtain By using this equation, we can perform a numerical study on the evolution of the potential versus the scalar field. The result is shown in figure 3. At the initial times, where inflation happens, the potential is large, and also the friction term 3Hφ. It is worth checking if inflation ends in this setup. The inflation ends when one of the slow-roll parameters reaches unity. In this regard, we study the first slow-roll parameter versus the e-folds number N . The result is shown in figure 4. As the figure shows, in this model inflation ends after about 60 e-folds (which is also corresponding to the minimum of the potential) and the graceful exit of inflation towards the matter-dominated era can be achieved. Note that, since the scale factor in the intermediate inflation evolves faster than the one in the power-law inflation and slower than the one in the exponential inflation, the same situation happens for the slow-roll parameter . In this model also, it is possible to have the seeds for the observed dark matter. In fact, by considering the conservation of the energy-momentum tensor (as ∇ µ T µν = 0), we obtain By considering the constraint (2) and H 2 = κ 2 3 ρ, form equations (6) and (46) we have After the end of the inflation, at the moment we reach the minimum of the potential, we have V = F −1 = constant. At that point, the slope of the potential is zero and so is the right-hand side of the equation (46). In this case, we have and therefore where C is a constant with dimension [M ] 4 . In this way, the seeds of the dark matter are obtained and the term C a 3 determines the amount of dark matter in the intermediate MDBI model. Note that, although in the inflation era (where the right-hand side of the equation (47) is not zero) it is possible to have a term like C a 3 , this term is diluted away quickly. However, after the end of inflation, this term has an important role. Also, by considering the gravitational particle production at the end of inflation, it is possible to get the observed radiation and baryons in the universe [97]. Although the potential is a function of time and increases again, the gravitational created particles dominate the potential term near the minimum. Then, by increasing the time, the potential term (and also, F −1 √ 1 − F term) dominates the particles and probably becomes the dark energy component leading to the late-time acceleration of the universe.
In the next section, we study the observational viability of this model with the Planck2018 data.

Comparing with the Planck2018 Observational Data
When an inflation model is constructed, it is important to check if its results are consistent with observational data. The observational data give some constraints on the perturbation parameters such as the scalar spectral index, the tensor spectral index, and the tensor-to-scalar ratio. Also, the observational data sets constraints on the amplitudes of the non-gaussianity in the equilateral configuration. Therefore, by studying these parameters in an inflation model and comparing the results with observational data, we can explore the viability of the model. The constraint on the scalar spectral index, from Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data, based on ΛCDM+r + dns d ln k model, is n s = 0.9658 ± 0.0038. This dataset gives the constraint on the tensorto-scalar ratio as r < 0.072. Also, Planck2018 TT, TE, EE +lowE+lensing+BK14+BAO+LIGO and Virgo2016 constraint on the tensor spectral index is as −0.62 < n T < 0.53. By using these data, we can obtain some constraints on the intermediate parameters β and b. By substituting equations (40), (41) and (43) for the intermediate inflation in the equations (19), (22), (23), we obtain the tensor-to-scalar ratio, the scalar spectral index and the tensor spectral index in terms of the parameters β and b. Then, we perform numerical analysis on these parameters and compare the results with observational data. Our numerical analysis shows that, for 0 < b ≤ 10, depending on the values of b, the scalar spectral index in the intermediate MDBI model is consistent with observational data if 0.345 < β < 0.387. This is shown in the left-upper panel of figure 5. The tensor-to-scalar ratio in this model is consistent with observational data if, depending on the values of b, 0.341 < β < 1. This is shown in the right-upper panel of figure 5. Also, we have found that the range 0.044 < β < 1 leads to the observationally viable values of the tensor spectral index in  figure 5. We have also studied the tensor-to-scalar ratio versus the scalar spectral index in the background of Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data at 68% CL and 95% CL. The results are shown in figure 6. From this figure, we have obtained some constraints on the parameter β that are summarized in table 1. The tensor-to-scalar ratio versus the tensor spectral index in the background of Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data (figure 7) is another studied case that gives some more constraints presented in table 1. Another important property in studying an inflation model is the non-gaussian feature of the primordial perturbations. As we mentioned earlier, the amplitude of the non-gaussianity is related to the sound speed of the perturbations and in this way, it is related to the model's parameter space. In this regard, studying the non-gaussian features of the perturbations in the inflation models helps us get more information about the models. We study the amplitude of the nongaussianity in equilateral configurations. By considering the constraints on the model's parameters obtained from the observational viability of the scalar and tensor spectral index and the tensor-toscalar ratio from Planck2018, we analyze the non-gaussianity in the intermediate MDBI inflation numerically. Here, we consider the non-gaussianity parameter f equil , which is related to the sound speed by equation (34). By using equations (10), (37) and (38) and also assuming F −1 = V , we relate the square of the sound speed of the primordial perturbations to the model's parameters. The square of the sound speed is also related to the tensor-to-scalar ratio via equation (23). In this regard, the constraints on r set some constraints on c 2 s that is shown in figure 8. The planck2018 combined temperature and polarization analysis gives the constraint on the sound speed of the DBI model as c DBI s ≥ 0.086, at 95% CL. Also, this data gives the constraint on the sound speed of the general P (X, φ) model (where, X = − 1 2 ∂ ν φ ∂ ν φ) as c s ≥ 0.021, at 95% CL. According to our numerical analysis, shown in figure      equation (34). The Planck2018 data gives the constraint on the equilateral configuration of the non-gaussianity as f equil = −26 ± 47. By using this constraint, we can find the range of the parameter space leading to the observationally viable values of the equilateral non-gaussianity in the intermediate MDBI model. The result is shown in the left panel of figure 9. Also, the right panels of figure 9 show the prediction of the model for the equilateral amplitude of the non-gaussianity, where we have used the observationally viable ranges of the parameters b and β which are obtained from the comparing of the tensor-to-scalar ratio and scalar spectral index with Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data at 95% CL. Note that, we have used the same ranges of the parameters to plot figures 6 and 7. However, plot 6 is r − n s behavior and plot 7 is r − n T behavior. From the slope of the r − n s plot, we see that it is possible to have larger values of r, but these larger values are corresponding to smaller values of n s , which are not observationally viable. Therefore, we haven't shown those parts of the plot which are not consistent with observational data. However, the situation is different for r − n T . As figure 7 shows, if we consider the larger values of r, the parameter n T is still consistent with observational data.
Up to here, we have numerically studied the perturbation parameters n s , n T and r and also the non-gaussianity parameter f equil N L , to find some constraints on the model's parameters. According to the obtained results, for every parameter, there are some ranges in the parameter space which make the model observationally viable. However, we are interested in the case where all parameters are consistent with Planck2018 data in the same range. By this, we mean that it is interesting to find a range for the parameters β and b where we have an instabilities-free and observationally viable intermediate MDB model. Our data analysis to obtain such ranges shows that by 0 < b ≤ 10 and 0.345 < β < 0.387, it is possible to have an instabilities-free intermediate MDBI model that gives the observationally viable perturbations. Also, with these ranges, the values of the equilateral amplitude of the non-gaussianity are consistent with observational data.  In Ref. [72], we have studied the power-law DBI and power-law MDBI models, where a = a 0 t n , with details. We have obtained the main inflation and perturbation parameters in both models and performed a numerical analysis on those parameters. According to our analysis in Ref. [72], the result of the numerical analysis on the perturbation parameters of the power-law DBI model is not consistent with Planck2018 observational data. Then, we have considered the power-law MDBI model. We have shown that the power-law MDBI model is an instabilities-free model. We have also shown that the scalar spectral index and the tensor-to-scalar ratio in the power-law MDBI model are consistent with Planck2018 TT, TE, EE+lowE +lensing data at 95% CL. However, these perturbation parameters in the power-law MDBI model are not consistent with Planck2018 TT, TE, EE+lowE+lensing+BK14+BAO data, where the combination of the BICEP2/Keck Array 2014 and Planck2018 data is considered. Note that, in the case of constant sound speed, it is possible to get observational consistency, however, our attention is on the varying sound speed case. In this paper, We have shown that if we consider an intermediate MDBI model, the perturbation parameters n s and r are consistent with Planck2018 TT, TE, EE+lowE+lensing+BK14+BAO data at both 68% CL and 95% CL. This is an interesting advantage of the intermediate MDBI model over the power-law MDBI model. We have also explored the non-gaussian feature of the primordial perturbation in the power-law DBI and power-law MDBI models, in Ref. [72]. In that paper, it has been shown that the amplitude of the primordial non-gaussianity in the power-law DBI model is too large to be consistent with Planck2018 observational data. Also, we have shown that the prediction of the power-law MDBI model for the amplitude of the equilateral non-gaussianity is very small (of the order of 10 −4 ). However, as we have seen in the current paper, the equilateral non-gaussianity in the intermediate MBI model is in the range −16.7 < f equil < −12.5. This is a result that is consistent with planck2018 data.
It seems that, by considering both perturbation and non-gaussianity parameters, the intermediate MDBI model is consistent with Planck2018 data and therefore is more favorable.

Summary and Conclusion
Recently, it has been shown that to have a ghost and gradient instabilities-free mimetic gravity model, one can consider a DBI-like term in the action of the mimetic gravity and adopt a power-law scale factor. In this paper, we have considered a MDBI model with intermediate scale factor as a = a 0 exp(bt β ). In this regard, we have studied the intermediate inflation in the MDBI model. We have shown that, with the intermediate MDBI model, it is possible to have a mimetic gravity model that is free of ghost and gradient instabilities in some ranges of the intermediate parameters b and β. This means that, in those ranges of the parameters, we have 0 < c 2 s ≤ 1 and W s > 0. To seek the observational viability of the models, we have studied the perturbation and nongaussianity parameters of this model and compare the results with observational data. From Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 data, we have the value of the scalar spectral index as n s = 0.9658 ± 0.0038. This implies that, for 0 < b ≤ 10, the constraint on the β is as 0.345 < β < 0.387. Planck2018 TT, TE, EE+lowE+lensing +BAO +BK14 constraint on the tensor-to-scalar ratio is r < 0.072, leading to 0.341 < β < 1 for 0 < b ≤ 10. From Planck2018 TT, TE, EE +lowE+lensing+BK14+BAO+LIGO and Virgo2016 data, the constraint on the tensor spectral index is −0.62 < n T < 0.53. In this regard, we find that for 0 < b ≤ 10, the range 0.044 < β < 1 leads to the observationally viable values of the tensor spectral index in the intermediate MDBI model. We have also studied r − n s and r − n T behaviors of the intermediate MDBI model in comparison to the observational data at 68% CL and 95% CL and found some constraints summarized in the tables.
Another important aspect of the inflation models is the non-gaussian feature of the primordial perturbations. In this paper, we have studied the non-gaussianity in the equilateral configuration. We have considered the k 1 = k 2 = k 3 limit, where the equilateral configuration has a peak. By studying this configuration of the non-gaussianity and considering the observationally viable ranges of the parameters b and β, we have predicted the amplitudes of the non-gaussianity in our intermediate MDBI model.
As a summary, we have shown that our proposed intermediate MDBI model with 0 < b ≤ 10 and 0.345 < β < 0.387, is instabilities-free and gives the observationally viable perturbation and non-gaussianity parameters.