Anisotropic power-law inflation for models of non-canonical scalar fields non-minimally coupled to a two-form field

In this paper, we investigate the validity of the so-called cosmic no-hair conjecture in the framework of anisotropic inflation models of non-canonical scalar fields non-minimally coupled to a two-form field. In particular, we focus on two typical {\it k}-inflation and Dirac-Born-Infeld inflation models, in which we find a set of exact anisotropic power-law inflationary solutions. Interestingly, these solutions are shown to be stable and attractive during an inflationary phase using the dynamical system analysis. The obtained results indicate that the non-minimal coupling between the scalar and two-form fields acts as a non-trivial source of generating stable spatial anisotropies during the inflationary phase and therefore violates the prediction of the cosmic no-hair conjecture, even when the scalar field is of non-canonical forms. In connection with the Planck 2018 data, tensor-to-scalar ratios of these anisotropic solutions are investigated. As a result, it appears that the tensor-to-scalar ratio of the anisotropic power-law inflationary solution of {\it k}-inflation model turns out to be more highly consistent with the Planck 2018 data than that of Dirac-Born-Infeld model.


I. INTRODUCTION
Cosmic inflation was firstly introduced as a solution to the longstanding puzzles in the standard Big Bang cosmology, such as the horizon, flatness, and primordial monopole problems [1][2][3].More interestingly, a rapid expansion during the inflationary phase will stretch the primordial density perturbations, which are created from quantum fluctuations, from microscopic scales to galactic scales, and therefore the large-scale structure of the present universe can be produced from these primordial density perturbations [4].For a recent interesting review on cosmic inflation, see Ref. [5].
Notably, observational data of the cosmic microwave background radiation (CMB) from the Wilkinson Microwave Anisotropy Probe (WMAP) [6] and the Planck satellite [7,8] have been consistent very well with theoretical predictions of the standard inflationary models, whose underlying assumption is the cosmological principle, which states that our universe is spatially homogeneous and isotropic on large scales.In cosmology, there exists a unique spacetime called the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, which has both spatial homogeneity and isotropy.This is a reason why the FLRW metric has been widely used as a background spacetime in standard inflationary models [9].
However, some unavoidable anomalies in the CMB map, such as the cold spot and hemispheric asymmetry, have presented challenges to the standard inflationary models based on the cosmological principle [10].In other words, if the cosmological principle is valid during the inflationary phase, the probability of the existence of the CMB anomalies is very small.It is worth noting that some other interesting observational evidences against the validity of cosmological principle have been summarized in a recent interesting review [11].Therefore, violating the cosmological principle during the inflationary phase might be a reasonable resolution.And one of the simplest ways to archive this violation is replacing the standard FLRW spacetime with the Bianchi spacetimes, which are spatially homogeneous but anisotropic [12].Consequently, we will end up with an anisotropic inflation [10].Very interestingly, many theoretical predictions for anisotropic inflation had been worked out even when the observed anomalies were not detected [13,14].It turns out that anisotropic inflation has been considered seriously by many people, e.g., see Ref. [15] and references therein.
Besides the cosmological principle, there has existed the so-called cosmic no-hair conjecture proposed by Hawking et al. a long time ago, which is also about the homogeneity and isotropy of universe's spacetime [16].Basically, this conjecture implies that the late-time state of any accelerating universe is independent of its initial state.Eventually, the universe becomes isotropic and homogeneous as an attractor of cosmic evolution.Physicists and cosmologists have long been challenged by the cosmic no-hair conjecture.However, they have made notable progress in providing partial proofs for this conjecture.Remarkably, the first proof, specifically for the Bianchi spacetimes, which are homogeneous but anisotropic metrics, in the presence of a positive cosmological constant Λ, has been proposed by Wald using an approach based on energy conditions [17].It should be noted that the Wald's proof deals only with background anisotropies.Remarkably, Starobinsky pointed out independently that the cosmological constant is the best "isotropizer", i.e., it is capable of eliminating or extending over very large scales of all types of inhomogeneities [18].Furthermore, he concluded that the cosmic no-hair conjecture should be valid locally, i.e., inside the future de Sitter event horizon.More interestingly, an inhomogeneous time-independent tensor (we can regard it as hairs) existing outside of the future de Sitter event horizon may have an arbitrarily large amplitude.In order to archive this conclusion, Starobinsky considered an inhomogeneous metric, whose scale factors depend on not only a cosmic time but also three spatial coordinates.This metric is nothing but an extension of de Sitter one with small inhomogeneous perturbations.In a follow-up study, Starobinsky and his colleagues showed that inhomogeneous time-independent hairs also exist in a power-law inflation model of scalar field [19].It is worth noting that similar conclusions have been archived later by other people in different scenarios, in which background spacetimes are of exact anisotropic and/or inhomogeneous forms like the so-called Tolman-Bondi spacetime [20][21][22].According to these papers, one could state that the cosmic no-hair conjecture is invalid globally in the view of inhomogeneous hairs (a.k.a.constant perturbation modes in the super-Hubble regime) and may therefore be valid locally.Interestingly, these hairs could even be observable after the end of inflation and we could expect to see them directly in fluctuations of CMB temperature and polarization.
Recently, various cosmological models have been employed to test the validity of the cosmic no-hair conjecture, such that the higher curvature models [23,24], the Lorentz Chern-Simons theory [25], and the Horndeski models [26].It turns out that some claimed counterexamples from these models no longer hold due to their instability during the inflationary phase as indicated in Ref. [27].From 2009 to 2010, a vivid counterexample to the cosmic no-hair conjecture was successfully constructed by Kanno, Soda, and Watanabe (KSW) within a supergravity-motivated model involving a non-minimal coupling between the U (1) gauge and scalar fields such as f 2 (ϕ)F µν F µν [28,29].As a result, the KSW model admits homogeneous and anisotropic Bianchi type I solutions, whose spatial anisotropies turn out to be stable during the inflationary phase, in contrast to the prediction of the cosmic no-hair conjecture.Subsequently, several non-trivial extensions of the KSW model have been proposed by considering non-canonical scalar fields instead of canonical one, such as the Dirac-Born-Infeld (DBI) model [30][31][32][33], the generalized ghost condensate model [31], supersymmetric DBI model [34], and the k-inflation model [35].As a result, the cosmic no-hair conjecture is always violated in these non-canonical extensions.For CMB imprints of non-canonical anisotropic inflation, see Refs.[34][35][36], while that of canonical one can be found in Refs.[37][38][39][40][41][42].All these results indicate that the nonminimal coupling, f 2 (ϕ)F µν F µν , has played a leading role in order to generate stable spatial anisotropies during the inflationary phase.Interesting reviews on cosmological implications of the anisotropic inflation based on the vector field can be found in Refs.[43,44].
One might ask if spatial anisotropies can also be caused by other mechanisms.Interestingly, people have figured out that a non-minimal coupling between scalar and two-form fields such as f 2 (ϕ)H µνρ H µνρ could play a similar role as the coupling f 2 (ϕ)F µν F µν [45][46][47][48][49].It is well known that a two-form field B µν can be found in string theory [50].In Ref. [45], the authors have pointed out that the anisotropy induced by the two-form field corresponds to the prolate type, i.e., the expansion of the Universe slows down in the (y, z) plane, in contrast to the oblate type stemming from the vector field (a.k.a.one-form field).In Ref. [46], the authors have derived the corresponding observational constraints of two anisotropic inflation models, one for the vector field and the other for the two-form field.Very interestingly, some significant gaps between observational predictions of these two types of anisotropic inflation have been identified accordingly.The authors have come to the main conclusion that the precise measurements of g * as well as the TB correlation will clarify which anisotropic model is favored over the other [46].In Ref. [47], the authors have investigated anisotropic hairs in the presence of both one-form and two-form fields.As a result, they have obtained an important conclusion that there always exists one stable anisotropic fixed point in this model.In Ref. [48], the authors have indicated that the five-dimensional (5D) two-form field can be shown to be equivalent to a 5D gauge field via a Routh transformation.Interestingly, the cosmic no-hair conjecture has been shown to be broken down in this 5D model.In Ref. [49], the authors have considered a more general scenario, which involves not only one-and two-form fields but also a three-form one.It should be noted that all these studies of anisotropic inflation in the presence of the two-form field have only been done for the canonical field.Other physical and cosmological implications of two-form field can be found in Refs.[51][52][53][54][55][56][57][58][59][60][61].
Motivated by studies of anisotropic inflation models of non-canonical scalar fields [30,31,34,35] and anisotropic inflation models of two-form fields [45][46][47][48][49], we would like to investigate in the present paper the validity of the cosmic no-hair conjecture in a mixed scenario, in which a non-canonical scalar field is allowed to non-minimally coupled to the two-form field.Specifically, we focus on two well-known types of non-canonical scalar fields, one is from the k-inflation model [62,63] and the other is from the string-inspired DBI inflation model [64][65][66].As a result, we are able to find exact anisotropic power-law solutions for both non-canonical scalar fields.Furthermore, we confirm, using the dynamical system method, that these solutions are indeed stable and attractive during the inflationary phase and therefore act as additional counterexamples to the cosmic no-hair conjecture.In connection with the Planck 2018 data [7,8] as well as future detections like the CMB-S4 project [67], we will derive a tensor-to-scalar ratio for a general case of non-canonical scalar field non-minimally coupled to the two-form field.Then we will focus on the two mentioned models, i.e., k-inflation and DBI ones, to see whether the corresponding ratios are consistent with the Planck 2018 data or not.Very interestingly, the obtained tensor-to-scalar ratio of the k-inflation two-form field model turns out to be highly consistent with the Planck 2018 data.
As a result, this paper will be organized as follows: (i) An introduction of our study has been written in Sect.I. (ii) A general action of studied models will be presented in Sect.II.(iii) In Sect.III, we derive the corresponding set of anisotropic power-law solutions for k-inflation model and investigate its stability.(iv) In Sect.IV, we extend our analysis to the DBI model.(v) In Sect.V, we compare the anisotropic parameter |Σ/H| derived in this paper with that obtained in the previous papers for heuristic reasons.(vi) In Sect.VI, tensor-to-scalar ratios of the obtained anisotropic power-law inflationary solutions will be investigated.(vii) Finally, concluding remarks will be written in Sect.VII.Additional calculations will be presented in the Appendix.

