BRST symmetry for Regge-Teitelboim based minisuperspace models

The Einstein-Hilbert action in the context of Higher derivative theories is considered for finding out their BRST symmetries. Being a constraint system, the model is transformed in the minisuperspace language with the FRLW background and the gauge symmetries are explored. Exploiting the first order formalism developed by Banerjee et. al. the diffeomorphism symmetry is extracted. From the general form of the gauge transformations of the field, the analogous BRST transformations are calculated. The effective Lagrangian is constructed by considering two gauge fixing conditions. Further, the BRST (conserved) charge is computed which plays an important role in defining the physical states from the total Hilbert space of states. The finite field dependent BRST (FFBRST) formulation is also studied in this context where the Jacobian for functional measure is illustrated specifically.


Introduction
In field theory, higher derivative terms are often considered to describe renormalisation of the corresponding theory. In [1] Stelle showed that higher derivative gravity models can be renormalised. Being a simple model, the Einstein -Hilbert(EH) action, when expressed with respect to the metric components shows the higher derivative nature. A careful examination of the action shows that one can take out a total derivative term which is known as the Gibbons-Hawking term [2]. This makes the Einstein-Hilbert theories a first order one. Although, in usual theories the surface terms are neglected, but in gravitation they serve as important information about entropy of the system. Therefore, removing the surface term trivially for the gravity models may lead us to losing important information about the thermodynamics of the system. This motivates us to keep the total derivative term intact to the action thereby making it a higher derivative theory.
In the minisuperspace version one considers finite degrees of freedom of a model to study the full theory. With Friedmann-Rebertson-Lamaitre-Walker (FRLW) background we construct the minisuperspace version of the EH action which is analogous to the minisuperspace version of the Regge-Teitelboim (RT) model [3]. Hamiltonian analysis of this minisuperspace model has already been performed in different ways and quantization was further explored in [4][5][6][7][8]. Particularly, in [4] the same model with the higher derivative term was analysed without neglecting the surface term. There the authors have explored the reparammetrisation symmetries of the model exploiting the first order formalism developed earlier in [9]. To our knowledge, the BRST analysis for this model is still lacking in the literature which can be important from the view of quantization of the gravity models. Therefore, this paper will serve the purpose visualizing the role of BRST symmetries in the cosmological models which will serve important for the purely higher derivative cosmological models.
The BRST symmetry and the associated concept of BRST cohomology provide the most used covariant quantization method for constrained systems such as gauge [10], string and gravity theories [11,12]. There appear geometrical constraints, for instance in topological solitons such as nonlinear sigma model, CP(N) model, Skyrmion model and chiral bag model, which can be rigorously treated in the Hamiltonian quantization scheme in extended phase spaces [13]. The extended phase space includes ghost and antighost fields of the theories. This symmetry guarantees the quantization, renormalizability, unitarity and other aspects of the first-class constrained theories. The derivation of Slavnov-Taylor identities utilize the BRST symmetry transformation. The generalization of BRST transformation by allowing the parameter finite and field-dependent, known as FFBRST symmetry [14], are studied extensively for diverse area [15][16][17][18][19][20].
In the present paper, the minisuperspace version of the EH action is considered to study its BRST symmetries. For that we first perform the Hamiltonian formulation exploiting the first order formalism developed in [9]. Dirac's method for constraint analysis is followed to find out the full constraint structure of the model [11,21,22]. The theory is found to have only one primary first class constraint and one secondary first class constraint. Existence of the primary first class constraint is in conformity with one independent gauge symmetry which is in the course of calculation identified as reparametrisation invariance. To canonically quantize the model two gauges were proposed in [4]. We considered the same gauge fixing conditions to find out their implications in the BRST symmetries. We found out the nilpotent symmetries which are in line with the standard characteristics of the BRST symmetries. Also the BRST charge is calculated using Noether charge definition. Further the BRST symmetry is generalized by integrating the transformation parameter after making it field-dependent. By doing so we obtain a new symmetry of the theory. The main feature of this symmetry is that this leads to a non-trivial Jacobian for functional measure under change of variables, however the effective action remains unchanged. We compute the Jacobian systematically for an arbitrary and for some specific choices. We observe that such analysis amounts certain change in the effective action. Since this change appears in the BRST exact part, so it is not harmful at all for the theory on physical ground.
The paper is organized as follows. In section II we give a very brief introduction to the RT model. In the section III we study the gauge symmetries using the first order formalism. Section IV is our main contribution which we devoted to explore the BRST symmetries. In this section we calculated the BRST transformation using the same gauges as considered in [4]. Also, we calculated the BRST charge and the Jacobian of the path integral under FFBRST symmetry in section V. Finally, we conclude in section VI.

