Dynamical suppression of spacetime torsion

A surprising feature of our present four dimensional universe is that its evolution appears to be governed solely by spacetime curvature without any noticeable effect of spacetime torsion. In the present paper, we give a possible explanation of this enigma through"cosmological evolution"of spacetime torsion in the backdrop of a higher dimensional braneworld scenario. Our results reveal that the torsion field may be significant at early phase of our universe, but it gradually decreases with expansion of the universe. This leads to a negligible footprint of torsion in our present visible universe. We also show that at early epoch, when the amplitude of the torsion field might not be suppressed, our universe underwent through an inflationary stage having a graceful exit within a finite time. Further in order to link the model with observational constraints, we determine the spectral index for curvature perturbation ($n_s$) and tensor to scalar ratio ($r$) in the present context, which match with the results of $Planck$ 2015 \cite{Planck}.


Introduction
A surprising feature of the present universe is that its large scale behaviour appears to be controlled by one type of geometrical deformation only, namely curvature; while we notice practically no effect of another type of deformation, namely torsion. The most straightforward way of including torsion is to add an antisymmetric component to the connection Γ α µν , which is the essence of the so-called Einstein-Cartan theory [2]. Once torsion enters into the theory in this manner, it can in principle couple with all matter fields having non zero spin. From dimensional argument, it can be easily shown that such interaction terms in general are of dimension 5, and are suppressed by the Planck mass (M p ), just as in the case of graviton couplings. But there has been no experimental evidence of the footprint of spacetime torsion on the present universe. An example is the Gravity Probe B experiment which was designed to estimate the precession of a gyroscope to observe any signature of spacetime torsion [3]. However, all such probes, within the limit of their experimental precision, have consistently produced negative results and thereby disfavored the presence of the torsion in the spacetime geometry of our (3 + 1) dimensional visible universe [4,5,6]. Therefore the apparent torsion free universe indicates that the torsion field, if exists, must be severely suppressed at the present scale of the universe. Thus the question that naturally arises is : why are the effects of spacetime torsion are less perceptible than the spacetime curvature ? There is no satisfactory answer to this in the domain of four dimensional classical gravity models. The proposals to remove torsion by quantum effects in 4 dimensional spacetime have appeared much before in [7], where the authors showed the invisibility of spacetime torsion (on the present energy scale of our universe) through the consideration of quantum corrected (caused by vacuum effects) gravitational action with torsion. There was also attempts to seek an answer to this in the context of higher dimensional braneworld models [8,9,10,11,12,13,14]. In particular, in Randall-Sundrum (RS) scenario [10] which involves one extra compact spacelike dimension with S 1 /Z 2 orbifolding along the extra dimension proposed a possible explanation for this suppression of spacetime torsion in four dimension. This kind of scenario postulates gravity in the five-dimensional bulk, whereas our four-dimensional universe is confined to one of the two 3-branes located at the two orbifold fixed points along the compact dimension. However it has been already shown that a rank-2 antisymmetric tensor field, generally known as Kalb-Ramond (KR) field (B MN ), can act as a source of spacetime torsion where the torsion is identified with rank-3 antisymmetric field strength tensor H MN L having a relation with B MN as H MN L = ∂ [M B N L] [15]. In the RS like scenario where both the gravity and the KR field propagate in the bulk, the exponential warping nature of spacetime geometry causes the KR field (or equivalently the torsion) to be diluted on the visible 3-brane [8,16,17,18]. Also there is a recent work on spacetime torsion with antisymmetric tensor fields in higher curvature gravity model in the context of both four dimensional and five dimensional spacetime [19], where the authors showed that due to the effect of higher curvature term(s), the amplitude of torsion field gets suppressed in the course of the universe evolution. However in the background of cosmological evolution, the suppression of spacetime torsion (on our present universe) sourced by Kalb-Ramond field still awaits a proper understanding. Furthermore one of our authors showed earlier that the amplitude of KR field may be significant and can play a relevant role in the early phase of the universe. This motivates to explore whether the "dynamical evolution" of KR field (from early universe) actually leads to a negligible footprint of torsion on the present universe in the backdrop of braneworld scenario. We also want to explore the "cosmological evolution" of KR field from very early universe to examine whether the universe underwent through an inflationary expansion [20,21,22,23,24,25,26,27]. In particular, the questions that we address in the present paper are : • How does the Kalb-Ramond field evolve from early era of our universe? Does this evolution lead to an explanation of why the effect of torsion is so much weaker than that of curvature on the present visible brane?
• In such circumstance, does the four dimensional universe undergo an accelerating expansion at early epoch? If such an inflationary scenario is allowed, then what is the dependence of the duration of inflation on the KR field energy density? Moreover what are the values of the spectral index (n s ) and tensor to scalar ratio (r) in the present context?
The present paper serves a natural explanation of the above questions in the backdrop of Randall-Sundrum scenario. However the warped RS geometry in its original form is intrinsically unstable due to intervening bulk gravity. A popular way of stabilizing the interbrane separation (also known as modulus or radion) is via Goldberger-Wise (GW) mechanism [28,29] which proposes the existence of a bulk stabilizing scalar field. Some variants of RS model and its modulus stabilization are discussed in [14,30,31,32,33,34,35,36]. Following the GW mechanism, here we propose a dynamical stabilization method of the extra dimensional modulus field (coupled to the KR field through the effective field equations). Our paper is organized as follows: the model is described in section-II, while section-III is reserved for presenting the cosmological field equations and their possible solutions from the perspective of four dimensional effective theory. Their implications and possible consequences are discussed in the remaining part of the paper.

