Anisotropic constant-roll inflation for the Dirac-Born-Infeld model

In this paper, we study a non-canonical extension of a supergravity-motivated model acting as a vivid counterexample to the cosmic no-hair conjecture due to its unusual coupling between scalar and electromagnetic fields. In particular, a canonical scalar field is replaced by the string-inspired Dirac-Born-Infeld one in this extension. As a result, exact anisotropic inflationary solutions for this Dirac-Born-Infeld model are figured out under a constant-roll condition. Furthermore, numerical calculations are performed to verify that these anisotropic constant-roll solutions are indeed attractive during their inflationary phase.


I. INTRODUCTION
Cosmic inflation has played as one of the leading paradigms in modern cosmology [1]. This is due to the fact that many theoretical predictions of cosmic inflation have been well confirmed by the cosmic microwave background radiation (CMB) observations of the Wilkinson Microwave Anisotropy Probe (WMAP) [2] as well as the Planck [3] satellites.
It appears that all standard inflationary models have been investigated under the cosmological principle [4], which states that our universe is just simply homogeneous and isotropic on large scales as described by the Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime [5]. However, testing the validity of this principle seems to be not straightforward [6]. Recently, some anomalies of the CMB temperature such as the hemispherical asymmetry and the cold spot, which have been detected by the WMAP and confirmed by the Planck, have challenged the validity of the cosmological principle [7]. Consequently, this issue has received a lot of attention. In particular, a number of mechanisms have been proposed to explain the origin of these anomalies [7]. For example, some papers have appeared with an interesting idea that the mentioned anomalies could just be instrumental rather than cosmological [8]. In particular, the investigations in these papers have pointed out, using either the WMAP or Planck data, that the asymmetric beams could lead to the CMB statistical anisotropy. However, an independent analysis in Ref. [9] has been performed to show that the asymmetric beams seem to be unimportant. It is worth noting that a recent study has suggested that the Hubble tension might be a smoking gun for a breakdown of the FLRW cosmology [10].
All of these motivations lead us to think of a modification of the cosmological principle. One of reasonable possibilities is that the isotropy of spacetime of the early time universe should be slightly broken down [7,11]. Consequently, the spacetime of the early time universe would be described by the Bianchi metrics [11,12], which are homogeneous but anisotropic, rather than the FLRW one. If the early time universe was anisotropic, one could ask a related question that would the late time universe still be anisotropic or not. It appears that answering this question is not straightforward task at all. Theoretically, there has existed an important hint based on the so-called cosmic no-hair conjecture proposed by Hawking and his colleagues a few decades ago [13]. According to this conjecture, a late time state of the universe should obey the cosmological principle, i.e., should be homogeneous and isotropic, regardless of any initial states of the universe, which might have hairs, i.e., inhomogeneity and/or anisotropy. However, it is very complicated to draw a complete proof to this conjecture. Indeed, enormous efforts have been made to prove this conjecture for specific scenarios since the first partial proof by Wald for the Bianchi metrics in the presence of the cosmological constant Λ [14][15][16][17]. It is important to note that if the cosmic no-hair conjecture was valid, it would only be valid locally, according to the investigations in Refs. [18,19]. Along with the theoretical hint, it is worth noting that there has been a recent observational claim that the current universe might be anisotropic [20]. If this observation was the case, it would be an important issue directly related to validity of the cosmic no-hair conjecture.
Along with the proofs, many papers have examined the validity of the cosmic no-hair conjecture within specific inflationary models to see whether it is broken down or not [21][22][23][24][25][26]. It turns out that many papers among them have claimed to admit counterexamples to the conjecture. However, most of these claimed counterexamples have been pointed out to be invalid due to their instability [24]. Up to now, an anisotropic inflation found in a supergravitymotivated model by Kanno, Soda, and Watanabe (KSW) has been regarded as the first counterexample to the cosmic no-hair conjecture [25,26]. In particular, this anisotropic inflation has been shown to be stable and attractive due to the existence of supergravity-motivated coupling between scalar and vector fields, f 2 (φ)F µν F µν [26]. Consequently, this model has received a lot of attention. As a result, a number of its extensions have been proposed in order to see whether other counterexamples to the conjecture appear or not [27][28][29][30][31][32][33][34]. It motivated a revisiting of the Wald's no-hair theorem for inflationary settings [35]. Additionally, CMB imprints of the anisotropic inflation within the context of the KSW model as well as its non-canonical extensions have been investigated in Refs. [36][37][38][39][40], while primordial gravitational waves in anisotropic inflation have been studied in Ref. [41]. For other cosmological aspects of the KSW anisotropic inflation, see interesting reviews in Ref. [42].
Among non-trivial extensions of the KSW model, there are two interesting approaches, which we are currently interested in. The first approach is non-canonical extensions [31][32][33], in which a canonical scalar field is replaced by non-canonical ones such as the Dirac-Born-Infeld (DBI) field [31], whose origin can be realized in the D3-brane theory [43]. For interesting cosmological aspects of the DBI inflation, see Refs. [44][45][46][47][48]. As a result, these non-canonical models have been shown to admit stable and attractive anisotropic power-law inflationary solutions, which do violate the validity of the cosmic no-hair conjecture [31][32][33]. This result really confirms the leading role of the coupling f 2 (φ)F µν F µν in breaking down the validity of the cosmic no-hair conjecture. The last approach is a new type of anisotropic inflation, which has been found within the KSW model under a constant-roll condition [34].
Basically, the constant-roll condition is based on an assumption that one of slow-roll parameters, η = −φ/(Hφ), should be constant during an inflationary phase [49]. Interestingly, this constant-roll inflation is regarded as an interpolation between the well-known slow-roll inflation and a novel ultra-slow-roll inflation proposed in Ref. [50].
Hence, it would be interesting if we were able to figure out anisotropic inflation under the constant-roll condition within the KSW model. It is worth noting that the study in Ref. [34] has pointed out that it is possible to have such anisotropic inflation for the KSW model of canonical scalar field. Furthermore, it has been shown that the anisotropic constant-roll inflation is really attractive [34]. All of these motivations lead us to propose to investigate in this paper whether the DBI extension of the KSW model [31] admits anisotropic constant-roll inflation. As a result, we will point out the existence of anisotropic DBI constant-roll inflationary solutions. Furthermore, we will numerically confirm the attractive property of these solutions.
As a result, the present paper will be organized as follows: (i) A brief introduction has been written in Sec. I. (ii) Basic setup of the DBI model will be presented in Sec. II. (iii) Anisotropic constant-roll inflationary solutions will be solved in Sec. III. (iv) Attractor property of the obtained solutions will be shown in Sec. IV. (v) Finally, concluding remarks will be given in Sec. V.