Minisuperspace version of the Einstein-Hilbert action
In this section we convert the Einstein-Hilbert action to the minisuperspace version keeping the FRLW background. As mentioned earlier, this has a direct representation as the Regge-Teitelboim model [see [4,7]]. As the Universe has degrees of freedom, limiting the degrees of freedom can make life easier. So we reduce the dimension of spacetime which is the underlying concept for minisuperspace. Let us take a d-dimensional brane Σ evolving in a N dimensional bulk spacetime with fixed Minkowski metric η µν . The world volume swept out by the brane is a d + 1 dimensional manifold m defined by the embedding x µ = X µ (ξ a ) where x µ are the local coordinates of the background spacetime and ξ a are local coordinates for m. The theory is given by the action functional Clearly, in this action integral β is should be having the dimension [L] 1−d . Here g is the determinant of the metric, Λ is the cosmological constant and R, the Ricci scalar.
Next, we take a homogeneous and isotropic universe in 4-dimensions and embed it locally in a 5-dimensional Minkowski spacetime. The reason behind taking this choice of dimension is worth mentioning. Algebraically, the embedding of a 4 dimensional space in a 5D Minkowski space is correct but on the physical ground it has a limitation only to the vacuum solutions with cosmological constant. This is in line with the findings of Kasner who first pointed out that most curved solutions in 4D can not be embedded in 5 dimensions [25]. On the other hand, we can embed a all 4 dimensions with vacuum solutions to 6 or higher dimensions. This helps us to understand the dynamics of 4D space with the help of extra dimension. Now we decompose the Lagrangian in the action (1) in the ADM formalism [23]. In (4+1) dimensions if we take the normal vector (N = ṫ 2 −ȧ 2 is the lapse function and the over dot ' . ' is derivative with respect to τ ) the induced metric on the world volume is, We will consider the above metric for the analysis of the system. The system is parametrise by τ . For a spherically symmetric system only t(τ ) and a(τ ) are the dynamical variables and other variables can be integrate out. In this situation, the Ricci scalar is found out to be From (1) the Lagrangian takes the form in 4 dimensions as [4,7] 1 L(a,ȧ,ä,ṫ,ẗ) = aṫ Note that the Lagrangian (5) contains higher derivative terms of the field a. However we can write it as [7] It is worthy to note that the boundary term in this case plays no role in the equations of motion. However being a gravitational system, it is better to keep it at the first place as it carries important information about entropy. Also this will enable us to understand a HD gravity theory. For the Hamiltonian analysis of the system we will consider the HD version i.e. (5).