The model
We consider a five dimensional compactified warped geometry two brane model with spacetime torsion in the bulk. In the present context, the source of torsion is taken as rank-2 antisymmetric Kalb-Ramond (KR) field B MN (where latin indices run from 0 to 4). Torsion can be identified with rank-3 antisymmetric field strength tensor H MN L which is related to the KR field as The spacetime is S 1 /Z 2 orbifolded along the extra dimension, where the orbifolded fixed points are identified with two 3-branes. Considering φ as extra dimensional angular coordinate, two branes are located at φ = 0 (hidden brane) and at φ = π (visible brane) respectively while the latter one is identified with the visible universe. One of the crucial aspects of this braneworld scenario is to stabilize the distance between the branes (known as modulus or radion). For this purpose, one needs to generate a suitable modulus potential with a stable minima and in order to do this, here we consider a massive scalar field in the five dimensional bulk. Therefore the action of the model is given by, where R (5) is the five dimensional Ricci scalar formed by the metric G MN , 1 2κ 2 = M 3 (M is the 5 dimensional Planck mass), Λ(< 0) is the bulk cosmological constant and V h , V v are the brane tensions on hidden, visible brane respectively. Ψ is the stabilizing scalar field with m denoting its mass. The KR field action is represented by the last term in the above action. Considering a negligible backreaction of the KR field (B MN ) and the scalar field (Ψ) on the background spacetime, the solution of metric G MN turns out to be same as well known RS model i.e where k = −Λ 24M 3 and r c is the interbrane separation. With this metric, the scalar field equation of motion in the bulk is following, where Ψ is taken as the function of φ only. Considering non-zero value of Ψ on the branes, the above eqn.
(3) has the general solution, where σ = 4 + m 2 k 2 . Further the integrations constants A and B are obtained from the boundary conditions, Ψ(0) = v h and Ψ(π) = v v as follows, Using the five dimensional spacetime metric (see eqn. (2)), different components of stress tensor of the stabilizing scalar field (Ψ) can be obtained as, Putting the bulk scalar field solution (eqn. 3) in the expression of T φφ (Φ) and T µν (Φ) and using the form of A and B in terms of v v and v h (eqn. 5), one can show that the ratio of corresponding component of stress tensor between bulk scalar field and bulk cosmological constant varies as where T φφ (Λ) and T µν (Λ) are different components of stress tensor for the bulk cosmological constant. Similarly the Lagrangian density for Kalb-Ramond field leads to the ratio of KR field stress tensor with the bulk cosmological constant as ∼ H MN L H MN L /M 5 . Thus the stress tensor for the bulk scalar field as well as for the KR field is less than that of the bulk cosmological constant for v 2 v /M 3 and H MN L H MN L /M 5 less than unity. These conditions allow us to neglect the backreaction of the stabilizing scalar field and the KR field (on the background spacetime) in comparison to bulk cosmological constant. In order to introduce the radion field, we consider a fluctuation of branes around the stable configuration (r c ). So, the interbrane separation can be considered as a field (T (x)) and here, for simplicity, we assume that this new field depends only on the brane coordinates. The corresponding metric ansatz is, Correspondingly the introduction of radion field leads to the bulk scalar field (Ψ) solution as follows, where A and B are given by the following expressions, Having these set-up, now we proceed to obtain the effective four dimensional action leading to a viable physical description of our visible universe. In the following few subsections, we individually determine the explicit form of 4D effective action for various parts of the original five dimensional action (eqn.(1)).

