---
abstract: 'Energy momentum tensor of a conformally coupled quantum scalar field in five dimensional warped cosmological spacetimes is studied. We look at situations where the four dimensional part represents a cosmological thick brane and the scale of the extra dimension is time dependent. Renormalization of the components of the energy momentum tensor is achieved using adiabatic regularization method. The resulting energy and pressure densities explicitly show the effects of warping and the dynamic extra dimension on the created matter. We discussed how the created matter may accumulate to form thick branes along the extra dimension.'
author:
- 'Suman Ghosh[^1]${}^{}$'
title: Renormalized stress tensor of a quantized massless scalar field in warped cosmological braneworld background
---

Introduction
============

Since Kaluza-Klein [@Kaluza; @Klein:1926tv], a large variety of cosmological models with extra dimensions have been proposed over the years. The braneworld models [@Randall:1999ee; @Randall:1999vf; @Gogberashvili:1998vx; @ArkaniHamed:1999hk] where our world is considered as a four dimensional sub-manifold (3-brane) embedded in five dimensions are popular today largely because of their potential in proposing achievable experimental signatures of extra dimensions. The warped type assumes a curved higher dimensional spacetime and the line element on the 3–brane is [*scaled*]{} by a [*warp factor*]{}, thus rendering the higher dimensional metric non-factorisable. Brane-world models seem to provide a viable resolution of the long-standing hierarchy problem in high energy physics. How a quantized field behave in such background is therefore an important question. Analysis of quantum fields in a higher dimensional spacetime with a Kaluza–Klein-like extra dimension has been reported in [@RandjbarDaemi:2000cr; @Saharian:2005xf; @Garriga:1989jx; @Nojiri:2003vk; @Mak:1999tf; @Huang:1989zy]. Gravitational particle production in braneworld cosmology and its implications have been studied in [@Bambi:2007nz]. However, investigations on quantum fields have mostly done in time-independent braneworld backgrounds which do not take into account the expansion of the 3-brane as such. Interestingly large class of warped cosmological i.e. time dependent bulk solutions in presence of various bulk fields were found e.g. in [@Ghosh:2008vc]. In [@Ghosh:2008zs], the author have studied how such various kind of warped and dynamic extra dimensions effect the particle creation rate of a massless bulk scalar field.

Quantum field theory in curved spacetime [@PT; @MW; @Wald-QFTCS; @Buchbinder; @Fulling; @BD] is an useful tool to study particle creation by fields due to evolving geometry of spacetime. It has also been instrumental to analyse the inhomogeneities in the cosmic microwave background and the large-scale structure of the Universe [@Liddle:2000cg; @Hawking:1974sw]. The concept of the so-called adiabatic vacuum gives us a notion of particles in curved space that comes closest to the definition of field quanta in flat space. One important feature of the energy momentum tensor (EMT) of the created particles in curved spacetime are the quadratic and logarithmic ultra-violet (UV) divergences in addition to the expected quartic divergence. Various renormalization methods were developed to tame these infinities. The one we are going to use in this article is adiabatic regularization [@parker:2012; @Parker:1974qw; @Zeldovich:1971mw]. In this method, finite expressions are obtained from the formal one containing UV divergences by subtracting mode by mode (under the integral sign) each term in the adiabatic expansion of the integrand which contains at least one UV divergent part for arbitrary values of the parameters in the theory.

A study of quantum scalar field and renormalization in the context of warped spacetimes where there is a dynamic extra dimension is our objective here. Recent discovery of Higgs particle also motivates us with the following investigation. The plan of the article is as follows. In Section II, we review the conformally coupled scalar field equations and the adiabatic regularization method to find the renormalized EMT (REMT). Here we have essentially used the algorithm introduced in [@Zeldovich:1971mw]. We write the equations to determine particle number density in terms of, which may be called, Zeldovich-Starobinsky (ZS) variables that were introduced in [@Zeldovich:1971mw]. The equations for these ZS variables can be solved as adiabatic series where each term in the series are essentially functions of metric functions and their time derivatives. In Section III, we derive the REMT using the adiabatic subtraction method. We emphasize on how components of REMT depend on metric functions to figure out the distinguishing roles of the warping factor and the time-dependent scale factors as such. Apart from the adiabatic regularization we discuss how this created matter density may be localised along the extra dimension, thus providing a notion of a [*physical*]{} thick brane. Finally, in Section IV, we conclude with comments and future plans.