Hamiltonian analysis
tmp In this section we will perform the Hamiltonian formulation to extract the inherent symmetries of the higher derivative system (5). Particularly, we will follow the method of first order formalism developed in [9]. This is a different approach than the well known Ostrogradski method. We give results rather shortly as details can be gathered from [4].
To begin with, we convert the HD Lagrangian (5) into a first order Lagrangian by defining the time derivative of the fields as new fields. Thus in the equivalent first order formalism, we define the new fields as,ȧ These definitions introduce new constraints in the system given by Enforcing the constraints (8) through the Lagrange multipliers λ a , and λ t as we get the first order Lagrangian, known as the auxiliary Lagrangian. The Euler Lagrange equation of motion, obtained from the first order Lagrangian (9), by varying w.r.t. a, A, t, T, λ a and λ t , are respectively given by To perform the Hamiltonian formulation we notice that the phase space is spanned by q µ = a, t, A, T, λ a , λ t and their corresponding momenta as Π qµ = Π a , Π t , Π A , Π T , Π λa , Π λt with µ = 0, 1, 2, 3, 4, 5.
Here the momenta is defined by This is the point of departure of our Hamiltonian formulation from the Ostrogradsky convention of [7].
Using the definition of momenta we find out the following set of primary constraints, To obtain the full set of primary first class constraints we consider the constraint combination Computing the Poisson brackets we find that only Φ ′ 3 gives the zero Poisson brackets with all the constraints. The nonzero Poisson brackets between the newly defined primary set of constraints Now, expression for the canonical Hamiltonian is given by while the total Hamiltonian becomes Here Λ 1 , Λ 2 , Λ 3 , Λ 4 , Λ 5 , Λ 6 are undetermined Lagrange multipliers. Demanding the time evolution of the constraints Φ 1 , Φ 5 , Φ 6 to be zero ({Φ i , H T } ≈ 0) the following Lagrange multipliers get fixed whereas, conservation of Φ 2 gives the following condition between Λ 4 and Λ 6 Time preservation of the constraint Φ ′ 3 gives rise to the following secondary constraint Likewise, Φ 4 yields the following secondary constraint Time preservation of Ψ 1 trivially gives 0 = 0. A similar analysis involving Ψ 2 yields, on exploiting (22), , .
The iterative procedure is thus closed and no more secondary constraints or other relations are generated.
Finally, after manipulation of the constraints we arrive to the following set of constraints Here {F 1 , F 2 } are the pair of first class constraints and {S 1 , S 2 } are the pair of send class constraints. F 1 is the primary first class constraint which is consistent with the fact that there is only one undetermined multiplier in the theory. For removal of the second class constraints the Poisson brackets will be replaced by the Dirac brackets which is defined by where we have taken with ǫ 12 = 1 and i, j = 1, 2.
We calculate the Dirac brackets between the basic fields which are given below (only the nonzero brackets are listed) , , .
Due to the Dirac brackets all the second class constraints become strongly zero. So in this theory there is only two first class constraints. Now we proceed to show the gauge symmetry of the system.

Construction of the gauge generator
To construct the BRST symmetries, it is important to find out the gauge symmetries of the theory. So, we construct the gauge generator which is the linear combination of all the first class constraints given by, Here {Φ a } is the whole set of constraints and ǫ a , a = 1, 2 are the gauge parameters. However not all the gauge parameters ǫ a are independent. The independent gauge symmetries can be found out by applying the algorithm developed in [24]. In addition, one should careful about finding out the independent symmetries since it is a HD theory. For the HD theories the commutativity of gauge variation and time translation i.e.
where q n,α denotes the α-th order time derivative of q, can play very important role in finding out the independent gauge parameters [9]. Now, to find out the independent gauge parameters we use [24] which gives us So here ǫ 2 may be chosen as the independent gauge parameter. This means that there is only one independent gauge transformation which essentially is in conformation with the fact that there is only one independent primary first class constraints.
The gauge transformations of the fields are given by This completes our analysis for finding out the gauge symmetries of the system. We have seen that the RT model has one independent primary first class constraint which is consistent with the fact that there is only one independent gauge symmetry.

BRST quantization
As we have seen that the RT model possesses a redundant degrees of freedom and has been identified as the diffeomorphism symmetry in [4]. The principle of gauge invariance is essential to constructing a workable RT model. But it is not feasible to perform a perturbative calculation without first fixing the gauge i.e. adding terms to the Lagrangian density of the action principle which break the gauge symmetry to suppress these unphysical degrees of freedom. Since there are two first class constraints, so we take two well motivated gauge fixing term and The first gauge is the very well known cosmic gauge and the second gauge was proposed in [4]. Here α is some number with the condition α = H. With these choice of gauges the first class constraints become second class. It is also known that such gauge fixing condition will induce ghost terms for the model which play s an important role in the proof of unitarity of the theory. Now we perform the BRST symmetry analysis of this model.