Effective action for 5D Einstein-Hilbert term
With the metric in eqn.(6), a Kaluza-Klein reduction for the five dimensional Einstein-Hilbert action reduces to four dimensional effective action as, where R (4) is the four dimensional Ricci scalar formed by the on-brane metric g µν . As it is evident that T (x) is not canonical and thus we redefine the field by the following transformation : In terms of the canonical radion field ξ(x), S (1) ef f takes the following form,

Effective action for bulk scalar field (Ψ): Radion potential
Plugging the bulk scalar field solution (see eqn. (7)) back into the five dimensional scalar field action S scalar = d 4 xdφ and integrating over the extra dimensional coordinate φ yields the effective action as follows, k . Further to derive the above expression, we use the relation between T (x) and ξ(x) as shown in eqn. (9). However it may be noticed that the integrand of eqn.(11) acts as a potential term for the radion field. Afterwards we denote this potential by V (ξ) i.e Eqn. (12) clearly indicates that the potential V (ξ) goes to zero in absence of the bulk scalar field (i.e v h = v v = 0). Therefore as mentioned earlier, the potential term for the radion field is generated entirely due to the presence of the bulk scalar field Ψ [28,29]. The potential in eqn. (12) has a minimum at and a maxima at respectively. Moreover V (ξ) goes to zero as ξ = 0. In figure [1] we give a plot of V (ξ) against ξ. With the expression of < ξ >, we determine the squared mass of the radion field as

Effective action for KR field action
Recall that the 5D KR field action is given by,  (6)) and keeping B 4µ = 0, the above action turns out to be, The Kaluza-Klein decomposition for the KR field can be written as, where B (n) µν (x) and χ (n) (x, φ) represent the nth mode of on-brane KR field and extra dimensional KR wave function respectively. It may be mentioned that the wave function χ (n) is considered to be a function of brane coordinates also (apart from the coordinate φ), this is because our motive is to investigate whether the "dynamical evolution" of KR field leads to its invisibility on our present universe. Substituting the decomposition in the 5-dimensional action S H and integrating over the extra dimension, the four dimensional effective action turns out to be: provided χ (n) (x, φ) satisfies the following equation of motion, along with the normalization condition as, where m n denotes the mass of nth KK mode. As we will see later that χ (n) (x, π) is important to determine the coupling between the KR field and various Standard Model fields on the visible brane. Further eqn. (20) clearly demonstrates that the dynamical evolution of χ (n) (x, φ) is coupled with the modulus (or radion) field T (x).
Eqns. (10), (11) and (19) immediately lead to the final form of the four dimensional effective action as follows : where the radion potential V (ξ) is explicitly shown in eqn. (12). At this stage it deserves mentioning that the zeroth Kaluza-Klein (KK) mode of the field strength tensor (i.e H (0) µνλ ) can be identified with spacetime torsion and thus from now on, we deal with the zeroth mode of the KR field for which m n=0 = 0. With this lowest KK mode, the four dimensional effective action turns out to be, However due to the presence of the potential V (ξ), the radion field acquires a certain dynamics governed by the effective field equations. In this scenario, our motivation is to investigate whether the dynamics of the radion field can trigger such a evolution on the KR wave function χ (0) (x, φ), that will lead to the fact that the effect of KR field (or equivalently the torsion field) evolves to get suppressed with the expansion of our universe. Further in [18], it was shown that the energy density of B (0) µν may be dominant and can have a significant role at early phase of the universe. Therefore it is crucial to explore the dynamical evolution of KR field from very early era of the universe where it is also an intriguing part to examine whether the early universe passes through an inflationary period or not. Motivated by this idea, we try to solve the cosmological Freidmann equations obtained from the four dimensional effective action S ef f . This is demonstrated in the next section. At this stage it deserves mentioning that the energy scale, or the compactification scale, of the five dimensional bulk is ∼ Planck scale. However as mentioned earlier that here we are interested on inflation on our 4D visible universe, where the energy scale (or the inverse of the duration of inflation) comes with ∼ 10 10 GeV which is consistent with the Planck observations as has been described later. Thus the 4D inflationary energy scale is lesser compared to the 5D bulk scale and we can consider the 4D effective action where the extra dimensional component of 5D metric i.e the modulus appears as radion field. Thus the approach here is motivated by the calculation of the effective action proposed by Goldberger and Wise in [28,29].