Quantum field coupled to a spacetime with an extra dimension
============================================================

Let us consider the background line element (using conformal time) to be generically of the form [@Ghosh:2008vc] $$ds^2 = e^{2f(\sigma)} a^2(\eta)[- d\eta^2 + dx^2 + dy^2 + dz^2] + b^2(\eta) d\sigma^2 ,\label{eq:metric}$$ where $ e^{2f(\sigma)}$ is the warp factor, $a(\eta)$ and $b(\eta)$ are the scale factors associated with the ordinary space (${\vec x}$) and the extra dimension ($\sigma$) respectively. $\eta$ denotes conformal time.

The Lagrangian density of a scalar field $\psi(\eta,\vec{x},\sigma)$ coupled to the background geometry is given by [@BD] = - ,      A = 0,1,2,3,4. In case of conformal coupling in $d$ spacetime dimensions, the coupling constant $\xi = \f{d-2}{4(d-1)}$. The corresponding momenta is given by = = a\^2 e\^[2f]{} . Eq. (\[eq:lag\]) leads to the following Klein-Gordon equation, $$\begin{aligned}
& & -\ \frac{e^{-2f}}{a^2} \frac{\partial^2 \psi}{\partial \eta^2} + \frac{e^{-2f}}{a^2} \left(\frac{\partial^2 \psi}{\partial x^2} + \frac{\partial^2 \psi}{\partial y^2} + \frac{\partial^2 \psi}{\partial z^2}\right) + \frac{1}{b^2} \frac{\partial^2 \psi}{\partial \sigma^2}\nonumber \\ 
& & -\ \frac{e^{-2f}}{a^2} \left (\frac{2\dot a}{a} + \frac{\dot b}{b} \right )\frac{\partial \psi}{\partial \eta} + \frac{4f'}{b^2} \frac{\partial \psi}{\partial \sigma} - (m^2 + \xi R)\psi = 0 ,\label{eq:genel} \end{aligned}$$ where $m$ is the mass of the scalar particle, $\xi$ is the conformal coupling constant ($\xi = \frac{3}{16}$ in five dimensions) $R$ is the five-dimensional curvature scalar. A dot ($^.$) denotes differentiation w.r.t $\eta$ and a prime ($\prime$) denotes differentiation w.r.t $\sigma$. For simplicity we choose to investigate a massless field ($m = 0$). To quantize the field, we impose the following commutation relations between the field and its conjugate momenta as usual = 0, = 0, = i \^3([x]{}-[x’]{})(-’). We introduce the creation and annihilation operators using the following mode decomposition as = d\^4q ,      q {,k\_} where $\Phi_q(x^A)$ are the eigen-solutions of the field equation and $d^4q/4\pi^2$ is the integral measure in five dimensional momentum space. This decomposition implies = \[a\_q\^, a\_[q’]{}\^\] = 0, \[a\_q,a\_[q’]{}\^\] = \^4(q-q’). The vacuum state is defined as $a_q|0\> =0$ for all $q$.

To separate variables in the field equation we use the following ansatz: $$\Phi_q (\eta,{\bf x},\sigma) = \frac{1}{N} \phi_q (\eta)F({\bf x})G(\sigma) .
 \label{eq:separating}$$ where $N = e^fa\sq{b}$. The normalization condition for $\psi$ or $\Phi_q$ gives the Wronskian relation, $$\dot \phi_q^*\phi_q - \dot\phi_q\phi_q^* = i  .\label{eq:wronskian}$$ Let us set, $$\frac{1}{F({\bf x})}\left\{ \frac {d^2F({\bf x}) }{dx^2}+ \frac {d^2F({\bf x})}{dy^2} + \frac {d^2F({\bf x})}{dz^2} \right\} = -{\bf k}^2 ,\label{eq:Fx}$$ and $$e^{2f} \left\{\frac {G''(\sigma)}{G(\sigma)} + 2 f'\frac{G'(\sigma)}{G(\sigma)} + \left(\frac{f''}{2} + \frac{3f'^2}{4}\right) \right\} = - k_\sigma^2 .
\label{eq:Gsigma}$$

The above two considerations imply, $$\ddot \phi_q(\eta) + \left[ \left({\bf k}^2 + \frac{a^2}{b^2} k_\sigma^2 \right) + \frac{\ddot a}{8a} - \frac{\ddot b}{8b} + \frac{\dot b^2}{4b^2} - \frac{\dot a\dot b}{4a b} \right]\phi_q(\eta) = 0 .
\label{eq:Phit}$$

One can write Eq.(\[eq:Phit\]). as, $$\ddot \phi_q(\eta) + \left[ \Omega_q^2(\eta) + Q(\eta) \right]\phi_q(\eta) = 0 ,\label{eq:Phit1}$$ where, $$\Omega_q^2(\eta) = \left({\bf k}^2 + \frac{a^2}{b^2} k_\sigma^2 \right)\label{eq:Omgforms}$$ and $$Q(\eta) =\frac{\ddot a}{8a} - \frac{\ddot b}{8b} + \frac{\dot b^2}{4b^2} - \frac{\dot a\dot b}{4a b}.\label{eq:Qforms}$$ Note that for an observer on the brane this scalar field (with a specific $k_\s$), which is massless in the bulk, appears to be a scalar field with mass ‘$k_\s$’ [@Ghosh:2008zs] with an effective on brane scale factor $a(\eta)/b(\eta)$. Thus different $k_\s$ modes, if excited, may be interpreted as particles with different masses on the brane as such. Moving ahead, we note that, Eq. \[eq:Phit1\] admits WKB solutions of the form [@parker:2012; @Zeldovich:1971mw], $$\phi_q(\eta) = \frac{\alpha_q(\eta)}{\sqrt{2\Omega_q}} e_-  + \frac{\beta_q(\eta)}{\sqrt{2\Omega_q}} e_+ , ~~~~\mbox{where}~~~~ e_{\pm} = e^{\pm i\int^\eta \Omega_q d\eta}
\label{eq:WKBsol}$$ with a further restriction, $$\dot\phi_q(\eta) = -i\Omega_q \left[\frac{\alpha_q(\eta)}{\sqrt{2\Omega_q}} e_- - \frac{\beta_q(\eta)}{\sqrt{2\Omega_q}} e_+\right] ,\label{eq:WKBrestriction}$$ where $\alpha_q$ and $\beta_q$ are the Bogoliubov coefficients.

Putting Eqs.(\[eq:WKBsol\]) and (\[eq:WKBrestriction\]) in Eqs.(\[eq:Phit1\]) and using the condition (\[eq:wronskian\]) we get, $$\begin{aligned}
\dot \alpha_q& = &\frac{1}{2} \left( \frac{\dot \Omega_q}{ \Omega_q} - i\frac{Q}{\Omega_q}\right)\beta_q\ e_+ - i \frac{Q}{2 \Omega_q}\alpha_q  ,\nonumber\\
                                             \label {eq:alphabetadot}       \\
\dot \beta_q& = &\frac{1}{2} \left( \frac{\dot \Omega_q}{ \Omega_q} + i\frac{Q}{\Omega_q}\right)\alpha_q\ e_- + i \frac{Q}{2 \Omega_q}\beta_q ,\nonumber\end{aligned}$$ $$\mbox{and,} \hspace{2cm} |\alpha_q|^2 - |\beta_q|^2 = 1 .\label{eq:modsqrdiff}$$

Corresponding initial conditions are, $$\alpha_q (\eta_0) = 1 \hspace{0.5cm} \mbox{and}\hspace{0.5cm} \beta_q(\eta_0) = 0. \label{eq:initial1}$$ where, $\eta_0$ is a suitably chosen initial time of a chosen adiabatic vacuum[^2]. Then the number of particles thus created in mode $q$ is given by, $$N_q = \lim_{\eta\rightarrow \infty} |\beta_q|^2 .\label{eq:partno}$$

Following Zel’dovich and Starobinsky [@Zeldovich:1971mw], Eqs.(\[eq:alphabetadot\]) and (\[eq:modsqrdiff\]), with the initial conditions (\[eq:initial1\]), can be recast in the form, $$\begin{aligned}
\dot s_q & = & \frac{\dot \Omega_q}{2 \Omega_q}u_q + \frac{Q}{2 \Omega_q}\t_q \la{eq:s}, \\
\dot u_q & = & \frac{\dot \Omega_q}{\Omega_q}(1 + 2s_q) - \left[\frac{Q}{\Omega_q} + 2 \Omega_q \right]\t_q \la{eq:u},\\ 
\dot \t_q & = & \frac{Q}{\Omega_q}(1 + 2s_q) + \left[\frac{Q}{\Omega_q} + 2 \Omega_q \right]u_q \la{eq:t},\end{aligned}$$ where, s\_q = |\_q|\^2,    u\_q = 2 Re,    \_q = 2i Im.are the ZS variables which obey the following initial conditions, $$s_q(\eta_0) = u_q(\eta_0) = \t_q(\eta_0) = 0. \label{eq:sut-conds}$$ Thus $s_q(\eta)$ represents the particle number density and is related to the energy density in the $\eta \ra \infty$ or the adiabatic limit. Similarly, $u_q(\eta)$ and $\t_q(\eta)$ are essentially related to the pressure and current density components of the EMT of the created particles[^3].

To get the number of particles created in mode $q$ one has to solve the first order system of differential equations (\[eq:s\]-\[eq:t\]) with initial conditions (\[eq:sut-conds\]) and determine $s_q$ when $\eta\rightarrow \infty $. These equations can be evolved numerically using standard codes, in cases where one is unable to find an analytic solution [@Garriga:1989jx; @Ghosh:2016epo]. In the following, we solve Eqs. \[eq:s\]-\[eq:t\] analytically using the adiabatic approximation. In the adiabatic or quasi classical regime $|\dot\O_q|<<|\O_q|^2$. This means if we expand the ZS variables in adiabatic series as s\_q = \_[r=1]{}\^s\_q\^[(r)]{}, u\_q = \_[r=1]{}\^u\_q\^[(r)]{}, \_q = \_[r=1]{}\^\_q\^[(r)]{}, where, ‘$r$’ denotes the adiabatic order[^4], in the adiabatic limit e.g. $s_q^{(r)}$ is suppressed as $|\dot\O_q / \O_q^2|^r$ and $s_q^{(r+1)}/s_q^{(r)} << 1$. If we replace Eq. \[eq:sut-series\] in Eqs. \[eq:s\], \[eq:u\] and \[eq:t\], we get the following recursion relations among the terms of different adiabatic orders of $s_q$, $u_q$ and $\t_q$, \[eq:s\]    s\_q\^[(r)]{} &=& d( u\_q\^[(r)]{} + \_q\^[(r-1)]{}),\
\[eq:u\]    \_q\^[(r+1)]{} &=& s\_q\^[(r)]{} - u\_q\^[(r)]{} - \_q\^[(r-1)]{},\
\[eq:t\]    u\_q\^[(r+2)]{} &=& - (2s\_q\^[(r)]{}+u\_q\^[(r)]{}), where $r = 2,4,..$ and \_q\^[(1)]{} = ,    s\_q\^[(1)]{}=u\_q\^[(1)]{}=\_q\^[(2)]{}=0 s\_q\^[(r)]{}=u\_q\^[(r)]{}=\_q\^[(r+1)]{}=0      r. Using these results back in Eqs \[eq:s-ad\], \[eq:u-ad\] and \[eq:t-ad\], we get, u\_q\^[(2)]{} = - = (-2 + - 2Q) s\_q\^[(2)]{} = . The higher order terms can be derived in a straightforward manner, upto arbitrary order, from Eqs \[eq:s-ad\], \[eq:u-ad\] and \[eq:t-ad\] (see Appendix \[app:sut\]).

Energy momentum tensor
======================

The energy momentum tensor of a conformally coupled field is given by [@Zeldovich:1971mw; @Flanagan:1996gw] T\_[AB]{} = \_A\^\* \_B- + (G\_[AB]{} + g\_[AB]{}- \_A\_B)||\^2 where $G_{AB}$ is the Einstein tensor and $\square = g^{AB}\na_A\na_B$. For metric (\[eq:metric\]) with separation of variables given by Eq. (\[eq:separating\]), the vacuum expectation value of the 00-component of the stress tensor (with respect to the adiabatic vacuum ) i.e. the energy density created during the evolution of warped cosmological braneworlds is given by, T\_[00]{} = &&\
= && |G|\^2

where we have defined = ( - ),   \_0Ø\_q\^2() = k\^2 + K\_0\^2(k\_;), where K\_0\^2(k\_;) = {|G’|\^2 - \_|G|\^2 - ()|G|\^2 - \_\^2 |G|\^2}

Using equations \[eq:WKBsol\], \[eq:WKBrestriction\] and \[eq:sut-def\], Eq. \[eq:T00\] can be written as T\_[00]{} = |G|\^2 && Note that the $1/2$-factors appearing in Eq. \[eq:T00\_sut\] (and in the following equations) represent the vacuum contributions (which will later be subtracted for regularization). The above expression for energy density does contain features of a dynamic warped extra dimension. In the limit $a(\eta)=b(\eta)$ and $f(\s)=0$, we have $\g=0$ and $\O_q^2 = {}_0\O_q^2 = k^2 + k_\s^2$, which would make all the terms in Eq. \[eq:T00\_sut\] vanish to recover the standard result that no particles would be created from a massless scalar field in a conformally flat spacetime.

Similarly, the other components of the energy momentum tensor can be written, e.g. the pressure density along the on-brane spatial directions, as T\_[ii]{} = && |G|\^2\
= && |G|\^2

where, \_iØ\_q\^2() = k\^2 + K\_i\^2(k\_;) and K\_i\^2(k\_;) = {|G’|\^2 + \_|G|\^2 + ()|G|\^2 - \_\^2 |G|\^2}

The pressure density along the extra spatial direction ‘$\s$’ is given by T\_ = &&\
= &&

where, \_Ø\_q\^2() = k\^2 + K\_\^2(k\_;) and K\_\^2(k\_;) = {|G’|\^2 + \_|G|\^2 + |G|\^2}.

The only off-diagonal element of EMT, which is a distinguishing feature of such warped spacetimes is given by T\_[0]{} &=& |G|\^2\
&=& |G|\^2 where K\_[0]{}\^2(k\_;) = (\_|G|\^2 - 3f’|G|\^2)

After getting the expression of bulk EMT we note that these components of bulk EMT are functions of both time and the extra dimension as expected. In the following we explore what these dependencies physically mean.

Localisation of matter density
------------------------------

The functional dependence of the EMT components on time and extra dimension can be studied as a two dimensional surface in $T,t,\s$ space. So, these components can be analyzed graphically. However, these surfaces, in general, will have complicated shapes. Further, it would be difficult to carry out the integration over $k_\s$ to find how $\< T_{AB}\> $ (or its renormalized version) varies along $\s$. Note that, the $\s$-dependence in $\< T_{00}\> $, for example, comes from an over all factor of $|G|/e^f$ and ${}_0\O_q(\s)$ (which carries $K_0^2(k_\s,\s)$). For simplicity, let us look how these factors depend on $\s$, for few fixed values of $k_\s$. This will let us compare, e.g. how much matter is created at different location along the extra dimension due to the bulk scalar field. We can do this even before the adiabatic subtraction as regularization will only effect Z-S variables which do not depend on $\s$.

Firstly, we look at what happens to the factors $K_A^2(k_\s,\s)$ at the location of the brane. Eq. \[eq:Gsigma\], near $\s=0$, in presence of a growing warp factor, reduces to G”() + ( + k\_\^2)G() = 0, i.e. $G(\s)$ behaves like a plane wave around $\s=0$. Eq. \[eq:Gnear0\] further implies that, near $\s=0$ K\_0\^2 = K\_i\^2 = k\_\^2 + ,      K\_\^2 = k\_\^2 + . Thus near $\s=0$, ${}_A\O_q(\s)$ may be considered as independent of $\s$. Similar $\s$-dependence (or independence) is found for $\< T_{ii}\> $ as the pre-factor is same. In $\< T_{\s\s}\> $, there is an extra factor of $e^{-2f}$ present which does not change any qualitative behavior of the pre-factor in any significant manner near $\s=0$. This simple analysis tells us it is enough to analyze the pre-factor in order to understand the spatial variation of EMT components.

In the following, we have graphically presented the behaviour of the pre-factor in the integrand of $\< T_{00}\> $ i.e. $|G|^2e^{-2f}$ as a function of $\s$ (for suitable lowest values of $k_\s$). In order to do that, we have solved Eq. \[eq:Gsigma\] numerically (with initial conditions $G(0)=1$ and $G'(0)=0.1$) for two different functional forms of the warp factor, namely $f(\s)=\pm \log(\cosh \k\s)$, where $1/\k$ represents a length scale along the extra dimension and the ‘+’ sign (‘-’ sign) represents the so-called growing (decaying) warp factor in thick brane models. We have set $\k=1$ for all numerical computations.

![Left: Variation of the overall pre-factor in presence of decaying warp factor for $k_\s = 0, 0.5, 1$. Right: Variation of the overall pre-factor in presence of growing warp factor for $k_\s = 0, 10, 20$.[]{data-label="fig:case1and2"}](case1.pdf "fig:"){width="2.5in"} ![Left: Variation of the overall pre-factor in presence of decaying warp factor for $k_\s = 0, 0.5, 1$. Right: Variation of the overall pre-factor in presence of growing warp factor for $k_\s = 0, 10, 20$.[]{data-label="fig:case1and2"}](case2.pdf "fig:"){width="2.5in"}

The numerical values of $k_\s$ are chosen[^5] such that the curves show one, three and five minima or maxima in Fig. \[fig:case1and2\]. The left graph with the decaying warp factor implies that most of the energy density is created away from the brane location ($\s=0$). Further, density increases indefinitely in the $\s \ra \pm \infty$ limit. Thus this scenario is not physically interesting as it may only be interpreted as creation of two [*physical*]{} branes, formed due to accumulation of infinite amount of matter, separated by infinite distance. On the other hand, the right graph with the growing warp factor implies that most of the created energy density is localised near $\s=0$. Thus a thick brane is formed with localised matter with masses $k_\s$. Thickness of this brane may be quantified by the ‘width’ of the envelope of all the pre-factor curves for all $k_\s$ modes (the zero-mode curve itself looks like such an envelope!). This scenario is physically more appealing. Therefore we carry on our analysis on braneworlds with growing warp factor only. Note that the central maxima for non-zero $k_\s$ modes lie on or near $\s=0$, however, the zero mode seems to be little off. The location of this peak depends on the initial conditions used. To verify how sensitive the location of central maxima is to initial conditions, we have retraced these curves with significantly different initial conditions namely with $G'(0)=1$ (left) and $G'(0)=-1$ (right) in Fig. \[fig:case21and22\].

![Left: The overall pre-factor in presence of growing warp factor for $G(0)=1$ and $G'(0)=1$. Right: The overall pre-factor in presence of growing warp factor for $G(0)=1$ and $G'(0)=-1$.[]{data-label="fig:case21and22"}](case21.pdf "fig:"){width="2.5in"} ![Left: The overall pre-factor in presence of growing warp factor for $G(0)=1$ and $G'(0)=1$. Right: The overall pre-factor in presence of growing warp factor for $G(0)=1$ and $G'(0)=-1$.[]{data-label="fig:case21and22"}](case22.pdf "fig:"){width="2.5in"}

Fig. \[fig:case21and22\] clearly shows that though the zero mode can peak away from the so-called brane location, localisation of the non-zero modes are insensitive to the initial conditions. So we may conclude that localisation of massive particles on the brane is a generic feature of braneworld with a growing warp factor [^6].

Adiabatic regularization of bulk EMT
------------------------------------

Let us now turn to the regularization of the components of $\< T_{AB}\> $ which obviously contains ultra-violet divergences, coming from the large modes, apart from the usual Minkowski space divergence. According to the adiabatic subtraction algorithm, the REMT is found after subtracting the vacuum contributions and the terms upto the fourth adiabatic order from the adiabatic series expansion of EMT. Note that as the adiabatic regularization is done mode by mode, components of REMT of different adiabatic order and of each mode totally decouple from each other and are covariantly conserved by definition.

In our case, let us define the on-brane [*conformal*]{} Hubble parameter to be $h = \dot{a}/{a}$. Note that, as long as $a k_{\s}/b << h$, the adiabatic corrections (which go as powers of $|\dot\O_q / \O_q^2|$) are much smaller than the zeroth order ones even for the $k = 0$ mode making the adiabatic expansion legitimate for all $k$-modes. On the other hand, for $k_{\s}=0$ mode, the adiabatic expansion fails for modes with $k < h$ (i.e. for on-brane superhorizon modes). However, the expansion is still valid for the mode functions satisfying $k >> h$ (i.e. for on-brane subhorizon modes). Thus the momentum integrals can be performed in the interval ($k_* , \infty$), for some $k_*>>h$. Since the divergences are essentially caused by large modes we expect the adiabatic regularization to be valid in this case. Subtracting terms upto fourth adiabatic order from the diverging energy density, we get the renormalized energy density to be, T\_[00]{} \_[Ren]{} &=& |G|\^2\
& & - |G|\^2 . In the last line we have approximated the renormalized quantity by neglecting terms of order higher than the sixth in the adiabatic series (see Appendix \[app:sut\]). Similarly, the other REMT components can be written as

T\_[ii]{} \_[Ren]{} - && |G|\^2

T\_ \_[Ren]{} - &&

T\_[0]{} \_[Ren]{} && |G|\^2 Using the expressions of $s_q^{(6)}$, $u_q^{(6)}$ and $\t^{(5)}$ as given in Appendix \[app:sut\], one gets analytic expressions for all the non-vanishing components of the REMT of a massless bulk quantized scalar field in cosmological braneworld background.

The projected renormalised energy density on the cosmological brane at $\s=0$ (in case of growing warp factor) is given by T\_[00]{} \_[Ren]{}\^[(brane)]{} & & - . Similarly, the renormalized pressure densities are T\_[ii]{} \_[Ren]{}\^[(brane)]{} & & - . In the conformally flat limit (which implies $k_i=k_\s$) the above expressions leads to usual five dimensional renormalised EMT for massless scalar field. To get physically interesting outputs from the above complicated expressions one needs to look at specific examples of viable cosmological braneworld scenarios as such. A detailed study of such various models will be reported elsewhere.

Discussion
==========

In this article we have used the adiabatic regularization method to study the effect of particle creation due to a massless bulk scalar field in a warped cosmological braneworld scenario. Our main goal was to uncover the distinguishing effects of the key components of the geometry of such models, i.e. the warping factor, the cosmological expansion factor and the dynamic scale of the extra dimension. It is found that the warping factor plays key role in deciding how much of the created matter would be localised at different locations along the extra dimension, thus providing us a guiding principle to locate where a four dimensional thick brane may be formed. In our case, it is found that presence of a decaying warp factor is not suitable for such matter brane to form. On the other hand, a growing warp factor helps matter to accumulate near the minimum of the warp factor which may be considered as the four dimensional thick brane we live in. This further suggests that as many parallel branewrolds will be formed as the number of minima in the warp factor.

Adiabatic regularization is performed on the EMT by subtracting diverging terms upto fourth adiabatic order from the formal vacuum EMT. This gives us a finite bulk REMT which carry the effect of dynamic nature of the extra dimension. It would be an interesting exercise to compare the matter density projected on the brane with the background on-brane critical energy density to test the viability of various models. However, it will depend on the exact dynamic nature of $a(\eta)$ and $b(\eta)$. Note that, for brane dynamics the ratio $a(\eta)/b(\eta)$ plays the role of effective cosmological scale factor. A large class of cosmological braneworld solutions with various combinations of $a(\eta)$ and $b(\eta)$ were found in [@Ghosh:2008vc]. Based on the formalism developed here, a detailed quantitative study of how adiabatically regularized EMT components depend on these expansion factors and effects the braneworld cosmology will be reported in a separate article.

Recently renormalization of EMT of spin 1/2 field is achieved under adiabatic regularization scheme [@Landete:2013axa; @Landete:2013lpa; @delRio:2014cha; @Ghosh:2015mva; @Ghosh:2016epo]. Few studies on fermionic fields in non-factorisable geometries, e.g. in Randall-Sundrum braneworld background has been reported in literature [@Ichinose:2002kg; @Grossman:1999ra; @Koley:2004at]. It will be interesting to carry out an analysis on quantum spinor fields in warped cosmological background, as we have done for a bulk scalar field in this article. We hope to report on that in a future communication.

Acknowledgements {#acknowledgements .unnumbered}
================

This research is supported by a start-up grant awarded by University Grant Commission, India, with grant number no. F.30-420/2018(BSR). Author is thankful to Prof. Sayan Kar for insightful discussions and the ‘Visitors Program’ of the Centre for Theoretical Studies at the Indian Institute of Technology Kharagpur where the initial part of this work was done.

Useful geometric quantities
===========================

Affine connections: \^0\_[00]{} = \^0\_[ii]{} = \^i\_[i0]{} = ,    \^i\_[i]{} = \^\_[0]{} = f’;      i=1,2,3. \^\_[00]{} = - \^\_[ii]{} = f’ ,     \^0\_ =

   R = (6 + 4 + 2 ) - ,

$s_q$, $u_q$ and $\t_q$ of different adiabatic orders
=====================================================

\_q\^[(3)]{} = u\_q\^[(4)]{} &=& s\_q\^[(4)]{} &=& \_q\^[(5)]{} &=& u\_q\^[(6)]{} &=& s\_q\^[(6)]{} &=&

[^1]: Electronic address : [*suman.ghosh@igntu.ac.in*]{}

[^2]: e.g. one may choose $a(\eta_0) = b(\eta_0) = 1$.

[^3]: Note that, in case of EMT of a spin-$\f{1}{2}$ field, corresponding axial anomaly or the axial current in $m\ra 0$ limit is equal to $\t_q$ integrated over all the modes.

[^4]: An adiabatic order essentially equals the number of time derivative of $a(\eta)$ contained in the term.

[^5]: Allowed values of $k_\s$ are discrete in case of finite extra dimension [@Ghosh:2008zs].

[^6]: Note that, earlier, [*classical*]{} localisation of massive particles in such spacetimes has been demonstrated through geodesic confinement in [@Ghosh:2009ig]