Gauge fixing and Faddeev-Popov action
For simplicity we rename the gauge parameter ǫ 2 as ǫ. From (32) we can see that the other gauge parameter becomes −Λ 3 ǫ −ǫ. The above gauge conditions can be incorporated at quantum level by adding the following term to the classical action: here both the first class constraints are incorporated by a single Lagrange multiplier λ as there is only one independent gauge symmetry. To construct the ghost term in the effective Lagrangian we first compute the gauge variation of the gauges ϕ 1 and ϕ 2 . They are given by So we can easily calculate the ghost term for the effective Lagrangian as where c is a ghost field andc is an anti-ghost field. These ghost and anti-ghost fields have a geometrical interpretation Maurer-Cartan forms on diffeomorphism. By adding the gauge fixing and ghost terms to the classical Lagrangian, the effective Lagrangian is then given by Following the structure of the gauge transformations (33 -36) of the fields, we are able to construct the BRST transformations as which leaves the above effective action invariant. Here η is the Grassmann parameter of transformation. It is easy to check that the above transformation is nilpotent, i.e. δ 2 b = 0. This transformations can be used to compute the Ward identities which will yield the relation between different Green's function.
Here we are also able to define the anti-BRST transformation just by replacing the ghost by anti-ghost field and vice-versa.

BRST charge calculation
In this subsection we derive the total BRST charge corresponding to the above BRST symmetry. Utilizing the Noether's formula, we calculate the BRST charges for the different fields as , Total BRST charge, which is a hermitian operator, is given by .
This operator operator Q 2 which implements the BRST symmetry in quantum Hilbert space should be nilpotent which is evident from the above expression. In order to get probabilistic interpretation of the model we must project out all the physical states in the positive definite Hilbert space. Now, this charge (Q) annihilates the physical states of the total Hilbert space as follows: which helps in defining the physical states in total Hilbert space of the theory.
Similarly, it is also possible to compute the anti-BRST symmetry of this model where the role of ghost fields will be changed by anti-ghost field. The conserved charge for anti-BRST symmetry (Q) must also satisfy the Kugo-Ojima condition: The solution of this equation will give us the quantum mechanical wave function of the universe. The BRST charge obtained in (50)can be mapped easily to the Wheeler DeWitt equation obtained in [4]. Here, in this paper, in the BRST approach we have not solved the first class constraints and hence the quantization is actually done in the extended phase space. The corresponding form of the WDW potential will also depend upon the variables of extended phase space. Here we did only the Lagrangian formulation to construct the BRST charges. Whereas. there exist another process namely the BFV [26] formulation by which one can construct the Hamiltonian of the model with ghost fields in the extended phase space. In this approach the ghost field and the lagrange multipliers are treated as dynamical variable and is a gauge independent way. The Hamiltonian thus obtained enables one to directly write down the Wheeler DeWitt equation and the Wheeler DeWitt potential will be straight forward.

Finite field-dependent BRST transformation
The purpose of this section is to study the extended BRST symmetry, known as FFBRST transformation, by making the transformation parameter finite and field dependent. Such transformations have been studied in various contexts with various important motivations. We try to build such formulation first time in RT based minisuperspace models. We achieve this goal by deriving first the methodology of FFBRST transformation, originally advocated in Ref. [14], in a much elegant way. Then, we discuss its illustration part.
We start with the consideration of the fields of RT model, written collectively as φ, as a function of parameter κ : 0 ≤ κ ≤ 1 in such a manner that the original fields and finitely transformed fields are described by its extremum values. Specifically, φ(τ, κ = 0) = φ(τ ) defines the original (nontransformed) fields, however, φ(τ, κ = 1) = φ ′ (τ ) refers the finite field-dependent BRST transformed fields. Now, fields transform under an infinitesimal field-dependent BRST transformation as [14], Here s b φ refers the Slavnov variation. To get the FFBRST transformation, we first integrate the above equation w.r.t. κ from 0 to κ, leading to which at boundary (κ = 1) yields the FFBRST transformation [14], Here Θ[φ] is a finite field-dependent parameter related to infinitesimal version by Θ The remarkable feature of this FFBRST transformation is that this leaves the Faddeev-Popov action invariant, however, the functional measure is not, leading to a non-trivial Jacobian. We compute the Jacobian of path integral measure under such FFBRST transformation by following the similar procedure (up to some good extent) as discussed in [14]. The Jacobian of functional measure under an infinitesimal change is given by which further reads where ± signs are considered according to the nature of the fields φ. For instance, + is used for bosonic fields and − for Fermionic ones. Upon Taylor expansion, the relation (56) yields, This further simplifies to Upon integration w.r.t κ with its limiting values, this leads to following expression: By exponentiating the above expression, we get the explicit form of the Jacobian of functional measure as follows: Due to this Jacobian under FFBRST transformation, an effective action S[φ] of the generating functional of the theory undergoes to following change where φ ′ refers the transformed fields collectively. Here we draw a conclusion that under the whole procedure, the effective action of the theory gets a precise modification in their original expression by an extra piece. We will notice that under such an analysis the theory does not change on physical ground. The resulting action depends on a arbitrary value of Θ ′ . In the next subsection, we illustrate this result by considering an specific value of Θ ′ .