Effective cosmological equations and their possible solutions
In order to obtain the effective field equations, first we determine the energymomentum tensor for ξ(x) and B µν (x) as, and respectively. The on-brane metric ansatz that fits our purpose is the flat FRW metric i.e where a(t) is the scale factor of the visible universe. However before presenting the field equations, we want to emphasize that due to antisymmetric nature, H µνλ has four independent components on the visible 3-brane, they can be expressed as, With these independent components along with the metric shown in eqn. (26) The above set of equations has the following solution, Using this solution, one easily obtains total energy density and pressure for the where the fields are taken to be homogeneous in space and an overdot denotes d dt ). As a result, the diagonal Friedmann equations take the following form, where H =ȧ a is known as Hubble parameter. Further, the effective field equations for the zeroth mode of KR field (B (0) µν ) and the radion field (ξ) are given by, respectively, where V (ξ) is explicitly shown in eqn. (12). However the only information that we get from eqn. (31) is that the non-zero component of H  29) and (30) immediately lead to an expression as 2Ḣ = −ξ 2 − 1 2 h 4 h 4 . Equating these two expressions ofḢ and using the radion field equation, one finally lands with the following time evolution for h 4 h 4 as, Solving the above differential equation, we obtain where h 0 is an integration constant which is restricted to take only positive values in order to get a real solution of h 4 . Recall that the term 1 2 h 4 h 4 represents the energy density contributed from the KR field i.e ρ KR = 1 2 h 4 h 4 . Therefore eqn. (33) clearly indicates that the energy density of the KR field (zeroth mode) decreases monotonically as the universe expands with time. This leads to a negligible footprint of spacetime torsion on our present visible universe. However at the same time eqn. (33) also demonstrates that the energy density of the KR field should play an important role at early phase of the universe (when a(t) is small compared to the present one). Therefore in order to understand the dynamical suppression of the KR field, it is crucial to determine the time evolution of h 4 from very early universe where it is also important to examine whether the universe undergoes through an inflationary stage or not. To investigate these phenomena, we need to solve the scale factor during initial era. Using the above form of h 4 h 4 (see eqn. (33)), there remain two independent effective field equations, These two equations are sufficient to determine the two unknowns namely the scale factor (a(t)) and the radion field (ξ(t)). As mentioned earlier, we are interested to solve eqns. (34), (35) during early universe and for this purpose, the potential energy of the radion field is considered to be greater than that of the kinetic energy (known as slow-roll approximation) i.e.
Under this approximation, eqn.(34) and eqn. (35) are simplified to, and 3Hξ + ∂V ∂ξ = 0 (38) respectively. Using the explicit form of V (ξ) (see eqn. (12)), we solve the above two equations for ξ(t), a(t) as, and where D = v h vv k 24M 3 σ/2 and ξ 0 , C are integration constants with a 0 = C exp [−ξ 2 0 /8]. Further g 1 (t) has the following form, where 2F 1 symbolizes the hypergeometric function. Similarly the form of g 2 (t) is given by, It may be noticed from eqn. (39) and eqn.(40) that for Ψ → 0 (or v v = 0), the solution of the radion field and the Hubble parameter become and H ∝ 1 a 3 respectively. This is expected because in the absence of bulk scalar field (Ψ), the potential V (ξ) (see eqn. (12)) goes to zero and thus the radion field has no dynamics which in turn makes the variation of the Hubble parameter as H ∝ 1 a 3 (solely due to the KR field having equation of state parameter = 1). Further eqn. (39) clearly indicates that ξ(t) decreases with time. Comparison of eqn. (13) and eqn. (39) reveals that the radion field reaches at its vacuum expectation value (vev) asymptotically (within the slow roll approximation) at large time (t ≫ t 0 ) i.e.
This vev of radion field leads to the stabilized interbrane separation (between Planck and TeV branes) as, where m 2 is squared mass of the stabilizing scalar field Ψ.

