Generalized Sparling-Thirring form in the Brans-Dicke theory

The definition of the Sparling-Thirring form is extended to the Brans-Dicke theory. By writing the Brans-Dicke field equations in a formally Maxwell-like form, a superpotential and a corresponding pseudo energy-momentum form are defined. The general energy expression provided by the superpotential in the Jordan frame is discussed in relation to the corresponding expression in the Einstein frame. In order to substantiate its formal definition, the generalized Sparling-Thirring form is used to calculate the energy for the spherically symmetric vacuum solution in the Brans-Dicke theory.


I. INTRODUCTION
According to the equivalence principle, the gravitational field, as described by Einstein as the geometry of spacetime, can locally be transformed away at a given point so that a curved spacetime looks locally like a flat Minkowskian at an infinitesimal level. Conversely, an accelerated frame in flat spacetime can emulate gravity [1]. Consequently, one can argue that the equivalence principle renders any definition of a local energy density to be nontensorial. This line of thought lead to the definition of pseudo energy-momentum forms that depend on a chosen frame. In particular, for an asymptotically flat spacetime, a working definition of energy-momentum can be provided. A well-known example of such a definition of pseudo-energy momentum form and the superenergy, in the framework of the general relativity theory (GR), is the Sparling-Thirring 2-form and the associated pseudo-energy momentum forms [1][2][3]. In GR, the Sparling-Thirring form is introduced to define the total energy and it is related to the definition of energy by Arnowitt-Deser-Misner (known as the ADM energy in short). The Sparling-Thirring form can be related to other definitions of energy, such as the well-known Landau-Lifschitz energy [2,[4][5][6].
The constructions of a pseudotensor using field equations of gravity and the definitions involving the Noether symmetries have previously been discussed in connection with the quasilocal quantities in [7][8][9][10][11]. In a completely different approach, but using a similar mathematical notation and set up, the Bell tensor is derived from a variational principle in the context of a particular tensor-tensor model of gravitation [12]. More recently, conserved Noether quantities were also studied in the context of the motion of spinning particles in the BD theory [13]. A general formulation of generalized Bianchi identities relative to conformally related frames and corresponding conservation laws in various modified gravity models have been discussed in [14].
The scalar-tensor theory of gravitation proposed by Brans and Dicke [15] almost half a century ago is still one of the popular theories of gravity which received rekindled interest in the applications in particular to the cosmological dynamics as possible alternatives to dark energy and dark matter [16,17]. The dynamical equivalence between f (R)-type modified gravity models to the scalar tensor theories, introduced in [18,19] for particular modified gravity models, also increased the interest to these theories [20][21][22]. For a treatment of this equivalence in the differential form calculus see [23]. Any modified theory of gravity should compete with the GR on the observational front at cosmological and solar system scale.
On the theoretical side, a prediction of such a modified theory is to be compared to those of GR. The conserved currents for BD theory has been studied by applying the Lagrange-Noether machinery [24][25][26]. More recently, using the exterior algebra of differential forms, and in a more general context of Noether's theorem, a new conserved charge is defined as an application in [27]. In an earlier work, the Landau-Lifshitz superpotential is extended to BD theory [28].
In the approach making use of the field equations, the definition of a pseudo-tensor can be applied to any geometrical theory that follow from any Lagrangian once the field equations formulated in terms of the tensor-valued differential forms. For example, a generalized Sparling-Thirring form was previously defined [29] for the dimensionally continued Euler-Poincare form where L EP can be considered as natural generalization of Einstein-Hilbert Lagrangian leading to second order field equations (in the notation is introduced in the following section).
In this generalization, an exact differential is separated out from the field equations following from (1) using the definition of the curvature 2-forms as in the case of the Einstein field equations.
In the present work, the Sparling-Thirring construction is extended to the BD scalartensor theory. The BD scalar field is incorporated into the definition of Sparling-Thirring form, as a consequence of the vanishing torsion constraint on the connection 1-form obtained by the constrained first order formalism to the BD theory Lagrangian.
The organization of the paper is as follows. In the next section, after briefly discussing the notation of the tensor-valued forms and the algebra of differential forms, the BD field equations are obtained by the variation of the BD action via a first order formalism. The use of the modern language of differential forms has several advantages over tensorial treatment, in that the differential forms are natural integrands on manifolds and allows one to use Stokes' theorem to convert volume integrals to flux integrals. Moreover, the use of differential forms also allow one to define energy for asymptotically flat spacetimes without introducing an approximation scheme. The definition of energy superpotential is therefore requires the BD field equations to be formulated relative to an orthonormal coframe. Due to this, we express the well known BD field equations relative to an orthonormal coframe in this section. In section 3, following the construction of Sparling-Thirring form in GR, the BD field equations in this form is used to define a superpotential and a corresponding energy momentum form which involves the nonminimally coupled dynamical BD scalar field as well. In accordance with the equivalence principle, the resulting expressions for energy and momentum depends on a chosen frame as in GR. Since the definition of BD pseudoenergy-momentum form is particularly suitable to discuss the mathematical expressions in Jordan frame and conformally related Einstein frames, we obtain relevant expressions in the Einstein frame as well in section 4. In section 5, the energy definition is used to define the total mass of the spherically symmetric and asymptotically flat vacuum BD solution. The paper ends with general remarks regarding to some possible energy definitions in related modified gravity theories in line with the discussion below.

II. BD FIELD EQUATIONS
The notation used for the algebra of differential forms is adopted from [30]. The calculations in this paper will be carried out relative to a set of orthonormal basis coframe 1-forms {e a } for which the metric reads g = η ab e a ⊗ e b with η ab = diag(− + ++). The set of basis frame fields is {X a } and the abbreviation i Xa ≡ i a is used for the contraction operator with respect to the basis frame field X a . * denotes the Hodge dual operator acting on the basis forms and * 1 = e 0 ∧ e 1 ∧ e 2 ∧ e 3 is the oriented volume element. The abbreviations of the form e a ∧ · · · ∧ e c ∧ · · · ≡ e a··· c··· for the exterior products of basis 1-forms are extensively used for the convenience of the notation. The first structure equations of Maurer-Cartan with the vanishing torsion two form, i. e., Θ a = 0. D is the covariant exterior derivative operator, acting on tensor-valued forms, and a suitable definition and its relation to covariant derivative can be found for example in [1]. The curvature 2-form Ω a b with Ω a b = 1 2 R a bcd e cd satisfies the second structure equation of the Maurer-Cartan The discussion below will be based exclusively on the original BD Lagrangian among other scalar-tensor theories. The approach to the definition of a total energy in BD theory below involves surface integrals of pseudo-tensorial quantities, therefore it is convenient to use the algebra of differential forms defined on pseudo-Riemannian manifolds that are already built for integration.
Expressed in terms of differential forms, the total Lagrangian 4-form for the original BD Lagrangian with matter fields included has the form where the gravitational part in the so-called Jordan frame is given by and the matter part L matter [g, ψ] is assumed to be independent of the BD scalar and the spinor fields. The gravitational coupling constant G in GR is replaced by a dynamical scalar field φ −1 with a corresponding kinetic term for the scalar field. ω is the free BD parameter and GR limit is recovered in the ω → ∞ limit. For the vacuum solutions or for when the matter energy momentum tensor has a non vanishing trace, however, things are more complicated and the above limit may not give [31][32][33] the GR limit, unless the arbitrary parameters are fixed by some physical grounds such as post Newtonian extension for asymptotically flat solutions [34].
In the general framework of first order formalism for gravity, the independent gravitational variables can be taken as the set of basis coframe 1-forms {e a } and the connection 1-forms {ω a b }. The local Lorentz invariance forbids any gravitational action to have explicit dependence on {ω a b } and the first order derivatives of de a and dω a b enters into a gravitational Lagrangian with local Lorentz symmetry by the tensorial quantities such as Θ a and Ω a b . The minimal coupling prescription for the matter fields implies that de a and dω a b occurs only in the gravitational sector in the presence of the matter fields. On the other hand, the BD scalar field φ couples to the metric nonminimally. Because of the nonminimal coupling, the BD scalar field is dynamical even in the absence of the kinetic term for it.
An essential ingredient of the discussion is the formulation of the BD equations in terms of exterior forms which then allows one to write the conservation laws in differential and integral forms as well easily with the help of Stokes' theorem with considerable technical advantage over the methods using tensor components. Hence, a derivation of BD equations in the desired form using exterior algebra of forms using the Lagrange multiplier method is presented below.
The vanishing torsion constraint can be implemented into the variational procedure by introducing Lagrange multiplier 4-form term to the original Lagrangian form L BD , where the Lagrange multiplier 2-form λ a is a vector- for the extended gravitational part then has the explicit form The total variational derivative of L ext. with respect to independent variables can be found as up to an omitted boundary term. The energy-momentum 3-form of the scalar field * T a Assuming that there are no spinor fields to couple to the connection, the independent connection equations read These equations can be considered as the defining equation for the Lagrange multiplier 2forms λ a and they can uniquely be solved for the Lagrange multiplier 2-forms by calculating its contractions by taking the constraint Θ a = 0 into account. One finds that Consequently, using the expression (12) for the Lagrange multiplier form in the metric field equations induced by the coframe variational derivative δL ext. /δe a ≡ − * E a = 0 in (9) read where E a ≡ E a b e b is vector-valued 1-form. In the presence of the matter fields, the BD field equations take the form * E a − 8π where * T a [ψ] stands for the energy-momentum form for the matter field ψ derived from the variational derivative of L matter [g, ψ] with respect to the basis coframe 1-forms. * T a [ψ] depends on the metric tensor as well and therefore in the notation used here, it may also be appropriate to state the dependence by writing it in the form * T a [ψ, e a ]. However, the metric dependence is surpassed for simplicity in the case that the metric dependence is obvious from the context in the discussions relative to a given coframe.
Although it is customary to write the BD equations (13) by dividing it with φ, the Lagrange multiplier and the ωφ −1 * T a [φ] terms are considered to be on the "geometry" side of the BD field equations. As a consequence of the diffeomorphism invariance of the BD Lagrangian, it follows from the corresponding Noether identity that D * E a = 0 [35].
Alternatively, in a more straightforward way, and with the help of the identities it is possible to arrive at the identity as expected. The right-hand side vanishes identically provided that the field equation for the BD scalar is satisfied since the terms on the right-hand side inside the bracket are proportional to the field equations for the BD scalar given below. This consideration is in line with the well-known case of the Einstein-massless scalar field equations, R a = (i a dφ)dφ, from which the field equations for φ follows from the Bianchi identity. Consequently, one has D * T a [ψ] = 0, from which one can derive geodesic postulate for test particles by introducing the matter energy momentum of ideal fluid [30]. As it is well-known, this is not the case in the conformally related Einstein frame, where scalar field couples nonminimally to a matter field.
On the other hand, the field equation for the BD scalar that follows from the variational derivative δL BD /δφ = 0 is given by The BD scalar couples to the matter energy momentum through the last term in (17). Then, by combining the scalar field equation with the trace of the metric equations, the equation for the BD scalar reduces to for which the trace T [ψ] ≡ T a a [ψ] of the matter energy-momentum tensor act as the source term. As pointed out in the invariance of L BD discussion above, the reduced field equation (18) follows from Bianchi identity for the BD field equations together with its trace.

MOMENTUM FORMS
The definitions of the superpotential for the BD theory can be given in terms of the original Sparling-Thirring 2-form in GR. Thus it is appropriate to recall the construction of Sparling-Thirring form [1][2][3]. The Einstein 3-form is defined in terms of the following where 1-form G a can be defined in terms of the Einstein tensor components as G a ≡ G a b e b . The definition (19) is suitable to separate out an exact differential out of Einstein tensor regardless of the presence of any matter Lagrangian. Explicitly, by inserting the second structure equation (4) into the right-hand side of (19), probably technically in the most straightforward way, one finds * G a = d * F a + * t a , with the Sparling-Thirring 2-form * F a = − 1 2 ω bc ∧ * e abc , and the gravitational energy-momentum (pseudo-tensor) 3-form * t a ≡ 1 2 (ω bc ∧ ω a d ∧ * e bcd − ω bd ∧ ω d c ∧ * e abc ).
The classical definitions of the Sparling-Thirring superpotential and the corresponding energy-momentum pseudo-tensor above can be related to well-known expressions for other superpotentials [2]. In particular, the components t ab of the pseudo-tensorial object * t a = t a b * e b are symmetrical relative to a coordinate basis [1]. The definitions (21) and (22) then help to split * E a , in an equally straightforward way as in the GR case above, in the following Maxwell-like form for BD theory * E a = d * F a + * T a = 8π with the generalized definitions of the superpotential * F a ≡ φ * F a − * (dφ ∧ e a ), (24) and the corresponding energy-momentum form for the superpotential and the pseudo-energy-momentum forms respectively for the BD vacuum case. A conserved pseudo 4-current then follows from the equations (23) as the simple consequence of the differential identity d 2 ≡ 0. The expression on the right-hand side of (25) is pseudo tensorial and it involves connection 1-forms which can be transformed away at a given point in accordance with the equivalence principle.
The formulation of the field equations in the form of a conservation law in terms of differential forms allows one to rewrite the conserved quantities as flux integrals by making use of Stokes's theorem. By defining total energy-momentum pseudo form * τ a ≡ 8π a conserved 4-current P a BD then can be expressed as the following flux integral with ∂U as the boundary of a three-dimensional submanifold U ⊂ M on a pseudo-Riemannian manifold M.
The BD superpotential * F a can be put in a simplified form analogous to the Sparling- with a modified connection 1-form denoted by Λ bc = −Λ cb that is given explicitly by the following form By definition, 1-form Λ bc incorporates the constraint on the connection 1-form resulting from nonminimal coupling of the BD scalar in a peculiar way.
The BD superpotential * F a defined in (28) naturally involves BD scalar which carries dynamical degrees of freedom and the scalar contribution to the superenergy definition results from the vanishing torsion constraint term on the independent connection. Consequently, the dynamical coupling constant in BD theory also incorporated into the total energy by definition through the Lagrange multiplier term and at the same time the construction is guided by the mathematical structure of the field equations formulated in terms of differential forms. By construction, the conserved quantities require the field equations to be satisfied.
In parallel to the original definition of the total energy in terms of Sparling-Thirring form for an asymptotically flat geometry in GR, one can define the total energy as the temporal component of conserved 4-current E BD ≡ P 0 BD as where the flux integral is over two dimensional sphere with infinite radius denoted by S 2 ∞ evaluated at a constant value of t. As for the other classical pseudo energy-momentum and superpotential forms, the 4-momentum definition (27) depends on the frame in which it is computed [7]. In the same way as the superpotentials are calculated in GR for an asymptotically flat spacetime, the BD superpotential is to be calculated in a coordinate system that is most nearly Minkowskian as well.
In contrast to the previous energy-momentum pseudo-tensors for scalar-tensor theories in the literature, the use of exterior forms in defining conserved currents in BD theory is technically straightforward which amounts to a judicious arrangement of the terms in the field equations. This arrangement of the field equations not only separates the higher order terms but also singles out the leading terms that are linear in the derivatives of the field variables.

IV. JORDAN FRAME VS. EINSTEIN FRAME
The definition of the total energy in BD theory above involves no approximation scheme exercised in other approaches, for example linearization of the field equations around a suitable background solution [37,38]. The corresponding definition in the Einstein frame below allow one to discuss the energy definition relative to Einstein frame in a straightforward manner.
It is a well-known fact that the conformal scaling of the metric components g µν (relative to a coordinate coframe with Greek indices) by brings the BD Lagrangian into the Einstein frame with a new scalar field coupled nonminimally to matter fields. The conformal transformation (31) is equivalent to the scaling of the coframe basis 1-forms as [39] e a →ẽ a = φ 1/2 e a .
Furthermore, the interior product operators of the frames are related by i Xa = φ 1/2 iX a , whereas the invariant volume forms are related by * 1 = φ −2 * 1. Consequently, the connection 1-forms transform as The conformal transformation (32) brings BD Lagrangian into the so-called Einstein frame with a minimally-coupled massless scalar field α as up to an omitted closed form. The scalar field α is related to φ by α = (ω + 3/2) 1/2 ln φ.
The metric field equations that follow from the scaled Lagrangian can be obtained as where the conformally scaled metric is defined asg = e −2bα g and the constant b is defined in terms of the BD parameter ω as b = 2 2ω + 3 In the Einstein frame the gravitational part assumes the familiar form while coupling to the matter becomes nonminimal while in the Jordan frame coupling to the matter is through a dynamical BD scalar field. As a consequence, contrary to a matter energy-momentum forms relative to the Jordan frame, a matter energy-momentum is not covariantly constant in the Einstein frame. On the other hand, the construction of the expression for the BD superpotential is not modified by a matter field Lagrangian coupled to the BD Lagrangian.
The familiar Sparling-Thirring superpotential form can be adopted for the field equations (36) defined relative to the Einstein frame. The above definition of BD superpotential (28) facilitates the expressions in the Jordan frame. By applying the conformal transformation used above to the generalized Sparling-Thirring form * F a in the Jordan frame * F a = − whereω ab is given by (33). Then one immediately finds Finally, by using the defining relations for the conformal transformations (32), one ends up with the remarkable result * F a = φ 1/2 * F a , whereF a is the Sparling-Thirring form in the Einstein frame of the formF a = −1/2ω bc ∧ * ẽ abc similar to (21). This equation relates the Sparling-Thirring form in the Einstein frame and the generalized BD superpotential defined in the corresponding Jordan frame.

VACUUM METRIC
In testing the validity of the formal definition of a conserved charge that makes explicit use of the field equations, an exact solution to the field equations provides a valuable tool in explicit applications. In the discussion below, a particular spherically symmetric, static, vacuum solution to the BD theory will be taken into account. In the so-called isotropic coordinates, the metric for the spherically symmetric vacuum solution to the BD theory has the form [15] with the explicit forms of the metric functions given by The metric functions depend only on the function r defined by r 2 = δ ij x i x j in terms of the the spatial coordinates {x i }, for i, j = 1, 2, 3. α 0 , β 0 , B are integration constants and the remaining constants C and λ are related further by As a function of r, the BD scalar is given by The The explicit expressions for the Levi-Civita connection 1-forms relative to the orthonormal coframe defined by e 0 = f dt and e k = h dx k are given by where ′ stands for the derivative with respect to r. The expansion of connection 1-forms into the associated coordinate basis coframe forms in (46)-(47), rather then the orthonormal coframe is convenient in order to evaluate the flux integral in what follows.
Moreover, with the assumption φ = φ(r), one finds that the modified connection 1-form Λ jk has the explicit expression Consequently, the expression in (48) can be used in the definition (27) to calculate the total energy E BD as where ǫ 0ijk = −ǫ 0ijk and ǫ ijk stands for the permutation symbol on the spatial submanifold t = constant. By considering (48), the total energy becomes The integration measure above can be put into a form convenient for the spherical symmetry by recalling the basic geometrical formula for S 2 on Euclidean space R 3 . On a three dimensional Euclidian space R 3 , the unit S 2 can simply be defined by setting r = 1. In terms of the Cartesian coordinates {x i }, one has r 2 = δ ij x i x j . By differentiating this expression and taking the Hodge dual one finds where ⋆ is the Hodge dual in R 3 . Restriction of the expression on the right-hand side to S 2 r gives the volume element on S 2 r which is conveniently taken to be r 2 dΩ with dΩ standing for an infinitesimal element of the solid angle. Eventually, the energy expression becomes For φ ∼ φ 0 =constant, the expression on the right-hand side of Eq. (52) becomes proportional to the corresponding expression in GR (see, for example, the expression given in [1]) up to a constant.
Evaluating (52) for the Brans-I solution defined by Eqs. (41)-(45) above, one finds The expression on the right-hand side contains five parameters, but due to the relation (44), only four of them are free.
Let us compare this expression with corresponding expression of GR. By considering the linear expansion of the static vacuum BD solution given in (41) and the scalar field (45), and subsequently matching it to a static Newtonian source, the parameters B, C, λ, φ 0 , β 0 can be expressed in terms of the corresponding GR parameters together with the BD parameter ω as [15,41] Note that in our notation we have set G = 1. Putting these into (53), we have found that the Einstein limit of the total energy is given by It is possible to obtain this result in a more straightforward way as well. The choice C = 0, λ = 1 brings the Brans-I solution to the corresponding one in GR, which yields the result in the Einstein limit defined by ω → ∞. In the expression above, B E is a constant and the particular choice B E = M/2 yields the correct GR limit.

VI. CONCLUDING REMARKS
Both Abbott-Deser [36] and Deser-Tekin [37] charge definitions make essential use of the field equations as well as the Killing symmetries of the background solution. For example, the Deser-Tekin charges are constructed by linearizing the field equations around a background solution (flat Minkowski space in our case above). Subsequently, a conserved quantity is expressed as a flux integral in the background by making use of the linearized equations.
The construction of the generalized Sparling-Thirring superpotential is achieved in the same spirit as the construction of a Deser-Tekin charge and the use of exterior forms and Stokes's theorem allows one to obtain flux integrals on background spacetime in a practically useful form see, for example, the construction in [38] by making use the Killing vector fields of the background spacetime. As for the BD example above, the energy in the Deser-Tekin approach is simply obtained by linearizing the generalized superpotential * F a .
As another important technical point, note that the construction of the BD superpotential is insensitive to a possible potential term for the scalar field. The definition also covers the special case where the kinetic term for the BD scalar is absent and therefore the construction can be applied to more general f (R) models since they are known to have scalar-tensor equivalent models with a potential term for the scalar field. To facilitate a comparison with these theories, let us consider the simplest modification of the Einstein-Hilbert Lagrangian with f (R) assumed to be an arbitrary differentiable and nonlinear function of the scalar curvature R. Particular forms of the function f are studied with different motivations, for example, arising from cosmological applications.
The metric field equations for vacuum that follow from the modified Lagrangian above take the form [43] f ′ * G a − D * (df ′ ∧ e a ) + 1 2 (Rf ′ − f ) * e a = 0.
By comparing the field equations with the BD field equation (13) above one can see that the field redefinition f ′ ≡ φ which entails the potential term V (φ) by the Legendre trans- for the scalar field φ. Consequently, in terms of the scalar field φ, the scalar-tensor equivalent equations for (59) become φ * G a − D * (dφ ∧ e a ) + 1 2 V (φ) * e a = 0, which are in fact the vacuum BD equations with ω = 0 and with the potential term for scalar field. The form of the field equations (59) is not very common in the vast literature involving various f (R) models and reader is referred to [44] in relating (59) to the more familiar coordinate expression. In either form, the metric equations allow one to define gravitational energy in parallel to the BD case presented above immediately. By using the form of the field equation given in (59), one finds the total energy by the flux integral For the simple case f (R) = R the above formula gives back the Sparling-Thirring form up to a constant multiple. For the next simple case with f (R) = R + αR 2 the leading term is the second term on the right-hand side and the above formula reduces up to a constant multiple to the energy definition of Deser-Tekin for a particular subcase of general quadratic curvature gravity in four dimensions [37]. It follows by definition that the energy has the expression for the R 2 term is of the form for asymptotically flat solutions [45]. Further scrutiny of the generalized Sparling-Thirring forms for modified gravity models in relation to the definitions of the Deser-Tekin charges will be taken up elsewhere.