II. DIRAC-BORN-INFELD MODEL
In this paper, we would like to study a non-canonical extension of the KSW model [31], in which a canonical scalar field is replaced by the non-canonical Dirac-Born-Infeld (DBI) field [43] as follows where the reduced Planck mass M p is set to be one for convenience. In addition, V (φ) is the potential of scalar field φ and F µν ≡ ∂ µ A ν − ∂ ν A µ is the field strength of the vector field A µ . It is noted that the notation γ defined as [43] γ ≡ plays as the Lorentz factor characterizing the motion of the D3-brane [43]. It is clear that γ ≥ 1 for non-negative f (φ). The last ingredient in the action (2.1) is a non-trivial coupling function h(φ), whose existence breaks down the conformal invariance of the Maxwell field. It should be noted that the existence of f (φ) in the definition of γ leads us to modify the original KSW coupling f 2 (φ)F µν F µν to h 2 (φ)F µν F µν in order to avoid any misunderstanding [31]. It is clear that if we take the canonical limit, i.e., f → 0, then the DBI model will reduce to the original KSW model. As a result, a stable and attractive Bianchi type I power-law inflation has been found in this DBI model [31]. This implies that the non-canonical property of the DBI field does not affect on the breaking of the cosmic no-hair conjecture within the KSW model. In other words, the coupling between the scalar and vector field, h 2 (φ)F µν F µν , does play the leading role in violating the cosmic no-hair conjecture. It is worth noting that the CMB imprints of non-canonical anisotropic model have been investigated recently in Ref. [40].
As a result, the corresponding field equations of the DBI can be derived as where V ≡ dV /dφ and so on. In this paper, similar to the previous work in Ref. [34], we will consider the following Bianchi type I metric, which is homogeneous but anisotropic, here a(t) acts as an isotropic scale factor, while b(t) plays as a deviation from isotropy. In addition, we will assume that the scalar field is homogeneous, i.e., φ = φ(t), while the vector field A µ takes the following configuration as A µ = (0, A x (t), 0, 0) in order to be compatible with the Bianchi type I metric. As a result, a non-trivial solution to the field equation of vector field (2.5) can be solved to bė where p A is an integration constant. Thanks to this solution, the corresponding field equation for the scalar field (2.4) now reduces toφ where the corresponding Lorentz factor now becomes (2.9)  . Once again, it appears that if we take the canonical limit γ → 1, or equivalently f → 0, then all these field equations will reduce to that derived in the KSW model for canonical scalar field [34]. Additionally, if a(t) and b(t) both take exponential forms as a(t) = exp[α(t)] and b(t) = exp[σ(t)] then all these field equations will become that derived in Ref. [31] for the DBI field.