Beginning of inflation
After obtaining the solution of a(t) (in eqn. (40)), we can now examine whether this form of scale factor corresponds to an accelerating era of the early universe (i.e. t t 0 ) or not. In order to check this, we expand a(t) in the form of Taylor series (about t = t 0 ) and retain the terms only up to first order in t − t 0 : where a 0 is the value of the scale factor at t = t 0 and related to the integration constant C as, a 0 = C exp [−ξ 2 0 /8]. Eqn.(45) leads to the acceleration of the universe at t → t 0 as follows: It may be noticed that for the condition the early universe undergoes through an accelerating stage while for 3 ,ä(t → t 0 ) becomes less than zero. At this stage, it deserves mentioning that the parameters v v and h 0 controls the strength of the radion field and the KR field energy density respectively. Therefore the interplay between the radion field and the KR field fixes whether the early universe evolves through an accelerating stage or not. However in order to solve the flatness and horizon problems (for a review, we refer to [21,22]), the universe must passes through an accelerating stage at early epoch and from this requirement, here we stick to the condition shown in eqn.(47).

End of inflation and reheating
In the previous section, we show that the very early universe expands with an acceleration and this accelerating stage is termed as the inflationary epoch. In this section, we check whether such acceleration of the scale factor has an end in a finite time or not. The end point of an inflationary era is defined by, We now examine whether this condition is consistent with the field equations shown in eqn. (37) and eqn. (38). Near the end of inflation, one can safely neglect the term proportional to 1/a 6 and thus eqn.(37) takes the following form (at end regime of inflation): Differentiating both sides of this equation with respect to t, we get the time derivative of the Hubble parameter as follows, where we use the equation of the radion field (3Hξ + V ′ (ξ) = 0). Plugging back the expressions of H 2 andḢ into eqn.(48), one gets the following condition on radion field, where t f is the time when the radion field acquires the value 2 √ 2 (in Planckian unit). Eqn. (50) clearly indicates that the inflationary era of the universe continues as long as the radion field remains greater than ξ f (= 2 √ 2). Correspondingly the duration of inflation (i.e. t f − t 0 ) can be calculated from the solution of ξ ( t) as follows, Simplifying the above expression, we obtain σ/2 and σ = m 2 4k 2 . Therefore it is clear that the inflation comes to an end in a finite time. In order to estimate the duration of inflation explicitly, one needs the value of the parameters h 0 , ξ 0 and v v , which can be determined from the expressions of spectral index and tensor to scalar ratio as discussed in the next section. However before moving to the next section, here we discuss the reheating in the present context and the possible effects of KR field (or equivalently the spacetime torsion) on it. Needless to say that reheating describes the production of Standard Model matter at the end of the period of accelerated expansion. For this purpose, we consider an example where the radion field (i.e the inflaton) is coupled to another scalar field ζ, given by the interaction Lagrangian, where g is a dimensionless coupling constant and λ is a mass scale. With this interaction Lagrangian, the decay rate of the inflaton into ζ particles becomes recall that m rad is the mass of the radion field (see eqn. (15)). Generally the energy loss of the inflaton due to the production of ζ particles is taken into account by adding a damping term to the inflaton equation of motion as, Eqn.(54) clearly indicates that the radion field losses energy due to the expansion of the universe and due to transfer to the ζ particles, accounted by the damping terms 3Hξ and Γξ respectively. As a result the production of ζ particles becomes effective when the Hubble parameter becomes less or comparable to Γ, otherwise the energy loss into particles is negligible compared to the energy loss due to the expansion of space as occurred during early phase of the inflation. Therefore the time scale t h (let us call it the reheating time) after when the production of ζ becomes effective is given by With the solution of scale factor (see eqn.(40)), the above equation turns out to be, where g 1 and g 2 are shown in eqns.(41) and (42) respectively. Recall, h 0 represents the energy density of the KR field during early universe and the presence of h 0 in the above expression entails that the KR field indeed affects the reheating time t h . In order to understand the effect of KR field more clearly, h is the reheating time in absence of KR field (h 0 = 0) i.e Thus δt is the deviation of reheating time from t h + δt, we get the following expression of δt where we use eqn.(57) and retain up to the term first order in δt. Clearly δt becomes zero as h 0 → 0, as expected. Using the explicit expressions of g 1 , g 2 along with the condition m k < 1 (i.e ratio of bulk scalar field mass to bulk curvature is less than unity, which is also consistent with Planck observations as described in the next section), we determine the termg 1 −g 2 (sitting in the denominator of eqn.(58)) as follows: Thus the termg 1 −g 2 is positive. As a result, eqn.(58) immediately leads to the condition δt < 0 which in turn makes t h < t (0) h . Thereby the presence of Kalb-Ramond field makes the reheating time lesser in comparison to the case when the KR field is absent. However, this is expected because the KR field corresponds to a deceleration of the universe i.e due to the appearance of KR field the Hubble parameter (H(t)) decreases with a faster rate by which H(t) reaches to Γ more quickly relative to the situation where the KR field is absent.