II. GENERAL ACTION
Let us begin by introducing a general action of non-canonical extension of the KSW model [31,36], where P (ϕ, X) is an arbitrary function of scalar field ϕ and its kinetic term defined as X ≡ −(1/2)∂ µ ϕ∂ µ ϕ [62].In addition, f (ϕ) is the gauge kinetic function depending only on ϕ, while the rank-2 tensor is the field strength of the vector field A µ .It is noted that we have set the reduced Planck mass as one, i.e., M p = 1, just for convenience.Specific forms of P (ϕ, X) have been considered in previous papers [30,31,34,35].
In this paper, we would like to investigate a modification of the above action, in which the one-form field A µ is replaced by a two-form field B µν .As a result, the corresponding action is given by where Without the loss of generality, one can assume that the (y, z) plane aligns with the direction of the two-form field.Consequently, we can express B µν in the following form [45-47] where v B (t) is a function of cosmic time t.In this paper, we will study action (2) for two typical non-canonical scalar fields, one is from the k-inflation model [35,62] and the other is from the string-inspired DBI inflation model [30,[64][65][66].

III. K-INFLATION CASE
In this section, we focus on the k-inflation model [35,62], in which the function P (ϕ, X) takes the following form P (ϕ, X) = K(ϕ)X + L(ϕ)X 2 .Thus, the action (2) becomes here K(ϕ) and L(ϕ) are arbitrary functions of ϕ.Following Refs.[45][46][47], we will adopt the Bianchi type I metric, as the background spacetime for the cosmic evolution.Here, N (t) represents the lapse function that allows us to derive the Friedmann constraint equation, while α represents the average expansion measured in terms of the number of e-foldings and σ corresponds to the spatial anisotropy.By substituting the background metric into the action (5), we obtain the following expression, where an overdot stands for a derivative with respect to the cosmic time t.The equation of motion for the two-form field can be solved to give a non-trivial solution of v B (t) as where p B is a constant of integration [45][46][47].By varying the above action (7) with respect to N , α, σ, as well as ϕ, we are able to derive the following background equations of motion, after setting N = 1, respectively.Here, the subscript in K ϕ , L ϕ and f ϕ indicates a derivative with respect to the field ϕ, i.e., K ϕ ≡ ∂K/∂ϕ.To figure out analytical solutions for these derived field equations, we will consider the following ansatz as used in many previous papers [29,30,35,[45][46][47], along with the exponential functions of the scalar field, where k 0 , l 0 , f 0 , ζ, η, ξ, ϕ 0 , λ, κ, and ρ are all constant.For convenience, we will introduce the following new parameters to aid in our analysis, Thus, the field Eqs. ( 9), ( 10), (11), and ( 12) can be simplified to the following algebraic equations, In addition to these equations, we have the following constraint equations that ensure all terms in the field equations are proportional to t −2 , Apparently, the condition for κ = 0 gives us a constant value of K(ϕ) = k 0 for the ansatz (13).
A. Anisotropic power-law inflation We will focus on seeking anisotropic solutions.From Eq. ( 25) we have On the other hand, we have from Eqs. ( 21) and ( 22) that thanks to the constraint equation ( 23) as well as the relation shown in Eq. ( 26).Substituting η, u, and w defined above into either Eq. ( 19) or Eq. ( 20), we can find non-trivial solutions of ζ such as Now, we rewrite Eq. ( 26) as follows This equation implies an important point that the ratio −2ρ/λ will mainly determine the value of ζ since η should be much smaller than ζ.Therefore, the constraint for the existence of inflation, ζ ≫ 1, will imply that ζ ≃ −2ρ/λ ≫ 1, or equivalently |ρ| ≫ |λ|.Consequently, we observe that only the solution, is suitable for describing the inflationary phase.Note that the positivity of ρ will imply the negativity of λ, or vice versa.In this paper, we will prefer the choice, in which ρ is assumed to be positive and λ will be negative.Then, the corresponding expression for η turns out to be while the corresponding anisotropy parameter reads The real values of ζ and η require that k 0 must satisfy the following constraint, However, the smallness of η implies, according to Eq. ( 32), the following constraint, which leads to an approximated value of k 0 such as due to the constraint, ζ ≃ −2ρ/λ.It is clear that this value is consistent with the inequality (34).Consequently, the value for ζ and η can be approximate as follows It is clear that the negativity of η is consistent with the positivity of w, according to Eq. ( 21).However, the negativity of η in this model of two-form field is indeed in contrast to the positivity of η required in models of vector field [29,30,35].This result is consistent with the previous investigation for canonical scalar field [47].
As a result, the following approximated value of the anisotropy parameter turns out to be Of course the absolute value of this ratio is much smaller than one, i.e., |Σ/H| ≪ 1 as expected.In fact, to be consistent with the cosmological observation the absolute value of the anisotropy parameter, i.e., |Σ/H|, must be much smaller than one [28,29].For heuristic reasons, we will compare |Σ/H| derived in our current model with that obtained in the KSW model [29] and the k-inflation model [35] of vector field, as well as with that derived in Ref. [47] for the canonical scalar field coupled to the two-form field, using specific values of field parameters such as |λ| = 0.1 and ρ = 50.As a result, we obtain the corresponding values as |Σ/H| KSW ≃ 0.0004, |Σ/H| k-one-form ≃ 0.0005, |Σ/H| canonical-two-form ≃ 0.0008, and |Σ/H| k-two-form ≃ 0.0003 < |Σ/H| canonical-two-form .It turns out that the anisotropy induced by the two-form field non-minimally coupled to the k-inflation field is the smallest one among these four values.

B. Stability analysis
To investigate the stability of the anisotropic power-law solution defined above we will introduce the corresponding dynamical variables [47], along with two auxiliary variables [35], As a result, the corresponding dynamical system of autonomous equations are defined to be Here, α = αdt is understood as a new time coordinate [29,30,35,47].In the above equations, there exist α, σ, and φ, which can be determined from the field equations ( 10), (11), and (12) with the help of the Friedmann constraint equation (9), which now becomes as As a result, explicit expressions of autonomous equations can be defined to be Now, we are going to find out anisotropic fixed points with x ̸ = 0 to this dynamical system by solving the following equations, First, the equation dω κ /dα = 0 gives κ = 0 and then ω κ = 1, consistent with the power-law solution derived in the previous subsection.As a result, the equation dω λ /dα = 0 gives while the equation dz/dα = 0 implies Then, we can obtain the expression for y from these two equations, On the other hand, by using two equations, dx/dα = 0 and dz/dα = 0, we can obtain the expression for z 2 as Thanks to these results, we are able to obtain non-trivial solutions for x from either equation dy/dα = 0 or dz/dα = 0, It is clear that the solution x + is absolutely equivalent to the anisotropic power-law one found in the previous subsection due to the result x + = η/ζ − .This indicates that the stability of the anisotropic fixed point is also that of the anisotropic power-law solution.
To examine the stability of the anisotropic fixed point, we perturb the autonomous equations around this fixed point as follows Then, we take exponential perturbations, to obtain the corresponding matrix equation, It is well known that this matrix equation admits non-trivial solutions if and only if which can be expanded to be a polynomial equation for ω as follows where Apparently, the constraint that λ < 0 and ρ > 0 will indicate the fact that all the coefficients a i (i = 0 − 3) are positive.Hence, the polynomial equation ( 73) only admits negative roots ω < 0. Therefore, the anisotropic fixed point turns out to be stable against perturbations.In addition, the numerical result shown in Fig. 1 confirms that the anisotropic fixed point is clearly attractive since trajectories with different initial conditions all converge to the anisotropic fixed point (displayed as the black point) rather than an isotropic fixed point with x = z = 0.For now, we can conclude that the anisotropic power-law solution is stable and attractive during the inflationary phase.And this model provides us one more counterexample to the cosmic no-hair theorem.

IV. DIRAC-BORN-INFELD CASE
In this section, we extend our analysis to the DBI model, whose origin comes from string theory [64][65][66].Anisotropic power-law inflation for the DBI model in the presence of one-form field has been firstly studied in Ref. [30] then revisited in Refs.[31][32][33].An action of the Dirac-Born-Infeld model in the presence of the two-form field is given by where we have renamed the gauge kinetic function as f (ϕ) → h(ϕ) to avoid any misunderstanding.It turns out that the action (79) will become that of canonical scalar field considered in Refs.[45][46][47] if the limit f (ϕ) → 0 is taken.Similar to the previous case, we will use the Bianchi type I metric shown in Eq. ( 6) to define the corresponding form of the action (79), As a result, the corresponding field equations of this model are given by here we have used the following definition, as well as the solution of v B (t), It is noted that γ ≡ 1/ 1 + f (ϕ)∂ µ ϕ∂ µ ϕ is called the Lorentz factor characterizing the motion of the D3-brane [64,66].The positivity of f (ϕ) implies that γ > 1.The model of canonical scalar field studied in Ref. [47] will correspond to γ = 1.

A. Anisotropic power-law solution
With the above setting, to find exact anisotropic power-law solutions to the DBI model we continue to use the ansatz shown in Eq. ( 13) along with exponential functions such as [30] V where V 0 , f 0 , h 0 , λ, τ , and ρ are all constant.We also defined new variables [30], It turns out that the corresponding γ is given by Obviously, if τ ξ = 2 then γ will become a constant γ 0 defined as Now, we are going to find out ζ and η from the field equations because their values will determine whether the inflationary phase exists or not.Similar to the previous analysis, we will define a set of the algebraic equations from the field equations ( 81), (82), (83), and (84) to be with the help of the corresponding constraints, As a result, the relation between λ and τ can be expressed as follows Hence, the positivity of λ will imply the negativity of τ and vice versa.After some algebra, a set of non-trivial solutions is found to be As a result, the scale factors for this model are given by Apparently, the positive scale factors require that both λ and ρ are positive definite.Therefore, τ should be negative according to Eq. ( 102).As a result, the corresponding anisotropy parameter is given by It is straightforward to verify that these solutions will recover that defined in Ref. [47] in the canonical limit γ 0 → 1.
For the anisotropic power-law inflation with ρ ≫ λ, the ratio Σ/H is indeed much smaller than one as expected.It appears that Σ/H DBI-two-form ≃ 0.0005 < Σ/H canonical-two-form ≃ 0.0008 provided λ = 0.1, ρ = 50, and γ 0 = 1.5.This result indicates that the larger γ 0 is, the smaller Σ/H is.Therefore, the anisotropy parameter Σ/H seems to be reduced in the context of non-canonical scalar field models, according to this result as well as that in the previous case.

B. Stability analysis
The stability of the obtained anisotropic power-law solution of the DBI model will be investigated in this subsection.To do this task, we firstly introduce the corresponding dimensionless dynamical variables as Consequently, the field equations ( 81), ( 82), (83), and (84) will be converted into the corresponding dynamical system of autonomous equations given by Here, we have defined an auxiliary variable as γ = 1/γ for the completeness of dynamical system [30,66].Now, we are going to seek anisotropic fixed points with x ̸ = 0 to this dynamical system by solving a set of equations, Combining both equations, dx/dα = 0 and dz/dα = 0, will yield an equation of z 2 , On the other hand, the equation dγ/dα = 0 gives Finally, the equation dz/dα = 0 implies Thanks to these useful results, an important expression for y can be revealed to be Now, we are about to obtain an equation of x by substituting Eq. ( 116) and Eq. ( 119) into the equation for dz/dα = 0, Solving this equation gives a non-trivial solution for x, Very interestingly, the equation dy/dα = 0 can be rewritten as which clearly implies that τ = −λ, or equivalently γ = γ0 with γ0 = 1/γ 0 .This result is indeed consistent with the power-law solution found above.It is also consistent with Ref. [30], in which the DBI field is non-minimally coupled to the vector field.Consequently, a complete set of anisotropic fixed point solution could be represented in terms of λ, ρ, and γ0 as follows It is clear that the positivity of z 2 will be ensured if the following inequality, is satisfied.Furthermore, this anisotropic fixed point is exactly equivalent to the anisotropic power-law inflation solution found in the previous subsection.Therefore, we will examine the stability of the anisotropic fixed point solution, which will tell us the stability of the corresponding anisotropic power-law solution.
To do this, we will first take approximated values for x, y, and z 2 , during the inflationary phase, Taking exponential perturbations, we are able to define the following perturbed equations, which can be written in a matrix equation given by Consequently, the corresponding equation for ω is defined as where It is noted that the positivity of the last coefficient b 0 is ensured by the positivity of z 2 .It is crystal that Eq. ( 134) with all positive coefficients will always admit negative roots of ω.Thereby, the anisotropic fixed point of the DBI case is indeed stable during the inflationary phase.Consistently, the numerical result displayed in Fig. 2 points out that this fixed point is attractive.It is safe to conclude that the anisotropic power-law of the DBI case is stable and disfavors the cosmic no-hair conjecture.
FIG. 2: Attractive behavior of the anisotropic fixed point of the DBI case (displayed as a black point) for the parameters chosen as λ = 0.1, ρ = 50, and γ 0 = 1.5.

V. COMPARISONS OF |Σ/H| FOR DIFFERENT MODELS A. Two-form field models
For heuristic reasons, we would like to compare the anisotropy parameter, |Σ/H|, for three different anisotropic power-law inflation models of two-form field.As a result, the expression of |Σ/H| for the canonical scalar field can be found in Eq. ( 48) of Ref. [47].According to Fig. 3, it appears that |Σ/H| varies at different rates for the canonical, k-inflation, and DBI scalar fields.Specifically, Fig. 3a contains three curves of |Σ/H| for the fixed values ρ = 50 and γ 0 = 1.5, along with |λ| running from 0.03 to 0.5.On the other hand, Fig. 3b is plotted for the fixed values |λ| = 0.1 and γ 0 = 1.5, along with ρ varying from 15 to 100.Note that the ranges for |λ| and ρ have been chosen to accommodate the stability conditions for the canonical model [47] (λ > ρ 2 + 2 − ρ ≃ 0.02 for a fixed value of ρ = 50 and ρ > 1/λ − λ/2 = 9.95 for a fixed value of λ = 0.1) as well as the DBI model (126) (λ > γ2 0 ρ 2 + 2γ 0 /γ 0 − ρ ≃ 0.03 for fixed values ρ = 50 and γ 0 = 1.5 and ρ > 1/(γ 0 λ) − λ/2 = 14.95 for fixed values λ = 0.1 and γ 0 = 1.5).Also note that, these ranges are automatically consistent with the stability of the k-inflation model.It is crystal, according to these two figures, that the values of |Σ/H| for non-canonical scalar fields, i.e., the k-inflation and DBI ones, are always smaller than that of the canonical scalar one during the inflationary phase with ρ ≫ |λ|.More interestingly, the value of |Σ/H| for the k-inflation case turns out to be smallest among three types of scalar fields.This result would be useful to investigate the observational constraints of anisotropic inflation due to the two-form field.

B. One-form field vs. two-form field
One might ask about more comparisons among six anisotropy parameters |Σ/H| derived in two types of anisotropic power-law inflation, one is due to the one-form field (a.k.a.vector field) and the other is due to the two-form field, for three different kinds of scalar fields.In this subsection, therefore, we would like to make such comparisons.First, we should note that the derivations of |Σ/H| for one-form anisotropic inflation have been shown in Eq. (2.34) of Ref. [29], Eq. (3.34) of Ref. [30], and Eq.(3.21) of Ref. [35] for the canonical, DBI, and k-inflation scalar fields, respectively.In particular, Fig. 4a contains six curves of |Σ/H| for the fixed values ρ = 50 and γ 0 = 1.5, along with |λ| ∈ [0.06, 0.5].On the other hand, Fig. 4b is plotted for the fixed values |λ| = 0.1 and γ 0 = 1.5, along with ρ ∈ [30,100].Once again, the value ranges of |λ| and ρ have been chosen to accommodate the stability conditions for all six models we have been interested in.According to two Figs.4a and 4b, we observe that the |Σ/H| of k-inflation two-form anisotropic inflation tends to be one of the smallest ones among that of six models when ρ increases or |λ| decreases during the inflationary phase.More interestingly, it turns out that the values of |Σ/H| of canonical and DBI two-form anisotropic inflation are always larger than that of canonical and DBI one-form anisotropic inflation, respectively.In contrast, |Σ/H| of k-inflation two-form anisotropic inflation is always smaller than that of k-inflation one-form anisotropic inflation.All these remarkable points would be useful when judging which anisotropic inflation model is cosmologically viable in the light of observational data.

VI. TENSOR-TO-SCALAR RATIO
In connection with the Planck 2018 data [7,8] as well as future detections like the CMB-S4 project [67], we would like to investigate in this section the corresponding tensor-to-scalar ratio of the present non-canonical anisotropic inflationary models, following our previous works [34][35][36] for the one-form field as well as other works in Refs.[37][38][39][40][41][42] for the one-form field and in Refs.[45,46] for the two-form field.
As discussed earlier, the anisotropy deviation σ should be much smaller than the isotropy parameter α in order to be consistent with the observations of WMAP and Planck.Therefore, it is reasonable to regard the background metric as the spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, rather than the Bianchi type I metric, for simplicity as done in Refs.[34][35][36][37][38][39][40][41][42] for the one-form field as well as Refs.[45,46] for the two-form field.
It is important to mention that when the statistical isotropy of CMB is broken down the scalar power spectrum will be modified, according to Ref. [69], to be where g * is a non-trivial constant characterizing a deviation from the spatial isotropy and is expected to be smaller than one, i.e., |g * | < 1.In addition, P ζ(0) k is nothing but the isotropic scalar power spectrum, while θ is the angle between the wave number k and the privileged direction V close to the ecliptic poles [69].On the observation side, it is worth noting that small values of g * have been identified via several analyses in Refs.[70][71][72][73][74]. On the theoretical side, the analytical formula of g * has been derived within the context of the anisotropic inflation of one-form field in several published papers, e.g., in Refs.[37,46] for the canonical scalar field and in Refs.[34][35][36] for the non-canonical scalar field.More interestingly, the analytical formula of g * for the anisotropic inflation of two-form field has been determined in Refs.[45,46].However, this formula is limited only to the canonical scalar field.We will therefore derive in this section a more general formula of g * for non-canonical scalar fields.

A. Scalar perturbations
Following our previous works done in Refs.[34,36] as well as the seminal paper [63], the metric of scalar perturbations is given by with η ≡ a −1 dt is a conformal time.In order to compute the scalar power spectrum, we adopt the spatially flat gauge with Ψ = 0, which has been used in our previous paper [36] here H is the Hubble parameter.It is noted that this gauge has also been used in many previous papers on anisotropic inflation, e.g., see Refs.[42,45,46].It is also noted that we have used modified notations such as η and ζ, since ζ and η have been used for the power-law solution in the previous sections.
In the absence of non-minimal coupling between the scalar and two-form fields, the isotropic scalar power spectrum P ζ(0) k can be obtained [35,36] which is identical to that firstly obtained in Ref. [63].Here, the notation ' * ' implies the pivot scale (a.k.a.horizon-exit scale), where the spacetime of universe can be approximated as the de Sitter one with a * ≃ −(H * η * ) −1 , while the superscript (0) stands for the de Sitter background.In addition, ϵ is the slow-roll parameter defined as ϵ ≡ − Ḣ/H 2 , while c s is the speed of sound, whose definition is given by [63] where p and ρ the pressure and energy density parameters, respectively.
In the presence of non-minimal coupling between the non-canonical scalar field and the two-form field, the full scalar power spectrum can be calculated, following our previous works [34,36] as well as Refs.[45,46], to be where the correction term δP ζ k is given by (see the Appendix A for derivations) Here, E yz ≡ h/a 3 H 0yz = h/a 3 B ′ yz is the background value of the two-form field [45,46], a ≡ e α = −(ηH) −1 is a scale factor of the background de Sitter spacetime, ′ represents a derivative with respect to the conformal time η, i.e., B ′ yz = dB yz /dη, and h is nothing but the gauge kinetic function.It is noted that the key element to derive the above formula is the standard Bunch-Davies (BD) vacuum state for the non-canonical scalar field, whose definition has been shown in Ref. [75] as This BD vacuum state has also been used in our previous papers [34,36] to compute the CMB imprints of the non-canonical anisotropic inflation based on the one-form field.As a result, the full scalar power spectrum can be expressed as with N csk ≃ 60 is the e-fold number.Therefore, the corresponding g * for non-canonical scalar fields can be defined as where for the canonical scalar field [45,46].Subsequently, the scalar spectral index can be shown as where η ≡ ε/(ϵH) and s ≡ ċs /(c s H) [63,65].Given that the average value of cos 2 θ is ⟨cos 2 θ⟩ = 1/3 [45,46], Eq. (150) now becomes It is important to note that g 0 * and g * have been shown to be negative definite for the one-form field [34][35][36][37]46], while they turn out to be positive definite for the two-form field, according to Eqs. ( 148) and (149).This difference is a very interesting point, which could be useful to distinguish two types of anisotropic inflation, one is due to the one-form field and the other is due to the two-form field, in the light of observational constraints [70][71][72][73][74].

B. Tensor perturbations
In this subsection, we consider tensor perturbations for a model of non-canonical scalar fields non-minimally coupled to the two-form field.It is important to note that tensor perturbations for models of a canonical scalar field nonminimally coupled to either the one-form or two-form field have been investigated in the literature, e.g., see Refs.[45,46].Additionally, tensor perturbations for a model of non-canonical scalar fields non-minimally coupled to the one-form field have been presented in our previous papers [34,36].According to these papers, the metric of tensor perturbations can be expressed as with h ij is the traceless (δ ij h ij = 0) and transverse (∂ i h ij = 0) tensor perturbations satisfying the condition As a result, h ij has two degrees of freedom, h + and h × , which are denoted as polarizations.Since the gravitational sector is not affected by the scalar field as well as the non-minimal coupling between the non-canonical scalar and two-form fields [45,46], then the tensor power spectrum for anisotropic inflation remains identical to that for isotropic inflation, which can be found in Ref. [63], i.e., Consequently, the tensor spectral index reads C. Full tensor-to-scalar ratio Based on the results derived above, we now end up with the full tensor-to-scalar ratio for the two-form case as follows which can be reduced to if we take ⟨cos 2 θ⟩ = 1/3 [45,46].Apparently, this tensor-to-scalar ratio will recover that derived in Ref. [46], in the canonical limit c s → 1.

D. Application to the anisotropic power-law solutions
The tensor-to-scalar ratio is perhaps the most significant parameter of any inflationary model, by which we can judge how viable it is in the light of Planck 2018 data [7,8] or future detections like the CMB-S4 project [67].Therefore, we will consider in this subsection whether the obtained anisotropic power-law inflationary solutions are consistent with the Planck 2018 data.It is noted that the results derived in the previous subsections are general, which could be valid for many types of non-canonical scalar field.
We first consider the k-inflation model, which was firstly mentioned in the one-form case [35] and has been investigated in Sect.III, with the corresponding P (ϕ, X) is given by As a result, the corresponding speed of sound turns out to be Furthermore, c 2 s becomes, for the anisotropic power-law inflationary solution derived in Sect.III, as along with the corresponding slow-roll parameter given by Next, we consider the DBI inflation model with the corresponding form of P (ϕ, X) given by It turns out that the corresponding speed of sound is defined as It is apparent that γ = γ 0 for the power-law inflation as pointed out in Sect.IV, i.e., c 2 s = γ −2 0 .In addition, the corresponding slow-roll parameter reads It is clear that ϵ will not be expressed in terms of c s in this DBI case, unlike the k-inflation case.Due to the fact that η = s = 0 for the power-law solutions, the corresponding scalar spectral index for both models reads where N csk ≃ 60 will be chosen as usual.
In order to compare with the Planck 2018 data [7,8], we are going to plot the n s −r nc diagram.First, we will choose g 0 * = +0.03,according to the analytical formula shown in Eq. ( 149) as well as the recent observational constraints of g 0 * published in Refs.[71][72][73][74] for both k-inflation and DBI models.In addition, we consider a range for the speed of sound as 0 < c s ≤ 0.035 for the k-inflation model.For the DBI inflation model, we choose a range of the ratio λ/ρ as 0 < λ/ρ ≤ 0.055 along with the fixed value for speed of sound is c s ≃ 0.077 (equivalent to γ = γ 0 = 13).According to Fig. 5, it is obvious that the k-inflation model turns out to be more viable, in the light of the Planck 2018 data, than the DBI one.Indeed, the anisotropic power-law inflationary of the k-inflation model is highly consistent with the Planck 2018 data [7,8], similar to the solution found in Ref. [35] for the one-form field.More interestingly, it appears that the tensor-to-scalar ratio of k-inflation one-form field model [35] turns out to be smaller than that of k-inflation two-form field model (see Fig. 6 for the detailed comparison).It should be noted that we have used a more precise formula, c 2 s ≃ −λ/(48ρ), instead of an approximated formula shown in Eq. (5.15) of Ref. [35] in order to plot the n s − r nc diagram of k-inflation one-form field model as displayed in Fig. 6.Remarkably, the tensor-to-scalar ratios of both k-inflation one-form and two-form models seem to be relevant to a forecast of the CMB-S4 project [67].

VII. CONCLUSIONS
We have studied the anisotropic power-law inflation in the presence of the two-form field non-minimally coupled to non-canonical scalar fields.For this purpose, we have focused on examining two typical non-canonical forms of scalar field, which come from the k-inflation [35,62,63] and DBI inflation [30-33, 64, 66].As a result, the anisotropic inflationary solutions derived in these models become stable and attractive as confirmed by the stability analysis based on the dynamical system method.Our present work together with the previous ones by other people in Refs.[45][46][47][48][49] demonstrate that the cosmic no-hair conjecture proposed by Hawking et al. [16] is extensively violated due to the existence of non-minimal coupling between the scalar and two-form fields such as h 2 (ϕ)H µνρ H µνρ .More interestingly, we have pointed out, by taking some simple comparisons, that the non-canonical property of scalar field seems to reduce the magnitude of anisotropy parameter |Σ/H|.This interesting point would be useful when investigating observational constraints of the anisotropic inflation based on two-form field.Additionally, we have made comparisons among six anisotropy parameters |Σ/H| derived in two types of anisotropic power-law inflation, one is due to the one-form field and the other is due to the two-form field, for three different kinds of scalar fields.As a result, the |Σ/H| of k-inflation two-form inflation tends to be one of the smallest ones among six anisotropy parameters.This additional interesting point would also be useful when judging which anisotropic inflation model is cosmologically viable in the light of observational data.In connection with the Planck 2018 data [7,8] and future detections like the CMB-S4 project [67], we have investigated the tensor-to-scalar ratio of the obtained anisotropic power-law solutions.It turns out that the tensor-to-scalar ratio of the anisotropic power-law solution of k-inflation model is more consistent with the Planck 2018 data than that of DBI model.More interestingly, additional analysis has shown that the tensor-to-scalar ratio of k-inflation one-form field model is more viable than that of k-inflation two-form field model in the light of Planck 2018 data [7,8].Furthermore, these tensor-to-scalar ratios seem to be relevant to a forecast of the CMB-S4 project [67].We hope that our study would contribute a new perspective on counterexamples to the cosmic no-hair theorem.It should be noted that the research on the CMB imprints of anisotropic inflation model, in which the canonical scalar field is non-minimally coupled to the two-form field, has been investigated in Refs.[45,46].Our next step would therefore be a study of the effect of non-canonical scalar fields on the CMB imprints of two-form anisotropic inflation like what we have done in Ref. [36].Besides, we will continue to examine whether anisotropic inflation could be caused by other non-minimal couplings between scalar and other fields such as the Yang-Mills one [68].We leave these issues for our future works.
here we have only kept the second-order terms of perturbations.Next, we apply the Fourier transform to δE ij (x, η) such as [45,46] δE ij (x, η) = d 3 k (2π) 3/2 e ik•x δE ij (k, η), (A4) As a result, we can obtain a solution in the super-Hubble regime (|k η| ≪ 1) as [45,46]  ) is the polarization tensor, whose explicit components can be found in Ref. [46].Consequently, the interacting Hamiltonian H ζ can be defined to be [45,46]  which will be identical to that defined in the case of canonical scalar field [45,46] if we take limit c s → 1. Thanks to the expression (A5), which can be written explicitly as [46] δE yz (η we are able to define the following formula, (A12) The integral in the above equation can be calculated as follows [36,45] η −1/(csk)

FIG. 1 :
FIG. 1: Numerical results demonstrate the attractive behavior of the obtained anisotropic fixed point (displayed as a black point).Here, the parameters have been chosen as λ = −0.1 and ρ = 50.
FIG. 4: Comparisons among six anisotropy parameters |Σ/H| derived in two types of anisotropic power-law inflation, one is due to the one-form field and the other is due to the two-form field, for three different kinds of scalar fields.The left figure is plotted for the fixed values, ρ = 50 and γ 0 = 1.5, while the right figure is plotted for the fixed values, |λ| = 0.1 and γ 0 = 1.5.
FIG.6:The n s − r nc diagram of anisotropic power-law inflationary solution of k-inflation one-form field model (lower curve) vs. that of k-inflation two-form field model (upper curve).