A. Constant-roll conditions
We recall that in the case of canonical inflation, γ = 1, Eq. (2.8) reduces tö It is noted that the slow-roll inflation is characterized by some slow-roll parameters, which are normally constrained by the observational data of WMAP or Planck [2,3]. One of them is defined as [49] Usually, the value of η should be much smaller than one within the context of slow-roll inflation. In the constant-roll inflation [49], η has been assumed to be constant during the inflationary phase, i.e., withβ is a constant. Note thatβ is related to other definitions of constant-roll constant in Refs. [34,49,54] such aŝ β = −β = 3 + α with α and β are also constants. So α → −3 will lead toβ → 0. As a result, along with the power-law inflation [51] other types of inflation corresponding to different potentials can be obtained under this constant-roll assumption [49]. More interestingly, the similar result also happens in the KSW model [34]. It has been shown in Ref. [65] that the slow-roll parameter η could be modified for the DBI scalar field. To see how to modify η, let us rewrite Eq. (2.8) as By comparing this equation with that of canonical scalar field (3.1), we can arrive at the modified definition of η for the DBI model such as which is nothing but that introduced in a recent paper [65]. However, it turns out that η DBI = η if γ = constant. It should be noted that the condition γ = constant is necessary for obtaining the corresponding power-law inflation to the DBI model as shown in Refs. [31,47]. Hence, one can expect the same thing for the DBI constant-roll inflation, which involves the power-law inflation as its subclass. Hence, we will assume from now on the constant Lorentz factor condition, with γ 0 ≥ 1 is a constant. Consequently, the constant-roll condition for the DBI model is now given by Additionally, we will impose the constant anisotropy condition from Ref. [34], with n is also a constant, for seeking anisotropic constant-roll solutions. It appears that |n| 1 is a sufficient constraint for viable anisotropic constant-roll inflation.
In the following steps, we will employ the method presented in Ref. [34] with some modifications necessary for the DBI model. More specifically, the strategy of seeking solutions is as follows [34]: (i) First, we use three conditions shown in Eqs. (3.6), (3.7), and (3.8) to derive the corresponding forms of V (φ), h(φ), and f (φ). (ii) Then, we put back the derived forms of V (φ), h(φ), and f (φ) to solve numerically the field equations (2.8), (2.10), (2.11), and (2.12) for different initial conditions along with some values of γ 0 . The attractor properties such as η DBI (t) →β, H b (t)/H a (t) → n, and γ(t) → γ 0 after some e-folds during the inflationary phase, will therefore be confirmed by numerical calculations accordingly.