Spectral index, tensor to scalar ratio and number of e-foldings
In order to test the broad inflationary paradigm as well as particular models against precision observations [1], we need to calculate the value of spectral index(n s ) and tensor to scalar ratio (r) and for this purpose, here we define a dimensionless parameter (known as slow roll parameter) as, Recall the slow roll equation, H 2 = 1 3 V (ξ) + h0 6a 6 . Differentiating both sides of this equation with respect to time, we get where we use the field equation for radion field. These expressions ofḢ and H 2 lead to the slow roll parameter ǫ as follows, The spectral index and tensor to scalar ratio are defined by, With the expression of ǫ obtained in eqn.(61), r and n s turn out to be, and where U 1 and U 2 have the following expressions: respectively. It may be observed that the spectral index and tensor to scalar ratio depend on the parameters v v , h 0 and ξ 0 . To fix these parameters, we use the observational results of Planck 2018 ( combining with BICEP-2 Keck-Array data ) [1] which put a constraint on n s and r as n s = 0.9649 ± 0.0042 and r < 0.064 respectively. Here we take, It may be mentioned that these values of v v and h 0 are consistent with the condition that is necessary for neglecting the backreaction of the bulk scalar field and the KR field on the background five dimensional spacetime. Using eqns.  2] and [3] clearly demonstrate that for 34 < ξ 0 < 38 (in Planckian unit), both the observable quantities n s and r remain within the constraints provided by P lanck 2018 [1]. Further with the estimated values of v v , h 0 and ξ 0 , the duration of inflation (t f − t 0 , see eqn.(51)) comes as 10 −10 (Gev) −1 if the ratio m/k (bulk scalar field mass to bulk curvature ratio) is taken as 0.2 [28]. We also determine the number of e-foldings, defined by N = △t 0 Hdt (△ t = t f − t 0 , duration of inflation), numerically and lands with N ≃ 58 (with ξ 0 = 36, in Planckian unit). In table [1], we now summarize our results: Table [1] clearly indicates that the present model may well explain the infla-Parameters Estimated values n s 0.969 r 0.100 t f − t 0 10 −10 (GeV) −1 N 58 Table 1: Estimated values of various quantities for κv v = √ h0 M 2 ≃ 10 −7 and ξ 0 = 36 tionary scenario of the universe in terms of the observable quantities n s and r as per the results of P lanck 2018. Using the solutions of ξ(t), a(t) (see eqns. (39), (40)) along with the estimated values of the parameters (v v , h 0 , ξ 0 ), we give the plots for the interbrane separation (T (t), see figure [4]) and the deceleration parameter (q = −ä/a, see figure [5]) against a dimensionless time variablet = t t f N .  The equation for the zeroth mode of KR wave function (χ (0) (t, φ)) follows from eqn. (20) and given by, As we may notice that the dynamics of the interbrane separation controls the evolution of χ (0) (t, φ). It may be mentioned that the overlap of χ (0) (t, φ) with the brane φ = π (i.e. χ (0) (t, π)) regulates the coupling strengths between KR field and various Standard Model fields on the visible brane. These interaction terms play the key role to determine the observable signatures of KR field on our universe and thus we are interested to solve eqn.(64) in the vicinity of φ = π (i.e. near the visible brane). Near the regime of φ ≃ π, eqn.(64) can be written as, where χ With this expression, eqn.(65) turns out to be, As it is evident that the left and right hand side of eqn.(66) are functions of time and φ alone respectively. Therefore both sides of eqn.(66) can be separately equated with a constant as follows: and 1 f 2  Fig. [6] clearly demonstrates that in the regime φ ≃ π, the KR wave function monotonically decreases with time and the decaying time scale (t ≃ 25) is less than the exit time of the inflation (t = 55). This may explain why the present universe carries practically no observable signatures of the rank two antisymmetric Kalb-Ramond field (or equivalently the torsion field). Thereby as a whole, the solution of χ where f 1 (t) is obtained in figure [6]. Using this solution as a boundary condition, we solve eqn.(64) (evolution of KR wave function in the whole bulk) numerically as plotted in Fig.[7].
Fig. [7] reveals that the zeroth mode of KR wave function χ (0) (t, φ) decreases with time in the whole five dimensional bulk i.e for 0 ≤ φ ≤ π. However for a fixed t, χ (0) (t, φ) has different values (in Planckian unit) on hidden and visible brane and such hierarchial nature of χ (0) (t, φ) (between the two branes) is controlled by the constant b. For T (t) =< T >, the zeroth mode of KR wave function acquires a constant value throughout the bulk and given by where we use the normalization condition as shown in eqn. (21). This result is also in agreement with [8]. Using the above expression of χ (0) (t, φ) T =<T > , we obtain the coupling strengths of Kalb-Ramond field with U (1) gauge field and fermion field on the visible brane as follows [8]: and where M p = M 3 /k. For k < T >≃ 12 (required for solving the gauge hierarchy problem), e −kπ<T > becomes of the order 10 −16 . Thereby eqns.(70), (71) clearly indicate that the interaction strengths of KR field to the matter fields are heavily suppressed over the usual gravity-matter coupling strength 1/M p . This may well serve as an explanation why the large scale behaviour of our present universe is solely governed by gravity and carries practically no observable footprints of antisymmetric Kalb-Ramond field.