Jacobian calculation
For illustrating the above results discussed in the last section, we compute the Jacobian of functional measure explicitly. To compute the Jacobian under a specific FFBRST transformation, we first construct an infinitesimal field-dependent parameter with following form: where F [ϕ] is an arbitrary gauge-fixing condition with ϕ ≡ T, A, a, α. In construction of this parameter, we take care of the ghost number of Θ ′ , which must be −1.
Now, the expression (60) together with (62) yields This corresponds to the Jacobian for path integral measure under FFBRST transformation with parameter in (62). With this Jacobian, the functional integral (under FFBRST transformation) changes to which is nothing but the expression for generating functional of minisuperspace model in an arbitrary gauge. This assures that the FFBRST transformation maps the minisuperspace model in one specific gauge to any other arbitrary gauge. Since the gauge choice from one to another does not change the theory on physical ground as gauge-fixed action is BRST exact. Thus, we conclude that the FFBRST transformation with any particular parameter does not change the theory numerically.

Conclusions
In this paper, we have analysed the BRST symmetries of the Einstein-Hilbert(EH) action in the minisuperspace representation. Einsitein's theory includes first class constraints for which BRST symmetries play significant role during quantisation. We constrain the generalised EH action, to 4+1 dimension. So the bulk we considered is 5 dimensional while the gravity is induced on the 4D hypersurface. The cosmological constant has been considered here for the reason that five dimensional theory cannot have a vacuum solution [25]. For the homogeneous and isotropic FRLW background, the Lagrangian which was obtained appear to be as pseudo higher derivative in nature. The Lagrangian we derived with respect to the minisuperspace variables contains a total derivative term in addition to the first order term. Without avoiding the surface term we considered the higher derivative version as existence of the surface term is directly linked to the thermodynamics of the system . Dirac's constraint analysis showed that there are one primary first class constraint and one secondary first class constraint. We have seen that the model is consist of only one gauge symmetry i.e. the reparametrisation invariance. The gauge symmetry is due to the existence of the primary first class constraints. Looking at the gauge structure of the system, analogous BRST symmetries were constructed. We considered the cosmic gauge and the gauge proposed in [4] to construct the effective Lagrangian. As there is only one gauge symmetry these two constraints were added to the effective Lagrangian by a single Lagrange multiplier. The BRST transformations for all the fields in the extended phase space were calculated. Further we have calculated the BRST charges utilizing Noether theorem. The charge constructed so annihilate the physical states. The finite field-dependent BRST (FFBRST) transformations were also analyse for the theory. We have computed the Jacobian for the path integral under FFBRST transformation with an arbitrary parameter which amounts a precise change in effective action. The results were also illustrated with a particular construction of parameter which maps the action in a specific gauge to an arbitrary gauge.
The aim to construct the BRST symmetries for a higher derivative system has been completed in the previous sections. Utility of the first order formalism to extract the gauge symmetries for the higher derivative system was invoked in the present paper. In fact, apart from just a mathematical tool the method described here can be helpful for other higher dimensional models like the Randall-Sundurum background or more likely the Dvali-Gabadadze-Porrati extensions [27]. As the Randall-Sundurum(RS) model deals with the extra dimensions, one can intriguingly find connection of the model discussed in this paper. The induced metric 3 which is four dimensional here can give the results as obtained here only after consideration of the normal vector 2. So in the embedded space the choice of normal vector is crucial. Consequently, a proper normal vector will lead us to the metric discussed in [27]. One can further use these result to quantize the system and use the more tools like renormalization from quantum field theories and can see the outcomes. The analysis is left here with the scope for future consideration and as a step towards quantization in the path integral method for better understanding of the more complicated higher derivative actions.