B. Solutions
As a result, we get, from Eqs. (2.11), (2.12), and (3.8), an equation of H a such as which can be rewritten, by regarding H a as a function of φ, as follows It turns out that Eq. (3.10) is a quadratic equation ofφ and can therefore be easily solved to give the corresponding non-trivial solutions, Furthermore, by differentiating these solutions with respect to t, we get along with the following important relation between two parameters n andβ [34], n =β 6 . (3.14) Here, C 1 and C 2 are integration constants. It should be noted that this solution is different from that derived for an isotropic inflation with n = 0 in Ref. [49]. Similarly, we are able to figure out, thanks to the relation (3.14), a generalized solution for H a in the DBI model such as It turns out that different values of C 1 and C 2 will lead to different types of inflation. It is interesting to note that if we set either C 1 = 0 or C 2 = 0, the corresponding anisotropic power-law inflation [31] can be derived accordingly. Note again that the condition γ = γ 0 has been shown to be attractive in the DBI model for anisotropic power-law inflation [31]. Furthermore, following Refs. [34,49] we will focus on three special solutions derived from the general one (3.15), where M is an integration constant determining the amplitude of the power spectrum of the curvature perturbation [49]. As a result, the first solution (3.16) is obtained by setting C 1 = M and C 2 = 0, while the second one (3.17) is derived by setting C 1 = C 2 = M/2. Additionally, the last solution (3.18) corresponds to C 1 = M/2 and C 2 = −M/2.
It is noted that an unlisted solution corresponding to C 1 = 0 and C 2 = M , i.e., H a = M exp − 6γ0β 6+β φ , can be obtained from the solution in Eq. (3.16) by redefining φ → −φ [49]. It is clear thatβ should be positive definite for the first solution shown in Eq. (3.16), whileβ can be negative definite for the last two solutions shown in Eqs. (3.17) and (3.18) [49]. However, the relation (3.14) implies that |β| should be much smaller than one such that |H b | H a in order to be consistent with the observations. It appears that these solutions will recover that derived in Ref. [34] for the canonical scalar field once taking a limit γ 0 → 1. However, these solutions as well as the solutions found in Ref. [34] cannot reduce to that derived in Ref. [49] for an isotropic inflation under a limit n → 0. This is because these solutions have been derived using a special relation (3.14), which was firstly adopted in Ref. [34] for the anisotropic inflation.  and, up to a constant, respectively. From the definition of γ shown in Eq. (2.9), we obtain the corresponding definition of f (φ) as with the help of the solution (3.19). In conclusion, for the solution (3.16) under the constant-roll condition, the corresponding functions f (φ) and h(φ) as well as the corresponding potential V (φ) all appear to be exponential functions of φ. As a result, Eq. (3.19) can be solved to give the following solution, which leads the scale factors a(t) and b(t) to the following power-law expansion, up to a constant, 26) respectively. Hence, it turns out that 0 <β 1 is a sufficient condition in order to have anisotropic inflation. This power-law solution is nothing but that derived in Ref. [31] for the DBI model. As a result, substituting the second solution (3.17) into Eq. (3.11) indicates thatβ must be negative definite in order to makeφ real definite. However, it is confirmed in Ref. [34] that this solution is not an attractor for the canonical scalar field with γ 0 = 1. Additionally, this is also the case for the DBI model. Therefore, we will ignore this solution from now on. It should be noted that the "cosh" type solution of Hubble constant H has been shown to be physical for the isotropic inflation [49]. This is indeed a different point between the isotropic constant-roll inflation found in Ref. [49] and the anisotropic constant-roll inflation found in Ref. [34] as well as in this paper.
Consequently, the corresponding definitions of V (φ), h(φ), and f (φ) turn out to be respectively. Here, we have only consideredβ > 0. Forβ < 0, it appears that the solution (3.18) is similar, up to a field redefinition, to the solution (3.17) [49], and therefore is neglected. It turns out that these solutions will be reduced to that derived for the canonical scalar field in Ref. [34] by taking the limit γ 0 → 1. As a result, once Eq. (3.27) is solved to give the solution of φ, we can easily determine the related scale factors a(t) and b(t) as well as the functions f (φ), h(φ), and V (φ). As a result, a non-trivial solution of Eq. (3.27) for the upper " + " sign is given by which leads to the corresponding scale factors given by It appears that a(t), b(t), and φ(t) are all real definite for M 2β2 t 2 > 1. However, it turns out that only the region Mβt > 1 is suitable for inflationary solutions since the other region Mβt < −1 leads toȧ(t) < 0. Moreover, it is clear that inflation will happen if 0 <β 1. Hence, M should be positive definite in this case. Additionally, it is noted that this solution in the canonical limit, γ 0 = 1, has been claimed to suffer a strong coupling problem around φ = 0, i.e., h(φ 0) < 1 [34]. Consequently, this problem leads to a requirement that the inflation needs to finish before φ becomes zero. It is straightforward to check that this result is also the case for γ 0 > 1 (see Fig. 1 for a demonstration). On the other hand, a non-trivial solution of Eq. (3.27) for the lower " − " sign can be defined to be which implies the corresponding scale factors given by It appears that a(t), b(t), and φ(t) are all real definite for −1 < Mβt < 1. Sinceȧ(t) < 0 for 0 < Mβt < 1, we will only be interested in the region −1 < Mβt < 0. Similar to the above solutions, it is clear that inflation will happen if 0 <β 1. Hence, M should be negative definite in this case. Interestingly, the corresponding h(φ) does not meet the strong coupling problem in this case. However, the corresponding potential V (φ) will be negative definite around φ = 0 provided that 0 <β 1 (see Fig. 2 for a demonstration). Consequently, the corresponding inflation also needs to finish before φ approaches zero.
In conclusion, we have completed the first step of the proof that is working out the corresponding forms of V (φ), h(φ), and f (φ) compatible with the constant-roll condition, η DBI (t) =β, and the associated ones, H b (t)/H a (t) = n =β/6 and γ(t) = γ 0 . Up to now, we have figured out three different inflationary solutions, that are: • solution I: is described by a set of Eqs.  Here, h0 is an integration constant. This plot shows that h(φ) is always larger than one. In addition, it also shows that V (φ) is negative definite when φ approaches zero (the negativity region of V (φ) is displayed as the dotted line).