Conclusion
We consider a five dimensional braneworld model with spacetime torsion caused by a rank-2 antisymmetric Kalb-Ramond (KR) field in the bulk. The extra spatial dimension is S 1 /Z 2 orbifolded where the orbifolded fixed points are identified with hidden and visible brane respectively. A massive scalar field is also considered in the bulk in order to generate a stable potential term for the radion field (ξ) -required for stabilizing the interbrane separation. We determine the explicit form of the radion potential (V (ξ)) as shown in eqn. (12). It may be observed that V (ξ) goes to zero in absence of the bulk scalar field, as expected [28]. However the presence of the potential V (ξ) activates a dynamics to the radion field governed by the effective field equations. In this scenario, we want to investigate whether the dynamics of the radion field can trigger such a dynamical evolution on the KR field, that may lead to an explanation of why the effect of torsion is so much weaker than that of curvature on the present visible universe. Motivated by these ideas, we solve the cosmological field equations from the perspective of four dimensional effective theory. Our findings are as follows: • We find that the Kalb-Ramond energy density (ρ KR ) on our visible universe depends on the on-brane scale factor a(t) as ρ KR ∝ 1/a 6 (see eqn. (33)). As we may observe that ρ KR monotonically decreases as the universe expands with time, which leads to a negligible footprint of the KR field on the present universe. However eqn. (33) also entails that the energy density of the KR field may be significant in early universe. This points us to explore the dynamical evolution of the KR field from very early phase of the universe. For this purpose, we solve the coupled Freidmann equations for the radion field (ξ(t)) and the scale factor (a(t)) during initial era and the solutions are given in eqns. (39) and (40)) respectively. It is demonstrated in Fig.[4] that the interbrane separation increases with time and saturates at a constant value (< T >) asymptotically. It is also found that without any fine tuning of the parameters, the asymptotic value of the modulus can address the gauge hierarchy problem. On the other hand, the solution of the scale factor corresponds to an accelerating expansion of the early universe and the rate of expansion depends on the parameters v v and h 0 (with v v and h 0 controls the energy density of the bulk scalar field and the KR field respectively). At this stage, it deserves mentioning that in absence of the bulk scalar field (Ψ), the radion field becomes constant while the Hubble parameter varies as H ∝ 1/a 3 . This is expected because for Ψ → 0 (or v v = 0), the potential V (ξ) goes to zero and thus the radion field has no dynamics which in turn makes the variation of the Hubble parameter as H ∝ 1/a 3 (solely due to the KR field having equation of state parameter = 1). The duration of inflation (t f − t 0 ) is obtained in eqn.(51) which reveals that the accelerating phase of the universe terminates within a finite time. Further we also discuss the possible effects of the KR field on the reheating in the present context. We explained that the presence of Kalb-Ramond field makes the reheating time (the time interval after which the production of new particles becomes effective) lesser in comparison to the case when the KR field is absent. However, this is expected because the KR field corresponds to a deceleration of the universe i.e due to the appearance of KR field the Hubble parameter (H(t)) decreases with a faster rate by which H(t) reaches to Γ (the decay amplitude) more quickly relative to the situation where the KR field is absent.
• In order to test the model with the observations of P lanck 2018 (combining with BICEP-2 Keck-Array data), it is crucial to calculate the spectral index of curvature perturbation (n s ) and tensor to scalar ratio (r), which are defined in terms of the slow-roll parameter (ǫ). Using these definitions, the expressions of n s and r are explicitly determined in the present context and as a result, we find that for suitable values of the parameters (v v , h 0 , ξ 0 ), n s and r remain within the constraints provided by P lanck 2018 [1] (see Table[ 1]). Moreover the duration of inflation comes as 10 −10 (GeV) −1 if the ratio m/k (bulk scalar field mass to bulk curvature ratio) is taken as 0.2 [28].
• However the overlap of the zeroth mode KR wave function (χ (0) (t, φ)) with the visible brane actually fixes the coupling strengths of KR field with various Standard Model fields on the brane. Keeping this in mind, we solve χ (0) (t, φ) on the visible brane, numerically, as plotted in Fig.[6]. It is clearly demonstrated that at φ = π the KR wave function monotonically decreases with time and the decaying time scale is less than the exit time of the inflation. Further we also determine the numerical solution for the KR wave function in the whole bulk (see Fig.[7]), which reveals that the effect of χ (0) (t, φ) decreases with time in the full five dimensional bulk i.e. for 0 ≤ φ ≤ π. However it may be mentioned that the dynamics of χ (0) (t, φ) is controlled by the evolution of the radion field and it turns out that for T (t) =< T >, χ (0) (t, φ) acquires a constant value throughout the bulk as obtained in eqn.(69). Consequently we determine the coupling strengths of KR field with various matter fields on our present visible universe. As a result, such interaction strengths come with a heavily suppressed factor over the usual gravity-matter coupling 1/M p . This may provide a natural explanation why the large scale behaviour of our present universe is solely governed by gravity and carries practically no observable footprints of spacetime torsion.
• The second rank antisymmetric Kalb-Ramond field is related to a pseudoscalar field, known as axion field (Z(x)) given by H µνα = ǫ µναβ ∂ β Z. It may be mentioned that there exist some dark matter models where the axion field was considered as a possible candidate to solve the mystery of dark matter [37,38,39]. However an experimental program named ABRACADABRA is designed to search for axion dark matter and the first results of ABRACADABRA is recently published in [40] where the authors, through estimating the axion-photon coupling, have found no evidence for axion-like cosmic dark matter with 95 percentage C.L. This is consistent with the results of our present paper i.e the present universe carries no evidence of axion-like dark matter coming from Kalb-Ramond field.

Appendix-2
The field equation for the zeroth mode Kalb-Ramond field is given by, where g is the determinant of the on-brane metric. Using the FRW metric ansatz, one obtains √ −g = a 3 (t), where a(t) is the scale factor of the universe. Thus eqn.(72) takes the following form,   Here the third term survives, which ensures that H (0)123 is independent of y.