IV. ATTRACTOR PROPERTY
As mentioned above, we will prove in this section that the conditions (3.6), (3.7), and (3.8) are really attractive by solving numerically the field equations (2.8), (2.10), (2.11), and (2.12) for three inflationary solutions I, II, and III with three different values γ 0 = 1, 1.5, and 2. Once numerical solutions are solved, we will plot the corresponding time evolution of η DBI (t), H b (t)/H a (t), and γ(t) to see whether they become constant during the inflationary phase.
It is noted that 0 <β 1 for all three inflationary solutions. Additionally, it is also noted that M > 0 for both solutions I and II, while M < 0 for the solution III. Therefore, we will chooseβ = 0.1, |M | = 10 −4 , along with initial conditions for scalar field as φ(0) = 15 andφ(0) = 0, for all three solutions I, II, and III. By choosing that, it is easy to compare these solutions to each other. It is noted that these values are similar to that used in the canonical case [34] since we would like to see the effect of γ 0 .
The time evolution of H b (t)/H a (t) is depicted in Fig. 3. It is clear that this ratio converges to the assumed constant, n =β/6 0.0167, after 10−20 e-folds. On the other hand, the value of γ 0 does affect on the evolution of H b (t)/H a (t). In particular, the increase of the value of γ 0 tends to delay the convergence of H b (t)/H a (t). The time evolution of γ(t) = 1/ 1 − f (φ)φ 2 is displayed in Fig. 4. It is noted that γ(t) should always be larger than or equal to one. For γ 0 = 1, γ(t) should remain as one constantly as time evolves. As a result, these properties can be easily confirmed by the plots in Fig. 4. Additionally, these plots also clearly indicate that γ(t) → γ 0 after 10 − 20 e-folds. More interestingly, it is shown that the value of γ 0 also affects on the evolution of γ(t) in the same way it does for H b (t)/H a (t). These results confirm that the condition γ(t) = γ 0 as chosen in Eq. (3.6) is really attractive.
Next, we plot in Fig. 5 the phase space of H b (t)/H a (t) and η DBI (t). We see that for three different values of γ 0 , the trajectories all converge to the same fixed point β , n (0.1, 0.0167). Finally, we would like to discuss the difference between three solutions I, II, and III. According to the plots, it turns out that the behavior of these solutions are not significantly different from each other. The only gap between the first two solutions I and II and the last solution III is the e-fold number as displayed in Fig. 4. In particular, the inflationary phases corresponding to the first two solutions I and II still survive when the e-fold number reaches to 60. However, this is not the case for the last solution III. Indeed, the corresponding inflationary phase of this solution stops before the e-fold number reaches to 60 due to the negativity of the corresponding potential V (φ) happening when φ approaches zero as mentioned above. One additional interesting point should be noted that the value of γ 0 does affect on the e-fold number of the inflationary solution III. As a result, it turns out that the larger γ 0 is, the larger e-fold number the inflationary solution III has. This result can be explained according to the Fig. 6. In particular, this figure clearly displays that the potential V (φ) tends to become negative sooner when γ 0 decreases to one.

V. CONCLUSIONS
In this paper, we have extended our DBI model proposed in Ref. [31], which is a non-canonical extension of the KSW model [25], to a novel type of inflation called the constant-roll inflation [49]. Basically, our paper is just a non-trivial extension of the paper [34], in which anisotropic constant-roll inflationary solutions have been found to the KSW model. As a result, we have figured out three exact anisotropic constant-roll inflationary solutions named as I, II, and III, to the DBI model. As a result, these solutions have been achieved under the constant-roll condition [49,65], η DBI (t) =β, along with the compatible ones, H b (t)/H a (t) = n =β/6 for the Bianchi type I spacetime [34] and γ(t) = γ 0 for the DBI field. While the solution I is nothing but the well-known power-law type [31], the other ones II and III turn out to be novel ones. Numerical calculations have been performed to confirm that these conditions are really attractive. These results together with that investigated in Ref. [34] indicate that the cosmic no-hair conjecture is extensively violated in the KSW model as well as in its non-canonical extension, the DBI model, even when the constant-roll condition is applied. It would be interesting if ones could find anisotropic constant-roll inflationary solutions to other non-canonical models such as the Galileon inflation and the k-inflation [32,33,63] as well as to multi-field models [28,68]. Additionally, it is important to investigate the CMB imprints of the anisotropic constant-roll inflation. For example, one might ask whether this model satisfies the observational constraints on the spectral index and tensor-to-scalar ratio. However, it turns out that these observational issues are not straightforward to answer within not only the framework of constant-roll inflation [54,55,65] but also the framework of anisotropic inflation [38][39][40]. We will therefore leave them to our further studies, whose results will be presented in a sequel to this paper. We hope that our present paper would be useful to studies of anisotropic inflation as well as constant-roll inflation.